信号与系统实验2015解读.doc

信号与系统实验报告 13级通信四班 赵豆 学号实验一 连续时间信号的采样 一、 实验目的 1、进一步加深对采样定理和连续信号傅立叶变换的理解。 2、了解MATLAB提供的在相邻点间内插的几种方法。 3、使用sinc函数并且研究时域中的混叠问题 二、 实验步骤 1．复习采样定理和采样信号的频谱 采样定理 如果采样频率大于有限带宽信号带宽的两倍，即 （1） 则该信号可以由它的采样值重构。否则就会在中产生混叠。该有限带宽模拟信号的被称为乃魁斯特频率。 必须注意，在被采样以后，表示的最高模拟频率为Hz（或）。 2．熟悉如何用MATLAB语言实现模拟信号表示 严格地说，除了用符号处理工具箱(Symbolics)外，不可能用MATLAB来分析模拟信号。然而如果用时间增量足够小的很密的网格对采样，就可得到一根平滑的曲线和足够长的最大时间来显示所有的模态。这样就可以进行近似分析。令是栅网的间隔且，则 （2） 可以用一个数组来仿真一个模拟信号。不要混淆采样周期和栅网间隔，因为后者是MATLAB中严格地用来表示模拟信号的。类似地，付利叶变换关系也可根据（2）近似为： （3） 现在，如果（也就是）是有限长度的。则公式（3）与离散付利叶变换关系相似，因而可以用同样的方式以MATLAB来实现，以便分析采样现象。 3．根据提供的例子程序，按照要求编写实验用程序； 三、实验内容 （1）通过例一熟悉用MATLAB语言实现描绘连续信号的频谱的过程，并在MATLAB语言环境中验证例1的结果； 例1 令，求出并绘制其付利叶变换。 解：根据傅立叶变换公式有 （4） 因为是一个实偶信号，所以它是一个实值函数。为了用数值方法估计，必须先把用一个栅格序列来近似。 利用，注意可以用一个在（或等效地[-5,5]毫秒）之间的有限长度信号来近似。 模拟过程及结果如下： Exp(-5)=0 信号的时间范围选取-0.005—0.005 w=2*pi*2000 ws=4*pi*2000 fs=2000t=0.00005Ts=0.00025 % 模拟信号 Dt=0.00005; t=-0.005:Dt:0.005; xa=exp(-1000*abs(t)); %连续时间傅立叶变换 Wmax=2*pi*2000; K=500; k=0:1:K; W=k*Wmax/K; Xa=xa*exp(-j*t*W)*Dt; Xa=real(Xa); W=[-fliplr(W),W(2:501)];%频率从-Wmax to Wmax Xa=[fliplr(Xa),Xa(2:501)];%Xa 介于 -Wmax和 Wmax之间 subplot(1,1,1) subplot(2,1,1);plot(t*1000,xa); xlabel(t 毫秒); ylabel(xa(t)); title(模拟信号) subplot(2,1,2); plot(W/(2*pi*1000),Xa*1000); xlabel(频率（单位：Hz)); ylabel(Xa(jW)*1000) title(连续时间傅立叶变换) （2）仿照例2用MATLAB语言实现对连续信号 的采样；并验证采样定理。 例2 为了研究采样对频域各量的影响，这里用两个不同的采样频率对例1中的进行采样。 a.以样本/秒采样得到。求并画出。 b.以样本/秒采样得到。求并画出。 解：a.因为的带宽是2kHz，奈魁斯特频率为4000样本/秒。它比所给的采样频率低，因此混叠将（几乎）不存在。 % 模拟信号 Dt=0.00005; t=-0.005:Dt:0.005; xa=exp(-1000*abs(t)); %离散时间信号 Ts=0.0002;n=-25:1:25;x=exp(-1000*abs(n*Ts)); %离散时间傅立叶变换 K=500; k=0:1:K; w=pi*k/K; X=x*exp(-j*n*w); X=real(X); w=[-fliplr(w),w(2:K+1)]; X=[fliplr(X),X(2:K+1)]; subplot(1,1,1) subplot(2,1,1);plot(t*1000,xa); xla