边际谱是否正确呢？

nkdtxf

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;
tspan=0.4;
plot(t,z)

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

%%%%%%%%%%%%%求边际谱
figure;
fs=1000; %fs为采样频率;
N=1000;     %采样点数
[A,fa,tt]=hhspectrum(z);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E);
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs*1/tspan;  
end
H=size(E,1);
f=(0:H-1)/H*(fs/2);
plot(f,bjp);
上面的程序得到的边际谱是否正确呢，我综合论坛的程序做的。谢谢！

zhangnan3509

做出结果来了吗 把效果图贴出来