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

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）MATLAB 数据处理部分要求有正确的运行结果及结果分析；（3）对带通滤波器要求每位同学有自己的理解

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

4、.6工作计划及进度：2016 年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.AD 转换电路25.A/D 采集卡的应用2二 设计内容3三 设计原理3四 MATLAB 软件介绍4五 设计过程51.信号的产生52.对连续输入信号进行采

5、样，进行频谱分析7.3.切比雪夫带通滤波器84.对信号进行滤波10六 结果分析12七 心得体会及总结13参考文献13附录：总程序代码15.一 基于 USB 总线 A/D 卡的报告1.USB 总线介绍通用串行总线标准USB是 1995 年微软、康柏、 DEC、IBM 等公司为解决传统总线不足的问题而推出的一种新型通信标准。该标准具有高带宽、安装方便、易扩展诸多优点，逐渐成为现代数据传输的主要手段。基于USB 数据采集卡充拥有 USB 总线的上述优点，有效优化解决传统数据采集卡的安装麻烦、扩展不便等弊端。 USB 总线为通用串行总线， USB 接口位于 PS/2 接口和串并口之间，允许外设在开机状

6、态下热插拔，最多可串接下来127 个外设，传输速率可达480Mb/S ， P 它可以向低压设备提供5 伏电源，同时可以减少PC 机 I/O 接口数量。USB是基于通用连接技术，实现外设的简单快速连接，达到方便用户、降低成本、扩展 PC 连接外设范围的目的。2.USB 接口电路设计图 1 电路设计如图 1 所示， R3 是上拉电阻器，它可使USB 口的 D+ 端上拉到 DS2490S 的VB 端，表示 USB主机系统是高速设备， 同时这个上拉电阻器告诉主机有USB设备插入。该上拉电阻器的设置对适配器的影响很大，它的负载值和1-Wire 网络的总长决定 1-Wire 总线电压上升到5 V 的速度。

7、经过实验测试选择R3 的阻值为.27±lO%。 R1、 R2 为 USB 数据线保护电阻器。 L、 L2 具有禁止高频干扰并且减弱 EMI 辐射的功能。 LF33CV 为 3.3 V 电压稳压器， 与周围元件 C1、C2 组成强上拉部分，给 EEPROM或温度传感器等器件提供额外的电源。3.USB 接口的数据采集系统的设计实现数据采集系统使用采集卡进行数据采集 ,然后经过 A/D 转换器供计算机加工处理。基于 USB 接口的数据采集与频谱分析系统本系统结构由硬件部分和软件部分组成 ,硬件部分主要有计算机、 I/O 接口设备。计算机作为硬件平台的核心可采用台式机 ,系统采用的 I/O

8、设备为 A/D 数据采集卡 ,该采集卡是一种基于 USB总线数据采集产品 ,可与带 USB 接口的各种台式计算机、笔记本电脑、工控机连接构成高性能的数据采集测量系统 .整个系统主要由 4 部分组成： USB 接口芯片及外围电路、控制电路、数据缓冲电路和AD 转换电路。 USB 接口芯片选择了Cypress 公司的EZ-USB 2131Q ，该芯片内嵌8051 控制器，因此整个系统以EZ-USB 控制器为核心，由EZ-USB 经控制电路实现对AD 转换电路和数据缓冲电路的控制，模拟信号转换后的数据送入数据缓冲器，当数据缓冲器存满之后，通知 EZ-USB 控制器，由主机取出数据。整个系统框图如图2

9、 所示。图 2 系统框图4.AD 转换电路A/D 转换芯片采用了MAXIM 公司的 MAX122 ，芯片是 12b 的高速的 A/D 转换器。在完全转换模式下，他的转换时间可以达到26μs ，采样率为 333kS/s。.MAX122 有 5 种工作模式，在数据采集系统中，采用了模式2 即连续转换模式。在这种模式下，每次转换需要 1314 个时钟脉冲节拍， 转换可以不间断地进行，但是需要提供开始转换使能信号，并且要保证使能信号和时钟信号同步，读信号和片选始终处于有效状态。 数据输出使能信号一直有效， 在转换结束时产生新的数据。5.A/D 采集卡的应用数据采集就是把来自各种传感器的信

