资源描述：
《基于matlab的iir滤波器的设计及应用(信号去噪)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、.※※※※※※※※※※※※※※※※※※※※※※※※2010级数字信号处理课程设计数字信号处理课程设计报告书课题名称基于MATLAB的IIR滤波器的设计及应用（信号去噪）姓名学号院、系、部电气工程系专业电子信息工程指导教师2013年6月28日...基于MATLAB的IIR滤波器的设计及应用（信号去噪）一、实验目的1．学会MATLAB的使用，掌握MATLAB的程序设计方法。2．掌握数字信号处理的基本概念、基本理论和基本方法。3.在MATLAB环境下产生噪声信号。4．掌握MATLAB设计IIR数字滤波器的方法。5．学会用MATLAB对信号进行分析和处理。二、实验原理数字滤波器的设计：巴特沃斯
2、(Butterworth)滤波器的幅度平方函数用下式表示：(2.1)式中，N为滤波器的阶数，幅度下降的速度与N有关，N越大，通带越平坦，过渡带越窄，总的频响特性与理想低通滤波器的误差越小。切比雪夫(Chebyshev)滤波器的幅频特性在通带或者阻带有等波纹特性，可以提高选择性，其幅度平方函数用下式表示：(2.2)式中，为小于1的正数，表示通带幅度波动的程度，越大，波动幅度也越大，Ωp称为通带截止频率。椭圆(Ellipse)滤波器的通带和阻带呈现等波纹幅频特性时，通带和阻带波纹幅度越小，过渡带就越宽。所以椭圆滤波器的阶数由通带边界频率、阻带边界频率、通带边界衰减、阻带边界衰减共同决定。三
3、、主要实验仪器及材料微型计算机、MATLAB6.5教学版四、实验内容1．噪声信号的频谱分析。2．设计数字滤波器和画出频率响应：低通滤波器性能指标，fp=1000Hz，fs=1800Hz，As=100dB，Ap=1dB；在MATLAB中，可以利用函数butte、cheby1和ellip设计IIR滤波器；最后，利用MATLAB中的函数freqz画出各滤波器的频率响应。3．用滤波器对噪声信号进行滤波。4．比较几种滤波器去滤波前后噪声信号的波形及频谱。...五、程序设计1、噪声与信号的产生分析：N=1000;Fs=5000;T=1/Fs;Tp=N*T;t=-(N-1)*T:T:(N-1)*T;
4、k=-N-1:N-1;f=k/Tp;fm1=100;fc1=500;nt=cos(2*pi*fc1*t);subplot(3,3,1);plot(t,nt);axis([-Tp/8,Tp/8,-2,2]);title('噪声波形');xlabel('t(s)');ylabel('n(t)')fst=fft(nt,N);fst=fftshift(fst);k=-N/2:N/2-1;f=k/Tp;abs(fst);subplot(3,3,2);stem(f,abs(fst),'.');grid;axis([-Fs/2,Fs/2,0,1.2]);title('噪声频谱');xlabel('f
5、(Hz)');ylabel('幅度(dB)')yt=cos(2*pi*fm1*t);subplot(3,3,3);plot(t,yt);xlabel('t(s)');ylabel('y(t)');axis([-Tp/4,Tp/4,-2,2]);title('信号波形');fmt=fft(yt,N);fmt=fftshift(fmt);subplot(3,3,4);plot(f,fmt);axis([-Fs/2,Fs/2,0,1.2]);title('信号频谱');xlabel('f(Hz)');ylabel('幅度(dB)')xt=cos(2*pi*fm1*t)+cos(2*pi*fc
6、1*t);subplot(3,3,5);plot(t,xt);xlabel('t(s)');ylabel('x(t)');axis([-Tp/8,Tp/8,-2,2]);title('叠加波形');fxt=fft(xt,N);fxt=fftshift(fxt);subplot(3,3,6);plot(f,fxt);axis([-Fs/2,Fs/2,0,1.2]);title('叠加频谱');xlabel('f(Hz)');ylabel('幅度(dB)')2、巴特沃斯低通滤波器：fp=100;fs=1800;Ap=1;As=100;...wp=2*pi*fp;ws=2*pi*fs;T=8
7、000;Fs=1/T;Wp=wp/T;Ws=ws/T;[N,wc]=buttord(Wp,Ws,Ap,As,'s');[B,A]=butter(N,wc,'s');[Bz,Az]=impinvar(B,A);[h,w]=freqz(Bz,Az);subplot(3,3,7);plot(w/pi,20*log10(abs(h)));grid;xlabel('频率(KHZ)');ylabel('幅度(dB)');title('滤波器频率响应');z
