fft.m

用MATLAB实现的DSP函数

N=128;
t=linspace(0,3,N);
f=2*exp(-3*t); %2e(-3t)
Ts=t(2)-t(1);
Ws=2*pi/Ts;
F=fft(f);
Fp=F(1:N/2+1)*Ts;
W=Ws*(0:N/2)/N;
Fa=2./(3+j*W);
plot(W,abs(Fa),'red',W,abs(Fp),'+');
xlabel('频率,弧度/秒'),ylabel('|F(w)|');