数字信号处理习题答案西安电子第7章课件.ppt

《数字信号处理习题答案西安电子第7章课件.ppt》由会员分享，可在线阅读，更多相关《数字信号处理习题答案西安电子第7章课件.ppt（103页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、有限脉冲响应(FIR)数字滤波器的设计第章解：（1）由所给h(n)的取值可知，h(n)满足h(n)=h(N 1n)，所以FIR 滤波器具有A 类线性相位特性：由于N=6 为偶数(情况2)，所以幅度特性关于=点奇对称。（2）由题中h(n)值可知，h(n)满足h(n)=h(N 1n)，所以FIR 滤波器具有B 类线性相位特性：由于7为奇数(情况3)，所以幅度特性关于=0,2 三点奇对称。有限脉冲响应(FIR)数字滤波器的设计第章2 已知第一类线性相位FIR 滤波器的单位脉冲响应长度为16，其16个频域幅度采样值中的前9个为：Hg(0)=12，Hg(1)=8.34，Hg(2)=3.79，Hg(3)H

2、g(8)=0 根据第一类线性相位FIR 滤波器幅度特性Hg()的特点，求其余7个频域幅度采样值。解：因为N=16 是偶数（情况2），所以FIR 滤波器幅度特性Hg()关于=点奇对称，即Hg(2)=Hg()。其N 点采样关于k=N/2 点奇对称，即Hg(N k)=Hg(k)k=1,2,15综上所述,可知其余7个频域幅度采样值：Hg(15)=Hg(1)=8.34，Hg(14)=Hg(2)=3.79，Hg(13)Hg(9)=0有限脉冲响应(FIR)数字滤波器的设计第章3 设FIR 滤波器的系统函数为求出该滤波器的单位脉冲响应h(n)，判断是否具有线性相位，求出其幅度特性函数和相位特性函数。解:对FI

3、R 数字滤波器，其系统函数为有限脉冲响应(FIR)数字滤波器的设计第章由h(n)的取值可知h(n)满足：h(n)=h(N 1n)N=5所以，该FIR 滤波器具有第一类线性相位特性。频率响应函数H(ej)为所以其单位脉冲响应为有限脉冲响应(FIR)数字滤波器的设计第章有限脉冲响应(FIR)数字滤波器的设计第章幅度特性函数为 相位特性函数为4 用矩形窗设计线性相位低通FIR 滤波器，要求过渡带宽度不超过/8 rad。希望逼近的理想低通滤波器频率响应函数Hd(ej)为有限脉冲响应(FIR)数字滤波器的设计第章（1）求出理想低通滤波器的单位脉冲响应hd(n)；（2）求出加矩形窗设计的低通FIR 滤波器

4、的单位脉冲响应h(n)表达式，确定 与N 之间的关系；（3）简述N 取奇数或偶数对滤波特性的影响。解:（1）有限脉冲响应(FIR)数字滤波器的设计第章（2）为了满足线性相位条件，要求，N 为矩形窗函数长度。因为要求过渡带宽度 rad,所以要求,求解得到N32。加矩形窗函数,得到h(n)：有限脉冲响应(FIR)数字滤波器的设计第章（3）N 取奇数时，幅度特性函数Hg()关于=0,2 三点偶对称，可实现各类幅频特性；N 取偶数时，Hg()关于=奇对称，即Hg()=0，所以不能实现高通、带阻和点阻滤波特性。5 用矩形窗设计一线性相位高通滤波器，要求过渡带宽度不超过/10 rad。希望逼近的理想高通滤

