北交大 数字信号研讨实验2

1、数字信号处理课程研究性学习报告姓名 陶威东 学号 14212129 同组成员 张欣悦 14212138 指导教师 刘留 时间 2015 年6月1日 数字滤波器设计专题研讨【目的】(1) 掌握IIR和FIR数字滤波器的设计方法及各自的特点。(2) 掌握各种窗函数的时频特性及对滤波器设计的影响。(3) 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。【研讨题目】 基本题 1分析矩形窗、汉纳窗、哈明窗、布莱克曼窗、凯泽窗的频域特性，并进行比较。【题目分析】不同窗函数的主瓣宽度，主瓣高度以及过渡带长度和阻带波动都有明显差别。下面我们将通过对时域长度相同的窗函数进行分析来比较不同窗函数的频

2、域特性。取时域长度为50【仿真结果】【结果分析】各种窗有何特点?通带衰减由小到大的排列为：，；阻带衰减由小到大的排列为：，；过渡带由小到大的排列为：，；在设计滤波器的时候，我们希望得到最大的通带衰减，最小的阻带衰减及最小的过渡带。但在实践中，这三者是不可能同时满足的。所以，我们要根据实际需求，选择合适的窗函数，以达到最佳的滤波效果。最后我们还可以看出，主瓣宽度与主瓣幅值的乘积是定值，即主瓣宽度的增加必将导致主瓣幅值的降低。最后我们会发现凯泽窗幅度响应曲线和矩形窗的完全吻合。这是因为凯泽窗是是一个可以变幻形状的窗函数，当系统没有给值时，默认的为矩形窗，所以是完全符合的。【自主学习内容】几种不同窗

3、函数的设计【阅读文献】1 陈后金.数字信号处理M.北京：高等教育出版社.2008.11【发现问题】 (专题研讨或相关知识点学习中发现的问题)：对于不同窗函数的选择上，我们需要根据需要在主瓣宽度和过渡带长短之间做出选择【问题探究】在谱分析中如何选择窗函数，在滤波器设计中如何选择窗函数？提高FIR滤波器阻带衰减是以增加过渡带宽度为代价的。实际中，在满足阻带衰减的前提下，尽可能地选择主瓣宽度较小的窗函数。凯泽窗是可调的，它通过改变窗函数的形状来控制窗函数旁瓣的大小，设计中可根据滤波器的衰减指标来确定窗函数的形状。【仿真程序】M=50;L=512; fsam=1000;wsam=2*pi*fsam;

