| 查看: 581 | 回复: 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，能帮忙解决下吗？不胜感激 |
回复此楼» 猜你喜欢
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有523人回复
-
物理学I论文润色/翻译怎么收费?
已经有84人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
