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

回复 #15 dengbangfei 的帖子

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

[ 本帖最后由 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

回复 #18 zhlong 的帖子

这能说明什么?

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

回复 #1 dengbangfei 的帖子

中间少了一步:

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

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

回复 #18 zhlong 的帖子

你看看这效果,很郁闷的

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

回复 #20 zhangnan3509 的帖子

请问中间少了哪一步呢?

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

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

回复 #22 dengbangfei 的帖子

对不起,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);
%%%%%%%%%%%%%求边际谱
=hhspectrum(y1);
=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

回复 #23 zhangnan3509 的帖子

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

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

回复 #24 破凰 的帖子

谢谢你!
页: 1 [2]
查看完整版本: 边际谱 问题求助