问题描述:
Matlab中使用solve命令出现sym问题
在相机成像的模型中
需要用matlab解算一个方程组:
S=solve('Xw1*r11+Yw1*r12+A(1)*Xw1*r31+A(1)*Yw1*r32+t1+A(1)*t3=0',...
'Xw1*r21+Yw1*r22+B(1)*Xw1*r31+B(1)*Yw1*r32+t2+B(1)*t3=0',...
'Xw2*r11+Yw2*r12+A(2)*Xw2*r31+A(2)*Yw2*r32+t1+A(2)*t3=0',...
'Xw2*r21+Yw2*r22+B(2)*Xw2*r31+B(2)*Yw2*r32+t2+B(2)*t3=0',...
'Xw3*r11+Yw3*r12+A(3)*Xw3*r31+A(3)*Yw3*r32+t1+A(3)*t3=0',...
'Xw3*r21+Yw3*r22+B(3)*Xw3*r31+B(3)*Yw3*r32+t2+B(3)*t3=0',...
'Xw4*r11+Yw4*r12+A(4)*Xw4*r31+A(4)*Yw4*r32+t1+A(4)*t3=0',...
'Xw4*r21+Yw4*r22+B(4)*Xw4*r31+B(4)*Yw4*r32+t2+B(4)*t3=0',...
'r11^2+r21^2+r31^2=1',...
'r21^2+r22^2+r32^2=1',...
'r11*r12+r21*r22+r31*r32=0');
其中向量A=[0.2252 0.2453 -0.1968 -0.1834]
B=[0.1665 -0.1506 -0.1757 0.0798]
Xw1=0;Xw2=0;Xw3=1;Xw4=1;
Yw1=0.6;Yw2=0;Yw3=0;Yw4=0.6; 这些都为已知量
需要通过8个一次线性方程和3个二次方程共11个方程来解算9个未知量:r11 r12 r21 r22 r31 r32 t1 t2 t3
而直接使用solve命令解算11个方程 CPU%100很长时间都没算出来
如果只用前8个一次线性方程解算 出现了 带sym的解
在相机成像的模型中
需要用matlab解算一个方程组:
S=solve('Xw1*r11+Yw1*r12+A(1)*Xw1*r31+A(1)*Yw1*r32+t1+A(1)*t3=0',...
'Xw1*r21+Yw1*r22+B(1)*Xw1*r31+B(1)*Yw1*r32+t2+B(1)*t3=0',...
'Xw2*r11+Yw2*r12+A(2)*Xw2*r31+A(2)*Yw2*r32+t1+A(2)*t3=0',...
'Xw2*r21+Yw2*r22+B(2)*Xw2*r31+B(2)*Yw2*r32+t2+B(2)*t3=0',...
'Xw3*r11+Yw3*r12+A(3)*Xw3*r31+A(3)*Yw3*r32+t1+A(3)*t3=0',...
'Xw3*r21+Yw3*r22+B(3)*Xw3*r31+B(3)*Yw3*r32+t2+B(3)*t3=0',...
'Xw4*r11+Yw4*r12+A(4)*Xw4*r31+A(4)*Yw4*r32+t1+A(4)*t3=0',...
'Xw4*r21+Yw4*r22+B(4)*Xw4*r31+B(4)*Yw4*r32+t2+B(4)*t3=0',...
'r11^2+r21^2+r31^2=1',...
'r21^2+r22^2+r32^2=1',...
'r11*r12+r21*r22+r31*r32=0');
其中向量A=[0.2252 0.2453 -0.1968 -0.1834]
B=[0.1665 -0.1506 -0.1757 0.0798]
Xw1=0;Xw2=0;Xw3=1;Xw4=1;
Yw1=0.6;Yw2=0;Yw3=0;Yw4=0.6; 这些都为已知量
需要通过8个一次线性方程和3个二次方程共11个方程来解算9个未知量:r11 r12 r21 r22 r31 r32 t1 t2 t3
而直接使用solve命令解算11个方程 CPU%100很长时间都没算出来
如果只用前8个一次线性方程解算 出现了 带sym的解
问题解答:
我来补答展开全文阅读