temp.m

来自「西交大刘树棠写的《数字信号处理》那本书的所有源代码」· M 代码 · 共 13 行

Dt=0.00005;t=-0.005:Dt:0.005;xa=exp(-1000*abs(t));
Ts=0.0002;n=-25:1:25;x=exp(-1000*abs(n*Ts));
k=0:50; N=51
w=pi*k/(N-1);
X=dft(x,N);X=real(X);
subplot(2,1,1);plot(t*1000,xa);
xlabel('ms');ylabel('x1(n)');
title('discrete time signal');hold on;
stem(n*Ts*1000,x);hold off
subplot(2,1,2);stem(w/pi,X);
xlabel('w');ylabel('X');
title('DTFT');