w.m

来自「用MATLAB实现模拟立方体的重构」· M 代码 · 共 17 行

function w=W(x1,x2)

%计算绝对二次曲线的像的函数
% x1,x2分别是三正交方向消影点坐标矩阵

A(1,:)=[x1(1,1)*x1(1,2) x1(2,1)*x1(1,2)+x1(1,1)*x1(2,2) x1(2,1)*x1(2,2) x1(1,1)+x1(1,2) x1(2,1)+x1(2,2) 1];
A(2,:)=[x1(1,1)*x1(1,3) x1(2,1)*x1(1,3)+x1(1,1)*x1(2,3) x1(2,1)*x1(2,3) x1(1,1)+x1(1,3) x1(2,1)+x1(2,3) 1];
A(3,:)=[x1(1,2)*x1(1,3) x1(2,2)*x1(1,3)+x1(1,2)*x1(2,3) x1(2,2)*x1(2,3) x1(1,2)+x1(1,3) x1(2,2)+x1(2,3) 1];
A(4,:)=[x2(1,1)*x2(1,2) x2(2,1)*x2(1,2)+x2(1,1)*x2(2,2) x2(2,1)*x2(2,2) x2(1,1)+x2(1,2) x2(2,1)+x2(2,2) 1];
A(5,:)=[x2(1,1)*x2(1,3) x2(2,1)*x2(1,3)+x2(1,1)*x2(2,3) x2(2,1)*x2(2,3) x2(1,1)+x2(1,3) x2(2,1)+x2(2,3) 1];
A(6,:)=[x2(1,2)*x2(1,3) x2(2,2)*x2(1,3)+x2(1,2)*x2(2,3) x2(2,2)*x2(2,3) x2(1,2)+x2(1,3) x2(2,2)+x2(2,3) 1];
[U,D,V]=svd(A);
c=V(:,6)/V(6,6);
w=[c(1) c(2) c(4);...%绝对二次曲线的像
   c(2) c(3) c(5);...
   c(4) c(5) c(6)];