clc;
clear?all;


%x=X097_DE_time（1:6000）;
t=0:0.001:1-0.001;
x=2+cos（3*t）.*cos（2*t+2*t.*t）+0.2*randn（1length（t））;


%x=（1+0.5*cos（9*pi*t））.*cos（200*pi*t+2*cos（10*pi*t））+3*cos（20*pi*t.*t+6*pi*t）;
c?=?x;
N?=?length（x）;

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

while（1）

??si?=?c;
??a?=?1;
??
???while（1）
????h?=?si;
??????maxV?=?[];
??????minV=?[];
???%?找出所有的極大值和極小值點
??????for?i?=?2:?N?-?1
?????????if?h?（i?-?1）??h?（i?+?1）
????????????maxV=?[maxV?i];?		
?????????end
?????????if?h?（i?-?1）?>?h?（i）?&?h?（i）?????????????minV?=?[minV?i];?		
?????????end?????????
??????end
???%?檢查是不是余量
??????if?（length?（maxV）?+?length?（minV））??????????break;
??????end
??%?端點的處理
??????lenmax=length（maxV）;
??????lenmin=length（minV）;
??????%左端點
??????if?h（1）>0
??????????if（maxV（1）??????????????yleft_max=h（maxV（1））;
??????????????yleft_min=-h（1）;
??????????else
??????????????yleft_max=h（1）;
??????????????yleft_min=h（minV（1））;
??????????end
??????else
??????????if?（maxV（1）??????????????yleft_max=h（maxV（1））;
??????????????yleft_min=h（1）;
??????????else
??????????????yleft_max=-h（1）;
??????????????yleft_min=h（minV（1））;
??????????end
??????end
??????%右端點
??????if?h（N）>0
??????????if（maxV（lenmax）?????????????yright_max=h（N）;
?????????????yright_min=h（minV（lenmin））;
??????????else
??????????????yright_max=h（maxV（lenmax））;
??????????????yright_min=-h（N）;
??????????end
??????else
??????????if（maxV（lenmax）??????????????yright_max=-h（N）;
??????????????yright_min=h（minV（lenmin））;
??????????else
??????????????yright_max=h（maxV（lenmax））;
??????????????yright_min=h（N）;
??????????end
??????end
??????%得到極大值點和極小值點的包絡，采用的是樣條差值進行擬合
??????maxE=spline（[1?maxV?N][yleft_max?h（maxV）?yright_max]1:N）;
??????minE=spline（[1?minV?N][yleft_min?h（minV）?yright_min]1:N）;
????mm?=?（maxE?+?minE）/2;
????aa?=?abs（maxE?-?minE）/2;
????h?=?h-mm;
????si?=?h./aa;
????a?=?a.*aa;????
????ai?=?aa;
????B?=?length（ai）;
????C?=?ones（1B）;
????bb?=?norm（ai-C）;
????if（bb?????????break;
????end?????
???end?
???pf?=?a.*si;
???PF?=?[PF;?pf];
???bbb?=?length?（maxV）?+?length?（minV）;
?%?檢查是不是余量
??????if?（length?（maxV）?+?length?（minV））??????????break;
??????end?
??c?=?c-pf;
end
line=size（PF1）;?
NN?=?length（PF（1:））;
n?=?linspace（01NN）;
figure（1）;
subplot（line+211）plot（nx）ylabel（‘X（t）‘）;
subplot（line+212）plot（nPF（1:））ylabel（‘PF_1（t）‘）;
subplot（line+213）plot（nPF（2:））ylabel（‘PF_2（t）‘）;
subplot（line+214）plot（nPF（3:））ylabel（‘PF_3（t）‘）;
subplot（line+215）plot（nc）ylabel（‘u（t）‘）;
%subplot（line+216）plot（nPF（5:））ylabel（‘PF_5（t）‘）;
%subplot（line+217）plot（nc）ylabel（‘u（t）‘）;
%subplot（line+118）plot（nPF（7:）*1000）ylabel（‘PF_7（t）‘）;
%subplot（line+119）plot（nPF（8:）*1000）ylabel（‘PF_8（t）‘）;
%subplot（line+1110）plot（nPF（9:）*1000）ylabel（‘PF_9（t）‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3034??2015-04-18?16:02??LMD1.m

-----------?---------??----------?-----??----

?????????????????3034????????????????????1