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

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

24小時(shí)熱門(mén)版塊排行榜

返回列表

當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

cuphunter

金蟲(chóng) (知名作家)

應(yīng)助: 2 (幼兒園)
金幣: 1195.5
散金: 1440
紅花: 13
沙發(fā): 26
帖子: 6312
在線(xiàn): 612.1小時(shí)
蟲(chóng)號(hào): 794813
注冊(cè): 2009-06-15
性別: GG
專(zhuān)業(yè): 化學(xué)反應(yīng)工程

[求助] matlab擬合動(dòng)力學(xué)參數(shù)

本人為matlab菜鳥(niǎo)，希望有高手相助，給出matlab代碼，本人萬(wàn)分感謝。這里給出的數(shù)據(jù)只是一個(gè)溫度下的，我想先學(xué)習(xí)一下，看自己能不能弄?jiǎng)e的溫度，再擬合活化能。這些等之后再交流。

數(shù)據(jù)：
初始值： L1(0) L2(0) L3(0) L4 (0) L5(0) L6(0)
X0 4.29 32.77 15.50 21.41 5.80 20.23

實(shí)驗(yàn)值：
編號(hào) L1 L2 L3 L4 L5 L6
1 6.95 36.54 13.13 17.56 3.86 21.96
2 6.47 34.98 11.95 17.09 4.07 25.45
3 9.45 40.82 11.87 11.80 3.03 23.03
4 8.62 50.21 10.01 7.49 1.71 21.96
5 10.33 46.08 9.11 9.63 2.55 22.3
6 9.13 40.62 11.38 12.01 3.29 23.56
7 7.26 37.41 13.64 16.65 3.46 21.57
8 7.61 37.20 13.22 15.21 3.46 23.3
9 8.07 37.86 12.28 14.67 3.58 23.54
10 7.32 36.66 14.39 14.58 3.72 23.54

動(dòng)力學(xué)模型為：

回復(fù)此樓

珍愛(ài)生命，遠(yuǎn)離木蟲(chóng).............................

1樓 2011-12-22 15:50:29

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

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線(xiàn): 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

1stOpt求解很方便，只是你的數(shù)據(jù)擬合結(jié)果不是太好�。�

CODE:

Variable x,a1,a2,a3,a4,a5,a6;

ODEFunction a1'=k31*a3;

            a2'=k42*a4;

            a3'=-(k31+k34+k35)*a3;

            a4'=-((k42+k45)*a4-k34*a3);

            a5'=-(k56*a5-k35*a3-k45*a4);

            a6'=k56*a5;

Data;

x        a1        a2        a3        a4        a5        a6

0        4.29        32.77        15.50        21.41        5.80        20.23

1        6.95        36.54        13.13        17.56        3.86        21.96

2        6.47        34.98        11.95        17.09        4.07        25.45

3        9.45        40.82        11.87        11.80        3.03        23.03

4        8.62        50.21        10.01        7.49        1.71        21.96

5        10.33        46.08        9.11        9.63        2.55        22.3

6        9.13        40.62        11.38        12.01        3.29        23.56

7        7.26        37.41        13.64        16.65        3.46        21.57

8        7.61        37.20        13.22        15.21        3.46        23.3

9        8.07        37.86        12.28        14.67        3.58        23.54

10        7.32        36.66        14.39        14.58        3.72        23.54

參數(shù) 最佳估算
-------------------- -------------
k31 0.0419690845909743
k42 0.0662882491041956
k34 0.413188420544208
k35 -0.42242795564006
k45 0.375861403372349
k56 0.0940392220979913

贊一下

回復(fù)此樓

3樓2011-12-22 18:27:57

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

查看全部 5 個(gè)回答

dbb627

榮譽(yù)版主 (著名寫(xiě)手)

專(zhuān)家經(jīng)驗(yàn): +4
計(jì)算強(qiáng)帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線(xiàn): 1794.8小時(shí)
蟲(chóng)號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專(zhuān)業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

X 和a分別是什么變量

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2樓2011-12-22 16:37:49

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

zl030091126

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 15 (小學(xué)生)
金幣: 303.8
散金: 605
紅花: 2
帖子: 524
在線(xiàn): 141.5小時(shí)
蟲(chóng)號(hào): 1444718
注冊(cè): 2011-10-16
性別: GG
專(zhuān)業(yè): 化學(xué)反應(yīng)工程

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1
cuphunter(金幣+2): 謝謝參與 2011-12-24 12:56:23

實(shí)驗(yàn)編號(hào)是間歇式的結(jié)果？或是一次反應(yīng)結(jié)束后的實(shí)驗(yàn)結(jié)果？a和X是什么變量，樓主都說(shuō)清楚了啊。a是組成？1 2 3...為時(shí)間？

贊一下

回復(fù)此樓

A--B,虛線(xiàn)變成實(shí)線(xiàn)則是你的軌跡

4樓2011-12-23 04:42:00

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

zt_chem

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 12 (小學(xué)生)
金幣: 2147.8
紅花: 6
帖子: 872
在線(xiàn): 981.6小時(shí)
蟲(chóng)號(hào): 286920
注冊(cè): 2006-10-15
性別: GG
專(zhuān)業(yè): 催化化學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1
cuphunter(金幣+2): 謝謝參與 2011-12-24 12:56:35

用lsqnonlin編程優(yōu)化吧

贊一下

回復(fù)此樓

5樓2011-12-23 11:55:31

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

查看全部 5 個(gè)回答

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 268求調(diào)劑 +9	簡(jiǎn)單點(diǎn)0 2026-03-17	9/450	2026-03-21 15:37 by lature00
[考研] 初試 317 +3	半拉月丙 2026-03-20	3/150	2026-03-21 13:25 by zhukairuo
[考研] 279分求調(diào)劑一志愿211 +14	chaojifeixia 2026-03-19	15/750	2026-03-21 13:24 by zhukairuo
[考研] 346求調(diào)劑[0856] +4	WayneLim327 2026-03-16	7/350	2026-03-21 04:02 by JourneyLucky
[考研] 310求調(diào)劑 +3	baibai1314 2026-03-16	3/150	2026-03-21 03:56 by JourneyLucky
[考研] 一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 317求調(diào)劑 +8	申子申申 2026-03-19	13/650	2026-03-21 00:09 by 劉國(guó)森
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4	葵梓衛(wèi)隊(duì) 2026-03-18	6/300	2026-03-20 23:02 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡(jiǎn)木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無(wú)懈可擊111
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 機(jī)械專(zhuān)碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky