91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 6KB
    文件類型: .rar
    金幣: 2
    下載: 0 次
    發布日期: 2021-05-19
  • 語言: 其他
  • 標簽: LMD代碼??

資源簡介

很簡單的局部均值分解代碼,可能不夠完善,需要修改能夠使用

資源截圖

代碼片段和文件信息

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


評論

共有 條評論

相關資源