专业课程设计格式.doc
课 程 设 计 报 告课程名称: 学生姓名: 学 号: 专业班级: 指导教师: 完成时间: 201?年?月?日 评阅意见: 评阅教师 日期 报告成绩: XXXXXXXX的设计二号黑体占三行1.课程设计目的宋体小三,段前段后空0.5行通过对常用数字滤波器的设计和实现,掌握数字信号处理的工作原理及设计方法;熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法,掌握利用数字滤波器对信号进行滤波的方法,掌握数字滤波器的计算机仿真方法,并能够对设计结果加以分析宋体,小四,1.5倍行距。2.课题要求采用双线性变换法设计一IIR数字带通滤波器,抽样频率为,性能要求为:通带范围从250Hz到400Hz,在此两频率处衰减不大于, 在150Hz和500Hz频率处衰减不小于,采用巴特沃思型滤波器。3.设计原理3.1双线性变换原理宋体,四号占一行先用非线性频率压缩方法,将整个频率轴上的频率范围压缩到所有的公式都用公式编辑器编写之间,再用转换到平面上。也就是说,第一步先将整个平面压缩映射到平面的一条横带里;第二步再通过变换关系将此横带变换到整个平面上去。这样就使平面与Z平面建立了一一对应的单值关系,消除了多值变换性,也就消除了频谱混叠现象,映射关系如图1所示。图1双线性变换的映射关系五号黑体为了将平面的整个虚轴压缩到平面轴上段,可以通过以下的变换实现所有的公式都用公式编辑器编写 (1标号右对齐)表1列出了将给定的对数字滤波器接转换为对低通模拟滤波器的和转换公式。表1 双线性变换和频率预校正的计算公式五号黑体变换类型变换关系频率预校正备注低通变换表中字体用五号,宋体.其中,T为抽样周期,f为模拟频率中心频率其中,分别为数字带通滤波器通带的上下边界角频率,或数字带阻滤波器阻带的上下边界角频率。4.设计思路4.1 设计步骤为了将平面的整个虚轴压缩到平面轴上段,可以通过以下的变换实现(1)(2)(4)4.2 设计过程(1)对应的各数字的频率(2) 求D1,E1。采用归一化原型的低通滤波器作为变换原型,则低通到带阻的变换中所需常数D1,E1分别为(3)求。它是满足数字带阻滤波器要求的归一化原型模拟低通滤波器的阻带起始频率,所以有:应该取=min(),这是因为在较小的频率处如果满足衰减要求,则在较大的频率处一定满足衰减要求。所以=。(4)求N。N是归一化原型模拟低通滤波器的阶数。将代入巴特沃思低通滤波器频率响应的公式可求得。按程序运行结果得N=3。(5)求,根据求得的N,查表可求得。(6)求数字带阻滤波器的系统函数,也可按程序求出,.再用分式,所以有5实验程序及结果5.1 实验程序w2=2.*pi.*(fp2./fsa);wp1=2*pi*fp1*T;wp2=2*pi*fp2*T;ws2=2.*pi.*(fs2./fsa);cosw0=(sin(w1+w2)./(sin(w2)+sin(w1);w0=acos(cosw0);bw=wp2-wp1;Wp=(cosw0-cos(w2)./sin(w2);Ws=(cosw0-cos(ws2)./sin(ws2);N,omgn=cheb1ord(Wp,Ws,Rp,As,'s')5.2 实验结果图2 数字带阻滤波器的幅频响应图及相频特性图6.心得体会在之前数字信号与处理的学习以及完成课后的作业的过程中,已经使用过MATLAB,对其有了一些基础的了解和认识。通过这次练习是我进一步了解了信号的产生、采样及频谱分析的方法。 以及其中产生信号和绘制信号的基本命令和一些基础编程语言。让我感受到只有在了解课本知识的前提下,才能更好的应用这个工具;并且熟练的应用MATLAB也可以很好的加深我对课程的理解,方便我的思维。这次设计使我了解了MATLAB的使用方法,学会分析滤波器的优劣和性能,提高了分析和动手实践能力。同时我相信,进一步加强对MATLAB的学习与研究对我今后的学习将会起到很大的帮助。在课程设计中也出现了很多问题,主要有以下几个方面:(1)对IIR滤波器的设计原理掌握不牢,在设定通带阻带频率上碰到点麻烦。在经过很多的尝试以后,我们采用了巴特沃思低通滤波器,并把阻带通带设定在了一个合理的位置上。(2)对Matlab相关函数不太了解,很多函数都是现学的。对于Matlab里面的各种滤波器的函数也不时很了解,经过了大量的实验之后,才有了本实验程序与结果。(3)在结果屡试都达不到要求的情况下,很急躁导致程序出现细微的错误都察觉不出来。后经过多次检查及翻阅书籍才达到要求。7.参考资料1 丁玉美,高西全.数字信号处理.西安:西安电子科技大学出版社,20102 A.V.奥本海姆,R.W.谢弗.数字信号处理.北京:科学出版社,2009