EMD 算法MATLAB 程序

数据包里有原始数据，打开就可以看到效果。先用采集卡采集信号，得到信号后做EMD去趋势处理，去直流处理，然后时域分析波形（脉冲宽度等），频域分析（峰值频率，低频，中心频率等）；

function [allmode] = eemd（Y NoiseLevel NE numImf varargin）
% fast EMD/EEMD/CEEMD code:
%  Copyright （C） RCADA National Central University; 2013
%  ***** For Educational and Academic Purposes Only ********* 
%  Author : Yung-Hung Wang RCADA NCU
%
% Please Cite: 
% Y. H. Wang C. H. Yeh H. W. V. Young K. Hu and M. T. Lo
% “On the computational complexity of the empirical mode decomposition
% algorithm“
% Physica A: Statistical Mechanics and its Applications vol. 400Issue 15 pp. 159-167 2014
%
% Usage: 
% （A） [allmode]=eemd（Y NoiseLevel NE numImf）
% （B） [allmode]=eemd（xNoiseLevelNEnumImfrunCEEMDmaxSifttypeSplinetoModifyBCrandTypeseedNocheckInput）;
%---------------------------------------------------------------
% INPUTS : 
%        Y: input signal
%        NoiseLevel: noise level
%        NE: ensemble number; if NE=1 and NoiseLevel = 0 then run EMD [2]
%        numImf: number of prescribed imf; if it is less than zeroit will be log2（n） n=data length 
% 
%-------------Additional Input Properties-----------------------------------------
% runCEEMD: default=0; 0: run EEMD [3]; 1: CEEMD [4] add anti-phase noise so input signal is reconstructed; 
% maxSift:default = 10; sifting iteration number; 
% typeSpline: default=2; 1: clamped spline; 2: not a knot spline; 3: natural cubic spline;
% toModifyBC: default=1; 0: None ; 1: modified linear extrapolation; 2: Mirror Boundary
% randType: fefault=2; 1: uniformly distributed random noise; 2: gaussian white noise
% seedNo: default=1; random seed number in generating white noise; The seed number must be an integer between 0 and 2^32 - 1
% checkSignal: default=1;  0=> dont check  the input signal;  1: check there is any  NaN or Infinity in input;
%
% OUTPUTS :
%  allmode[m][t]: returned imfs; m=imf index t=time index
%
% References
%    [1] Y. H. Wang C. H. Yeh H. W. V. Young K. Hu and M. T. Lo 
%        “On the computational complexity of the empirical mode decomposition algorithm“
%        Physica A: Statistical Mechanics and its Applications vol. 400 Issue 15 pp. 159-167 2014
%    [2] N. E. Huang Z. Shen S. R. Long M. L. Wu H. H. Shih Q. Zheng N. C. Yen C. C. Tung and H. H. Liu 
%        “The Empirical Mode Decomposition and Hilbert spectrum for nonlinear and non-stationary time series analysis“
%        Proc. Roy. Soc. London A Vol. 454 pp. 903995 1998.
%    [3] Z. Wu and N. E. Huang “Ensemble Empirical Mode Decomposition: A
%        Noise-Assisted Data Analysis Method“ Advances in Adaptive Data Analysis vol. 1 pp. 1-41 2009.
%    [4] J. R. Yeh J. S. Shieh and N. E. Huang “Complementary ensemble empirical
%        mode decomposition: A novel noise enhanced data analysis method“ 
%        Advances in Adaptive Data Analysis vol. 2 Issue 2 pp. 135156 2010.
%
%


fprintf（‘Copyright （C） RCADA NCU Taiwan.\n‘）;

allmode = [];
%verifyCode = 20041017; % For verification of emd

[YNoiseLevelNE

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

     文件     112697  2014-04-01 17:02  FEEMD\BFVL.mat

     文件       7870  2014-09-04 17:56  FEEMD\eemd.m

     文件      22140  2014-09-04 16:52  FEEMD\emd.mexa64

     文件      16896  2014-08-01 13:27  FEEMD\emd.mexw32

     文件      21504  2014-09-04 17:26  FEEMD\emd.mexw64

     文件       2551  2015-06-15 09:38  FEEMD\example_eemd.asv

     文件       2541  2015-06-22 11:00  FEEMD\example_eemd.m

     文件     756829  2014-04-01 17:01  FEEMD\ref.pdf

     文件      53820  2015-06-06 22:59  FEEMD\直探头 BFA086 20150605 202719 衰减200m.txt

     文件      53908  2015-06-16 22:16  FEEMD\直探头 EDH004 20150616 220903 50m.txt

     文件      53947  2015-06-16 22:10  FEEMD\直探头 EDH004 20150616 221000.txt

     文件      54282  2015-06-18 22:01  FEEMD\直探头 EDH004 20150618 215701.txt

     文件      54324  2015-06-18 18:10  FEEMD\直探头 EDH019 20150618 180538.txt

     目录          0  2015-06-22 10:56  FEEMD

----------- ---------  ---------- -----  ----

              1213309                    14