《数字信号处理》实验指导书.doc

《《数字信号处理》实验指导书.doc》由会员分享，可在线阅读，更多相关《《数字信号处理》实验指导书.doc（57页珍藏版）》请在三一文库上搜索。

1、数字信号处理实验指导书实验一 MATLAB数学软件的使用一、实验目的1、掌握 Matlab 软件使用的基本方法；2、熟悉 Matlab 的离散时间序列的产生。二、实验仪器与软件 1. PC机 1台 2. MATLAB7.0环境 三、实验原理1、窗口简介MATLAB环境是一种为数值计算、数据分析和图形显示服务的交互式的环境。MATLAB有3种窗口，即：命令窗口（The Command Window）、m-文件编辑窗口（The Edit Window）和图形窗口（The Figure Window），而Simulink另外又有Simulink模型编辑窗口。1）命令窗口（The Command Wi

2、ndow）当MATLAB启动后，出现的最大的窗口就是命令窗口。用户可以在提示符“”后面输入交互的命令，这些命令就立即被执行。在MATLAB中，一连串命令可以放置在一个文件中，不必把它们直接在命令窗口内输入。在命令窗口中输入该文件名，这一连串命令就被执行了。因为这样的文件都是以“.m”为后缀，所以称为m-文件。2）m文件编辑窗口（The Edit Window）我们可以用m文件编辑窗口来产生新的m文件，或者编辑已经存在的m-文件。在MATLAB主界面上选择菜单“File/New/M-file”就打开了一个新的m-文件编辑窗口；选择菜单“File/Open”就可以打开一个已经存在的m-文件，并且可

3、以在这个窗口中编辑这个m-文件。2、序列的产生可通过命令窗口键入不带.m的程序名并按回车来执行。也可从命令窗口中的File菜单上，选择Open，选定所要的M文件，再使用Tools菜单下的Run命令，可运行该M文件。四、实验内容：编制程序产生下述4种信号，并绘出其图形。1、产生单位抽样序列 在MATLAB中可以利用zeros()函数实现。如果在时间轴上延迟了k个单位，得到即： 2单位阶跃序列 在MATLAB中可以利用ones()函数实现。 3正弦序列在MATLAB中4复指数序列在MATLAB中五、实验要求 利用所学知识，完成上述各项实验内容，并将实验过程和实验步骤和结果写在报告中。53实验二 数

4、字信号处理中常见算法的MATLAB实现一、实验目的1、了解序列移位与周期延拓运算；2、了解序列翻褶与序列累加运算；3、了解序列卷积和相关运算。二、实验仪器与软件1. PC机 1台2. MATLAB7.0环境 三、实验原理1、序列移位与周期延拓运算；（1）序列移位的数学描述为：。序列移位的MATLAB实现：y=x;ny=nx-m。（2）序列周期延拓的数学描述为：。其中M表示延拓周期。序列周期延拓的MATLAB实现：ny=nxs:nxf;y=x(mod(ny,M)+1)。试编制序列移位与周期延拓运算的MATLAB程序。2、序列翻褶与序列累加运算；（1）序列翻褶的数学描述为：。序列翻褶的MATLAB

5、可由函数fliplr来实现。该函数的功能将行向量左右翻转，其调用格式：y=fliplr（x）。（2）序列累加的数学描述为：。序列累加的MATLAB可由函数cumsum来实现。其调用格式：y=cumsum（x）。试编制序列翻褶与序列累加的MATLAB程序。3、序列卷积和相关运算。 （1）两序列卷积运算的数学描述为：。两序列卷积运算的MATLAB实现：y=conv(x1,x2)。序列x1(n)和x2(n)必须长度有限。 请通过MATLAB计算下列卷积，并图示各序列及其卷积结果。（2）两序列相关运算的数学描述为：。两序列相关运算的MATLAB实现：y=xcorr(x1,x2)。说明，若两序列具有相同

