基于matlab的切比雪夫及巴特沃斯低通、高通滤波器的设计.docx

精品文档巴特沃斯低通、切比雪夫低通、高通IIR滤波器设计05941401 1120141454焦奥一、设计思路IIR滤波器可以分为低通、高通、带通、带阻等不同类型的滤波器，而以系统函数类型又有巴特沃斯、切比雪夫等滤波器。其中巴特沃斯较为简单，切比雪夫较为复杂；低阶比高阶简单，但却有着不够良好的滤波特性。在满足特定的指标最低要求下，低阶、巴特沃斯滤波器能更大程度地节省运算量以及复杂程度。滤波器在不同域内分为数字域和模拟域。其中数字域运用最广泛。在设计过程中，一般是导出模拟域的滤波器，之后通过频率转换变为数字域滤波器，实现模拟域到数字域的传递。在针对高通、带通、带阻的滤波器上，可以又低通到他们的变换公式来进行较为方便的转换。综上，IIR滤波器的设计思路是，先得到一个满足指标的尽可能简单的低通模拟滤波器，之后用频域变换转换到数字域。转换方法有双线性变换法、冲激响应不变法等。虽然方法不同，但具体过程有很多相似之处。首先将数字滤波器的指标转换为模拟滤波器的指标，之后根据指标设计模拟滤波器，再通过变换，将模拟滤波器变换为数字滤波器，是设计IIR滤波器的最基本框架。以下先讨论较为简单的巴特沃斯低通滤波器。二、巴特沃斯低通滤波假设需要一个指标为04hz内衰减小于3db、大于60hz时衰减不小于30db的滤波器。其中抽样频率为400hz。以双线性变换方法来设计。首先将滤波器转换到模拟指标。T=1fs=1400c=2fc=8c=cT=0.02s=2fs=120s=sT=0.3根据双线性变换=2Ttan2得到c=25.14s=407.62这就得到了模拟域的指标。由巴特沃斯的方程2=Haj2=11+c2N20lgHaj=-10lg1+c2N20lgHajc-320lgHajs-30dB得到-10lg1+cc2N-3-10lg1+sc2N-30当N取大于最小值的整数时，解出N=2，因此为二阶巴特沃斯低通滤波器。查表，可知二阶巴特沃斯低通滤波器的系统函数为Has=c2s2+2cs+c2代入求得数值Has=632.02s2+32.45s+632.02这样便得到了模拟滤波器的函数，利用双线性变换的公式转换到数字域s=2T 1-z-11+z-1得到数字域滤波器Hz=0.00095+0.0019z-1+0.00095z-21-1.91z-1+0.91z-2至此，设计完成。接下来用matlab实现上述变换设计。b=0 0 632.02;a=1 35.45 632.02;bz,az=bilinear(b,a,400)w=0:500*pi/500;H,w=freqz(bz,az);f=400*w/(2*pi);subplot(121);plot(f,abs(H),k);xlabel(频率);ylabel(幅值);axis(0 50 0 1);subplot(122);plot(f,20*log10(abs(H),k);xlabel(频率);ylabel(DB值);axis(0 100 -70 2);hold onplot(0 4,-3 -3,r);hold onplot(4 4,-180 -3,r);hold onplot(60 60,-180 -48.5,b);hold onplot(0 60,-48.5 -48.5,b);plot(0 200,-30 -30,g);text(6,-3,6hz,3db);text(60,-45,60hz,大于30db);text(60,-30,30db);得到结果bz =0.0009 0.0019 0.0009az =1.0000 -1.9114 0.9152bz为分母，az为分子，从左到右分别为1、负一次方、负二次方。验证表示结果正确。因为N取了2，取为整数，大于计算的值，所以在60hz时，衰减超出指标。以下讨论切比雪夫滤波器。三、切比雪夫滤波器切比雪夫滤波器是另一种滤波器，使用较为广泛的是低通滤波器，也能够通过变换转化为其他滤波类型的滤波器。以下设计切比雪夫低通滤波器，仍用双线性变换法。频率脉冲响应法与双线性变换法的区别仅仅在于两个域的转换方面，在设计模拟滤波器的方面，是没有影响的。频率脉冲法计算更加复杂，而且有混叠效应；双线性法通过预畸减小了误差，并且没有混叠。给出指标：030hz内衰减小于3db、大于90hz时衰减不小于40db的滤波器。其中抽样频率为200hz。T=1fs=1200c=2fc=60c=cT=0.3s=2fs=180s=sT=0.9根据双线性变换=2Ttan2得到c=203.81s=2525.50得到模拟域指标。求=100.11-1=100.3-1=1求NNarcch1100.12-1arcchsc=arcch104-1arcch2525.50203.81=arcch99.99arcch12.39=5.33.2=1.66N=2查表,得到N=2的归一化系统函数Hans=0.50s2+0.64s+0.71代入s=sc=s203.81并令新的s=2T 1-z-11+z-1得到HLPz=0.5+z-1+0.5z-25.82-6.26z-1+3.30z-2=0.086+0.17z-1+0.086z-21-1.08z-1+0.57z-2Matlab仿真Wc=203.81;b=0 0 0.5*Wc2;a=1 0.64*Wc 0.71*Wc2;bz,az=bilinear(b,a,200)w=0:500*pi/500;H,w=freqz(bz,az);f=200*w/(2*pi);subplot(121);plot(f,abs(H),k);xlabel(频率);ylabel(幅值);axis(0 50 0 1);subplot(122);plot(f,20*log10(abs(H),k);xlabel(频率);ylabel(DB值);axis(0 95 -50 1);hold onplot(0 30,-3 -3,r);hold onplot(30 30,-180 -3,r);hold ontext(30,-3,30hz3db);text(40,-48,90hz,40db);满足指标要求。bz = 0.0859 0.1719 0.0859az = 1.0000 -1.0801 0.5682数值正确。以下设计切比雪夫高通滤波器四、切比雪夫高通滤波器求高通滤波器的常用方法是模拟低通转化为数字高通。指标：通带截止频率3khz，衰减不大于3db；阻带上限截止频率2khz，衰减不小于14db。抽样频率10khz。T=1fs=110kc=2fc=6kc=cT=0.6s=2fs=4ks=sT=0.4根据双线性变换，并令c=1C1=ctanc2=1.38求低通原型st=C1cotc2=1.89求阶次N=100.11-1=100.3-1=1Narcch1100.12-1arcchstc=arcch4.91arcch1.89=2.271.25=1.97则N=2，系统函数（因为c为1，故系统函数即原型函数）Hans=0.50s2+0.64s+0.71代入变换式s=C1 1+z-11-z-1=1.38z+1z-1得到HHPz=0.5-z-1+0.5z-23.50+2.39z-1+1.73z-2Matlab仿真程序f=1.38/2;b=0 0 0.5;a=1 0.64 0.71;bz=0.5 -1 0.5;az=3.5 2.39 1.73;w=0:500*pi/500;H,w=freqz(bz,az);f=10000*w/(2*pi);subplot(121);plot(f,abs(H),k);xlabel(频率);ylabel(幅值);subplot(122);plot(f,20*log10(abs(H),k);xlabel(频率);ylabel(DB值);hold onplot(0 3000,-3 -3,r);hold onplot(3000 3000,-180 -3,r);hold onplot(2000 2000,-180 -15.9,b);hold onplot(0 2000,-15.9 -15.9,b);plot(0 5000,-14 -14,g);text(3000,-3,3000hz,3db);text(2100,-21,2000hz,大于14db);text(1000,-10,14db);完成指标。五、总结本次作业分别设计了巴特沃斯低通、切比雪夫低通和高通滤波器，所用到的方法是双线性变换法，模拟域转数字域。滤波器的设计思路大体为：将数字指标转化为模拟指标，设计模拟滤波