| 查看: 1331 | 回复: 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 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有18人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有12人回复
-
有没有仍没收到信息的
已经有3人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有8人回复
-
基金不中，共勉
已经有11人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有9人回复
-
2026年叶企孙基金
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有13人回复
-
基金系统什么内容也没有
已经有9人回复
-
怎么查啊
已经有6人回复
2楼2016-06-07 23:30:32
30