mvut实验五rv.docVIP

实验五 连续信号的采样与重构 实验目的 验证采样定理； 熟悉信号的抽样与恢复过程； 通过实验观察欠采样时信号频谱的混跌现象； 掌握采样前后信号频谱的变化，加深对采样定理的理解； 掌握采样频率的确定方法。 实验原理 一个有限频宽的连续信号f（t），其最高频率为Wm，经过等间隔抽样后，只要抽样频率Ws不小于信号最高频率Wm的2倍，既满足Ws=2Wm,就能从抽样信号fs（t）中恢复原信号，得到f0（t）。俩信号相比没有失真，只有幅度和相位的差异。一般把低的抽样频率Wsmin=2Wm成为乃奎斯特抽样频率。 实验内容 正弦信号的采样。 源程序： clf; t=0:0.0005:1; f=13; xa=cos(2*pi*f*t); subplot(2,1,1); plot(t,xa);grid xlabel(ê±??￡?mses);ylabel(·ù?μ); title(á?D?ê±??D?o?x_{a}(t)); axis([0 1 -1.2 1.2]); subplot(2,1,2); T=0.1; n=0:T:1; xs=cos(2*pi*f*n); k=0:length(n)-1; stem(k,xs);grid; xlabel(ê±??￡?msec);ylabel(·ù?μ); title(á?D?ê±??D?o?x[n]); axis([0 (length(n)-1) -1.2 1.2]); 实验结果及图。 采样的性质。 源程序： clf; t=0:0.005:10; xa=2*t.*exp(-t); subplot(2,2,1); plot(t,xa);grid; xlabel(ê±??,msec);ylabel(·ù?μ); title(á?D?ê±??D?o?x_{a}(t)); subplot(2,2,2); wa=0:10/511:10; ha=freqs(2,[1 2 1],wa); plot(wa/(2*pi),abs(ha));grid; xlabel(?μ?ê￡?kHz);ylabel(·ù?μ); title(|X_{a}(j\-mega)|); axis([0 5/pi 0 2]); subplot(2,2,3); T=1; n=0:T:10; xs=2*n.*exp(-n); k=0:length(n)-1; stem(k,xs);grid; xlabel(ê±??,n);ylabel(·ù?μ); title(à?é￠ê±??D?o?x[n]); subplot(2,2,4); wd=0:pi/255:pi; hd=freqz(xs,1,wd); plot(wd/(T*pi),T*abs(hd));grid; xlabel(?μ?ê￡?kHz);ylabel(·ù?μ); title(|X(e^{j\-mega})|); axis([0 1/T 0 2]); 实验结果及图。 模拟低通滤波器设计。 源程序： clf; Fp=3500;Fs=4500; Wp=2*pi*Fp;Ws=2*pi*Fs; [N,Wn]=buttord(Wp,Ws,0.5,30,s); [b a]=butter(N,Wn,s); wa=0:(3*Ws)/511:3*Ws; h=freqs(b,a,wa); plot(wa/(2*pi),20*log(abs(h)));grid; xlabel(Frequency,Hz);ylabel(Gain,dB); title(Gain response); axis([0 3*Fs -60 5]); 实验结果及图。 时域欠采样。 源程序： clf; n=0:49; m=0:50*3-1; x=sin(2*pi*0.042*m); y=x([1:3:length(x)]); subplot(2,1,1); stem(n,x(1:50));axis([0 50 -1.2 1.2]); title(ê?è?DòáD); subplot(2,1,2); stem(n,y);axis([0 50 -1.2 1.2]); title(ê?3?DòáD); 实验结果及图。 程序设计实验 源程序： clf; t=0:0.0005:1; f=150; xa=3*sin(2*pi*f*t); subplot(2,2,1); plot(t,xa);grid; xlabel(ê±??￡?mses);ylabel(·ù?μ); title(á?D?ê±??D?o?x1_{a}(t)); axis([0 2.5/150 -3.2 3.2]); subplot(2,2,3); T=0.0008; n=0:T:1; xs=3*sin(2*pi*5120*n); stem(n,xs);grid;