版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

1064126551

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 916
帖子: 15
在線: 25.4小時
蟲號: 3099587
注冊: 2014-03-29
專業(yè): 分離過程

[求助] 用遺傳算法和模擬退火來進(jìn)行動力學(xué)參數(shù)的擬合

用matlab中的lsqnonlin函數(shù)、ode45函數(shù)進(jìn)行動力學(xué)擬合，擬合指前因子、活化能、級數(shù)，一直未果。老師說讓自己用遺傳算法和模擬退火組合算，自己一點都沒接觸過。
我的問題是：原料是：甲醇、CO、氧氣。主反應(yīng)：甲醇、CO、氧氣合成DMC，副反應(yīng)：CO和氧氣合成CO2
任務(wù)：擬合出生成DMC和CO2的動力學(xué)參數(shù)（指前因子、活化能、級數(shù)）
動力學(xué)方程為dy1/dw=f(y1、y2)和dy2/dw=f(y1、y2),具體形式見附圖。
整理后的實驗數(shù)據(jù)：N0(原料初始的摩爾數(shù))、P（壓力）、t(溫度)、x10(甲醇初始的摩爾分?jǐn)?shù))、x20(CO初始的摩爾分?jǐn)?shù))、x30(O2初始的摩爾分?jǐn)?shù))、y1(DMC產(chǎn)物的摩爾分?jǐn)?shù))、y2(O2產(chǎn)物的摩爾分?jǐn)?shù))、W=0.375(催化劑的質(zhì)量)。(每次的實驗條件不同，原料經(jīng)過固定床上催化劑的質(zhì)量0.375g，反應(yīng)后得到y(tǒng)1和y2)
目標(biāo)函數(shù)：（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
我下一步想用遺傳算法和模擬退火，該怎么做呢？？請大神給些指導(dǎo)

圖片1.jpg

圖2.jpg

回復(fù)此樓

» 猜你喜歡

焦慮已經(jīng)有10人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有3人回復(fù)
297求調(diào)劑已經(jīng)有3人回復(fù)
302求調(diào)劑已經(jīng)有8人回復(fù)
環(huán)境工程調(diào)劑已經(jīng)有7人回復(fù)
268求調(diào)劑已經(jīng)有6人回復(fù)
311求調(diào)劑已經(jīng)有11人回復(fù)
工科材料085601 279求調(diào)劑已經(jīng)有4人回復(fù)
301求調(diào)劑已經(jīng)有9人回復(fù)

1樓 2015-03-09 16:30:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

何fox

新蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 29.9
散金: 1879
紅花: 9
帖子: 1020
在線: 765小時
蟲號: 4033976
注冊: 2015-08-23
專業(yè): 高分子科學(xué)

現(xiàn)在這個算得咋樣了

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

2樓2016-06-07 23:30:32