声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3450|回复: 25

[HHT] 边际谱 问题求助

[复制链接]
发表于 2007-6-11 22:57 | 显示全部楼层 |阅读模式

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

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

x
clear;
fs=1200;    %fs为采样频率;
N=1200;     %采样点数
t=1/fs:1/fs:N/fs;
y1=1*sin(800*pi*t+900);
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y1);
[E,tt1]=toimage(A,fa,tt,length(tt));
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs;  
end
f=(1:N-2)/N*fs;
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');

其中信号明明频率应该是400的,为何图形得到是800呢?

[ 本帖最后由 zhlong 于 2007-6-13 17:46 编辑 ]

评分

2

查看全部评分

回复
分享到:

使用道具 举报

发表于 2007-6-11 23:08 | 显示全部楼层

回复 #1 dengbangfei 的帖子

哥们 请问这个信号频率是多少?你真的知道吗?应该是400
 楼主| 发表于 2007-6-11 23:12 | 显示全部楼层

回复 #2 zhangnan3509 的帖子

对啊,就是400,可是图形出来的是800,不是有问题吗?
发表于 2007-6-11 23:15 | 显示全部楼层

回复 #3 dengbangfei 的帖子

你看清楚了 这个正弦的频率是多少?
 楼主| 发表于 2007-6-11 23:23 | 显示全部楼层

回复 #4 zhangnan3509 的帖子

我先问你sin(2*pi*t)这个正弦信号频率应该是多少呢?难道不是1吗?
发表于 2007-6-11 23:24 | 显示全部楼层

回复 #5 dengbangfei 的帖子

800*pi=2*pi*f,这个频率是多少
发表于 2007-6-11 23:25 | 显示全部楼层
f=(1:N-2)/N*fs;

不是fs吧,应该是fs/2
发表于 2007-6-11 23:28 | 显示全部楼层

回复 #7 zhlong 的帖子

这没错!求边际谱的频率序列,如果除以2就对应不上了。这出自破皇的边际谱

[ 本帖最后由 zhangnan3509 于 2007-6-11 23:30 编辑 ]
 楼主| 发表于 2007-6-11 23:31 | 显示全部楼层

回复 #6 zhangnan3509 的帖子

我们的意思是一样的,难道是程序有问题?请指教
发表于 2007-6-11 23:32 | 显示全部楼层

回复 #9 dengbangfei 的帖子

我们的意思是不一样的!你连频率都算错了。
 楼主| 发表于 2007-6-11 23:33 | 显示全部楼层

回复 #6 zhangnan3509 的帖子

难道这个频率不是400?
发表于 2007-6-11 23:34 | 显示全部楼层

回复 #11 dengbangfei 的帖子

你不是说800吗
 楼主| 发表于 2007-6-11 23:35 | 显示全部楼层

回复 #12 zhangnan3509 的帖子

看来是你一直搞错我的意思了,我就是认为应该是400,而程序运行出来的图形显示的是800.
发表于 2007-6-11 23:37 | 显示全部楼层

回复 #13 dengbangfei 的帖子

不会吧!我这电脑上没有MATLAB,没:@L 办法试验
 楼主| 发表于 2007-6-11 23:37 | 显示全部楼层

回复 #14 zhangnan3509 的帖子

真的啊,我郁闷了好几天那!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-30 12:03 , Processed in 0.079017 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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