边际谱问题求助

zhangnan3509 · 发表于 2007-6-11 23:40

明天我去试验室看看！破凰的应该没什么问题

[ 本帖最后由 zhlong 于 2007-6-12 16:32 编辑 ]

dengbangfei · 发表于 2007-6-11 23:42

还有就是,为什么只是任意改变信号的幅值大小，出来的频谱图，频率都不一样，纳闷ing。

zhlong · 发表于 2007-6-11 23:50

胡劲松基于EMD和HT的旋转机械振动信号时频分析振动、测试与诊断 2004。6 vol24.N2

zhangnan3509 · 发表于 2007-6-11 23:53

这能说明什么？

zhangnan3509 · 发表于 2007-6-12 09:18

中间少了一步：

[ 本帖最后由 zhangnan3509 于 2007-6-12 09:34 编辑 ]

zhangnan3509 · 发表于 2007-6-12 10:29

你看看这效果，很郁闷的

dengbangfei · 发表于 2007-6-12 12:36

请问中间少了哪一步呢?

[ 本帖最后由 zhlong 于 2007-6-12 15:42 编辑 ]

zhangnan3509 · 发表于 2007-6-12 17:08

对不起，dengbangfei 是我们疏忽了！这么多天了也没注意到问题，今天和zhlong讨论了很久，破凰也做出了修改，中间那部分是按照zhlong的意见，把fs/2,谢谢的细心和对论坛的支持，为我的武断向你道歉:handshake

破凰 · 发表于 2007-6-12 17:10

本帖最后由 wdhd 于 2016-9-6 14:45 编辑

原帖由 dengbangfei 于 2007-6-12 12:36 发表
请问中间少了哪一步呢?

帮你改了下程序
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=(0:N-3)/N*fs/2;
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
另外那篇贴子的程序有点小错。

dengbangfei · 发表于 2007-6-12 18:46

呵呵,版主们太客气了,真的很谢谢你们最近一段时间内不厌其烦的解答我的问题,谢谢了!

dengbangfei · 发表于 2007-6-12 18:48

谢谢你!

账号		自动登录	找回密码
密码			我要加入

[HHT] 边际谱问题求助

回复 #15 dengbangfei 的帖子

回复 #18 zhlong 的帖子

回复 #1 dengbangfei 的帖子

回复 #18 zhlong 的帖子

回复 #20 zhangnan3509 的帖子

回复 #22 dengbangfei 的帖子

评分

回复 #23 zhangnan3509 的帖子

回复 #24 破凰的帖子

[HHT] 边际谱 问题求助

回复 #15 dengbangfei 的帖子

回复 #18 zhlong 的帖子

回复 #1 dengbangfei 的帖子

回复 #18 zhlong 的帖子

回复 #20 zhangnan3509 的帖子

回复 #22 dengbangfei 的帖子

评分

回复 #23 zhangnan3509 的帖子

回复 #24 破凰 的帖子

[HHT] 边际谱问题求助

回复 #24 破凰的帖子