双通带切比雪夫带通滤波器的设计matlab

1、中北大学课程设计说明书学生姓名：学号：学 院：信息与通信工程学院专 业：电子信息工程题 目：信息处理信息实践：信号的产牛与带通滤波器的设计指导教师：职称：2016年1月4日中北大学课程设计任务书2015/2016 学年第 一 学期学院：信息与通信工程学院专业：电子信息工程学生姓名：学生姓名：学号：学生姓名：学号：_课程设计题目：信程处理信息实践：信号的产牛与带通滤波器的设计起迄日期：2016年1月4日2016年1月22日课程设计地点：学院楼201实验室指导教师：王玉下达任务书日期：2016 年1月4日课程设计任务书1 .设计目的：(1)通过实践对滤波器的知识有更好的认识和巩固;(2)熟悉滤波器

2、的设计方法及作用；(3) 了解付氏变换的应用情况。2 .设计内容和要求(包括原始数据、技术参数、条件、设计要求等) ：(1)产生三种不同频率的正弦信号并叠加，通过 fft变换得到频域图；(2)设计在200-300hz, 500-600hz的双通带带通滤波器，并进行频谱分析；(3)对合成信号进行滤波，对滤波前后信号进行对比；(4)要求3位同学完成。3 .设计工作任务及工作量的要求包括课程设计计算说明书(论文)、图纸、实物样品等：(1)要求设计组的每个成员都要了解设计的要求和思路；(2) matla嗷据处理部分要求有正确的运行结果及结果分析；(3)对带通滤波器要求每位同学有自己的理解；(4)每位同

3、学针对上述内容撰写设计说明书(每人 1份)。课程设计任务书4 .主要参考文献：1桂志国，杨民，陈友兴，郝利华.数字信号处理原理及应用m.北京：国防工业 出版社，2012,12宋宇飞，潘子宇，魏蛆.数字信号处理实验与学习指导m.北京：清华大学出版 社，2012.83赵谦.通信系统中matlan础与仿真应用m.西安：西安电子科技大学出版社， 2010.34王彬，于丹，汪洋.matla吸字信号处理m.北京：机械工业出版社，2010.55刘卫国.matlab®序设计与应用（第二版）m.北京：高等教育出版社，20065 .设计成果形式及要求:课程设计说明书程序运行结果6 .工作计划及进度：20

4、16 年1月4日1月6日：查资料，了解基于 usb总线或pci总线a/d卡的通用结构以 及a/d采集卡的应用；1月7日1月20日：在老师指导下完成课程设计内容；1月20日1月22日：撰写课程设计说明书、答辩。负责人审查意见:签字：年 月 日一基于usb总线a/d卡的报告11 .usb总线介绍12 .usb接口电路设计13 .usb接口的数据采集系统的设计实现 14 .a/d转换电路25 .a/d 采集卡的应用 22 设计内容 33 设计原理 3四matla瞅件介绍45 设计过程 51. 信号的产生52. 对连续输入信号进行采样，进行频谱分析 73. 切比雪夫带通滤波器84. 对信号进行滤波10

5、6 结果分析 127 心得体会及总结13参考文献 13附录：总程序代码15一基于usb总线a/d卡的报告1 .usb总线介绍通用串行总线标准usb是1995年微软、康柏、dec旧m等公司为解决传统 总线不足的问题而推出的一种新型通信标准。该标准具有高带宽、安装方便、易 扩展诸多优点，逐渐成为现代数据传输的主要手段。基于usbt据采集卡充拥有 usb总线的上述优点，有效优化解决传统数据采集卡的安装麻烦、扩展不便等弊 端。usb总线为通用串行总线，us眼口位于ps/2接口和用并口之间，允许外设 在开机状态下热插拔，最多可用接下来 127个外设，传输速率可达480mb/s, p 它可以向低压设备提供

6、5伏电源，同时可以减少pc机i/o接口数量。usb是基 于通用连接技术，实现外设的简单快速连接，达到方便用户、降低成本、扩展 pc连接外设范围的目的。2 .usb接口电路设计图1电路设计如图1所示，r3是上拉电阻器，它可使 usb口的d+®上拉到ds2490s勺vb 端，表示usb主机系统是高速设备，同时这个上拉电阻器告诉主机有 us股备插 入。该上拉电阻器的设置对适配器的影响很大，它的负载值和 1-wire网络的总 长决定1-wire总线电压上升到5 v的速度。经过实验测试选择 r3的阻值为27 q±lo% r1、r2为usbr据线保护电阻器。l、l2具有禁止高频干扰并且