6、的长度M，则相关序列y(n)的长度为2M-1。若两序列的长度不同，则短者将自动添0。请通过MATLAB计算下列相关，并图示各序列及其相关结果。四、实验步骤与内容1、认真复习本实验中的序列运算，阅读本实验原理与方法。2、编制实验原理中要求的程序。3、调通并运行实验程序，完成实验内容五、实验要求1、列出本次实验编写的所有文件及各项实验结果的曲线，加注必要的说明。2、总结实验体会及实验中存在的问题。 实验三 信号、系统及系统响应一、 实验目的1、熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解；2、熟悉时域离散系统的时域特性；3、利用卷积方法观察分析系统的时域特性；4、掌握序列傅立

7、叶变换的计算机实现方法，利用序列的傅立叶变换对连续信号、离散信号及系统响应进行频域分析。二、 实验原理及方法采样是连续信号数字处理的第一个关键环节。对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生变化以及信号信息不丢失的条件，而且可以加深对傅立叶变换、Z变换和序列傅立叶变换之间关系式的理解。对一个连续信号进行理想采样的过程可用下式表示：其中为的理想采样，p(t)为周期脉冲，即 的傅立叶变换为 上式表明为的周期延拓。其延拓周期为采样角频率（）。只有满足采样定理时，才不会发生频率混叠失真。在实验时可以用序列的傅立叶变换来计算。公式如下：离散信号和系统在时域均可用序列来表示。为了在实验中观

8、察分析各种序列的频域特性，通常对在0，2上进行M点采样来观察分析。对长度为N的有限长序列x(n)，有： 其中，k=0,1,M-1三、 实验内容及步骤1、认真复习采样理论，离散信号与系统，线性卷积，序列的傅立叶变换及性质等有关内容，阅读本实验原理与方法。2、编制实验用主程序及相应子程序。3、调通并运行实验程序，完成下列实验内容：分析采样序列的特性，产生采样信号序列，使A=444.128，a=50,。（的无失真采样频率约为1000Hz）。a. 取采样频率=1kHz，即T=1ms。观察所得采样的幅频特性|和原图中的幅频特性曲线在折叠频率附近有无明显差别。应当注意，实验中所得频谱是用序列的傅立叶变换公

9、式求得的，所以在频率度量上存在关系：，为数字频率，为模拟频率。b. 改变采样频率，=300Hz，观察 |的变化，并做记录（打印曲线）；进一步降低采样频率，=200Hz，观察频谱混叠是否明显存在，说明原因，并记录（打印）这时的 |曲线。 a步实验结果如下图所示：由图形可知，当采样频率为1000Hz时，采样序列在折叠频率附近处，即=处无明显频谱混叠。b步实验结果如下图所示：由图可知，当采样频率进一步降低时，主瓣宽度逐渐变宽，频率混叠现象也逐渐严重。存在较明显的失真现象。时域离散信号、系统和系统响应分析。a. 观察信号和系统的时域和频域特性；利用线性卷积求信号通过系统的响应y(n)，比较所求响应y(

10、n)和的时域及频域特性，注意它们之间有无差别。绘图说明，并用所学理论解释所得结果。实验结果如下图所示：b.观察系统对信号的响应特性。利用线性卷积求系统响应y(n)，并判断y(n)图形及其非零值序列长度是否与理论结果一致，对，说出一种定性判断y(n)图形正确与否的方法。调用序列傅立叶变换数值计算子程序，求得，观察|特性曲线，定性判断结果的正确性。改变的长度，取N=5，重复该实验。注意参数变化的影响，说明变化前后的差异，并解释所得结果。实验结果如下图所示：N=10：N=5：欲判断结果正确与否，可以先对其进行运算，算出其卷积，再与图形对照。当N=10时，峰值较高，且峰值很窄，变换之后图形频带主值部分

