| 查看: 464 | 回复: 0 | ||
[求助]
求助，方程组求解出数值解后，第一个解是下个循环的一个参数
|
|
求解方程组循环的问题，程序是这样的 clear; x0=[300 300 300 300 300 0.5]; x1 = fsolve(@myfun,x0) function F = myfun(x); Tf=x(1); Tb=x(2); Tri=x(3); Tro=x(4); Tco=x(5); es=x(6); eq1=3.14*Df*l*hbf*(Tb-Tf)-(kb*Df*l/ab*(Tri-Tb)); eq2=2*3.14*l*kr*(Tro-Tri)/log(Dro/Dri)-(kb*l*Df/ab*(Tri-Tb)); eq3=G*0.34*l*abs*ref*tao-2*3.14*l*kr*(Tro-Tri)/log(Dro/Dri)-3.14*Dro*hro*l*(Tro^2+Tco^2)*(Tro+Tco)*(Tro-Tco); eq4=(3.14*b*Dco*aco*l*(Tco^4-Tsky^4)+3.14*Dco*l*hca*(Tco-Ta))-G*0.34*l*0.05*ref-3.14*Dro*hro*l*(Tro^2+Tco^2)*(Tro+Tco)*(Tro-Tco); eq5=2*m*Cf*(Tf-Tfi)+ql-G*0.34*l*abs*ref*tao+G*0.34*l*0.05*ref; eq6=es-2*m*Cf*(Tf-Tfi)/(G*0.34*l); F=[eq1 eq2 eq3 eq4 eq5 eq6] 有五个方程，给定一个Tfi 初始值，求解五个出数值解后的Tf下一个循环的Tfi，不知这一步如何解决，求教大神，不胜感激！ |
回复此楼» 猜你喜欢
十年后又回来了，论文投稿求助
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
fileCode有新解读？
已经有17人回复
-
这样的filecode谁见过
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
