k1=0:999; k2=-999:-1; m1=length(k1); m2=length(k2); f=[fs0*k2/m2,fs0*k1/m1]; w=[-2*pi*k2/m2,2*pi*k1/m1]; fx1=eval(fy); FX1=fx1*exp(-j*[1:length(fx1)]'*w); figure subplot(2,1,1),plot(t,fx1,'r') title('原信号'), xlabel('时间t (s)') axis([min(t),max(t),min(fx1),max(fx1)]) subplot(2,1,2),plot(f,abs(FX1),'r') title('原信号幅度频谱') , xlabel('频率f (Hz)') axis([-100,100,0,max(abs(FX1))+5]) Ts=1/fs; t1=-tp:Ts:tp; f1=[fs*k2/m2,fs*k1/m1]; t=t1; fz=eval(fy); FZ=fz*exp(-j*[1:length(fz)]'*w); figure subplot(2,1,1),stem(t,fz,'.'), title('取样信号') , xlabel('时间t (s)') line([min(t),max(t)],[0,0]) subplot(2,1,2),plot(f1,abs(FZ),'m') title('取样信号幅度频谱') , xlabel('频率f (Hz)') 2.恢复程序： function fh=huifu(fz,fs) T=1/fs; dt=T/10; tp=; t=-tp:dt:tp; n=-tp/T:tp/T; TMN=ones(length(n),1)*t-n'*T*ones(1,length(t)); fh=fz*sinc(fs*TMN); k1=0:999; k2=-999:-1; m1=length(k1); m2=length(k2); w=[-2*pi*k2/m2,2*pi*k1/m1]; FH=fh*exp(-j*[1:length(fh)]'*w); figure subplot(2,1,1),plot(t,fh,'g'), st1=sprintf('由取样频率fs=%d',fs); st2='恢复后的信号'; st=[st1,st2]; title(st) , xlabel('时间t (s)') axis([min(t),max(t),min(fh),max(fh)]) line([min(t),max(t)],[0,0]) f=[10*fs*k2/m2,10*fs*k1/m1]; subplot(2,1,2),plot(f,abs(FH),'g') title('恢复后信号的频谱') , xlabel('频率f (Hz)') axis([-100,100,0,max(abs(FH))+2]); 3.三种采样实现程序及其显示结果： ①欠采样：f1='cos(2*pi*80*t)+2*sin(2*pi*30*t)+cos(2*pi*40*t-pi/3)'; fs2=caiyang(f1,140); fr2=huifu(fs2,140); 图 fs=140Hz原信号波形及频谱 图 fs=140Hz取样信号波形及频谱 图 fs=140Hz恢复后信号波形及频谱 ② 临界采样：f1='cos(2*pi*80*t)+2*sin(2*pi*30*t)+cos(2*pi*40*t-pi/3)'; fs2=caiyang(f1,160); fr2=huifu(fs2,160); 图 fs=160Hz原信号波形及频谱 图 fs=160H取样信号波形及频谱 图 fs=160Hz恢复后波形及频谱 ③过采样：f1='cos(2*pi*80*t)+2*sin(2*pi*30*t)+cos(2*pi*40*t-pi/3)'; fs2=caiyang(f1,200); fr2=huifu(fs2,200);
实验六 抽样定理的MATLAB仿真
