也许你可以给更好的初始猜测x0解决方案比我可以…
x0 = 10 *兰德(7 - 1);
选择= optimset (“MaxFunEvals”,1000000,“麦克斯特”,1000000);
规范(乐趣(x0))
x = fsolve (x0, @fun选项)
规范(有趣的(x))
函数F (x) =乐趣
Pt = x (1);
Tt = x (2);
rhot = x (3);
Zt型= x (4);
Vt = x (5);
= x (6);
xt = x (7);
di = 0.04;% k
α= 12.67;%的学位
Ti = 291.65;% K
π= 90;%的酒吧
Vi = 43.5;% m / s
Ci = 443.8;% m / s
米= 0.098;
子= 1;
rhoi = 60.17;% ^ 3公斤/米
k = 1.32;
R = 8.3114;% kj / kmol.K
Ai = 0.0013;% m ^ 2
%方程:
F (1) = Pt / rhot - Tt Zt型* R *;
F (2) =√k * R * Zt型* Tt) - Vt;
F (3) = (k * R * Zt型* Tt)——(Vi ^ 2 + (2 * k) / (k - 1)) *(子* * Ti-Zt * R * Tt));
F(4) =((子* Ti) / (Zt型* Tt))——(rhoi / rhot) ^ (k - 1);
F (5) = (rhoi / rhot) ^ (k - 1) -(π/ Pt) ^ ((k - 1) / k);
F (6) = Ai * (rhoi / rhot) * (Vi / Vt)——在;
F (7) = (Ai在^ ^ 0.5 - 0.5)/((π^ 0.5)*罐内(α))- xt;
结束
