《简单低通滤波器设计及matlab仿真》由会员分享，可在线阅读，更多相关《简单低通滤波器设计及matlab仿真（10页珍藏版）》请在装配图网上搜索。
1、东北大学研究生考试试卷评分考试科目：课程编号： 阅 卷 人: 考试日期：姓 名： xl学 号： 注 意 事 项1. 考前研究生将上述项目填写清楚.2. 字迹要清楚,保持卷面清洁.3. 交卷时请将本试卷和题签一起上交.4. 课程考试后二周内授课教师完成评卷工作,公共课成绩单与试卷交研究生院培养办公室,专业课成绩单与试卷交各学院,各学院把成绩单交研究生院培养办公室.东北大学研究生院培养办公室数字滤波器设计技术指标：通带最大衰减： =3dB，通带边界频率： 100Hz 阻带最小衰减： 20dB 阻带边界频率： 200Hz采样频率：Fs=200Hz目标：1、根据性能指标设计一个巴特沃斯低通模拟滤波器。
2、2、通过双线性变换将该模拟滤波器转变为数字滤波器。原理：一、模拟滤波器设计每一个滤波器的频率范围将直接取决于应用目的，因此必然是千差万别。为了使设计规范化，需要将滤波器的频率参数作归一化处理。设所给的实际频率为（或），归一化后的频率为，对低通模拟滤波器令，则，。令归一化复数变量为，则。所以巴特沃思模拟低通滤波器的设计可按以下三个步骤来进行。(1)将实际频率规一化(2)求c和N这样c和N可求。根据滤波器设计要求3dB，则C1，这样巴特沃思滤波器的设计就只剩一个参数N，这时(3)确定因为，根据上面公式有由解得，1，2，2N这样可得求得后，用代替变量，即得实际需要得。二、双线性变换法双线性变换法是将
3、平面压缩变换到某一中介平面的一条横带里，再通过标准变换关系将此带变换到整个z平面上去，这样就使平面与平面之间建立一一对应的单值关系，消除了多值变换性。为了将平面的轴压缩到平面的轴上的到一段上，可以通过以下的正切变换来实现：这样当由经变化到时，由经过0变化到，也映射到了整个 轴。将这个关系延拓到整个平面和平面，则可以得到再将平面通过标准变换关系映射到平面，即令得到同样对z求解，得到这样的变换叫做双线性变换。为了验证这种映射具有s平面的虚轴映射到z平面单位圆上的特性，考虑 ,，得除了使s平面的虚轴映射到单位圆上之外，s平面的左半部分映射到单位圆的内部，s平面的右半部分映射到单位圆的外部。设计过程一
4、、计算滤波器阶数N和截止频率c。根据公式：代入数据，计算可得：ps=0.5Kps=0.99885N=3.3151所以取N=4，即滤波器为四阶滤波器。计算3dB截止频率，根据公式 计算可得c=112.6096二、计算系统函数已知四阶归一化低通巴特沃斯模拟滤波器系统函数为：将S用S/c替代，求出系统函数：Ha(s) = b/ (s4+a3s3+ a2s2+ a1s+ a0)其中：b=2.5063*10(11)a3=1.8489*10(5)a2=1.7092*10(6)a1=9.2560*10(8)a0=2.5063*10(11)三、程序实现模拟滤波器（MATLAB）clear; close all
5、fp=100;fs=200;Rp=3;As=20; %滤波器指标N,fc=buttord(fp,fs,Rp,As,s) %计算阶数N和3dB截至频率fcB,A=butter(N,fc,s); %设计低通巴特沃斯模拟滤波器hf,f=freqs(B,A); %计算模拟滤波器频率响应plot(f,20*log10(abs(hf), r)grid; xlabel(频率Hz);ylabel(幅度dB)title(模拟低通滤波器);axis(0,250,-25,5)line(0,250,-3,-3);line(100,100,-25,5);line(0,250,-20,-20);line(200,200,
6、-25,5);响应曲线如下四、通过双线性法变换将模拟滤波器转变为数字滤波器首先根据公式将 H (s)转换为H(z)。计算可得：H(z)分子系数为：Mz = 0.0197 0.0786 0.1179 0.0786 0.0197H(z)分母系数为：Nz = 1.0000 -1.7271 0.0871 1.8682 -0.9138 1-1.7271 Z-1+0.0871 Z-2+1.8682 Z-3-0.9138 Z-40.0197+0.0786 Z-1+0.1179 Z-2+0.0786 Z-3+0.0197 Z-4H(z) = 五、程序实现数字滤波器（MATLAB）T=0.005;M=2.506
7、3*10(11); N=1,1.8489*10(5) ,1.7092*10(6) ,9.2560*10(8) ,2.5063*10(11)Mz,Nz=bilinear(M,N,1/T); %对模拟滤波器双线性变换Mz,Nzh1,w1=freqz(Mz,Nz); %数字滤波器的幅频响应figureplot(w1/pi,20*log10(abs(h1), r); grid;xlabel(/);ylabel(幅度(dB); title(数字低通滤波器);axis(0,1.1,-160,20)响应曲线如下六、设计级联型滤波器根据1-1.7271 Z-1+0.0871 Z-2+1.8682 Z-3-0.
8、9138 Z-40.0197+0.0786 Z-1+0.1179 Z-2+0.0786 Z-3+0.0197 Z-4H(z) =编写matlab程序求级联型滤波器系数。程序如下：a=0.0197,0.0786,0.1179,0.0786,0.0197;b=1.0000,-1.7271,0.0871,1.8682,-0.9138;sos,g=tf2sos(a, b) %求级联型结构系数结果如下：sos = 1.0000 2.3550 1.4405 1.0000 0.3905 -0.6025 1.0000 1.6348 0.6942 1.0000 -2.1176 1.5166g = 0.0197所
9、以可由级联型结构系数写出系统函数：H(z) = 0.0197( )( )1+0.3905 Z-1-0.6025 Z-21+1.6348 Z-1+0.6942Z-21+2.355Z-1+1.4405Z-21-2.1176 Z-1+1.5166 Z-2级联型结构图如下图所示：X(n)Y(n)0.01972.3550.6942-0.39052.11761.6348-1.5166Z-1Z-1Z-11.44050.6025Z-1结果分析1、模拟滤波器性能分析观察响应曲线，在通带边界频率100Hz处，幅度最大衰减为3 dB；在阻带边界频率200Hz处，幅度衰减达到20dB。因此模拟滤波器设计符合要求。2、
10、双线性变换所得数字低通滤波器分析双线性变换优点：是频率坐标变换是线性的，即，如果不考虑频率混叠现象，用这种方法设计的数字滤波器会很好的重现原模拟滤波器的频率特性。另一个优点是数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，时域逼近性好。双线性变换法的缺点：数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变。例如，一个模拟微分器，它的幅度与频率是线性关系，但通过双线性变换后，就不可能得到数字微分器。另外，一个线性相位的模拟滤波器经双线性变换后，滤波器就不再有线性相位特性。虽然双线性变换有这样的缺点，但它目前仍是使用得最普遍、最有成效的一种设计工具。这是因为大多数滤波器都具有分段常数的频响特性，如低通、高通、带通和带阻等，它们在通带内要求逼近一个衰减为零的常数特性，在阻带部分要求逼近一个衰减为的常数特性，这种特性的滤波器通过双线性变换后，虽然频率发生了非线性变化，但其幅频特性仍保持分段常数的特性。 双线性变换比脉冲响应法的设计计算更直接和简单。由于s与z之间的简单代数关系，所以从模拟传递函数可直接通过代数置换得到数字滤波器的传递函数。置换过程：频响： 这些都比脉冲响应不变法的部分分式分解便捷得多。
简单低通滤波器设计及matlab仿真
