有关于FFT变换的相频特性问题
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),'.');%对比二者相频特性图像
结果出现的图像并不相符,甚至相差比较明显,想请教如何解决这个问题。 一来你的f是一个衰减函数,是一个连续频率成分的信号,而FFT只是对平稳信号,
二来只用FFT来求相位一般都是不对的. 原帖由 yangzj 于 2007-4-14 14:42 发表
一来你的f是一个衰减函数,是一个连续频率成分的信号,而FFT只是对平稳信号,
二来只用FFT来求相位一般都是不对的.
感谢楼上的朋友。
请问有什么办法解决么?我实际要处理的信号比这个复杂的多,笔头去推算傅立叶变换太不实际了。
页:
[1]