基于Matlab的椭圆带阻IIR数字滤波器设计(数字频率变换).doc

PAGE 2 - 数字信号处理课程设计 数字频带变换的椭圆IIR带阻滤波器设计 南华大学 电气工程学院 电子091班 苗润武 20094470138 一 手工计算完成椭圆IIR数字带阻滤波器初始设计 1、设计要求： 滤波器的设计指标要求为 ⑴通带下截止频率, ⑵阻带下截止频率, ⑶阻带上截止频率, ⑷通带上截止频率, ⑸通带最大衰减, ⑹阻带最小衰减。 其中为我学号的后两位。我的学号：20094470138，所以=38。 由此计算性能指标（保留小数点后1位）得： ⑴通带下截止频率：==1.3 rad/sample ⑵阻带下截止频率==2.0 rad/sample = 3 \* GB2 ⑶阻带上截止频率：==4.7 rad/sample = 4 \* GB2 ⑷通带下截止频率：==5.3 rad/sample = 5 \* GB2 ⑸通带最大衰减： = 6 \* GB2 ⑹阻带最小衰减： 由得到性能指标在MATLAB中的常用形式（保留小数点后1位）： = 1 \* GB2 ⑴ rad /sample， = 2 \* GB2 ⑵ rad/sample， = 3 \* GB2 ⑶ rad/sample， = 4 \* GB2 ⑷ rad/sample 2、数字边界频率转换成模拟边界频率： 转变换关系为： 其中，令T=2s计算模拟边界频率（保留小数点后1位）得： = 1 \* GB2 ⑴0.8 rad/s = 2 \* GB2 ⑵1.6 rad/s = 3 \* GB2 ⑶-1.0 rad/s = 4 \* GB2 ⑷-0.5 rad/s 另外： = 5 \* GB2 ⑸ = 6 \* GB2 ⑹ 3、将带阻滤波器的设计要求转换为低通原型滤波器的设计要求 = 1 \* GB2 ⑴=1 = 2 \* GB2 ⑵=2.0 = 3 \* GB2 ⑶ = 4 \* GB2 ⑷ 将和的几何平均值作为频率归一化的基准频率，即： =2.0 定义频率的选择性因数为两个截止频率之比： 4、设计低通原型滤波器转移函数H(s) 椭圆滤波器的幅频响应函数公式为： 其中，为波纹系数，表示波纹情况；为截止频率，为椭圆函数。这里： ==0.5 此外，这时通带截止频率的归一化值和阻带截止频率的归一化值分别为： =0.7071 在以为归一化基准频率的条件下，通带、阻带阶值归一化频率互为倒数。 令： ， ， 分别计算得： 由入阶数计算公式 计算得 最低阶数N=5 由N=5，，=2.0查表可确定归一化转移函数（参数保留小数点后15位）得： 5、将椭圆模拟低通原型滤波器转换得数字滤波器转移函数 由转换公式：（参数保留小数点后15位）得： 6、将数字椭圆低通滤波器变换为带阻滤波器 利用数字频带转换公式： 其中：， 其中：， 由数字低通滤波器转换为数字带阻滤波器，其阶数增加为2N，即为10阶的滤波器 则初步运算所得滤波器转移函数的系数为： den=[1,0,-2.12658637235093,0,2.92411213439205,0,-2.37562859929732,0,1.21306969101691,0,-0.312305136646181] num=[0.0207958528193616,0,0.0543457659047958,0,0.0861892398331069,0,0.0861892398331069,0,0.0543457659047958,0,0.0207958528193616] 二 滤波器的不同结构对性能指标的影响 在理想状态下，对于同一个传递函数几乎对应着无数种等效结构，然而这些结构却并不一定都能实现。在无限参数字长的情况下，所有能实现传递函数的结构之间，其表现完全相同。 然而，在实际中，由于参数字长有限的限制，各实现结构的表现并不相同。下面我们就将对比直接型（包括直接I、II型）和级联型两种结构在本例中对性能指标的影响。 在MATLAB中可以利用FDATOOL工具箱构建不同类型的数字滤波器。在此为了使对比效果明显，我们不妨先将将上述初步设计的椭圆数字IIR带阻滤波器的设计参数的字长（即转移函数中分子、分母各项前的系数）进行保留小数点后3位的进一步的缩减。缩减后的参数如下： den=[1,0,-2.126,0,2.924,0,-2.375,0,1.213,0,-0.312] num=[0.020,0,0.0543,0,0.086,0,0.0861,0,0.054,0,0.021] 将上述参数输入FDATOOL中的filter coefficients工具中如图1所示。 图1 filter coefficients工具工作界面 1、利用直接型结构构建数字滤波器 直接型结构滤波器的典型频率响应表达式为： 对应的差分