10、号数据实时地、准确地测量或汇集起来 ,用计算机进行实时处理或记录存储,实时完成测试和控制功能。数据采集系统结构通过微机的标准接口连接各种功能模块、仪器仪表和传感器,组成测量系统。数据采集监测技术已成为日益重要的检测技术，广泛应用于需要监控湿度、 温度、压力等工农业中。 数据采集过程是工业控制等系统中的重要过程，通常主要采用一些功能独立的单片机系统来采集，是测控系统不可缺少的一部分，因为数据采集的性能特点直接影响到整个系统的稳定性。在生产过程中这一系统可对现场的工艺参数进行监视、采集以及记录，从而提高产品质量，提供信息和手段，降低生产成本。而在科学研究中应用该系统可以获得大量动态然后研究瞬间物理

11、过程，也成为获取科学奥秘重要手段之一。总之，无论在哪个领域，数据采集系统应用的越及时，工作效率就愈高，获得经济利益就愈大，这是毋庸置疑的。二 设计内容产生三种不同频率的正弦信号并叠加，通过 FFT变换得到频域图； 并设计在200-300Hz 的带通滤波器， 并进行频谱分析； 再对合成信号进行滤波， 对滤波前后信号进行对比。.三 设计原理1.利用 MATLAB 软件产生三个连续信号，并对三个连续信号进行合成。2.对信号进行抽样，进行频谱分析。（1）时域采样（奈奎斯特采样）定理：为了避免产生混叠现象，能从抽样信号无失真地恢复出原信号，抽样频率必须大于或等于信号频谱最高频率的两倍。本设计中信号最高频

12、率是500Hz，抽样频率采用2000Hz。（2）频谱分析：频谱分析是指对信号进行频域谱的分析，观察其频域的各个分量的功率大小， 其理论基础是傅立叶变换， 现在一般采用数字的方法，也就是将时域信号数字化后做FFT，可以得到频域的波形。3.带通滤波器滤波的设计现代生活中， 为了滤除谐波干扰， 获得所需要的高精度的模拟信号，经常要用到滤波器对信号进行滤波。典型的模拟滤波器有巴特沃斯(Butterworth) 滤波器、切比雪夫 (Chebyshev)滤波器和椭圆 (Ellipse)滤波器等。其中，巴特沃斯滤波器又叫最平坦响应滤波器，顾名思义，它的响应最为平坦，通带内没有波纹，其频率响应在通带和阻带中都

13、是单调的，且在靠近零频处最平坦， 而在趋向阻带时衰减单调增大，巴特沃斯响应能够最大化滤波器的通带平坦度。该响应非常平坦，非常接近 DC信号，然后慢慢衰减至截止频率点为-3dB ，最终逼近 -20ndB/decade的衰减率，其中 n为滤波器的阶数。切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。本次课程设计将使用切比雪夫 1型带通滤波器对信号进行滤波。滤波器的结构框图如下图3所示：.连续混合输出采样带通滤波器（ 2000HZ ）信号图 3 滤波器的结构框图四 MATLAB 软件

14、介绍MATLAB 软件是由美国Mathworks 公司推出的用于数值计算和图形处理的科学计算系统环境。在MATLAB 环境下，用户可以集成地进行程序设计、数值计算、图形绘制、输入输出、文件管理等各项操作。MATLAB 提供了一个人机交互的数学系统环境， 该系统的基本数据结构是矩阵，在生成矩陈对象时， 不要求作明确的维数说明。与利用 C 语言或 FORTRAN语言作数值计算的程序设计相比，利用 MATLAB 可以节省大量的编程时间。 MATLAB 因其提供计算、 图形、文字处理的统一环境而深受欢迎。MTALAB 系统由五个主要部分组成，下面分别加以介绍。(1)MATALB 语言是高层次的矩阵数组

15、语言。具有条件控制、 函数调用、数据结构、输入输出、面向对象等程序语言特性。利用它既可以进行小规模编程，完成算法设计和算法实验的基本任务，也可以进行大规模编程， 开发复杂的应用程序。(2)MATLAB 工作环境这是对MATLAB 提供给用户使用的管理功能的总称。包括管理工作空间中的变量据输入输出的方式和方法，以及开发、调试、管理 M文件的各种工具。(3)MATLAB 图形系统的基础，包括完成2D 和 3D 数据图示、图像处理、动画生成、图形显示等功能的高层MATLAB 命令，也包括用户对图形图像等对象进行特性控制的低层MATLAB 命令，以及开发 GUI 应用程序的各种工具。.(4)MATLA

