matlab解三元方程组,请问如何用Matlab求解带参数的三元三次方程组

本帖最后由 zounidefeng 于 2012-10-24 11:28 编辑

zqyadam 发表于 2012-10-24 10:33

y1、y2、y3都是两排结果，这上下两排结果只差一个符号，请问是每个y1、y2、y3都对应着两个值吗？

计算结果如下(红色字体是我改的，区分两排结果)：

y1 =

(d*a*e-f*c*e+e*b*f-e*b*d-d*a*f-(-e*c^2*a*f^2+e*c^2*b*f^2-e^2*c^2*d*a+e^2*c^2*a*f-e^2*c^2*b*f+e^2*c^2*b*d+a*f^2*c^2*d+a^2*f^2*c*e-a^2*f^2*d*c-d*c^2*b*f^2-d^2*c^2*a*f+d^2*c^2*a*e-d^2*c^2*e*b+d^2*c^2*b*f-b^2*f^2*e*c+b^2*f^2*c*d+e^2*a^2*d*c-e^2*a^2*f*c+e*b^2*c*d^2-b^2*f*c*d^2+e^2*b^2*c*f-e^2*b^2*c*d+d^2*a^2*c*f-d^2*a^2*c*e+d*a^2*b*f^2-d^2*a^2*b*f+d^2*a^2*b*e-d^2*a*b^2*e-d*a*b^2*f^2+d^2*a*b^2*f-e*b*a^2*f^2+e^2*b*a^2*f-e^2*b*a^2*d-e^2*b^2*a*f+e^2*b^2*a*d+e*b^2*a*f^2)^(1/2)+d*c*f)/(e*a-e*c-a*f+d*c+b*f-b*d)

(d*a*e-f*c*e+e*b*f-e*b*d-d*a*f+(-e*c^2*a*f^2+e*c^2*b*f^2-e^2*c^2*d*a+e^2*c^2*a*f-e^2*c^2*b*f+e^2*c^2*b*d+a*f^2*c^2*d+a^2*f^2*c*e-a^2*f^2*d*c-d*c^2*b*f^2-d^2*c^2*a*f+d^2*c^2*a*e-d^2*c^2*e*b+d^2*c^2*b*f-b^2*f^2*e*c+b^2*f^2*c*d+e^2*a^2*d*c-e^2*a^2*f*c+e*b^2*c*d^2-b^2*f*c*d^2+e^2*b^2*c*f-e^2*b^2*c*d+d^2*a^2*c*f-d^2*a^2*c*e+d*a^2*b*f^2-d^2*a^2*b*f+d^2*a^2*b*e-d^2*a*b^2*e-d*a*b^2*f^2+d^2*a*b^2*f-e*b*a^2*f^2+e^2*b*a^2*f-e^2*b*a^2*d-e^2*b^2*a*f+e^2*b^2*a*d+e*b^2*a*f^2)^(1/2)+d*c*f)/(e*a-e*c-a*f+d*c+b*f-b*d)

y2 =

(-e*c^2*a*f^2+e*c^2*b*f^2-e^2*c^2*d*a+e^2*c^2*a*f-e^2*c^2*b*f+e^2*c^2*b*d+a*f^2*c^2*d+a^2*f^2*c*e-a^2*f^2*d*c-d*c^2*b*f^2-d^2*c^2*a*f+d^2*c^2*a*e-d^2*c^2*e*b+d^2*c^2*b*f-b^2*f^2*e*c+b^2*f^2*c*d+e^2*a^2*d*c-e^2*a^2*f*c+e*b^2*c*d^2-b^2*f*c*d^2+e^2*b^2*c*f-e^2*b^2*c*d+d^2*a^2*c*f-d^2*a^2*c*e+d*a^2*b*f^2-d^2*a^2*b*f+d^2*a^2*b*e-d^2*a*b^2*e-d*a*b^2*f^2+d^2*a*b^2*f-e*b*a^2*f^2+e^2*b*a^2*f-e^2*b*a^2*d-e^2*b^2*a*f+e^2*b^2*a*d+e*b^2*a*f^2)^(1/2)/(e*a-e*c-a*f+d*c+b*f-b*d)