7、减弱 emi辐射的功能。lf33cv为3.3 v电压稳压器，与周围元件 c1、c2组成强上拉 部分，给eeprom温度传感器等器件提供额外的电源。3 .usb接口的数据采集系统的设计实现数据采集系统使用采集卡进行数据采集，然后经过a/d转换器供计算机加工 处理。基于usbr口的数据采集与频谱分析系统本系统结构由硬件部分和软件部 分组成，硬件部分主要有计算机、i/o接口设备。计算机作为硬件平台的核心可 采用台式机，系统采用的i/o设备为a/d数据采集卡，该采集卡是一种基于 usb 总线数据采集产品，可与带usb8口的各种台式计算机、笔记本电脑、工控机连 接构成高性能的数据采集测量系统.整个系统主

8、要由4部分组成：usb接口芯片 及外围电路、控制电路、数据缓冲电路和a/d转换电路。usb接口芯片选择了cypress公司的ez-usb2131q该芯片内嵌8051控制器，因此整个系统以ez-usb 控制器为核心，由ez-usb®空制电品&实现对a/ d转换电路和数据缓冲电路的控 制，模拟信号转换后的数据送入数据缓冲器，当数据缓冲器存满之后，通知 ez-usb空制器，由主机取出数据。整个系统框图如图2所示。图2系统框图4 .a/d转换电路a/d 转换芯片采用了 maxi泌司的max122芯片是12b的高速的a/d转换器。 在完全转换模式下，他的转换时间可以达到2.6&

9、mu;s,采样率为333ks/s°max122 有5种工作模式，在数据采集系统中，采用了模式2即连续转换模式。在这种模 式下，每次转换需要1314个时钟脉冲节拍，转换可以不间断地进行，但是需 要提供开始转换使能信号，并且要保证使能信号和时钟信号同步， 读信号和片选 始终处于有效状态。数据输出使能信号一直有效，在转换结束时产生新的数据。5 .a/d采集卡的应用数据采集就是把来自各种传感器的信号数据实时地、准确地测量或汇集起来， 用计算机进行实时处理或记录存储，实时完成测试和控制功能。数据采集系统结 构通过微机的标准接口连接各种功能模块、仪器仪表和传感器，组成测量系统。数据采集监测技术已

10、成为日益重要的检测技术， 广泛应用于需要监控湿度、温度、 压力等工农业中。数据采集过程是工业控制等系统中的重要过程，通常主要采用一些功能独立的单片机系统来采集， 是测控系统不可缺少的一部分， 因为数据采集的性能特点直接影响到整个系统的稳定性。 在生产过程中这一系统可对现场的工艺参数进行监视、采集以及记录，从而提高产品质量，提供信息和手段，降低生产成本。而在科学研究中应用该系统可以获得大量动态然后研究瞬间物理过程，也成为获取科学奥秘重要手段之一。总之，无论在哪个领域，数据采集系统 应用的越及时，工作效率就愈高，获得经济利益就愈大，这是毋庸置疑的。二 设计内容产生三种不同频率的正弦信号并叠加，通过

11、fft变换得到频域图；并设计在200-300hz 的带通滤波器，并进行频谱分析；再对合成信号进行滤波，对滤波前 后信号进行对比。三 设计原理1. 利用matla歆件产生三个连续信号，并对三个连续信号进行合成。2. 对信号进行抽样，进行频谱分析。( 1)时域采样(奈奎斯特采样)定理：为了避免产生混叠现象，能从抽样信号无失真地恢复出原信号，抽样频率必须大于或等于信号频谱最高频率的两倍。本设计中信号最高频率是500hz,抽样频率采用2000hz( 2)频谱分析：频谱分析是指对信号进行频域谱的分析，观察其频域的各个分量的功率大小， 其理论基础是傅立叶变换， 现在一般采用数字的方法， 也就 是将时域信号

12、数字化后做fft,可以得到频域的波形。3. 带通滤波器滤波的设计现代生活中， 为了滤除谐波干扰， 获得所需要的高精度的模拟信号， 经常要用到滤波器对信号进行滤波。典型的模拟滤波器有巴特沃斯(butterworth) 滤波器、 切比雪夫 (chebyshev) 滤波器和椭圆 (ellipse) 滤波器等。 其中， 巴特沃斯滤波器又叫最平坦响应滤波器，顾名思义，它的响应最为平坦，通带内没有波纹，其频率响应在通带和阻带中都是单调的， 且在靠近零频处最平坦， 而在趋向阻带时衰减单调增大， 巴特沃斯响应能够最大化滤波器的通带平坦度。 该响应非常平坦，非常接近dc信号，然后慢慢衰减至截止频率点为-3db,

