syms x1 x2 y1 y2 m1 m2 n1 n2 t a b;
m=[cos(t), -sin(t);sin(t), cos(t)]
k=[x1,y1;x2,y2];
kk=[m1,n1;m2,n2];
f=k-m*kk-[a,b];
ssd=sum(sum(f.*f))
diff(ssd,a)
diff(ssd,b)
diff(ssd,t)
solve([diff(ssd,a), diff(ssd,b), diff(ssd,t)],[a, b, t])
本来写以上代码可以解决该问题的。但是分析一下,发现,中点平移到中点,然后旋转到两直线方向重合就解决了。基于我们写的主成分分析PCA入门上一篇。