11、比较集中；N=5时情况与之相反。*卷积定理的验证。将实验中的信号换为，使a=0.4，A=1，T=1，重复实验a，打印|曲线；对主程序做简单修改，按式（10.3.9）计算，并绘出|曲线，与前面直接对y(n)进行傅立叶变换所得幅频特性曲线进行比较，验证时域卷积定理。实验所得结果如下：（*的内容对黄海学院的同学不做要求）四、实验要求 利用所学知识，完成上述各项实验内容，并将实验过程和实验步骤和结果写在报告中。实验四 用FFT作谱分析一、 实验目的1、 进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质)。2、熟悉FFT算法原理和

12、FFT子程序的应用。3、学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。二、 实验步骤1、 复习DFT的定义、性质和用DFT作谱分析的有关内容。2、 复习FFT算法原理与编程思想，并对照DITFFT运算流图和程序框图，读懂本实验提供的FFT子程序。3、 编制信号产生子程序，产生以下典型信号供谱分析： ， ， ， 其它n ， ， 0 ， 其它n 应当注意，如果给出的是连续信号，则首先要根据其最高频率确定采样速率以及由频率分辨率选择采样点数N，然后对其进行软件采样(即计算，)，产生对应序列。对信号，频率分辨率的选择要以能分辨开其中的

13、三个频率对应的谱线为准则。对周期序列，最好截取周期的整数倍进行谱分析，否则有可能产生较大的分析误差。4、 编写主程序下图给出了主程序框图，供参考。本实验提供FFT子程序和通用绘图子程序。开始读入长度N调用信号产生子程序产生实验信号调用绘图子程序(函数)绘制时间序列波形图调用FFT子程序(函数)计算信号的DFT 调用绘图子程序(函数)绘制曲线结束主程序框图三、 实验结果直接运行程序，按照实验内容及程序提示键入18，分别对及、进行谱分析。输出的波形及其8点DFT和16点DFT，的16点、32点和64点采样序列及其DFT。1、及其8点和16点DFT2、及其8点和16点DFT3、及其8点和16点DFT

14、4、的8点和16点波形及其DFT5、的8点和16点波形及其DFT6、的16点、32点和64点采样序列波形及其DFT选7时，计算并图示和及其DFT。程序自动计算并绘图验证DFT的共轭对称性。当N=16时，。即为的共轭对称分量，而是的共轭反对称分量。根据DFT的共轭对称性，应有以下结果：的8点和16点波形及其DFT绘出和的模。它们正是图中16点的和。选8时，计算并图示和及其DFT。程序自动计算并绘图验证DFT的共轭对称性的第二种形式：如果，则，。其中，。的8点和16点DFT程序计算结果如下：及，正好与图中的16点及相同。及，正好与图中16点的及相同。四、 实验总结本实验主要是求、的DFT变换。其中

15、是直接给出了离散序列，而、则是经过、运算得到的、。离散傅立叶变换可以看作是在时的Z变换，即表明的N点DFT是的Z变换在单位圆上的N点等间隔采样。离散傅立叶变换也可以看作在时的傅立叶变换，即表明可以看作的傅立叶变换在区间上的N点等间隔采样。五、 思考题1、在N=8时，和的幅频特性会相同吗？为什么？N=16呢？2、如果周期信号的周期预先不知道，如何用FFT进行谱分析？ 实验五 用双线性变换法设计IIR数字滤波器一、 实验目的 1、熟悉用双线性变换法设计IIR数字滤波器的原理与方法； 2、掌握数字滤波器的计算机仿真方法； 3、通过观察对实际心电图信号的滤波作用，获得数字滤波的感性知识。二、 实验内容

