《数字信号处理课后答案第6章.ppt》由会员分享，可在线阅读，更多相关《数字信号处理课后答案第6章.ppt（96页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、教材第教材第6章习题与上机题解答章习题与上机题解答1 设计一个巴特沃斯低通滤波器，要求通带截止频率fp=6 kHz，通带最大衰减ap=3 dB，阻带截止频率fs=12kHz，阻带最小衰减as=25 dB。求出滤波器归一化系统函数G(p)以及实际滤波器的Ha(s)。解解:（1）求阶数N。将ksp和sp值代入N的计算公式,得所以取N=5（实际应用中，根据具体要求，也可能取N=4，指标稍微差一点，但阶数低一阶，使系统实现电路得到简化）。（2）求归一化系统函数G(p)。由阶数N=5直接查教材第157页表6.2.1,得到五阶巴特沃斯归一化低通滤波器系统函数G(p)为或当然，也可以先按教材（）式计算出极点
2、：再由教材（）式写出G(p)表达式为最后代入pk值并进行分母展开，便可得到与查表相同的结果。（3）去归一化（即LP-LP频率变换），由归一化系统函数G(p)得到实际滤波器系统函数Ha(s)。由于本题中ap=3 dB，即c=p=26103 rad/s，因此对分母因式形式，则有如上结果中，c的值未代入相乘，这样使读者能清楚地看到去归一化后，3 dB截止频率对归一化系统函数的改变作用。2 设计一个切比雪夫低通滤波器，要求通带截止频率fp=3 kHz，通带最大衰减p=0.2 dB，阻带截止频率fs=12 kHz,阻带最小衰减s=50 dB。求出滤波器归一化系统函数G(p)和实际的Ha(s)。解解:（1
3、）确定滤波器技术指标。p=0.2 dB,p=2fp=6103 rad/ss=50 dB,s=2fs=24103 rad/s p=1,(4)求阶数N和。为了满足指标要求，取N=4。(3)求归一化系统函数G(p)其中，极点pk由教材（）式求出如下：(4)将G(p)去归一化，求得实际滤波器系统函数Ha(s)：其中，sk=ppk=6103pk,k=1,2,3,4。因为p4=p1*,p3=p2*,所以，s4=s1*,s3=s2*。将两对共轭极点对应的因子相乘，得到分母为二阶因子的形式，其系数全为实数。也可得到分母多项式形式，请读者自己计算。3 设计一个巴特沃斯高通滤波器，要求其通带截止频率fp=20 k
4、Hz，阻带截止频率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)求相应的归一化低通滤波器技术指标要求：套用图中高通到低通频率转换公式，p=1,s=p/s，得到 p=1,ap=3 dB as=15 dB(3)设计相应的归一化低通G(p)。题目要求采用巴特沃斯类型，故所以，取N=3，查教材中表，得到三阶巴特沃斯归一化低通G(p)为(4)频率变换。将G(p)变换成实际高通滤波器系统函数H(s)：式中c=2fc=220
5、103=4104 rad/s 4.已知模拟滤波器的系统函数Ha(s)如下：(1)(2)式中a、b为常数，设Ha(s)因果稳定，试采用脉冲响应不变法将其转换成数字滤波器H(z)。解解:该题所给Ha(s)正是模拟滤波器二阶基本节的两种典型形式。所以，求解该题具有代表性，解该题的过程，就是导出这两种典型形式的Ha(s)的脉冲响应不变法转换公式。设采样周期为T。(1)Ha(s)的极点为s1=a+jb,s2=ajb将Ha(s)部分分式展开（用待定系数法）：比较分子各项系数可知，A1、A2应满足方程：解之得，A1=1/2，A2=1/2，所以套用教材（）式，得到按照题目要求，上面的H(z)表达式就可作为该题
6、的答案。但在工程实际中，一般用无复数乘法器的二阶基本节结构来实现。由于两个极点共轭对称，所以将H(z)的两项通分并化简整理，可得这样，如果遇到将用脉冲响应不变法转换成数字滤波器时，直接套用上面的公式即可，且对应结构图中无复数乘法器，便于工程实际中实现。(2)Ha(s)的极点为 s1=a+jb,s2=ajb将Ha(s)部分分式展开：套用教材（）式，得到 通分并化简整理，得到5 已知模拟滤波器的系统函数如下：(1)(2)试采用脉冲响应不变法和双线性变换法将其转换为数字滤波器。设T=2 s。解解:.用脉冲响应不变法(1)方法一 直接按脉冲响应不变法设计公式，Ha(s)的极点为将T=2代入上式，得方法
7、二 直接套用4题（2）所得公式。为了套用公式，先对Ha(s)的分母配方，将Ha(s)化成4题中的标准形式：c为一常数由于所以对比可知，套用公式，得(2)或通分合并两项得 用双线性变换法（1）(2)6 设ha(t)表示一模拟滤波器的单位冲激响应，即用脉冲响应不变法，将此模拟滤波器转换成数字滤波器（用h(n)表示单位脉冲响应，即 h(n)=ha(nT)）。确定系统函数H(z)，并把T作为参数，证明：T为任何值时，数字滤波器是稳定的，并说明数字滤波器近似为低通滤波器还是高通滤波器。解解:模拟滤波器系统函数为Ha(s)的极点s1=0.9，故数字滤波器的系统函数应为H(z)的极点为z1=e0.9T,|z
8、1|=e0.9T题6解图所以，T0时，|z1|slsu，所以不满足教材（）式。按照教材（）式,增大sl，则采用修正后的设计巴特沃斯模拟带通滤波器。(3)将带通指标转换成归一化低通指标。套用图中带通到低通频率转换公式，求归一化低通边界频率：(4)设计模拟归一化低通G(p)：取N=3。查教材表，得到归一化低通系统函数G(p)：(5)频率变换，将G(p)转换成模拟带通Ha(s)：(6)用双线性变换公式将Ha(s)转换成H(z)：以上繁杂的设计过程和计算，可以用下面几行程序ex612.m实现。程序运行结果如题12解图所示。得到的系统函数系数为B=0.0234 0 0.0703 0 0.0703 0 0
9、.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)；%调用自编绘图函数myplot绘制带通滤波器的损耗函数曲线题12解图13*设计巴特沃斯数字带通滤
10、波器，要求通带范围为0.25 radw 0.45 rad,通带最大衰减为3 dB，阻带范围为0w0.15 rad和0.55 radw rad，阻带最小衰减为40 dB。调用MATLAB工具箱函数buttord和butter设计，并显示数字滤波器系统函数H(z)的系数，绘制数字滤波器的损耗函数和相频特性曲线。这种设计对应于脉冲响应不变法还是双线性变换法？解解：调用函数buttord和butter设计巴特沃斯数字带通滤波器程序ex613.m如下：%程序ex613.mwp=0.25，0.45；ws=0.15，0.55；rp=3；rs=40；N，wc=buttord(wp，ws，rp，rs)；B，A=
11、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函数buttord和butter是采用双线性变换法来设计巴特沃斯数字滤波器的。数字滤波器的损耗函数和相频特性曲线如题13*解图所示。题13*解图14*.
12、设计一个工作于采样频率80 kHz的巴特沃斯低通数字滤波器，要求通带边界频率为4 kHz，通带最大衰减为0.5 dB，阻带边界频率为20 kHz，阻带最小衰减为45 dB。调用MATLAB工具箱函数buttord和butter设计，并显示数字滤波器系统函数H(z)的系数，绘制损耗函数和相频特性曲线。题14*解图解解：本题以模拟频率给定滤波器指标，所以，程序中先要计算出对应的数字边界频率，然后调用MATLAB工具箱函数buttord和butter来设计数字滤波器。设计程序为ex614.m。%程序ex614.mFs=80000；T=1/Fs；wp=2*pi*4000/Fs；ws=2*pi*2000
13、0/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 2.7188 1.3066 0.2425数字滤波器的损耗函数和相频特性曲线如题14*解图所示。由图可见，滤波器通带截止频率大于0.1（对应的模拟频率分别为4 kHz），阻带截止频率为0.5（对应的模拟频率分别为20 kHz），完全
14、满足设计要求。15*设计一个工作于采样频率80 kHz的切比雪夫型低通数字滤波器，滤波器指标要求与题14*相同。调用MATLAB工具箱函数cheb1ord和cheby1设计，并显示数字滤波器系统函数H(z)的系数，绘制损耗函数和相频特性曲线。与题14*的设计结果比较，简述巴特沃斯滤波器和切比雪夫型滤波器的特点。解解：本题除了调用的MATLAB工具箱函数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=4
15、5;%数字滤波器指标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)的系数：B=0.0023 0.0069 0.0069 0.0023A=1.0000 2.5419 2.2355 0.6753数字滤波器的损耗函数和相频特性曲线如题15*解图所示。由图可见，完全满足设计要求。巴特沃斯滤波器和切比雪夫型滤波器的特点见教材第179页。题15*解图16*设计一个工作于采样频率2500
16、kHz的椭圆高通数字滤波器，要求通带边界频率为325 kHz，通带最大衰减为1 dB，阻带边界频率为225 kHz，阻带最小衰减为40 dB。调用MATLAB工具箱函数ellipord和ellip设计，并显示数字滤波器系统函数H(z)的系数，绘制损耗函数和相频特性曲线。解解：本题求解程序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，hig
17、h)；clf；mpplot(B，A，rs)；%调用本书绘图函数mpplot绘图程序运行结果：阶数N=5，数字滤波器系统函数H(z)的系数：B=0.2784 1.2102 2.2656 2.2656 1.2102 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，阻带最小衰减为
18、50 dB。调用MATLAB工具箱函数ellipord和ellip设计，并显示数字滤波器系统函数H(z)的系数，绘制损耗函数和相频特性曲线。解解：本题求解程序ex617.m如下：%程序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)；%
19、调用本书绘图函数mpplot绘图程序运行结果：阶数N=4，2N阶数字带通滤波器系统函数H(z)的系数：B=0.0043 0.0184 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，阻带边界
20、频率为1050 kHz和1400 Hz，阻带最小衰减为40 dB。调用MATLAB工具箱函数ellipord和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
21、)clf；mpplot(B，A，rs)；%调用本书绘图函数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设计过渡模拟低通滤波器，再调用脉冲响应不变法设计函数im
22、pinvar，将过渡模拟低通滤波器转换成低通数字滤波器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=butter(N，
23、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点的衰减mmin=20*log10(m(end)/max(m)%D
24、F在=点的衰减%以下绘制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*解图