动态时间归整.m

来自「线性预测这一术语是维纳于1947年首次提出的」· M 代码 · 共 27 行

n=size(t,1);
 m=size(r,1);
 dd=zeros(n,m);
 for i=1:n
     for j=1:m
         dd(i,j)=sum((t(i,:)-r(j,:)).^2);
     end
 end
 D=ones(n,m)*realmax;
 D(1,1)=dd(1,1);
 for i=2:n
     for j=1:m
         D1=D(i-1,j);
         if j>1
             D2=D(i-1,j-1);
         else
             D2=realmax;
         end
         if j>2
             D3=D(i-1,j-2);
         else
             D3=realmax;
         end
         D(i,j)=dd(i,j)+min([D1,D2,D3]);
     end
 end
 dist=D(n,m)