Question

matlab solve 解这个方程组得出的结果不对

G5*G10DivideByG5*X10+G17*X17=G18*X18

(G18-G5)*X8+G6*X6=G13*X13

G6=G5*(1-G10DivideByG5)

其中：G10DivideByG5=0.9271 ； X10=0.3030；G17=1； X17=0.6；X18=0.35；
X8=0.35； X6=0.9485； G13=1；X13=0.6；

用solve执行如下命令：
[G5,G18,G6]=solve('G5*G10DivideByG5*X10+G17*X17=G18*X18','(G18-G5)*X8+G6*X6=G13*X13','G6=G5*(1-G10DivideByG5)','G5','G18','G6');

得到：
G18=(G13*X13*X18-G17*X17*X8)/(X18*X6-X18*X8-G10DivideByG5*X18*X6+G10DivideByG5*X10*X8);
G5=(G17*X17*X6-G17*X17*X8+G13*G10DivideByG5*X10*X13-G10DivideByG5*G17*X17*X6)/(X18*X6-X18*X8-G10DivideByG5*X18*X6+G10DivideByG5*X10*X8);
G6=(G13*X13*X18-G17*X17*X8-G13*G10DivideByG5*X13*X18+G10DivideByG5*G17*X17*X8)/(X18*X6-X18*X8-G10DivideByG5*X18*X6+G10DivideByG5*X10*X8);

计算结果得到：
G5=1.7143
G18=0；
G6=0；
这个结果明显不能符合原来的方程组. 到底哪里错了呢？

Answer 1

%% 这样就好了，matlab不能把解一一对应

clear;clc
G10DivideByG5=0.9271;
X10=0.3030;G17=1;
X17=0.6;X18=0.35;
X8=0.35;X6=0.9485;G13=1;
X13=0.6;
syms G5 G6 G18
eq1=G5*G10DivideByG5*X10+G17*X17-G18*X18;
eq2=(G18-G5)*X8+G6*X6-G13*X13;
eq3=G6-G5*(1-G10DivideByG5);
D=solve(eq1,eq2,eq3);
g5=D.G5
g6=D.G6
g18=D.G18