数字低通滤波器的设计.docx

《数字低通滤波器的设计.docx》由会员分享，可在线阅读，更多相关《数字低通滤波器的设计.docx（7页珍藏版）》请在第一文库网上搜索。

1、数字低通滤波器的设计摘要：日常生活中数字低通滤波器使用非常广泛。根据技术指标不同，低通滤波器的设计方法和思路也是不同的。利用MAT1AB强大的运算能力和仿真功能，可以用编程的方法或MAT1AB提供的FDAtOoI和Sin1U1ink实现各种低通滤波器的设计和仿真。关键词：数字低通滤波器；MAT1AB；无限长冲激响应滤波器；有限长冲激响应滤波器Abstract：inthedai1y1ife,thedigita11ow-passfi1teriswide1yused.Basedonthedifferenttechnica1targets,thedesignmethodsandthoughtsofth

2、edigita11ow-passfi1teraredifferenttoo.Foundonthepowerfu1ca1cu1atingabi1ityandemu1ationfaci1ityoftheMat1ab,wecanusetheprogrammab1ewaysortheFDAtooi,simu1inkprovidedbytheMat1abtoachievethekindsofdesignandsimu1ationofthe1ow-passfi1ter.Keywords:thedigita11ow-passfi1ter;MAT1AB;IIR;FIR（一）实验目的：掌握R数字低通滤波器的设计

3、方法。（二）实验原理：1、滤波器的分类滤波器分两大类：经典滤波器和现代滤波器。经典滤波器是假定输入信号（）中的有用成分和希望取出的成分各自占有不同的频带。这样，当无（）通过一个线性系统（即滤波器）后可讲欲去除的成分有效的去除。现代滤波器理论研究的主要内容是从含有噪声的数据记录（又称时间序列）中估计出信号的某些特征或信号本身。经典滤波器分为低通、高通、带通、带阻滤波器。每一种又有模拟滤波器（AF）和数字滤波器（DF）。对数字滤波器，又有R滤波器和FIR滤波器。IIRDF的转移函数是：MSbrZTH(Z)=X(2)k1+2JQFk=1FIRDF的转移函数是:N-IH（Z）=Eh（n）Z-n=0FI

4、R滤波器可以对给定的频率特性直接进行设计，而IIR滤波器目前最通用的方法是利用已经很成熟的模拟滤波器的设计方法进行设计。2、滤波器的技术要求低通滤波器：p：通带截止频率（又称通带上限频率）“一阻带下限截止频率0。：通带允许的最大衰减、阻带允许的最小衰减(戊，a5的单位dB)Qp：通带上限角频率Qs：阻带下限角频率3、IIR数字滤波器的设计步骤：D按一定规则将给出的数字滤波器的技术指标转换魏模拟低通滤波器的技术指标。2)根据转换后的技术指标设计模拟低通滤波器G(S)；3)再按一定的规则将G(S)转换成(2)。4)若是高通、带通或带阻数字滤波器则将它们的技术指标先转化为低通模拟滤波器的技术指标，然

5、后按上述步骤2)设计出低通GW,再将G(S)转换为所需的“(2)。4.几种不同类型的滤波器的介绍：因为我们设计的滤波器的冲击响应一般都为实数，所以有G(S)G*(S)=G(S)G(F)IEQ=IG(JQS)广这样，如果我们能由,。，%,OS求出IG(JQ)I,那么就容易得到所需要的G(S)O不同类型的叵。)的表达式，代表了几种不同类型的滤波器。(1)巴特沃思(BUtterWOrth)滤波器：IG(JQ)=919C为待定常数，N为待定的滤波器阶次。(2)切比雪夫I型(ChebySheV-1)滤波器：5.巴特沃思模拟低通滤波器的设计由于每一个滤波器的频率范围将直接取决于设计者的应用目的，因此必然是

