数字低通滤波器课设
实验四IIR数字滤波器的设计(1)(2)课案
实验四 IIR 数字滤波器的设计及网络结构一、实验目的1．了解IIR 数字滤波器的网络结构。
2.掌握模拟滤波器、IIR 数字滤波器的设计原理和步骤。
3．学习编写数字滤波器的设计程序的方法。
二、实验内容数字滤波器：是数字信号处理技术的重要内容。
它的主要功能是对数字信号进行处理，保留数字信号中的有用成分，去除信号中的无用成分。
1．数字滤波器的分类滤波器的种类很多，分类方法也不同。
（1）按处理的信号划分：模拟滤波器、数字滤波器 （2）按频域特性划分；低通、高通、带通、带阻。
（3）按时域特性划分：FIR 、IIR2．IIR 数字滤波器的传递函数及特点数字滤波器是具有一定传输特性的数字信号处理装置。
它的输入和输出均为离散的数字信号，借助数字器件或一定的数值计算方法，对输入信号进行处理，改变输入信号的波形或频谱，达到保留信号中有用成分去除无用成分的目的。
如果加上A/D 、D/A 转换，则可以用于处理模拟信号。
设IIR 滤波器的输入序列为x(n)，则IIR 滤波器的输入序列x(n)与输出序列y(n)之间的关系可以用下面的方程式表示：1()()()M Ni j i j y n b x n i a y n j ===-+-∑∑（5-1）其中，j a 和i b 是滤波器的系数，其中j a 中至少有一个非零。
与之相对应的差分方程为：10111....()()()1....MM NN b b z b z Y z H Z X z a z a z ----++==++ （5-2）由传递函数可以发现无限长单位冲激响应滤波器有如下特点： （1） 单位冲激响应h(n)是无限长的。
（2） 系统传递函数H(z)在有限z 平面上有极点存在。
（3） 结构上存在着输出到输入的反馈，也就是结构上是递归型的。
3．IIR 滤波器的结构IIR 滤波器包括直接型、级联型和并联型三种结构：① 直接型：优点是简单、直观。
但由于系数bm 、a k 与零、极点对应关系不明显，一个bm 或a k 的改变会影响H(z)所有零点或极点的分布，所以一方面，bm 、a k 对滤波器性能的控制关系不直接，调整困难；另一方面，零、极点分布对系数变化的灵敏度高，对有限字长效应敏感，易引起不稳定现象和较大误差。
(完整版)fir低通滤波器设计(完整版)
电子科技大学信息与软件工程学院学院标准实验报告（实验）课程名称数字信号处理电子科技大学教务处制表电 子 科 技 大 学实 验 报 告学生姓名： 学 号： 指导教师： 实验地点： 实验时间：14-18一、实验室名称：计算机学院机房 二、实验项目名称：fir 低通滤波器的设计 三、实验学时： 四、实验原理：1. FIR 滤波器FIR 滤波器是指在有限范围内系统的单位脉冲响应h[k]仅有非零值的滤波器。
M 阶FIR 滤波器的系统函数H(z)为()[]Mkk H z h k z-==∑其中H(z)是kz-的M 阶多项式，在有限的z 平面内H(z)有M 个零点，在z平面原点z=0有M 个极点.FIR 滤波器的频率响应()j H e Ω为 0()[]Mj jk k H e h k e Ω-Ω==∑它的另外一种表示方法为()()()j j j H e H e e φΩΩΩ=其中()j H e Ω和()φΩ分别为系统的幅度响应和相位响应。
若系统的相位响应()φΩ满足下面的条件()φαΩ=-Ω即系统的群延迟是一个与Ω没有关系的常数α，称为系统H(z)具有严格线性相位。
由于严格线性相位条件在数学层面上处理起来较为困难，因此在FIR 滤波器设计中一般使用广义线性相位。
如果一个离散系统的频率响应()j H e Ω可以表示为()()()j j H e A e αβΩ-Ω+=Ω其中α和β是与Ω无关联的常数，()A Ω是可正可负的实函数，则称系统是广义线性相位的。
如果M 阶FIR 滤波器的单位脉冲响应h[k]是实数，则可以证明系统是线性相位的充要条件为[][]h k h M k =±-当h[k]满足h[k]=h[M-k],称h[k]偶对称。
当h[k]满足h[k]=-h[M-k],称h[k]奇对称。
按阶数h[k]又可分为M 奇数和M 偶数，所以线性相位的FIR 滤波器可以有四种类型。
2. 窗函数法设计FIR 滤波器窗函数设计法又称为傅里叶级数法。
FIR数字低通滤波器的(汉宁)窗函数法设计
)(9cos 15.0)(12cos 15.0)(1919n R n n R N n n w ⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛--=ππ2.3进行语音信号的采集（1）按“开始”-“程序”－“附件”－“娱乐”－“录音机”的顺序操作打开W ｉｎdo ｗｓ系统中的录音机软件。
如图１所示。
图１ wind ｏｗs 录音机（2)用麦克风录入自己的声音信号并保存成wav 文件。
如图２所示。
图2 保存文件保存的文件按照要求如下:① 音信号文件保存的文件名为“yuxue ｊiao ．wav ”。
②语音信号的属性为“８．00０ＫＨz,8位,单声道 7ＫＢ／秒” ，其它选项为默认。
2.４语音信号的分析将“y ｕx ｕejia ｏ．wav ”语音文件复制到计算机装有Matlab 软件的磁盘中相应图3语音信号的截取处理图在图3中，其中第一个图为原始语音信号；第二个图是截短后的信号图。
图4频谱分析图其中第二个图是信号的ＦFT 结果，其横坐标的具体值是Ｘ(k)中的序号ｋ；第三个图是确定滤波频率范围的参考图，其横坐标的具体值应当是遵循D ＦT 定义式和频率分辨率求得的:∑-===10)()]([)(N n k N W n x n x DFT k X π当k 等于0时， 020j kn Njk knNe eW ==⋅-=π，从数字角频率上看，对应的正好是0=ω即直流的位置,也就是说,在取滤波频段时，当将主要能量（即红色框的部分)保留，其余频段部分的信号滤除。
)]([)(n x DFT k X =相当于是信号)(n x 的实际频谱)]([)(n x DFT ej X w =采样,而)(n x 又是连续时间语音信号)(t x 的采样。
)(k X 的每两个相邻取值之间的频率间隔大小对应到语音信号)(n x 的频谱中去，其频率间隔大小正好是采样结果的长度采样速率===∆L f f f s det f ∆称频率分辨率，其中Hz f s 8000=，10000=L ,ｐ２=sum(s2.^２)-suｍ(s1．^２);SNR1=１0*loｇ10（p１/ｐ２);p3=suｍ（ｓ4.^2)/8０0０;p4=sum(s3.^2)/８０00-sum（s4.^2)／8000;SNR2＝10＊loｇ１０(p３/ｐ4);２．６噪声叠加图5 语音信号与加噪声后语音信号对比图五为语音信号与加噪声后语音信号对。
数字信号处理课程设计-等波纹数字FIR低通滤波器
设计题目：等波纹数字FIR低通滤波器2．对课程设计成果的要求〔包括图表（或实物）等硬件要求〕：滤波器的初始设计通过手工计算完成；在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种以上合适的滤波器结构进行分析）；在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响；以上各项要有理论分析和推导、原程序以及表示计算结果的图表；课程设计结束时提交设计说明书。
3．主要参考文献：[1]高息全丁美玉.《数字信号处理》[M].西安：西安电子科技大学出版社，2008.8[2]陈怀琛.《数字信号处理教程——MATLAB释义与实现》[M].北京：电子工业出版社，2004.12[3]张德丰.《详解MATLAB数字信号处理》[M].北京：电子工业出版社，2010.6[4]飞思科技产品研发中心.《MATLAB7辅助信号处理技术与应用》[M].北京：电子工业出版社，2005.34．课程设计工作进度计划：序号起迄日期工作内容接到题目，搜集资料1 2016.12.26-2016.12.31整理资料，构思设计方案2 2016.12.31-2016.1.3手工计算进行滤波器的初步设计3 2016.1.3-2016.1.5完善初步设计，学习Matlab软件操作4 2016.1.5-2016.1.7通过Matlab软件分析设计内容，逐步落实课题目标5 2016.1.8-2016.1.9上交课程设计，并做细节修改并完成设计6 2016.1.10-2016.1.13主指导教师日期：年月日1.前言数字滤波器(digital filter)是由数字乘法器、加法器和延时单元组成的一种装置，在通信、图像、语音、雷达等许多领域都有着十分广泛的应用。
在数字信号处理中，数字滤波占有极其重要的地位。
目前对数字滤波器的设计有多种方法。
其中Matlab软件已成为设计数字滤波器的强有力工具。
传统的数字滤波器设计过程复杂、计算工作量大、滤波特性调整困难，但利用Matlab信号处理工具箱可以快速有效地实现由软件组成的常规数字滤波器的设计、分析和仿真，极大地减轻了工作量，有利于滤波器设计的最优化。
用Matlab设计的低通数字滤波器
指导教师评定成绩：审定成绩：重庆邮电大学自动化学院综合课程设计报告低通数字滤波器设计单位（二级学院）：自动化学生姓名：专业：自动化班级：0810902学号：指导教师：设计时间：2012 年11月重庆邮电大学自动化学院制摘要本次设计以AT89S52单片机为控制器件，利用AD/DA转换芯片构成低通数字滤波器。
首先利用MATLAB根据要求设计低通滤波器模型并在该平台进行初步仿真得到较为理想的效果。
接着利用proteus软件搭建硬件电路，并对编写的程序进行检测；程序检测无误后，再用Protel画出电路图。
首先利用ADC0809将模拟信号转换成数字信号输入单片机，正弦模拟信号由AD9850模块产生，是一般正常的低频正弦波；另一个噪声信号是通过单片机产生的高频脉冲信号。
然后通过滤波程序将采集的数据进行算法滤波，滤波后的数据通过DAC0832输出，将DA0832的输出信号通过一级运放（UA741）转换成电压输出。
最后将运放输出的结果接到示波器上进行显示。
关键字：AT89C51、DDS、ADC0809、DAC0832、UA741设计题目题目与要求低通滤波器是允许低频讯号通过，而不允许高频讯号通过的滤波器。
本题目要求用单片机（型号任选）制作一个低通数字滤波器，要求：1.根据确定的采样频率，用MATLAB设计低通滤波器模型以满足不同截止频率需求，并进行仿真；2.用单片机设计并制作一个低通数字滤波器，实现数据采集、滤波处理、信号输出等功能（包括软、硬件）；3.滤波器能与PC机通信，设计通信协议，实现截止频率参数的调节、原始信号采样数据与滤波后数据的上传；4.利用高级语言开发工具，设计PC机端的程序，可以调节滤波器参数，接收滤波器上传的原始信号采样数据与滤波后数据，并能绘图、保存；5.设计几种不同频率、类型的噪声，利用滤波器进行试验，记录实验情况并进行分析、总结；6.按规定的统一格式设计报告。
目录第一章 MATLAB仿真 (4)1.1方案论证 (4)1.2设计方法及原理介绍 (4)1.3 matlab编程实现 (5)1.4 Simulink仿真结构图 (5)第二章硬件设计方案 (6)2.1硬件功能整体设计 (6)2.2单片机最小系统 (7)2.3 ADC0809和74HC74 (7)2.4 DAC0832和UA741 (8)2.5 串口部分 (9)2.6 DDS部分 (9)第三章软件设计 (10)3.1 单片机软件设计 (10)3.2 PC机通信软件设计 (12)第四章心得体会 (13)参考文献 (14)附件1：方波源程序 (15)附件2：DDS正弦信号模块 (15)附件3：主单片机程序 (19)附录4：protues原理图 (21)附录5：protel电路图 (22)附录6：DDS原理图 (22)第一章 MATLAB 仿真1.1方案论证数字滤波器，按功能分类，可以分成低通、高通、带通、带阻、全通滤波器；按实现的网络结构或单位抽样响应分类：可以分成无限脉冲响应滤波器（IIR 滤波器）、有限脉冲响应滤波器（FIR 滤波器）。
巴特沃斯数字低通滤波器课程设计
巴特沃斯数字低通滤波器课程设计目录1.题目.......................................................................................... .22.要求 (2)3.设计原理 (2)3.1 数字滤波器基本概念 (2)3.2 数字滤波器工作原理 (2)3.3 巴特沃斯滤波器设计原理 (3)3.4脉冲响应不法 (5)3.5实验所用MATLAB函数说明 (7)4.设计思路 (9)5、实验内容 (9)5.1实验程序 (9)5.2实验结果分析 (13)6.设计总结 (13)7.参考文献 (14)一、题目：巴特沃斯数字低通滤波器二、要求：利用脉冲响应不变法设计巴特沃斯数字低通滤波器，通带截止频率100HZ,阻带截止频率120Hz,采样频率1000HZ，通带最大衰减为0.5HZ，阻带最小衰减为10HZ，画出幅频、相频相应曲线。
并假设一个信号x(t)=sin(2*pi*f1*t)+sin(2*pi*f2*t),其中f1=50HZ,f2=200HZ。
用此信号验证滤波器设计的正确性。
三、设计原理1、数字滤波器的基本概念所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤波器除某些频率成分的数字器件或程序，因此，数字滤波的概念和模拟滤波相同，只是的形式和实现滤波方法不同。
正因为数字滤波通过数值运算实现滤波，所以数字滤波处理精度高、稳定、体积小、质量轻、灵活、不存在阻抗匹配问题，可以实验模拟滤波器无法实现的特殊滤波功能。
如果要处理的是模拟信号，可通过A\DC和D\AC,在信号形式上进行匹配转换，同样可以使用数字滤波器对模拟信号进行滤波。
2、数字滤波器的工作原理数字滤波器是一个离散时间系统，输入x(n)是一个时间序列，输出y(n)也是一个时间序列。
如数字滤波器的系统函数为H(Z),其脉冲响应为h(n),则在时间域内存在下列关系y(n)=x(n) h(n)在Z域内，输入输出存在下列关系Y(Z)=H(Z)X(Z)式中，X(Z),Y(Z)分别为输入x(n)和输出y(n)的Z 变换。
基于切比雪夫1的IIR低通滤波器
课程设计课程设计名称：DSP课程设计专业班级：电信1005班学生姓名：学号：指导教师：李相国课程设计时间：2013年6月数字信号处理专业课程设计任务书说明：本表由指导教师填写，由教研室主任审核后下达给选题学生，装订在设计（论文）首页填表说明1．“课题性质”一栏：A．工程设计；B．工程技术研究；C．软件工程（如CAI课题等）；D．文献型综述；E．其它。
2．“课题来源”一栏：A．自然科学基金与部、省、市级以上科研课题；B．企、事业单位委托课题；C．校、院（系、部）级基金课题；D．自拟课题。
1 需求分析数字滤波器是指完成信号滤波处理功能的，用有限精度算法实现的离散时间线性非时变系统，其输入是一组数字量，其输出是经过变换的另一组数字量。
因此，它本身即可以是用数字硬件装配成的一台完成给定运算的专用数字计算机，也可以是将所需运算编成程序，让通用计算机来执行。
在本次课程设计中，我做的是基于双线性变换法设计 Chebyshev-I 型 IIR 数字低通滤波器，看到这个题目，我们很快就能联系到数字信号处理中的 IIR 数字滤波器的设计，根据以前学习的东西我们知道，要想设计一个数字低通滤波器，我们可以把所给的数字指标转换为模拟参数，通过设计一个模拟低通的滤波器，对设计好的模拟低通滤波器进行数字化就可以得到一个数字低通滤波器。
在本次实验中，我们首先将所给的滤波器数字指标转换为模拟指标，利用Chebyshev-I型函数设计一个 Chebyshev-I 型低通模拟滤波器，然后对所设计好的 Chebyshev-I型低通模拟滤波器进行参数分析，接下来利用双线性变化法将此模拟低通滤波器转变为数字低通滤波器。
2 概要设计3 运行环境PC 机，windows2000及其以上，matlab 软件。
4 开发工具和编程语言MATLAB 软件，编程语言为DSP 设计语言。
5 详细设计第一步：将数字低通滤波器的性能指标转化成模拟低通滤波器指标Chebyshev-I 型IIR 数字低通滤波器的性能指标：p ω=0.25π，s ω=0.4π，p δ=0.01，s δ=0.001程序模块： %转化成模拟量 rp=-20*log10(1-lp); rs=-20*log10(ls); Omgp=tan(wp/2) Omgs=tan(ws/2) Fs=0.5;第二步：利用切比学夫逼近法设计模拟滤波器切比雪夫1低通滤波器的原理：幅度平方函数为： 其中ε为小于1的正数，表示通带内幅度波动的程度，p Ω称为通带截止频率。
低通滤波器的设计
电子综合课程设计课程名称电子技术综合课程设计院（系）机电学院专业电子科学与科技班级10电子（一）班学号201010330140姓名黄玉燕指导老师杨波课程设计题目低通滤波器的设计程设计的内容及要求：一、设计说明设计一个低通滤波器。
低通滤波器在工业现场主要用于信号的滤波，提高有效信号的信噪比。
实际环境下的有效信号一般是传感器输出信号或通信传输的信号。
目前随着计算机技术的快速发展，诞生了很多方便的设计软件，例如：TI公司的FilterPro,就是一款很好的滤波器设计软件，本次设计建议使用。
二、技术指标1．截止频率10Hz，通带增益20dB；截止带增益-30dB。
三、设计要求1．设计建议采用TI公司的FilterPro滤波器设计软件。
2．输入信号为正弦波，要求通过示波器观测输入输出的相位差。
3．要求将输出正弦波整形为方波并设计电路测量输出方波信号频率。
4．可以考虑频率输出信号与频率测量电路之间增加光电隔离。
5．画出电路原理图（元器件标准化，电路图规范化）。
四、实验要求1．根据技术指标制定实验方案。
2．部分仿真或实际验证所设计的电路。
3．进行实验数据处理和分析。
五、推荐参考资料1.姚福安. 电子电路设计与实践[M]济南：山东科学技术出版社，2001年2.阎石. 数字电子技术基础. [M]北京：高等教育出版社，2006年3.刘贵栋主编.电子电路的Multisim仿真实践[M].哈尔滨：哈尔滨工业大学出版社，2008年4.童诗白、华成英主编.模拟电子技术基础,[M]北京：高等教育出版社，2007年.一、概述滤波器（filter），是一种用来消除干扰的器件，它的主要作用是让有用信号尽可能无衰减的通过，对无用信号尽可能大的衰减。
其功能就是得到一个特定频率或消除一个特定频率。
低通滤波器是一个可以通过低频信号，但是对于高频信号会产生抑制或衰减的部件。
理想滤波器电路的频响在通带内应具有一定幅值和线性相移，而在阻带内其幅值应为零。
巴特沃斯低通滤波器课程设计
电路基础课程设计巴特沃斯低通滤波器设计目标:通带边界频率ωc=4396rad/s (f c=700Hz);通带最大衰减αmax=3dB;阻带边界频率ωs=26376rad/s(f s=4200Hz); 阻带最小衰减αmin=30dB;1.设计步骤⑴设计电压转移函数①将给定的电压衰减技术指标进行频率归一化选取归一化角频率ωr=ωc，这样通带边界频率Ωc=ωc/ ωr=1,阻带边界频率Ωs=ωs/ ωr=ωs/ωc。
②根据归一化的技术指标求出电压转移函数巴特沃斯低通滤波器的阶数n=Log(100.1αmin−1) 2Log(Ωs)带入数据求得n=1.93 取整得n=2由a k=2sin(2k−1)π2n，b k=1和H(s)=U out(s)U in(s)=∏A ks2+a k s+b kn2k=1可得到电压转移函数H(s)=U out(s)U in(s)=1s2+√2s+1将转移函数进行反归一化，即另s=sωc 得到实际转移函数H(s)=U out(s)U in(s)=1s243962+√2s4396+1⑵转移函数的实现选取下图作为实现转移函数的具体电路：列节点方程求解转移函数节点1 U1(1R1+1R2+s∗C1)−1R1U in−1R2−s∗C1∗U2=0节点2 (1R2+s∗C2)U2−1R2U1=0又有U out=U3解得H(s)=U outU in=11+(R2+R2)s∗C2+C1C2R1R2s2对比解得的电压转移函数和推得的电压转移函数里各项的系数并且令R1= R2,C1=1μF,可以得到C1=11000000F=1μFR1=250000√21099Ω=321.705ΩR2=250000√21099Ω==321.705ΩC2=12000000F=0.5μF因实验室没有0.5μF的电容因此取C2=0.47μF2.计算机仿真⑴软件环境：Multisim 10⑵电路图：⑶仿真结果：①700Hz下的波形图②4200Hz下的波形图③波特图◎700Hz下衰减2.673dB◎4200Hz下衰减30.491dB3．实验室实际操作因实验室没有0.5μF的电容和321.705Ω的电阻，因此取C2=0.47μFR1=R2=330Ω实际连电路时，选取集成电路块的第1、2、3引脚分别作为放大器的输出端、负端和正端，第4和11引脚作为供电端，C2一端连接电压源的接地线。
数字信号Butterworth滤波器设计方法
一. 设计任务 1. 设计目的：（1）熟悉和巩固模拟滤波器的设计方法和原理（2）掌握Butterworth 滤波器设计方法 （3）实现滤波器设计的有关经典算法（4）熟练掌握使用高级语言程序设计各种要求的数字滤波器 （5）熟练掌握冲激响应不变法。
2.设计技术指标：（1）按要求设计Butterworth 型数字低通滤波器， （2）性能指标如下：① 通带截止频率πω2.0=p ； ② 通带最大衰减αp=3dB ； ③ 阻带起始频率πω3.0=s ； ④ 阻带最小衰减αs=20dB ；二．数字滤波器的设计根据给定技术指标和数字信号处理理论相关知识进行计算；（1）计算模拟低通滤波器的阶次：N=lg )lg(/11011010/10/12c s ΩΩ--δδ 一般情况下计算结果N 为小数，应向上取整数，对于本任务书给出的技术指标，N=6。
（2） 设计模拟低通滤波器：查参考书或教材表格（不同参考书可能形式略有不同）可得归一化原型滤波器系统函数，根据归一化原型滤波器系统函数，得模拟滤波器系统函数如下： csp p G s G Ω==|)()(注：此处也可以不采用查表法，直接求G(S)，即Ha(s)极点，得出模拟滤波器系统函数。
（3） 根据冲激响应不变法，将模拟低通滤波器转化为数字低通滤波器冲激响应不变法原理： 冲激响应不变法是使数字滤波器的单位冲激响应序列()h n 模仿模拟滤波器的单位冲激响应()a h t ，将模拟滤波器的单位冲激响应加以等间隔抽样，使()h n 正好等于()a h t 的抽样值，即满足：()()a h n h nT =式中：T 为抽样周期。
冲激不变法把稳定的()a H s 转换为稳定的()H z 。
由此方法可得到一阶系统的最基本的转换关系为：1111aTs s a e z --⇒+-（5）本题设计结果：H(z)=0.0034118*(21212756.00318.110332.10328.21----+-++z z z z*21214127.01427.119996.11----+-++z z z z *21217358.04041.119679.09676.11----+-++z z z z ) 三、实验代码：wp=0.2*pi; %性能指标： 通带截止频率0.2*pi ws=0.3*pi; %阻带起始频率0.3*pi Rp=3; %通带最大衰减3dBAs=20; %阻带最小衰减αs=20dB T=1;Rip=10^(-Rp/20); Atn=10^(-As/20); OmgP=wp*T; OmgS=ws*T;[N,OmgC]=buttord(OmgP,OmgS,Rp,As,'s'); %选取滤波器阶数N=6 [cs,ds]=butter(N,OmgC,'s');[b,a]=impinvar(cs,ds,T); %用冲击不变法进行转换 [db,mag,pha,grd,w]=freqz_m(b,a);subplot(2,2,1); plot(w/pi,mag); title('Magnitude Response'); xlabel('w(/pi)'); ylabel('|H(jw)|'); axis([0,1,0,1.1]); subplot(2,2,2); plot(w/pi,db); title('Magnitude in dB'); xlabel('w(/pi)'); ylabel('dB'); axis([0,1,-40,5]);subplot(2,2,3); plot(w/pi,pha/pi); title('Phase Response'); xlabel('w(/pi)'); ylabel('pha(/pi)'); axis([0,1,-1,1]);subplot(2,2,4); plot(w/pi,grd); title('Group delay'); xlabel('w(/pi)'); ylabel('Sample'); axis([0,1,0,12]);function[db,mag,pha,grd,w]=freqz_m(b,a)%滤波器的幅值响应(相对、绝对)、相位响应及群延迟 %Usage:[db,mag,pha,grd,w]=freqz_m(b,a) %a 系统函数 H(z)的分母项(对 FIR ，a=1)[H,w]=freqz(b,a,500); mag=abs(H); db=20*log10(mag/max(mag)); pha=angle(H); grd=grpdelay(b,a,w);matlab 仿真结果Magnitude Responsew(/pi)|H (j w )|Magnitude in dBw(/pi)d BPhase Responsew(/pi)p h a (/p i )Group delayw(/pi)S a m p l e四、激响应不变法的频率混叠失真和优缺点分析数字滤波器的频率响应与模拟滤波器的频率响应间的关系为：12()()j ak H e H j k TTωπ∞=-∞Ω-=∑上式表明，数字滤波器的频率响应是模拟滤波器的周期延拓，根据奈奎斯特抽样定理，只有当模拟滤波器的频率响应是严格限带的，且带限于折叠频率[/2,/2]s s -ΩΩ以内时，才能使数字滤波器的频率响应在折叠频率以内，重现模拟滤波器的频率响应而不产生混叠失真。
课程设计-低通滤波器设计（含matlab程序）
课程设计-低通滤波器设计（含matlab程序）2010/2011学年第2 学期学院：信息与通信⼯程学院专业：电⼦信息科学与技术学⽣姓名：学号：课程设计题⽬：低通滤波器设计起迄⽇期： 6 ⽉13 ⽇～6⽉24⽇课程设计地点：指导教师：系主任：下达任务书⽇期: 2011 年 6 ⽉12 ⽇课程设计任务书课程设计任务书⽬录1 设计⽬的及要 (5) 1.1设计⽬的 (5)1.2设计内容和要求 (5)2 设计原理 (5)2.1 FIR滤波器 (5) 2.2窗函数 (6)2.3矩形窗 (7)3 设计过程 (8)3.1 设计流程图 (8)3.2 产⽣原始信号并分析频谱 (8)3.3 使⽤矩形窗设计不同特性的数字滤波器 (10)3.4 信号滤波处理 (11)4 实验结果及分析 (12)5 课程设计⼼得体会 (12)6 参考⽂献 (13)附录： (14)低通滤波器的设计1 设计⽬的及要求1.1设计⽬的设计⼀种低通滤波器并对信号进⾏滤波。
低通滤波器的作⽤是滤去信号中的中频和⾼频成分，增强低频成分。
要求做到:1.了解MATLAB 的信号处理技术；2.使⽤MATLAB 设计低通滤波器，掌握其滤波处理技术；3.对滤波前和滤波后的波形进⾏时域和频域⽐较。
1.2设计内容和要求1.熟悉有关采样，频谱分析的理论知识，对信号作频谱分析；2.熟悉有关滤波器设计理论知识，选择合适的滤波器技术指标，设计低通滤波器对信号进⾏滤波，对⽐分析滤波前后信号的频谱；3.实现信号频谱分析和滤波等有关MATLAB 函数；2 设计原理本次课程设计，我们主要是基于矩形窗的FIR 滤波器来设计⼀个低通滤波器。
2.1 FIR 滤波器FIR 滤波器即有限抽样响应因果系统，其单位抽样响应h(n)是有限长的；极点皆位于z=0处；结构上不存在输出到输⼊的反馈，是⾮递归型的。
其系统函数表⽰为：()()n -1-N 0n z n h z H ∑==普通的FIR 滤波器系统的差分⽅程为：()()()i n x i h n y 1N 0i -=∑-=式中：N 为FIR 滤波器的抽头数；x(n)为第n 时刻的输⼊样本；h(i)为FIR 滤波器第 i 级抽头系数。
电子线路课程设计：低通滤波器
电子线路课程设计课程名称电子线路课程设计院（系）电子信息工程学院专业通信工程班级学号姓名设计题目低通滤波器一．设计任务和要求：设计一个低通滤波器。
设计要求：（禁止使用集成模块）①截止频率：100KHz②通带增益：20dB③截止带增益：-30dB二．设计设备：低通滤波器在工业现场主要用于信号的滤波，提高有效信号的信噪比。
实际环境下的有效信号一般是传感器输出信号或通信传输的信号。
目前随着计算机技术的快速发展，诞生了很多方便的设计软件。
此次课程设计的模拟仿真，我选择使用Filter Wiz RRO.Filter Wiz Pro是一款很好的滤波器设计软件。
三．概述滤波器（filter），是一种用来消除干扰的器件，它的主要作用是让有用信号尽可能无衰减的通过，对无用信号尽可能大的衰减。
其功能就是得到一个特定频率或消除一个特定频率。
本设计为低通滤波器，低通滤波器是容许低于截止频率的信号通过，但高于截止频率的信号不能通过的电子滤波装置。
通过采用Filter Wiz RRO滤波器设计软件，通过输入截止频率100K Hz，通带增益20dB；截止带增益-30dB后经过自动分析处理后直接算出滤波器的性能及所有滤波器原件的值得到低通滤波器电路图。
四．方案论证：本设计的方案为通过采用Filter Wiz Pro滤波器设计软件设计出符合条件的低通滤波器。
Filter Wiz Pro是一款功能强大的滤波器设计软件，其能够帮助用户设计软件，并可进行低通、高通、带通和带阻滤波器等设置。
尽管低通滤波器在现代电子学领域的地位越来越重要，但其设计及定型工作仍是冗长乏味且耗时巨大的。
不过现在有了Filter Wiz Pro，用户就可以比较迅速地设计、优化和仿真一套完整的多级有源滤波器解决方案。
软件使用了精选的TI运算放大器和TI供应商合作伙伴提供的无源组件，因此起可帮助用户设计出最佳的滤波器，并且，软件还可通过对比带宽、电流、成本和其他参数对增益带宽进行评估，为用户的设计选择最佳的运算放大器。
设计数字低通滤波器(用matlab实现)
DSP 设计滤波器报告姓名：张胜男 班级：07级电信（1）班 学号：078319120一·低通滤波器的设计（一）实验目的：掌握IIR 数字低通滤波器的设计方法。
（二）实验原理：1、滤波器的分类滤波器分两大类：经典滤波器和现代滤波器。
经典滤波器是假定输入信号)(n x 中的有用成分和希望取出的成分各自占有不同的频带。
这样，当)(n x 通过一个线性系统（即滤波器）后可讲欲去除的成分有效的去除。
现代滤波器理论研究的主要内容是从含有噪声的数据记录（又称时间序列）中估计出信号的某些特征或信号本身。
经典滤波器分为低通、高通、带通、带阻滤波器。
每一种又有模拟滤波器（AF ）和数字滤波器（DF ）。
对数字滤波器，又有IIR 滤波器和FIR 滤波器。
IIR DF 的转移函数是：∑∑=-=-+==Nk k k M r r rz a z b z X z Y z H 101)()()(FIR DF 的转移函数是：∑-=-=10)()(N n nz n h z HFIR 滤波器可以对给定的频率特性直接进行设计，而IIR 滤波器目前最通用的方法是利用已经很成熟的模拟滤波器的设计方法进行设计。
2、滤波器的技术要求低通滤波器：p ω：通带截止频率（又称通带上限频率） s ω：阻带下限截止频率 p α：通带允许的最大衰减s α：阻带允许的最小衰减 （p α，s α的单位dB ） p Ω：通带上限角频率 s Ω：阻带下限角频率 （s p p T ω=Ω，s s s T ω=Ω）即 C p p F ωπ2=Ω C s s F ωπ2=Ω3、IIR 数字滤波器的设计步骤：1）按一定规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标。
2）根据转换后的技术指标设计模拟低通滤波器)(s G ；3）再按一定的规则将)(s G 转换成)(z H 。
4）若是高通、带通或带阻数字滤波器则将它们的技术指标先转化为低通模拟滤波器的技术指标，然后按上述步骤2）设计出低通)(s G ，再将)(s G 转换为所需的)(z H 。
课程设计---IIR数字低通滤波器的MATLAB实现
IIR数字低通滤波器的MATLAB实现目录错误！未找到引用源。
一、引言 (4)（一）设计目的 (4)（二）设计思想.........................................错误！未定义书签。
二、IIR数字滤波器的MATLAB实现三、实验步骤 ..................................................错误！未定义书签。
四、实验结果及截图 ......................................错误！未定义书签。
结果.............................................................错误！未定义书签。
参考文献.....................................................错误！未定义书签。
附录.............................................................错误！未定义书签。
错误！未找到引用源。
一、引言设计目的（1）熟悉Matlab在信号处理方面的应用（2）深刻了解频域上的频谱特性。
（3）了解数字滤波器的设计。
（4）学会比较滤波前后语音信号的差异。
（5）利用GUI工具做一个简洁大方的图形用户界面。
二、IIR数字低通滤波器的MATLAB实现（一）IIR数字滤波器的传递函数及特点设IIR滤波器的输入序列为x(n)，则IIR滤波器的输入序列x(n)与输出序列y(n)之间的关系可以用下面的方程式表示：)()()(1j n x a i n x b n y Nj j M i i -+-=∑∑==其中，j a 和i b 是滤波器的系数，其中j a 中至少有一个非零。
与之相对应的差分方程为：NN MM z a z a z b z b b z X z Y Z H ------++==....1....)()()(11110由传递函数可以发现无限常单位冲激响应滤波器有如下特点： （1） 单位冲激响应h(n)是无限长的。
低通滤波器matlab课程设计
低通滤波器matlab课程设计一、课程目标知识目标：1. 理解低通滤波器的原理和数学模型；2. 掌握使用MATLAB进行低通滤波器的设计与实现；3. 学会分析低通滤波器的频率特性及其在实际应用中的作用。
技能目标：1. 能够运用MATLAB软件设计不同类型的低通滤波器；2. 能够通过调整滤波器参数，优化滤波效果；3. 能够利用所学的知识解决实际信号处理问题。
情感态度价值观目标：1. 培养学生对信号处理领域的兴趣，激发学习热情；2. 培养学生严谨的科学态度，注重实验数据的真实性；3. 培养学生的团队协作意识，提高沟通与交流能力。
课程性质：本课程为电子信息工程及相关专业的高年级学生设计，结合了理论知识与实践操作。
课程以低通滤波器的设计为主线，通过MATLAB软件的运用，使学生更好地理解信号处理技术。
学生特点：学生已具备一定的电路基础和MATLAB编程能力，具有较强的学习能力和动手实践能力。
教学要求：课程要求学生在掌握低通滤波器理论知识的基础上，通过MATLAB 软件实现滤波器的设计与优化。
教学过程中注重培养学生的实际操作能力和问题解决能力，将理论知识与实际应用紧密结合。
通过课程学习，使学生能够达到上述课程目标，为后续相关课程打下坚实基础。
二、教学内容1. 理论知识：- 低通滤波器的基本概念与分类；- 低通滤波器的传递函数与频率特性分析；- MATLAB环境下滤波器设计的基本原理。
2. 实践操作：- 使用MATLAB设计不同类型的低通滤波器（如巴特沃斯、切比雪夫等）；- 分析和调整滤波器参数，观察滤波效果的变化；- 对实际信号进行处理，验证滤波器的性能。
3. 教学大纲：- 第一周：低通滤波器的基本概念与分类，介绍相关教材章节；- 第二周：低通滤波器的传递函数与频率特性分析，结合教材相关内容；- 第三周：MATLAB环境下滤波器设计原理，结合教材实例进行讲解；- 第四周：实践操作，指导学生使用MATLAB设计低通滤波器，并进行性能分析；- 第五周：课程总结与作业布置，巩固所学知识。
(完整word版)用双线性变换法设计IIR数字低通滤波器课程设计
V=课程设计报告书姓名:班级：学号:时间:设计题目用双线性变换法设计IIR数字低通滤波器设计要求1. 通过实验加深对双线性变换法设计IIR滤波器基本方法的了解.2. 了解MATLAB有关双线性变换法的子函数。
3．掌握用双线性变换法设计数字低通滤波器的方法。
本次课程设计是采用双线性变换法基于MATLAB设计一个IIR数字低通滤波器, 其中要求通带截止频率为ωp=0.25π;通带最大衰减Rp=1dB;阻带最小衰减As=15dB;阻带截止频率ωs=0.4π;滤波器采样频率Fs=100Hz.设计过程摘要: 根据IIR滤波器的特点, 在MATLAB坏境下用双线性变换法设计IIR数字滤波器。
利用MATLAB设计滤波器, 可以随时对比设计要求和滤波器特性调整参数, 直观简便, 极大的减轻了工作量, 有利于滤波器设计的最优化。
1.关键词：双线性变换法 , 数字滤波器 , MATLAB , IIR2.设计原理与步骤1.1设计原理滤波器的种类很多, 从功能上可分为低通、高通、带通和带阻滤波器, 每一种又有模拟滤波器和数字滤波器两种形式。
如果滤波器的输人和输出都是离散时间信号, 则该滤波器的冲击响应也必然是离散的, 这种滤波器称之为数字滤波器。
数字滤波器是一种用来过滤时间离散信号的数字系统, 通过对抽样数据进行数学处理来达到频域滤波的目的。
数字滤波器也是具有一定传输选择特性的数字信号处理装置, 其输入、输出均为数字信号, 实质上是一个由有限精度算法实现的线性时不变离散系统。
IIR数字滤波器采用递归型结构, 即结构上带有反馈环路。
IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成, 可以组合成直接型、正准型、级联型、并联型四种结构形式, 都具有反馈回路。
数字滤波器根据其冲激响应函数的时域特性, 可分为两种, 即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。
IIR 数字滤波器的特征是, 具有无限持续时间冲激响应, 需要用递归模型来实现, 其差分方程为：（1-1）（1-2）设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数H(z), 使其频率响应H(z)满足所希望得到的频域指标, 即符合给定的通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、题目：巴特沃斯数字低通滤波器二、要求：利用脉冲响应不变法设计巴特沃斯数字低通滤波器，通带截止频率100HZ,采样频率1000HZ ，通带最大衰减为0.5HZ ，阻带最小衰减为10HZ ，画出幅频、相频相应相应曲线。
并假设一个信号x(t)=sin(2*pi*f1*t)+sin(2*pi*f2*t),其中f1=50HZ,f2=200HZ 。
用此信号验证滤波器设计的正确性。
三、设计原理1、数字滤波器的基本概念所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤波器除某些频率成分的数字器件或程序，因此，数字滤波的概念和模拟滤波相同，只是的形式和实现滤波方法不同。
正因为数字滤波通过数值运算实现滤波，所以数字滤波处理精度高、稳定、体积小、质量轻、灵活、不存在阻抗匹配问题，可以实验模拟滤波器无法实现的特殊滤波功能。
如果要处理的是模拟信号，可通过A\DC 和D\AC,在信号形式上进行匹配转换，同样可以使用数字滤波器对模拟信号进行滤波。
2、数字滤波器的工作原理数字滤波器是一个离散时间系统，输入x(n)是一个时间序列，输出y(n)也是一个时间序列。
如数字滤波器的系统函数为H(Z),其脉冲响应为h(n),则在时间域内存在下列关系y(n)=x(n) h(n)在Z 域内，输入输出存在下列关系Y(Z)=H(Z)X(Z)式中，X(Z),Y(Z)分别为输入x(n)和输出y(n)的Z 变换。
同样在频率域内，输入和输出存在下列关系Y(jw)=X(jw)H(jw)式中，H(jw)为数字滤波器的频率特性，X(jw)和Y(jw)分别为x(n)和y(n)的频谱。
w 为数字角频率，单位rad 。
通常设计H(jw)在某些频段的响应值为1，在某些频段的响应为0.X(jw)和H(jw)的乘积在频率响应为1的那些频段的值仍为X(jw)，即在这些频段的振幅可以无阻碍地通过滤波器，这些频带为通带。
X(jw)和H(jw)的乘积在频段响应为0的那些频段的值不管X(jw)大小如何均为零，即在这些频段里的振幅不能通过滤波器，这些频带称为阻带。
一个合适的数字滤波器系统函数H(Z)可以根据需要输入x(n)的频率特性，经数字滤波器处理后的信号y(n)保留信号x(n)中的有用频率成分，去除无用频率成分。
3、巴特沃斯滤波器设计原理(1) 基本性质巴特沃斯滤波器以巴特沃斯函数来近似滤波器的系统函数。
巴特沃斯滤波器是根据幅频特性在通频带内具有最平坦特性定义的滤波器。
巴特沃思滤波器的低通模平方函数表示1()ΩΩ+=Ωc N /22a 11)(j H N=1,2，…… （2-6）下面归纳了巴特沃斯滤波器的主要特征a 对所有的N ，()1a j H 20=Ω=Ω。
b 对所有的N ，()707.0a j 2c =ΩΩH =Ω即()dB 3a lg 20j H c =Ω=ΩΩ c ()Ωj H a 2是Ω的单调下降函数。
d ()Ωj H a 2随着阶次N 的增大而更接近于理想低通滤波器。
如下图2所示，可以看出滤波器的幅频特性随着滤波器阶次N 的增加而变得越来越好，在截止频率Ωc 处的函数值始终为1/2的情况下，通带内有更多的频带区的值接近于1；在阻带内更迅速的趋近于零。
图2 巴特沃思低通滤波平方幅频特性函数（2）系统函数设巴特沃斯的系统函数为H a （s ），则：（3）设计过程巴特沃思低通滤波技术指标关系式为a p >-20log|H a (j Ω)|，Ω<ΩPa s <-20log|H a (j Ω)|，Ω>Ωs其中：Ωp 为通带边界频率，Ωs 为阻带边界频率。
代入式1.4.1可得:经过化简整理可得：取满足上式的最小整数N作为滤波器的阶数。
再将N代入可得：或查表求得归一化传输函数H(s),令s/Ωc代替归一化原型滤波器系统函数中的s,即得到实际滤波器传输函数。
4、脉冲响应不变法所谓脉冲响应不变法就是数字滤波器的脉冲响应序列h(n)等于模拟滤波器的响应ha(t)的采样值，即h(n)=ha(t)|t=nT=ha(nT)式中，T为采样周期。
因此数字滤波器的系统函数H(Z)可由下式求得H(z)=Z[h(n)]=Z[ha(nT)]Z[-]表示[-]的内容进行变换，变换的内容请参考相应的数字信号处理材料。
如果已经获得了满足性能指标的模拟滤波器的传递函数Ha(s) ,求与之对应的数字滤波器的传递函数H(z)的方法是：（1）、求模拟滤波器的单位脉冲响应ha(t)。
式中，L[Ha(s)]表示对Ha(s)的Laplace.逆变换。
Laplace变换内容请参考高等数学的积分变换或信号处理教材。
（2）、求模拟滤波器单位冲激响应ha(t)的采样值，即数字滤波器冲激响应序列h(n)。
（3）、对数字滤波器的冲激h(n)响应进行z变换，得到传递函数H(z)。
由上述方法推论出更直接地由模拟滤波器系统函数Ha(s)求出数字滤波器系统函数H(z)的步骤是：（1）利用部分分式展开将模拟滤波器的传递函数H(z)展开成Ha(s)= Rk\(S-Pk)在MA TLAB中这步可通过residue函数实现若调用residue函数的形式为[b,a]=residue(R,P,K)形式。
若为[R,P,K]=residue(a,b)则为上面调用形式的反过程。
(2)将模拟极点Pk变换为数字极点e^pkT即得到数字系统的传递函数H(z)= Rk\(1-e^pkT*z*(-1))式中T为采样间隔。
（3）将上式转换为传递函数形式，可采用[R,P,K]=residue(b,a)。
对于上面的步骤，中已经提供了冲激响应不变法设计数字滤波器的函数，调用格式为[bz,az]=impinvar(b,a[ ,Fs],Fp)式中，b,a为模拟滤波器分子和分母多项式系数向量；Fs为采样频率（所滤波数据），单位Hz，缺省时为1Hz，为预畸变频率（prewarped frequency），是一个“匹配”频率，在该频率上，频率响应在变换前后和模拟频率可精确匹配。
一般设计中不考虑。
bz,az分别为数字滤波器分子和分母多项式系数向量。
5、实验所用MATLAB函数说数。
（1）[N,wc]=buttord(wp,ws,RP,As,’s’)该格式用于计算巴特沃斯模拟滤波器的阶数N和3db截止频率wc。
Wp、ws和wc是实际模拟角频率(rad\s)。
Rp和As为通带最大衰减和最小衰减。
（2）[Z,P,k]=buttap(N)该格式用于计算N阶巴特沃斯归一化模拟低通原型滤波器系统函数的零、极点和增益因子，返回长度为N的列向量Z和P，分别给出N个零点和极点的位置，K表示滤波器增益。
（3）Y=filter(b,a,x)式中b表示系统传递函数的分子多项式的系数矩阵；a表示系统传递函数的分母多项式的系数矩阵；x表示输入序列；filter表示输出序列。
IIR函数实现的直接形式。
(4) [b，a]=butter（N，wc，‘ftype’）计算N阶巴特沃斯数字滤波器系统函数分子、分母多项式的系数向量b、a。
说明：调用参数N和wc分别为巴特沃斯数字滤波器的阶数和3dB截止频率的归一化值，一般是调用buttord格式（1）计算N和wc。
系数b、a是按照z-1的升幂排列。
(5) [B，A]=butter（N，Ωc，‘ftype’，‘s’）计算巴特沃斯模拟滤波器系统函数的分子、分母多项式系数向量。
说明：调用参数N和Ωc分别为巴特沃斯模拟滤波器的阶数和3dB截止频率(实际角频率)，可调用buttord（2）格式计算N和Ωc。
系数B、A按s的正降幂排列。
tfype为滤波器的类型：◇ftype=high时，高通；Ωc只有1个值。
◇ftype=stop时，带阻；Ωc=[Ωcl,Ωcu]，分别为带阻滤波器的通带3dB下截止频率和上截止频率。
◇ ftype缺省时：若Ωc只有1个值，则默认为低通；若Ωc有2个值，则默认为带通；其通带频率区间Ωcl<Ω <Ωcu。
(6)[H,w]=freqz(b,a,N)b和a分别为离散系统的系统函数分子、分母多项式的系数向量，返回量H则包含了离散系统频响在 0～pi范围内N个频率等分点的值（其中N为正整数），w则包含了范围内N 个频率等分点。
调用默认的N时，其值是512。
可以先调用freqz()函数计算系统的频率响应，然后利用abs()和angle()函数及plot()函数，绘制出系统的频响曲线。
(7)lp2lp函数 [bt,at]=lp2lp(b,a,w0)该函数用于实现由低通模拟原型滤波器至低通滤波器的频率变换，可以用传递函数和状态空间进行转换，但无论哪种形式，其输入必须是模拟滤波器原型。
（8）[bz,az]=impinvar（b,a,fs)把具有[b,a]模拟滤波器传递函数模型转换为采样频率为fs的数字滤波器的传递函数模型[bz,az]，如果在函数中没有确定频率fs时，函数默认为1Hz.四、设计思路设定信号↓模拟低通滤波器原型→频率变换→模拟离散化→IIR数字滤波器→输出信号五、设计内容1.MATLAB程序设计Wp=2*pi*100; Ws=2*pi*150; %滤波器截止频率Rp=0.5; Rs=10; %通带最大衰减和阻带最小衰减Fs=1000; %采样频率Nn=128; %调用freqz所用的频率点数[N,Wc]=buttord(Wp,Ws,Rp,Rs,'s'); %模拟滤波器的最小阶数[z,p,k]=buttap(N); %设计模拟低通原型Butterworth滤波器[Bap,Aap]=zp2tf(z,p,k); %将零点极点增益形式转换为传递函数形式[b,a]=lp2lp(Bap,Aap,Wc) %进行频率转换[bz,az]=impinvar(b,a,Fs); %运用脉冲响应不变法得到数字滤波器的传递函数figure(1)[H,W]=freqz(bz,az,Nn,Fs); %绘制数字滤波器的幅值特性和相频特性subplot(2,1,1)plot(W,20*log10(abs(H)));xlabel('频率');ylabel('幅度');grid on;subplot(2,1,2);plot(W,180\pi*unwrap(angle(H)));xlabel('频率');ylabel('幅度');grid on;figure(2)f1=50; f2=200; %输入信号的频率N=100; %数据长度dt=1\Fs; n=0:N-1; t=n*dt; %采样间隔和时间序列x=sin(2*pi*f1*t)+sin(2*pi*f2*t); %滤波器输入信号subplot(2,1,1); %绘制输入信号plot(t,x);title('输入信号'); %用filter函数对输入信号滤波y1=filter(bz,az,x); subplot(2,1,2); %绘出输出波形xlabel('时间')title('输出信号');2.实验结果分析实验得到的两幅图如上所示，在第一幅图中，小于100处衰减小于3Hz，而在大于120Hz 处衰减大于15dB，满足滤波器设计指标。