北理工大学数字信号处理实验报告

1、1 DFT用于分析信号频谱。一.实验的目的1.加深对DFT原理的理解。2.通过DFT分析信号的频谱。3.深刻理解DFT分析信号频谱的原理，分析实现过程中的现象和解决方法。二。实验设备和环境计算机和MATLAB软件环境。三、实验原理1.DFT与DTFT的关系频率范围内的N个等距点上的有限序列的离散时间傅立叶变换的N个采样值可以用下面的公式表示:从上式可以看出，序列x(n)的n点DFT X(k)实际上是x(n)序列的DTFT在n个等距频点上的样本X(k)。2.用离散傅立叶变换寻找DTFT方法1:从X (k)中回收的方法:方法二:但是在实际的MATLAB计算中，上述插值运算不一定是最好的方式。由于D

2、FT是DTFT的采样值，并且两个相邻频率采样点之间的距离是2/N，如果我们增加数据的长度N，DFT谱将更精细，其包络将更接近DTFT的结果，这样我们就可以使用DFT来近似DTFT。如果没有更多的数据，您可以通过填充零来增加数据长度。3.用DFT分析连续时间信号的频谱。用计算机分析连续时间信号的频谱，首先要对连续时间信号进行离散化。这里需要两个操作:采样和截断。DFT分析连续非周期信号频谱的步骤总结如下:(1)确定时域采样间隔t，得到离散序列x(n)；(2)确定截距长度m，得到m个点的离散序列的，其中是窗函数。(3)确定频域采样点数n，并要求。(4)通过FFT计算离散序列的N点DFT，得到。(5

3、)根据公式(2-6)，计算采样点的近似值。用上述方法计算的频谱，我们需要注意以下三个问题:(1)频谱混叠。如果采样定理的条件不满足，频谱就会出现混叠误差。对于具有无限频谱的信号，应考虑覆盖大部分主要频率成分的范围。(2)栅栏效应和光谱分辨率。用DFT计算频谱，结果只有N个谱特征值，样本值之间的频谱未知。比如透过栅栏观察频谱，就叫“栅栏效应”。频谱分辨率与记录长度成反比，所以要提高频谱分类率，就要增加记录时间。(3)频谱泄漏。截断信号会将窗函数的频谱引入信号频谱，表现为频谱泄漏。解决这个问题的主要方法是使用旁瓣较小的窗函数，频谱泄漏和窗函数都会引起误差。因此，需要合理选择采样间隔和截取长度，必要

