ex3p17.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 12 行

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('t in msec');ylabel('xa(t)')
title('Analog Signal')
subplot(2,1,2);plot(W/(2*pi*1000),Xa*1000);
xlabel('Frequency in KHz');ylabel('Xa(jw)*1000')
title('Continuous -time Fourier Transform')