《数字信号处理课后答案第6章课件.ppt》由会员分享,可在线阅读,更多相关《数字信号处理课后答案第6章课件.ppt(95页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、3ssp3p2 12 10226 10 将ksp和sp值代入N的计算公式, 得lg17.7944.15lg2N 所以取N=5(实际应用中, 根据具体要求, 也可能取N=4, 指标稍微差一点, 但阶数低一阶, 使系统实现电路得到简化)。 (2) 求归一化系统函数G(p)。 由阶数N=5直接查教材第157页表6.2.1, 得到五阶巴特沃斯归一化低通滤波器系统函数G(p)为54321( )3.2365.23615.23613.23611G pppppp或221( )(0.6181)(1.6181)(1)G pppppp当然, 也可以先按教材(6.2.13)式计算出极点: 121j22e 0,1,2,
2、3,4kNkpk再由教材(6.2.12)式写出G(p)表达式为401( )()kkG ppp最后代入pk值并进行分母展开, 便可得到与查表相同的结果。(3) 去归一化(即LP-LP频率变换), 由归一化系统函数G(p)得到实际滤波器系统函数Ha(s)。由于本题中ap=3 dB, 即c=p=26103 rad/s, 因此caa5c54233245ccccc( )( )|3.23615.23615.23613.2361spHsHpsssss对分母因式形式, 则有caa52222ccccc( )( )|(0.6180)(1.6180)()spcHsHpsssss如上结果中,c的值未代入相乘, 这样使
3、读者能清楚地看到去归一化后,3 dB截止频率对归一化系统函数的改变作用。 2 设计一个切比雪夫低通滤波器, 要求通带截止频率fp=3 kHz,通带最大衰减p=0.2 dB,阻带截止频率fs=12 kHz, 阻带最小衰减s=50 dB。 求出滤波器归一化系统函数G(p)和实际的Ha(s)。 解解: (1) 确定滤波器技术指标。 p=0.2 dB, p=2fp=6103 rad/ss=50 dB,s=2fs=24103 rad/s p=1,ssp4(4) 求阶数N和。sp1s0.110.1arch arch 1011456.65101arch 1456.653.8659arch 4aakNkN为了
4、满足指标要求, 取N=4。 p0.11010.2171a (3) 求归一化系统函数G(p)411111( )2()1.7368()NNkkkkQ ppppp其中, 极点pk由教材(6.2.46)式求出如下: (21)(21)ch sinjch cos1,2,3,4221111arsharsh0.558040.2171kkkpkNNN 1234ch0.5580sinj ch0.5580cos0.4438j1.07158833ch0.5580sinj ch0.5580cos1.0715j0.44388855ch0.5580sinj ch0.5580cos1.0715j0.44388877ch0.5
5、580sinj ch0.5580cos0.4438j1.0788pppp 15(4) 将G(p)去归一化, 求得实际滤波器系统函数Ha(s): pa44p44p11( )( )|1.7368()1.7368()sppkkkkHsQ pspss其中, sk=ppk=6103pk, k=1, 2, 3, 4。 因为p4=p1*, p3=p2*, 所以, s4=s1*, s3=s2*。 将两对共轭极点对应的因子相乘, 得到分母为二阶因子的形式, 其系数全为实数。 16222211222482487.2687 10( )(2Re | )(2Re| )7.2687 1016(1.6731 104.779
6、1 10 )(4.0394 104.7790 10 )aHsss ssss ssssss也可得到分母多项式形式, 请读者自己计算。 3 设计一个巴特沃斯高通滤波器, 要求其通带截止频率fp=20 kHz, 阻带截止频率fs=10 kHz, fp处最大衰减为3 dB, 阻带最小衰减as=15 dB。 求出该高通滤波器的系统函数Ha(s)。解解: (1) 确定高通滤波器技术指标要求: p=20 kHz, ap=3 dB fs=10 kHz, as=15 dB (2) 求相应的归一化低通滤波器技术指标要求: 套用图5.1.5中高通到低通频率转换公式, p=1, s=p/s, 得到 p=1, ap=3
7、 dB , 2spsas=15 dB(3) 设计相应的归一化低通G(p)。 题目要求采用巴特沃斯类型, 故ps0.1sp0.1ssppspsp1010.181012lglg0.182.47lglg2aakkN 所以, 取N=3, 查教材中表6.2.1, 得到三阶巴特沃斯归一化低通G(p)为321( )221G pppp(4) 频率变换。 将G(p)变换成实际高通滤波器系统函数H(s): c33223ccc( )( )|22pssH sG psss式中c=2fc=220103=4104 rad/s 4. 已知模拟滤波器的系统函数Ha(s)如下: a22( )()saHssab(1)(2)a22(
8、 )()bHssab式中a、 b为常数, 设Ha(s)因果稳定, 试采用脉冲响应不变法将其转换成数字滤波器H(z)。 解解: 该题所给Ha(s)正是模拟滤波器二阶基本节的两种典型形式。 所以, 求解该题具有代表性, 解该题的过程, 就是导出这两种典型形式的Ha(s)的脉冲响应不变法转换公式。 设采样周期为T。 (1)a22( )()saHssabHa(s)的极点为s1=a+jb, s2=ajb将Ha(s)部分分式展开(用待定系数法): 12a2212122122121 22 122( )()()()()()()AAsaHssabssssA ssA sssabAA sAsA ssab比较分子各项
9、系数可知, A1、 A2应满足方程: 121 22 11AAAsA sa解之得, A1=1/2, A2=1/2, 所以a1/21/2( )(j )(j )Hssabsab 套用教材(6.3.4)式, 得到21(j )1(j )111/21/2( )1 e1 e1 ekkTsab Tab TkAH zzzz 按照题目要求, 上面的H(z)表达式就可作为该题的答案。 但在工程实际中, 一般用无复数乘法器的二阶基本节结构来实现。 由于两个极点共轭对称, 所以将H(z)的两项通分并化简整理, 可得11221ecos()( )1 2ecos()eaTaTaTzbTH zbT zz这样, 如果遇到将a22
10、( )()saHssab用脉冲响应不变法转换成数字滤波器时, 直接套用上面的公式即可, 且对应结构图中无复数乘法器, 便于工程实际中实现。 a22( )()bHssab(2)Ha(s)的极点为 s1=a+jb, s2=ajb将Ha(s)部分分式展开: ajj22( )(j )(j )Hssabsab 套用教材(6.3.4)式, 得到 (j )1(j )1jj22( )1 e1 eab Tab TH zzz 通分并化简整理, 得到1122esin()( )1 2ecos()eaTaTaTzbTH zbT zz5 已知模拟滤波器的系统函数如下: a21( )1Hsss (1)(2)a21( )23
11、1Hsss试采用脉冲响应不变法和双线性变换法将其转换为数字滤波器。 设T=2 s。 解解: . 用脉冲响应不变法(1)a21( )1Hsss 方法一 直接按脉冲响应不变法设计公式, Ha(s)的极点为121313j,j2222ss a1313jj22221133jj33( )1313jj232233jj33( )1 e1 eTTHsssH zzz 将T=2代入上式, 得1 j 311 j 3111112233jj33( )1 e1 e2 3e sin331 2ecos 3eH zzzzzz 方法二 直接套用4题(2)所得公式。 为了套用公式, 先对Ha(s)的分母配方, 将Ha(s)化成4题中
12、的标准形式: a22( )()bHscsabc为一常数由于2222131312422ssss 所以a222312 32( )131322Hssss 对比可知, , 套用公式, 得13,22ab112221111222 3esin()( )31 2ecos()e2 3e sin331 2ecos 3eaTaTaTTzbTH zzbTzzzz(2)a2111( )123112Hsssss21112112111( )|1 e1 e111 e1 eTTTH zzzzz或通分合并两项得12112132(ee )( )1 (ee )ezH zzz 用双线性变换法(1) 11a22111,21111 21
13、2111 21221( )( )|11111(1) (1)(1)(1)(1)12 3zsTTzH zHszzzzzzzzzzzz(2)11a21111111 21 221 21211( )( )|1123111(1)2(1)3(1)(1) 12 62zszH zHszzzzzzzzzzz6 设ha(t)表示一模拟滤波器的单位冲激响应, 即0.9ae 0( )0 0tth tt用脉冲响应不变法, 将此模拟滤波器转换成数字滤波器(用h(n)表示单位脉冲响应, 即 h(n)=ha(nT))。 确定系统函数H(z), 并把T作为参数, 证明: T为任何值时, 数字滤波器是稳定的, 并说明数字滤波器近似
14、为低通滤波器还是高通滤波器。 解解: 模拟滤波器系统函数为0.9a01( )eed0.9tstHstsHa(s)的极点s1=0.9, 故数字滤波器的系统函数应为110.9111( )1 e1 es TTH zzzH(z)的极点为z1=e0.9T, |z1|=e0.9T题6解图所以, T0时, |z1|slsu, 所以不满足教材(6.2.56)式。 按照教材(6.2.57)式, 增大sl, 则plpuslsu0.3538=0.30221.1708 采用修正后的设计巴特沃斯模拟带通滤波器。 (3) 将带通指标转换成归一化低通指标。 套用图5.1.5中带通到低通频率转换公式, sl220slpssl
15、W1, B求归一化低通边界频率: p1,2220slssl0.35380.30221.97440.3022 0.4399WBps3 dB,15 dBaa(4) 设计模拟归一化低通G(p):ps0.10.3sp0.11.5ssppspsp101101 0.18031011011.9744lglg 0.18032.5183lglg1.9744kkNaa 取N=3。 查教材表6.2.1, 得到归一化低通系统函数G(p):321( )221G pppp(5) 频率变换, 将G(p)转换成模拟带通Ha(s): 220wa33w22 322222223300w0ww365432( )( )|()2()2(
16、)0.0850.87981.44840.70760.51240.11010.0443spsBHsG pB ssssBss Bs Bsssssss(6) 用双线性变换公式将Ha(s)转换成H(z): 11a21115123415561234561( )( )|(0.0181 1.7764 100.05434.44090.05432.7756 100.0181)(1 2.2723.51513.26852.31290.96280.278)zsTzH zHszzzzzzzzzzzz以上繁杂的设计过程和计算, 可以用下面几行程序ex612.m实现。 程序运行结果如题12解图所示。 得到的系统函数系数为B
17、 = 0.0234 0 0.0703 0 0.0703 0 0.0234A= 1.0000 2.2100 3.2972 2.9932 2.0758 0.8495 0.2406与手算结果有差别, 这一般是由手算过程中可能产生的计算误差造成的。%程序ex612.m wp=0.25, 0.45; ws=0.15, 0.55; Rp=3; As=15; %设置带通数字滤波器指标参数N, wc=buttord(wp, ws, Rp, As); %计算带通滤波器阶数N和3 dB截止频率WcB, A=butter(N, wc); %计算带通滤波器系统函数分子分母多项式系数向量A和Bmyplot(B, A);
18、 %调用自编绘图函数myplot绘制带通滤波器的损耗函数曲线题12解图13* 设计巴特沃斯数字带通滤波器, 要求通带范围为0.25 rad 0.45 rad, 通带最大衰减为3 dB, 阻带范围为00.15 rad和0.55 rad rad, 阻带最小衰减为40 dB。 调用MATLAB工具箱函数buttord和butter设计, 并显示数字滤波器系统函数H(z)的系数, 绘制数字滤波器的损耗函数和相频特性曲线。 这种设计对应于脉冲响应不变法还是双线性变换法?解解: 调用函数buttord和butter设计巴特沃斯数字带通滤波器程序ex613.m如下: %程序ex613.mwp=0.25, 0
19、.45; ws=0.15, 0.55; rp=3; rs=40; N, wc=buttord(wp, ws, rp, rs); B, A=butter(N, wc)clf; mpplot(B, A, rs)程序运行结果: 数字滤波器系统函数H(z)的系数: B=0.000100.000700.002200.00360 0.003500.002200.000700.0001A =1.00005.3093 16.291334.729756.939974.512280.0136 71.117071.117052.640832.227016.16966.4618 1.98310.4218 0.0524
20、函数buttord和butter是采用双线性变换法来设计巴特沃斯数字滤波器的。 数字滤波器的损耗函数和相频特性曲线如题13*解图所示。题13*解图14*. 设计一个工作于采样频率80 kHz的巴特沃斯低通数字滤波器, 要求通带边界频率为4 kHz, 通带最大衰减为0.5 dB, 阻带边界频率为20 kHz, 阻带最小衰减为45 dB。 调用MATLAB工具箱函数buttord和butter设计, 并显示数字滤波器系统函数H(z)的系数, 绘制损耗函数和相频特性曲线。 题14*解图解解: 本题以模拟频率给定滤波器指标, 所以, 程序中先要计算出对应的数字边界频率, 然后调用MATLAB工具箱函数
21、buttord和butter来设计数字滤波器。 设计程序为ex614.m。% 程序ex614.mFs=80000; T=1/Fs; wp=2*pi*4000/Fs; ws=2*pi*20000/Fs; rp=0.5; rs=45; N, wc=buttord(wp/pi, ws/pi, rp, rs); B, A=butter(N, wc); clf; mpplot(B, A, rs); %调用本书绘图函数mpplot绘图程序运行结果: 阶数N=4, 数字滤波器系统函数H(z)的系数: B= 0.0028 0.0111 0.0166 0.0111 0.0028A= 1.0000 2.6103
22、2.7188 1.3066 0.2425数字滤波器的损耗函数和相频特性曲线如题14*解图所示。 由图可见, 滤波器通带截止频率大于0.1(对应的模拟频率分别为4 kHz), 阻带截止频率为0.5(对应的模拟频率分别为20 kHz), 完全满足设计要求。 15* 设计一个工作于采样频率80 kHz的切比雪夫型低通数字滤波器, 滤波器指标要求与题14*相同。 调用MATLAB工具箱函数cheb1ord和cheby1设计, 并显示数字滤波器系统函数H(z)的系数, 绘制损耗函数和相频特性曲线。 与题14*的设计结果比较, 简述巴特沃斯滤波器和切比雪夫型滤波器的特点。 解解: 本题除了调用的MATLA
23、B工具箱函数cheb1ord和cheby1与题14*不同以外, 程序与14*题完全相同。 本题求解程序ex615.m如下: % 程序ex615.mFs=80000; T=1/Fs; wp=2*pi*4000/Fs; ws=2*pi*20000/Fs; rp=0.5; rs=45; %数字滤波器指标N, wp=cheb1ord(wp/pi, ws/pi, rp, rs); B, A=cheby1(N, rp, wp); clf; mpplot(B, A,rs); %调用本书绘图函数mpplot绘图程序运行结果: 阶数N=3, 比题14*设计的巴特沃斯滤波器低1阶。 数字滤波器系统函数H(z)的系
24、数: B=0.0023 0.0069 0.0069 0.0023A=1.0000 2.5419 2.2355 0.6753数字滤波器的损耗函数和相频特性曲线如题15*解图所示。 由图可见, 完全满足设计要求。 巴特沃斯滤波器和切比雪夫型滤波器的特点见教材第179页。题15*解图16* 设计一个工作于采样频率2500 kHz的椭圆高通数字滤波器, 要求通带边界频率为325 kHz, 通带最大衰减为1 dB, 阻带边界频率为225 kHz, 阻带最小衰减为40 dB。 调用MATLAB工具箱函数ellipord和ellip设计, 并显示数字滤波器系统函数H(z)的系数, 绘制损耗函数和相频特性曲线
25、。 解解: 本题求解程序ex616.m如下: % 程序ex616.mFs=2500000; fp=325000; rp=1; fs=225000; rs=40; %滤波器指标wp=2*fp/Fs; ws=2*fs/Fs; %将边界频率转换为数字频率N, wpo=ellipord(wp, ws, rp, rs); B, A=ellip(N, rp, rs, wpo, high); clf; mpplot(B, A, rs); %调用本书绘图函数mpplot绘图程序运行结果: 阶数N=5, 数字滤波器系统函数H(z)的系数: B=0.2784 1.2102 2.2656 2.2656 1.2102
26、 0.2784A=1.0000 2.1041 2.5264 1.4351 0.4757 0.0329数字滤波器的损耗函数和相频特性曲线如题16*解图所示。 由图可见, 完全满足设计要求。题16*解图17*. 设计一个工作于采样频率5 MHz的椭圆带通数字滤波器, 要求通带边界频率为560 kHz和780 kHz, 通带最大衰减为0.5 dB, 阻带边界频率为375 kHz和1 MHz, 阻带最小衰减为50 dB。 调用MATLAB工具箱函数ellipord和ellip设计, 并显示数字滤波器系统函数H(z)的系数, 绘制损耗函数和相频特性曲线。 解解: 本题求解程序ex617.m如下: % 程
27、序ex617.mfpl=560000; fpu=780000; fsl=375000; fsu=1000000; Fs=5000000; %滤波器指标wp=2*fpl/Fs,2*fpu/Fs;ws=2*fsl/Fs, 2*fsu/Fs; rp=0.5; rs=50; %将边界频率转换为数字频率N, wpo=ellipord(wp, ws, rp, rs); B, A=ellip(N, rp, rs, wpo); clf; mpplot(B, A, rs); %调用本书绘图函数mpplot绘图程序运行结果: 阶数N=4, 2N阶数字带通滤波器系统函数H(z)的系数: B =0.0043 0.01
28、84 0.0415 0.0638 0.0734 0.0638 0.0415 0.01840.0043A =1.0000 5.1091 13.4242 22.3290 25.6190 20.5716 11.39363.9943 0.7205数字滤波器的损耗函数和相频特性曲线如题17*解图所示。 由图可见, 完全满足设计要求。题17*解图18*. 设计一个工作于采样频率5 kHz的椭圆带阻数字滤波器, 要求通带边界频率为500 Hz和2125 Hz, 通带最大衰减为1 dB, 阻带边界频率为1050 kHz和1400 Hz, 阻带最小衰减为40 dB。 调用MATLAB工具箱函数ellipord和
29、ellip设计, 并显示数字滤波器系统函数H(z)的系数, 绘制损耗函数和相频特性曲线。 题18*解图解解: 本题求解程序ex618.m如下: % 程序ex618.mfpl=500; fpu=2125; fsl=1050; fsu=1400; Fs=5000; rp=1; rs=40; %滤波器指标wp=2*fpl/Fs, 2*fpu/Fs; ws=2*fsl/Fs, 2*fsu/Fs; %将边界频率转换为数字频率N, wpo=ellipord(wp, ws, rp, rs)B, A=ellip(N, rp, rs, wpo, stop)clf; mpplot(B, A,rs);%调用本书绘图
30、函数mpplot绘图程序运行结果: 阶数N=3, 2N阶数字带阻滤波器系统函数H(z)的系数: B =0.0748 0.0557 0.1618 0.0897 0.1618 0.0557 0.0748A =1.0000 0.2604 1.2316 0.0633 1.0458 0.0040 0.3412数字滤波器的损耗函数和相频特性曲线如题18*解图所示。 19*. 用脉冲响应不变法设计一个巴特沃斯低通数字滤波器, 指标要求与题14*的相同。编写程序先调用MATLAB工具箱函数buttord和butter设计过渡模拟低通滤波器, 再调用脉冲响应不变法设计函数impinvar, 将过渡模拟低通滤波器
31、转换成低通数字滤波器H(z), 并显示过渡模拟低通滤波器和数字滤波器系统函数的系数, 绘制损耗函数和相频特性曲线。 请归纳本题的设计步骤和所用的计算公式, 并比较本题与题14*的设计结果, 观察双线性变换法的频率非线性失真和脉冲响应不变法的频谱混叠失真。 解解: 本题求解程序ex619.m如下: % 程序ex619.mFs=80000; T=1/Fs; fp=4000; fs=20000; rp=0.5; rs=45; %相应的模拟滤波器指标wp=2*pi*fp; ws=2*pi*fs; %将边界频率转换为角频率N, wc=buttord(wp, ws, rp, rs, s); B, A=bu
32、tter(N, wc, s); Bz, Az=impinvar(B, A, Fs) %调用转换函数impinvar将AF转换成DF%以下计算AF和DF的频响特性fk=0: 10: Fs/2; omega=2*pi*fk; %对AF频响函数在0, Fs/2上以间隔10 Hz采样Hs=freqs(B, A, omega); ms=abs(Hs); ps=angle(Hs); H, W=freqz(Bz, Az, 1000); %对DF频响函数在0, Fs/2采样1000点m=abs(H); p=angle(H); msmin=20*log10(ms(end)/max(ms) %AF在f=Fs/2点
33、的衰减mmin=20*log10(m(end)/max(m) %DF在 = 点的衰减%以下绘制AF和DF的损耗函数和相频特性曲线(省略)程序运行结果: 阶数N=4, N阶数字低通滤波器系统函数H(z)的系数: Bz =0.0000 0.0043 0.0128 0.0024 0Az =1.0000 2.8902 3.2452 1.6605 0.3250模拟滤波器的损耗函数和相频特性曲线如题19*解图(a)和(b)所示, 数字滤波器的损耗函数和相频特性曲线如题19*解图 (c)和(d)所示。 由图可见, 脉冲响应不变法设计的数字滤波器的频响特性基本模拟了模拟滤波器的频响形状, 但存在频谱混叠失真。 模拟滤波器的损耗函数在f=Fs/2点的衰减为msmin =69.0823 dB, 而数字滤波器的损耗函数在 = 点的衰减为mmin =63.4990 dB, 这就是频谱混叠失真引起了-5.5832 dB的衰减误差。 题14*是用双线性变换法设计的, 不存在频谱混叠失真, 但存在频率非线性失真, 所以数字滤波器的频响曲线形状与模拟滤波器的频响形状差别较大, 而且, 频率越高, 频率非线性失真越严重。 本题的设计步骤和所用的计算公式请读者在教材6.3节查找。 题19*解图
限制150内