[信号实验2.docVIP

实验2 信号的抽样与重构 姓 名 ： 王 林 学 号 ： 班 级 ： 电气5班 实验目的： 学习在Matlab软件下信号的抽样与重构过程的计算机仿真； 学习掌握抽样定理，通过不同抽样频率的取样与重构，加深对定理的理解； 学习Matlab软件在信号处理中的应用。 实验内容： 设研究的信号如下： (1)分别采用1KHz 、2KHz和5KHz进行采样，作出原信号和采样信号的图形； (2)利用重构公式对采样信号进行重构并画出图形，然后与原信号进行比较并给出结论； (3)分析信号间的频谱关系并解释重构的差异。 程序如下： %构造模拟信号 Dt=0.001;t=-0.01:Dt:0.01; xa=exp(-1000*abs(t)); %定义两种采样周期并对原信号进行采样 Ts1=0.001;n1=-5:1:5; Ts2=0.0005;n2=-25:1:25; Ts3=0.0002;n3=-30:1:30; x1=exp(-500*abs(n1*Ts1)); x2=exp(-500*abs(n2*Ts2)); x3=exp(-500*abs(n3*Ts3)); %将原信号和采样信号以图形方式表示出来 figure(Name,模拟信号的取样。,NumberTitle,off); subplot(4,1,1);plot(t*500,xa); subplot(4,1,2);stem(n1*Ts1*500,x1); subplot(4,1,3);stem(n2*Ts2*500,x2); subplot(4,1,4);stem(n3*Ts3*500,x3); 运行结果如下： 结果分析:由运行结果得到的图形可得，在所取采样频率1KHZ，2KHZ和5KHZ中，利用采样频率5KHZ得到的采样信号的图形与原信号图形越接近。 (2)程序如下： Fs1=1/Ts1;Fs2=1/Ts2;Fs3=1/Ts3;n1Ts1=n1*Ts1;n2Ts2=n2*Ts2;n3Ts3=n3*Ts3; xa1=x1*sinc(Fs1*(ones(length(n1Ts1),1)*t-n1Ts1*ones(1,length(t)))); error1=max(abs(xa-xa1)); xa2=x2*sinc(Fs2*(ones(length(n2Ts2),1)*t-n2Ts2*ones(1,length(t)))); error2=max(abs(xa-xa2)); xa3=x3*sinc(Fs3*(ones(length(n3Ts3),1)*t-n3Ts3*ones(1,length(t)))); error3=max(abs(xa-xa3)); %显示重构后的信号波形图 figure(Name,信号的重构,NumberTitle,off); subplot(4,1,1);plot(t*500,xa); subplot(4,1,2);plot(t*500,xa1); subplot(4,1,3);plot(t*500,xa2); subplot(4,1,4);plot(t*500,xa3); 运行结果如下： 结果分析：实验所得信号图形如上，通过原信号图形与重构得到的信号图形比较可得，重构后得到的信号图形与与原信号图形存在着一定的差异，由图可得采用采样频率2KHZ得到的重构信号图形与原信号图形比较接近。 程序如下： Wmax=2*pi*770; 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)]; w=pi*k/K; X1=x1*exp(-j*n1*w); X2=x2*exp(-j*n2*w); X3=x3*exp(-j*n3*w); X1=real(X1); X2=real(X2); X3=real(X3); w=[-fliplr(w),w(2:K+1)]; X1=[fliplr(X1),X1(2:K+1)]; X2=[fliplr(X2),X2(2:K+1)]; X3=[fliplr(X3),X3(2:K+1)]; %作出信号的频谱图 figure(Name,信号的频谱。,NumberTitle,off); subplot(4,1,1);plot(W/(2*pi*500),Xa*500,blue) subplot(4,1,2);plot(w/pi,X1,green