基于MATLAB的IIR滤波器的设计和应用（信号去噪).docVIP

word专业整理 学习资料 整理分享 ※※※※※※※※※※※ ※※※※※※※※※ ※※ ※※ ※※ ※※※※※※※※※ 2010级数字信号处理 课程设计 数字信号处理课程设计报告书 课题名称 基于MATLAB的IIR滤波器的设计及应用（信号去噪） 姓 名 学 号 院、系、部 电气工程系 专 业 电子信息工程 指导教师 2013年 6 月28日 学习资料 整理分享 基于MATLAB的IIR滤波器的设计及应用（信号去噪） 一、实验目的 1．学会MATLAB的使用，掌握MATLAB的程序设计方法。 2．掌握数字信号处理的基本概念、基本理论和基本方法。 3. 在MATLAB环境下产生噪声信号。 4．掌握MATLAB设计IIR数字滤波器的方法。 5．学会用MATLAB对信号进行分析和处理。 二、实验原理 数字滤波器的设计： 巴特沃斯(Butterworth)滤波器的幅度平方函数用下式表示： (2.1) 式中，N为滤波器的阶数，幅度下降的速度与N有关，N越大，通带越平坦，过渡带越窄，总的频响特性与理想低通滤波器的误差越小。 切比雪夫(Chebyshev)滤波器的幅频特性在通带或者阻带有等波纹特性，可以提高选择性，其幅度平方函数用下式表示： (2.2) 式中，为小于1的正数，表示通带幅度波动的程度，越大，波动幅度也越大，Ωp称为通带截止频率。 椭圆(Ellipse)滤波器的通带和阻带呈现等波纹幅频特性时，通带和阻带波纹幅度越小，过渡带就越宽。所以椭圆滤波器的阶数由通带边界频率、阻带边界频率、通带边界衰减、阻带边界衰减共同决定。 三、主要实验仪器及材料 微型计算机、MATLAB6.5教学版 四、实验内容 1．噪声信号的频谱分析。 2．设计数字滤波器和画出频率响应： 低通滤波器性能指标，fp=1000Hz，fs=1800 Hz， 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;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(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*fc1*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