《实验一信号、系统及系统响应.docVIP

实验一 信号、系统及系统响应 实验目的： 1、熟悉常用离散时间信号的产生方法； 2、利用卷积方法观察分析系统的时域特性； 3、熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解； 4、熟悉MATLAB编程。 二、实验内容： 产生如下序列，并绘出离散图： （1），≤≤5； （2），值可由键盘输入，即编写一个产生矩形序列的函数； （3），≤≤； （4），≤≤实部和虚部分开显示。n=0:1:15 n=0:1:15 x4=exp((pi/4)*n*i) x4r=real(x4) subplot(1,2,1) stem(n,x4r,.) x4i=imag(x4) subplot(1,2,2) stem(n,x4i,.) 设线性时不变系统的单位取样响应为，设输入序列为：，求系统的输出。要求显示、和3个序列的图形。 n:0:1:9 h=0.9.^n subplot(3,2,1) title(h(n)) stem(h) x=ones(1,10) subplot(3,2,2) title(x(n)) stem(x) subplot(3,2,3) title(y(n)) conv(x,h)=y stem(y) 已知连续时间信号为，求： （1）的傅立叶变换； （2）取样频率为5000Hz，绘出； （3）取样频率为2000Hz，绘出； （4）重构两种取样频率下恢复出的时域连续信号。 提示：严格说来，在MATLAB中一般不能精确分析模拟信号，但以充分小的时间间隔取样时，可产生平滑的图形，当采样足够长的时间时，就可以实现近似地分析。本连续时间信号的傅立叶变换为（请查《信号与系统》）。在精度下，为的带限信号，因此若取，采样的信号就能够很好地逼近原信号。 4、分析采用序列的特性，验证时域采样定理。 （1）对连续信号： 进行采样，可得到采样序列 ， 0≤n＜50 其中A=444.128，，。采样频率有3种，分别是，，，显示不同采样频率下的幅频特性，观察频谱混叠现象。 三、实验原理： 由公式y(n)=x(n)*h(n)即输入序列与系统单位取样响应的时域卷积为系统输出，利用Matlab中的线性卷积函数conv即可实现。 根据傅立叶变换公式 ，可得到：，据此可以通过将积分转化为累加和来求，即得到时域模拟信号的傅立叶变换。 对时域信号进行采样，当取样频率大于其最高频率两倍时，可以无失真的恢复出原始信号，当取样频率小于其最高频率两倍时，会产生失真。 信号重建可由公式 得到，可以利用Ｍatlab中的函数sinc(x)=sin(pi*x)/（pi*x）来实现。 四、实验报告要求： 按实验步骤附上程序和实验得到的图形，并对结果进行分析和解释； 请将卷积conv()得到的结果与理论计算结果作一个大致比较； 总结实验中的主要结论； 附录： 绘离散时间信号图时，用stem函数。 函数fliplr的功能是翻转。 函数axis的功能是调整显示图形的坐标范围。 3、实验内容（3）中模拟信号的傅立叶变换参考程序： 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)]; Xa=[fliplr(Xa),Xa(2:501)]; figure(1) subplot(2,1,1) plot(t*1000,xa,r); xlabel(t in msec.); ylabel(x(t)); subplot(2,1,2) plot(W/(2*pi*1000),Xa*1000,r) xlabel(Frequency in KHZ); ylabel(X(W));