Hilbert-Huang变换（HHT）代码

优秀论文及配套源码。Hilbert-Huang变换(HHT)是一种新的非平稳信号处理技术，该方法由经验模态 分解(EMD)与Hilbert谱分析两部分组成。任意的非平稳信号首先经过EMD方法处理后被分解为一系列具有不同特征尺度的数据序列，每一个序列称为一个固有模态函数(IMF)，然后对每个IMF分量进行Hilbert谱分析得到相应分量的Hilbert谱，汇总所有Hilbert谱就得到了原信号的谱图。该方法从本质上讲是对非平稳信号进行平稳化处理，将信号中真实存在的不同尺度波动或趋势逐级分解出来，最终用瞬时频率和能量来表征原信号的频率含量。 本文研究了基于HHT的暂态电能质量扰动检测方法，介绍了HHT的基本原理和利用HHT检测电能质量多扰动信号的实现方法。仿真试验表明该方法可以实时检测扰动的起止时刻，持续时间和扰动幅度，适用于电能质量多扰动的监测和辨识系统。

function [imfortnbits] = emd（varargin）

[xtsdsd2tolMODE_COMPLEXndirsdisplay_siftingsdtsd2trimfknbitNbItMAXITERATIONSFIXEFIXE_HMAXMODESINTERPmask] = init（varargin{:}）;

if display_sifting
  fig_h = figure;
end


%main loop : requires at least 3 extrema to proceed
while ~stop_EMD（rMODE_COMPLEXndirs） && （k < MAXMODES+1 || MAXMODES == 0） && ~any（mask）

  % current mode
  m = r;

  % mode at previous iteration
  mp = m;

  %computation of mean and stopping criterion
  if FIXE
    [stop_siftmoyenne] = stop_sifting_fixe（tmINTERPMODE_COMPLEXndirs）;
  elseif FIXE_H
    stop_count = 0;
    [stop_siftmoyenne] = stop_sifting_fixe_h（tmINTERPstop_countFIXE_HMODE_COMPLEXndirs）;
  else
    [stop_siftmoyenne] = stop_sifting（mtsdsd2tolINTERPMODE_COMPLEXndirs）;
  end

  % in case the current mode is so small that machine precision can cause
  % spurious extrema to appear
  if （max（abs（m））） < （1e-10）*（max（abs（x）））
    if ~stop_sift
      warning（‘emd:warning‘‘forced stop of EMD : too small amplitude‘）
    else
      disp（‘forced stop of EMD : too small amplitude‘）
    end
    break
  end


  % sifting loop
  while ~stop_sift && nbit
    if（~MODE_COMPLEX && nbit>MAXITERATIONS/5 && mod（nbitfloor（MAXITERATIONS/10））==0 && ~FIXE && nbit > 100）
      disp（[‘mode ‘int2str（k）‘ iteration ‘int2str（nbit）]）
      if exist（‘s‘‘var‘）
        disp（[‘stop parameter mean value : ‘num2str（s）]）
      end
      [imiM] = extr（m）;
      disp（[int2str（sum（m（im） > 0））‘ minima > 0; ‘int2str（sum（m（iM） < 0））‘ maxima < 0.‘]）
    end

    %sifting
    m = m - moyenne;

    %computation of mean and stopping criterion
    if FIXE
      [stop_siftmoyenne] = stop_sifting_fixe（tmINTERPMODE_COMPLEXndirs）;
    elseif FIXE_H
      [stop_siftmoyennestop_count] = stop_sifting_fixe_h（tmINTERPstop_countFIXE_HMODE_COMPLEXndirs）;
    else
      [stop_siftmoyennes] = stop_sifting（mtsdsd2tolINTERPMODE_COMPLEXndirs）;
    end

    % display
    if display_sifting && ~MODE_COMPLEX
      NBSYM = 2;
      [indminindmax] = extr（mp）;
      [tmintmaxmminmmax] = boundary_conditions（indminindmaxtmpmpNBSYM）;
      envminp = interp1（tminmmintINTERP）;
      envmaxp = interp1（tmaxmmaxtINTERP）;
      envmoyp = （envminp+envmaxp）/2;
      if FIXE || FIXE_H
        display_emd_fixe（tmmprenvminpenvmaxpenvmoypnbitkdisplay_sifting）
      else
        sxp=2*（abs（envmoyp））./（abs（envmaxp-envminp））;
        sp = mean（sxp）;
        display_emd（tmmprenvminpenvmaxpenvmoypsspsxpsdtsd2tnbitkdisplay_siftingstop_sift）
      end
    end

    mp = m;
    nbit=nbit+1;
    NbIt=NbIt+1;

    if（nbit==（MAXITERATIONS-1） && ~FIXE && nbit > 100）
      if exist（‘s‘‘var‘）
        warning（‘emd:warning‘[‘forced stop of sifting : too many iterations... mode ‘int2str（k）‘. stop parameter mean value : ‘num2str（s）]）
      else
        warning（‘emd:warning‘[‘forced stop of sifting : too many iterations... mode ‘int2str（k）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      17417  2010-04-26 23:21  复杂谐波\emd.m

     文件       1133  2010-04-26 23:58  复杂谐波\emdyy.m

     文件       1916  2010-06-16 20:10  复杂谐波\hhspectrum.m

     文件       2256  2010-06-06 15:34  复杂谐波\IMFeg.asv

     文件       1626  2010-06-16 20:10  复杂谐波\IMFeg.m

     文件       1942  2010-06-16 20:10  复杂谐波\toimage.m

     文件      17417  2010-04-26 23:21  复杂暂降（幅值）\emd.m

     文件       1133  2010-04-26 23:58  复杂暂降（幅值）\emdyy.m

     文件       1915  2010-06-16 20:09  复杂暂降（幅值）\hhspectrum.m

     文件       2292  2010-06-12 18:13  复杂暂降（幅值）\IMFeg.asv

     文件       1692  2010-06-18 12:25  复杂暂降（幅值）\IMFeg.m

     文件       1942  2010-06-16 20:09  复杂暂降（幅值）\toimage.m

     文件      17417  2010-04-26 23:21  复杂暂降（时刻）\emd.m

     文件       1133  2010-04-26 23:58  复杂暂降（时刻）\emdyy.m

     文件       1915  2010-06-16 20:09  复杂暂降（时刻）\hhspectrum.m

     文件       2292  2010-06-12 18:13  复杂暂降（时刻）\IMFeg.asv

     文件       1632  2010-06-17 21:09  复杂暂降（时刻）\IMFeg.m

     文件       1942  2010-06-16 20:09  复杂暂降（时刻）\toimage.m

     文件      12200  2010-06-06 19:07  复杂中断\3435567788.fig

     文件      17417  2010-04-26 23:21  复杂中断\emd.m

     文件       1133  2010-04-26 23:58  复杂中断\emdyy.m

     文件       2850  2010-04-26 23:06  复杂中断\hhspectrum.m

     文件       2299  2010-06-13 00:03  复杂中断\IMFeg.asv

     文件       2378  2010-06-13 00:04  复杂中断\IMFeg.m

     文件       3003  2007-06-27 22:12  复杂中断\toimage.m

     文件      17417  2010-04-26 23:21  间谐波\emd.m

     文件       1133  2010-04-26 23:58  间谐波\emdyy.m

     文件       1916  2010-06-16 20:07  间谐波\hhspectrum.m

     文件       2262  2010-06-07 09:41  间谐波\IMFeg.asv

     文件       2262  2010-06-07 10:12  间谐波\IMFeg.m

............此处省略24个文件信息