数字信号处理实验五

实验五：抽样定理 一、实验目的 1、了解用 MATLAB 语言进行时域、频域抽样及信号重建的方法。 2、进一步加深对时域、频域抽样定理的基本原理的理解。 3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和内插公式的编程方法。 二、实验内容及步骤 1、阅读并输入实验原理中介绍的例题程序，观察输出的数据和图形，结合基本原理理解每一条语句的含义。 2、已知一个连续时间信号 f(t)=sinc(t)，取最高有限带宽频率 fm=1Hz。 （1）分别显示原连续信号波形和 Fs=fm、Fs=2fm、Fs=3fm三种情况下抽样信号的波形； dt=0.1;f0=1;T0=1/f0; fm=1;Tm=1/fm; t=-2:dt:2; f=sinc(t); subplot(4,1,1);plot(t,f); 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'); axis([min(n),max(n),1.1*min(f),1.1*max(f)]); 课程名称 数字信号处理 实验成绩 指导教师 实 验 报 告 院系 信息工程学院 班级 学号 姓名 日期 end -2-1 .5-1-0 .500 .511 .5200 .51原连续信 号 和 抽 样 信 号-2-1 .5-1-0 .500 .511 .5200 .51-2-1 .5-1-0 .500 .511 .5200 .51-2-1 .5-1-0 .500 .511 .5200 .51 （2）求解原连续信号和抽样信号的幅度谱； dt=0.1;f0=1;T0=1/f0; fm=1;Tm=1/fm; t=-2:dt:2; N=length(t); f=sinc(t); wm=2*pi*fm; k=0:N-1; w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt; subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3; if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2; N=length(n); f=sinc(n); wm=2*pi*fs; k=0:N-1; w=k*wm/N; F=f*exp(-j*n'*w)*Ts; subplot(4,1,i+1);plot(w/(2*pi),abs(F)); axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end 00 .511 .522 .533 .540 .20 .40 .60 .811 .200 .511 .522 .533 .54012-2-1 .5-1-0 .500 .511 .5200 .51-2-1 .5-1-0 .500 .511 .5200 .51 （3）用时域卷积的方法（内插公式）重建信号。 dt=0.01;f0=1;T0=1/f0; fm=1;Tm=1/fm; t=0:dt:3*T0; x =sinc(t); su bplot(4,1,1);plot(t,x ); ax is([min(t),max (t),1.1*min(x ),1.1*max (x )]); title('用时域卷积重建抽样信号'); for...