设计的巴特沃斯低通滤波器指标与等波纹低通滤波器相同 ，通带截至频率4kHZ，阻带截至频率为5kHZ，通带衰减不大于0.5DB，阻带衰减不小于40DB，采样频为20kHZ。
设计巴特沃斯低通滤波器程序如下：
fp = 4000; fs = 5000; Fs = 20000; wp = 2*pi*fp/Fs; ws = 2*pi*fs/Fs; wp = tan(wp/2);
ws = tan(ws/2); % 通带最大衰减为0.5dB，阻带最小衰减为40dB [N, wn] = buttord(wp, ws, 0.5, 40, 's'); % 模拟低通滤波器极零点 [z, p, k] = buttap(N); % 由极零点获得转移函数参数 [b, a] = zp2tf(z, p, k); % 由原型滤波器获得实际低通滤波器 [B, A] = lp2lp(b, a, wp); [bz, az] = bilinear(B, A, .5); [h, w] = freqz(bz, az, 256, Fs); subplot(2,3,2); plot(w, abs(h))
grid;title('巴特沃斯低通滤波器')
11
巴特沃斯低通滤波器的基于Matlab的幅频响应曲线 如下：
七、结论心得
经过这几天紧张的课程设计，我感觉收获特别大，首先，
12
我学会了如何使用MATLAB软件；其次，又由于对基础知识掌握不牢固，又重新翻开了数字信号处理课本，复习研究其中的原理，在熟悉理论知识环节和编程环境的前提下，开始思考如何利用编程实现我们的课设要求；最后，通过查找资料和独立处理一些问题，锻炼了自己的能力，对我今后的学习和工作起到了不可或缺的帮助。
参考文献：
[1]高西全 丁美玉.《数字信号处理》[M].西安：西安电子科技大学出版社，2008.8
[2]陈怀琛.《数字信号处理教程——MATLAB释义与实现》[M].北京：电子工业出版社，2004.12
[3]张德丰.《详解MATLAB数字信号处理》[M].北京：电子工业出版社，2010.6
[4]飞思科技产品研发中心.《MATLAB7辅助信号处理技术与应用》[M].北京：电子工业出版社，2005.3
13