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

数字信号处理课程设计 数字频带变换的椭圆IIR带阻滤波器设计 南华大学 电气工程学院 电子091班 苗润武 20094470138 一 手工计算完成椭圆IIR数字带阻滤波器初始设计 1、设计要求： 滤波器的设计指标要求为 ⑴通带下截止频率, ⑵阻带下截止频率, ⑶阻带上截止频率, ⑷通带上截止频率, ⑸通带最大衰减, ⑹阻带最小衰减。 其中为我学号的后两位。我的学号：20094470138，所以=38。 由此计算性能指标（保留小数点后1位）得： ⑴通带下截止频率：==1.3 rad/sample ⑵阻带下截止频率==2.0 rad/sample ⑶阻带上截止频率：==4.7 rad/sample ⑷通带下截止频率：==5.3 rad/sample ⑸通带最大衰减： ⑹阻带最小衰减： 由得到性能指标在MATLAB中的常用形式（保留小数点后1位）： ⑴ rad /sample，⑵ rad/sample，⑶ rad/sample，⑷ rad/sample 2、数字边界频率转换成模拟边界频率： 转变换关系为： 其中，令T=2s计算模拟边界频率（保留小数点后1位）得： ⑴0.8 rad/s ⑵1.6 rad/s ⑶-1.0 rad/s ⑷-0.5 rad/s 另外： ⑸ ⑹ 3、将带阻滤波器的设计要求转换为低通原型滤波器的设计要求 ⑴=1 ⑵=2.0 ⑶ ⑷ 将和的几何平均值作为频率归一化的基准频率，即： =2.0 定义频率的选择性因数为两个截止频率之比： 4、设计低通原型滤波器转移函数H(s) 椭圆滤波器的幅频响应函数公式为： 其中，为波纹系数，表示波纹情况；为截止频率，为椭圆函数。这里： ==0.5 此外，这时通带截止频率的归一化值和阻带截止频率的归一化值分别为： =0.7071 在以为归一化基准频率的条件下，通带、阻带阶值归一化频率互为倒数。 令： ， ， 分别计算得： 由入阶数计算公式 计算得 最低阶数N=5 由N=5，，=2.02N，即为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、利用直接型结构构建数字滤波器 直接型结构滤波器的典型频率响应表达式为： 对应的差分方程表达式为： 直接型的结构流图如图2所示： 图2 直接型I型结构流图 选择filter structure选项框中的 Direct-Form I选项，点击窗口下方的Import Filter按钮，构建直接1型结构的椭圆数字IIR带阻滤波器，结果如图3所示。 图3 Direct-Form I型结构的滤波器幅频响应图 读图可以得Direct-Form I结构的滤波器技术指标（，，，1所示： 表1 Direct-Form I结构滤波器对性能指标的影响 性能指标） 初始设计指标 Direct-Form I △ 0.2 0.2033251 0.0033251 0.3 0.2816907 -0.0183093 0.7 0.7183093 0.0183093 0.8 0.7966749 -0.0033251 60 58.423