4、w=(-wsam/2+(0:L-1)*wsam/L)/fsam;w1=ones(1,M);X1=fftshift(fft(w1,L);plot(w,abs(X1),g); hold on;xlabel(Normalized frequency); legend(square); title(幅度响应);grid on;w2=hanning(M); X2=fftshift(fft(w2,L); subplot(2,2,1),plot(w,abs(X2),r); hold on;xlabel(Normalized frequency); legend(hanning); title(幅度响应);g

5、rid on;w3=hamming(M); X3=fftshift(fft(w3,L); subplot(2,2,2),plot(w,abs(X3),y); hold on;xlabel(Normalized frequency); legend(hamming); title(幅度响应);grid on;w4=blackman(M); X4=fftshift(fft(w4,L); subplot(2,2,3),plot(w,abs(X4),b); hold on;xlabel(Normalized frequency); legend(blackman); title(幅度响应);grid

6、on;w5=kaiser(M); X5=fftshift(fft(w5,L); subplot(2,2,4),plot(w,abs(X5),-);xlabel(Normalized frequency); legend(kaiser); title(幅度响应);grid on;【研讨题目】 基本题 2（M5-5）在用窗口法设计FIR滤波器时，由于理想滤波器的频幅响应在截频处发生突变，使得设计出的滤波器的频幅响应发生振荡，这个现象被称为Gibbs现象。解决这个问题的一个方案是本书中介绍的用逐步衰减的窗函数。另一个方案是使理想滤波器过渡带为渐变的，如下图所示具有线性过渡带的理想低通滤波器的频率响应

7、，试用窗口法设计逼近该频率响应的FIR滤波器。 题2图【设计步骤】渐变的窗选为hamming窗。取wp=0.55pi,ws=0.45pi,As=25db,Ap=1db。设hamming窗的长度为M,矩形窗的长度为M1。M1取不同的值作图，其中M=7【单位脉冲响应证明】试证该滤波器的单位脉冲响应为其中：，证明如图【仿真结果】M=8M=16M=32M=64【结果分析】1. 通过逐步衰减的窗函数法和渐变过渡带两种方法设计的滤波器都能较好的设计较小Gibbs现象的滤波器。通过理论分析我们可以知道，FIR滤波器的波动是由于窗函数的主瓣面积和旁瓣面积决定的，所以增加采样点并不能改变通阻带的波动。而对于滤波

8、器的渐变法，则可以通过改变点数很好的减少波动的范围。【自主学习内容】逐步减变法的设计【阅读文献】1 陈后金.数字信号处理M.北京：高等教育出版社.2008.11【发现问题】 (专题研讨或相关知识点学习中发现的问题)：【问题探究】通过实验讨论如何控制滤波器的阻带衰减【仿真程序】wp=0.55*pi;ws=0.45*pi;Ap=1;As=25;N=ceil(7*pi/(wp-ws);N=mod(N+1,2)+N;M=N-1;w=hamming(N);wc=(wp+ws)/2;k=0:M;hd=(wc/pi)*sinc(wc*(k-0.5*M)/pi);h=hd.*w;omega=linspace(

9、0,pi,512);mag=freqz(h,1,omega);magdb=abs(mag);plot(omega/pi,magdb,b);grid;w=ws-wp;M1=8;k2=-M1:M1;wc=(wp+ws)/2;hd=sinc(w*k2/2).*(sin(wc*k2)./(k2.*pi);hd(M1+1)=wc/pi;omega2=linspace(0,pi,512);mag2=freqz(hd,1,omega2);magdb2=abs(mag2);hold on;plot(omega2/pi,magdb2,r);legend(逐步衰减,过渡带渐变 )grid on 【研讨题目】 基本

10、题 3试用频率取样法设计频率响应逼近的FIR数字微分器。分别采用III型IV型线性相位滤波器，并将所得结果进行比较。【题目分析】可利用III型滤波器的频率响应函数和频率取样法设计FIR滤波器【设计步骤】【仿真结果】III型IV型【自主学习内容】【阅读文献】1 陈后金.数字信号处理M.北京：高等教育出版社.2008.11【发现问题】 (专题研讨或相关知识点学习中发现的问题)：在0点处，所设计的滤波器效果并不好，有一定的衰减【问题探究】未发现好的解决办法【仿真程序】III型M=64;Wp1=-pi;Wp2=pi;m=0:M/2;Wm=2*pi*m./(M+1);Ad=double(Wm=Wp1)&

11、(Wm=Wp1)&(Wm=Wp1)&(Wm=Wp1)&(Wm=Wp1)&(Wm=Wp2);Ad(mtr1)=0.29;Ad(mtr2)=0.30;Hd=Ad.*exp(-j*0.5*M*Wm);Hd=Hd conj(fliplr(Hd(2:M/2+1);h=real(ifft(Hd);w=linspace(0.1,pi,1000);H=freqz(h,1,w);plot(w/pi,20*log10(abs(H);grid;（3）矩形窗clear allFsamp = 10e3;fp = 4.2e3;fs = 3.8e3;Ap = 1;As = 50;Wp = 2*pi*fp/Fsamp; %W

12、s = 2*pi*fs/Fsamp;N=ceil(1.8*pi/(Wp-Ws);N=mod(N+1,2)+N;M=N-1;w=ones(1,N);fprintf(N=%.0fn,N);Wc=(Wp+Ws)/2;k=0:M;hd=-(Wc/pi)*sinc(Wc*(k-0.5*M)/pi);hd(0.5*M+1)=hd(0.5*M+1)+1;h=hd.*w;omega=linspace(0,pi,512);mag=freqz(h,1,omega);magdb=20*log10(abs(mag);plot(omega/pi,magdb); grid;xlabel(Normalized freque

13、ncy);ylabel(Gain,dB);hann窗wp=0.84*pi;ws=0.76*pi;Ap=1;As=50;N=ceil(6.4*pi/(wp-ws);N=mod(N+1,2)+N;M=N-1;w=hann(N);wc=(wp+ws)/2;k=0:M;hd=-(Wc/pi)*sinc(Wc*(k-0.5*M)/pi);hd(0.5*M+1)=hd(0.5*M+1)+1;h=hd.*w;omega=linspace(0,pi,512);mag=freqz(h,1,omega);magdb=20*log10(abs(mag);plot(omega/pi,magdb,b);grid;w=w

14、s-wp;N对于其他窗函数只需改变窗函数写法即可得到（4）Fp=0.84;Fs=0.76;ds=0.001;dp=ds; f=Fs Fp;a=0 1;dev=ds dp;M,fo,ao,w = remezord(f,a,dev);h = remez(M,fo,ao,w);w=linspace(0,pi,1000);mag=freqz(h,1,w);plot(w/pi,20*log10(abs(mag);xlabel(Normalized frequency);ylabel(Gain, dB);grid;相位谱程序Wp=0.84*pi; Ws=0.76*pi; Ap=1; As=50;T=2;F

15、s=1/T; wp=2*tan(Wp/2)/T;ws=2*tan(Ws/2)/T;wp1=1/wp;ws1=1/ws;N,wc=buttord(wp1,ws1,Ap,As,s);num,den=butter(N,wc,s);numa,dena=lp2hp(num,den,1);numd,dend=bilinear(numa,dena,Fs);w=linspace(0,pi,1024);h=freqz(numd,dend,w); plot(w/pi, angle(h); grid;xlabel(Normalized frequency);ylabel(Gain,dB);title( BW);Wp=0.84*pi;Ws=0.76*pi;Ap=1;As=50;N=ceil(1.8*pi/(Wp-Ws);N=mod(N+1,2)+N;M=N-1;