5、波器频率响应函数Hd(ej)为有限脉冲响应(FIR)数字滤波器的设计第章（1）求出该理想高通的单位脉冲响应hd(n)；（2）求出加矩形窗设计的高通FIR 滤波器的单位脉冲响应h(n)表达式，确定 与N 的关系；（3）N 的取值有什么限制？为什么？解:（1）直接用IFTHd(ej)计算：有限脉冲响应(FIR)数字滤波器的设计第章有限脉冲响应(FIR)数字滤波器的设计第章hd(n)表达式中第2项正好是截止频率为c的理想低通滤波器的单位脉冲响应。而(n)对应于一个线性相位全通滤波器：Hdap(ej)=ej即高通滤波器可由全通滤波器减去低通滤波器实现。（2）用N 表示h(n)的长度，则h(n)=hd(

6、n)RN(n)=有限脉冲响应(FIR)数字滤波器的设计第章为了满足线性相位条件：h(n)=h(N 1n)要求满足（3）N 必须取奇数。因为N 为偶数时（情况2），H(ej)=0，不能实现高通。根据题中对过渡带宽度的要求，N 应满足：,即N40。取N=41。6 理想带通特性为 有限脉冲响应(FIR)数字滤波器的设计第章（1）求出该理想带通的单位脉冲响应hd(n)；（2）写出用升余弦窗设计的滤波器的h(n)表达式，确定N 与 之间的关系；（3）要求过渡带宽度不超过/16 rad。N 的取值是否有限制？为什么？解:（1）有限脉冲响应(FIR)数字滤波器的设计第章上式第一项和第二项分别为截止频率c+B

7、 和c的理想低通滤波器的单位脉冲响应。所以，上面hd(n)的表达式说明，带通滤波器可由两个低通滤波器相减实现。（2）h(n)=hd(n)w(n)为了满足线性相位条件，与N 应满足有限脉冲响应(FIR)数字滤波器的设计第章实质上，即使不要求具有线性相位，与N 也应满足该关系，只有这样，才能截取hd(n)的主要能量部分，使引起的逼近误差最小。（3）N 取奇数和偶数时，均可实现带通滤波器。但升余弦窗设计的滤波器过渡带为8/N，所以，要求，即要求N128。7 试完成下面两题：（1）设低通滤波器的单位脉冲响应与频率响应函数分别为h(n)和H(ej)，另一个滤波器的单位脉冲响应为h1(n)，它与h(n)的

8、关系是h1(n)=(1)nh(n)。试证明滤波器h1(n)是一个高通滤波器。有限脉冲响应(FIR)数字滤波器的设计第章（2）设低通滤波器的单位脉冲响应与频率响应函数分别为h(n)和H(ej)，截止频率为c,另一个滤波器的单位脉冲响应为h2(n)，它与h(n)的关系是h2(n)=2h(n)cos0n，且c0(c)。试证明滤波器h2(n)是一个带通滤波器。解:（1）由题意可知对h1(n)进行傅里叶变换,得到有限脉冲响应(FIR)数字滤波器的设计第章上式说明H1(ej)就是H(ej)平移 的结果。由于H(ej)为低通滤波器，通带位于以=0 为中心的附近邻域，因而H1(ej)的通带位于以=为中心的附近

9、，即h1(n)是一个高通滤波器。有限脉冲响应(FIR)数字滤波器的设计第章这一证明结论又为我们提供了一种设计高通滤波器的方法（设高通滤波器通带为 c,）：设计一个截止频率为c的低通滤波器hLp(n)。对hLp(n)乘以cos(n)即可得到高通滤波器hHp(n)cos(n)=(1)nhLp(n)。(2)与(1)同样道理，代入h2(n)=2h(n)cos0n,可得有限脉冲响应(FIR)数字滤波器的设计第章因为低通滤波器H(ej)通带中心位于=2k，且H2(ej)为H(ej)左右平移0，所以H2(ej)的通带中心位于=2k0处，所以h2(n)具有带通特性。这一结论又为我们提供了一种设计带通滤波器的方

10、法。8 题8图中h1(n)和h2(n)是偶对称序列，N=8，设 H1(k)=DFT h1(n)k=0，1，N 1 H2(k)=DFT h2(n)k=0，1，N 1（1）试确定H1(k)与 H2(k)的具体关系式。|H1(k)|=|H2(k)|是否成立？为什么？（2）用h1(n)和h2(n)分别构成的低通滤波器是否具有线性相位？群延时为多少？有限脉冲响应(FIR)数字滤波器的设计第章题8图有限脉冲响应(FIR)数字滤波器的设计第章解：（1）由题8图可以看出h2(n)与h1(n)是循环移位关系：h2(n)=h1(n+4)8R8(n)由DFT 的循环移位性质可得(2)由题8图可知，h1(n)和h2(

11、n)均满足线性相位条件：h1(n)=h1(N 1n)h2(n)=h2(N 1n)有限脉冲响应(FIR)数字滤波器的设计第章所以，用h1(n)和h2(n)构成的低通滤波器具有线性相位。直接计算FT h1(n)和h2(n)也可以得到同样的结论。设 所以，群延时为有限脉冲响应(FIR)数字滤波器的设计第章9 对下面的每一种滤波器指标，选择满足FIRDF 设计要求的窗函数类型和长度。（1）阻带衰减为20 dB，过渡带宽度为1 kHz，采样频率为12 kHz；（2）阻带衰减为50 dB，过渡带宽度为2 kHz，采样频率为20 kHz；（3）阻带衰减为50 dB，过渡带宽度为500 Hz，采样频率为5 k

