HHT中平稳度的定义

rogen · 发表于 2012-5-29 15:56

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

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

rogen · 发表于 2012-5-29 16:02

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

% 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

复制代码

rogen · 发表于 2012-5-29 17:00

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

静强 · 发表于 2013-4-17 14:39

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

ben1984 · 发表于 2013-11-6 14:44

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

账号		自动登录	找回密码
密码			我要加入

[HHT] HHT中平稳度的定义

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。