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

巴特沃斯低通、切比雪夫低通、高通IIR滤波器设1120141454焦奥一、设计思路IIR滤波器可以分为低通、高通、带通、带阻等不同类型的滤波器，而以系统函数类型又有巴特沃斯、切比雪夫等滤波器。其中巴特沃斯较为简单，切比雪夫较为复杂；低阶比高阶简单，但却有着不够良好的滤波特性。在满足特定的指标最低要求下，低阶、巴特沃斯滤波器能更大程度地节省运算量以及复杂程度。滤波器在不同域内分为数字域和模拟域。其中数字域运用最广泛。在设计过程中，一般是导出模拟域的滤波器，之后通过频率转换变为数字域滤波器，实现模拟域到数字域的传递。在针对高通、带通、带阻的滤波器上，可以又低通到他们的变换公式来进行较为方便的转换。综上，IIR滤波器的设计思路是，先得到一个满足指标的尽可能简单的低通模拟滤波器，之后用频域变换转换到数字域。转换方法有双线性变换法、冲激响应不变法等。虽然方法不同，但具体过程有很多相似之处。首先将数字滤波器的指标转换为模拟滤波器的指标，之后根据指标设计模拟滤波器，再通过变换，将模拟滤波器变换为数字滤波器，是设计IIR滤波器的最基本框架。以下先讨论较为简单的巴特沃斯低通滤波器。二、巴特沃斯低通滤波假设需要一个指标为0~4hz内衰减小于3db、大于60hz时衰减不小于30db的滤波器。其中抽样频率为400hz。以双线性变换方法来设计。首先将滤波器转换到模拟指标。根据双线性变换得到这就得到了模拟域的指标。由巴特沃斯的方程得到当N取大于最小值的整数时，解出N=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时，衰减超出指标。以下讨论切比雪夫滤波器。三、切比雪夫滤波器切比雪夫滤波器是另一种滤波器，使用较为广泛的是低通滤波器，也能够通过变换转化为其他滤波类型的滤波器。以下设计切比雪夫低通滤波器，仍用双线性变换法。频率脉冲响应法与双线性变换法的区别仅仅在于两个域的转换方面，在设计模拟滤波器的方面，是没有影响的。频率脉冲法计算更加复杂，而且有混叠效应；双线性法通过预畸减小了误差，并且没有混叠。给出指标：0~30hz内衰减小于3db、大于90hz时衰减不小于40db的滤波器。其中抽样频率为200hz。根据双线性变换得到得到模拟域指标。求求NN=2查表,得到N=2的归一化系统函数代入并令新的得到Matlab仿真Wc=203.81;b=[0 0 0.5*Wc^2];a=[1 0.64*Wc 0.71*Wc^2];[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 -