12、Hz。解：我们知道，根据阻带最小衰减选择窗函数类型，根据过渡带宽度计算窗函数长度。为了观察方便，重写出教材第211 页中表7.2.2。有限脉冲响应(FIR)数字滤波器的设计第章有限脉冲响应(FIR)数字滤波器的设计第章结合本题要求和教材表7.2.2，选择结果如下：（1）矩形窗满足本题要求。过渡带宽度1 kHz 对应的数字频率为B=200/12 000=/60，精确过渡带满足：1.8/N/60，所以要求N1.860=108。（2）选 哈 明 窗，过 渡 带 宽 度1 kHz 对 应 的 数 字 频 率 为B=4000/20 000=/5，精 确 过 渡 带 满 足：6.6/N/5，所 以要求N6

13、.65=33。（3）选哈明窗，过渡带宽度1 kHz 对应的数字频率为B=1000/5000=/5，精确过渡带满足：6.6/N/5，所以要求N6.65=33。10 利用矩形窗、升余弦窗、改进升余弦窗和布莱克曼窗设计线性相位FIR 低通滤波器。要求希望逼近的理想低通滤波器通带截止频率c=/4 rad，N=21。求出分别对应的单位脉冲响应。有限脉冲响应(FIR)数字滤波器的设计第章解:(1)希望逼近的理想低通滤波器频响函数Hd(ej)为 其中,a=(N 1)/2=10。(2)由Hd(ej)求得hd(n)：有限脉冲响应(FIR)数字滤波器的设计第章(3)加窗得到FIR 滤波器单位脉冲响应h(n)：升余

