matlab 解方程问题

a =-5.0000 -2.0000 1.0000 4.0000
1.5000 -4.8000 0.2000 1.3000
1.0000 0.6000 -2.5000 -0.5000
0.9000 -3.0000 -1.5000 -7.0000
syms x y z ；
equation1=abs(x*a(1,3))+abs(a(2,3))+a(3,3)+abs(a(4,3))-(a(1,1)+abs(1/x*a(2,1))+abs(1/x*a(3,1))+abs(1/x*a(4,1)))；

equation2=abs(a(1,2)*y)+a(2,2)+abs(a(3,2))+abs(a(4,2))- (a(1,1)+abs(1/y*a(2,1))+abs(1/y*a(3,1))+abs(1/y*a(4,1)))；

equation3=abs(z*a(1,4))+abs(a(2,4))+abs(a(3,4))+a(4,4)-(a(1,1)+abs(1/z*a(2,1))+abs(1/z*a(3,1))+abs(1/z*a(4,1)))；
问题：用
x=solve(eval(equation1))；
y=solve(eval(equation2))；
z=solve(eval(equation3))；
和用
[x y z]=solve(eval(equation1),eval(equation2),eval(equation3))
求解得到的答案不同。 后一个解法中含有新引入的符号变量。 不知为何？求教各位高手，非诚勿扰~~~谢谢

clear
clc
close all
a =[-5.0000 -2.0000 1.0000 4.0000
1.5000 -4.8000 0.2000 1.3000
1.0000 0.6000 -2.5000 -0.5000
0.9000 -3.0000 -1.5000 -7.0000];
syms x y z
equation1=abs(x*a(1,3))+abs(a(2,3))+a(3,3)+abs(a(4,3))-(a(1,1)+abs(1/x*a(2,1))+abs(1/x*a(3,1))+abs(1/x*a(4,1)));
equation2=abs(a(1,2)*y)+a(2,2)+abs(a(3,2))+abs(a(4,2))- (a(1,1)+abs(1/y*a(2,1))+abs(1/y*a(3,1))+abs(1/y*a(4,1)));
equation3=abs(z*a(1,4))+abs(a(2,4))+abs(a(3,4))+a(4,4)-(a(1,1)+abs(1/z*a(2,1))+abs(1/z*a(3,1))+abs(1/z*a(4,1)));
x=solve(eval(equation1));
y=solve(eval(equation2));
z=solve(eval(equation3));
double(x)
double(y)
double(z)
[x0 y0 z0]=solve(equation1,equation2,equation3);
double(x0)
double(y0)
double(z0)
你可以看看结果是不是一样的

举个例子，希望有所帮助。 代码 % by lyqmath clc; clear all; close all; syms a [x,y]=solve('a*x*x-b*y=10,c*y-4*x=3') subs(x, a, 4) subs(y, a, 4) 结果     x =    (8*b + 3*a*c + 4*(4*b^2 + 3*a*b*c + 10*a*c^2)^(1/2))/(4*a*c) - 3/4  (8*b + 3*a*c - 4*(4*b^2 + 3*a*b*c + 10*a*c^2)^(1/2))/(4*a*c) - 3/4     y =    (8*b + 3*a*c + 4*(4*b^2 + 3*a*b*c + 10*a*c^2)^(1/2))/(a*c^2)  (8*b + 3*a*c - 4*(4*b^2 + 3*a*b*c + 10*a*c^2)^(1/2))/(a*c^2)     ans =    (8*b + 12*c + 4*(4*b^2 + 12*b*c + 40*c^2)^(1/2))/(16*c) - 3/4  (8*b + 12*c - 4*(4*b^2 + 12*b*c + 40*c^2)^(1/2))/(16*c) - 3/4     ans =    (8*b + 12*c + 4*(4*b^2 + 12*b*c + 40*c^2)^(1/2))/(4*c^2)  (8*b + 12*c - 4*(4*b^2 + 12*b*c + 40*c^2)^(1/2))/(4*c^2)   >>

如以上回答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！

点此我要举报以上问答信息