IIR数字滤波器专业课程设计
IIR数字带通滤波器设计 一、设计内容 1、设计任务 利用双线性变换法基于MATLAB设计一个IIR带通滤波器。 2、设计要求 其中带通中心频率为ωp00.5π,;通带截止频率ωp10.4π,ωp20.6π;通带最大衰减αp3dB;阻带最小衰减αs15dB;阻带截止频率ωs20.7π。 3、设计分析 数字滤波器是对数字信号实现滤波线性时不变系统。数字滤波实质上是一个运算过程,实现对信号运算处理。输入数字信号(数字序列)经过特定运算转变为输出数字序列,所以,数字滤波器本质上是一个完成特定运算数字计算过程,也能够了解为是一台计算机。描述离散系统输出和输入关系卷积和差分方程只是给数字信号滤波器提供运算规则,使其根据这个规则完成对输入数据处理。时域离散系统频域特征 其中、分别是数字滤波器输出序列和输入序列频域特征(或称为频谱特征),是数字滤波器单位取样响应频谱,又称为数字滤波器频域响应。输入序列频谱经过滤波后,所以,只要根据输入信号频谱特点和处理信号目标, 合适选择,使得滤波后满足设计要求,这就是数字滤波器滤波原理。 数字滤波器依据其冲激响应函数时域特征,可分为两种,即无限长冲激响应IIR数字滤波器和有限长冲激响应FIR数字滤波器。IIR 数字滤波器特征是,含有没有限连续时间冲激响应,需要用递归模型 来实现,其差分方程为 系统函数为 设计IIR滤波器任务就是寻求一个物理上可实现系统函数Hz,使其频率响应Hz满足所期望得到频域指标,即符合给定通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数。 二、设计方法 1、设计步骤 1 依据任务,确定性能指标在设计带通滤波器之前,首先依据工程实际需要确定滤波器技术指标 带通滤波器阻带边界频率相关中心频率ωp0几何对称,所以ws1wp0- ws2-wp00.3π通带截止频率wc10.4π,wc20.6π;阻带截止频率wr10.3π,wr20.7π;阻带最小衰减αs3dB和通带最大衰减αp15dB; 2 用Ω2/T*tanw/2对带通数字滤波器Hz数字边界频率预畸变,得到带通模拟滤波器Hs边界频率关键是通带截止频率ωp1,ωp2;阻带截止频率ωs1,ωs2转换。 为了计算简便,对双线性变换法通常T2s 通带截止频率wc12/T*tanwp1/2tan0.4π/20.7265 wc22/T*tanwp2/2tan0.6π/21.3764 阻带截止频率wr12/T*tanws1/2tan0.3π/20.5095 wr22/T*tanws2/2tan0.7π/21.9626 阻带最小衰减αs3dB和通带最大衰减αp15dB; 3 利用低通到带通频率变换公式λΩ2-Ω02/B*Ω将模拟带通滤波器指标转换为模拟低通滤波器指标。 Bwc2-wc10.6499 normwr1wr12-w02/B*wr12.236 normwr2wr22-w02/B*wr22.236 normwc1wc12-w02/B*wc11 normwc2wc22-w02/B*wc21 得出,normwc1,normwr2.236 模拟低通滤波器指标normwc1,normwr2.236,αp3dB,αs15dB 4 设计模拟低通原型滤波器。用模拟低通滤波器设计方法得到模拟低通滤波器传输函数Has;借助巴特沃斯Butterworth滤波器、切比雪夫Chebyshev滤波器、椭圆Cauer滤波器、贝塞尔Bessel滤波器等。 5 调用lp2bp函数将模拟低通滤波器转化为模拟带通滤波器。 6 利用双线性变换法将模拟带通滤波器Has转换成数字带通滤波器Hz。 2、程序步骤框图 开始 ↓ 读入数字滤波器技术指标 ↓ 将指标转换成归一化模拟低通滤波器指标 ↓ 设计归一化模拟低通滤波器阶数N和3db截止频率 ↓ 模拟域频率变换,将GP变换成模拟带通滤波器Hs ↓ 用双线性变换法将Hs转换成数字带通滤波器Hz ↓ 输入信号后显示相关结果 ↓ 结束 三、MATLAB程序及运行 MATLAB 程序以下 clear wp00.5*pi;wp10.4*pi;wp20.6*pi; Ap3;ws20.7*pi;As15;T2; 数字带通滤波器技术指标 ws1wp0-ws2-wp0; 计算带通滤波器阻带下截止频率 wc12/T*tanwp1/2;wc22/T*tanwp2/2; wr12/T*tanws1/2;wr22/T*tanws2/2; w02/T*tanwp0/2; 频率预畸变 Bwc2-wc1; 带通滤波器通带宽度 normwr1wr12-w02/B*wr1; normwr2wr22-w02/B*wr2; normwc1wc12-w02/B*wc1; normwc2wc22-w02/B*wc2; 带通到低通频率变换 if absnormwr1absnormwr2 normwrabsnormwr2 else normwrabsnormwr1 end normwc1; 将指标转换成归一化模拟低通滤波器指标 Nbuttordnormwc,normwr,Ap,As,s; 设计归一化模拟低通滤波器阶数N和3db截止频率 [bLP,aLP]butterN,normwc,s; 计算对应模拟滤波器系统函数Gp [bBP,aBP]lp2bpbLP,aLP,w0,B; 模拟域频率变换,将GP变换成模拟带通滤波器Hs [b,a]bilinearbBP,aBP,0.5; 用双线性变换法将Hs转换成数字带通滤波器Hz wlinspace 0,2*pi,500; hfreqzb,a,w; subplot2,1,2; plotw,absh; grid on xlabelwrad ylabel|Hjw| title频谱函数 subplot2,2,1; plotw,20*log10absh; axis[0,2*pi,-120,20]; grid on xlabelwrad ylabel20*lg|Hjw|db titl