Question

用matlab解三元二次方程组

如图所示，方程组里面有个三元二次方程，然后用如图中命令解这个方程组，结果没能得到数值解。图中的y的解是解析解吗？i是什么？我如何才能得到数值解。另外，这组方程肯定有解的。

Answer 1

图中的y的解是解析解，i是虚数。

如要得到数值解，可以这样做

>> fun=@(x)[0.132611776*(x(1)-8.104)+0.895881152*(x(2)-5.017)+0.146428341*(5-5.418),(x(1)-8.104)^2+(x(2)-5.017)^2+(5-5.418)^2];

>> [x,h]=fsolve(fun,rand(1,2))

执行得到如下结果

从运算结果来看，上述值并不是真正的数值解。

这说明该方程组无实数解，有复数解。

Answer 2

x1=-226.11; y1=238.802; z1=-1577; x2=223.1317; y2=201.5533; z2=-1577; x3=-277.2719; y3=-202.0605; z3=-1577; syms t1 t2 t3 eq1=(t1*x1-t2*x2)^2+(t1*y1-t2*y2)^2+(t1*z1-t2*z2)^2-532.8689^2; eq2=(t1*x1-t3*x3)^2+(t1*y1-t3*y3)^2+(t1*z1-t3*z3)^2-532.8689^2; eq3=(t3*x3-t2*x2)^2+(t3*y3-t2*y2)^2+(t3*z3-t2*z2)^2-753.5904^2; [t1 t2 t3]=solve(eq1,eq2,eq3,'t1','t2','t3') 楼主我用solve可以解出解来，你参考一下