端点极值延拓

%参考程序G.Rillng的 emd分解   %Bivariate Empirical Mode
%Decomposition，得出imf，，其中t0由信号具体决定其长度%图3-3（对应11行）与图3-7其中（对应12行）
clc
clear all
dt=0.001;
 tic;
t=0:dt:1;
%dt=1;
%x = sin（2*pi*10*t） + sin（2*pi*20*t）;% + sin（t）.*exp（5*t）/10;
%x=cos（100*pi*t）+2*cos（50*pi*t+0.5*sin（10*pi*t））;
%x=cos（25*pi*t）+2*cos（12.5*pi*t+0.125*sin（2.5*pi*t））;图3-3
x=5*sin（2*pi*10*t）+5*sin（2*pi*35*t）;%图3-7
%plot（x）
figure
%x = sin（2*pi*10*t） + sin（2*pi*50*t） +
%sin（2*pi*100*t）;%+sin（t）.*exp（10*t）/50;
%[x fs nbit]=wavread（‘D:\Desktop\plot_hht\Hum.wav‘）;
pp=0;
%for ff=0.04:0.0025:0.1
%x=0.5*sin（2*pi*0.1*t）+sin（2*pi*ff*t）;
pp=pp+1;
num=length（x）;
z=max（x）;
   h=x;
   imf=[];
   n=1;
   xrest=x;
   T=0:0.001:2;
while （max（x）>z*0.025）  %x始终是剩余的函数
xa=diff（h）;
xb=xa（1:num-2）;
xc=xa（2:num-1）;%微分后的xa至多有num-1个点此处为 参考Gabriel Rilling，ON EMPIRICAL MODE DECOMPOSITIONAND ITS ALGORITHMS
nummin=[];
nummax=[];
nummin=find（xb.*xc<0 & xb<0）;
nummax=find（xb.*xc<0 & xc<0）;
if length（nummin）<=0 || length（nummax