|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
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 |
|