版塊導(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)版塊排行榜

返回列表

cuphunter

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

應(yīng)助: 2 (幼兒園)
金幣: 1195.5
散金: 1440
紅花: 13
沙發(fā): 26
帖子: 6312
在線: 612.1小時(shí)
蟲(chóng)號(hào): 794813
注冊(cè): 2009-06-15
性別: GG
專業(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的回帖

dbb627

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

專家經(jīng)驗(yàn): +4
計(jì)算強(qiáng)帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲(chóng)號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(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的回帖

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 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的回帖

zl030091126

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

應(yīng)助: 15 (小學(xué)生)
金幣: 303.8
散金: 605
紅花: 2
帖子: 524
在線: 141.5小時(shí)
蟲(chóng)號(hào): 1444718
注冊(cè): 2011-10-16
性別: GG
專業(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,虛線變成實(shí)線則是你的軌跡

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
在線: 981.6小時(shí)
蟲(chóng)號(hào): 286920
注冊(cè): 2006-10-15
性別: GG
專業(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的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 cuphunter 的主題更新

返回列表

普通表情龍兔虎貓

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

[考研] 281求調(diào)劑（0805） +6	煙汐憶海 2026-03-16	16/800	2026-03-19 00:59 by 煙汐憶海
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開(kāi)心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 307求調(diào)劑 +5	超級(jí)伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629