二阶IIR滤波器的设计

第1
章 二阶IIR 滤波器的设计思路与步骤
IIR 数字滤波器可用一个n 阶差分方程
y(n)=Σb r x(n-r)+Σa k y(n-k),
或用它的Z 域系统函数:
对照模拟滤波器的传递函数:
不难看出,数字滤波器与模拟滤波器的设计思路相仿,其设计实质也是寻找一组系数{b,a},去逼近所要求的频率响应,使其在性能上满足预定的技术要求;不同的是模拟滤波器的设计是在S 平面上用数学逼近法去寻找近似的所需特性H(S),而数字滤波器则是在Z 平面寻找合适的H(z)。

IIR 数字滤波器的单位响应是无限长的,而模拟滤波器一般都具有无限长的单位脉冲响应,因此与模拟滤波器相匹配。

已经证明,冲击响应不变法和双线性变换法能较好地担当此任,则在此基础上,数字滤波器的设计就可首先归结为模拟滤波器的设计了。

数字滤波器的设计步骤如图所示。

第2章二阶IIR滤波器的设计方法
二阶IIR滤波器的设计方法有多种,可归纳为下述两种。

2.1 传统设计方法
根据前述设计思路,首先设计一个模拟原型滤波器(截止频率为1rad/s的低通滤波器),然后在模拟域(S平面)进行频率变换,将模拟原形滤波器转换成所需类型(指定截止频率的低通、高通、带通、带阻)的模拟滤波器,再将其数字离散化,从S平面映射至Z平面,得到所需技术指标的数字滤波器。

上述过程中,也可先将模拟原型离散化,得到数字原型滤波器,继而在数字域(Z平面)进行频率变换,得到所需类型的数字滤波器。

模拟滤波器到数字滤波器的转换可在时域进行也可在频域实现,时域转换的关键是要使数字滤波器与模拟滤波器时域响应的采样值相等,以保持其瞬态特性不变,常用的是冲击响应不变法。

频域变换法必须使得数字滤波器在
-π≤ω≤π范围内的幅频特性与模拟滤波器在-π/T≤Ω≤π/T 范围内的幅
频特性一致,即保证S平面与z平面上幅频特性的一一单值对应关系,常用的是双线性变换法。

2.2 计算机辅助设计方法
传统设计方法思路清晰,步骤详尽,可参阅公式、手册循章而行。

但由于计算繁琐,手工计算大多只能用来进行简单低阶选频滤波器(如LP,HP,BP及BS等)的设计。

计算机辅助设计方法是集电路理论、网络图论、数值分析、矩阵运算、元件建模、优化技术、高级计算机语言等多交叉学科于一身的新领域,它把计算机的快速、高精度、大存储容量、严格的逻辑判断和优良的数据处理能力与人的思维创造能力充分结合起来,极大地简化了数字滤波器的设计过程。

在优秀科技应用软件MATLAB的信号处理工具箱中,提供了一整套模拟、数字滤波器的设计命令和运算函数,方便准确,简单易行,使得设计人员除了可按上述传统设计步骤快速地进行较复杂高阶选频滤波器的计算、分析外,还可通过原型变换法直接进行各种典型数字滤波器设计,即应用MATLAB设计工具从模拟原型直接变换成满足原定频域指标要求的数字滤波器。

此外,MATLAB软件还提供了强大的数字滤波器优化设计功能。

即预先确定一种最佳设计准则,然后直接求得在该准则下滤波器系统的传递函数{b,a}。

这种最优设计方法可方便地用于任意幅频特性要求的多带通复杂滤波器系统的设计。

第3章各种设计方法的MATLAB实现
下面以一带通数字滤波器的设计为例,说明基于MATLAB的三种实现方法。

3.1 基于模拟低通原型的MATLAB实现
%通过模拟低通原型滤波器进行数字带通设计程序：
fp=480 %模拟低通通带上限频率
fs=520 %模拟低通阻带下限频率
wp=2*pi*fp %模拟低通通带上限角频率
ws=2*pi*fs %模拟低通阻带下限角频率
rp=3 %通带波动
rs=20 %阻带衰减
%巴特沃斯模拟低通原型滤波器设计
[n,wn]=buttord(wp,WS,rp ,rs,'s')
[z,p,k]=buttap(n) %模拟低通原型零、极点系数
[b1,a1]=zp2tf(z,p,k) %零、极点系数转换为传递函数
%巴特沃斯模拟低通原型滤波器频率响应
[hl,w1]=freqs(b1,a1)
mag1=abs(h1)
%模拟低通原型滤波器幅频特性曲线
subplot(221);semilogx(w1,mag1)
fw=40 %模拟带通滤波器带宽频率
bw=2*pi*fw %模拟带通滤波器带宽角频率
% 由模拟低通原型变换为模拟带通滤波器
[b2,a2]=lp2bp(b1,a1,wn,bw) %模拟带通滤波函数系数
%巴特沃斯模拟带通滤波器频率响应
[h2,w2]= freqs(b2,a2)
mag2 =abs(1l2)
%模拟带通滤波器幅频特性曲线(db)
subplot(222);plot(20*log10(mag2))
%数字带通滤波器直接设计程序
fp= [480,520];fs=[450,550] %模拟通带、阻带频率
rp=3;rs=20 %通带波动、阻带衰减
fo=10000 %采样频率
%频率指标变换
wp=2*pi*fp/f0 %数字通带频率
ws=2*pi*fs/fo %数字阻带频率
%切比雪夫1型数字带通滤波器直接设计
[n,wn]=cbeblord(wp/pi,ws/pi,rp,rs)
[b,a]=chebyl(n,rp,wn) %数字带通滤波器系数
%切比雪夫1型数字带通滤波器频率响应
[h,w]= freqz(b,a,128,10000)
mag=abs(h;pha=unwrqp(angle(h))
subplot(321);plot(w,mag) %幅频曲线
subplot(322);plot(w,pha) %相频曲线
hi=impz(b,a) %冲击响应
subplot(324);plot(hi) %响应曲线
n=0:500;t=n/fc
x1=2*square(2*pi*500*t) %500Hz方波信号subplot(323);plot(t,x1) %500Hz方波波形
yi=conv(hi,x1) %时域卷积输出
subplot(326);plot(yi) %卷积输出波形
y1=filter(b,a,x1) %数字滤波函数输出
subplot(325); stem(y1) %数字滤波器输出波形
第4章结束语
基于MATLAB的信号处理工具箱为数字滤波器设计带来了全新的实现手段,设计快捷方便,仿真波形直观。

上述三种设计方案均可实现设计指标,但以直接原型变换法最为简便。

实际应用中,数字滤波器也可以对连续时间信号进行处理,但需要先对连续信号进行A/D变换,经数字滤波后,再经D/A转换得到所需要的连
续信号。

参考文献
[1] 倪养华,等．数字信号处理与实现．上海：上海交通大学出版社,1998；
[2]施阳．等．MATLAB语言工具箱．西安：西北工业大学出版社,1999。