声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1066|回复: 3

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

[复制链接]
发表于 2008-1-19 17:15 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
各位大虾们,请教一个概念的问题:
对于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 编辑 ]
回复
分享到:

使用道具 举报

发表于 2008-1-19 20:35 | 显示全部楼层
Pyy =abs(Y)/512
发表于 2008-1-21 10:38 | 显示全部楼层
我也上来问一句,楼主的程序是基于一个给定生成的数据做分析的,如果我手头的数据存放在一个名为S的变量里,S为行或者列向量。程序该如何修改下。
我是新手,请大家指导下。

[ 本帖最后由 zhlong 于 2008-1-21 16:20 编辑 ]
发表于 2008-1-21 14:56 | 显示全部楼层
Pyy =abs(fft(S))/512
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-14 18:46 , Processed in 0.060468 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表