首页 理论教育 线性相位有限冲激响应系统的网络结构优化

线性相位有限冲激响应系统的网络结构优化

时间:2023-06-26 理论教育 版权反馈
【摘要】:下面从式出发推导线性相位有限冲激响应滤波器结构。利用线性相位系统的对称条件,在实现上可以减少乘法次数,从而简化了结构。图4.4.5 线性相位FIR系统的零点分布4)zi是实数又在单位圆上,则四点合一,即z4点。由上面分析可知,线性相位有限冲激响应系统的系统函数H,用级联形式结构实现时,可以分别由一阶、二阶、四阶子系统级联而成,并且每一子系统都是线性相位的,由此级联而成的整个系统也必保持线性相位的特性。

线性相位有限冲激响应系统的网络结构优化

在许多实际应用,如图像处理中,要求数字滤波器具有线性相位。具有线性相频特性的滤波器系统函数H(ejω)为

978-7-111-42877-0-Chapter04-29.jpg

式中,θω)=-τω为滤波器的相位响应函数,是ω的线性函数(即线性相频特性);978-7-111-42877-0-Chapter04-30.jpg是滤波器的幅度响应函数,是ω的实函数。

如果有限冲激响应滤波器的单位抽样响应h[n]为实数(0≤nN-1),且满足

h[n]=±h[N-1-n] (0≤nN-1) (4.4.2)

式中,“+”号代表偶对称,“-”号代表奇对称,则此时的有限冲激响应滤波器是线性相位的。式(4.4.2)说明h[n]对(N-1)/2是偶对称或奇对称的。下面从式(4.4.2)出发推导线性相位有限冲激响应滤波器结构。

h[n]=h[N-1-n],N为偶数时,其频率特性为

978-7-111-42877-0-Chapter04-31.jpg

对式(4.4.3)中第二个求和项令m=N-1-n,则得

978-7-111-42877-0-Chapter04-32.jpg

z=ejω代入式(4.4.4),得

978-7-111-42877-0-Chapter04-33.jpg

式中

978-7-111-42877-0-Chapter04-34.jpg

由式(4.4.5)可以得出其中相移常数α=(N-1)/2,在给定序列长度N时,α是一常数,因而系统的群延时

978-7-111-42877-0-Chapter04-35.jpg

τ为常数,说明该系统不存在线性失真,也就是说信号的不同频率分量通过系统时其延时相等。

从式(4.4.4)可知h[n]偶对称,且N为偶数时的网络结构形式,如图4.4.1所示。从图可以看出所需的乘法次数是N/2(注意:图4.3.1的结构需要N次乘法)。

h[n]=h[N-1-n],N取奇数,其频率特性为

978-7-111-42877-0-Chapter04-36.jpg

对式(4.4.9)中最后一个求和项令m=N-1-n,则得

978-7-111-42877-0-Chapter04-37.jpg

978-7-111-42877-0-Chapter04-38.jpg

图4.4.1 偶对称h[n](N=偶数)及其结构流图

将式(4.4.2)代入式(4.4.10),并化简得

978-7-111-42877-0-Chapter04-39.jpg

根据式(4.4.11)画出的直接结构形式如图4.4.2所示。由图可以看出所需的乘法次数为(N+1)/2,而不像图4.3.1所示的结构那样需要N次乘法。

978-7-111-42877-0-Chapter04-40.jpg

图4.4.2 偶对称h[n](N=奇数)及其结构流图

同理可以证明,h[n]满足奇对称h[n]=-h[N-1-n]条件时,对应的频率特性为

978-7-111-42877-0-Chapter04-41.jpg

式中

978-7-111-42877-0-Chapter04-42.jpg

式(4.4.12)与式(4.4.5)相比,相移常数中增加了常数相移π/2,但它的群延时仍等于α,因而还是属于线性相位系统。

满足式(4.4.1)所示的奇对称条件的h[n]如图4.4.3a和图4.4.4a所示,其中图4.4.3a是N=偶数的情况,图4.4.4a是N=奇数的情况。它的直接形式结构图与具有偶对称条件的h[n]的结构相似,只是将后者的延时相加改换成了延时相减。另外,在N为奇数时,满足978-7-111-42877-0-Chapter04-43.jpg,因而978-7-111-42877-0-Chapter04-44.jpg,即978-7-111-42877-0-Chapter04-45.jpg项不存在。具体流图结构如图4.4.3b和图4.4.4b所示。

利用线性相位系统的对称条件,在实现上可以减少乘法次数,从而简化了结构。线性相位有限冲激响应系统除去直接形式外,同样存在级联形式的结构,和一般的有限冲激响应系统相似,对多项式Hz)进行因式分解就可以得到相应的级联形式,但考虑到线性相位h[n]的对称性,Hz)的零点分布还必须满足一定的条件。

978-7-111-42877-0-Chapter04-46.jpg

图4.4.3 奇对称h[n](N=偶数)及其结构流图

线性相位有限冲激响应滤波器的零点分布特点如下:

h[n]=h[N-1-n]时

978-7-111-42877-0-Chapter04-47.jpg

图4.4.4 奇对称h[n](N=奇数)及其结构流图

978-7-111-42877-0-Chapter04-48.jpg

m=N-1-n,则

978-7-111-42877-0-Chapter04-49.jpg

h[n]=-h[N-1-n]时,同理可得

Hz)=-z-(N-1)Hz-1) (4.4.17)

由式(4.4.16)、式(4.4.17)可知,线性相位有限冲激响应滤波器的系统函数具有如下的零点分布特点。