13、最终逼近-20ndb/decade的衰减率，其中n为滤波器的阶数。切比雪夫滤波器在过渡带比巴 特沃斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。切比雪夫滤波器 和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。本次课程设计将使用切比雪夫1型带通滤波器对信号进行滤波。滤波器的结构框 图如下图3所示：图3滤波器的结构框图四matla琳件介绍matla欧件是由美国mathworks公司推出的用于数值计算和图形处理的科 学计算系统环境。在matla解境下，用户可以集成地进行程序设计、 数值计算、 图形绘制、输入输出、文件管理等各项操作。matlab!供了一个人机交互的数学系统

14、环境，该系统的基本数据结构是矩阵，在生成矩陈对象时，不要求作明确 的维数说明。与利用 c语言或fortran言作数值计算的程序设计相比，利用 matlabt以节省大量白编程时间。matlab3其提供计算、图形、文字处理的统 一环境而深受欢迎。mtala系统由五个主要部分组成，下面分别加以介绍。(1)matalb®言是高层次的矩阵/数组语言。具有条件控制、函数调用、数 据结构、输入输出、面向对象等程序语言特性。利用它既可以进行小规模编程， 完成算法设计和算法实验的基本任务， 也可以进行大规模编程，开发复杂的应用 程序。(2)matlab工作环境这是对matlabg供给用户使用的管理功能

15、的总称。包 括管理工作空间中的变量据输入输出的方式和方法，以及开发、调试、管理 m 文件的各种工具。(3)matlab图形系统的基石包括完成 2d和3d数据图示、图像处理、动画 生成、图形显示等功能的高层 matla脸令，也包括用户对图形图像等对象进行 特性控制的低层matlabt令，以及开发gui应用程序的各种工具。（4）matlab数学函数库这是对matla改用的各种数学算法的总称。包括各 种初等函数的算法，也包括矩阵运算、矩阵分析等高层次数学算法。（5）matlab应用程序接口（api）这是matla的用户提供的一个函数库，使得 用户能够在matlam境中使用c程序或fortran程序，

16、包括从matla呻调用于 程序（动态链接），读写matt件的功能。可以看出 matla是一个功能十分强大 的系统，是集数值计算、图形管理、程序开发为一体的环境。五 设计过程1.信号的产生首先产生一个连续输入信号，包含低频（f=30hz ） ，中频（f=250hz）高频（f=500hz）分量，并实现三个信号的合成。1）程序代码f1=30;f2=250;f3=500;t=(1:100)/2000;x1=sin(2*pi*t*f1);figure(1);subplot(2,1,1);plot(x1);xlabel('t');ylabel('x1(t)');title(

17、'低频信号');grid;x2=sin(2*pi*t*f2);subplot(2,1,2);plot(x2);xlabel('t');ylabel('x2(t)');title(' 中频信号 ');grid;x3=sin(2*pi*t*f3);figure(2);subplot(2,1,1);plot(x3);xlabel('t');ylabel('x3(t)');title(' 高频信号 ');%产生三种频率%绘制x1（t） 的图形%绘制x2（t） 的图形%绘制x3（t） 的图形2

18、4grid；x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3);%实现三种信号的叠加subplot(2,1,2);plot(x);%绘制 x(t)的图形xlabel('t');ylabel('x(t)');title('合成信号');grid;(2)程序运行结果如图4、图5所示低频信号 10 5 营0081010203040 扣 6470 so 91j 1mt中频信号 1 0 5哥00 5图4低频中、频信号波形图高频信方 10 5s 0-0 5=101020304050 go 70 bq 901004合成

