实验五双线性变换法设计IIR数字滤波器.doc

PAGE  PAGE 12 实验五、双线性变换法设计IIR 数字滤波器 一、实验目的： 1、熟悉用双线性变换法设计IIR 数字滤波器的原理与方法。 2、掌握数字滤波器的计算机仿真方法。 3、熟悉Batterworth滤波器设计方法及特点 二、实验原理 （一）、IIR数字滤波器的设计步骤： ① 按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标； ② 根据模拟滤波器技术指标设计为响应的模拟低通滤波器； ③ 跟据脉冲响应不变法和双线性不变法把模拟滤波器转换为数字滤波器； ④ 如果要设计的滤波器是高通、带通或带阻滤波器，则首先把它们的技术指标转化为模拟低通滤波器的技术指标，设计为数字低通滤波器，最后通过频率转换的方法来得到所要的滤波器。 在MATLAB中，经典法设计数字滤波器主要采用以下步骤：模拟滤波器原型 buttap,cheb1ap 频率变换 模拟离散化 bilinear,impinvar IIR数字滤波器 IIR数字滤波器设计步骤 （二）、用模拟滤波器设计数字滤波器的方法 1、冲激响应不变法： 冲激响应不变法是从时域出发，要求数字滤波器的冲激响应h(n) 对应于模拟滤波器h(t) 的等间隔抽样。 优点：时域逼近良好；保持线性关系。 缺点：频域响应混叠。只适用于限带低通滤波器和带通滤波器 2、双线性变换法 优点：克服了频域混叠 缺点：高频时会引起畸变 1）冲激响应不变法impinvar 格式：[BZ,AZ]= impinvar（B,A,Fs） 功能：把具有[B,A]模拟滤波器传递函数模型转换为采样频率为Fs的数字滤波器的传递函数模型[BZ,AZ]，Fs默认值为1。 例：一个4阶的Butterworth模拟低通滤波器的系统函数如下： 试用冲激响应不变法求出Butterworth模拟低通数字滤波器的系统函数。 num=1; den=[1,sqrt(5),2,sqrt(2),1]; [num1,den1]=impinvar(num,den) 2）双线性变换法bilinear 格式一：[Zd,Pd,Kd]= bilinear(Z,P,K,Fs) 功能：把模拟滤波器的零极点模型转换成数字滤波器的零极点模型，Fs是采样频率 格式二：[numd,dend]= bilinear(num,den,Fs) 功能：把模拟滤波器的传递函数模型转换为数字滤波器的传递函数模型。 例：一个三阶的模拟Butterworth模拟低通滤波器的系统函数如下： ，试用双线性变换法求出数字Butterworth数字低通滤波器的系统函数。 num=1; den=[1,sqrt(3),sqrt(2),1]; [num1,den1]=bilinear(num,den,1) IIR数字滤波器的频率变换实现 步骤： 按一定的规则将数字滤波器的技术指标转换为模拟低通滤波器的技术指标 根据转换后的技术指标使用滤波器阶数函数，确定滤波器的最小阶数N和截止频率Wc 利用最小阶数N产生模拟低通原型 利用截止频率Wc把模拟低通滤波器原型转换为模拟低通、高通、带通、带阻滤波器 利用冲激响应不变法或双线性变换法把模拟滤波器转换为数字滤波器 表一 IIR滤波器阶次估计 函数名功能说明buttord计算Butterworth滤波器的阶次及截止频率cheb1ord计算ChebyshevⅠ滤波器的阶次cheb2ord计算ChebyshevⅡ滤波器的阶次ellipord计算椭圆滤波器的最小阶次表二 模拟低通滤波器原型设计 函数名功能说明buttapButterworth模拟低通滤波器原型设计cleb1apChebyshevⅠ模拟低通滤波器原型设计cheb2apChebyshevⅡ模拟低通滤波器原型设计ellipap椭圆模拟低通滤波器原型设计表三 模拟滤波器变换函数 函数名功能说明lp2bp模拟低通转换为带通lp2bs模拟低通转换为带阻lp2hp模拟低通转换为高通lp2lp改变模拟低通的截止频率（三）、数字滤波器的设计 数字滤波器的设计参数 滤波器的4个重要的通带、阻带参数为： ：通带截止频率（Hz） ：阻带起始频率（Hz） ：通带内波动（dB），即通带内所允许的最大衰减； ：阻带内最小衰减 设采样速率（即奈奎斯特速率）为，将上述参数中的频率参数转化为归一化角频率参数： ：通带截止角频率（rad/s） ，； ：阻带起始角频率（rad/s） ， 通过以上参数就可以进行离散滤波器的设计。 2、巴特沃斯滤波器设计 巴特沃斯滤波器阶数的选择： 在已知设计参数，，，之后，可利用“buttord”命令可求出所需要的滤波器的阶数和3dB截止频率，其格式为： [n，Wn]=butto