-(-e*c^2*a*f^2+e*c^2*b*f^2-e^2*c^2*d*a+e^2*c^2*a*f-e^2*c^2*b*f+e^2*c^2*b*d+a*f^2*c^2*d+a^2*f^2*c*e-a^2*f^2*d*c-d*c^2*b*f^2-d^2*c^2*a*f+d^2*c^2*a*e-d^2*c^2*e*b+d^2*c^2*b*f-b^2*f^2*e*c+b^2*f^2*c*d+e^2*a^2*d*c-e^2*a^2*f*c+e*b^2*c*d^2-b^2*f*c*d^2+e^2*b^2*c*f-e^2*b^2*c*d+d^2*a^2*c*f-d^2*a^2*c*e+d*a^2*b*f^2-d^2*a^2*b*f+d^2*a^2*b*e-d^2*a*b^2*e-d*a*b^2*f^2+d^2*a*b^2*f-e*b*a^2*f^2+e^2*b*a^2*f-e^2*b*a^2*d-e^2*b^2*a*f+e^2*b^2*a*d+e*b^2*a*f^2)^(1/2)/(e*a-e*c-a*f+d*c+b*f-b*d)

y3 =

-(e*b*a-e*b*c+d*a*c-d*a*b-a*f*c+b*f*c+(-e*c^2*a*f^2+e*c^2*b*f^2-e^2*c^2*d*a+e^2*c^2*a*f-e^2*c^2*b*f+e^2*c^2*b*d+a*f^2*c^2*d+a^2*f^2*c*e-a^2*f^2*d*c-d*c^2*b*f^2-d^2*c^2*a*f+d^2*c^2*a*e-d^2*c^2*e*b+d^2*c^2*b*f-b^2*f^2*e*c+b^2*f^2*c*d+e^2*a^2*d*c-e^2*a^2*f*c+e*b^2*c*d^2-b^2*f*c*d^2+e^2*b^2*c*f-e^2*b^2*c*d+d^2*a^2*c*f-d^2*a^2*c*e+d*a^2*b*f^2-d^2*a^2*b*f+d^2*a^2*b*e-d^2*a*b^2*e-d*a*b^2*f^2+d^2*a*b^2*f-e*b*a^2*f^2+e^2*b*a^2*f-e^2*b*a^2*d-e^2*b^2*a*f+e^2*b^2*a*d+e*b^2*a*f^2)^(1/2))/(e*a-e*c-a*f+d*c+b*f-b*d)

-(e*b*a-e*b*c+d*a*c-d*a*b-a*f*c+b*f*c-(-e*c^2*a*f^2+e*c^2*b*f^2-e^2*c^2*d*a+e^2*c^2*a*f-e^2*c^2*b*f+e^2*c^2*b*d+a*f^2*c^2*d+a^2*f^2*c*e-a^2*f^2*d*c-d*c^2*b*f^2-d^2*c^2*a*f+d^2*c^2*a*e-d^2*c^2*e*b+d^2*c^2*b*f-b^2*f^2*e*c+b^2*f^2*c*d+e^2*a^2*d*c-e^2*a^2*f*c+e*b^2*c*d^2-b^2*f*c*d^2+e^2*b^2*c*f-e^2*b^2*c*d+d^2*a^2*c*f-d^2*a^2*c*e+d*a^2*b*f^2-d^2*a^2*b*f+d^2*a^2*b*e-d^2*a*b^2*e-d*a*b^2*f^2+d^2*a*b^2*f-e*b*a^2*f^2+e^2*b*a^2*f-e^2*b*a^2*d-e^2*b^2*a*f+e^2*b^2*a*d+e*b^2*a*f^2)^(1/2))/(e*a-e*c-a*f+d*c+b*f-b*d)