《数字滤波器研究报告课程论文》由会员分享，可在线阅读，更多相关《数字滤波器研究报告课程论文（16页珍藏版）》请在金锄头文库上搜索。
1、- .研 究 生 课 程 论 文(2010-2011学年第一学期)数字滤波器的研究研究生：提交日期： 研究生签名：学 号学 院电子与信息学院课程编号S0809074课程名称现代数字信号处理学位类别硕士任课教师教师评语：成绩评定： 分 任课教师签名： 年 月 日. word.zl- .中文题目：数字滤波器的研究研究生摘要：无限冲激响应IIR数字滤波器和有限冲激响应FIR数字滤波器广泛应用于数字信号处理系统中。FIR 滤波器具有很好的线性相位特性，因此线性相位特性在工程实际中如数据通信、图像处理等应用领域具有比拟重要的意义。IIR 数字滤波器具构造简单和实现容易的特点, 这对于相位要求不是特别严格
2、的场合有很重要的意义。利用MATLAB 实现IIR 低通数字滤波器的设计,并对被干扰的信号进展滤波,到达了预期结果。关键词: IIR；FIR；线性相位；数字滤波器；MATLAB1引言数字滤波器是现代数字信号处理的不可或缺的一局部，数字滤波器按其响应形式可分为无限脉冲响应IIR滤波器和有限脉冲响应FIR滤波器和两种。当系统无严格相位要求时，IIR 数字滤波器与FIR 数字滤波器相比，可用较低的阶数获得较高的选择性1。对于一样的设计指标，FIR 滤波器所要求的阶数比IIR滤波器高510 倍，而且信号的延迟也较大2。MATLAB 作为一种矩阵运算为根底的交互式程序语言，是进展科学研究常用且必不可少的
3、工具。MATLAB着重针对科学计算、工程计算和绘图的需求。它用解释方式工作，键入程序立即得出结果，人机交互性能好，为科学人员所乐于承受。MATLAB 提供了丰富的函数，其中fir1 函数实现了加窗线性相位FIR 滤波器设计的经典方法，fir1 主要用于常用的标准通带滤波器设计，包括低通、带通、高通和带阻数字滤波器3。一般来讲，数字滤波器的设计要经过三步：确定设计指标、模拟逼近和数字转换。通常在设计滤波器之前，先根据具体的应用确定一些技术指标，然后就可以根据数学知识和滤波器的根本原理提出一个滤波器的模型来逼近给定的指标，逼近的结果通常是得到以差分方程或脉冲响应描述的滤波器，最后可以根据这个描述用
4、硬件或软件实现，至此完成一个滤波器设计的全过程3。2.数字滤波器的设计2.1模拟椭圆滤波器对加噪语音信号的滤波：设计模拟滤波器的一般方法是：首先根据技术指标确定滤波器的传输函数H(s,然后综合电路网络实现该传递函数。而设计滤波器传输函数的关键是找到逼近函数，在各种滤波器响应中，椭圆函数响应最优越4-6。首先，翻开的是一个摩托车引擎信号的波形文件，它是一个频率围大致分布在5000-6000Hz的一个低频语音信号。对于该语音信号，施加高斯白噪声，设计一个1dB截止频率6000Hz 最小阻带衰减为100dB的椭圆低通滤波器。其原始信号波形、加噪声后信号波形与滤波后信号波形分别图1和图2：图1.加噪声
5、后语音信号波形与频谱图2.原始语音信号与滤波后语音信号波形与频谱比照 通过matlab仿真可以发现，原始语音信号根本可以从高斯白噪声中恢复，但是滤波后的语音信号比原始语音信号有明显的波形延迟，这是模拟椭圆低通滤波器缺乏的地方。2.2 IIR滤波器的构造与设计这里采用的是直接II型构造实现IIR 滤波器，如式1所示： 1此构造便于准确地实现滤波器的零点、极点，也便于调整滤波器的频率响应性能7；另一个优点是所需的存储单元较少，在硬件实现时甚至还可以用一个二阶节进展时分复用，进一步降低了对现场可编程逻辑阵列硬件资源的要求8。IIR滤波器的系数在做强噪音下的语音增强，需要通过共振峰对语音进展端点识别。
6、前两个共振峰对区别不同语音有非常重要的作用9；1 kHz以下根本为第一共振峰的围,对语音感知、语意识别作用比拟重要的第二共振峰根本在1 kHz之外9。因此对带噪语音信号参加预处理环节,即进展数字滤波(通带下限1 kHz,Butterworth滤波器高通,阻带衰减3 dB，语音信号采样频率为8000 Hz)。利用MATLAB提供的FDATool便可直接得到直接II型各个子系由于硬件当中对小数进展运算消耗资源比拟大，需要在计算精度与速度之间优化选择10。下面用双线性变换设计一个0.4dB截止频率为10KHz且在30KHz处有最小阻带衰减为50dB的数字巴特沃兹低通滤波器，其抽样率为100KHz；设
7、计的具体步骤：1：先设计一个符合以上指标的模拟巴特沃兹低通滤波器，作为与要设计的数字滤波器的参照滤波器。先用手算计算出滤波器所需的阶数是Ns=7, 与Ns,wn=buttord(wp,ws,Rp,Rs,s);得出的结果一样。2： 双线性变换公式： 2其中，是模拟的频率，是数字的归一化频率，T是积分步长。用公式2预畸所求数字滤波器的的数字频率指标数字滤波器通带截止频率0.2，阻带截止频率0.6，得到一个等价的模拟低通滤波器的频率指标。3：得到原型模拟低通滤波器的指标后作频率响应，从而得模拟的传输函数H(s)，然后对传输函数进展双线性变换，从而得到数字滤波器的传输函数H(z).结果分析：图3.模拟
8、原型与数字滤波器幅度相应比拟图4. 模拟原型与数字滤波器相位相应比拟由图1可以看出模拟原型低通滤波器的通带截止频率为10KHz对应着数字低通滤波器的截止频率0.2 rad与双线性变换结果一致。由图2可以看出，相位响应就有很大差异了，这是由双线性变换的频率畸变所致。2.3 FIR滤波器的构造与设计对于IIR滤波器必须明确要求推导出来的传输函数是稳定的，而另一方面，对于FIR滤波器的设计，由于滤波器的传输函数是以1 /z的多项式表示的，所以FIR数字滤波器都是稳定的11。利用傅里叶加窗级数法设计线性相位低通FIR滤波器没有使用到fir1函数设计指标：通带截止频率4rad/s，阻带截止频率6rad/
9、s，最大通带衰减为0.2dB,最小阻带衰减为42dB，抽样率为18rad/s.图5.由截短得到的理想FIR低通滤波器的冲击响应图6.各种窗函数设计的低通滤波器结果分析：Hamming窗函数的长度是30通带波纹是Rp = 0.0394dB最小阻带衰减时As = 52dBHanning窗函数的长度是29通带波纹是Rp= 0.0711dB最小阻带衰减时As = 44dBBlackman窗函数的长度是52通带波纹是Rp= 0.0027dB最小阻带衰减时As = 75dB均满足滤波器的性能指标。结论在电子系统中，具有严格的线性相位特性的FIR 数字滤波器被广泛使用，是信号处理的根本构成之一。利用matl
10、ab的强大运算功能，基于matlab信号处理工具箱signal processing toolbox的数字滤波器设计法可以快速有效的设计由软件组成的常规数字滤波器，设计方便、快捷，极大的减轻了工作量。在设计过程中可以比照滤波器特性，随时更改程序参数，以到达滤波器设计的最优化。利用matlab设计数字滤波器在信号处理软件和微机保护中，有着广泛的应用前景。本文通过调用simulink中的功能模块构成数字滤波器的仿真框图，编写matlab语言进展仿真，从结果可以看出各种滤波器的性能并反映实际的情况。参考文献1 惠琼，孟令军，王宏涛.基于Nios 的IIR数字滤波器的设计J.通信技术,2009,42(
11、08)2 王秀敏,汪毓铎,洋,世华. 通信系统中FIR数字滤波器的设计研究J. 通信技术, 2009,(09) . 3梁辰 基于MATLAB的FIR数字滤波器的设计J.机械设计与制造，2010，124 史燕,小雪. 基于Matlab和Multisim的综合性实验椭圆滤波器设计与仿真J. 北华航天工业学院学报, 2008,(S1) . 5 王靖,永全. 数字椭圆滤波器的Matlab设计与实现J. 现代电子技术, 2007,(06) . 6 抒珍,童子权,任丽军,小红. DDS波形合成技术中低通椭圆滤波器的设计J. 理工大学学报, 2004,(05) .7 晓琳. 循环码理论及其译码算法研究设计距
12、离为n 的二元BCH 码构造及其B-M 迭代译码算法实现D.：理工大学,2008.8 王宇. BCH 编码在GPS 探空仪中的应用J.信息平安与通信,2010，(07)9 亚梅, 杜红棉, 志杰. 基于MATLAB一种IIR数字带通滤波器的设计与仿真J. 微计算机信息, 2007,(13)10世华,王秀敏,豪威;基于DSP Builder和FPGA的IIR滤波器设计J,2010,(12)11LAWRENCE.R.R FIR Digital Filter Design Techniques UsingWeighted Chebyshev ApproximationPROCEEDINGS OF T
13、HE IEEE, VOL. 63, NO. 4, APRIL 1975.附仿真源程序：程序1：y1,Fs,bits=wavread(engine.wav,2048);%Fs是音频采样率，bits是每个采样点包含的bit数目sound(y1,Fs,bits);hold onY1=fft(y1,4096);f=Fs*(0:4095)/4096;y2=y1+0.2*randn(2048,1);figure(1)subplot(211)plot(y2);xlabel (时间序号n); ylabel (振幅);legend(叠加噪声后信号波形,1)axis (0 ,2100, -1 , 1);Y2=ff
14、t(y2,4096);subplot(212)plot(f,fftshift(abs(Y2);xlabel (信号频率Hz); ylabel (振幅);legend(叠加噪声后信号频谱,1)axis (4000,7000, 0 , 150);Ap=1;%通带波纹As=100;%最小阻带衰减fp=200;fs=300;Ft=20000;wp=2*pi*fp/Ft;%归一化通带截止角频率ws=2*pi*fs/Ft;%归一化阻带截止角频率n,wn=ellipord(wp,ws,Ap,As);num,den=ellip(n,Ap,As,wn);y3=filter(num,den,y2);Y3=fft(y3,4096);figure(2)subplot(211)plot(y1,r-);hold on