16、及原理 1、用双线性变换法设计一个巴特沃斯低通IIR数字滤波器。设计指标参数为：在通带内截止频率低于0.2时，最大衰减小于1dB；在阻带内0.3,频率区间上，最小衰减大于15dB。 2、以0.02为采样间隔，打印出数字滤波器在频率区间0，/2上的幅频响应特性曲线。 3、用所设计的滤波器对实际心电图信号采样序列进行仿真滤波处理，并分别打印出滤波前后的心电图信号波形图，观察总结滤波作用与效果。教材例中已求出满足本实验要求的数字滤波系统函数：,式中 A=0.09036, 实验所得结果如下所示：对频率的压缩符合下列公式： 这样的变换叫做双线性变换。用双线性变换法来设计数字滤波器，由于从s面映射到s1面

17、具有非线性频率压缩的特点，因此不可能产生频率混叠现象，而且转换成的H(z)是因果稳定的，这是双线性变换法的最大优点。其缺点是w与之间的非线性关系直接影响数字滤波器频香逼真的模仿模拟滤波器的频响。数字滤波器的输入和输出均为数字信号，通过一定的运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分。数字滤波器可以通过模拟其网络传输函数进行实现。如图中所示，滤波器对其高于截止频率的频段产生很高的衰减，所得信号较之原信号剔除了高频的成分。三、 思考题用双线性变换法设计数字滤波器过程中，变换公式 中T的取值，对设计结果有无影响？为什么？实验六 用窗函数法设计FIR数字滤波器一、 实验目的1、掌握

18、用窗函数法设计FIR数字滤波器的原理和方法。2、熟悉线性相位FIR数字滤波器的特性。3、了解各种窗函数对滤波特性的影响。二、 实验原理与方法如果所希望的滤波器的理想频率响应函数为（），则其对应的单位脉冲响应为 （41）窗函数设计法的基本原理时用有限长单位脉冲响应序列h(n)逼近。用窗函数将截断，并进行加权处理，得到： （42）就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数为 （43）式中，N为所选窗函数的长度。由书本第七章可知，用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。

19、 这样选定窗函数类型和长度N后，求出单位脉冲响应，并按式（43）求出。是否满足要求，要进行验算。一般在h(n)尾部加零使长度满足2的整数次幂，以便用FFT计算。如果要观察细节，补零点数增多即可。如果补满足要求，则要重新选择窗函数类型和长度N，再次验算，直至满足要求。如果要求线性相位特性，则h(n)还必须满足： 根据上式中的正负号和长度N的奇偶性又将线性相位FIR滤波器分为四类。要根据所设计的滤波器特性正确选择其中一类。例如，要设计线性相位低通特性，可选择一类，而不能选择一类。三、 实验内容1、用升余弦窗设计一线性相位低通FIR数字滤波器，截止频率。窗口长度N15，33。要求在两种窗口长度情况下

20、，分别求出h(n),打印出相应的幅频特性和相频特性曲线，观察3dB带宽和20dB带宽。总结窗口长度N对滤波特性的影响。 设计低通FIR数字滤波器时，一般以理想低通滤波特性为逼近函数，即 其中 2、 n33，用四种窗函数设计线性相位低通滤波器。绘制相应的幅频特性曲线，观察3dB和20dB带宽以及阻带最小衰减，比较四种窗函数对滤波器特性的影响。四、 实验结果运行程序，根据实验内容要求和程序提示选择要进行的实验参数：1、窗函数长度（the length of the window : N=）默认值为 N15，选择范围为1064。2、逼近理想低通滤波器截止频率（：）默认值为，选择范围为0.5, 。3、

21、窗函数类型（请选择窗函数类型：）输入1：（boxcar）矩形窗。输入2：hamming 窗。输入3：hanning窗。输入4：blackman窗。4）、输出图形： 矩形窗（N15） 矩形窗（N33） Hamming窗（N15） Hamming窗（N33） Hanning窗（N15） Hanning窗（N33） Blackman窗（N15） Blackman窗（N33）五、 实验总结及心得用窗函数法设计FIR滤波器的主要特点：设FT为希望逼近的频响特性函数,H() =FTh(n)为用窗函数法设计的实际滤波器的频响函数。通常取H()相应的理想频响特性作为。知识要点如下：1、希望逼近的理想滤波器频响

