example3_15.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

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');