1)若z=zi是零点,Hzi)=0,则它的倒数z=zi-1也必定是零点。这是因为

Hzi-1)=±zN-1)iHzi)=0

2)由于h[n]是实数,故Hz)的零点必须以共轭对出现,所以z=ziz=(zi-1也一定是零点。

这种互为倒数的共轭对有以下几种可能的情况(见图4.4.5):

1)zi不在实轴上又不在单位圆上,此时必然是四个互为倒数的两组共轭对,四个零点z1978-7-111-42877-0-Chapter04-50.jpgz1978-7-111-42877-0-Chapter04-51.jpg同时出现。(www.xing528.com)

2)zi是实数又不在单位圆上,其共轭也就是zi本身,则两个零点z2978-7-111-42877-0-Chapter04-52.jpg同时出现。

3)zi是单位圆上的复零点,其共轭倒数就是其本身,则两个零点z3z3同时出现。

978-7-111-42877-0-Chapter04-53.jpg

图4.4.5 线性相位FIR系统的零点分布

4)zi是实数又在单位圆上,则四点合一,即z4点。

按照这四种零点分布情况就可以知道,线性相位级联形式的结构分别由四阶、二阶和一阶子系统串联而成。位于单位圆内的复零点z1及其共轭、镜像反演对称的零点978-7-111-42877-0-Chapter04-54.jpgz1978-7-111-42877-0-Chapter04-55.jpg可以构成一个四阶实系数的子系统,即有

H1z)=(1-z1z-1)(1-z1z-1)(1-z-1/z1)(1-z-1/z1

978-7-111-42877-0-Chapter04-56.jpg,则

978-7-111-42877-0-Chapter04-57.jpg

978-7-111-42877-0-Chapter04-58.jpg

则得

H1z)=1+bz-1+cz-2+bz-3+z-4 (4.4.18)

从式(4.4.18)可以看出,H1z)本身也是线性相位的,图4.4.6画出了H1z)的线性相位直接形式结构流图。对式(4.4.18)所示的四阶子系统也可以用二阶子系统级联来实现,如图4.4.7所示。

978-7-111-42877-0-Chapter04-59.jpg

图4.4.6 四阶线性相位子系统

978-7-111-42877-0-Chapter04-60.jpg

图4.4.7 二阶子系统的级联

对于位于实轴上而又不在单位圆上的实零点z2=r2,构成的子系统是二阶的,即

978-7-111-42877-0-Chapter04-61.jpg

式(4.4.19)所描述的子系统显然是线性相位的,其网络结构如图4.4.8所示。

978-7-111-42877-0-Chapter04-62.jpg

图4.4.8 位于实轴上而又不在单位圆上的实零点组成的子系统流图

位于单位圆上而又不在实轴上的复零点978-7-111-42877-0-Chapter04-63.jpg978-7-111-42877-0-Chapter04-64.jpg所构成的子系统为

978-7-111-42877-0-Chapter04-65.jpg

978-7-111-42877-0-Chapter04-66.jpg

其网络结构与图4.4.8相同,只是把系数978-7-111-42877-0-Chapter04-67.jpg换成-2cosθ3即可。

位于单位圆上又在实轴上的实零点z4=±1,则构成的一阶子系统可表示为

H4z)=1∓z-1 (4.4.21)由此可见,当N为偶数时,Hz)有奇数个零点,其中含z=±1的零点必为奇数个;当N为奇数时Hz)有偶数个零点,如果含有z=±1的零点,则必是偶数,否则就没有z=±1的零点。

由上面分析可知,线性相位有限冲激响应系统的系统函数Hz),用级联形式结构实现时,可以分别由一阶、二阶、四阶子系统级联而成,并且每一子系统都是线性相位的,由此级联而成的整个系统也必保持线性相位的特性。

【例4.4.1】 已知一有限冲激响应系统的系统函数为

Hz)=1-2.05z-1+3.2025z-2-1.05z-3-1.05z-4+3.2025z-5-2.05z-6+z-7

求用级联形式实现的结构流图及零点分布。

由于Hz)中的N=8,因而必定有一个零点在z=±1处,由于最高阶z-7的系数为正,因而在z=-1点处有一个零点,由此可得H4z)=1+z-1Hz)除以H4z)后余下的多项式为

978-7-111-42877-0-Chapter04-68.jpg

H′z)的最高阶次N-1=6,因而有两个零点落在单位圆上,即978-7-111-42877-0-Chapter04-69.jpg978-7-111-42877-0-Chapter04-70.jpg,故对应的二阶子系统的系统函数为

978-7-111-42877-0-Chapter04-71.jpg

又设余下的四阶子系统为

H1z)=1+bz-1+cz-2+bz-3+z-4

因为H′z)=H1zH3z)和H1z)、H3z)展开式比较系数间的关系,可得

978-7-111-42877-0-Chapter04-72.jpg

所以有-2cosθ3=-1,解得978-7-111-42877-0-Chapter04-73.jpg。再利用

978-7-111-42877-0-Chapter04-74.jpg

解上述联立方程,得

978-7-111-42877-0-Chapter04-75.jpg

所以七个零点分别为

978-7-111-42877-0-Chapter04-76.jpg

这样,Hz)的零点分布如图4.4.9所示,级联形式的结构流图如图4.4.10所示。

978-7-111-42877-0-Chapter04-77.jpg

4.4.9 例4.4.1的零点分布图

978-7-111-42877-0-Chapter04-78.jpg

图4.4.10 例4.4.1的级联形式结构流图

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