fft_exp.asv

matlab从入门到精通第6章书后源码对初学者很有帮助

% fft_exp.m
% 连续信号傅立叶变换比较
N=128;
% 估计的时间点序列
t=linspace(0,3,N);
f=2*exp(-3*t);
% 采样时间
Ts=t(2)-t(1);
% 采用频率
Ws=2*pi/Ts;
% 计算傅立叶变换结果
F=fft(f);
% 将傅立叶变换的半值翻转
Fc=fftshift(F)*Ts;
% 频率轴
W=Ws*(-N/2:(N/2)-1)/N;
% 理论分析傅立叶变换
Fa=2./(3+j*W);
% 绘值