IIR数字滤波器可视化平台设计与实现.docVIP

IIR数字滤波器可视化平台设计与实现 　　【摘要】本设计利用Matlab语言对IIR数字滤波器设计建立一个软件仿真平台,系统由巴特沃斯滤波器、切比雪夫Ⅰ型滤波器、切比雪夫Ⅱ型滤波器、椭圆滤波器等组成,设计结果由时域单位脉冲响应图形及频谱图形表示,可任意修改各系统参数以分析研究滤波器设计结果,形象直观。本设计开发目的是给数字滤波器学习分析与研究提供一个软件平台,可作为《数字信号处理》、《信号与系统》、《自动控制原理》、《通信原理》等课程教学辅助。 　　【关键词】IIR数字滤波器 MatLab语言 仿真 多媒体教学 　　【中图分类号】G41 【文献标识码】A 【文章编号】1009-9646(2008)10(a)-0170-04 　　 　　滤波器的功能就是允许某一部分频率的信号顺利通过,而另外一部分频率的信号则受到较大的抑制,它实质上是一个选频电路。滤波器中,把信号能够通过的频率范围,称为通频带或通带;信号受到很大衰减或完全被抑制的频率范围称为阻带;通带和阻带之间的分界频率称为截止频率;理想滤波器在通带内的电压增益为常数,在阻带内的电压增益为零;实际滤波器的通带和阻带之间存在一定频率范围称过渡带,在通带内的电压增益不是常数而是一定的波动范围称通带波动,在阻带内的电压增益不是零而是一定的波动范围称阻带波动。滤波器指标为：通带边缘频率为ωp,阻带边缘频率为ωs,通带电压增益波动为Rp,阻带电压增益波动为As。由于滤波器设计数学模型复杂,不易理解学习,本文采用Matlab程序,根据给定滤波器参数进行设计,得出设计结果并绘出频域和时域特性图,这种图形化、可视化,形象直观便于理解和学习研究,可提高教学效率。本设计已应用于教学实践中,得到学生普遍好评。 　　 　　1 系统设计与仿真 　　本设计开发目的是给滤波器设计的学习分析与研究提供一个软件平台,也是为“信号与系统”、“数字信号处理”、“自动控制原理”、“通信原理”等课程教学提供多媒体教学辅助平台,方便学生对系统的各种特征进行观察,分析比较,提高教学效率。因此,本设计采用菜单形式,把各滤波器分别独立模块,方便模块间进行比较观察分析。 　　1.1 巴特沃斯低通/高通数字滤波器 　　巴特沃斯低通滤波器原型的平方幅度响应函数为:,式中ωc为低通滤波器的截止频率,N为滤波器的阶数。Butterworth低通滤波器特点是通带内具有最大平坦的幅频特性,且随频率增大平滑单调下降;阶数N越高,特性越接近理想滤波器,过渡带越窄;传递函数无零点,极点等距离分布在以|s|=ωc为半径的圆周上。 　　程序算法：输入滤波器设计参数即通带边界数字频率wp、阻带边界数字频率ws、通带最大衰减Rp、阻带最小衰减As、采样频率fs;数字频率转化为模拟角频率op=wp*Fs、os=ws*Fs;计算最小阶数和截止频率[N,oc]=buttord(op,os,Rp,Rs,s);计算低通滤波器系统函数分子分母系数:低通[b,a]=butter(N,oc,s),高通[b,a]=butter(N,oc,high,s);计算归一化模拟滤波器零点、极点和增益[z,p,k]=tf2zp(b,a);模拟滤波器转化为数字滤波器[bz,az]=impinvar(b,a,fs)或者[bz,az]=bilinear(b,a,fs);确定显示频谱特性的角频率范围w=0:0.01:pi;计算频率响应H=freqz(bz,az,w); 计算单位脉冲响应[h,t]=impz(bz,az); 计算单位阶跃响应[h1]=dstep(bz,az,t);用plot函数绘出增益系数频谱特性图axes(handles.axes1);plot(w,abs(H),-r);分贝系数频谱特性图axes(handles.axes2);plot(w,20*log10(abs(H)),-r);单位脉冲响应特性图axes(handles.axes3);plot(t,real(h));单位阶跃响应特性图axes(handles.axes4);plot(t,real(h1));程序运行仿真结果如图1所示。 　　仿真结果分析：分别绘出增益系数频谱特性图和分贝系数频谱特性图,目的是让两种频谱特性图放在一起进行对比,两种指标的表现能力及优缺点清晰可见。由巴特沃斯低通滤波器的幅频特性曲线图可看出在通带内幅频特性比较平坦,且随着频率增大单调下降;传递函数无零点,极点分布在圆周上;由图1所示给出设计参数,设计结果频谱曲线可看出通带最大衰减Rp=0.99dB,通带指标有富余;阻带最小衰减,As=14.99,能满足指标要求。任意改变设计参数,仿真结果与滤波器设计理论相符合。 　　1.2 巴特沃斯带通/带阻IIR数字滤波器 　　程序算法：输入设计参数wp1、wp2、ws1、ws2、