数字信号处理实验六时域采样与信号的重建.docVIP

实验目的： 1.了解用MATLAB语言进行时域抽样与信号重建的方法 2.进一步加深对时域信号抽样与恢复的基本原理的理解 3.掌握采样频率的确定方法和内插公式的编程方法。 实验内容 1认真阅读并输入实验原理与方法中介绍的例子，观察输出波形曲线，理解每一条语句的含义。. 2．已知一个连续时间信号f(t)=sinc(t)。取最高有限带宽频率fm=1Hz。 （1）分别显示原连续时间信号波形和Fm=fm、Fm=2fm、Fm=3fm三种情况下抽样信号的波形。 实验程序： dt=0.1; f0=1; T0=1/f0; fm=f0; Tm=1/fm; t=-2:dt:2; f=sinc(t); subplot(4,1,1),plot(t,f,k); axis([min(t) max(t) 1.1*min(f) 1.1*max(f)]); title(原连续信号和抽样信号); for i=1:3; fs=i*fm; Ts=1/fs; n=-2:Ts:2; f=sinc(n); subplot(4,1,i+1),stem(n,f,filled,k); axis([min(n) max(n) 1.1*min(f) 1.1*max(f)]); end 实验截图： 求解原连续信号波形和抽样信号所对应的幅度谱。 实验程序： dt=0.1;t=-4:dt:4; N=length(t);f=sinc(t);Tm=1;fm=1/Tm; wm=2*pi*fm;k=1:N; w1=k*wm/N; F1=f*exp(-j*t*w1)*dt; subplot(4,1,1),plot(w1/(2*pi),abs(F1));grid axis([0 max(4*fm) 1.1*min(F1) 1.1*max(F1)]); for i=1:3; if i= 2 c=0 ,else c=0.2,end fs=(4-i+c)*fm; Ts=1/fs; n=-4:Ts:4; f=sinc(n); N=length(n); wm=2*pi*fs; k=1:N; w=k*wm/N; F=f*exp(-j*n*w)*Ts; subplot(4,1,5-i),plot(w/(2*pi),abs(F),k);grid axis([0 max(4*fm) 1.1*min(F) 1.1*max(F)]); end 实验截图： 用时域卷积方法（内插公式）重建信号。 实验程序： f0=1; T0=1/f0; dt=0.01; fm=f0; Tm=1/fm; t=-4:dt:4*f0; x=sinc(t); subplot(4,1,1),plot(t,x); axis([min(t),max(t),1.1*min(x),1.1*max(x)]); title(用时域卷积重建抽样信号) for i=1:3; fs=i*fm; Ts=1/fs; n=-4:(4*f0)/Ts; t1=-4:Ts:4*f0; x1=sinc(n/fs); T_N=ones(length(n),1)*t1-n*Ts*ones(1,length(t1)) xa=x1*sinc(fs*pi*T_N); subplot(4,1,i+1),plot(t1,xa); axis([min(t1),max(t1),1.1*min(xa),1.1*max(xa)]); end 实验截图：