Hilbert能量谱求助？

Mrchen

[A,fa,tt]=hhspectrum(imf);
[E,tt1]=toimage(A,fa,tt,length(tt));
for k=1:size(E,1)
    bjp(k)=sum(E(k,:)*1/fs;  
end
这是求得边际谱，那么Hilbert能量谱是不是sum((E(k,:).^2)*1/fs;