基于MATLAB的IIR滤波器的设计及应用(信号去噪)

※※※※※※※※※※※※※※※※※※※※※※※※※※※※※※※※※※※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=1800Hz，As=100dB，Ap=1dB；在MATLAB中，可以利用函数butte、cheby1和ellip设计IIR滤波器；最后，利用MATLAB中的函数freqz画出各滤波器的频率响应。3．用滤波器对噪声信号进行滤波。4．比较几种滤波器去滤波前后噪声信号的波形及频谱。图6-1噪声波形及频谱图6-2信号波形及原理图6-3叠加波形及频谱三种滤波器的通带边界频率=100HZ，阻带截止频率=1800HZ时，分别设计的巴特沃斯滤波器频率响应见图6-4，切比雪夫滤波器的频率响应见图6-5，椭圆滤波器的频率响应见图6-6。图6-4巴特沃斯滤波器的频率响应图6-5切比雪夫滤波器的频率响应图6-6椭圆滤波器的频率响应经过巴特沃斯滤波器后的信号波形和频谱见图6-7，经过切比雪夫滤波器后的信号波形和频谱见图6-8，经过椭圆滤波器后的信号波形和频谱见图6-9。图6-7经巴特沃斯滤波器后的波形和频谱图6-8经切比雪夫滤波器后的波形和频谱图6-9经椭圆滤波器后的波形和频谱经过三种滤波器后的波形稍有区别，经过巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器的波形的起点位置不同。频谱的区别不大，都是实现了滤除高频的功能。设计总结模拟滤波器从滤波特性上分类，可以分为低通、高通、带通、带阻等滤波器。这次小学期设计的是模拟低通滤波器，巴特沃斯模拟滤波器具有单调下降的幅频特性；切比雪夫的幅频特性在通带或者阻带有等波纹特性，可以提高选择性；椭圆滤波器的选择性相对前两种是最好的但通带和阻带内均呈现等波纹幅频特性，相位特性的非线性也稍严重。设计滤波器时，先设计低通滤波器再通过频率变换将低通滤波器转换成希望类型的滤波器。这次设计让我学会了MATLAB的使用方法，学会分析滤波器的优劣和性能，提高了分析和动手实践能力。同时我相信，进一步加强对MATLAB的学习与研究对我今后的学习将会起到很大的帮助。八、参考文献[1]丁玉美,高西全.数字信号处