椭圆IIR高通数字滤波器.docxVIP

南华大学电气工程学院数字信号处理课程设计设计题目：椭圆高通IIR数字滤波器专 业：10级电子02班学生姓名: 王 超学 号:20104470245起迄日期:2013年12月26日—2014年1月11日指导教师: 陈忠泽总体设计思路根据rand函数产生的0-1之间随机值计算数字高通滤波器的技术指标。将数字高通滤波器的边界频率转换成模拟高通滤波器的边界频率。将模拟高通滤波器的的技术指标转换成归一化的椭圆模拟低通滤波器的技术指标。设计归一化椭圆低通滤波器原型(s)。通过频率转换将(s)转换成过渡的模拟高通滤波器。采用双线性变换法将过渡的转换成数字高通滤波器。手工计算完成椭圆IIR数字高通滤波器的初始设计数字高通滤波器的设计指标通带截止频率, 阻带截止频率,通带最大衰减,阻带最小衰减其中，，为两个用MATLAB的rand函数随机产生的0-1之间的值，且有。如图1.1.1所示，图1.1.1（1）通带截止频率：==1.7181（2）阻带截止频率：==0.8749（3）通带最大衰减：（4）阻带最小衰减：将数字高通滤波器设计指标转换成模拟高通滤波器设计指标频率转换关系为：=，为了计算方便，令T=2s 计算得:1.1594rad/s0.4677 rad/s转换成归一化椭圆滤波器设计指标归一化边界频率椭圆滤波器的归一化低通原型一般选择关于通带边界频率归一化的低通系数。在虚轴上低通到高通滤波器的映射关系为如下频率变换公式综合上述，得：归一化通带边界频率归一化阻带边界频率确定椭圆滤波器最小阶数定义频率的选择性因数k为：令: =-0.1175-0.1973i=-0.1677-0.1046i==椭圆滤波器的最小阶次N可表示为：经计算，最小阶数为N=4。设计归一化椭圆低通滤波器原型(s)椭圆滤波器的幅频响应函数公式为：其中，是与通带衰减有关的参数，为截止频率，为雅可比椭圆函数。这里：将和的几何平均值作为频率归一化的基准频率，即：这时通带截止频率的归一化值和阻带截止频率的归一化值分别为：在以为归一化基准频率的条件下，通带、阻带阶值归一化频率互为倒数。已知：最低阶数 N=4通带最大衰减阻带最大衰减通带边界频率查表可得归一化传递函数：(p)=将归一化椭圆低通滤波器转换为模拟高通滤波器从低通到高通滤波器的映射关系为因此，模拟高通滤波器的传递函数为模拟高通滤波器转换为数字高通滤波器从Z平面上的一点映射到S平面一点双线性变换为因此,数字传输函数和原模拟函数之间的关系为检验检验上述变量计算结果为了检验上述变量是否计算正确，我编写了从第一步到第三步所有变量计算的M文件脚本。M文件和计算结果如图1.1和图1.2所示：图1.1 M文件脚本图2.2运行结果结论：经检验，手工计算与计算机计算得出的结果完全相符。检验最终得出的数字高通滤波器传递函数根据数字信号处理教材上关于椭圆滤波器的一些介绍，我编写了直接设计椭圆数字滤波器的M文件，源程序和运行结果如图2.3和图2.4所示：图2.3图2.4结论：通过matlab函数直接设计椭圆数字高通滤波器检验，结果与上述手工计算相符。滤波器的不同结构对性能指标的影响利用直接型结构构建数字滤波器在MATLAB中可以利用FDATOOL工具箱构建不同类型的数字滤波器。在此为了使对比效果明显，我们不妨先将将上述初步设计的椭圆数字IIR带阻滤波器的设计参数的字长（即转移函数中分子、分母各项前的系数）进行保留小数点后3位的进一步的缩减。缩减后的参数如下:den=[ 0.051, -0.165, 0.235, -0.165, 0.051 ]num=[ 1.000, 1.117, 1.347, 0.750, 0.267 ]图4.1.1 filter coefficients工具工作界面直接型的结构流图如图4.1.2所示：x(n) 0.051 y(n) -0.165-1.117 0.235 -1.347 -0.165-0.750 0.051 -0.267图4.1.2直接型I型结构流图选择filter structure选项框中的 Direct-Form I选项，点击窗口下方的Import Filter按钮，构建直接1型结构的椭圆高通滤波器，结果如图4.1.3所示：图4.1.3Direct-Form I型结构的滤波器幅频响应图读图可以得Direct-Form I结构的滤波器技术指标（ws,wp单位为；,,单位为dB）如表1所示：表4.1.1Direct-Form I结构滤波器对性能指标的影响性能指标初始设计指标Direct-Form I