Question

用matlab解方程组，结果不正确！向大家求教，多谢！

Answer 1

clc;clear
syms f d;
F1 = d+1-d/f-0.6;
F2 = 0.34*f/(0.34+f)+0.6-2*d;
[d0,f0] = solve(F1,F2,'d','f');
fprintf('解为\nd1=%s\nf1=%s\n或\nd2=%s\nf2=%s\n',(d0(1,1)),f0(1,1),d0(2,1),f0(2,1))
% 校验误差,可以知道误差计算为0
err1 = simplify((subs(F1,[d,f],([d0(1,1),f0(1,1)]))));
err2 = simplify((subs(F2,[d,f],([d0(1,1),f0(1,1)]))));
err3 = simplify((subs(F1,[d,f],([d0(2,1),f0(2,1)]))));
err4 = simplify((subs(F2,[d,f],([d0(2,1),f0(2,1)]))));
fprintf('解的计算误差分别为\n%s\n%s\n%s\n%s\n',err1,err2,err3,err4)

结果为

解为
d1=109/670 - 25549^(1/2)/670
f1=2/15 - 25549^(1/2)/435
或
d2=25549^(1/2)/670 + 109/670
f2=25549^(1/2)/435 + 2/15
解的计算误差分别为
0
0
0
0

追问

你好，你上面的结果是对的，我验证过了。可是我的程序错在哪里?（d、f数据值交换了）

追答

你应该用[d,f] = solve();
因为matlab对未知数的顺序有一套自己的规则，一般来说x优先，然后是y、z然后字母按照排序来算