22、函数的表达式。因为数字滤波器一般要求设计成线性相位特性，所以必须满足上述线性相位FIR滤波器的频域特点。2、熟悉各种常用窗函数的技术指标和加窗后对滤波特性的影响，这样才能根据设计指标正确选择窗函数类型及其长度N。3、检验设计结果：窗函数法的设计结果单位脉冲响应。而检验一般在频域进行。所以要计算检验3dB截止频率和阻带最小衰减，其计算量相当大，必须用计算机进行。4、熟悉窗函数设计法的特点：设计过程简单.方便实用。但边界频率不易精确控制。所以设计完以后，必须检验结果。六、 思考题1、如果给定通带截止频率和阻带截止频率以及阻带最小衰减，如何用窗函数法设计线性相位低通滤波器？写出设计步骤。2、如果要求

23、用窗函数法设计带通滤波器，且给定上、下边带截止频率为和，试求理想带通的单位脉冲响应。附录 MATLAB下的数字信号处理实现示例实验一MATLAB数学软件的使用1、绘制单位样本序列clf;n=-10:20;u=zeros(1,10) 1 zeros(1,20);stem(n,u);xlabel(时间序号n);ylabel(振幅);title(单位样本序列);axis(-10 20 0 1.2);2、绘制长度为20的单位样本序列clf;n=1:20;u=ones(1,20);stem(n,u);xlabel(时间序号n);ylabel(振幅);title(单位样本序列);axis(-10 20 0

24、 1.2);3、绘制复指数序列clf;c=-(1/12)+(pi/6)*i;k=2;n=0:40;x=k*exp(c*n);subplot(2,1,1);stem(n,real(x);xlabel(时间序号n);ylabel(振幅);title(实部);subplot(2,1,2);stem(n,imag(x);xlabel(时间序号n);ylabel(振幅);title(虚部);4、产生一个正弦序列clf;n=0:40;f=0.1;phase=0;A=1.5;arg=2*pi*f*n-phase;x=A*cos(arg);stem(n,x);axis(0 40 -2 2);grid;titl

25、e(正弦序列);xlabel(时间序号n);ylabel(振幅);实验二数字信号处理中常见算法的MATLAB实现1、序列移位与周期延拓运算N=24;M=8;m=3;n=0:N-1;x1=(0.8).n; %生成指数序列x2=(n=0)&(nM); %形成指数序列x=x1.*x2; %截取操作形成新序列x(n)xm=zeros(1,N);for k=m+1:m+M xm(k)=x(k-m); %产生序列移位x(n-3)end;xc=x(mod(n,M)+1); %产生x(n)的周期延拓x(n)8xcm=x(mod(n-m,M)+1); %产生移位序列x(n-3)的周期延拓x(n-3)8subpl

26、ot(4,1,1),stem(n,x,.);ylabel(x(n);subplot(4,1,2),stem(n,xm,.);ylabel(x(n-3);subplot(4,1,3),stem(n,xc,.);ylabel(x(n)_8);subplot(4,1,4),stem(n,xcm,.);ylabel(x(n-3)_8);2、序列翻褶与叠加n=0:N-1; %x(n)序列的时间序列x=3*exp(-0.2*n); %x(n)序列大小y=fliplr(x); %x(n)序列翻褶n1=-fliplr(n); %时间序列的翻褶（翻褶点为原点）n2=fliplr(-(n-3); %在指定位置m=

27、3处的时间序列的翻褶subplot(2,1,1);stem(n,x); %绘制x(n)序列xlabel(n),ylabel(x(n);subplot(2,1,2);stem(n1,y); %绘制y(n)=x(-n)序列xlabel(n),ylabel(y(n)=x(-n);s=cumsum(x); %求累加序列figure(2);subplot(2,1,1);stem(n2,y); %绘制点为3处的翻褶序列y(n)=x(-n+3)xlabel(n),ylabel(y(n)=x(-n+3);subplot(2,1,2);stem(n,s); %绘制s(n)序列xlabel(n),ylabel(s

28、(n);3、 序列的卷积与相关运算（1） 卷积程序： Nx=20;Nh=10;m=5; %设定Nx,Nh和位移值mn=0:Nx-1;x1=(0.9).n; %产生x1(n)x2=zeros(1,Nx+m);for k=m+1:m+Nx %产生x2(n)=x1(n-m) x2(k)=x1(k-m); %完成位移endnh=0:Nh-1;h1=ones(1,Nh); %产生h1(n)h2=h1;y1=conv(x1,h1); %计算y1(n)=x1(n)*h1(n)subplot(3,2,1);stem(n,x1); %绘制x1(n)xlabel(n),ylabel(x1(n);subplot(3

29、,2,3);stem(nh,h1); %绘制h1(n) xlabel(n),ylabel(h1(n);subplot(3,2,5);stem(y1); %绘制y1(n)xlabel(n),ylabel(y1(n);y2=conv(x2,h2); %计算y2(n)=x2(n)*h2(n)subplot(3,2,2);stem(x2); %绘制x2(n) xlabel(n),ylabel(x2(n);subplot(3,2,4);stem(nh,h2); %绘制h2(n) xlabel(n),ylabel(h2(n);subplot(3,2,6);stem(y2); %绘制y2(n) xlabel

30、(n),ylabel(y2(n);(2)相关运算Nx=20;Nh=10;m=5; %设定Nx,Nh和位移值mn=0:Nx-1;x1=(0.9).n; %产生x1(n)x2=zeros(1,Nx+m);for k=m+1:m+Nx %产生x2(n)=x1(n-m) x2(k)=x1(k-m); %完成位移endnh=0:Nh-1;h1=ones(1,Nh); %产生h1(n)h2=h1;y1=xcorr(x1,h1); %计算相关subplot(3,2,1);stem(n,x1); %绘制x1(n)xlabel(n),ylabel(x1(n);subplot(3,2,3);stem(nh,h1)

31、; %绘制h1(n) xlabel(n),ylabel(h1(n);subplot(3,2,5);stem(y1); %绘制y1(n)xlabel(n),ylabel(y1(n);y2=xcorr(x2,h2); %计算相关subplot(3,2,2);stem(x2); %绘制x2(n) xlabel(n),ylabel(x2(n);subplot(3,2,4);stem(nh,h2); %绘制h2(n) xlabel(n),ylabel(h2(n);subplot(3,2,6);stem(y2); %绘制y2(n) xlabel(n),ylabel(y2(n); 卷积运算结果图 相关运算结

32、果图 实验三 信号、系统和系统响应1、理想采样信号序列 分析采样序列的特性:a,bf=input(采样率fs=);subplot(1,2,1);n=0:49;A=444.128;a=50*1.414*pi;w0=50*1.414*pi;xa=A*exp(-a)*n/f).*sin(w0*n/f);stem(n,xa,.);xlabel(n);ylabel(x(n);title(fs取1000Hz);subplot(1,2,2);k=-200:200;c=xa*(exp(-j*pi/100).(n*k);y=abs(c);plot(k/100,y);xlabel(w/);ylabel(Xa(ej

33、w);title(|X(ejw)|);时域离散信号、系统和响应分析a:观察信号xb(n)和系统ha(n)的时域和频域特性ha(n)的时域和频域特性：子程序1（本文档后所附）xb(n)的时域和频域特性：子程序2（本文档后所附）信号xb(n)和系统ha(n)的响应分析：程序如下n=1:4;ha=zeros(1,4);ha(1)=1; ha(2)=2.5; ha(3)=2.5; ha(4)=1;close all;%subplot(5,1,1);stem(ha);title(系统)n=1:1;x=zeros(1);x(1)=1;%subplot(5,1,2);stem(x);title(输入信号);

34、y=conv(x,ha);%subplot(5,1,3);stem(y);title(输出信号);subplot(3,1,1);stem(y);title(输出信号);n=1:4;k=-25:25;X=y*(exp(-j*pi/12.5).(n*k);magX=abs(X);%subplot(5,1,4);stem(magX);title(y(n)的幅度谱);subplot(3,1,2); stem(magX);title(y(n)的幅度谱);angX=angle(X);%subplot(5,1,5);stem(angX);title(y(n)的相位谱);subplot(3,1,3);stem

35、(angX);title(y(n)的相位谱);b：观察系统hb(n)对信号xc(n)的响应特性。N=10时：n=1:50;x=sign(sign(10-n)+1);close all;n=1:50;hb=sign(sign(10-n)+1);y=conv(x,hb);subplot(3,1,1);stem(y);title(输出信号);n=1:99;k=-25:25;X=y*(exp(-j*pi/12.5).(n*k);magX=abs(X);subplot(3,1,2); stem(magX);title(y(n)的幅度谱);angX=angle(X);subplot(3,1,3);stem

36、(angX);title(y(n)的相位谱);N=5时：n=1:50;x=sign(sign(5-n)+1);close all;n=1:50;hb=sign(sign(10-n)+1);y=conv(x,hb);subplot(3,1,1);stem(y);title(输出信号);n=1:99;k=-25:25;X=y*(exp(-j*pi/12.5).(n*k);magX=abs(X);subplot(3,1,2); stem(magX);title(y(n)的幅度谱);angX=angle(X);subplot(3,1,3);stem(angX);title(y(n)的相位谱);*卷积定

37、理验证n=1:50;hb=zeros(1,50);hb(1)=1; hb(2)=2.5; hb(3)=2.5; hb(4)=1;close all;subplot(2,1,1);stem(hb);title(系统h(n);k=-25:25;Hb=hb*(exp(-j*pi/12.5).(n*k);magHb=abs(Hb);subplot(2,1,2); plot(magHb);title(hb(n)的幅度谱);figure;n=1:50;A=1;a=0.4;T=1;w0=2.0734;x=A*exp(-a*n*T).*sin(w0*n*T);subplot(2,1,1);stem(x);ti

38、tle(输入信号x(n);k=-25:25;X=x*(exp(-j*pi/12.5).(n*k);magX=abs(X);subplot(2,1,2); plot(magX);title(x(n)的幅度谱);figure;y=conv(x,hb);subplot(2,1,1);stem(y);title(输出信号y(n);n=1:99; k=-25:25;Y=y*(exp(-j*pi/12.5).(n*k);magY=abs(Y);subplot(2,1,2); plot(magY);title(y(n)的幅度谱);figure;%验证结果XHb=X.*Hb;Subplot(2,1,1);pl

39、ot(abs(XHb);title(x(n)的幅度谱与hb(n)的幅度谱相乘);Subplot(2,1,2);plot(abs(Y);title(y(n)的幅度谱);参考子程序1、单位脉冲序列的幅度谱及相位谱n=1:50;x=zeros(1,50);x(1)=1;close all;subplot(3,1,1);stem(x);title(dan wei chong ji xin hao xu lie);k=-25:25;X=x*(exp(-j*pi/12.5).(n*k);magX=abs(X);subplot(3,1,2);stem(magX);title(fu du pu);angX=a

40、ngle(X);subplot(3,1,3);stem(angX);title(xiang wei pu);2、特定冲激串n=1:50;x=zeros(1,50);x(1)=1;x(2)=2.5;x(3)=2.5;x(4)=1;close all;subplot(3,1,1);stem(x);title(hb(n)k=-25:25;X=x*(exp(-j*pi/12.5).(n*k);magX=abs(X);subplot(3,1,2);stem(magX);title(hb(n)的幅度谱);angX=angle(X);subplot(3,1,3);stem(angX);title(hb(n)的相位谱);3、矩形序列的