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

mzy

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)

[A,f,tt]=hhspectrum(imf);
[im,tt]=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

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

form

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

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

mzy

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

我爱振动

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

                               
登录/注册后可看大图

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

帖子打不开怎么回事？