请高手将记录的一段波形进行傅立叶分析并画图

laoxuzi

楼上历害,怎么看出来有直流分量的.
只是有些程序看不明白,
N2=N/2;
n2=1:N2;
df=fs/N;
F=(n2-1)*df;
plot(F,abs(S(n2)));
这个与
Y=fft(x,1024);
p=Y.*conj(Y)/1024;
f=fs*(0:512)/1024;
plot(f,p(1:513))
有什么不同?

[ 本帖最后由 laoxuzi 于 2009-9-26 11:16 编辑 ]

ChaChing

高人songzy41做的频谱与我的相同了!
楼主something is wrong!! LZ给的资料与图有问题!?

songzy41

原帖由 laoxuzi 于 2009-9-26 11:07 发表
楼上历害,怎么看出来有直流分量的.
只是有些程序看不明白,
N2=N/2;
n2=1:N2;
df=fs/N;
F=(n2-1)*df;
plot(F,abs(S(n2)));
这个与
Y=fft(x,1024);
p=Y.*conj(Y)/1024;
f=fs*(0:512)/1024;
plot(f,p(1:513) ...

从信号的波形图上可看到信号的起伏不是在0附近，而是在-0.015附近，直观就告诉我们信号的均值在-0.015附近，FFT后必然有很大的直流分量；同时LZ在121.jpg图上也看到在0频率上有一个很大的值，这就表明直流分量很大。
N2=N/2;
n2=1:N2;
df=fs/N;
F=(n2-1)*df;
plot(F,abs(S(n2)));
这一段程序主要目的是要计算FFT变换后的频率刻度F，这样在plot中才能把频率标注出来。当然也可以改为
f=fs*(0:12056)/24112;
plot(f,abs(S(1:12057)));
这时要先知道数据长为24112。

laoxuzi

谢谢songzy41,ChaChing帮助,很是受益,以后还要跟二位多学习.:@P