%參考程序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