16、B 数学函数库这是对MATLAB 使用的各种数学算法的总称。包括各种初等函数的算法，也包括矩阵运算、矩阵分析等高层次数学算法。(5)MATLAB 应用程序接口 (API)这是 MATLAB 为用户提供的一个函数库，使得用户能够在 MATLAB 环境中使用 C 程序或 Fortran 程序，包括从 MATLAB 中调用于程序 (动态链接 )，读写 MAT 文件的功能。可以看出MATLAB 是一个功能十分强大的系统，是集数值计算、图形管理、程序开发为一体的环境。五 设计过程1.信号的产生首先产生一个连续输入信号，包含低频（f=30Hz ），中频（ f=250Hz ） ,高频（ f=500Hz ）分

17、量，并实现三个信号的合成。（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);%绘制 x1(t)的图形xlabel('t');ylabel('x1(t)');title(' 低频信号 ');grid;x2=sin(2*pi*t*f2);subplot(2,1,2);plot(x2);%绘制 x2(t)的图形xlabel('t');ylabel('x2(t)');tit

18、le(' 中频信号 ');.grid;x3=sin(2*pi*t*f3);figure(2);subplot(2,1,1);plot(x3);%绘制 x3(t)的图形xlabel('t');ylabel('x3(t)');title(' 高频信号 ');grid;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)');ti

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

20、的抽样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;（ 2）程序运行结果如图6、图 7 所示图 6 连续信号抽样结果波形图图 7 滤波前的频谱分析波形图（ 3）结果分析图6、图 7是对包含低频、中频、高频的连续信

21、号抽样，抽样点数取为512。从图中看出抽样的模型是对图6中的连续信号波形的采样，并且满足时域采样定.理。 MATLAB的优点是直接可以采用FFT进行快速傅里叶变换。本设计就是直接利用 FFT进行频谱分析，得到了 20HZ、250HZ、500HZ的频谱。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;%数字滤波器的归一化阻带截

22、止频率wp=2*fpl/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

23、);%求解第二个通带数字滤波器的频率响应H=h1+h2;W=w1+w2;%实现两个通带滤波器.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

24、;xlabel('频率 /HZ');ylabel('angle');title(' 切比雪夫 1 型带通滤波器的相频特性');（ 2）程序运行结果如图 8、图 9所示图 8 切比雪夫 1型双通带带通滤波器的幅频特性.图 9 切比雪夫 1型双通带带通滤波器的相频特性（ 3）分析设定 IIR 切比雪夫1 型带通数字滤波器指标：通带范围为:200-300Hz和500-600Hz, 第一个通带阻带上限为： 350HZ,阻带下限为 150Hz。第二个通带阻带上限为： 650HZ,阻带下限为 450Hz。通带最大衰减Ap=0.5dB ，阻带最小衰减为As=2

25、0dB，采样频率为 fsa=2000Hz 。其中， N, wpo=cheblord(wp,ws, Rp, As)该格式用于计算切比雪夫1 滤波器的阶数 N 和通带中心频率wpo 。wp 和 ws 分别为滤波器的通带边界频率和阻带截止频率的归一化值，要求 0<=wp 和 ws<=1,表示模拟频率 Fs/2。Rp 和 As 分别为通带最大衰减和阻带最小衰减。ws<wp 时，为高通滤波器。 B, A=cheby1(N, Rp ，wpo) ，该格式用于计算切比雪夫1 滤波器系统函数的分子分母多项式系数向量B 和 A。调用参数 N 和 wpo 分别为切比雪夫 1 滤波器的阶数N 和通带

26、中心频率。值得我们注意的是，切比雪夫滤波器分为 1 型，2 型滤波器， 1 型滤波器为通带内是等波纹的，阻带是单调递减的；2型滤波器为通带是单调递减的，而阻带是等波纹的。.4.对信号进行滤波（ 1）程序代码 y1=filter(B,A,x); y2=filter(C,D,x); y=y1+y2; figure(7) plot(t,y); xlabel('时间 (t)');ylabel('滤波后的信号 y');title(' 滤波后的信号 ');grid on;axis(0 0.05 -2 2);Y=fft(y,512);%对滤出来的信号进行频谱分

