求教FFT或功率谱问题

zhujyabc

小弟最近采集到一组数据，但是画不出FFT或功率谱图，我的采集频率为25000HZ，请哪位高手指教一下 数据在附件里 我写的程序如下： load tt.txt; y=fft(tt,1024); p=y.*conj(y)/1024; ff=10000*(0:511)/1024; plot(ff,p(1:512))

songzy41

采样频率为25000HZ，为什么
ff=10000*(0:511)/1024;
是否应改为
ff=25000*(0:511)/1024;

zhujyabc

我改了 发现波形是这样的，感觉不对呀

zhujyabc

我改了 发现波形是这样的，感觉不对呀

songzy41

由于数据中的直流分量太大所造成的，在FFT分析前先除去直流分量：
load tt.txt;
tt=tt-mean(tt)
y=fft(tt,1024);
p=y.*conj(y)/1024;
ff=25000*(0:511)/1024;
plot(ff,p(1:512))
得图为

[ 本帖最后由 songzy41 于 2008-3-18 10:25 编辑 ]

zhujyabc

谢谢了 我知道怎么做了，但是还有个问题 就是直流分量怎么形成的哇？

wanyeqing2003

直流分量是信号的均值，也就是动力系统的平衡位置。

形成的原因可能有两个：1、信号本身的均值不为零；2、测试系统的零位点漂移。