EMD中这个函数怎么理解？

dingdingysu

function [zjsy,begin]=sig_match(sig)
[indmax indmin]=extr(sig);
l=length(sig);
lmax=length(indmax);
lmin=length(indmin);
if lmax>=lmin
    k=lmin;
else
    k=lmax;
end
for i=1:(k-1)
    tx(i)=(indmax(1)*indmin(i+1)-indmin(1)*indmax(i+1))/(indmax(1)-indmin(1));
end
n=1:l;
y=interp1(n,sig,tx,'linear');
for i=1:length(y)
    e(i)=abs(sig(indmax(i+1))-sig(indmax(1)))+abs(sig(indmin(i+1))-sig(indmin(1)))+abs(y(i)-sig(1));
end
[ea,indea]=sort(e);
inde=indea(1);
zjsy=tx(inde)-1;
zjsy=floor(zjsy);
if inde<=1
    begin=1;
else
    begin=tx(inde-1);
end
begin=floor(begin);
end

杨德昌

这个应该是波形匹配函数吧！