HHT中平稳度的定义

rogen

在论坛中搜了好久，发现，只有对程序的一些基本问题提问的，而对于在huang的经典文章中提到的平稳度的概念则无人问津，提出这个概念，希望大家探讨！

rogen

我对其进行了分析，但做出来的结果好像不太合适！具体代码如下：

1. % Marginal hilbert spectrum
   
2. tl = length(im)*Ts;
   
3. for k = 1:size(im,1)
   
4.     mhs(k) = sum(im(k,:))*Ts;                % marginal hilbert spectrum
   
5.     mmhs(k) = sum(im(k,:))*Ts*1/tl;          % Mean marginal hilbert spectrum
   
6.     %mhset(k) = sum(im(k,:).^2)*1/Ts*1/tl;   % Energy
   
7.     mhset(k) = sum(im(k,:).^2*Ts);           % Accumulate Energy
   
8.     mhsef(k) = sum(im(k,:).^2*fs/NN/2);      % Instantaneous Energy
   
9.    
   
10.     Ds(k)=sum((1-(im(k,:)/mmhs(k))).^2)/tl;  %Degrade of stationarity
    
11.    
    
12. end
    
13. 
    

复制代码

rogen

平稳度从公式来理解，应该是：每个时间点的幅值远离平均值的程度并对时间球平均得到的结果，其中幅值采用的是im用的是toimage出来的结果，我的理解也应该是对的，平均值是对边际谱对时间的平均，也是对的，程序应该是没有问题的，对于同样的程序对x=cos(2*pi*t/0.1)+cos(2*pi*t/0.2);做平稳度，结果如下图：可以结果就是不让人满意！

静强

您好，我最近也在用hht计算平稳度，遇到的问题和你一样，苦恼很多天也没结果，不知道您怎么解决的谢谢

ben1984

DS(K)=...是否还得乘以Ts？