hht 的边际谱和能量谱计算

blusun_2008

我用rilling的HHT程序代码算出了各个IMF，然后求出hht谱图，用时域图表示出来，用这样的程序，不知对不对？
[imf,ort,nbits] = emd(pm);
   [A,f,tt] = hhspectrum(imf);
   subplot(2,2,1);
   plot(tt,f);

现在我想求它的边际谱和能量谱，原理是对hht谱的时间积分和hht谱的平方后频率积分，可是我用这个函数：ie=cumtrapz(f,A);来对hht谱积分，出现错误的，那应该怎样求它的边际谱和能量谱呢？谢谢指教了！

blusun_2008

http://forum.vibunion.com/forum/thread-50010-1-4.html
到这个网址上看看，有相关信息哦！