matlab solve 解这个方程组得出的结果不对G5*G10DivideByG5*X10+G17*X17=G18*X18(G18-G5)*X8+G6*X6=G13*X13G6=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
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/06 07:06:46
matlab solve 解这个方程组得出的结果不对G5*G10DivideByG5*X10+G17*X17=G18*X18(G18-G5)*X8+G6*X6=G13*X13G6=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
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;
这个结果明显不能符合原来的方程组.到底哪里错了呢?
matlab solve 解这个方程组得出的结果不对G5*G10DivideByG5*X10+G17*X17=G18*X18(G18-G5)*X8+G6*X6=G13*X13G6=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
%% 这样就好了,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