声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1827|回复: 7

[HHT] 求助:瞬时频率的算法是如何实现的?

[复制链接]
发表于 2009-4-25 20:31 | 显示全部楼层 |阅读模式

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

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

x
下面是用hilbert变换求瞬时频率的MATLAB程序,看了好久没懂,请教高手指导下这个算法是如何实现的。
% this is a function to calculate instantaneous based on EMD method
%
%   function omega = ifndq(vimf, dt)
%
% INPUT:   
%          vimf:        an IMF;
%          dt:          time interval of the imputted data    输入数据的时间间隔
% OUTPUT:
%          omega:       instantanesous frequency, which is 2*PI/T, where T
%          瞬时频率
%                       is the period of an oscillation
%
% References can be found in the "Reference" section.
%
% The code is prepared by Zhaohua Wu. For questions, please read the "Q&A" section or
% contact
%   zhwu@cola.iges.org
%
function omega = ifndq(vimf, dt)
Nnormal=5;
rangetop=0.90;
vlength = max( size(vimf) );
vlength_1 = vlength -1;
for i=1:vlength,
    abs_vimf(i)=vimf(i);
    if abs_vimf(i) < 0
        abs_vimf(i)=-vimf(i);
    end
end
for jj=1:Nnormal,
    [spmax, spmin, flag]=extrema(abs_vimf);
    dd=1:1:vlength;
    upper= spline(spmax(:,1),spmax(:,2),dd);
    for i=1:vlength,
        abs_vimf(i)=abs_vimf(i)/upper(i);
    end
end
for i=1:vlength,
    nvimf(i)=abs_vimf(i);
    if vimf(i) < 0;
        nvimf(i)=-abs_vimf(i);
    end
end
for i=1:vlength,
    dq(i)=sqrt(1-nvimf(i)*nvimf(i));
end
for i=2:vlength_1,
    devi(i)=nvimf(i+1)-nvimf(i-1);
    if devi(i)>0 & nvimf(i)<1
        dq(i)=-dq(i);
    end
end
rangebot=-rangetop;     
for i=2:(vlength-1),
    if nvimf(i)>rangebot & nvimf(i) < rangetop
        omgcos(i)=abs(nvimf(i+1)-nvimf(i-1))*0.5/sqrt(1-nvimf(i)*nvimf(i));
    else
        omgcos(i)=-9999;
    end
end
omgcos(1)=-9999;
omgcos(vlength)=-9999;
jj=1;
for i=1:vlength,
    if omgcos(i)>-1000
        ddd(jj)=i;
        temp(jj)=omgcos(i);
        jj=jj+1;
    end
end
temp2=spline(ddd,temp,dd);
omgcos=temp2;

for i=1:vlength,
    omega(i)=omgcos(i);
end
pi2=pi*2;
omega=omega/dt;
回复
分享到:

使用道具 举报

 楼主| 发表于 2009-4-26 13:09 | 显示全部楼层

高手何在?

请大家帮忙啊
 楼主| 发表于 2009-4-27 21:13 | 显示全部楼层

高手何在?

等的周期好长啊,高手何在?
发表于 2009-4-27 21:28 | 显示全部楼层
容我说句话, 高手那有空当家教!
 楼主| 发表于 2009-4-27 22:59 | 显示全部楼层

回复 地板 ChaChing 的帖子

你啥意思?你算个球啊
发表于 2009-4-28 11:15 | 显示全部楼层
楼主先勿生气! 楼主的帖子连续出现了三天, 楼主的急可想而知!
楼主的求知欲望值得嘉许, 但方式的确值得商榷!
若能先搜索并试下, 再询问重点(不是全部解说), 或许楼主能够得到更多!
 楼主| 发表于 2009-4-28 18:37 | 显示全部楼层

回复 6楼 ChaChing 的帖子

多谢,我知道了,我也是想把话朝简单方向说的,问题是这个程序其实就那么一两个公式,再抽象出来的话就没东西了。下来我自己再看看吧,不打扰大家了。
发表于 2010-5-18 15:49 | 显示全部楼层
上述程序中
for i=2:vlength_1,
    devi(i)=nvimf(i+1)-nvimf(i-1);
    if devi(i)>0 & nvimf(i)<1
        dq(i)=-dq(i);
    end
end
这段,为甚取负的平方根呢?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-13 18:46 , Processed in 0.065144 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表