边际谱 问题求助

zhangnan3509

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

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

dengbangfei

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

zhlong

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

zhangnan3509

这能说明什么？

zhangnan3509

中间少了一步：

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

zhangnan3509

你看看这效果，很郁闷的

dengbangfei

请问中间少了哪一步呢?

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

zhangnan3509

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

破凰

原帖由 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

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

dengbangfei

谢谢你!