% fft_exp.m
% 連續(xù)信號傅立葉變換比較
N=128;
% 估計的時間點序列
t=linspace(0,3,N);
f=2*exp(-3*t);
% 采樣時間
Ts=t(2)-t(1);
% 采用頻率
Ws=2*pi/Ts;
% 計算傅立葉變換結果
F=fft(f);
% 將傅立葉變換的半值翻轉(zhuǎn)
Fc=fftshift(F)*Ts;
% 頻率軸
W=Ws*(-N/2:(N/2)-1)/N;
% 理論分析傅立葉變換
Fa=2./(3+j*W);
% 繪值