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

问题描述:

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;
这个结果明显不能符合原来的方程组.到底哪里错了呢?
1个回答 分类:综合 2014-10-22

问题解答:

我来补答
%% 这样就好了,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
 
 
展开全文阅读
剩余:2000
下一页:老师第一十三题