[讲解]数字滤波器设计.docVIP

巴特沃斯滤波器的设计 1.1巴特沃斯滤波器的幅度平方函数及其特点 巴特沃斯模拟滤波器幅度平方函数的形式是 （5-6） 式中N为整数，是滤波器的阶次。=0时，=1时；当=时，=1/ ，所以又称为3dB截止频率。 1.2幅度平方函数的极点分布及的构成 将幅度平方函数写成s的函数 （5-7） 此式表明幅度平方函数有2N个极点，极点用下式表示 k=0,1,2，……（5-8） 这2N个极点分布在s平面半径为的圆上，角度间隔是π/N弧度。N=3时，极点间隔为π/3弧度或60度。极点对虚轴是对称的，且不会落在虚轴上。当N是奇数时，实轴上有极点；当N为偶数时，则实轴上没有极点。巴特沃斯滤波器的N个极点为 （5-9） 则的表达式即滤波器的系统函数为 （5-10） 1.3频率归一化问题 式（5-10）即为所求滤波器的系统函数，可看出与有关，即使滤波器的幅度衰减特性相同，只要不同，就不一样。为使设计统一，可将所以的频率归一化。这里采用对3dB截止频率归一化。 1.4设计步骤 总结以上讨论，低通巴特沃斯滤波器的设计步骤如下： 1）根据技术指标和，用式（5-17）求出N。 2）按照（5-13），求出归一化极点，代入（5-12），得到归一化传输函数。也可以直接查表. 3)将去归一化。将p=s/代入，得到实际的滤波器传输函数。如果技术指标没有给出3dB截止频率，可以按照式（5-18）或（5-19）求出。 2仿真过程及仿真图 2.1 MATLAB程序 Fs=1000; fp=100; fs=300; Ap=3; As=20; Wp=fp/(Fs/2);%计算归一化角频率 Ws=fs/(Fs/2); [N,Wc]=buttord(Wp,Ws,Ap,As);%计算阶数和截止频率 [b,a]=butter(N,Wc);%计算H(z)分子、分母多项式系数 [H,F]=freqz(b,a,500,Fs);%计算H(z)的幅频响应,freqz(b,a,计算点数,采样速率) subplot(2,2,2) plot(F,20*log10(abs(H))) xlabel(频率(Hz)); ylabel(幅度(dB)) axis([0 500 -30 3]); grid on subplot(2,2,1) plot(F,abs(H)); xlabel(频率(Hz)); ylabel(幅度); title(低通滤波器) grid on; subplot(2,2,3) pha=angle(H)*180/pi; plot(F,pha); xlabel(频率(Hz)); ylabel(相位(dB)) grid on 2.2仿真波形 3设计结论 首先，这个滤波器的设计的整体思路是先设计一个巴特沃斯低通滤波器，因为模拟滤波器的设计已经很成熟了，有严格的设计思路，现成的公式和图表以及曲线可以查阅。设计出模拟滤波器的转换成满足技术指标的数字滤波器。再学习了MATLAB编程，通过实践来学习软件的使用，根据提供的数据，查阅了基本的MATLAB的命令，中途在归一化这个方面遇到点小问题，不过经过反复的调试得到基本满意的波形。 4体会 总的来说，通过这次的课程设计我对MATLAB有了全面的认识，对数字信号处理的知识又有了深刻的理解，让我感受到只有在充分理解课本知识的前提下，才能更好的应用这个工具；并且熟练的应用 MATLAB也可以很好的加深我对课程的理解，方便我的思维。这次设计使我了解了MATLAB的使用方法，学会分析滤波器的优劣和性能，提高了分析和动手实践能力。同时我相信，进一步加强对 MATLAB 的学习与研究对我今后的学习将会起到很大的帮助！