example3_15b.m

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

Dt=0.00005;t=-0.005:Dt:0.005;xa=exp(-1000*abs(t));
N=length(xa);
Xa=Dt*dft(xa,N);
Xa=abs(real(Xa));
k=(1:N/4)*100/1000;
stem(k,Xa(1:N/4));
xlabel('Khz');