版塊導(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
專(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的回帖

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
在線: 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的回帖

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
專(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,虛線變成實(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
專(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的回帖

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

返回列表

普通表情龍兔虎貓

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

[考研] 生物學(xué)調(diào)劑招人！��！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋?xiě)?/a>
[考研] 294求調(diào)劑材料與化工專(zhuān)碩 +13	陌の森林 2026-03-18	13/650	2026-03-19 19:41 by maocaozhuxi
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	5/250	2026-03-19 18:03 by zcl123
[考博] 東華理工大學(xué)化材專(zhuān)業(yè)26屆碩士博士申請(qǐng) +8	zlingli 2026-03-13	8/400	2026-03-19 16:32 by 輕松不少隨
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +3	一瓶苯甲酸 2026-03-14	3/150	2026-03-19 15:17 by 盡舜堯1
[考研] 266求調(diào)劑 +5	陽(yáng)陽(yáng)哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽(yáng)陽(yáng)哇塞
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽(yáng)w
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 一志愿天大材料與化工（085600）總分338 +5	蔡大美女 2026-03-13	5/250	2026-03-19 10:44 by 是小劉呀～
[考研] 0703化學(xué)調(diào)劑 +8	妮妮ninicgb 2026-03-15	12/600	2026-03-19 05:55 by anny19840123
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無(wú)際的草原
[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無(wú)際的草原
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo