MATLAB抽样定理验证

目的 通过MATLAB编程实现对时域抽样定理的验证，加深抽样定理的理解。同时训练应用计算机分析问题的能力。 任务 连续信号ftcos8*pi*t2*sin40*pi*tcos24*pi*t，经过理想抽样后得到抽样信号fst，通过理想低通滤波器后重构信号ft。 方法 1、确定ft的最高频率fm。对于无限带宽信号，确定最高频率fm的方法设其频谱的模降到10-5左右时的频率为fm。 2、确定Nyquist抽样间隔TN。选定两个抽样时间TSTN，TSTN。 3、MATLAB的理想抽样为 n-200200;nTsn*Ts; 或 nTs-0.04Ts0.04 4、抽样信号通过理想低通滤波器的响应 理想低通滤波器的冲激响应为 系统响应为 由于 所以 MATLAB计算为 ftfs*Ts*wc/pi*sincwc/pi*oneslengthnTs,1*t-nTs*ones1,lengtht; 要求（画出6幅图） 当TSTN时 1、在一幅图中画原连续信号ft和抽样信号fSt。ft是包络线，fSt是离散信号。 2、画出重构的信号yt。 3、画出误差图，即 errorabsft-yt的波形。 当TSTN时同样可画出3幅图。 a wm40*pi; wc1.2*wm; 理想低通截止频率 Ts[0.02 0.03]; NlengthTs; for k1N; n-100100; nTsn*Tsk; fscos8*pi*nTs2*sin40*pi*nTscos24*pi*nTs.*unTspi-unTs-pi; t-0.250.0010.25; ftfs*Tsk*wc/pi*sincwc/pi*oneslengthnTs,1*t-nTs*ones1,lengtht; t1-0.250.0010.25; f1cos8*pi*t12*sin40*pi*t1cos24*pi*t1.*ut10.25-ut1-0.25; 在一副图中画原连续信号ft和样信号f_st。 figure3*k-2 plott1,f1,r,linewidth,2,hold on stemnTs,fs,grid on axis[-0.25 0.25 -4 4] line[-0.25 0.25],[0 0],color,k line[0 0],[-4 4],color,k xlabelnTs,ylabelfnTs; title[抽样信号Ts,num2strTsk,时的抽样信号fnTs] legend包络线,抽样信号,0 hold off 画重构的信号 figure3*k-1 plott,ft,grid on axis[-0.25 0.25 -4 4] line[-0.25 0.25],[0 0],color,k ine[0 0],[-4 4],color,k xlabelt,ylabelft; title由fnTs信号重建得到的信号 画出错信号 errorabsft-f1; figure3*k plott,error,grid on axis[-0.25 0.25 -0.1 4] line[-0.25 0.25],[0 0],color,k line[0 0],[-4 4],color,k xlabelt,ylabelerrort; title重建信号与原余弦信号的绝对误差 end 本文来自CSDN博客，转载请标明处出http//