声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1099|回复: 3

[编程技巧] [求助]关于fft帮助中的一点疑问

[复制链接]
发表于 2006-4-21 17:05 | 显示全部楼层 |阅读模式

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

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

x
在Matlab help里面查到fft的相关信息(如下所示),有些地方不太明白,关于1000Hz的那一句哪位能帮我翻译一下吗?另外在程序中还有几个问题不太明白,能解释一下吗(具体位置见下面)?谢谢!

A common use of Fourier transforms is to find the frequency components of a signal buried in a noisy time domain signal. Consider data sampled at 1000 Hz. Form a signal containing 50 Hz and 120 Hz and corrupt it with some zero-mean random noise:
  1. t = 0:0.001:0.6;
  2. x = sin(2*pi*50*t)+sin(2*pi*120*t);
  3. y = x + 2*randn(size(t));
  4. plot(1000*t(1:50),y(1:50))
  5. title('Signal Corrupted with Zero-Mean Random Noise')
  6. xlabel('time (milliseconds)')

  7. It is difficult to identify the frequency components by looking at the original signal. Converting to the frequency domain, the discrete Fourier transform of the noisy signal y is found by taking the 512-point fast Fourier transform (FFT):
  8. Y = fft(y,512);
  9. The power spectrum, a measurement of the power at various frequencies, is

  10. Pyy = Y.* conj(Y) / 512;
  11. 这一句是什么意思?

  12. Graph the first 257 points (the other 255 points are redundant) on a meaningful frequency axis:

  13. f = 1000*(0:256)/512;
复制代码

频率f的取值为什么只取一半呢?
  1. plot(f,Pyy(1:257))
  2. title('Frequency content of y')
  3. xlabel('frequency (Hz)')
复制代码
回复
分享到:

使用道具 举报

发表于 2006-4-22 21:56 | 显示全部楼层
fft变换之后除以采样频率的二分之一才是对应的幅值
取一般是因为fft后,前后两部分是对称的
 楼主| 发表于 2006-4-23 09:56 | 显示全部楼层
谢谢教授的解释。我编了一个小程序,得到的图形是对称的,不知道是不是程序的某些地方是错的?具体问题网址是:
http://forum.vibunion.com/thread-11170-1-1.html
发表于 2006-5-12 11:08 | 显示全部楼层
Pyy = Y.* conj(Y) / 512;
这一句是什么意思?我也不大明白
happy教授指点一下,谢谢!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-9-25 19:24 , Processed in 0.051611 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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