第4章 离散时间系统结构及滤波器设计4(FIR滤波器设计2)

1、n1、四大类FIR线性相位滤波器长度，对称特性，群延迟及特殊零点n2、求n3、求h(n)1jw)e(Hjdxxx)sin()2/Mn()2/Mn(sin)2/Mn()2/Mn(sinnhcccclp2n具有广义线性相位的理想高通滤波器的频率响应为：n其单位冲激响应为：高通滤波器cM/2j -cjhpe0)e (H)2/Mn()2/Mn(sin)2/Mn()2/Mn(sinnhchp34n所要求设计的高通滤波器的技术指标：n根据凯泽窗设计公式可得到：n由窗函数设计法的对称性可得：425. 0)/2(24M2.6,021. 0,5 . 0,35. 0,1)e (H1,)e (Hpsc21psp1j

2、hp1s2jhp5第I类线性相位高通滤波器67n通过验证计算可得此时的通带逼近误差为0.0213，高于预定指标0.021，n可以通过增大的方法而达到预定指标；n也可以保持不变,改变Mn将M增大至25，调整过渡带变窄而满足指标，但是此时为II类线性相位系统，在=处为零，不适合高通系统，n因此可以将M增大至26就会解决此问题。89n多通带滤波器1G2G3G4w1w2w3G1wcmbN1kk1kkmb2/M-jM/2)-(nM/2)-(nsin)GG()n(he)(逼近误差与跳变幅度成正比，在设计时需要进行指标的归一化，利用最苛刻的指标进行凯泽窗设计带通,带阻滤波器1011离散时间微分器n理想的连续

3、时间微分器的频率响应为：n其理想冲激响应为：2diffM/2)-(nM/2)-(nsinM/2)-(nM/2)-(nc)n(hos2/M-je )()(jwHdiff12n零相微分器000!1)1()sin(21)sin(21)sin(21)cos(221)n(h2dnnndxxxndwwnwdwwnwdwwnwjdwjwennnjwnjwHdiff)(13M=5,b=2.4M=10,b=2.414n频率抽样法从频域出发，把给定的理想频率抽样法从频域出发，把给定的理想频率响应频率响应Hd(ejw)加以等间隔抽样，即：加以等间隔抽样，即：n以此作为实际以此作为实际FIR数字滤波器的频率特数字滤波

4、器的频率特性的抽样值性的抽样值H(k)，即令：，即令： kHeHdkNwjwd2)( 1, 1 , 0,)()(2 NkeHkHkHkNwjwdd三、三、 频率抽样法频率抽样法15有限长序列的频域重构1,.2 ,1 ,0)(|)()(22NkenhzHkHnknNjezkNj)()(1)(102nRekHNnhNNkknNj101210121)(111)(1NkkNjNNkkNjNzekHNzzezkHN1010211010210)(1)(1)()(NkNnnkNjNnnNkknNjNnnezkHNzekHNznhzH161011)(1)(NkkNNzWkHNzzH kNwkHeHNkjw21

5、0212sin2sin1)(NjwewwNNw1021)2()1(211021102110102210)2(22sin2)2(sin1)(2sin2sin1)(2sin2sin1)(2sin2sin1)(2sin2sin1)(1)(1NkNkNwjNkNjNjwNkNkjNjwkjNkNkjNjwNkNkNkwjwNjNkNkwjjwNjwekNwNkNwNkHeeNkwkwNNkHeeNkwekwNNkHeeNkwwNNkHNkwwNeeNkHekHNeeH10)(12)(10)(12)(1)()1 (11)()1 (1)(NkzHkNjzHcNNNkzHkNjrzHcNNKKzrekHzr

6、NzrekHzrNzH 频率抽样法需要注意的问题n1、实现结构稳定性的考虑实现结构稳定性的考虑n极点在单位圆上，存在极点在单位圆上，存在系数量化误差系数量化误差时，时，极点会极点会移动移动,若若不能被梳状滤波器的零点所抵消，系统不能被梳状滤波器的零点所抵消，系统就就不能保持稳定不能保持稳定n为克服此缺点，可以将为克服此缺点，可以将所有零、极点都移到单位所有零、极点都移到单位圆内圆内某一靠近单位圆、半径为某一靠近单位圆、半径为r(小于或近似等于小于或近似等于1)的圆上（的圆上（r为正实数）。为正实数）。18n频率抽样点上滤波器实际频率响应严格地和理想频率响应数值相频率抽样点上滤波器实际频率响应严

7、格地和理想频率响应数值相等。等。n抽样点之间的频率响应则是由各抽样点的加权内插函数的延伸叠抽样点之间的频率响应则是由各抽样点的加权内插函数的延伸叠加而形成的，有一定的逼近误差。加而形成的，有一定的逼近误差。n逼近误差取决于理想频率响应曲线形状，采样点数等逼近误差取决于理想频率响应曲线形状，采样点数等.2、过渡带抽样优化设计、过渡带抽样优化设计19n一般过渡带取二、三点抽样值即可得到满意结果。一般过渡带取二、三点抽样值即可得到满意结果。n在低通设计中，不加过渡抽样点时，阻带最小衰在低通设计中，不加过渡抽样点时，阻带最小衰减为减为-2OdB，n一点过渡抽样的一点过渡抽样的最优最优设计，阻带最小衰减

8、可提高到设计，阻带最小衰减可提高到-4O到到-54dB左右，左右，n二点二点过渡抽样的最优设计可达过渡抽样的最优设计可达-6OdB到到-75dB左右，左右，n三点过渡抽样的最优设计则可达三点过渡抽样的最优设计则可达-80dB到到-95dB左右。左右。n频率抽样法的频率抽样法的优点优点是可以在是可以在频域直接设计频域直接设计，并且，并且适合于最优化设计适合于最优化设计 20n如果设计的是线性相位的如果设计的是线性相位的FIR数字滤波数字滤波器，则其抽样值的幅度和相位一定要满器，则其抽样值的幅度和相位一定要满足对称的约束条件。足对称的约束条件。n第第I类线性相位类线性相位FIR滤波器：滤波器： n

9、幅度函数应为偶对称的：幅度函数应为偶对称的：n如果抽样值也用幅值与相角表示，则：如果抽样值也用幅值与相角表示，则：n必须满足：必须满足：wNjjwewHeH)21()()()2()(wHwHkkjkjkNjeHekNHeHkH)2()()(2NkkNNk11221kNkHH3、线性相位的约束、线性相位的约束21n对对 进行频率抽样，就是在进行频率抽样，就是在z平面单位平面单位圆上的圆上的N个等间隔点上抽取出频率响应值。个等间隔点上抽取出频率响应值。n为保证线性相位特性，在单位圆上可以有为保证线性相位特性，在单位圆上可以有两种抽样方式：两种抽样方式：n第一种是第一个抽样点在第一种是第一个抽样点在

10、 处处(或或 )n第二种是第一个抽样点在第二种是第一个抽样点在 处处(或或 )n每种方式可分为每种方式可分为N是偶数、奇数两种是偶数、奇数两种)(jwdeH0w10jezNw/Njez4、频率抽样的两种方法、频率抽样的两种方法2223n第一种频率抽样第一种频率抽样(I型型)为：为：n其内插公式保持不变。其内插公式保持不变。10 ,)()()(2NkeHkHkHkNwjwdd1011)(1)(NkkNNzWkHNzzH24n第二种频率抽样第二种频率抽样(型型)为：为：nH(z)? 10 ,)()()()()(21010)21(2)21(22NkeenhenheHkHkHnkNjNnnNjNnnk

11、NjkNkNNwjwdd1,.2 , 1 , 0)(| )()(22NkenxzXkXnknNjezkNj)()(1)(102nRekXNnxNNkknNjnkNjNnnNjNnnkNjkNkNNwjweenxenxeXkX21010)21(2)21(22)()()()(102)(1NknkNjekXN10)21(2)(1)(NknkNjekHNnhnNjenx)(26n求求z变换，可得：变换，可得：101)21(21010)21(21)(1)(1)(NkkNjNNnnNknkNjzekHNzzekHNzHNje)0(HNje3) 1(NHNNje)12() 1 (HNzN1)(nx)(ny1

12、z1z1z27低通滤波器设计为例：设采样点数为N,截止频率为c , c 对应的样本点序号为kc ，N为偶数时 (1)/(1)/ 0,1,.,( )0 +1, 1 ,.,1j NkNcdccj NkNceekkHkkkNkkNkNN为奇数时 (1)/ 0,1,.,; ,1( )0 j NkNccdekkkNkNHkk为其他值频率抽样法设计任意滤波器n滤波器幅度取样值要求为n1）可以直接由频率内插公式求得频率响应，并由频率内插结构实现滤波器n2）也可以由H(K)求得h(n)，从而求取频率响应，并实现滤波器28elsekkkkH13,315,14,2 ,1 ,005.01)(频率抽样法设计低通滤波器

13、29按第一种频率抽样设计，N=303010Int42( )0Int151522ccNkH kNNk 515130sin 30sin 30sin230230230sin30sin30sin2230230jjkkkH eekk3421.6 10220.21321.5 10pppsswf3423.1 10220.41321.5 10stststsswf250dB过渡点优化n为改善频率特性，以满足指标要求，可在通带和阻带交界处安排一个或者几个不等于1的抽样值。n此处令：n优化后最小阻带衰减约为60dB左右31(6)0.1065H(5)0.5886H324.5 IIR与与FIR滤波器特性比较滤波器特性比

14、较n1、在相同的技术指标下，、在相同的技术指标下，IIR滤波器可用比滤波器可用比FIR滤波器较少的阶数来满足指标的要求，所滤波器较少的阶数来满足指标的要求，所用的存储单元少，运算次数少，较为经济。用的存储单元少，运算次数少，较为经济。n例如，用频率抽样法设计阻带衰减为例如，用频率抽样法设计阻带衰减为-20dB的的FIR滤波器，其阶数要滤波器，其阶数要33阶阶才能达到，而才能达到，而双线性变换法双线性变换法设计只需设计只需45阶阶的切贝雪夫的切贝雪夫IIR滤波器即可达到指滤波器即可达到指标要求，所以阶数要高标要求，所以阶数要高510倍左右。倍左右。n2、FIR滤波器可得到严格的线性滤波器可得到严

15、格的线性相位相位，IIR滤滤波器的频率选择性愈好，其相位的波器的频率选择性愈好，其相位的非线性失真非线性失真愈严重愈严重。33n 3、FIR滤波器主要采用非递归结构，因而无滤波器主要采用非递归结构，因而无论是从理论上还是从实际的有限精度的运算中论是从理论上还是从实际的有限精度的运算中它都是稳定的，有限精度运算的误差也较小。它都是稳定的，有限精度运算的误差也较小。nIIR滤波器必须采用递归结构，极点必须在滤波器必须采用递归结构，极点必须在z平面单平面单位圆内才能稳定，对于这种结构，有限精度处理有位圆内才能稳定，对于这种结构，有限精度处理有时会引起寄生振荡（零极限环）。时会引起寄生振荡（零极限环）

16、。n 4、对于、对于FIR滤波器，由于冲激响应是有限长滤波器，由于冲激响应是有限长的，因而可以用快速傅里叶变换算法，运算速的，因而可以用快速傅里叶变换算法，运算速度可以快得多。度可以快得多。IIR滤波器则不能这样运算。滤波器则不能这样运算。34n5、IIR滤波器可以利用模拟滤波器设计的现滤波器可以利用模拟滤波器设计的现成公式、数据和表格，因而计算工作量较小，成公式、数据和表格，因而计算工作量较小，对计算工具要求不高。对计算工具要求不高。nFIR窗函数法只给出窗函数的计算公式，但计算通窗函数法只给出窗函数的计算公式，但计算通带、阻带衰减仍无显式表达。带、阻带衰减仍无显式表达。nFIR滤波器优化设计，需要借助于计算机滤波器优化设计，需要借助于计算机 。n6、IIR滤波器主要是设计规