24小时热门版块排行榜    

CyRhmU.jpeg
查看: 396  |  回复: 0

ymgy

铜虫 (初入文坛)

[求助] 求助,方程组求解出数值解后,第一个解是下个循环的一个参数

求解方程组循环的问题,程序是这样的
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,不知这一步如何解决,求教大神,不胜感激!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ymgy 的主题更新
信息提示
请填处理意见