一种过渡带陡度可控的FIR半带滤波器设计方法与流程

一种过渡带陡度可控的fir半带滤波器设计方法
技术领域
1.本发明主要涉及滤波技术领域，具体涉及一种过渡带陡度可控的fir半带滤波器设计方法。


背景技术：

2.相量测量单元的广泛应用为实现智能电网的动态安全监控提供了数据基础，其动态行为是系统安全控制的关键。然而经分析，在系统动态情况下,传统的离散傅里叶变换算法会因频谱泄漏产生误差，引起振荡分析错误；另一方面，在广域测量系统中，子站pmu将相量上传主站时，有可能发生频率混叠，影响基于pmu的各种应用，数字滤波器的应用可以在一定程度解决上述问题。
3.最大平坦滤波器maxflat是一种重要的非递归有限脉冲响应fir滤波器，适用于需要高阻带衰减或平滑频率响应的情况。设计maxflat的fir滤波器的基本思想是使用一个数学证明的闭合解，该解在频带末端满足maxflat约束，并将其映射到传递函数中，计算滤波器的系数。然而，经典的设计涉及到用一个合适的闭型多项式来近似期望的频率响应，因为这种闭型解主要关注滤波器的平坦性，而不是精确的频率响应。目前已经提出了几种设计maxflat的fir半带滤波器的方法和实现技巧。这些滤波器具有精确的截止频率在频带的中间ω＝π/2，并允许计算效率的实现，因为几乎一半的脉冲响应ir系数是零。然而，它们的过渡带比其他滤波器相对更宽，只能通过增加滤波器的长度来缩小。
4.大多数流行的maxflat的fir半带滤波器是由拉格朗日半带多项式lhbp设计的，lhbp在z＝-1处具有最大的零数。这类滤波器有许多应用，如滤波器组、基于小波的压缩和多速率技术。然而，与现有的maxflat的fir半带滤波器相似，lhbp滤波器也没有任何独立参数。换句话说，没有直接控制lhbp滤波器的频率响应，以实现在阻带衰减和过渡带宽度之间的合理权衡。这是由于在z＝-1处的最大可能的零数是施加在半带封闭形式解上的，它没有留下自由度，因此没有独立的参数。
5.目前国内外的设计方法主要是实现窄带过渡滤波器，而关于实现窄带过渡带宽的半带滤波器的设计中，经常会出现宽频带问题，它不允许任意的自由度，即独立的参数，通过某种封闭的多项式来控制频率响应，因此，更需要具有可控频率特性的高精度滤波器，即在过渡带锐度和通带和阻带波纹之间进行合理的权衡。
6.为了解决目前fir滤波器设计中过渡带狭窄的问题，人们在频响控制方向上做了大量的工作。锐利过渡fir滤波器使用正弦频率函数来评估封闭形式的脉冲响应系数。该方法允许封闭参数的简单直接计算，但在阻带和通带会出现不可忽略的振幅畸变。基于智能元启发式的迭代方法使用多适应度函数结合加权误差函数，但是这种滤波器设计的缺点是需要进行复杂和大量的计算来调整频带的波纹和过渡频带的宽度，且存在半带滤波器设计的局限性。


技术实现要素：

