第05章离散时间系统的相位结构2

1、第5章 离散时间系统的相位、 结构与状态变量描述5.1 离散时间系统的相频响应5.2 FIR 系统的线性相位5.3 具有线性相位系统的零点分布5.4 全通系统和最小相位系统5.5 谱分解5.6 FIR 系统的结构5.1 离散时间系统的相频响应()()()jjjH eH ee 幅频响应相频响应() :( ):jH e 如果：( )k 我们称其为线性相位。若：( )k 也称线性相位()()()()()()xjjjjjkjkjY eH eX eX eeX ee()1jH e对输入 ，有( )x n假定：所以：( )()y nx nk输出是输入的简单移位，因此不会发生失真。( )k 例：令100( )

2、cos()cos(2()y nnknk则：波形没有发生失真00( )cos()cos(2)x nnn()jjkH ee具有线性相位例：令若：200( )cos(/4)cos(2)y nnn则：波形发生了失真00( )cos()cos(2)x nnn()()jjH ee 00/403/2( )3/2 -1001020-202-1001020-202-1001020-202( )x n1( )y n2( )y n如果令：0( )cos()x nAn0()1jAHe再令：00000( )cos()cos()y nnn 则：000( )() cos()jy nA H en 则：()()()jjjH e

3、H ee 由于：00000( )cos()cos()y nnn ( )( )p 定义：如果系统的相频响应不是线性的，那么系统的输出将不再是输入信号作线性移位后的组合，因此，输出将发生失真。 ( )( )gdd 定义：为系统的群延迟(Group Delay, GD)为系统的相位延迟(Phase Delay, PD)显然，若系统具有线性相位，则其GD为常数。00( )( )cos(),( ):Narrowband Signalacx nx nnx n若：则：0000( )()()cos()jagpy nH ex nn即：相位延迟 反映了载波信号的延迟， 而群延迟 反映了输出包络的延迟。 0()p0

4、()g思考：如何实现对信号的零相位滤波？若 要保证系统是因果的，又如何实现？5.2 FIR 系统的线性相位 在绝大部分信号处理的场合，人们都期盼系统具有线性相位，但是，如何实现线性相位？对于FIR 系统，如果保证：h(n)是实数，且( )(1)0,1,1h nh NnnN 则该系统具有线性相位。( )(1)h nh Nn :NN偶数奇数上述对称有四种情况：偶对称( )h n( )(1)h nh Nn :NN偶数奇数奇对称( )h n1. 为奇数N10()( )Njj nnH eh n e1(3) 2120(1)/21( )( )2NNNjj nj nnnNNh n eh n ehe1mNn 利

5、用 的对称性，并令 ，则有( )h n(3) 2(3) 2(1)00()( )(1)NNjjnj NmnmH eh n eh Nm e 1212NjNhe( )(1)h nh Nn h(n)偶对称: (3) 2(3) 2(1)00121(3) 2(1)20111(3)/22220()( )(1)121( )()2(NNjj nj NnnnNjNNjj nj NnnNNNNjjnjnnH eh n eh Nn eNheNh n eeheeh nee 1(3)/2201)2112cos22NNjnNhNNeh nnh(3) 2(1)/20112( )cos22Nj NmNNeh mmh()jH e

6、(1) 2nNm令：102( )121,2,(1) 22Nhna nNhnnN 令：(1) 2(1)/21112cos()22Nj NnNNehnnh()jH e实数(1) 2(1)/20()( )cos()Njj NnH eea nn最后有：(1) 20( )(1)2( )( )cos()NgnNHa nn 相位增益 所以，只要保证滤波器的系数偶对称，该滤波器必然具有线性相位。2(1)/211()2cos22Njj NnNH eehnn2. 为偶数N( )2 (),1,2,22NNb nhnn令：2(1)/211()( )cos() 2Njj NnH eeb nn则：( )(1)h nh N

7、n h(n)偶对称:3. 为奇数N1(1) 2221()( )sin()NNjjnH eec nn1( )21,2,(1) 22Nc nhnnN( )(1)/2/2N h(n)奇对称:( )- (1)h nh Nn 4. 为偶数N122211()( )sin2NNjjnH eed nn( )21,2,22Nd nhnnN( )(1)/2/2N 请掌握四种情况下线性相位表达式的推导方法。h(n)奇对称:( )- (1)h nh Nn 四种线性相位四种线性相位FIR滤波器特性滤波器特性四种线性相位四种线性相位FIR滤波器特性滤波器特性的线性组合，在 时， 易取得最大值，因此这一类滤波器易体现低通特

8、性，且是偶函数。通过频率移位，又可体现高通、带通、带阻特性。所以，经典的低通、高通、带通和带阻滤波器的 都是偶对称的。说明： 第1 、 2种情况 FIR 系统是cos()n0()jH e( )h n的线性组合，在 时， 的值为零，且是奇函数。这一类滤波器都是作为特殊形式的滤波器，如 Hilbert变换器、差分器等。第3 、 4种情况 FIR 系统是sin()n0()jH e 最好取为奇数，以便以中心点为对称。N思考：四类滤波器的对称点在何处5.3 具有线性相位系统的零点分布1100( )( )(1)NNnnnnH zh n zh Nn z ( )H z1()H z所以， 的零点也是 的零点,

9、反之亦然1mNn 令：1(1)0( )NNnnzh n z 则：1(1)0( )( )NNmmH zh m z (1)1()NzH z ( )H z 的零点分布: 零点分布可能有四种情况：不在实轴上，也不在单位圆上，应是两对共轭零点不在实轴上，但在单位圆上，也是一对共轭零点在实轴上，但不在单位圆上，无共轭，角度0， 在实轴上，且在单位圆上，无共轭，角度0，( )H z .0, ,1kkar1111( )111111kkkkjjkkkjjkkHzz r ez r ezezerr四个零点同时存在, 构成四阶系统.kz在单位圆内把该式展开，其系数也是对称的，是具有线性相位的子系统。.0, ,1,.k

10、kkbrz在实轴上111( )11mkkHzz rzr无共轭零点, 有镜象零点11( )11kkjjlH zz ez e.0 , ,1,.kkkcrz在单位圆上无镜象对称零点, 有共轭零点.kmlnnlmkzHzHzHzHzH)()()()()(一个具有线性相位的FIR数字滤波器的转移函数可表示为上述四类 FIR 子系统的级联，即：很容易证明，每一个子系统的系数都是对称的，因此它们都具有线性相位。.0, ,11,1kkkdrz)1 ()(1zzHn无镜象零点, 也无共轭零点.5.4 全通系统和最小相位系统 如果一个系统的幅频响应对所有的频率都等于1 (或一个常数), 即1| )(|japeH|

11、0则称系统 为全通系统。)(zHap( )kapHzz最简单的全通系统，纯延迟全通系统1111( )11apzHzz一阶全通系统：Pole:,Zero:1zz镜像对称212( )( )()apapapHzHz Hz22()()()jjjapapapHeHeHe二阶全通系统： 111111(1)(1 ()( )(1)(1)apzzHzzz一对位于单位圆内的共轭极点，一对共轭零点和极点以单位圆为镜像对称 。111( )11NkapkkkzHzz 高阶全通系统：高阶全通系统的另一种表示形式：1(1)111212( )1NNNNapNNaaza zzHza za za z 1()( )( )NapzA

12、 zHzA z 即：2*()()()1jjjapapapHeHeHe对该全通系统： 1 . 是IIR系统(不考虑纯延迟形式）； 2. 极点数和零点数相等； 3. 极点和零点是以单位圆镜像对称的； 4. 极点都在单位圆内,零点都在单位圆外； 5. 全通系统的群延迟始终为正值。全通系统的特点：IIR系统的 无限长，无法对称，即无法作到线性相位。在实际中，可以用一个全通系统和IIR系统相级联，在不改变幅频响应的情况下对相频响应做矫正，使其接近线性相位。( )h n全通系统的应用：全通系统还广泛应用在系统分析及一些特殊滤波器的设计方面（如功率互补IIR滤波器组）-101-1-0.500.51(a)00

13、.20.400.511.5(b)00.20.4-4-3-2-10(c)05101520-0.500.51(d) 一阶全通系统极零图幅频相频抽样响应-101-1-0.500.51(a)00.20.400.51(b)00.20.4-8-6-4-20(c)0102030-0.500.51(d) 三阶全通系统 一个离散系统，其极点必须在单位圆内，但对零点没有限制，如果： 所有的零点都在单位圆内： 最小相位系统；2. 所有的零点都在单位圆外： 最大相位系统；3. 单位圆内、外都有零点 ： 混合相位系统。最小相位系统 在具有相同幅频响应的因果的稳定的滤波器 集合中, 最小相位滤波器具有最小的相位偏移；最小

14、相位系统的性质：例：如下两个系统具有相同的幅频响应：1111( ),1,11zbbzH zabzaaz1211( ),1,11bzbzHzabzaaz那一个是最小相位系统1111( ),1,11zbbzH zabzaaz1211( ),1,11bzbzHzabzaaz幅频相频 在所有具有相同幅频响应的离散系统中, 最 小相位系统的 具有最小的延迟；( )h n令：20()( )0MnE Mh nM 累计能量有：22min00( )( )MMnnhnh n所以，最小相位系统的单位抽样响应又称最小延迟序列。思考： 具有线性相位的FIR系统是否是最小相位系统？/31 2/31 212/31/31/3

15、1/3122/31 2/31 232(1 0.5) (1 0.5)( )1 0.81(1 0.5)(1 0.5)(0.5)(0.5)( )1 0.81(0.5) (0.5)( )1 0.81jjjjjjjjezezH zzezezezezHzzezezHzz例. 三个系统：它们具有相同的幅频响应，试判断，那一个是最小相位系统？最大相位系统？混合相位系统？请注意：为保证系统具有相同的幅频响应（相同的定标）， 的表达式。123( ),( ),( )H zHzHz-101-101222-202-1012-202-10122200.250.50102000.250.5-20-1001001020-20

16、201020-20201020-20201020051015)(arg| )(|)(zHjezHzH( )ln( )ln |( )|arg( )H zH zH zjH z3. 设 为最小相位系统( )H z令：构成一对Hilbert变换)(jReH)(jIeHdeHeHdeHheHjRjIjjR)2cot()(21)()2cot()(21)0()(1则：和( )()( )jH zH eh n、复倒谱：Cepstrum 对于稳定因果系统，当且仅当其是最小相位 系统时, 该系统才有逆系统 （Inverse System)。 令：( )( )( )N zH zD z记：INV1( )( )( )(

17、)D zHzH zN z 的逆系统( )H z( )H zINV( )Hz( )y n( )x n( )x nDeconvolution(反卷积）System identification(系统辨识）5. 任一非最小相位的因果系统的转移函数均可由一个最小相位系统和一个全通系统 级联而成, 即：)()()(minzHzHzHap由于最小相位系统有着以上特殊的性质，因此有着广泛的应用，特别是在信号的建模与系统辨识方面。要理解，具有相同幅频响应的系统，它们所对应的转移函数可以是不相同的，区别就在于相位（或零点的位置）。那么，如何由一个最小相位系统得到具有相同幅频响应的最大相位、混合相位系统？5.5

18、谱分解（Spectral factorization)1( )( )()P zH z H z令：显然， 具有线性相位。将一个转移函数的极零点重新分配，得到两个转移函数，这一过程（或方法）就称为“谱分解”。最常用的是将具有线性相位系统的转移函数作分解，并且往往是分解成两个具有相同幅频响应的子系统。( )P z( )p n=1.0000 ，4.0500，8.1000 ，14.9956，27.7248，43.2996，51.1831，43.2996，27.7248，14.9956，8.1000，4.0500，1.0000例. 令显然，该系统具有线性相位，共有12个零点：0.8,1/0.8,0.6,1

19、/0.6,2/32/3/3/3,0.6,/0.6jjjjeeee 下图是对 作谱分解的结果，可以看出，分解后的两个系统具有相同的幅频响应。( )P z 谱分解的目的是想得到因果的、符合某种要求的系统，这在信号建模、特殊滤波器的设计中经常要用到。分解的一般方法是： 令一个系统是最小相位系统； 则另一个系统必然是最大相位系统。这样，两个系统都有着相同的幅频响应。 另外一种分解方法是得到两个混合系统，目的是保证它们都具有线性相位。0510020406000.250.500.51-2-101-101221200.250.500.510510012-101-101600.250.500.510510051015-2-101-1016( )P z( )P z( )H z1()H z5.6 FIR 系统的结构直接实现:12012( )( )MMY zX z bb zb zb z0( )MnnnH zb z一、 直接实现和级联实现1( )nMNbh n级联实现:LkkkkzzzH122110)()( )(1)h nh Nn :oddM:evenM乘法量减少一半二、 具有线