19、t信号0102030405。60706090 idoi图5高频、合成信号波形图(3)结果分析图4、5所示内容为按照设计要求作出的三种频率分量以及合成信号。种频率分量为正弦信号，但合成之后变成非正弦信号。2.对连续输入信号进行采样，进行频谱分析(1)程序代码%设定抽样频率为2000hz%直接求fft%频率横坐标的设定n=1:100;t=n/2000x=fft(x,512);w=(0:255)/256*1000;x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3);%叠加信号的抽样figure(3);stem(x);%绘制 x(n)的图形xlabel(

20、9;n');ylabel('x(n)');title('数字信号');grid；figure(4);plot(w,abs(x(1:256);% 绘制频谱图xlabel('hz');ylabel('频率响应幅度');title('滤波前的频谱图');grid;(2)程序运行结果如图6、图7所示图6连续信号抽样结果波形图滤波前的频谱图d o d 4 3 2 蚂里保整榔懦图7滤波前的频谱分析波形图00do 2304000670o00(3)结果分析图6、图7是对包含低频、中频、高频的连续信号抽样，抽样点数取为 51

21、2 从图中看出抽样的模型是对图6中的连续信号波形的采样，并且满足时域采样定 理。matlab优点是直接可以采用fft!行快速傅里叶变换。本设计就是直接利 用fft®行频谱分析，得至i了 20hz 250hz 500hzi勺频谱。3.iir切比雪夫带通滤波器(1)程序代码fsl=150;fsu=350;%第一个阻带截止频率fs1=450;fs2=650;%第二个阻带截止频率fpl=200;fpu=300;fp1=500;fp2=600;fs=2000;%第一个通带截止频率%第二个通带截止频率%抽样频率ws=2*fsl/fs,2*fsu/fs;%数字滤波器的归一化阻带截止频率wp=2*f

22、pl/fs,2*fpu/fs;%数字滤波器的归一化通带截止频率ws1=2*fs1/fs,2*fs2/fs;wp1=2*fp1/fs,2*fp2/fs;rs=20;rp=0.5;%衰减设置n,wpo=cheb1ord(wp,ws,rp,rs);%计算阶数和通带截止频率m,wp1=cheb1ord(wp1,ws1,rp,rs);b,a=cheby1(n,rp,wpo);c,d=cheby1(m,rp,wp1);h1,w1=freqz(b,a);h2,w2=freqz(c,d);h=h1+h2;w=w1+w2;%计算第一个通带分子分母的多项式系数%计算第二个通带分子分母的多项式系数 %求解第一个通带

23、数字滤波器的频率响应 %求解第二个通带数字滤波器的频率响应 %实现两个通带滤波器figure(5)plot(w*2000/(2*2*pi),abs(h);grid on;xlabel('频率/hz');ylabel('幅度');title('切比雪夫1型带通滤波器的幅频特性');set(gca,'ytickmode','manual','ytick',0 0.2 0.4 0.6 0.707 1);figure(6)plot(w/(2*pi),angle(h);grid on;xlabel('

24、频率/hz');ylabel('angle');title('切比雪夫1型带通滤波器的相频特性');(2)程序运行结果如图8、图9所示切比雪夫1型带逋漉波器的幅甄特性1111111x _一11例 0.707 0 &0.40 20hwvf)11.1j+一1- - - « i 一)10020030。400500 6cc 700600900 1000频率用工图8切比雪夫1型双通带带通滤波器的幅频特性切比雪夫1型带通谑波器的相频特性3ti i i i i ii i i00 10.20 30.40.5fl.60.70.80.9频率jhz图9切比雪

25、夫1型双通带带通滤波器的相频特性(3)分析设定iir切比雪夫1型带通数字滤波器指标：通带范围为：200-300hz和 500-600hz,第一个通带阻带上限为：350hz,阻带下限为150hz第二个通带阻带 上限为：650hz,阻带下限为450hz通带最大衰减ap=0.5db,阻带最小衰减为 as=20dr 采样频率为 fsa=2000hz。其中，n, wpo=cheblord(wp, ws, rp, as) 该格式用于计算切比雪夫1滤波器的阶数n和通带中心频率wpd wp和ws分别 为滤波器的通带边界频率和阻带截止频率的归一化值，要求0<=wp和ws<=1,表示模拟频率fs/2。

26、rp和as分别为通带最大衰减和阻带最小衰减。ws<wpm,为高通滤波器。b, a=cheby1(n, rp , wpo),该格式用于计算切比雪夫1滤波 器系统函数的分子分母多项式系数向量 b和a调用参数n和wpo分别为切比雪 夫1滤波器的阶数n和通带中心频率。值得我们注意的是，切比雪夫滤波器分为 1型，2型滤波器，1型滤波器为通带内是等波纹的，阻带是单调递减的；2型滤波器为通带是单调递减的，而阻带是等波纹的。4.对信号进行滤波(1)程序代码y1=filter(b,a,x);y2=filter(c,d,x);y=y1+y2;figure piot(t,y); xlabel('时间(

27、t)'); ylabel('滤波后的信号y'); title('滤波后的信号');grid on; axis(0 0.05 -2 2);y=fft(y,512);%对滤出来的信号进行频谱分析w=(0:255)/256*1000;%频率横坐标的设定figure(8) plot(w,abs(y(1:256);% 绘制频谱图xlabel('频率/hz'); ylabel('频率响应幅度'); title('滤波后的频谱图'); grid;(2)程序运行结果如图10、图11所示21 5涯'波后的信号5 0

28、6 o 410it妲s血度照05 0.5 ,04 0.04 0.5 03 035(0 02间 。时 020.5 j o 0.01 5 o 0 0.o图10滤波后信号时域波形图1403530o 5 0 5 0 5 0 2 2 112qo 3300 460率00 600700300900 o图11滤波后信号频谱图六结果分析设计过程中，首先产生连续输入信号，包含低频（f=30hz）,中频（f=250hz）, 高频（f=500hz）分量，然后对其进行采样，利用傅里叶变换进行频谱分析，并 由带通滤波器的参数设计带通滤波器对信号进行滤波处理，对应带通滤波器的通 带范围是（ 200, 300）和（400,5

29、00）,从运行结果图中可以看出，经过带通滤波 器滤波后信号对应的频率为原信号中的中频分量（f=250hz）和高频分量（f=500hz）。由上述结果显示，在误差允许的范围内实验结果与理论结果相同。出现误差的原因：在设计滤波器的参数时并不是十分的准确，在不同计算机上运行 matlab时会有一定的偶然误差，从而导致实验误差的存在。并且采用切比雪夫滤波器的 阻带衰减特性则有所改善。巴特沃斯带通滤波器的幅频特性随。的增加而单调下 降，当n较小时，阻带幅频特性下降较慢，要想使其幅频特性接近理想带通滤波 器，就必须增加滤波器的阶数，这就将导致模拟滤波器使用的原件增多，线路趋于复杂。如果阶次一定，巴特沃兹滤波

30、器在通带内幅度特性是单调下降的，则在靠近截止频率处，幅度变化范围比较大，或者说，为了使通带内的衰减足够小， 需要的阶次（n）很高，采用切比雪夫多项式逼近可以克服这一缺点。因为切比雪 夫滤波器在通带范围内是等幅起伏的， 所以同样的通带衰减，它的阶数比巴特沃兹滤波器要小。七 心得体会及总结这次课程设计， 老师给我们的时间比较长， 在课设的前半段时间， 我们设计了三个正弦信号的叠加，fft频谱分析，巴特沃斯带通滤波器及滤出正确的波形。在设计带通滤波器时我们遇到了麻烦， 我们先设计了模拟低通滤波器， 再转换为带通滤波器， 最后转换为数字带通滤波器， 我们所设计的滤波器满足了设计要求，但怎么也滤不出来波

31、。后来经过查资料，我们发现matlabt个很大的优点，可以不通过低通滤波器来设计带通滤波器， 可以直接根据参数设置来求得带通滤波器。 这使得我们的设计从头开始， 直接设计巴特沃斯带通滤波器， 成功把我们所需要的波给滤出来了。 在实验课上， 经过王老师的指导， 如果我们想要做得更完美，可以再设计一种到两种带通滤波器。我们小组最后决定再设计一个fir带通滤波器， iir 切比雪夫带通滤波器。由于任务比较多，我们把任务进行分块，而我主要设计切比雪夫带通滤波器。在设计时，我先利用w女t把模拟频率转换为数字角频率， 这涉及到归一化问题， 经过我们查阅资料， 在 iir 和 fir 滤波器设计中，频率是用

32、 fs/2 进行归一化，如果是弧度的话，就是用 pi 进行归一化。至于为什么要进行归一化，因为滤波器有各种类型，包括低通，高通，带通，带阻等，而且就一种类型而言，中心频率，截止频率又不同，但是设计表格不能有那么多， 一般只给出归一化原型滤波器的设计数据。 在设计出切比雪夫带通滤波器后，我们又设计出双通带的切比雪夫带通滤波器，最后出现了极好的滤波图。所以这次课程设计，真的使我们更好的理解了这学期学的数字信号处理的知识，也弥补了课堂上的不足，让我们学会了简单的matlab勺使用。参考文献1 桂志国，杨民，陈友兴，郝利华. 数字信号处理原理及应用 m. 北京：国防工业出版社， 2012,12 宋宇飞

33、，潘子宇，魏峘. 数字信号处理实验与学习指导 m. 北京：清华大学出版社， 2012.83赵谦.通信系统中matlab!础与仿真应用m.西安：西安电子科技大学出版社， 2010.34王彬，于丹，汪洋.matla吸字信号处理m.北京：机械工业出版社，2010.55刘卫国.matlab®序设计与应用（第二版）m.北京：高等教育出版社,2006附录：总程序代码% 信号的产生f1=30;f2=250;f3=500;t=(1:100)/2000;x1=sin(2*pi*t*f1);figure(1);subplot(2,1,1);plot(x1);xlabel('t');yla

34、bel('x1(t)');title(' 低频信号 ');grid;x2=sin(2*pi*t*f2);subplot(2,1,2);plot(x2);xlabel('t');ylabel('x2(t)');title(' 中频信号 ');grid;x3=sin(2*pi*t*f3);figure(2);subplot(2,1,1);plot(x3);xlabel('t');ylabel('x3(t)');title(' 高频信号 ');grid;x=sin(2*pi

35、*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3);% 绘制 x(t) 的图形% 绘制 x2(t) 的图形% 绘制 x3(t) 的图形subplot(2,1,2);plot(x);% 绘制 x(t) 的图形xlabel('t');ylabel('x(t)');title(' 合成信号 ');grid;% 对连续输入信号进行采样，进行频谱分析n=1:100;t=n/2000x=fft(x,512);w=(0:255)/256*1000;x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3)

36、;figure(3);stem(x);% 绘制 x(n) 的图形xlabel('n');ylabel('x(n)');title(' 数字信号 ');grid;figure(4);plot(w,abs(x(1:256);% 绘制频谱图xlabel('hz');ylabel(' 频率响应幅度');title(' 频谱图 ');grid;%fir 带通滤波器(汉明窗)wls=0.15*pi;wlp=0.2*pi;% 参数设置% 求两个过渡带的较小者截止频率取通带阻带边界频率的均值whp=0.3*pi;w

37、hs=0.35*pi;delta_w=min(wlp-wls),(whs-whp);wc1=(wls+wlp)/2;wc2=(whp+whs)/2; % %hamming 窗n2=ceil(6.6*pi/delta_w);%根据hamming窗精确过渡带宽6.6 u/n计算窗宽hn2=fir1(n2-1,wc1,wc2/pi,hamming(n2);h2,w2=freqz(hn2,1);% 绘图figure(5) %建立图形窗口subplot(3,1,1);n=0:n2-1;stem(n,hn2,'.');%绘制 hamming 窗单位脉冲响应axis(0,n2-1,-0.4,

38、0.4);%确定显示范围xlabel('n');ylabel('h(n)');grid on;title('hamming 窗单位脉冲响应h(n)');subplot(3,1,2);plot(w2*2000/(2*pi),20*log10(abs(h2);%绘制 hamming 窗幅频响应axis(0,1000,-150,5);xlabel(' 归一化角频率 ');ylabel(' 幅度(单位：分贝) ');grid on;title('hamming窗幅频响应');subplot(3,1,3);p

39、lot(w2/pi,180/pi*unwrap(angle(h2);%绘制 hamming 窗相频响应xlabel(' 归一化角频率 ');ylabel(' 单位：度 ');grid on;title('hamming窗相频相应');grid;y=filter(hn2,1,x);figure(6);subplot(2,1,1);plot(y);xlabel('t');ylabel('x(t)');title(' 连续信号 ');grid;y=fft(y,512);w=(0:255)/256*1000

40、;subplot(2,1,2);plot(w,abs(y(1:256);xlabel('hz');ylabel(' 频率响应幅度');title(' 频谱图 ');grid;% 巴特沃斯带通滤波器fp=200 300;fs=150 350;ap=2;as=30;fsa=2000;wp=fp/fsa*2;ws=fs/fsa*2;n,wn=buttord(wp,ws,ap,as);b,a=butter(n,wn);h,w=freqz(b,a,512);figure(7);subplot(2,1,1);plot(w*2000/(2*pi),abs(h)

41、;xlabel('hz');ylabel(' 频率响应幅度');title(' 带通滤波器');grid;% 绘制频谱图% 绘制带通频谱图subplot(2,1,2);plot(w/pi,angle(h);xlabel('hz');ylabel('angel');title(' 相位特性 ');grid;y=filter(b,a,x);figure(8);subplot(2,1,1);plot(y);xlabel('t');ylabel('x(t)');title(' 连续