| 查看: 554 | 回复: 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，能帮忙解决下吗？不胜感激 |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有158人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有26人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有34人回复
