HHT中平稳度的定义
在论坛中搜了好久，发现，只有对程序的一些基本问题提问的，而对于在huang的经典文章中提到的平稳度的概念则无人问津，提出这个概念，希望大家探讨！ 我对其进行了分析，但做出来的结果好像不太合适！具体代码如下：% Marginal hilbert spectrum
tl = length(im)*Ts;
for k = 1:size(im,1)
mhs(k) = sum(im(k,:))*Ts; % marginal hilbert spectrum
mmhs(k) = sum(im(k,:))*Ts*1/tl; % Mean marginal hilbert spectrum
%mhset(k) = sum(im(k,:).^2)*1/Ts*1/tl; % Energy
mhset(k) = sum(im(k,:).^2*Ts); % Accumulate Energy
mhsef(k) = sum(im(k,:).^2*fs/NN/2); % Instantaneous Energy
Ds(k)=sum((1-(im(k,:)/mmhs(k))).^2)/tl;%Degrade of stationarity
end
平稳度从公式来理解，应该是：每个时间点的幅值远离平均值的程度并对时间球平均得到的结果，其中幅值采用的是im用的是toimage出来的结果，我的理解也应该是对的，平均值是对边际谱对时间的平均，也是对的，程序应该是没有问题的，对于同样的程序对x=cos(2*pi*t/0.1)+cos(2*pi*t/0.2);做平稳度，结果如下图：可以结果就是不让人满意！
您好，我最近也在用hht计算平稳度，遇到的问题和你一样，苦恼很多天也没结果，不知道您怎么解决的谢谢 DS(K)=...是否还得乘以Ts？
页:
[1]