資源簡介
LMD局域均值分解matlab測試程序
寫的不完善,還存在沒有去處騎行波的問題
拿出來分享是希望有人能夠完善該程序

代碼片段和文件信息
function?[envmin?envmaxenvmoyindminindmaxindzer]?=?envelope(txINTERP)
%computes?envelopes?and?mean?with?various?interpolations
NBSYM?=?2; %?邊界延拓點數
DEF_INTERP?=?‘spline‘;
if?nargin?2
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
function?[tmintmaxxminxmax]?=?boundary_conditions(indminindmaxtxnbsym)
%?computes?the?boundary?conditions?for?interpolation?(mainly?mirror?symmetry)
lx?=?length(x);
%?判斷極值點個數
if?(length(indmin)?+?length(indmax)?3)
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?exte
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????6158??2010-12-14?19:35??lmd\envelope.m
?????文件???????5517??2011-01-15?21:17??lmd\example_lmd.asv
?????文件???????5514??2011-01-15?21:20??lmd\example_lmd.m
?????文件????????117??2010-12-15?21:35??lmd\example_shunpin.m
?????文件???????1246??2010-12-15?23:25??lmd\extrpoint.asv
?????文件???????1689??2011-01-15?15:19??lmd\hua_baoluo.m
?????文件????????687??2011-01-15?08:26??lmd\hua_fft1.m
?????文件???????1126??2011-01-13?16:29??lmd\hua_xihua.m
?????文件?????????58??2010-12-15?23:23??lmd\li
?????文件????????945??2010-12-19?11:39??lmd\lmd1.asv
?????文件????????952??2010-12-19?11:42??lmd\lmd1.m
?????文件???????1221??2011-01-15?22:32??lmd\lmd2.asv
?????文件???????1184??2011-01-15?22:34??lmd\lmd2.m
?????文件???????1788??2011-01-15?21:38??lmd\lmd2_example.m
?????文件???????1914??2010-12-16?10:45??lmd\move.asv
?????文件???????1914??2010-12-16?10:46??lmd\move.m
?????文件?????????96??2010-12-15?23:40??lmd\nengliang.m
?????文件????????185??2010-12-15?23:44??lmd\pos.m
?????文件???????1216??2010-12-15?23:43??lmd\position.m
?????文件???????1146??2010-12-16?11:00??lmd\shunpin.asv
?????文件???????1157??2010-12-16?11:04??lmd\shunpin.m
?????文件????????971??2010-12-16?10:56??lmd\smove.m
?????文件????????706??2010-12-17?08:13??lmd\test.m
?????文件???????1239??2011-01-15?22:21??lmd\testttt.m
?????文件????????416??2010-12-15?23:57??lmd\tiaozheng.m
?????文件???????6559??2011-01-15?21:22??lmd\zhaochun1.m
?????目錄??????????0??2011-03-14?11:00??lmd
-----------?---------??----------?-----??----
????????????????45721????????????????????27
............此處省略0個文件信息
評論
共有 條評論