数字信号处理 切比雪夫Ⅰ型带阻IIR数字滤波器的设计DOCWord文档下载推荐.docx
《数字信号处理 切比雪夫Ⅰ型带阻IIR数字滤波器的设计DOCWord文档下载推荐.docx》由会员分享,可在线阅读,更多相关《数字信号处理 切比雪夫Ⅰ型带阻IIR数字滤波器的设计DOCWord文档下载推荐.docx(14页珍藏版)》请在冰豆网上搜索。
学生姓名:
二零一四年一月=0.2551、0.5060、=0.6991、=0.8909课程设计任务书题目切比雪夫型带阻IIR数字滤波器的设计学生姓名学号专业班级设计内一、设计内容设计切比雪夫带阻滤波器,四个参数通过MATLAB的rand函容与要求数随机在01间产生,且满足。
(1)通带下截止频率=*rad
(2)阻带下截止频率=*rad(3)阻带上截止频率=*rad(4)通带上截止频率=*rad(5)通带最大衰减=1db(6)阻带最大衰减=60db二、设计要求
(1)滤波器的初始设计通过手工完成
(2)在计算机辅助计算的基础上分析滤波器结构对其性能指标的影响(至少选择两种以上合适的滤波器结构进行分析)(3)在计算机辅助计算的基础上分析滤波器参数的字长对其性能指标的影响(4)以上各项要有理论分析和推导,源程序以及表示计算结果的图表1.手工计算完成初始设计过程1.1整体设计思路
(1)
(2)(3)(6)(5)(4)图1-1-1整体设计流程图1.2计算过程
(1)通过MATLAB的rand函数随机在01间产生,且满足。
程序与运行结果如下:
rand(1,4)ans=0.25510.50600.69910.8909所以有:
通同带下截止频率=0.5*rad阻同带下截止频率=0.85*rad阻同带上截止频率=0.65*rad通同带上截止频率=0.35*rad通带最大衰减=1db阻带最大衰减=60db
(2)数字带阻指标转化为模拟带阻指标转化公式:
令T=2s,预畸变校正得到模拟边界频率:
=0.4236rad/s=5.7780rad/s=1db=1.0190rad/s=1.9558rad/s=80db而*=2.4476*=1.9930不满足修正=1.2515rad/s所以,转化后的模拟带阻的各项指标为:
=0.4236rad/s=5.7780rad/s=1db=1.2515rad/s=1.9558rad/s=60db(3)将模拟带阻指标转化为模拟低通指标带阻低通的映射为:
其中Bw为阻带宽度,=-=0.25*,取低通原型通带边界频率=1且=时的映射为=-带入上式得到:
=7.6026=1db=80db(4)设计切比雪夫模拟低通滤波器求滤波器阶数N和参数其中,=1965取符合条件的N的最小整数,取N=7;
=0.0.3493求归一化极点Pk其中,=0.375求出的四个归一化极点分别为:
-0.1395+0.9834j-0.3369+0.4073j-0.3369-0.4073j-0.1395-0.9834j求归一化系统函数带入后得到由于所求的是切比雪夫滤波器,故。
又由于要用到的是归一化的样本低通滤波器,故不必去归一化。
(5)将模拟低通滤波器转化为数字低通滤波器采用双线性变换的方法将将模拟低通滤波器转化为数字低通滤波器转换公式为所得数字低通滤波器系统函数为(6)在数字域内坐频带变换将数字低通滤波器转换为数字带阻滤波器转换过程如下:
样本低通DF截止频率则有其中,、分别为所求带阻滤波器上下通带的截止频率其中,为阻带中心频率。
求得转换公式为:
最终得到所求数字带阻滤波器的系统函数为:
2.结果分析2.1结合matlab工具进行幅频响应分析把相关参数输入进matlab程序里面,画出所成幅频响应:
图2-1-1matlab中运行的程序图2-1-2实际频率响应前面计算中求得N=7;
=0.3493,其中N影响过渡带的宽度,同时也影响通带内波动的疏密,N等于通带内最大值与最小值的总个数;
与通带内允许的波动幅度有关。
从图中也可看出,通带具有等波纹特性,在阻带内单调下降后再单调上升。
Fdatool中理想图形为:
图2-1-3matlab完成的理想频率响应由图2-1-2与图2-1-3对比可看出,实际计算的结果与matlab运行的结果存有一定的误差,但仍在可接受的范围内,满足设计的要求:
(1)通同带下截止频率=0.2551*rad
(2)阻同带下截止频率=0.5060*rad(3)阻同带上截止频率=0.6991*rad(4)通同带上截止频率=0.8909*rad(5)通带最大衰减=1db(6)阻带最大衰减=60db产生误差的原因可能是每次小数均取到小数点后的四位数,每次运算均有取舍,多次运算后造成误差积累。
2.2.分析结合滤波器结构对其性能指标的影响
(1)直接型画结构流图如下:
图2-2-1直接型结构流图图2-2-2直接型结构幅频响应表2-2-1直接型结构性能比较表性能指标初始指标直接型误差0.25510.2551270.000730.50600.49326610.01273390.69910.70996090.010866090.89090.8657710.02513(db)10.99852250.0014775(db)60600直接型结构由直接型结构经过延迟之路合并后得到,对于8阶差分方程只需要8个延迟单元,所需延迟单元最少,软件实现时可节省存储单元,硬件实现时可节省寄存器。
但同时此种结构也存在一定缺点:
虽然系统的零点由系数决定,系统的极点由系数决定。
但系数对滤波器的性能控制作用不明显,这是因为它们与系统函数的零极点关系不明显,因而调整困难;
这种结构极点对系数的变化过于灵敏,从而使系统频率响应对系统的变化过于灵敏,各个的有限字长误差都会影响到同一个极点,也就是对有限精度(有限字长)运算过于敏感,容易出现不稳定或产生较大误差。
为了克服直接性结构的不足,可通过对的变换,获得其他形式的IIR滤波器结构,下面用级联型结构对滤波器进行分析。
(2)级联型对的变换得到级联型系统函数:
图2-2-3级联型结构流图图2-2-4级联型型结构幅频响应表2-2-2级联结构性能比较表性能指标初始指标直接型误差0.25510.25476070.00033930.50600.4930420.0015580.69910.70996090.01086090.89090.878540.0155(db)10.99783790.0021621(db)60600由图2-2-4可看出:
调整、就能单独调整滤波器的第i对零点而不影响其他的零极点,调整、就能单独调整滤波器的第i对极点而不影响其他的零极点;
可以用零极点配对的方法,吧共轭零极点或相近的零极点组合在一个二阶滤波器中,可降低字长系数的敏感程度;
任意调整在级联结构中的位置是不改变整个系统函数的,这样可以通过交换级联顺序来减小有限字长所带来的系统误差;
在实际实现时,在各级联节之间应适当进行点电平的放大和缩小,以使变量值不会太大或太小。
不能太大是为了避免运算中产生溢出现象,不能太小是为了防止信号与噪声的比值太小。
2.3分析滤波器参数的字长对其性能指标的影响
(1)直接型分析字长影响A.滤波器系数量化字长均为16位,分子分母小数部分量化字长不同图2-3-1分子分母小数部分字长为2位图2-3-2分子分母小数部分字长为4位图2-3-3分子分母小数部分字长为8位图2-3-4分子分母小数部分字长为14位图中虚线为参考的理想字长下生成的滤波器的幅频响应曲线,实线为参数字长不同时的滤波器的幅频响应曲线,从图中可看出,当滤波器系数量化字长均为16位,随着分子分母小数部分量化字长增大,滤波器的幅频响应曲线与要求曲线逐渐重合,但速度较慢,当分子分母小数部分量化字长为14位时达到设计要求。
B.分子分母小数部分量化字长均为14位,滤波器系数量化字长不同图2-3-5滤波器系数量化字长为4位图2-3-6滤波器系数量化字长为12位图2-3-7滤波器系数量化字长为14位图2-3-8滤波器系数量化字长为16位从图中可看出,当分子分母小数部分量化字长均为14位时,随着滤波器系数量化字长增大,滤波器的幅频响应曲线与要求曲线逐渐重合,但速度较慢,滤波器系数量化字长为14位时才比较吻合,当滤波器系数量化字长为16位时达到设计要求。
(2)级联型分析字长影响A.滤波器系数量化字长均为16位,分子分母小数部分量化字长不同图2-3-9分子分母小数部分字长为2位图2-3-10分子分母小数部分字长为4位图2-3-11分子分母小数部分字长为8位图2-3-12分子分母小数部分字长为14位从图中可看出,当滤波器系数量化字长均为16位,随着分子分母小数部分量化字长增大,滤波器的幅频响应曲线与要求曲线逐渐重合,并且速度很快,当分子分母小数部分量化字长为4位时已基本与要求曲线逐渐重合。
B.分子分母小数部分量化字长均为14位,滤波器系数量化字长不同图2-3-13滤波器系数量化字长为4位图2-3-14滤波器系数量化字长为12位图2-3-15滤波器系数量化字长为14位图2-3-16滤波器系数量化字长为16位从图中可看出,当分子分母小数部分量化字长均为14位时,随着滤波器系数量化字长增大,滤波器的幅频响应曲线与要求曲线逐渐重合,但速度较慢,滤波器系数量化字长为14位时仍存在很大偏差,当滤波器系数量化字长为16位时达到设计要求。
分析总结:
滤波器参数量化字长越长,设计出的滤波器就越符合设计指标要求,误差越小,稳定性越好,滤波性能也越好,但总的来说,滤波器系数量化字长比分子分母小数部分量化字长对滤波器性能影响更大一些,并且,量化误差对直接型结构滤波器性能的影响比级联级结构大得多。
实验总结与体会此次课程设计让我学到了很多,最开始感到无从下手,后来自己看书,在错误的路上走了一遍又一遍,请教老师,查阅资料,最终于有了比较清晰的思路,但在完成的过程中又遇到了各种各样的问题,但我都一一解决了,整个课程设计下来,我对matlab有了整体的认识,对于数字和模拟滤波器的整个设计过程都有了比较完整的理解,虽然仿真结果并不是非常理想,但我的收获还是很多的。
感谢陈忠泽老师耐心的耐心讲解与分析,我自己也深刻的认识到了所学知识的不足与实际应用中的极大欠缺,以后的学习过程中,我会更加注意。
参考文献1高息全丁玉美.数字信号处理M.西安:
西安电子科技大学出版社,2008.082陈佩青数字信号处理教程M.北京:
清华大学出版社,2013.53从玉良数字信号处理原理及其MATLAB实现M.北京:
电子工业出版社,2009.7