4、时添加合适的窗口。对于连续时间周期信号，我们用计算机计算时，总是要对其进行截断，序列总是有限的，所以还是可以用上面的方法近似计算。四。实验能力1.已知x(n)=1，1，1，2，满足以下要求:(1)计算其DTFT，画出-pi，pi区间的波形。源代码:n = 0:3；x =2-1 1 1；w =-pi:0.01 * pi:pi；x = x * exp(-1j * n * w)；支线剧情(211)；plot(w，ABS(X)；xlabel( Omega/ pi )；标题(“量级”)；轴(-pi pi 0 4)；子情节(212)；plot(w，角度(X)/pi)；xlabel( Omega/ pi )

5、；轴(-pi pi-1 1)；标题(“阶段”)；实验结果:(2)计算4点DFT，并在(1)中绘制的图中显示结果。源代码:n = 0:3；x =2-1 1 1；w =-pi:0.01 * pi:pi；x = x * exp(-1j * n * w)；支线剧情(211)；plot(w，ABS(X)；xlabel( Omega/ pi )；标题(“量级”)；轴(-pi pi 0 4)；坚持住；w1 =-pi:0.5 * pi:0.5 * pi；x1 = x * exp(-1j * n * w1)；阀杆(w1，abs(X1)，填充)；子情节(212)；plot(w，角度(X)/pi)；坚持住；阀杆(w

6、1，角度(X1)/pi，“填充”)；xlabel( Omega/ pi )；轴(-pi pi-1 4)；标题(“阶段”)；实验结果:(3)对x(n)补零，计算64点DFT，并显示结果。源代码:n = 0:3；n1 = 0:63；x =2-1 1 1；x1=2 -1 1 1个零(1，60)；w =-pi:0.01 * pi:pi；w1 =-pi:pi/32:pi-pi/32；x = x * exp(-1j * n * w)；x1 = x1 * exp(-1j * n1 * w1)；支线剧情(211)；plot(w，ABS(X)；xlabel( Omega/ pi )；标题(“量级”)；轴(-pi

7、 pi 0 4)；坚持住；阀杆(w1，abs(X1)，填充)；子情节(212)；plot(w，角度(X)/pi)；坚持住；阀杆(w1，角度(X1)/pi，“填充”)；xlabel( Omega/ pi )；轴(-pi pi-1 4)；标题(“阶段”)；(4)根据实验结果，分析DTFT是否可以用DFT计算，如果可以，如何实现。DTFT可以通过DFT计算。从上图可以看出，n的值越大，DFT的包络越像DTFT的曲线。因此，我们可以直接得到答案，即当我们得到足够大的n时，我们可以用DFT来表示DTFT。2.调查顺序(1)在，DFT用于估计的频谱；零填充被扩展到长度为100点的序列。用DFT估计频谱。画

8、出相应的波形。源代码:n1 = 0:10；x1=(cos(0.48*pi。*n1)+cos(0.52*pi。* n1)；x1 = FFT(x1)；支线剧情(211)；xlabel( Omega/ pi )；标题(“量级”)；阀杆(n1，abs(X1)，填充)；n = 0:100；ni = 0:9；x=(cos(0.48*pi。*ni)+cos(0.52*pi。*ni)、零(1，91)；x = FFT(x)；子情节(212)；xlabel( Omega/ pi )；标题(“量级”)；阀杆(n/10，abs(X)，填充)；(2)在，用DFT估计x(n)的频谱，画出波形。源代码:n1 = 0:100

9、；x1=cos(0.48*pi。*n1)+cos(0.52*pi。* n1)；x1 = FFT(x1)；xlabel( Omega/ pi )；标题(“量级”)；阀杆(n1/100*2*pi，abs(X1)，填充)；(3)根据实验结果，分析如何提高光谱分辨率。当n增加并且n增加得足够大时，信号的频谱看起来更像原始实际信号的频谱。3.已知的信号，包括，和。从的表达式可以看出，它包含三个频率的正弦波。但从其时域波形(图E2-1)来看，似乎是一个正弦信号。利用DFT做频谱分析，确定合适的参数，使得到的频谱的频率分辨率满足要求。源代码:n = 0:0.1:4.9；x = 0.15 * sin(2 *

10、pi * n)+sin(4 * pi * n)-0.1 * sin(6 * pi * n)；x = FFT(x)；stem(n，abs(X)， filled )；4.用DFT近似分析连续时间信号的频谱(振幅谱)。分析不同采样间隔和截取长度的计算结果，最终确定合适的参数。源代码:t = 0:0.1:99.9；w =-pi:0.01:pi；x = exp(0.1 * t)；x = x * exp(-1j * t * w)；plot(w，ABS(X)；实验的主要结论、遇到的问题、解决方案、收获和经验。主要结论:1)可以用DFT代替DTFT，但是有一个条件，就是DFT中的值必须足够大，并且DFT的包络

11、能够代表DTFT的曲线。2)当时域信号离散时，我们可以增加采样时间，时间越长，离散信号的频谱越接近真实连续信号的频谱。遇到的问题:用MATLAB做DFT时，不能及时完成，信号出来时相位有延迟。后来发现时间采样不对。以及获得的经验:我们用MATLAB实际上看到了DFT和DTFT的区别和联系，课本上的知识不再只是想象的舞台。我们也可以用计算机软件来模拟，从其他方面给我们一个解决问题的思路和方法。实验IIR数字滤波器的设计一.实验的目的掌握用脉冲响应不变性方法和双线性变换方法设计IIR数字滤波器的原理和具体方法。了解数字滤波器和模拟滤波器的技术指标转换。掌握脉冲响应不变性方法和双线性变换方法在设计I

12、IR数字滤波器中的优缺点，以及它们的适用范围。二。实验设备和环境计算机和MATLAB软件环境三、实验原理IIR滤波器具有无限持续时间脉冲响应，而模拟滤波器通常具有无限持续时间的脉冲响应，因此，它与模拟滤波器相匹配。IIR滤波器设计的基本方法是首先设计一个合适的模拟滤波器，然后通过复值映射将模拟滤波器转换成数字滤波器。1.数字滤波器和模拟滤波器的一些指标通带、过渡带、通带响应容差、阻带容差、通带波动、阻带衰减等。2.模拟原型滤波器IIR滤波器设计方法数字滤波器是从现有的模拟滤波器获得的。我们称这些模拟滤波器为原型滤波器。常用的模拟原型滤波器包括巴特沃斯滤波器、切比雪夫(型和型)滤波器和椭圆滤波器

13、等。巴特沃斯过滤器切舍夫低通滤波器一、切比雪夫I型低通滤波器B.切比雪夫II低通滤波器椭圆滤波器(4)三种过滤器的比较给定相同的设计指标，椭圆滤波器所需阶数N最低，其次是切比雪夫滤波器和巴特沃兹滤波器。如果所需阶数相同，切比雪夫滤波器的过渡带比巴特沃兹滤波器陡，椭圆滤波器的过渡带比切比雪夫滤波器陡。但是从通带的相位响应来看，虽然椭圆滤波器提供了最好的幅度平方响应，但是通带的相位响应是非线性的，而巴特沃兹滤波器在通带中具有相当线性的相位，切比雪夫滤波器的相位特性介于两者之间。因此，在实际设计中，滤波器的选择应取决于实际用途和指标要求。3.从模拟滤波器到数字滤波器的转换模拟滤波器到数字滤波器的变换

14、是进一步从H a (s)得到H(z)，也就是从S平面到Z平面的变换，这就需要数字滤波器模仿模拟滤波器的特性。模拟滤波器到数字滤波器的转换方法有两种常用方法:脉冲响应不变量法和双线性变换法。3.1脉冲响应不变方法基本原理:从时域响应出发，数字滤波器的单位冲激响应h(n)模仿模拟滤波器的单位冲激响应ha(t)，h(n)等于ha(t)的采样值。3.2双线性变换法基本原理:从频率响应出发，直接将数字滤波器的频率响应近似为模拟滤波器的频率响应，进而得到数字滤波器的系统函数H(z)。利用MATLAB函数直接设计IIR数字滤波器我们介绍了用butter函数、cheby1函数、cheby2函数和ellip函数

15、设计模拟滤波器的方法。这些函数也可以直接用于设计IIR数字滤波器。由于双线性变换的优点，这些函数都采用双线性变换。这些函数的具体用法如下:b，a = butter(N，Wn)设计一个低通巴特沃兹数字滤波器。b，a = cheby1(N，Rp，Wn)设计低通Chebyshev 数字滤波器b，a = cheby2(N，As，Wn，)设计低通Chebyshev数字滤波器b，a = ellip(N，Rp，As，Wn)设计低通椭圆数字滤波器b，a = butter(N，Wn， ftype )设计巴特沃兹数字滤波器b，a = cheby1(N，Rp，Wn， ftype )设计切比雪夫型数字滤波器b，a =

16、 cheby2(N，As，Wn， ftype )设计切比雪夫型数字滤波器b，a= ellip设计的椭圆数字滤波器(n，rp，as，wn， ftype )“Ftype”表示数字滤波器的类型，可选选项有“high”、“low”和“stop”，分别表示高通、低通和带阻。上述函数中的参数n和Wn可以通过使用函数buttord、cheb1ord、cheb2ord和ellipord得到，在设计指标已知的情况下，可以给出滤波器的阶数n和截止频率Wn。这些函数的具体用法如下:N，Wn = buttord(Wp，Wst，Rp，As)根据数字滤波器的设计指标，给出数字巴特沃兹滤波器的阶数N和截止频率Wn；N，Wn

17、 = cheb1ord(Wp，Wst，Rp，As)根据数字滤波器的设计指标给出数字切比雪夫型滤波器的阶数N和截止频率Wn；N，Wn = cheb2ord(Wp，Wst，Rp，As)根据数字滤波器的设计指标，给出了数字切比雪夫型滤波器的阶数N和截止频率Wn；N，Wn = ellipord(Wp，Wst，Rp，As)根据数字滤波器的设计指标，给出数字椭圆滤波器的阶数N和截止频率Wn。其中Wp和Wst分别表示通带和阻带截止频率，Rp表示通带波动，As表示阻带衰减。四。实验能力将频率设置为，并设计一个数字低通滤波器，要求:要求分别设计Butterworth、Chebyshev I、Chebyshev

18、II和椭圆模拟原型滤波器，结合脉冲响应不变性方法和双线性变换方法分别设计。结合实验结果，讨论了用上述方法设计的数字滤波器能否满足给定指标的要求，分析了用脉冲响应不变性方法和双线性变换方法设计的IIR数字滤波器的优缺点和适用范围。巴特沃斯wp = 2 * pi * 1000ws = 2 * pi * 1500RP = 1；As = 15Fs = 10000OmegaP = wpOmegaS = wsep=sqrt(10(rp/10)-1)；纹波= sqrt(1/(1+EP * EP)；attn=1/(10(as/20)；N，OmegaC= but ord(OmegaP，OmegaS，Rp，aS，

19、 s )；z0，p0，k0= buttap(N)；p = P0 *c；z = z0 * OmegaCk=k0*omegacn；ba = k * real(poly(z)；aa = real(poly(p)；bd，ad=impinvar(ba，aa，Fs)；H，w=freqz(bd，ad，1000，整数)；H =(H(1:1:501)；w =(w(1:1:501)；mag = ABS(H)；db = 20 * log 10(mag+EPS)/max(mag)；pha =角度(H)；grd=grpdelay(bd，ad，w)；支线剧情(221)；绘图(w/pi，mag)；标题(“A”)；y lab

20、el( | H | )；网格打开；支线剧情(2，2，3)；绘图(w/pi，db)；标题( A(dB)；xlabel( f(pi)；y label( dB )；轴(0，1，-100，50)；网格打开；支线剧情(222)；plot(w/pi，PHA/pi)；标题(“阶段”)；xlabel(“”)；ylabel( Phase(pi)；轴(0，1，-1，1)；网格打开；支线剧情(224)；绘图(w/pi，grd)；标题( t )；xlabel( pi )；y label( x )；轴(0，1，0，20)；网格打开；切舍夫型wp = 0.2 * piws = 0.3 * piRP = 1；As = 15

21、Fs = 10000OmegaP = wp * FsOmegaS = ws * Fsep=sqrt(10(rp/10)-1)；纹波= sqrt(1/(1+EP * EP)；attn=1/(10(as/20)；N，OmegaC=cheb1ord(OmegaP，OmegaS，Rp，aS， s )；z0，p0，k0=cheb1ap(N，Rp)；p = P0 *c；z = z0 * OmegaCk=k0*omegacn；ba = k * real(poly(z)；aa = real(poly(p)；bd，ad=impinvar(ba，aa，Fs)；H，w=freqz(bd，ad，1000，整数)；H

22、=(H(1:1:501)；w =(w(1:1:501)；mag = ABS(H)；db = 20 * log 10(mag+EPS)/max(mag)；pha =角度(H)；grd=grpdelay(bd，ad，w)；支线剧情(2，2，1)；绘图(w/pi，mag)；标题( | A | )；xlabel(“”)；y label( | H | )；轴(0，1，0，1.1)；set(gca， XTickMode ， manual ， XTick ，0，0.4，0.6，1)；网格打开；set(gca， XTickMode ， manual ， YTick ，0，Attn，Ripple，1)；支线剧情

23、(2，2，3)；绘图(w/pi，db)；title( | A |(dB)；xlabel( Frep(pi)；y label( dB )；轴(0，1，-100，50)；网格打开；支线剧情(2，2，2)；plot(w/pi，PHA/pi)；标题(“阶段”)；xlabel(“”)；ylabel( Phase(pi)；轴(0，1，-1，1)；网格打开；支线剧情(2，2，4)；绘图(w/pi，grd)；标题(下岗)；xlabel( Frep(pi)；y label( x )；轴(0，1，0，20)；网格打开；切舍夫型wp = 0.2 * piws = 0.3 * piRP = 1；As = 15Fs =

24、 10000t = 1/Fs；OmegaP =(2/T)* tan(WP/2)；OmegaS =(2/T)* tan(ws/2)；ep=sqrt(10(rp/10)-1)；纹波= sqrt(1/(1+EP * EP)；attn=1/(10(as/20)；N，OmegaC=cheb2ord(OmegaP，OmegaS，Rp，aS， s )；z0，p0，k0=cheb2ap(N，As)；p = P0 *c；z = z0 * OmegaCk=k0*omegacn；ba0 = real(poly(z0)；ba0 = k0 * ba0aa0 =实数(poly(P0)；ba = real(poly(z)；

25、ba = k * baaa = real(poly(p)；bd，ad=双线性(ba，aa，Fs)；bd1，ad1=双线性(ba0，aa0，Fs/omega c)；子情节(1，1，1)H，w=freqz(bd，ad，1000，整数)；H =(H(1:1:501)；w =(w(1:1:501)；mag = ABS(H)；db = 20 * log 10(mag+EPS)/max(mag)；pha =角度(H)；grd=grpdelay(bd，ad，w)；支线剧情(2，2，1)；绘图(w/pi，mag)；标题(“A”)；xlabel(“”)；y label( | A | )；set(gca， XTi

26、ckMode ， manual ， XTick ，0，0.4，0.6，1)；网格打开；set(gca， XTickMode ， manual ， YTick ，0，Attn，Ripple，1)；支线剧情(2，2，3)；绘图(w/pi，db)；title( | A |(dB)；xlabel( freq(rad/s)；y label( dB )；网格打开；支线剧情(2，2，2)；plot(w/pi，PHA/pi)；标题(“阶段”)；xlabel(“”)；y label( pi )；网格打开；支线剧情(2，2，4)；绘图(w/pi，grd)；标题(下岗)；xlabel( freq(rad/s)；伊拉

27、贝尔(尼翁)；网格打开；椭圆模拟原型wp = 0.2 * piws = 0.3 * piRP = 1；As = 15Fs = 10000t = 1/Fs；OmegaP =(2/T)* tan(WP/2)；OmegaS =(2/T)* tan(ws/2)；ep=sqrt(10(rp/10)-1)；纹波= sqrt(1/(1+EP * EP)；attn=1/(10(as/20)；N，OmegaC=ellipord(OmegaP，OmegaS，Rp，aS， s )；z0，p0，k0=ellipap(N，Rp，As)；p = P0 *c；z = z0 * OmegaCk=k0*omegacn；ba0

28、= real(poly(z0)；ba0 = k0 * ba0aa0 =实数(poly(P0)；ba = real(poly(z)；ba = k * baaa = real(poly(p)；bd，ad=双线性(ba，aa，Fs)；bd1，ad1=双线性(ba0，aa0，Fs/omega c)；子情节(1，1，1)H，w=freqz(bd，ad，1000，整数)；H =(H(1:1:501)；w =(w(1:1:501)；mag = ABS(H)；db = 20 * log 10(mag+EPS)/max(mag)；pha =角度(H)；grd=grpdelay(bd，ad，w)；支线剧情(2，2

29、，1)；绘图(w/pi，mag)；标题(“A”)；xlabel(“”)；y label( | A | )；set(gca， XTickMode ， manual ， XTick ，0，0.4，0.6，1)；网格打开；set(gca， XTickMode ， manual ， YTick ，0，Attn，Ripple，1)；支线剧情(2，2，3)；绘图(w/pi，db)；title( | A |(dB)；xlabel( freq(rad/s)；y label( dB )；网格打开；支线剧情(2，2，2)；plot(w/pi，PHA/pi)；标题(“阶段”)；xlabel(“”)；y label(

30、 pi )；网格打开；支线剧情(2，2，4)；绘图(w/pi，grd)；标题(下岗)；xlabel( freq(rad/s)；伊拉贝尔(尼翁)；网格打开；用脉冲响应不变量法和双线性变换法设计IIR数字滤波器的优缺点和适用范围脉冲响应不变方法；优点:频率坐标变换是线性的。如果模拟滤波器的频率响应限于折叠频率，那么变换后得到的数字滤波器的频率响应可以不失真地反映源响应与频率的关系。缺点:频谱的周期性延伸效应，所以当模拟滤波器的频率响应不局限于折叠频率时，必然会产生频谱混叠失真。用法:当模拟滤波器的频率响应在折叠频率以上衰减较大时，用脉冲响应不变性方法设计的滤波器可以满足要求。双线性变换法:优点:消

31、除了脉冲响应不变方法引起的频谱混叠失真。缺点:频率之间存在严重的非线性。用法:变换后的连续时间系统的响应必须是分段常数，否则数字域的频率和连续时间域的频谱之间会有很大的失真。5.主要结论、遇到的问题和解决办法、收获和经验。主要结论:同一个系统用巴特沃斯滤波器，需要的阶数最少，也是最好看的。切比雪夫型在通带会有选择滤波，切比雪夫型在阻带会有一定的滤波效果，而椭圆滤波器的通带和阻带都会有一定的选择滤波效果。遇到的问题及解决方法:刚开始做切比雪夫type的时候，不知道他的功能是什么，怎么用。后来在实验教程里看了巴特沃斯的滤镜设计，慢慢就学会了。并获得经验:MATLAB功能强大，我用MATLAB学会了

32、几种分析方法。用4窗函数法设计FIR数字滤波器实验目的(1)掌握用窗函数法设计FIR数字滤波器的原理和方法。(2)熟悉线性相位FIR数字滤波器的特性。(3)了解各种窗函数对滤波特性的影响。二。实验设备和环境计算机和MATLAB软件环境。三、实验原理滤波器的理想频率响应函数为Hd(ej)，其对应的单位脉冲响应为窗口设计方法的基本原理是用有限单位脉冲响应序列h(n)强制。由于hd(n)通常是一个无限序列，并且是非因果的，所以使用窗口函数。W(n)被截断和加权:H(n)作为实际设计的FIR数字滤波器的单位冲激响应序列，其频率响应函数为用窗函数法设计的滤波器的性能取决于窗函数w(n)的类型和窗长N的取

33、值，在设计过程中应根据阻带衰减最小和过渡带宽度的要求选择合适的窗函数类型和窗长N。用窗函数法设计的滤波器的性能取决于窗函数(n)的类型和窗长N的取值，在设计过程中应根据阻带衰减最小和过渡带宽度的要求选择合适的窗函数类型和窗长N。各种类型的窗函数所能达到的阻带最小衰减和过渡带宽度可以在教科书中找到。选定窗函数类型和长度n后，得到单位冲激响应，从而得到。检查是否符合要求。一般在h(n)的尾部加零，使长度满足2的整数次幂，这样可以通过FFT计算。如果想观察细节，可以增加补零的次数。如果不符合要求，就要重新选择窗函数类型和长度N，再次检查，直到符合要求。如果需要线性相位特性，那么h(n)也必须满足:。

34、根据公式中的符号和长度n的奇偶性，线性相位FIR滤波器分为四类。应根据设计的滤波特性正确选择其中之一。四、实验公差和步骤1.根据以下技术规格设计一个数字低通FIR滤波器:矩形窗、汉宁窗、汉明窗、布莱克曼窗和凯泽窗被用来设计滤波器。结合实验结果，讨论了用上述方法设计的数字滤波器能否满足给定指标的要求。矩形窗口:Wp = 0.2 * piWst = 0.3 * piRp = 0.25As = 50tr _ width = Wst-Wp；n = ceil(1.8 * pi/tr _ width)+1；N = 0:(N-1)；WC =(Wp+Wst)/2；alpha =(N-1)/2；HD =(Wc/

35、pi)* sinc(Wc/pi)*(n-alpha)；w _ boxcar = boxcar(N)；h =高清。* w _ boxcar支线剧情(221)；词干(n，hd，填充)；轴紧；xlabel( n )；ylabel( HD(n)；Hr，w1=零相位(h)；支线剧情(222)；plot(w1/pi，Hr)；轴紧；xlabel( omega/ pi )；ylabel( H( omega)；支线剧情(223)；stem(n，h，填充)；轴紧；xlabel( n )；ylabel( h(n)；H，w=freqz(h，1)；支线剧情(224)；plot(w/pi，20 * log 10(ABS(

36、H)/max(ABS(H)；xlabel( omega/ pi )；y label( dB )；网格打开；从上图可以看出，后面的面积还是大于50dB，不符合设计要求。汉宁窗:Wp = 0.2 * piWst = 0.3 * piRp = 0.25As = 50tr _ width = Wst-Wp；n = ceil(6.2 * pi/tr _ width)+1；N = 0:(N-1)；WC =(Wp+Wst)/2；alpha =(N-1)/2；HD =(Wc/pi)* sinc(Wc/pi)*(n-alpha)；w _ hanning = hanning(N)；h =高清。* w _汉宁；支线

37、剧情(221)；词干(n，hd，填充)；轴紧；xlabel( n )；ylabel( HD(n)；Hr，w1=零相位(h)；支线剧情(222)；plot(w1/pi，Hr)；轴紧；xlabel( omega/ pi )；ylabel( H( omega)；支线剧情(223)；stem(n，h，填充)；轴紧；xlabel( n )；ylabel( h(n)；H，w=freqz(h，1)；支线剧情(224)；plot(w/pi，20 * log 10(ABS(H)/max(ABS(H)；xlabel( omega/ pi )；y label( dB )；网格打开；满足设计要求。海明创:Wp = 0

38、.2 * piWst = 0.3 * piRp = 0.25As = 50tr _ width = Wst-Wp；n = ceil(6.6 * pi/tr _ width)+1；N = 0:(N-1)；WC =(Wp+Wst)/2；alpha =(N-1)/2；HD =(Wc/pi)* sinc(Wc/pi)*(n-alpha)；w _ hamming = hamming(N)；h =高清。* w _海明；支线剧情(221)；词干(n，hd，填充)；轴紧；xlabel( n )；ylabel( HD(n)；Hr，w1=零相位(h)；支线剧情(222)；plot(w1/pi，Hr)；轴紧；xla

39、bel( omega/ pi )；ylabel( H( omega)；支线剧情(223)；stem(n，h，填充)；轴紧；xlabel( n )；ylabel( h(n)；H，w=freqz(h，1)；支线剧情(224)；plot(w/pi，20 * log 10(ABS(H)/max(ABS(H)；xlabel( omega/ pi )；y label( dB )；网格打开；满足设计要求。布莱克曼窗口:Wp = 0.2 * piWst = 0.3 * piRp = 0.25As = 50tr _ width = Wst-Wp；n = ceil(11 * pi/tr _ width)+1；N

40、= 0:(N-1)；WC =(Wp+Wst)/2；alpha =(N-1)/2；HD =(Wc/pi)* sinc(Wc/pi)*(n-alpha)；w _ blackman = blackman(N)；h =高清。* w _布莱克曼；支线剧情(221)；词干(n，hd，填充)；轴紧；xlabel( n )；ylabel( HD(n)；Hr，w1=零相位(h)；支线剧情(222)；plot(w1/pi，Hr)；轴紧；xlabel( omega/ pi )；ylabel( H( omega)；支线剧情(223)；stem(n，h，填充)；轴紧；xlabel( n )；ylabel( h(n)；H

41、，w=freqz(h，1)；支线剧情(224)；plot(w/pi，20 * log 10(ABS(H)/max(ABS(H)；xlabel( omega/ pi )；y label( dB )；网格打开；满足设计要求。凯瑟窗:Wp = 0.2 * piWst = 0.3 * piRp = 0.25As = 50tr _ width = Wst-Wp；n = ceil(As-7.95)/(2.285 * tr _ width)+1；N = 0:(N-1)；WC =(Wp+Wst)/2；alpha =(N-1)/2；HD =(Wc/pi)* sinc(Wc/pi)*(n-alpha)；w _ K

42、aiser = Kaiser(N)；h =高清。* w _凯泽；支线剧情(221)；词干(n，hd，填充)；轴紧；xlabel( n )；ylabel( HD(n)；Hr，w1=零相位(h)；支线剧情(222)；plot(w1/pi，Hr)；轴紧；xlabel( omega/ pi )；ylabel( H( omega)；支线剧情(223)；stem(n，h，填充)；轴紧；xlabel( n )；ylabel( h(n)；H，w=freqz(h，1)；支线剧情(224)；plot(w/pi，20 * log 10(ABS(H)/max(ABS(H)；xlabel( omega/ pi )；y

43、label( dB )；网格打开；不符合设计要求。2.按照以下技术规格设计数字带通FIR滤波器下阻带边缘:较低通带边缘:通带上限:上阻带边缘:Wp = 0.15 * piWst = 0.3 * piRP = 1；As = 60tr _ width = Wst-Wp；n = ceil(6.6 * pi/tr _ width)+1；N = 0:(N-1)；WC =(Wp+Wst)/2；alpha =(N-1)/2；HD =(Wc/pi)* sinc(Wc/pi)*(n-alpha)；w _ hamming = hamming(N)；h1 =高清。* w _海明；h=h1。* cos(0.5 * p

44、i * n)；支线剧情(221)；词干(n，hd，填充)；轴紧；xlabel( n )；ylabel( HD(n)；Hr，w1=零相位(h)；支线剧情(222)；plot(w1/pi，Hr)；轴紧；xlabel( omega/ pi )；ylabel( H( omega)；支线剧情(223)；stem(n，h，填充)；轴紧；xlabel( n )；ylabel( h(n)；H，w=freqz(h，1)；支线剧情(224)；plot(w/pi，20 * log 10(ABS(H)/max(ABS(H)；xlabel( omega/ pi )；y label( dB )；网格打开；3.FIR数字低

45、通滤波器采用频率采样设计方法设计，满足以下指标如果N = 20，则过渡区没有样本。源代码:N = 20alpha =(N-1)/2；l = 0:N-1；w1 =(2 * pi/N)* l；Hrs=1，1，1，零(1，15)，1，1；Hdr=1，1，0，0；wd1=0，0.25，0.25，1；k1 = 0:floor(N-1)/2)；k2 =地板(N-1)/2)+1:N-1；angH=-alpha*(2*pi)/N*k1，alpha *(2 * pi)/N *(N-k2)；h =小时。* exp(1j * angH)；h=ifft(H，N)；w =0:500* pi/500；H=freqz(h，

46、1，w)；Hr，wr=零相位(h)；支线剧情(221)；plot(wd1，Hdr，w1(1:11)/pi，Hrs(1:11)，o)；轴(0，1，-0.1，1.1)；xlabel( omega( pi)；ylabel( H(k)；支线剧情(222)；stem(l，h，填充)；轴(0，N-1，-0.1，0.3)；xlabel( n )；ylabel( h(n)；支线剧情(223)；plot(wr/pi，Hr，w1(1:11)/pi，Hrs(1:11)，o)；轴(0，1，-0.2，1.2)；xlabel( omega( pi)；ylabel( Hr(w)；支线剧情(224)；plot(w/pi，20

47、 * log 10(ABS(H)/max(ABS(H)；轴(0，1，-50，5)；网格；xlabel( omega( pi)；y label( dB )；取N = 40，过渡中有一个样本，T = 0.39。源代码:N = 40alpha =(N-1)/2；l = 0:N-1；w1 =(2 * pi/N)* l；Hrs=1，1，1，1，0.39，零(1，29)，0.39，1，1，1，1；Hdr=1，1，0.39，0，0；wd1=0，0.2，0.25，0.3，1；k1 = 0:floor(N-1)/2)；k2 =地板(N-1)/2)+1:N-1；angH=-alpha*(2*pi)/N*k1，al

48、pha *(2 * pi)/N *(N-k2)；h =小时。* exp(1j * angH)；h=ifft(H，N)；w =0:500* pi/500；H=freqz(h，1，w)；Hr，wr=零相位(h)；支线剧情(221)；plot(wd1，Hdr，w1(1:21)/pi，Hrs(1:21)，o)；轴(0，1，-0.1，1.1)；xlabel( omega( pi)；ylabel( H(k)；支线剧情(222)；stem(l，h，填充)；轴(0，N-1，-0.1，0.3)；xlabel( n )；ylabel( h(n)；支线剧情(223)；plot(wr/pi，Hr，w1(1:21)/p

49、i，Hr(1:21)，o)；轴(0，1，-0.2，1.2)；xlabel( omega( pi)；ylabel( Hr(w)；支线剧情(224)；plot(w/pi，20 * log 10(ABS(H)/max(ABS(H)；轴(0，1，-50，5)；网格；xlabel( omega( pi)；y label( dB )；N = 60，有两个样本在转换，T1 = 0.5925，T2 = 0.1099。N = 60alpha =(N-1)/2；l = 0:N-1；w1 =(2 * pi/N)* l；Hrs=1，1，1，1，1，1，0.5925，0.1099，0.1099，0.5925，1，1，1，1，1，1；Hdr=1，1，1，0.59