声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1239|回复: 7

[FFT] 我想知道我的程序为什么出现四个峰,不是只有60和90的频率么?

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

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

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

x
我想知道我的程序为什么出现四个峰,不是只有60和90的频率么?
clear
fs=200 ;
t=0:1/fs:0.6;
x=1.23*sin(2*pi*90*t+5)+1.46*sin(2*pi*60*t);
n=1024;
      y=fft(x,n);

n=0:length(y)-1;
f=fs*n/length(y);
plot(f,abs(y))

麻烦高手运行一下

[ 本帖最后由 江河之北 于 2008-6-25 17:14 编辑 ]
回复
分享到:

使用道具 举报

发表于 2008-6-25 18:02 | 显示全部楼层
只需取频率0-100之间,而100-200之间对应的是负频率。
 楼主| 发表于 2008-6-25 18:09 | 显示全部楼层

回复 2楼 的帖子

谢谢您老师 :handshake :handshake
 楼主| 发表于 2008-6-25 18:23 | 显示全部楼层

回复 2楼 的帖子

宋老师,如果我再计算这个频谱图的功率谱该如何计算呢,需要调用psd这个命令么?或者就是幅值的平方呢?
发表于 2008-6-25 19:55 | 显示全部楼层

回复 楼主 的帖子

对普通信号做傅里叶变换得到的谱图是对称的,要得到单边谱的话可以先对信号做希尔伯特变换使信号成为解析信号,这时再做傅里叶变换就是单边谱了。功率谱估计有多种方法,对应的matlab命令也有多个,你可以用psd,pwelch,periodogram等,具体格式用help看一下就行了
 楼主| 发表于 2008-6-26 09:31 | 显示全部楼层

回复 5楼 的帖子

我是想知道在频谱的基础上怎么求功率谱,psd好像是对顶楼的x来求功率谱的吧
 楼主| 发表于 2008-6-26 10:39 | 显示全部楼层

回复 5楼 的帖子

对信号做希尔伯特变换使信号成为解析信号,能否给出具体的程序,小弟不是学信号处理的呀
发表于 2008-6-26 12:49 | 显示全部楼层
用hilbert(x)对x做希尔伯特变换,并作为解析信号的虚部,x作为解析信号的实部,然后再做傅里叶变换就行了。呵呵,有些事情还是得自己动手啊。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-13 22:08 , Processed in 0.068113 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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