function?[envmin?envmaxenvmoyindminindmaxindzer]?=?envelope（txINTERP）
%computes?envelopes?and?mean?with?various?interpolations

NBSYM?=?2;?%?邊界延拓點數
DEF_INTERP?=?‘spline‘;


if?nargin?x?=?t;
t?=?1:length（x）;
INTERP?=?DEF_INTERP;
end

if?nargin?==?2
if?ischar（x）
INTERP?=?x;
x?=?t;
t?=?1:length（x）;
end
end

if?~ischar（INTERP）
error（‘interp?parameter?must?be?‘‘linear‘‘‘‘?‘‘cubic‘‘?or?‘‘spline‘‘‘）
end

if?~any（strcmpi（INTERP{‘linear‘‘cubic‘‘spline‘}））
error（‘interp?parameter?must?be?‘‘linear‘‘‘‘?‘‘cubic‘‘?or?‘‘spline‘‘‘）
end

if?min（[size（x）size（t）]）?>?1
error（‘x?and?t?must?be?vectors‘）
end
s?=?size（x）;
if?s（1）?>?1
x?=?x‘;
end
s?=?size（t）;
if?s（1）?>?1
t?=?t‘;
end
if?length（t）?~=?length（x）
error（‘x?and?t?must?have?the?same?length‘）
end

lx?=?length（x）;
[indminindmaxindzer]?=?extr（xt）;


%boundary?conditions?for?interpolation

[tmintmaxxminxmax]?=?boundary_conditions（indminindmaxtxNBSYM）;

%?definition?of?envelopes?from?interpolation

envmax?=?interp1（tmaxxmaxtINTERP）;?
envmin?=?interp1（tminxmintINTERP）;

if?nargout?>?2
envmoy?=?（envmax?+?envmin）/2;
end
end
function?[tmintmaxxminxmax]?=?boundary_conditions（indminindmaxtxnbsym）
%?computes?the?boundary?conditions?for?interpolation?（mainly?mirror?symmetry）


lx?=?length（x）;
%?判斷極值點個數
if?（length（indmin）?+?length（indmax）?error（‘not?enough?extrema‘）
end
%?插值的邊界條件
if?indmax（1）?if?x（1）?>?x（indmin（1））%?以第一個極大值為對稱中心
lmax?=?fliplr（indmax（2:min（endnbsym+1）））;
lmin?=?fliplr（indmin（1:min（endnbsym）））;
lsym?=?indmax（1）;
else%?如果第一個采樣值小于第一個極小值，則將認為該值是一個極小值，以該點為對稱中心
lmax?=?fliplr（indmax（1:min（endnbsym）））;
lmin?=?[fliplr（indmin（1:min（endnbsym-1）））1];
lsym?=?1;
end
else

if?x（1）?lmax?=?fliplr（indmax（1:min（endnbsym）））;
lmin?=?fliplr（indmin（2:min（endnbsym+1）））;
lsym?=?indmin（1）;
else%?如果第一個采樣值大于第一個極大值，則將認為該值是一個極大值，以該點為對稱中心
lmax?=?[fliplr（indmax（1:min（endnbsym-1）））1];
lmin?=?fliplr（indmin（1:min（endnbsym）））;
lsym?=?1;
end
end
%?序列末尾情況與序列開頭類似
if?indmax（end）?if?x（end）?rmax?=?fliplr（indmax（max（end-nbsym+11）:end））;
rmin?=?fliplr（indmin（max（end-nbsym1）:end-1））;
rsym?=?indmin（end）;
else
rmax?=?[lxfliplr（indmax（max（end-nbsym+21）:end））];
rmin?=?fliplr（indmin（max（end-nbsym+11）:end））;
rsym?=?lx;
end
else
if?x（end）?>?x（indmin（end））
rmax?=?fliplr（indmax（max（end-nbsym1）:end-1））;
rmin?=?fliplr（indmin（max（end-nbsym+11）:end））;
rsym?=?indmax（end）;
else
rmax?=?fliplr（indmax（max（end-nbsym+11）:end））;
rmin?=?[lxfliplr（indmin（max（end-nbsym+21）:end））];
rsym?=?lx;
end
end
%?將序列根據對稱中心，鏡像到兩邊
tlmin?=?2*t（lsym）-t（lmin）;
tlmax?=?2*t（lsym）-t（lmax）;
trmin?=?2*t（rsym）-t（rmin）;
trmax?=?2*t（rsym）-t（rmax）;

%?in?case?symmetrized?parts?do?not?extend?enough%?如果對稱的部分沒有足夠的極值點
if?tlmin（1）?>?t（1）?|?tlmax（1）?>?t（1）%?對折后的序列沒有超出原序列的范圍
if?lsym?==?indmax（1）
lmax?=?fliplr（indmax（1:min（endnbsym）））;
else

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

?????文件???????3957??2013-01-11?10:32??LMD\envelope.m

?????文件???????1705??2013-01-11?10:33??LMD\extr.m

?????文件?????????58??2010-12-15?23:23??LMD\link.m

?????文件????????843??2013-01-11?10:29??LMD\lmd1.m

?????文件???????1916??2013-01-14?17:17??LMD\move.m

?????文件????????100??2013-01-11?10:29??LMD\nengliang.m

?????文件????????189??2013-01-11?10:34??LMD\pos.m

?????文件????????991??2013-01-11?10:35??LMD\position.m

?????文件????????971??2010-12-16?10:56??LMD\smove.m

?????文件????????343??2013-01-11?10:35??LMD\Untitled8.m

?????文件????????386??2013-01-11?10:30??LMD\zhaochun1.m

?????目錄??????????0??2018-08-15?15:55??LMD

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

????????????????11459????????????????????12