mzy 发表于 2007-11-13 16:58

我想求EMD后的边际谱，我用了两个论坛的程序拼到一块儿，可就是不对

clear;
t=1:500;
t=t*1/2000;
x=sin(2*pi*20*t);
y=0.4*sin(2*pi*40*t+140);
z=x+y;
plot(t,z)

imf=emd(z);
emd_visu(z,1:length(z),imf)

=hhspectrum(imf);
=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))

%%%%%%%%%%%%%求边际谱
fs=1000; %fs为采样频率;
N=1000;   %采样点数
im=flipud(im);
for k=1:size(im,1)
bjp(k)=sum(im(k,:))*1/fs;
end

plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
------------------------------------------------------------
无法得到边际谱，各位请指点一下迷津，谢谢！！

zhlong 发表于 2007-11-13 18:35

边际谱程序用得不对，请检查！

form 发表于 2007-11-15 10:40

http://forum.vibunion.com/forum/viewthread.php?tid=37688
这是以前的帖子

[ 本帖最后由 form 于 2007-11-15 11:15 编辑 ]

mzy 发表于 2007-11-15 11:31

多谢FORM,那个贴子我以前看过，可用的时候就又忘了它了，晕

我爱振动 发表于 2014-2-9 18:38

form 发表于 2007-11-15 10:40 static/image/common/back.gif
http://www.chinavib.com/forum/viewthread.php?tid=37688
这是以前的帖子


帖子打不开怎么回事？
页: [1]
查看完整版本: 我想求EMD后的边际谱，我用了两个论坛的程序拼到一块儿，可就是不对