6、千差万别。为了使设计规范化，我们需要将滤波器的频率参数作归一化处理。设所给的实际频率为Q(或，),归一化后的频率为几，对低通模拟滤波器，令九=OgP显然，4=1=po又令归一化复数变量为p,p=以,显然P=JQ/Qp=sQp所以巴特沃思模拟低通滤波器的设计可按以下三个步骤来进行。将实际频率O规一化(2)求C和NC2=10%。-1N=IgJ柒二般这样C和N可求。若令0P=3dB,则C=I,这样巴特沃思滤波器的设计就只剩一个参数N,这时G()2=7T7=+g1Q)2N(3)确定G(S)因为P=JX,根据上面公式有G(p)G(-p)=+(pj)2N=+(T)N.2N1+J1)NP2N=O解得这样可得

7、Gk(P)=1(p-Pkp-PN+-k)1一2Z2/+N1、P-2pcos(占q万)+1求得G(P)后，用s/。，代替变量，即得实际需要得G(5)。6.用双线性Z变换法设计HR数字低通滤波器s平面到z平面的映射关系2z-1S=Ts2+1称为双线性Z变换，由此关系求出1 +(7；/2)5Z-1-(12)s.2Sin3/2)2 =)Tscos(62)2Q=tan(g/2)即TS=2arctan(7/2)2因为设计滤波器时系数方会被约掉，所以又有1+sZ=1s=ta32)-2arctan(三)实验内容题一：试用双线性Z变换法设计一低通数字滤波器，给定技术指标是2=100Hz,=300Hz,%=3c1

8、b,%=2。dB)抽样频率G=迹瓯提示：首先应该得到角频率力，然后再(1) 将数字滤波器的技术要求转换为模拟滤波器的技术要求。(2) 设计低通滤波器G(S)由依次求出4,4，再求出N,可得G(P)G(S)=G(p)工然后由P=E转换成G(S)(3)由G(S)求(z)源程序：fp=100;fs=300;%通带最大衰减%阻带最小衰减%抽样频率%数字角频率%角频率畸变ap=3;as=20;Fs1=IOOO;wp=2*pi*fpFs1;ws=2*pi*fsFs1;Fs=Fs1Fs1;T=1Fs;Op=2T*tan(wp2);Os=2T*tan(ws2);NWn=buttord(Op,Os,ap,as,

9、s,);zpk=buttap(N);ba=zp2tf(z,p,k);BA=1p21p(b,a,Op);BzAz=bi1inear(B,A,Fs*T2);Hw=freqz(Bz,Az,256,Fs*Fs1);p1ot(w,abs(H),r,);tit1e(低通滤波器);%G(p)的零点,极点,增益%G(p)的分子、分数系数%H(s)的分子、分数系数%H(z)的分子、分数系数%Hgridon;同理可以设计出低通滤波器和带通及带阻滤波器。MAT1AB参考程序和仿真内容functionx=ffts(mode,M)Nfft=2M;x=zeros(1,Nfft);%定义一个长度为Nfft的一维全0数组if

10、mode=1forn=O:Nfft-1x(n+1)=sin(2*pi*nNfft);endend%定义一个长度为Nfft的单周期正弦序列ifmode=2forn=O:Nfft-1x(n+1)=sin(4*pi*nNfft);endend%定义一个长度为Nfft的双周期正弦序列ifmode=3endifmode=4endifmode=5endifmode=6endforn=0:Nfft/2-1x(n+1)=sin(4*pi*n/Nfft);end%定义一个长度为Nfft/2的正弦序列，后面一半为0序列。forn=0:Nfft-1x(n+1)=square(2*pi*nNfft);endforn=

11、0:Nfft-1x(n+1)=square(2*pi*nNfft);endforn=0:Nfft/2-1x(n+1)=square(4*pi*nNfft);endn=0:Nfft-1;subp1ot(2,1,1);stem(n,x);axis(0Nfft-I1.1*min(x)1.1*max(x);x1abe1(,Points-,);y1abe1(x(n);X=abs(fft(x,Nfft);subp1ot(2,1,2);stem(n,X);axis(0Nfft-I1.1*min(X)1Pmax(X);x1abe1(,frequency-,);y1abe1(!X(k)!t);图2带通滤波器幅频特性曲线图3带通滤波器相频特性曲线图二0.1参考文献：程佩青.数字信号处理教程.清华大学出版社（第三版）2007.22MAT1AB及在电子信息课程中的应用（第三版）陈怀琛吴大正编著【3】张合全；张奎。基于MAT1AB的IIR数字滤波器设计，科技信息（学术研究）,2008/12