%%%局域均值分解
%%x原始信號

clc
clear
close?all
fs=2000;
t=0:1/fs:1;

%?x=15*（1+cos（40*pi*t））.*cos（600*pi*t）+5*（1+cos（40*pi*t））.*cos（200*t*pi）;
%x=cos（2*pi*30*t?+0.1*sin（2*pi*10*t））+sin（2*80*t）;
%x=5*（1+cos（40*pi*t））.*cos（150*pi*t）+5*（1+cos（40*pi*t））.*cos（100*t*pi）;
%plot（tx）;
%xlabel（‘時間t/s‘）
%ylabel（‘幅值‘）
x=xlsread（‘Book2.xlsx‘）;
figure（1）

plot（tx）;
xlabel（‘時間t/s‘）
ylabel（‘幅值‘）
figure（3）
y2=x;
L=length（y2）;
NFFT?=?2^nextpow2（L）;
Y?=?fft（y2NFFT）/L;
f?=?fs/2*linspace（01NFFT/2）;
plot（f2*abs（Y（1:NFFT/2）））
xlabel（‘頻率f/Hz‘）
ylabel（‘幅值‘）





%規定輸入信號為行向量
if?size（x1）~=1
????x?=?x‘;
end
c?=?x;
N?=?length（x）;


PF?=?[];
A1=[];
Si=[];
while（1）?%loop?1
????
????a?=?1;
????
????while（1）?%loop?2
????????h?=?c;
????????
????????maxVec?=?[];
????????minVec?=?[];
????????
????????%?look?for?max?and?min?point?找到h中的局部極大值和極小值
????????for?i?=?2:?N?-?1
????????????if?h?（i?-?1）??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??判斷是否滿足loop2的終止條件
????????if?（length?（maxVec）?+?length?（minVec））?????????????break;
????????end
????????
????????%?handle?end?point

????????%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（end）????????????????yright_max=h（N）;
????????????????yright_min=h（minVec（end））;
????????????else
?????????