27、析w=(0:255)/256*1000;%频率横坐标的设定figure(8)plot(w,abs(Y(1:256);%绘制频谱图xlabel('频率 /Hz');ylabel('频率响应幅度 ');title(' 滤波后的频谱图 ');grid;（ 2）程序运行结果如图10、图 11 所示.图 10 滤波后信号时域波形图图 11 滤波后信号频谱图六 结果分析设计过程中，首先产生连续输入信号， 包含低频（f=30Hz ），中频（ f=250Hz ）,高频（ f=500Hz ）分量，然后对其进行采样，利用傅里叶变换进行频谱分析，并由带通滤波器的参数设

28、计带通滤波器对信号进行滤波处理，对应带通滤波器的通.带范围是（ 200，300）和（400,500），从运行结果图中可以看出，经过带通滤波器滤波后信号对应的频率为原信号中的中频分量（f=250Hz ）和高频分量（ f=500Hz) 。由上述结果显示， 在误差允许的范围内实验结果与理论结果相同。出现误差的原因：在设计滤波器的参数时并不是十分的准确，在不同计算机上运行MATLAB 时会有一定的偶然误差， 从而导致实验误差的存在。 并且采用切比雪夫滤波器的阻带衰减特性则有所改善。巴特沃斯带通滤波器的幅频特性随的增加而单调下降，当N 较小时，阻带幅频特性下降较慢，要想使其幅频特性接近理想带通滤波器，

29、就必须增加滤波器的阶数， 这就将导致模拟滤波器使用的原件增多，线路趋于复杂。 如果阶次一定， 巴特沃兹滤波器在通带内幅度特性是单调下降的，则在靠近截止频率处，幅度变化范围比较大，或者说，为了使通带内的衰减足够小，需要的阶次(N)很高，采用切比雪夫多项式逼近可以克服这一缺点。因为切比雪夫滤波器在通带范围内是等幅起伏的，所以同样的通带衰减， 它的阶数比巴特沃兹滤波器要小。七 心得体会及总结这次课程设计， 老师给我们的时间比较长， 在课设的前半段时间， 我们设计了三个正弦信号的叠加， FFT频谱分析，巴特沃斯带通滤波器及滤出正确的波形。在设计带通滤波器时我们遇到了麻烦，我们先设计了模拟低通滤波器，再

30、转换为带通滤波器，最后转换为数字带通滤波器， 我们所设计的滤波器满足了设计要求，但怎么也滤不出来波。后来经过查资料，我们发现MATLAB 有个很大的优点，可以不通过低通滤波器来设计带通滤波器，可以直接根据参数设置来求得带通滤波器。这使得我们的设计从头开始，直接设计巴特沃斯带通滤波器，成功把我们.所需要的波给滤出来了。 在实验课上， 经过王老师的指导， 如果我们想要做得更完美，可以再设计一种到两种带通滤波器。我们小组最后决定再设计一个FIR带通滤波器， IIR 切比雪夫带通滤波器。由于任务比较多，我们把任务进行分块，而我主要设计切比雪夫带通滤波器。在设计时，我先利用W= T 把模拟频率转换为数字

31、角频率，这涉及到归一化问题，经过我们查阅资料，在IIR 和 FIR 滤波器设计中，频率是用fs/2 进行归一化，如果是弧度的话，就是用pi 进行归一化。至于为什么要进行归一化，因为滤波器有各种类型，包括低通，高通，带通，带阻等，而且就一种类型而言，中心频率，截止频率又不同，但是设计表格不能有那么多，一般只给出归一化原型滤波器的设计数据。在设计出切比雪夫带通滤波器后，我们又设计出双通带的切比雪夫带通滤波器，最后出现了极好的滤波图。所以这次课程设计，真的使我们更好的理解了这学期学的数字信号处理的知识，也弥补了课堂上的不足，让我们学会了简单的MATLAB 的使用。参考文献1 桂志国，杨民，陈友兴，郝

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

33、;subplot(2,1,1);plot(x1);%绘制 x(t)的图形xlabel('t');ylabel('x1(t)');title(' 低频信号 ');grid;x2=sin(2*pi*t*f2);subplot(2,1,2);plot(x2);%绘制 x2(t) 的图形xlabel('t');ylabel('x2(t)');title(' 中频信号 ');grid;x3=sin(2*pi*t*f3);figure(2);subplot(2,1,1);plot(x3);%绘制 x3(t) 的

34、图形xlabel('t');ylabel('x3(t)');title(' 高频信号 ');grid;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;%对连续输入信号进行采样，进行频谱分析n=1:100;t=n/2000X=fft(x,512);w=(0:255)/256*1000;x=sin

35、(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3);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;%参数设置wh

36、p=0.3*pi;whs=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 /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,'.');

37、%绘制 Hamming 窗单位脉冲响应axis(0,N2-1,-0.4,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('H

38、amming窗幅频响应 ');subplot(3,1,3);plot(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(' 连续信号 ')

39、;grid;Y=fft(y,512);w=(0:255)/256*1000;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(

40、7);subplot(2,1,1);plot(w*2000/(2*pi),abs(H);%绘制带通频谱图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(' 连续信号 ');.grid