Question

matlab求指数方程组

求教大神，怎么用MATLAB 实现下边这种方程组的求解？

Answer 1

He是常数吗如果是的话我的方法如下（因为没有具体数据，我就随便给了一下K、T和He的参数，你可以根据实际情况修改）。
主程序
clc;clear
K = [1;2;3;4];
T = [1;4;2;3];
He = 1;
Ct0 = rand(numel(K)+1,1);
options = optimset('TolFun',1e-16); % 精度要求
Ct = fsolve(@(Ct)eqs(Ct,K,T,He),Ct0,options);
C = Ct(1:end-1,1)
tmin = Ct(end,1)
函数程序

function F = eqs(Ct,K,T,He)
Ct = reshape(Ct,numel(Ct),1);
K = reshape(K,numel(K),1);
T = reshape(T,numel(T),1);
n = numel(K);
F = zeros(n+1,1);
F(1:n,1) = K.*(1-T/Ct(n+1,1).*(1-exp(-Ct(n+1,1)./T))) - Ct(1:n,1);
F(n+1,1) = 2*He/sum(Ct(1:n,1)) - Ct(n+1,1);
end
计算结果
C =
0.3803
0.2405
0.6647
0.6237
tmin =
1.0476