| 查看: 1332 | 回复: 1 | ||
[求助]
用遗传算法和模拟退火来进行动力学参数的拟合
|
|
用matlab中的lsqnonlin函数、ode45函数进行动力学拟合，拟合指前因子、活化能、级数，一直未果。老师说让自己用遗传算法和模拟退火组合算，自己一点都没接触过。 我的问题是：原料是：甲醇、CO、氧气。主反应：甲醇、CO、氧气合成DMC，副反应：CO和氧气合成CO2 任务：拟合出生成DMC和CO2的动力学参数（指前因子、活化能、级数） 动力学方程为dy1/dw=f(y1、y2)和dy2/dw=f(y1、y2),具体形式见附图。 整理后的实验数据：N0(原料初始的摩尔数)、P（压力）、t(温度)、x10(甲醇初始的摩尔分数)、x20(CO初始的摩尔分数)、x30(O2初始的摩尔分数)、y1(DMC产物的摩尔分数)、y2(O2产物的摩尔分数)、W=0.375(催化剂的质量)。(每次的实验条件不同，原料经过固定床上催化剂的质量0.375g，反应后得到y1和y2) 目标函数：（y1拟合-y1）的平方和+（y2拟合-y2）的平方和 最小 data=... [ 0.3091 0.1000 413.1500 0.6240 0.2500 0.1270 0.0150 0.0180 0.3260 0.1000 413.1500 0.5000 0.3340 0.1660 0.0160 0.0200 0.2796 0.1000 413.3500 0.5310 0.3490 0.1210 0.0150 0.0180 0.2856 0.1000 416.5500 0.5710 0.2870 0.1420 0.0170 0.0330 0.2801 0.1000 423.1500 0.5300 0.3480 0.1220 0.0200 0.0400 0.3260 0.1000 423.2500 0.5000 0.3340 0.1660 0.0210 0.0410 0.3091 0.1000 423.3500 0.6240 0.2500 0.1270 0.0210 0.0420 0.2853 0.1000 424.3500 0.5720 0.2860 0.1420 0.0210 0.0420 0.3190 0.1000 426.8500 0.5110 0.3450 0.1440 0.0210 0.0450 0.2853 0.1000 432.3500 0.5720 0.2860 0.1420 0.0250 0.0540 0.3962 0.1000 433.1500 0.5240 0.3520 0.1240 0.0230 0.0530 0.3190 0.1000 434.4500 0.5110 0.3450 0.1440 0.0250 0.0540 0.3085 0.2000 413.3500 0.6250 0.2500 0.1250 0.0180 0.0310 0.2793 0.2000 414.2500 0.5310 0.3490 0.1200 0.0180 0.0320 0.2804 0.2000 422.8500 0.5290 0.3480 0.1230 0.0250 0.0500 0.2856 0.2000 423.4500 0.5710 0.2860 0.1430 0.0260 0.0510 0.3083 0.2000 423.9500 0.6250 0.2500 0.1240 0.0260 0.0500 0.3085 0.2000 432.8500 0.6250 0.2500 0.1250 0.0300 0.0590 0.2788 0.2000 434.3500 0.5320 0.3500 0.1180 0.0300 0.0580 0.2368 0.3000 411.1500 0.5010 0.3340 0.1650 0.0230 0.0430 0.2600 0.3000 423.4500 0.5700 0.2870 0.1420 0.0280 0.0510 0.2597 0.3000 433.1500 0.5710 0.2880 0.1410 0.0360 0.0630 0.3207 0.3000 433.3500 0.5550 0.3340 0.1110 0.0340 0.0610 0.2344 0.3000 435.4500 0.5060 0.3370 0.1570 0.0380 0.0660 ]; % N0 P T x10（1初始量）x20 x30 y1 y2 我下一步想用遗传算法和模拟退火，该怎么做呢？？请大神给些指导  图片1.jpg  图2.jpg |
回复此楼» 猜你喜欢
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然评审意见
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有13人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
-
系统查不到
已经有11人回复
2楼2016-06-07 23:30:32
