clc;clear; wm=1; wc=wm; Ts=pi/wm; ws=2*pi/Ts; n=-100:100; nTs=n*Ts f=2*sinc(nTs/pi); Dt=0.0025;t=-15:Dt:15; fa=f*Ts*wc/pi*sinc((wc/pi)*(ones( length(nTs),1)*t-nTs'*ones(1, length(t)))); t1=-15:0.25:15; f1=2*sinc(t1/pi); subplot(121); stem(t1,f1); xlabel('kTs'); ylabel('f(kTs)'); title('sa(t)=sinc(t/pi)的临界抽样信号'); subplot(122); plot(t, fa) xlabel('t'); ylabel('fa(t)'); title('由sa(t)=sinc(t/pi)的临界抽样信号重构sa(t)'); grid;

这段代码是关于信号的抽样和重构的，其中采用了临界抽样，信号为基带信号，经过抽样后进行重构。抽样定理告诉我们，如果采样频率大于等于信号的最高频率的两倍，就可以完美地还原原始信号。在这段代码中，采样频率为 $\omega_c=\omega_m=\pi/T_s$，即采样间隔为 $T_s$，信号为 $s_a(t)=\operatorname{sinc}(t/\pi)$，经过采样后重构出来的信号为 $f_a(t)$，可以看到重构后的信号与原始信号非常接近。