| 查看: 612 | 回复: 0 | ||
一路向北1211新虫 (小有名气)
|
[求助]
1stOpt高版本帮忙求解动力学参数，谢谢了
|
|
r1=k1*(y(1)*101325).*(K1*y(5)*101325).^(n1)./((1+(K1*y(1)*101325).^(n1)+K2*y(5)*101325).^2); r2=k2*(y(2)*101325).^(m2)*(y(5)*101325).^(n2)./((1+K3*y(6)*101325).^(n2)); r3=k3*(y(1)*101325).^(m3)*(y(5)*101325).^(n3)./((1+K4*y(6)*101325).^(n3)); r4=k4*(y(2)*101325).^(m4)*(y(5)*101325).^(n4)./((1+K5*y(6)*101325).^(n4)); r5=k5*(y(3)*101325).^(m5)*(y(5)*101325).^(n5)./((1+K6*y(6)*101325).^(n5)); r6=k6*(y(4)*101325).^(m6)*(y(5)*101325).^(n6)./((1+K7*y(6)*101325).^(n6)); dydt=zeros(6,1); dydt(1)=-2*r1-r3; dydt(2)=r1-r2-r4; dydt(3)=r2-r5; dydt(4)=r3+2*r4+2*r5-r6; dydt(5)=-0.5*(r1+r2+r6)-1.5*r3-2.5*r4-2*r5; dydt(6)=r6; k1=k(1); k2=k(2); k3=k(3); k4=k(4); k5=k(5); k6=k(6); n1=k(7); m2=k(8); n2=k(9); m3=k(10); n3=k(11); m4=k(12); n4=k(13); m5=k(14); n5=k(15); m6=k(16); n6=k(17); K1=k(18); K2=k(19); K3=k(20); K4=k(21); K5=k(22); K6=k(23); K7=k(24); ExpData = ... [ 0.524933 0.015160 0.034247 0.049678 0.042493 0.048705 0.565928 0.016405 0.025675 0.035539 0.105374 0.036335 0.635708 0.017098 0.018350 0.023165 0.127280 0.025302 0.663704 0.015437 0.011760 0.016443 0.161281 0.019342 0.710197 0.012874 0.006357 0.010724 0.171001 0.013644]; yexp=ExpData(:,1:6); t=[20726.4,15544.8,12435.8,10363.2,7772.4]； @dingd，能帮忙解决下吗？不胜感激 |
回复此楼» 猜你喜欢
求助中医院临床中医中标的青年标书学习，谢谢
已经有6人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
物理学I论文润色/翻译怎么收费?
已经有89人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有0人回复