巴特沃斯高通数字滤波器设计资料

六、总结：
通过这次课程设计，使我对数字信号处理又有了更深一层的了解。巴特沃斯滤波器看似简单，但是在实际的设计过程之中却遭遇了很多意想不到的困难，许多问题并不是单纯地套用书本上的公式就可以解决了，往往需要通过思考、查阅广泛的资料才能解决。完成这个设计后，我对MATLAB已经有了初步的了解与掌握，同时也提高了自己的分析问题，解决问题的能力。可以说，对于知识的探索是永无止境的，我以后要多加努力，提高自己的专业素质。
3、巴特沃斯滤波器
巴特沃斯滤波器是电子滤波器的一种。巴特沃斯滤波器的特点是通频带的频率响应曲线最平滑。这种滤波器最先由英国工程师斯替芬·巴特沃斯（Stephen Butterworth）在1930年发表在英国《无线电工程》期刊的一篇论文中提出的。
【
巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。 在振幅的对数对角频率的波特图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。
-
2、设计过程：
（1） 求模拟低通滤波器指标。
（2）设计模拟低通滤波器
，取N=4，
查表得
（3）将模拟滤波器转换为数字滤波器
（4）(
（5）利用MATLAB设计巴特沃斯高通滤波器
源程序：
频率变换函数编程：
wp=*pi;ws=*pi;Ap=3;As=48;
Wp=tan(wp/2);Ws=tan(ws/2);
一阶巴特沃斯滤波器的衰减率为每倍频6分贝，每十倍频20分贝。二阶巴特沃斯滤波器的衰减率为每倍频12分贝、 三阶巴特沃斯滤波器的衰减率为每倍频18分贝、如此类推。巴特沃斯滤波器的振幅对角频率单调下降，并且也是唯一的无论阶数，振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器阶数越高，在阻频带振幅衰减速度越快。其他滤波器高阶的振幅对角频率图和低级数的振幅对角频率有不同的形状。
"
数字信号处理随着计算机技术信息技术的进步获得了飞速的发展。数字信号处理已广泛应用于科学研究和工程技术的各个领域，是新一代IT工程师必须掌握的信息处理技术。它在越来越多的应用领域中迅速替代传统的模拟信号处理技术，并且开辟出许多新的领域。
数字信号处理有很多深奥的数学概念，理论也相对抽象，而且是一门理论与实践密切结合的课程。我们通过课程设计深入掌握课程内容，深入理解与消化关于巴特沃斯滤波器的基本理论，锻炼我们独立解决问题的能力，培养我们的创新意识，加强我们的实践学习。
七、参考资料：
1、《数字信号处理》 西安电子科技大学出版社
2、《数字信号处理实践教程》 华中科技大学出版社
3、《信号与线性系统》 高等教育出版社
4、设计IIR数字滤波器的频率变换法
（1）、归一化模拟低通原型到数字高通滤波器的变换
设归一化模拟低通原型滤波器的系统函数是 ，p为模拟域内的拉普拉斯变量。在模拟域内从低通到高通变换，即以 代替p,有
进行反归一化后，即以 代入上式，得到反归一化后的高通滤波器的传输函数
；
—
再进双线性变换，得数字高通滤波器的系统函数为
Bw=W2-W1
W0=sqrt(W1*W2)
（3）低通模拟滤波器到高通模拟滤波器变换函数lp2hp
调用格式如下。
[bs,as]=lp2hp(b,a,Wh)
该函数将截止频率为1rad/s（归一化截止频率）的模拟低通原型滤波器变换成截止频率为 的高通滤波器。
-
（4）低通模拟滤波器到带阻模拟滤波器变换函数lp2bs
axis([0 1 -50 0]);grid;
xlabel('归一化频率');
ylabel('幅度/dB');
程序运行结果如下：
滤波器阶数N=4
.
分子系数b
+000 +000 +000
分母系数a
+000 +000 +000 +000
其幅度响应函数如图所示:
五、结果分析：
由最终所得到的运行结果以及该滤波器的幅度响应函数图像可以看出，我们所设计的滤波器完全符合最初的要求，是一个巴特沃斯高通滤波器，3dB数字截止频率为 ，阻带下边频 rad，阻带衰减为 。
（2）低通模拟滤波器到带通模拟滤波器的变换函数lp2bp
调用格式如下。
^
[bs,as]=lp2bp(b,a,W0,Bw)
该函数将截止频率为1rad/s（归一化截止频率）的模拟低通原型滤波器变换成具有指定带宽 和中心频率为 的带通滤波器。如果要求滤波器的低端截止频率为 ，高端截止频率为 ，则按下式可计算带宽 和中心频率为 的带通滤波器。
MATLAB信号处理工具箱函数lp2lp、lp2bp、lp2hp、lp2bs是模拟频率变换函数。
（1）低通模拟滤波器到低通模拟滤波器变换函数lp2lp
调用格式说明如下。
[bs,ቤተ መጻሕፍቲ ባይዱs]=lp2lp(b,a,Wl)
该函数将截止频率为1rad/s（归一化截止频率）的模拟低通原型滤波器变换成截止频率为Wl的低通滤波器。
二、设计原理：
1、数字滤波器
所谓数字滤波器，是指输入输出均为数字信号，通过数字运算处理改变输入信号所含频率成分的相对比例，或者滤除某些频率成分的数字器件或程序。正因为数字滤波通过数值运算实现滤波，所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊滤波功能。
（5）数字低通滤波器到数字高通变换
由截频 的数字低通滤波器系统函数 得到截频为 的数字高通滤波器系统函数 ，变换关系为
（6）数字低通滤波器到数字带通变换
由截频为 的数字低通滤波器系统函数 得到中心频率为 ，带通上、下截频分别为 的数字带通滤波器系统函数 ，变换关系为
-
5、MATLAB模拟频率变换函数
2、模拟滤波器
模拟滤波器的理论和设计方法已经发展的相当成熟，且有多种典型的模拟滤波器供我们选择，如巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器、贝塞尔滤波器等。这些滤波器有严格的设计公式、现成的曲线和图表供设计人员使用，而且所涉及的系统函数都满足电路实现条件。这些典型的滤波器各有特点：巴特沃斯滤波器具有单调下降的幅频特性；切比雪夫滤波器的幅频特性在通带或阻带有等波纹特性，可以提高选择性；贝塞尔滤波器通带内有较好的线性相位特性；椭圆滤波器的特性相对前三者是最好的，但通带和阻带内均呈现等波纹特性，相位特性的非线性也稍严重。设计时，根据具体要求选择滤波器的类型。
调用格式如下。
[bs,as]=lp2bs(b,a,Wh)
该函数将截止频率为1rad/s（归一化截止频率）的模拟低通原型滤波器变换成具有指定带宽 和中心频率为 的带阻滤波器。
三、设计步骤与流程框图：
四、设计过程：
1、设计要求：
设计一巴特沃斯高通数字滤波器，3dB数字截止频率为 ，阻带下边频 rad，阻带衰减为 。
[
巴特沃斯高通数字滤波器设计
要求：3dB数字截止频率为 ，阻带下边频 rad，阻带衰减为 。
一、课程设计目的：
数字信号处理（Digital Signal Processing DSP）是20世纪60年代以来，随着信息学科和计算机学科的高速发展而迅速发展起来的一门新兴学科。数字信号处理是利用计算机或其他专用处理设备，以数值计算的方法对信号进行采集、变换、滤波、压缩、传输、估值与识别等加工处理，借以达到提取信息和便于应用的目的的一种技术。
dB
应当注意，损耗函数 和幅频特性函数 只是滤波器幅频响应特性的两种描述方法。损耗函数的优点是对幅频响应 的取值非线性压缩，放大了小的幅度，从而可以同时观察通带和阻带的频响特性的变化情况。
模拟低通滤波器的设计指标参数有 、 、 和 。其中 和 分别称为通带边界频率和阻带截止频率， 称为通带最大衰减（即通带 中允许 的最大值）， 称为阻带最小衰减（即阻带 上允许 的最小值）， 和 的单位是dB。
；
直接由归一化低通原型变换为数字高通滤波器时，变换关系为
， ，
其中， 为模拟滤波器3dB截止频率； 为数字滤波器3dB截止频率。
（2）归一化模拟低通原型到数字带通滤波器的频率变换
直接寻求从模拟低通到数字带通之间的映射关系：
，
]
，
其中， 为模拟滤波器3dB截止频率； 分别是数字带通滤波器的中心频率、通带上边频、通带下边频和带通的上阻带下边频。
按照不同的分类方式，数字滤波器可以有很多种类型，但总起来可以分为两大类：经典滤波器和现代滤波器。经典滤波器的特点是其输入信号中有用的频率成分和希望滤除的成分分别占有不同的频带，通过一个合适的选频滤波器滤除干扰，得到纯净信号，达到滤波目的。但是，如果信号和干扰的频谱相互重叠，则经典滤波器无法有效滤除干扰，最大限度恢复信号，这就需要现代滤波器。现代滤波器是根据随机信号的一些统计特性，在某种最佳准则下，最大限度抑制干扰，同时最大限度恢复信号，达到最佳的滤波效果的目的。
[N,wn]=buttord(Wp,Ws,Ap,As,'s');
fprintf('滤波器阶数N=%.0f\n',N)
[b,a]=butter(N,1,'s');
~
[numa,dena]=lp2hp(b,a,Wp);
[numd,dend]=bilinear(numa,dena,;
disp('分子系数b');
fprintf('%.4e ',numd);
fprintf('\n');
disp('分母系数a');
fprintf('%.4e ',dend);
fprintf('\n');
,
w=linspace(0,pi,1024);
h=freqz(numd,dend,w);