Question

MATLAB求解方程组出错

当我用四个方程组，四个未知数来求解时，可以顺利得到结果。命令如：clc;clear all;close all;syms X Y W US=solve('X=Y*exp((76789-43.8*723.15-(169105-120.4*723.15)*X)/8.310/723.15)','W*(1-X)=exp(-0.011*(120-25))*(1-0.0182)','U*Y+W*X=0.0182','U+W=1');S.XS.YS.WS.U可是，当我想对其中一个值，比如120-25中的25用一个未知数T代替。并且给T赋值（T=0:5:120）就无法得到对应的XYWU，而我手动把这个25在0到120范围内改时，却可以求出每个XYWU。错误命令如下：clc;clear all;close all;syms X Y W UT=0:10:120S=solve('X=Y*exp((76789-43.8*723.15-(169105-120.4*723.15)*X)/8.310/723.15)','W*(1-X)=exp(-0.011*(120-T))*(1-0.0182)','U*Y+W*X=0.0182','U+W=1');S.XS.YS.WS.U

Answer 1

clc;clear all;close all;
syms X Y W U
% S=solve('X=Y*exp((76789-43.8*723.15-(169105-120.4*723.15)*X)/8.310/723.15)',...
% 'W*(1-X)=exp(-0.011*(120-25))*(1-0.0182)',...
% 'U*Y+W*X=0.0182','U+W=1');
T = 0:10:120;
result = zeros(length(T),4);
for i = 1:length(T)
Temp = T(i);
f1 = X-Y*exp((76789-43.8*723.15-(169105-120.4*723.15)*X)/8.310/723.15);
f2 = W*(1-X) - exp(-0.011*(120-Temp))*(1-0.0182);
f3 = U*Y+W*X - 0.0182;
f4 = U+W - 1;
S=solve(f1, f2, f3, f4);

result(i,1) = eval(S.X);
result(i,2) = eval(S.Y);
result(i,3) = eval(S.W);
result(i,4) = eval(S.U);
end
按上面代码求解，result存放你要的结果；
说明几点：solve是只能解单个方程组，方程组系数涉及向量或者矩阵无能为力；
原本可以将T看成一个符号，可以用solve求解其符号解，这样你可以把你的T(向量或者矩阵带入，用上述的eval函数求值)，但是这里解不出来符号解，可能是solve本身的问题；
所以这里没有办法，只能将T挨个值带入，利用循环求解。