验证采样定理.docVIP

目录 一、设计要求 - 2 - 二、原理及方法 - 2 - 1、时域采样定理的要点 - 2 - 2、信号的恢复与采样内插公式 - 3 - 三、实验步骤 - 3 - 四、MATLAB实现编程 - 4 - 五、实验结果与分析 - 12 - 1、时域波形与幅频特性曲线 - 12 - 2、结果分析与总结 - 13 - 六、总结与体会 - 15 - 七、参考文献 - 17 - 一、设计要求 利用Matlab对连续信号进行采样，对比采样频率，，三种情况下的频谱，并与理想频谱进行分析，进而从实验上验证时域采样定理。 二、原理及方法 1、时域采样定理的要点 （1）对模拟信号以间隔T进行时域等间隔理想采样，形成的采样信号的频谱是原模拟信号频谱以采样角频率（）为周期进行周期延拓。公式为： （2）采样频率必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。 理想采样信号和模拟信号之间的关系为： 对上式进行傅立叶变换，得到： 在上式的积分号内只有当时，才有非零值，因此： 2、信号的恢复与采样内插公式 内插公式表明模拟信号等于各采样点函数值乘以对应内插函数的总和，即只要采样频率高于信号频率的两倍，模拟信号就可用它的采样信号代表，而不会丢失任何信息。这种理想低通滤波器的模拟信号完全等于模拟信号，是一种无失真的恢复。 三、实验步骤 1、画出连续时间信号的时域波形及其幅频特性曲线，信号为其中A=100，a=10，=50*2*。 2、 对信号进行采样，得到采样序列 ，画出采样频率分别为200Hz，100 Hz，60 Hz时的采样序列波形。 3、对不同采样频率下的采样序列进行频谱分析，绘制其幅频曲线，对比各频率下采样序列和的幅频曲线。 4、由采样序列恢复出连续时间信号 ，画出其时域波形，与原连续时间信号的时域波形进行对比。 四、MATLAB实现编程 1、时域采样程序 clc; clear; close all; fs=500; fs1=200; fs2=100; fs3=60; %设置四重采样频率 t=0:1/fs:0.1; %采集信号长度为0.1s A=00; b=50*2*pi; a=10; %连续信号xa(t)的参数 xa=exp(-a*t).*sin(b*t); k=0:511; f=fs*k/512; %由wk=2pik/512=2pifT求得模拟频率f w=2*pi*k/512; Xa=xa*exp(-j*[1:length(xa)]*w);%近似模拟信号频谱 T1=1/fs1; t1=0:T1:0.1; %采集信号长度为0.1s x1=A*exp(-a.*t1).*sin(b*t1); %200Hz抽样序列x1(n) X1=x1*exp(-j*[1:length(x1)]*w);%x1(n)的512点DTFT T2=1/fs2; t2=0:T2:0.1; %采集信号长度为0.1s x2=A*exp(-a.*t2).*sin(b*t2); %100Hz抽样序列x2(n) X2=x2*exp(-j*[1:length(x2)]*w);%x2(n)的512点DTFT T3=1/fs3; t3=0:T3:0.1; %采集信号长度为0.1s x3=A*exp(-a.*t3).*sin(b*t3); %60Hz抽样序列x3(n) X3=x3*exp(-j*[1:length(x3)]*w);%x3(n)的512点DTFT figure(1); subplot(2,2,1); plot(t,xa); axis([0,max(t),min(xa),max(xa)]); title(模拟信号); xlabel(t(s)); %X-轴标注 ylabel(Xa(t)); %Y-轴标注 line([0,max(t)],[0,0]); %添加网格和绘X-轴线 subplot(2,2,2); plot(f,abs(Xa)/max(abs(Xa))); title(模拟信号的幅度频谱); axis([0,500,0,1]); xlabel(f(Hz)); ylabel(|Xa(jf)|); subplot(2,2,3); ste