有关于FFT变换的相频特性问题

xiaopeng66662

N=128;
t=linspace(0,3,N);
Ts=t(2)-t(1);
Ws=1/Ts;

f=exp(-3*t);%原函数
F=fft(f);%用FFT求傅立叶变换
Fc=fftshift(F)*Ts;

W=Ws*(-N/2:(N/2)-1)/N;
Fa=1./(3+j*W);%笔算出的傅立叶变换

plot(W,angle(Fa),W,angle(Fc),'.');%对比二者相频特性图像

结果出现的图像并不相符，甚至相差比较明显，想请教如何解决这个问题。

yangzj

一来你的f是一个衰减函数,是一个连续频率成分的信号,而FFT只是对平稳信号,
二来只用FFT来求相位一般都是不对的.

xiaopeng66662

原帖由 yangzj 于 2007-4-14 14:42 发表
一来你的f是一个衰减函数,是一个连续频率成分的信号,而FFT只是对平稳信号,
二来只用FFT来求相位一般都是不对的.

感谢楼上的朋友。
请问有什么办法解决么？我实际要处理的信号比这个复杂的多，笔头去推算傅立叶变换太不实际了。