毕业设计（论文）-数字滤波器的实现.docVIP

班级：测控一班 姓名：胡婷婷 学号：201006010125 指导老师：葛青老师 关键字：滤波器 无限长单位脉冲响应（IIR） 有限长单位脉冲响应（FIR） 设计 论文摘要：，　　 　　 有限冲激响应滤波器的传输函数为 在现代通信系统中，信号中经常混有各种复杂成分，所以很多信号的处理和分析都是基于滤波器而进行的。但是传统的数字滤波器设计过程复杂，滤波特性调整困难，利用MATALB信号处理工具箱（signal processing toolbox）可以快速有效的实现由软件组成的常规数字滤波器的设计、分析和仿真，极大地减轻了工作量，有利于滤波器设计的最优化。 脉冲响应不变法就是要求数字滤波器的脉冲响应序列h（n）与模拟滤波器的脉冲响应ha(t)采样值相等，根据模拟信号的拉普拉斯变换与离散序列z变换之间的关系，我们知道： 将Ha(s)映射为H(Z),即得到数字滤波器的系统函数H(Z): 1.双线性变换法设计IIR数字滤波器 脉冲响应不变法的主要缺点是会产生频谱混叠现象，为了克服这一缺点，可以采用双线性变换法，将整个频率轴上的频率范围压缩到s1平面上的-pi/T~pi/T之间，再用z=e^(s1*T)转换到Z平面上。 这里用正切变换实现频率压缩： 最终得到的转换关系为： 2.IIR数字巴特沃斯滤波器的MATALB实现 2.1 MATALB信号工具箱函数简介 MATALB信号工具箱提供了几个直接设计IIR数字滤波器的函数，直接调用这些函数就可以很方便地对滤波器进行设计。用MATALB工具箱函数涉及数字巴特沃斯滤波器的函数如下： 求最小阶N和3dB截止频率的函数buttord [N,wc]=buttord(wp,ws,Rp,As),此处wp,ws分别为数字滤波器的通带边界和阻带边界频率的归一化值。Rp和As分别指通带最大衰减和阻带最小衰减（dB）。 计算系统函数分子和分母多项式系数向量B和A的函数butter [B,A]=butter(N,wc,’ftype’),计算N阶巴特沃斯数字滤波器系统函数分子和分母多项式系数向量B和A,参数ftype用来区分滤波器的类型。 求模拟滤波器Ha(S)的频率响应函数freqs 脉冲响应不变法的转换函数impinvar和双线性变换法的转换函数bilinear [Bz,Az]=impinvar(B,A,Fs)和双线性变换法将分子和分母多项式的系数向量B和A的模拟滤波器的系统函数Ha(s)转换成数字滤波器的的系统函数H(Z)。 2.2 MATALB程序设计法示例 例1：设采样频率=4KHZ，分别用脉冲响应不变法，和双线性变换法设计一个三阶巴特沃斯低通滤波器，其3dB边界频率为1000Hz。 设计过程如下： [B,A]=butter(3,2*pi*1000,s); [Bz1,Az1]=impinvar(B,A,4000); [h1,w]=freqz(Bz1,Az1); [B,A]=butter(3,2/0.00025,s); [Bz2,Az2]=bilinear(B,A,4000); [h2,w]=freqz(Bz2,Az2); f=w/pi*2000; plot(f,abs(h1),-.,f,abs(h2),-);grid xlabel(频率/Hz);ylabel(幅值/dB) legend(脉冲响应不变法,双线性变换法); 下图给出了这两种设计方法得到的频响，虚线为脉冲响应法得到的结果，实线为双线性变换法的结果。 例2：设计一个巴特沃斯低通滤波器，性能指标如下：通带截止频率6KHz;通带截止频率最大衰减ap=3dB；通带最小衰减af=25dB;阻带截止频率12KHZ。假设输入该滤波器的信号为混合信号，低频信号f=5KHZ,高频信号f=20KHZ,取样间隔T=0.025us,要求从输入信号中滤除高频信号，仅留下低频信号。 设计过程如下： FL=5000,FH=20000;Fs=2*FH;n=0:49; input=cos(2*pi*FL*n/Fs)+cos(2*pi*FH*n/Fs);%绘制输入信号 subplot(2,2,3);plot(n,input);axis([0 49 -2 2]); ylabel(幅度);title(输入信号);%设计巴特沃斯低通滤波器 Wp=6000/FH;Ws=12000/FH; [N,Wn]=buttord(Wp,Ws,3,25); [b,a]=butter(N,Wn);%巴特沃斯滤波器的频响特性 [H,w]=freqz(b,a,FH,Fs); Hr=abs(H);Hphase=angle(H);%求频响 Hphase=unwrap(Hphase);%解卷绕 subplot(2,2,1);plot(w,Hr);grid on; sztitle