LMD分解程序

function [PFSI]=lmd4（x）
% 局域均值分析
% 局域均值分析
x=x‘;
c = x‘;
N = length（x）;

A= ones（1N）;
PF = [];
SI=[];
aii = 2*A;
while（1）

  si = c;
  a = 1;
  
   while（1）
    h = si;
      maxVec = [];
      minVec = [];
      
   % look for max and min point
      for i = 2: N - 1
         if h （i-1）< h （i）&h（i）>h（i+1）
            maxVec = [maxVec i]; 		
         end
         if h（i-1）>h （i）&h（i）            minVec = [minVec i]; 		
         end         
      end
      
   % check if it is residual
      if （length （maxVec） + length （minVec）） < 3
         break;
      end
           
  % handle end point 
      lenmax=length（maxVec）;
      lenmin=length（minVec）;
      %left end point
      if h（1）>0
          if（maxVec（1）              yleft_max=h（maxVec（1））;
              yleft_min=-h（1）;
          else
              yleft_max=h（1）;
              yleft_min=h（minVec（1））;
          end
      else
          if （maxVec（1）              yleft_max=h（maxVec（1））;
              yleft_min=h（1）;
          else
              yleft_max=-h（1）;
              yleft_min=h（minVec（1））;
          end
      end
      %right end point
      if h（N）>0
          if（maxVec（lenmax）