14、弦窗:有限脉冲响应(FIR)数字滤波器的设计第章 改进升余弦窗：布莱克曼窗:有限脉冲响应(FIR)数字滤波器的设计第章11 将技术要求改为设计线性相位高通滤波器，重复题10。解:方法一 将题10解答中的逼近理想低通滤波器(Hd(ej)、hd(n)改为如下理想高通滤波器即可。有限脉冲响应(FIR)数字滤波器的设计第章有限脉冲响应(FIR)数字滤波器的设计第章上式中(n10)对应于全通滤波器。上式说明，高通滤波器的单位脉冲响应等于全通滤波器的单位脉冲响应减去低通滤波器的单位脉冲响应。仿照10题，用矩形窗、升余弦窗、改进升余弦窗和布菜克曼窗对上面所求的hd(n)加窗即可。计算与绘图程序与题10解中类

15、同，只要将其中的h(n)用本题的高通h(n)替换即可。方法二 根据第7题（1）的证明结论设计。（1）先设计通带截止频率为/4 的低通滤波器。对四种窗函数所得FIR 低通滤波器单位脉冲响应为题9解中的hR(n)、hHn(n)、hHm(n)和hBl(n)。有限脉冲响应(FIR)数字滤波器的设计第章(2)对低通滤波器单位脉冲响应乘以cosn可得到高通滤波器单位脉冲响应：矩形窗：升余弦窗：有限脉冲响应(FIR)数字滤波器的设计第章 改进升余弦窗：布莱克曼窗：有限脉冲响应(FIR)数字滤波器的设计第章题12图12 利用窗函数（哈明窗）法设计一数字微分器，逼近题12图所示的理想微分器特性，并绘出其幅频特性

16、。有限脉冲响应(FIR)数字滤波器的设计第章解:(1)由于连续信号存在微分，而时域离散信号和数字信号的微分不存在，因而本题要求设计的数字微分器是指用数字滤波器近似实现模拟微分器，即用数字差分滤波器近似模拟微分器。下面先推导理想差分器的频率响应函数。设模拟微分器的输入和输出分别为x(t)和y(t)，即令x(t)=ejt，则y(t)=jket=jkx(t)对上式两边采样（时域离散化），得到有限脉冲响应(FIR)数字滤波器的设计第章其中=T。将x(nT)和y(nT)分别作为数字微分器的输入和输出序列，并用Hd(ej)表示数字理想微分器的频率响应函数，则即有限脉冲响应(FIR)数字滤波器的设计第章根据

17、题12图所给出的理想特性可知所以应取k=T，所以Hd(ej)=j取群延时=(N 1)/2，则逼近频率响应函数应为 Hd(ej)=jej=ej(/2)有限脉冲响应(FIR)数字滤波器的设计第章有限脉冲响应(FIR)数字滤波器的设计第章设FIR 滤波器h(n)长度为N，一般取=(N 1)/2。加窗后得到我们知道，微分器的幅度响应随频率增大线性上升，当频率=时达到最大值，所以只有N 为偶数的情况4才能满足全频带微分器的时域和频域要求。因为N 是偶数，=N/21/2=正整数1/2，上式中第一项为0，所以有限脉冲响应(FIR)数字滤波器的设计第章式就是用窗函数法设计的FIR 数字微分器的单位脉冲响应的通

18、用表达式，且具有奇对称特性h(n)=h(N 1n)。选定滤波器长度N 和窗函数类型，就可以直接按式得到设计结果。当然，也可以用频率采样法和等波纹最佳逼近法设计。本题要求的哈明窗函数：有限脉冲响应(FIR)数字滤波器的设计第章将式代入式得到h(n)的表达式：(2)对3种不同的长度N=20，40和41，用MATLAB 计算单位脉冲响应h(n)和幅频特性函数，并绘图的程序ex712.m 如下：有限脉冲响应(FIR)数字滤波器的设计第章%ex712.m：用哈明窗设计线性相位FIR 微分器clear all；close all；N1=20；n=0：N1 1；tou=(N1 1)/2；h1n=sin(n

19、tou)*pi)./(pi*(n-tou).2).*(hamming(N1)；N2=40；n=0：N2 1；tou=(N2 1)/2；h2n=sin(n tou)*pi)./(pi*(n-tou).2).*(hamming(N2)；N3=41；n=0：N3 1；tou=(N3 1)/2；h3n=sin(n tou)*pi)./(pi*(n tou).2).*(hamming(N3);h3n(N3 1)/2+1)=0；%因为该点分母为零，无定义，所以赋值0%以下为绘图部分（省略）有限脉冲响应(FIR)数字滤波器的设计第章程序运行结果即数字微分器的单位脉冲响应和幅频特性函数曲线如题12解图所示。由

20、图可见，当滤波器长度N 为偶数时，逼近效果好。但N=奇数时（本程序中N=41），逼近误差很大。这一结论与教材给出的理论一致（对第二类线性相位滤波器，N=奇数时不能实现高通滤波特性）。有限脉冲响应(FIR)数字滤波器的设计第章题12解图有限脉冲响应(FIR)数字滤波器的设计第章也可以采用调用等波纹最佳逼近法设计函数remez 来设计FIR 数字微分器的方法。hn=remez(N 1，f，m，defferentiator)设计N 1阶FIR 数字微分器，返回的单位脉冲响应向量hn具有奇对称特性。在大多数工程实际中，仅要求在频率区间0p上逼近理想微分器的频率响应特性，而在区间p 上频率响应特性不作要

21、求，或要求为零。对微分器设计，在区间p 上频率响应特性要求为零时，调用参数f=0,p/,(p+B)/,1，m=0,p/,0,0，其中B 为过渡带宽度（即无关区），p不能太靠近，B 也不能太小，否则设计可能失败。调用等波纹最佳逼近法设计函数remez 设计本题要求的FIR 数字微分器的程序ex712b.m如下：有限脉冲响应(FIR)数字滤波器的设计第章%ex712b.m：调用remez 函数设计FIR 微分器Wp=0.9；B=0.09；%设置微分器边界频率（关于 归一化）N=40；f=0，wp，wp+B，1；m=0，wp，0,0；hn=remez(N-1，f，m，defferentiator)；

