请各位大虾过来瞅瞅哦：HHT变换中的 瞬时能量谱

zxcv3339

灰太狼1985应该是对的，因为论文公式里有个dw,乘以2*pi是干什么，不明白

zxcv3339

受益良多，特此分享！
[A,fa,tt]=hhspectrum(cc);                                
[E,tt1]=toimage(A,fa,tt,length(tt));   
E1=flipud(E);   
for k=1:size(E1,2)
IE(k)=sum(E1(:,k).^2)*2*pi*fs/(2*N);
end
n=(0:N-3)/N*t;
figure(1)
plot(n,IE)