EMD中这个函数怎么理解？
function =sig_match(sig)=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
=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 这个应该是波形匹配函数吧！
页:
[1]