22、%调用remez 函数设计FIR 微分器%以下为绘图部分（省略）请读者运行该程序，观察设计效果。有限脉冲响应(FIR)数字滤波器的设计第章13.用窗函数法设计一个线性相位低通FIRDF，要求通带截止频率为/4 rad，过渡带宽度为8/51 rad，阻带最小衰减为45 dB。（1）选择合适的窗函数及其长度，求出h(n)的表达式。（2*）用MATLAB 画出损耗函数曲线和相频特性曲线。解：（1）根据教材7.2.2节所给步骤进行设计。根据对阻带衰减及过渡带的指标要求，选择窗函数的类型，并估计窗口长度N。由习题9中教材表7.2.2，本题应选择哈明窗。因为过渡带宽度Bt=8/51，所以窗口长度N 为N6

23、.6/Bt=42.075，取N=43。窗函数表达式为有限脉冲响应(FIR)数字滤波器的设计第章 构造希望逼近的频率响应函数Hd(ej):式中有限脉冲响应(FIR)数字滤波器的设计第章 求hd(n):加窗:有限脉冲响应(FIR)数字滤波器的设计第章（2）调用MATLAB 函数设计及绘图程序ex713.m 如下：%ex713.m：调用fir1 设计线性相位低通FIR 滤波器并绘图 wp=pi/4；Bt=8*pi/51；wc=wp+Bt/2；N=ceil(6.6*pi/Bt)；hmn=fir1(N-1，wc/pi，hamming(N)rs=60；a=1；mpplot(hmn，a，rs)%调用自编函数

24、mpplot 绘制损耗函数和相频特性曲线程序运行结果即损耗函数和相频特性曲线如题13解图所示，请读者运行程序查看h(n)的数据。有限脉冲响应(FIR)数字滤波器的设计第章题13解图有限脉冲响应(FIR)数字滤波器的设计第章14.要求用数字低通滤波器对模拟信号进行滤波，要求：通带截止频率为10 kHz，阻带截止频率为22 kHz，阻带最小衰减为75 dB，采样频率为Fs=50 kHz。用窗函数法设计数字低通滤波器。（1）选择合适的窗函数及其长度，求出h(n)的表达式。（2*）用MATLAB 画出损耗函数曲线和相频特性曲线。解：（1）根据教材7.2.2节所给步骤进行设计。根据对阻带衰减及过渡带的指

25、标要求，选择窗函数的类型，并估计窗口长度N。本题要求设计的FIRDF 指标：有限脉冲响应(FIR)数字滤波器的设计第章通带截止频率：阻带截止频率：阻带最小衰减：s=75 dB有限脉冲响应(FIR)数字滤波器的设计第章由习题9中教材表7.2.2可知，本题应选凯塞窗(=7.865)。窗口长度N10/Bt=10/(sp)=20.833，取N=21。窗函数表达式为，=7.865 构造希望逼近的频率响应函数Hd(ej):有限脉冲响应(FIR)数字滤波器的设计第章 求hd(n):加窗:有限脉冲响应(FIR)数字滤波器的设计第章（2）调用MATLAB 函数设计及绘图程序ex714.m 如下：%ex714.m

26、：调用fir1 设计线性相位低通FIR 滤波器并绘图Fs=50000；fp=10000；fs=22000；rs=75；wp=2*pi*fp/Fs；ws=2*pi*fs/Fs；Bt=ws-wp；wc=(wp+ws)/2；N=ceil(10*pi/Bt)；hmn=fir1(N-1，wc/pi，kaiser(N，7.865)；rs=100；a=1；mpplot(hmn，a，rs)%调用自编函数mpplot 绘制损耗函数和相频特性曲线程序运行结果即损耗函数和相频特性曲线如题14解图所示，请读者运行程序查看h(n)的数据。有限脉冲响应(FIR)数字滤波器的设计第章题14解图有限脉冲响应(FIR)数字滤波

27、器的设计第章15 利用频率采样法设计线性相位FIR 低通滤波器，给定N=21，通带截止频率c=0.15 rad。求出h(n)，为了改善其频率响应（过渡带宽度、阻带最小衰减），应采取什么措施？解:(1)确定希望逼近的理想低通滤波频率响应函数Hd(ej)：有限脉冲响应(FIR)数字滤波器的设计第章其中,a=(N 1)/2=10。采样：有限脉冲响应(FIR)数字滤波器的设计第章 求h(n)：有限脉冲响应(FIR)数字滤波器的设计第章因为所以有限脉冲响应(FIR)数字滤波器的设计第章损耗函数曲线绘图程序ex715.m 如下：%程序ex715.m N=21；n=0：N-1；hn=(1+2*cos(2*p

28、i*(n-10)/N)/N；rs=20；a=1；mpplot(hn，a，rs)%调用自编函数mpplot 绘制损耗函数和相频特性曲线运行程序绘制损耗函数曲线如题15解图所示，请读者运行程序查看hn的数据。为了改善阻带衰减和通带波纹，应加过渡带采样点，为了使边界频率更精确，过渡带更窄，应加大采样点数N。有限脉冲响应(FIR)数字滤波器的设计第章题15解图有限脉冲响应(FIR)数字滤波器的设计第章16 重复题15，但改为用矩形窗函数法设计。将设计结果与题15进行比较。解：直接调用fir1 设计，程序为ex716.m。%调用fir1 求解16题的程序ex716.m N=21；wc=0.15；hn=f

29、ir1(N-1，wc，boxcar(N)；%选用矩形窗函数(与上面求解中相同)rs=20；a=1；mpplot(hn，a，rs)%调用自编函数mpplot 绘制损耗函数和相频特性曲线运行程序绘制损耗函数曲线如题16解图所示。与题15解图比较，过渡带宽度相同，但矩形窗函数法设计的FIRDF阻带最小衰减约为20 dB，而15题设计结果约为16 dB。有限脉冲响应(FIR)数字滤波器的设计第章题16解图有限脉冲响应(FIR)数字滤波器的设计第章17 利用频率采样法设计线性相位FIR 低通滤波器，设N=16，给定希望逼近的滤波器的幅度采样值为解:由希望逼近的滤波器幅度采样Hdg(k)可构造出Hd(ej

30、)的采样Hd(k)：有限脉冲响应(FIR)数字滤波器的设计第章有限脉冲响应(FIR)数字滤波器的设计第章 18 利用频率采样法设计线性相位FIR 带通滤波器，设N=33，理想幅度特性Hd()如题18图所示。题18图有限脉冲响应(FIR)数字滤波器的设计第章解:由题18图可得到理想幅度采样值为有限脉冲响应(FIR)数字滤波器的设计第章有限脉冲响应(FIR)数字滤波器的设计第章19*设信号x(t)=s(t)+v(t)，其中v(t)是干扰，s(t)与v(t)的频谱不混叠，其幅度谱如题19*图所示。要求设计数字滤波器，将干扰滤除，指标是允许|s(f)|在0f15 kHz 频率范围中幅度失真为2%（1=

31、0.02）；f 20 kHz，衰减大于40 dB（2=0.01）；希望分别设计性价比最高的FIR 和IIR 两种滤波器进行滤除干扰。请选择合适的滤波器类型和设计方法进行设计，最后比较两种滤波器的幅频特性、相频特性和阶数。有限脉冲响应(FIR)数字滤波器的设计第章题19*图有限脉冲响应(FIR)数字滤波器的设计第章解：本题以模拟频率给定滤波器指标，所以，程序中先要计算出对应的数字边界频率，然后再调用MATLAB工具箱函数fir1 设计数字滤波器。由题意确定滤波器指标（边界频率以模拟频率给出）：fp=15 kHz，1=0.02，p=20 lg dB fp=20 kHz，2=0.01，s=40 dB

32、（1）确定相应的数字滤波器指标。根据信号带宽，取系统采样频率Fs=80 kHz。，1=0.02，p=20 lg dB，2=0.01，s=40 dB 有限脉冲响应(FIR)数字滤波器的设计第章（2）设计数字低通滤波器。为了设计性价比最高的FIR 和IIR 滤波器，IIR 滤波器选择椭圆滤波器，FIR 滤波器采用等波纹最佳逼近法设计。设计程序为ex719.m。%ex719.m：设计性价比最高的FIR 和IIR 滤波器Fs=80000；fp=15000；fs=20000；data1=0.02；rp=20*log10(1 data1)/(1+1)；data2=0.01；rs=40；wp=2*fp/Fs

33、；ws=2*fs/Fs；%计算数字边界频率（关于归一化）%椭圆DF 设计有限脉冲响应(FIR)数字滤波器的设计第章Ne，wpe=ellipord(wp，ws，rp，rs)；%调用ellipord 计算椭圆DF 阶数N 和通带截止频率wpBe，Ae=ellip(Ne，wpe，rs，wp)；%调用ellip 计算椭圆DF 系统函数系数向量Be 和Ae%用等波纹最佳逼近法设计FIRDF f=wp，ws；m=1,0；rip=data1,data2；Nr，fo，mo，w=remezord(f，m，rip)；hn=remez(Nr，fo，mo，w)；%以下为绘图部分（省略）有限脉冲响应(FIR)数字滤波器

34、的设计第章程序运行结果：椭圆DF 阶数Ne=5，损耗函数曲线和相频特性曲线如题图19*解图（a）所示。采用等波纹最佳逼近法设计的FIRDF 阶数Nr=29，损耗函数曲线和相频特性曲线如题19*解图(b)图所示。由图可见，IIRDF 阶数低得多，但相位特性存在非线性，FIRDF 具有线性相位特性。有限脉冲响应(FIR)数字滤波器的设计第章题19*解图有限脉冲响应(FIR)数字滤波器的设计第章20*.调用MATLAB 工具箱函数fir1 设计线性相位低通FIR 滤波器，要求希望逼近的理想低通滤波器通带截止频率c=/4 rad，滤波器长度N=21。分别选用矩形窗、Hanning 窗、Hamming

35、窗和Blackman 窗进行设计，绘制用每种窗函数设计的单位脉冲响应h(n)及其损耗函数曲线，并进行比较，观察各种窗函数的设计性能。有限脉冲响应(FIR)数字滤波器的设计第章题20*解图有限脉冲响应(FIR)数字滤波器的设计第章解：本题设计程序ex720.m 如下：%ex720.m：调用fir1 设计线性相位低通FIR 滤波器clear；close all；N=21；wc=1/4；n=0：20；hrn=fir1(N 1，wc，boxcar(N)；%用矩形窗函数设计hnn=fir1(N 1，wc，hanning(N)；%用hanning 窗设计hmn=fir1(N 1，wc，hamming(N)

36、；%用hamming 窗函数设计hbn=fir1(N 1，wc，blackman(N)；%用blackman 窗函数设计%以下为绘图部分（省略）有限脉冲响应(FIR)数字滤波器的设计第章程序运行结果：用矩形窗、Hanning 窗、Hamming 窗和Blackman 窗设计的单位脉冲响应h(n)及其损耗函数曲线如题20*解图所示。由图可见，滤波器长度N 固定时，矩形窗设计的滤波器过渡带最窄，阻带最小衰减也最小；blackman 窗设计的滤波器过渡带最宽，阻带最小衰减最大。21*.将要求改成设计线性相位高通FIR 滤波器，重作题20。有限脉冲响应(FIR)数字滤波器的设计第章解：本题的设计程序除

37、了在每个fir1 函数的调用参数中加入滤波器类型参数“high”外，与第20题的程序完全相同，请读者修改并运行程序，完成本题。22*.调用MATLAB 工具箱函数remezord 和remez 设计线性相位低通FIR 滤波器，实现对模拟信号的采样序列x(n)的数字低通滤波处理。指标要求：采样频率为16 kHz；通带截止频率为4.5 kHz，通带最小衰减为1 dB；阻带截止频率为6 kHz，阻带最小衰减为75 dB。列出h(n)的序列数据，并画出损耗函数曲线。解：本题设计程序ex722.m 如下：有限脉冲响应(FIR)数字滤波器的设计第章%ex722.m：调用remezord 和remez 设计

38、线性相位低通FIR滤波器Fs=16000；f=4500，6000；%采样频率，边界频率为模拟频率(Hz)m=1，0；rp=1；rs=75；dat1=(10(rp/20)-1)/(10(rp/20)+1)；dat2=10(-rs/20)；rip=dat1，dat2；M，fo，mo，w=remezord(f，m，rip，Fs)；M=M+1；%边界频率为模拟频率(Hz)时必须%加入采样频率Fs hn=remez(M，fo，mo，w)%以下为绘图部分（省略）有限脉冲响应(FIR)数字滤波器的设计第章程序运行结果：h（n）=0.0023 0.0026 0.0207 0.0131 0.0185 0.003

39、2 0.0278 0.0306 0.0176 0.0705 0.0402 0.1075 0.2927 0.6227 0.2927 0.1075 0.0402 0.0705 0.0176 0.0306 0.0278 0.0032 0.0185 0.0131 0.0207 0.0026 0.0023单位脉冲响应h(n)及其损耗函数曲线如题22*解图所示。有限脉冲响应(FIR)数字滤波器的设计第章题22*解图有限脉冲响应(FIR)数字滤波器的设计第章23*.调用MATLAB 工具箱函数remezord 和remez 设计线性相位高通FIR 滤波器，实现对模拟信号的采样序列x(n)的数字高通滤波处理。

40、指标要求：采样频率为16 kHz；通带截止频率为5.5 kHz，通带最小衰减为1dB；过渡带宽度小于等于3.5 kHz，阻带最小衰减为75 dB。列出h(n)的序列数据，并画出损耗函数曲线。解：滤波器的阻带截止频率fs=5500 3500=2000 Hz。本题设计程序ex723.m 如下：有限脉冲响应(FIR)数字滤波器的设计第章%ex723.m：调用remezord 和remez 设计线性相位高通FIR 滤波器Fs=16000；f=2000，5500；%采样频率，边界频率为模拟频率(Hz)m=0，1；rp=1；rs=75；dat1=(10(rp/20)-1)/(10(rp/20)+1)；da

41、t2=10(-rs/20)；rip=dat2，dat1；M，fo，mo，w=remezord(f，m，rip，Fs)；%边界频率为模拟频率(Hz)时必须加入采样频率Fs hn=remez(M，fo，mo，w)程序运行结果：滤波器长度为N=M+1=11，有限脉冲响应(FIR)数字滤波器的设计第章单位脉冲响应h(n)及其损耗函数曲线如题23*解图所示，请读者运行程序查看h(n)的数据。题23*解图有限脉冲响应(FIR)数字滤波器的设计第章24*.用窗函数法设计一个线性相位低通FIR 滤波器，要求通带截止频率为0.3 rad，阻带截止频率为0.5 rad，阻带最小衰减为40 dB。选择合适的窗函数及

42、其长度，求出并显示所设计的单位脉冲响应h(n)的数据，并画出损耗函数曲线和相频特性曲线，请检验设计结果。试不用fir1函数，直接按照窗函数设计法编程设计。解：直接按照窗函数设计法的设计程序ex724.m 如下：有限脉冲响应(FIR)数字滤波器的设计第章%ex724.m：直接按照窗函数设计法编程设计线性相位低通FIR 滤波器wp=0.3*pi；ws=0.5*pi；rs=40；%指标参数Bt=ws-wp；%过渡带宽度N=ceil(6.2*pi/Bt)；%选hanning 窗，求wn 长度N wc=(wp+ws)/2；r=(N-1)/2；%理想低通截止频率wc n=0：N-1；hdn=sin(wc*

43、(n-r)./(pi*(n-r)；%计算理想低通的hdnhdn(16)=wc/pi；%在n=(N-1)/2=15 点为0/0 型，直接赋值wn=0.5*(1-cos(2*pi*n/(N-1)；%求窗函数序列wn hn=hdn.*wn%加窗%以下为绘图部分（省略）有限脉冲响应(FIR)数字滤波器的设计第章题24*解图程序运行结果：单位脉冲响应h(n)及其损耗函数曲线如题24*解图所示，请读者运行程序查看h(n)的数据。有限脉冲响应(FIR)数字滤波器的设计第章25*.调用MATLAB 工具箱函数fir1 设计线性相位高通FIR滤波器。要求通带截止频率为0.6 rad，阻带截止频率为0.45，通带最大衰减为0.2 dB，阻带最小衰减为45 dB。显示所设计的单位脉冲响应h(n)的数据，并画出损耗函数曲线。解：本题设计程序ex725.m 如下：%ex725.m：调用fir1 设计线性相位高通FIR 滤波器wp=0.6*pi；ws=0.45*pi；rs=45；%指标参数wc=(wp+ws)/2；%理想低通截止频率wc Bt=wp-ws；%过渡带宽度N1=ceil(6.6*pi/Bt)；%hamming 窗w(n)长度