| 查看: 3719 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
一路向北1211新虫 (小有名气)
|
[求助]
1stOpt拟合动力学参数 已有2人参与
|
||
|
听说这个软件1stOpt可以拟合动力学参数，涉及到多参数的情况下比MATLAB更高效？ 实验数据ExpData = ... t y1 y2 y3 y4 y5 y6 [ 20276.4 0.524933 0.015160 0.034247 0.049678 0.042493 0.048705 15544.8 0.565928 0.016405 0.025675 0.035539 0.105374 0.036335 12435.8 0.635708 0.017098 0.018350 0.023165 0.127280 0.025302 10363.2 0.663704 0.015437 0.011760 0.016443 0.161281 0.019342 7772.40 0.710197 0.012874 0.006357 0.010724 0.171001 0.013644 15840.0 0.478691 0.015317 0.038807 0.056022 0.052035 0.049784 12672.0 0.544696 0.017613 0.031530 0.038851 0.095962 0.036315 10560.0 0.585210 0.018444 0.025214 0.028700 0.125622 0.029087 7920.00 0.629790 0.017852 0.015888 0.018250 0.166011 0.021326 6336.00 0.664716 0.016108 0.010087 0.013052 0.185704 0.015519 10758.0 0.510701 0.018743 0.039316 0.045492 0.078866 0.038200 8068.50 0.578239 0.020151 0.028095 0.028405 0.134897 0.024460 6454.80 0.622556 0.019838 0.019259 0.018934 0.166062 0.018166 5379.00 0.652604 0.018307 0.014439 0.014335 0.180215 0.014117 4610.60 0.667576 0.016048 0.009840 0.010565 0.200809 0.012332 10954.8 0.446400 0.017853 0.047938 0.056591 0.060466 0.045798 6572.90 0.561169 0.020556 0.033192 0.031908 0.132280 0.021964 4694.90 0.624082 0.020550 0.020227 0.017017 0.177808 0.014202 3651.60 0.659205 0.018635 0.012475 0.011012 0.200108 0.010419 2987.70 0.665006 0.015238 0.007869 0.010369 0.221306 0.009117 5575.80 0.514906 0.021258 0.042926 0.049267 0.068656 0.030418 4181.90 0.561172 0.021724 0.034264 0.037491 0.115704 0.019845 3345.50 0.641624 0.021579 0.021906 0.020604 0.146161 0.012826 2787.90 0.734110 0.021300 0.016030 0.015080 0.161670 0.009900 2389.60 0.850940 0.021100 0.011580 0.010600 0.172740 0.007470 ] 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; y0=[0.67,0,0,0,0.33,0]; %各种组分的初值 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); k0=[1.1*10^(-8),1.74*10^(-8),3.52*10^(-14),5*10^(-10),5*10^(-10),9.4*10^(-10),0.76,0.24,0.36,1,0.58,0.5,0.5,0.5,0.5,0.5,0.5,1.88*10^(-3),4.29*10^(-5),4.15*10^(-4),2.55*10^(-5),2.55*10^(-4),2.55*10^(-5),2.55*10^(-3)]; 还请帮忙，@dingd |
回复此楼» 收录本帖的淘帖专辑推荐
 常微分方程 |
» 猜你喜欢
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有71人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
-
3月 中国石油和化学工业联合会科学技术奖 马上开始，申获奖要点
已经有15人回复
-
中国石油和化学工业联合会科学技术奖的评审标准是什么？
已经有3人回复
-
经费如何使用，花不完怎么办？
已经有0人回复
一路向北1211
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 683.6
- 散金: 354
- 帖子: 161
- 在线: 36小时
- 虫号: 3433050
- 注册: 2014-09-22
- 性别: GG
- 专业: 化学反应工程
4楼2017-12-23 08:58:55
一路向北1211
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 683.6
- 散金: 354
- 帖子: 161
- 在线: 36小时
- 虫号: 3433050
- 注册: 2014-09-22
- 性别: GG
- 专业: 化学反应工程
送红花一朵 3楼2017-12-22 20:15:40
一路向北1211
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 683.6
- 散金: 354
- 帖子: 161
- 在线: 36小时
- 虫号: 3433050
- 注册: 2014-09-22
- 性别: GG
- 专业: 化学反应工程
6楼2017-12-24 22:21:06
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
7楼2017-12-25 13:34:17
