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

《数字信号处理》课程研究性学习报告姓名陶威东学号14212129同组成员张欣悦14212138指导教师刘留时间2015年6月1日数字滤波器设计专题研讨【目的】(1)掌握IIR和FIR数字滤波器的设计方法及各自的特点。(2)掌握各种窗函数的时频特性及对滤波器设计的影响。(3)培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。【研讨题目】基本题1．分析矩形窗、汉纳窗、哈明窗、布莱克曼窗、凯泽窗的频域特性，并进行比较。【题目分析】不同窗函数的主瓣宽度，主瓣高度以及过渡带长度和阻带波动都有明显差别。下面我们将通过对时域长度相同的窗函数进行分析来比较不同窗函数的频域特性。取时域长度为50【仿真结果】【结果分析】各种窗有何特点?通带衰减由小到大的排列为：ｂｌａｃｋｍａｎ，ｈａｍｍｉｎｇ，ｋａｉｓｅｒ；阻带衰减由小到大的排列为：ｋａｉｓｅｒ，ｈａｍｍｉｎｇ，ｂｌａｃｋｍａｎ；过渡带由小到大的排列为：ｋａｉｓｅｒ，ｈａｍｍｉｎｇ，ｂｌａｃｋｍａｎ；在设计滤波器的时候，我们希望得到最大的通带衰减，最小的阻带衰减及最小的过渡带。但在实践中，这三者是不可能同时满足的。所以，我们要根据实际需求，选择合适的窗函数，以达到最佳的滤波效果。最后我们还可以看出，主瓣宽度与主瓣幅值的乘积是定值，即主瓣宽度的增加必将导致主瓣幅值的降低。最后我们会发现凯泽窗幅度响应曲线和矩形窗的完全吻合。这是因为凯泽窗是是一个可以变幻形状的窗函数，当系统没有给值时，默认的为矩形窗，所以是完全符合的。【自主学习内容】几种不同窗函数的设计【阅读文献】[1]陈后金.数字信号处理[M].北京：高等教育出版社.2008.11【发现问题】(专题研讨或相关知识点学习中发现的问题)：对于不同窗函数的选择上，我们需要根据需要在主瓣宽度和过渡带长短之间做出选择【问题探究】在谱分析中如何选择窗函数，在滤波器设计中如何选择窗函数？提高FIR滤波器阻带衰减是以增加过渡带宽度为代价的。实际中，在满足阻带衰减的前提下，尽可能地选择主瓣宽度较小的窗函数。凯泽窗是可调的，它通过改变窗函数的形状来控制窗函数旁瓣的大小，设计中可根据滤波器的衰减指标来确定窗函数的形状。【仿真程序】M=50;L=512;fsam=1000;wsam=2*pi*fsam;w=(-wsam/2+(0:L-1)*wsam/L)/fsam;w1=ones(1,M);X1=fftshift(fft(w1,L));plot(w,abs(X1),'g');holdon;xlabel('Normalizedfrequency');legend('square');title('·幅度响应¦');gridon;w2=hanning(M)';X2=fftshift(fft(w2,L));subplot(2,2,1),plot(w,abs(X2),'r');holdon;xlabel('Normalizedfrequency');legend('hanning');title('·幅度响应');gridon;w3=hamming(M)';X3=fftshift(fft(w3,L));subplot(2,2,2),plot(w,abs(X3),'y');holdon;xlabel('Normalizedfrequency');legend('hamming');title('幅度响应¦');gridon;w4=blackman(M)';X4=fftshift(fft(w4,L));subplot(2,2,3),plot(w,abs(X4),'b');holdon;xlabel('Normalizedfrequency');legend('blackman');title('幅度响应');gridon;w5=kaiser(M)';X5=fftshift(fft(w5,L));subplot(2,2,4),plot(w,abs(X5),'--');xlabel('Normalizedfrequency');legend('kaiser');title('幅度响应');gridon;【研讨题目】基本题2．（M5-5）在用窗口法设计FIR滤波器时，由于理想滤波器的频幅响应在截频处发生突变，使得设计出的滤波器的频幅响应发生振荡，这个现象被称为Gibbs现象。解决这个问题的一个方案是本书中介绍的用逐步衰减的窗函数。另一个方案是使理想滤波器过渡带为渐变的，如下图所示具有线性过渡带的理想低通滤波器的频率响应，试用窗口法设计逼近该频率响应的FIR滤波器。pHL(ej)sps题2图【设计步骤】渐变的窗选为hamming窗。取wp=0.55pi,ws=0.45pi,As=25db,Ap=1db。设hamming窗的长度为M,矩形窗的长度为M1。M1取不同的值作图，其中M=7【单位脉冲响应证明】试证该滤波器的单位脉冲响应为0π)sin(Δ)2/Δsin(2,0,π][cckkkkkkkhL其中：psΔ，2/)(spc证明如图【仿真结果】M=8M=16M=32M=64【结果分析】1.通过逐步衰减的窗函数法和渐变过渡带两种方法设计的滤波器都能较好的设计较小Gibbs现象的滤波器。通过理论分析我们可以知道，FIR滤波器的波动是由于窗函数的主瓣面积和旁瓣面积决定的，所以增加采样点并不能改变通阻带的波动。而对于滤波器的渐变法，则可以通过改变点数很好的减少波动的范围。【自主学习内容】逐步减变法的设计【阅读文献】[1]陈后金.数字信号处理[M].北京：高等教育出版社.2008.11【发现问题】(专题研讨或相关知识点学习中发现的问题)：【问题探究】通过实验讨论如何控制滤波器的阻带衰减sA【仿真程序】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(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);holdon;plot(omega2/pi,magdb2,'r');legend('逐步衰减','过渡带渐变')gridon【研讨题目】基本题3．试用频率取样法设计频率响应逼近π,j)e(jDIFH的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)&(Wm=Wp2)]).*Wm;Hd=j*Ad.*exp(-j*0.5*M*Wm);Hd=[Hdconj(fliplr(Hd(2:M/2+1)))];h=real(ifft(Hd));w=linspace(-pi,pi,1000);H=freqz(h,[1],w);plot(w/pi,20*log10(abs(H)));axis([-11-5020]);grid;IV型M=63;Wp1=-pi;Wp2=pi;m=0:M/2;Wm=2*pi*m./(M+1);Ad=double([(Wm=Wp1)&(Wm=Wp2)]).*Wm;Hd=j*Ad.*exp(-j*0.5*M*Wm);Hd=[Hdconj(fliplr(Hd(2:(M+1)/2)))];h=real(ifft(Hd));w=linspace(-pi,pi,1000);H=freqz(h,[1],w);plot(w/pi,20*log10(abs(H)));grid;4．（1）（M5-6）利用频率取样法设计III型线性相位带通FIR滤波器，其通带截频分别为p1=0.3rad,p2=0.5rad（2）（M5-7）在通带和阻带间设置1个过渡点，探讨过渡点幅度的确定方法，研究过渡点幅度对设计结果的影响。【设计步骤】（1）学习利用频率取样法设计某I型线性相位带通FIR滤波器的方法；（2）确定理想滤波器的幅度函数；（3）确定理想滤波器的相位；（4）用实验法确定满足要求的滤波器的阶数；（5）用实验法确定过度点。【仿真结果】衰减大约为-16.87db(2)取过渡点t1=0.29,t2=0.3,此时可得实验最大衰减最大衰减大约为-52.9【结果分析】滤波器的阻带衰减，滤波器的阶数与设计结果的关系。知滤波器阶数越大，阻带衰减越大，但滤波器阶数越大，实际设计难度越大。【自主学习内容】III型FIR滤波器的设计方法【阅读文献】[1]陈后金.数字信号处理[M].北京：高等教育出版社.2008.11【发现问题】(专题研讨或相关知识点学习中发现的问题)：应用实验法寻找过渡点费时费力又不一定得出最好的过渡点【问题探究】如何确定过渡点的最佳幅度。由大及小，不断试验确定【仿真程序】（1）M=64;Wp1=0.3*pi;Wp2=0.5*pi;m=0:M/2;Wm=2*pi*m./(M+1);Ad=double([(Wm=Wp1)&(Wm=Wp2)]);Hd=Ad.*exp(-j*0.5*M*Wm);Hd=[Hdconj(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;（2）M=64;Wp1=0.3*pi;Wp2=0.5*pi;Wp=(Wp1+Wp2)/2;m=0:M/2;Wm=2*pi*m./(M+1);mtr1=ceil(Wp1*(M+1)/(2*pi));mtr2=floor(Wp2*(M+1)/(2*pi))+2;Ad=double([(Wm=Wp1)&(Wm=Wp2)]);Ad(mtr1)=0.29;Ad(mtr2)=0.30;Hd=Ad.*exp(-j*0.5*M*Wm);Hd=[Hdconj(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;【研讨题目】中等题5．Dhexian.wav是对频率为293.66,369.99,440Hz的D大调和弦以8000Hz抽样所得的数字音乐信号，试设计一数字滤波器从和弦中分离出440Hz的音符。要求：（1）设计IIR数字高通滤波器，通过实验研究sP,，sP,AA的选择对滤波效果及滤波器阶数的影响，给出滤波器指标选择的基本原则，确定你认为最合适的滤波器指标。（2）能否用IIR数字带通滤波器从和弦中分离出440Hz的音符？利用（1）确定的基本原则，给出数字带通滤波器的指标。设计IIR数字带通滤波器，并将结果与高通滤波器比较，给出你的结论。（3）用窗函数法设计FIR数字高通滤波器，分别利用矩形窗、汉纳窗、哈明窗、布莱