7.本发明要解决的技术问题就在于：针对现有技术存在的问题，本发明提供一种具有明显的过渡带宽控制的fir半带滤波器的设计方法。
8.为解决上述技术问题，本发明提出的技术方案为：
9.一种过渡带陡度可控的fir半带滤波器设计方法，包括步骤：
10.步骤一：构建ii型一般对称fir半带滤波器；
11.步骤二：构建一般拉格朗日半带多项式lhbp；
12.步骤三：引入过渡带陡度参数h0，构建广义拉格朗日半带多项式g-lhbp；
13.步骤四：构建g-lhbp的目标控制函数：过渡带陡度参数h0，得到fir半带滤波器的频率响应。
14.优选地，步骤一的具体过程为：
15.设一个ii型的一般对称fir半带滤波器h(z)，其脉冲响应hn为4k-2阶，表示为h(z)＝z-(2k-1)
qk(z)，传递函数
16.式(1)所示的maxflat的fir半带滤波器通过使用封闭形式的多项式，然后通过变换映射到滤波器函数中进行构建。
17.优选地，步骤二的具体过程为：
18.在qk(z)的设计中，当qk(z)在z＝-1处为零，即1 z-1
形式的项，设计4k-2阶maxflat的fir半带滤波器使用拉格朗日半带多项式lhbp，表达式为：
[0019][0020]
其中，l＝0,1,2,...,k，
[0021]
lhbp在z＝-1处有最大的零数，因此在ω＝π处有最大的平坦响应，即：
[0022][0023]
式(2)的lhbp滤波器不具有任何独立参数，并且对lhbp获得的滤波器的频率响应没有直接的控制。
[0024]
优选地，步骤三的具体过程为：
[0025]
当qk(z)在z＝-1处存在2(k-1)个零点，即：
[0026][0027]
将其条件应用于式(1)，并在重合点处使用拉格朗日插值，使qk(z)具有与式(2)相同的递归关系，用h0得到一个封闭的半带多项式，即广义拉格朗日半带多项式g-lhbp：
[0028][0029]
当过渡带陡度参数h0为：
[0030][0031]
广义拉格朗日半带多项式g-lhbp过渡带陡度参数h0是通过在式(6)上z＝-1处额外加一个零得到；对于闭式表达式，将式(6)映射到：
[0032][0033]
根据h0和k获得插值系数g
l
：
[0034][0035]
其中，l＝0,1,2,...,k，以及获得ir系数hn：
[0036][0037]
或者表示为：
[0038][0039]
当a《b或b《0时，
[0040]
式(9)代入式(10)时奇数索引系数为0，如式(11)所示，即h
2n-1
＝0；根据式(6)，滤波器的频率响应由k和h0表示：
[0041][0042]
根据式(12)，引入h0作为参数保持频率响应可控。
[0043]
优选地，步骤四的具体过程为：
[0044]
a)由递归关系计算q
k-1
(ω)
lhbp
和峰值频率
[0045]
b)选择确定过渡带陡度参数h
0,γ
。
[0046]
优选地，在步骤四之后，还包括对过渡带陡度参数h
0,γ
进行性能评估。
[0047]
优选地，性能评估的具体过程为：
[0048]
1)推导关系式，计算最大超调波推导不等式关系；
[0049]
2)设计两个g-lhbp滤波器；
[0050]
3)对比不同g-lhbp滤波器的相关参数。
[0051]
本发明还公开了一种过渡带陡度可控的fir半带滤波器设计系统，包括：
[0052]
第一程序模块，用于构建ii型一般对称fir半带滤波器；
[0053]
第二程序模块，用于构建一般拉格朗日半带多项式lhbp；
[0054]
第三程序模块，用于引入过渡带陡度参数h0，构建广义拉格朗日半带多项式g-lhbp；
[0055]
第四程序模块，用于构建g-lhbp的目标控制函数：过渡带陡度参数h0，得到fir半带滤波器的频率响应。
[0056]
本发明进一步公开了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序在被处理器运行时执行如上所述方法的步骤。
[0057]
本发明还公开了一种计算机设备，包括存储器和处理器，所述存储器上存储有计算机程序，所述计算机程序在被处理器运行时执行如上所述方法的步骤。
[0058]
与现有技术相比，本发明的优点在于：
[0059]
本发明针对现有半带滤波器设计中存在没有直接控制频率响应的独立参数的不足，提供了一种具有明显的过渡带宽控制的fir半带滤波器的设计方法。首先，构建一个广义拉格朗日半带多项式g-lhbp，对于给定阶滤波器，它的所有ir系数用第0个ir系数h0表示；然后，通过分析g-lhbp的线性递归关系，将h0参数化，直接控制该滤波器类型的过渡带陡度或宽度；最后，提供显式公式用于直接计算与选择所需滤波器特性相关的设计参数，通过在过渡带锐度和通带与阻带平坦度之间的合理权衡。此外，为最终滤波器的性能评估提供明确的公式，使得选择具有所需幅度响应特性的最佳滤波器时，即在过渡频带锐度和通带与阻带平坦度之间进行权衡，具有不同寻常的灵活性。
附图说明
[0060]
图1为本发明的方法在实施例的流程图。
具体实施方式
[0061]
以下结合说明书附图和具体实施例对本发明作进一步描述。
[0062]
如图1所示，本发明实施例的过渡带陡度可控的fir半带滤波器的设计方法具体步骤如下：
[0063]
步骤一：构建ii型一般对称fir半带滤波器；
[0064]
设一个ii型的一般对称fir半带滤波器h(z)，其脉冲响应hn为4k-2阶，滤波器的系统函数h(z)＝z-(2k-1)
qk(z)，传递函数为：
[0065][0066]
其中k为自然数，z是一个复变量，z-(2k-1)
是对信号进行z变换的一种运算因子；式(1)所示的maxflat的fir半带滤波器可通过使用合适的封闭形式的多项式，然后通过变换
映射到滤波器函数中进行设计。
[0067]
步骤二：构建一般拉格朗日半带多项式lhbp；
[0068]
在qk(z)的设计中，当qk(z)在z＝-1处为零，即1 z-1
形式的项，设计4k-2阶maxflat的fir半带滤波器使用拉格朗日半带多项式lhbp，表达式为：
[0069][0070]
其中，l＝0,1,2,...,k，
[0071]
lhbp在z＝-1处有最大的零数，因此，它在ω＝π处有最大的平坦响应，即：
[0072][0073]
其中为微分符号，表示对函数q求k阶偏导数；
[0074]
式(2)所述的lhbp滤波器不具有任何独立参数，并且对lhbp获得的滤波器的频率响应没有直接的控制。
[0075]
步骤三：引入过渡带陡度参数h0，构建广义拉格朗日半带多项式g-lhbp；
[0076]
当qk(z)在z＝-1处存在2(k-1)个零点，即：
[0077][0078]
将其条件应用于式(1)，并在重合点处使用拉格朗日插值，使qk(z)具有与式(2)相同的递归关系，可以用h0得到一个封闭的半带多项式，即广义拉格朗日半带多项式g-lhbp。
[0079][0080]
当过渡带陡度参数
[0081]
广义拉格朗日半带多项式g-lhbp过渡带陡度参数h0是通过在式(6)上z＝-1处额外加一个零得到。对于一般的闭式表达式，将式(6)映射到：
[0082][0083]
其中gk和g
k-l
为关于k的系数；
[0084]
根据h0和k获得插值系数g
l
：
[0085][0086]
其中，l＝0,1,2,...,k，以及获得ir系数hn：
[0087][0088]
或者表示为：
[0089][0090]
当a《b或b《0时，
[0091]
式(9)中代入式(10)时奇数索引系数为0，如式(11)所示，即h
2n-1
＝0。
[0092]
根据式(6)，滤波器的频率响应由k和h0表示：
[0093][0094]
根据式(12)，引入h0作为参数保持频率响应可控。
[0095]
步骤四：构建g-lhbp的目标控制函数：过渡带陡度参数h0；
[0096]
(a)由递归关系计算滤波器的频率响应q
k-1
(ω)
lhbp
和峰值频率
[0097]ak
(z)作为q
k-1
(z)
lhbp
到qk(z)的目标函数，滤波器的频率响应qk(ω)为：
[0098]
qk(ω)＝q
k-1
(ω)
lhbp
ak(ω)
ꢀꢀꢀ
(13)
[0099]
其中q
k-1
(ω)
lhbp
和ak(ω)分别为：
[0100][0101]
将ak(ω)简化为根据式(13)，对于给定k，由于q
k-1
(ω)是已知的maxflat低通函数，ak(z)用h0表征qk(ω)；从式(14)开始，ak(ω)经过ω＝0，ω＝π/2，ω＝π处的零值，在0≤ω≤π范围内，以π/2为中心呈反对称正弦波状。因此，有两个反对称峰值可用于识别ak(ω)在qk(ω)上h0的影响观测值。和是ak(ω)的两个峰值频率，可以通过在时，通过解可以得到：
[0102][0103]
将式(15)代入式(14)，得到和两个峰值：
[0104]
[0105]
其中，(双符号同序)表示和与的有关。由于q
k-1
(z)
lhbp
和ak(ω)作为0≤ω≤π的中心点，对ω＝π/2都是反对称的，因此滤波器的频率响应和峰值满足：
[0106][0107]
由式(13)可得，和满足：
[0108][0109]
(b)选择确定过渡带陡度参数h
0,γ
[0110]
用将式(16)代入式(13)，并对h0化简，可得
[0111][0112]
利用和作为该滤波器的两个过渡带边缘频率，即上下边缘频率，可以得到这类g-lhbp过滤器的斜率如下：
[0113][0114]
其中通过式(15)的积化和差变换得到。
[0115]
从式(20)可以看出，式(19)中所示的h0可用作陡度参数，以直接控制给定的k用和表示过渡带斜率。因此，令上边缘频率下边缘频率然后根据γ改写h0和slopek：
[0116][0117][0118]
为了根据式(21)确定h
0,γ
，使g-lhbp滤波器比maxflatg-lhbp滤波器有一个相对较窄的过渡带，γ必须在的范围内选择。其中，可根据式
(8)、式(13)、式(15)计算得到。对于给定的k，将代入式(21)计算得到的h
0,γ
与由式(8)计算得来的h0相等。而将h0与γ》1结合计算得到的qk(ω)在通带和阻带有更尖锐的过渡带和更大的畸变响应。
[0119]
本发明针对现有半带滤波器设计中存在没有直接控制频率响应的独立参数的不足，提供了一种具有明显的过渡带宽控制的fir半带滤波器的设计方法。首先，构建一个广义拉格朗日半带多项式g-lhbp，对于给定阶滤波器，它的所有ir系数用第0个ir系数h0表示；然后，通过分析g-lhbp的线性递归关系，将h0参数化，直接控制该滤波器类型的过渡带陡度或宽度；最后，提供显式公式用于直接计算与选择所需滤波器特性相关的设计参数，通过在过渡带锐度和通带与阻带平坦度之间的合理权衡。此外，为最终滤波器的性能评估提供明确的公式，使得选择具有所需幅度响应特性的最佳滤波器时，即在过渡频带锐度和通带与阻带平坦度之间进行权衡，具有不同寻常的灵活性。
[0120]
下面结合附图对本发明的具体过程进行详细说明：
[0121]
首先，设计过渡带陡度可控的fir半带滤波器，过程为：
[0122]
(a)给定k，构建g-lhbp滤波器一般形式qk(z)；
[0123]
给定k＝4的情况下，即滤波器的阶数为4k-2＝14，构建14阶g-lhbp的一般形式：
[0124]
(b)计算插值系数g
l
(l＝0,1,2,3,4)：
[0125][0126]
得到传递函数q4(z)，表示为：
[0127][0128]
(c)得到系数h
2n
和h
2n-1
，计算频率响应q4(ω)；
[0129]
半带滤波器的奇数索引系数为零，即h
2n-1
(n＝1,2,3)＝0；h
2n
(n＝1,2,3)可以通过将插值系数替换式(11)中的k＝4来获得。通过q4(ω)＝q3(ω)
lhbp
a4(ω)得到频率响应：
[0130]
(d)构建g-lhbp的目标控制函数：过渡带陡度参数h0；
[0131]
首先，计算q
k-1
(ω)
lhbp
和峰值频率
[0132]
得
[0133]
得
[0134]
选择确定过渡带陡度参数h
0,γ
；
[0135]
选择γ＝maxflat＝0.8592代入计算式中，得到过渡带陡度参数h
0,γ
＝h0＝-d
4,3
/2
14
＝-0.00122；选择γ＝1.0，代入确定过渡带陡度参数h
0,γ
＝-0.0584,因此得到滤波器的频率响应q4(ω)。
[0136]
γ＝maxflat和γ＝1.0情况下两个g-lhbp滤波器的g
l
和hn，如表1所示；
[0137]
表1 k＝4时两个g-lhbp滤波器的系数
[0138][0139]
然后评估过渡带陡度参数h
0,γ
性能，并与其他滤波器进行性能比较。选择h
0,γ
和h0在范围内，会分别在通带和阻带中发生过调和下调。
[0140]
具体地，对过渡带陡度参数h
0,γ
性能评估，过程如下：
[0141]
(1)推导关系式，计算最大超调波推导不等式关系；
[0142]
假定k＝3，且超调频率和欠调频率和欠调频率和分别在通带和阻带内达到峰值。在时，可以得到由k和h
0,γ
得到得到
[0143]
由于q
k-1
(z)
lhbp
和ak(ω)作为0≤ω≤π的中心点，对ω＝π/2都是反对称的，因此和满足性质和则得到关系式：
[0144][0145]
因此，最大超调波计算公式为：
[0146]
其中，当选择h
0,γ
＝0设计具有maxflat响应的g-lhbp滤波器时，根据式子和推导出不等式关系如下所示：
[0147][0148]
(2)设计两个g-lhbp滤波器；
[0149]
在γ＝1.0和γ
maxflat
＝0.8667两种情况下计算h
0,γ
并设计两个g-lhbp滤波器。结果表明，使用h
0,γ
可使g-lhbp滤波器具有窄过渡带，但在通带和阻带中会出现失真，如超调和欠调，因此，仍然存在由于使用h
0,γ
导致的这种不希望出现的失真是否在此类滤波器的设计可接受的限值范围内，具有容许的幅度失真但过渡带较窄。
[0150]
(3)对比不同的g-lhbp滤波器的相关参数；
[0151]
假定k＝2、3、5、6，在γ＝0.9、γ＝1.0和γ
maxflat
三种种情况下计算h
0,γ
设计不同的g-lhbp滤波器，对比相关参数，相关参数如下表2所示。
[0152]
表2 γ
maxflat
，γ＝0.9和γ＝1.0时的相关参数
[0153][0154]
结果表明：在给定的k条件下，随着γ的增加，过渡带斜率的陡度迅速增加，但振幅畸变略有增加。k增加越多，产生的影响越大。上述方法可导出具有可控频率特性的灵活fir半带滤波器，即在过渡带锐度和通带及阻带波纹之间进行合理的权衡。
[0155]
本发明采用不同的过渡带陡度参数设计了g-lhbp低通滤波器，验证了可控的过渡带陡度参数对于滤波器设计的有效性，通过控制过渡带陡度参数，g-lhbp低通滤波器在可接受范围内可以容忍失真，过渡带相对较窄，滤波器具有近似平坦的幅度响应，使得半带滤波器设计具有可控频率特性，与先进方法进行对比，能很好地在过渡带锐度和通带和阻带波纹之间进行合理的权衡，具有不同寻常的灵活性，为明确控制频率响应提供了一些独立参数的物理意义。
[0156]
本发明通过分析g-lhbp的线性递归关系，将h0作为过渡带陡度控制因子h0参数化，直接控制该滤波器类型的过渡带陡度(或宽度)，实现灵活选择具有所需幅度响应特性的最佳滤波器，是现有技术不具备的。
[0157]
本发明实施例还提供了一种过渡带陡度可控的fir半带滤波器设计系统，包括：
[0158]
第一程序模块，用于构建ii型一般对称fir半带滤波器；
[0159]
第二程序模块，用于构建一般拉格朗日半带多项式lhbp；
[0160]
第三程序模块，用于引入过渡带陡度参数h0，构建广义拉格朗日半带多项式g-lhbp；
[0161]
第四程序模块，用于构建g-lhbp的目标控制函数：过渡带陡度参数h0，得到fir半带滤波器的频率响应。
[0162]
本发明的设计系统与上述设计方法相对应，同样具有如上设计方法所述的优点。
[0163]
本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序在被处理器运行时执行如上所述方法的步骤。本发明实施例还提供了一种计算机设备，包括存储器和处理器，所述存储器上存储有计算机程序，所述计算机程序在被处理器运行时执行如上所述方法的步骤。
[0164]
以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。