请教!傅立叶幅值谱的一些问题

yexiao0771

各位大虾们，请教一个概念的问题：
对于fft帮助文档里的程序代码
t = 0:0.001:0.6;  x = sin(2*pi*50*t)+sin(2*pi*120*t);
y = x + 2*randn(size(t));
plot(1000*t(1:50),y(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')

Y = fft(y,512);
The power spectrum, a measurement of the power at various frequencies, is
Pyy = Y.* conj(Y) / 512;                 
f = 1000*(0:256)/512;
plot(f,Pyy(1:257))
title('Frequency content of y')
xlabel('frequency (Hz)')


Pyy = Y.* conj(Y) / 512;            问：  本句求得是是功率谱 ，如果要求的是傅立叶幅值谱，如何写？
是Pyy =abs（Y）还是Pyy =abs（Y）/512？还是Pyy =abs（Y）^2/512？
需要除以采样点数吗？ 请高手指点！论坛上还是没有给出什么时候该乘2阿

[ 本帖最后由 yexiao0771 于 2008-1-19 17:38 编辑 ]

VibrationMaster

Pyy =abs（Y）/512

gonghai88

我也上来问一句，楼主的程序是基于一个给定生成的数据做分析的，如果我手头的数据存放在一个名为S的变量里，S为行或者列向量。程序该如何修改下。
我是新手，请大家指导下。

[ 本帖最后由 zhlong 于 2008-1-21 16:20 编辑 ]

VibrationMaster

Pyy =abs（fft（S)）/512