用Matlab设的低通数字滤波器.docVIP

指导教师评定成绩： 审定成绩： 重 庆 邮 电 大 学 自 动 化 学 院 综合课程设计报告 低通低通滤波器是允许低频讯号通过，而不允许高频讯号通过的滤波器。低通要求： 用MATLAB低通，； 截止频率； 按规定的统一格式设计报告。 第一章 MATLAB仿真 2 1.1方案论证 2 1.2设计方法及原理介绍 2 1.3 matlab编程实现 3 1.4 Simulink仿真结构图 4 第二章 硬件设计方案 4 2.1硬件功能整体设计 4 2.2单片机最小系统 4 2.3 ADC0809和74HC74 5 2.4 DAC0832和UA741 6 2.5 串口部分 6 2.6 DDS部分 7 第三章 软件设计 8 3.1 单片机软件设计 8 3.2 PC机通信软件设计 9 第四章 心得体会 10 参考文献 10 附件1：方波源程序 11 附件2：DDS正弦信号模块 12 附件3：主单片机程序 15 附录4：protues原理图 16 附录5：protel电路图 16 附录6：DDS原理图 17 第一章 MATLAB仿真 1.1方案论证 数字滤波器，按功能低通、高通、带通、带阻、全通滤波器按实现的网络结构或单位抽样响应：无限脉冲响应滤波器（IIR滤波器）、有限脉冲响应滤波器（FIR滤波器）模拟滤波器设计已经有了相当成熟的技术和方法，有完整的设计公式，还有比较完整的图表可以查询，因此。IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，其设计工作量比较小，对计算工具的要求不高。：通带截止频率（又称通带上限频率） ：阻带下限截止频率 ：通带允许的最大衰减 ：阻带允许的最小衰减 （，的单位dB） ：通带上限角频率 ：阻带下限角频率 结合相关资料，总结IIR 巴特沃斯 .按一定规则将给出的数字滤波器的技术指标转换成模拟低通滤波器的技术指标。在这里要求在频率小于200Hz的通带内，幅度特性下降小于3dB；在频率大于480Hz的阻带内，衰减大于15dB，采样频率为1200Hz。 .根据转换后的技术指标设计模拟低通滤波器。因为我们设计的滤波器的冲击响应一般都为实数，所以有这样，如果我们能由，，，求出，那么就容易得到所需要的。巴特沃思(Butterworth)滤波器的为，其中C为待定常数，N为待定的滤波器阶次。总结后，巴特沃思滤波器的设计步骤分为以下五步： 确定技术指标，求得数字边缘频率： ， ， B、将数字低通滤波器的技术指标转换为模拟低通滤波器技术指标 用双线性变换法，预畸变公式为，可得 C、确定滤波器阶数N ， 取N＝1。 D、该一阶巴特沃斯归一化模拟低通滤波器的系统函数为： 反归一化后的模拟低通滤波器为： E、采用双线性变换法转换为数字低通滤波器的系统函数： 脉冲响应不变法设计的滤波器在附近存在哦率混叠失真，从而使其阻带内不满足指标要求，双线性变换法无频率混叠失真，满足要求，但其存在非线性频率失真，且频率越高，失真越明显。 1.3 matlab编程实现 源程序见附件一，在此介绍几个主要命令的使用方法。 （1）.[N,wc]=buttord(wp，ws，αp，αs)?用于计算巴特沃斯数字滤波器的阶数N和3dB截止频率wc。调用参数wp，ws分别为数字滤器的通带、阻带截止频率的归一化值，要求：0≤wp≤1，0≤ws≤1。1表示数字频率pi。 [z0，p0，k0]=buttap(N)?用于计算N阶巴特沃斯归一化(3dB截止频率Ωc=1)模拟低通原型滤波器系统函数的零、极点和增益因子。 [b,a]=zp2tf(z,p,k)零极点形式的传递函数的转换 [bt,at] = lp2lp(b,a,Wo) 低通滤波器截止频率变换 （5）. [H,F] = freqz(B,A,N,Fs)? 计算数字滤波器的频率响应的函数fp=200Hz;阻带下限截止频率fs=480Hz;通带最大衰减ap=3dB; 阻带最小衰减as=15dB; 抽样频率Fs=1200Hz。 运行程序得到图1.1。通过图形，我们可以看出滤波器的低频增益基本是常数1，高频衰减较快，通带上限截止频率、阻带下限截止频率等关键指标均满足设计要求。 1.4 Simulink仿真结构图 经过分析，Simulink仿真中应包含如下几个模块： 1.信号发生模块。采用两个离散的正弦信号发生器，一个用于产生能通过滤波器的理想信号，一个用于产生被滤波器滤除的噪声信号。 2.信号叠加模块。采用一个简单的加法器就能将噪声信号叠加到理想信号上去。 3.数字低通滤波器模块。我们采用simulink中自带的digital fil