试验报告试验一

1、实验一离散时间信号的 MATLAB 实现课程名称：数字信号处理实验日期：2016.4.21班级：数媒 1403姓名：李欣学号：一、实验目的掌握：(1) 序列的生成; (2) 离散序列的基本运算与波形变换的 MATLAB 实现;(3) 线性时不变系统求解的 MATLAB 实现; (4) 序列的傅里叶变换与信号重构MATLAB 实现二、实验要求按照实验材料内容，利用 MATLAB 编程实现相关任务。三、实验步骤（以实验材料第四部分“序列的傅里叶变换与信号重构 MATLAB 实现”为例， 给出实验步骤描述、相关代码等。）程序如下：Dt=0.00005;t=-0.005: Dt:0.005;xa=ex

2、p(-1000*abs(t);%模拟信号Wmax=2*pi*2000; K=500;k=0:1:K;W=k*Wmax/K;%连续时间傅立叶变换Xa=xa*exp(-1i*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,'.'); xlabel('t in msec'); ylabel('xa(t)');gtext('模拟信号');subplot(2,1,2

3、); plot(W/(2*pi*1000),Xa*1000,'.'); xlabel('Frequence in KHz'); ylabel('Xa(jw)*1000');gtext('连续时间傅立叶变换'); 程序运行结果：程序如下：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(-1i

4、*n'*w);X=real(X);w=-fliplr(w),w(2:501);X=fliplr(X),X(2:501);figure subplot(2,1,1); plot(t*1000,xa,'.');ylabel('xa1(t)');gtext('离散信号'); hold on ; stem(n*Ts*1000,x); hold off ; subplot(2,1,2); plot(w/pi,X,'.');ylabel('X1(jw)');gtext('离散时间傅立叶变换'); gte

5、xt('Ts=0.2msec')程序运行结果： 程序如下： Ts1=0.0002; Fs1=1/Ts1;n1=-25:1:25;nTs1=n1*Ts1;%离散时间信号x1=exp(-1000*abs(nTs1);Ts2=0.001;Fs2=1/Ts2;n2=-5:1:5;nTs2=n2*Ts2;x2=exp(-1000*abs(nTs2); Dt=0.00005;t=-0.005:Dt:0.005;%模拟信号重构xa1=x1*sinc(Fs1*(ones(length(nTs1),1)*t-nTs1'*ones(1,length(t); xa2=x2*sinc(Fs2*

6、(ones(length(nTs2),1)*t-nTs2'*ones(1,length(t); subplot(2,1,1);plot(t*1000,xa1,'.');ylabel('xa1(t)');title('从 x1(n)重构模拟信号 xa(t)'); hold on;stem(n1*Ts1*1000,x1); hold off;subplot(2,1,2);plot(t*1000,xa2,'.'); ylabel('xa2(t)');title('从 x2(n)重构模拟信号 xa(t)'); hold on;stem(n2*Ts2*1000,x2); hold off程序运行结果：四、实验体会(1) 掌握了 MATLAB 的基本操作以及学会