example3_15.m

来自「数字信号处理的比较全面的程序代码」· M 代码 · 共 11 行

Dt=0.00005;t=-0.005:Dt:0.005;xa=exp(-1000*abs(t));
Wmax=2*pi*2000;K=500;k=0:1:K;W=k*Wmax/K;
Xa=xa*exp(-j*t'*W)*Dt;
Xa=real(Xa);
W=[-fliplr(W),W(2:501)];
Xa=[fliplr(Xa),Xa(2:501)];
subplot(2,1,1);plot(t*1000,xa);
xlabel('ms');ylabel('xa(t)');
subplot(2,1,2);plot(W/(2*pi*1000),Xa*1000);
xlabel('kHz');ylabel('Xa(jw)*1000');