本帖產(chǎn)生 1 個計(jì)算強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

v1iv1ia1n

應(yīng)助: 0 (幼兒園)
金幣: 30721.4
散金: 18
紅花: 14
帖子: 5617
在線: 557.7小時(shí)
蟲號: 1193844
注冊: 2011-01-21
專業(yè): 化學(xué)環(huán)境污染與健康

[求助] matlab參數(shù)擬合

求助：
   不太會用matlab，請路過的大俠拔刀相助！！

   想求ku、kep、kem和km這幾個參數(shù)，微分方程如下：
   dCtot/dt = (ku*1.5)-(kep*Cp)-(kem*Cm)
   dCp/dt = (ku*1.5) -(km*Cp)-(kep*Cp)
   dCm/dt = (km*Cp)-(kem*Cm)

   已知數(shù)據(jù)為Ctot和t：
t       Ctot
0.5h
1       0.326
2       0.289
3       0.443
AVE       0.352666667
1h
1       0.362
2       0.38
3       0.417
AVE       0.386333333
2h
1       0.43
2       0.401
3
AVE       0.4155
4h
1       0.292
2       0.368
3       0.538
AVE       0.399333333
8h
1       0.507
2       0.492
3       0.513
AVE       0.504
12h
1       0.275
2       0.183
3       0.202
AVE       0.22
24h
1       0.111
2       0.122
3       0.094
AVE       0.109

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : data.xls

2012-04-13 16:06:47, 39.5 K

» 猜你喜歡

面向全球申請人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有172人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有2人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有9人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

matlab擬合方程參數(shù)時(shí)初值的選擇已經(jīng)有15人回復(fù)
matlab擬合求值已經(jīng)有17人回復(fù)
在matlab中擬合曲面的表達(dá)式已經(jīng)有29人回復(fù)
matlab擬合擬均相動力學(xué)參數(shù) 已經(jīng)有6人回復(fù)
需要用matlab擬合數(shù)據(jù)的外行人急需幫助已經(jīng)有10人回復(fù)
matlab非線性參數(shù)擬合問題已經(jīng)有7人回復(fù)
【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合已經(jīng)有17人回復(fù)
【求助】使用MATLAB擬合出按照公式的系數(shù) 已經(jīng)有19人回復(fù)
如何用MATLAB擬合指數(shù)函數(shù) 已經(jīng)有3人回復(fù)
【求助】matlab曲面擬合表達(dá)式已經(jīng)有10人回復(fù)
【求助】如何使得matlab擬合的曲線強(qiáng)制經(jīng)過一個點(diǎn) 已經(jīng)有8人回復(fù)
【求助】使用MATLAB怎么實(shí)現(xiàn)擬合力場參數(shù)的程序？【已解決】已經(jīng)有15人回復(fù)

1樓 2012-04-13 16:07:39

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

v1iv1ia1n

至尊木蟲 (知名作家)

應(yīng)助: 0 (幼兒園)
金幣: 30721.4
散金: 18
紅花: 14
帖子: 5617
在線: 557.7小時(shí)
蟲號: 1193844
注冊: 2011-01-21
專業(yè): 化學(xué)環(huán)境污染與健康

能提供寫思路也行

贊一下

回復(fù)此樓

2樓2012-04-13 21:37:20

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

逃之夭夭啦

銅蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 85
帖子: 36
在線: 20.8小時(shí)
蟲號: 1411106
注冊: 2011-09-21
專業(yè): 膠體與界面化學(xué)

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

感謝參與，應(yīng)助指數(shù) +1

不知道你是否學(xué)習(xí)過線性代數(shù)或者數(shù)值分析，你的就是三元一次方程組，求逆矩陣然后求解就可以了。這是matlab最擅長的，介意看一下matlab的寶典。

贊一下

回復(fù)此樓

3樓2012-04-14 09:59:06

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

逃之夭夭啦

銅蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 85
帖子: 36
在線: 20.8小時(shí)
蟲號: 1411106
注冊: 2011-09-21
專業(yè): 膠體與界面化學(xué)

不好意思，沒看清楚，你這個三方程，只有一組已知數(shù)dctot，而且還要求四個未知數(shù)的值，是不是有其他條件呢？

贊一下

回復(fù)此樓

4樓2012-04-14 10:05:00

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

逃之夭夭啦

銅蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 85
帖子: 36
在線: 20.8小時(shí)
蟲號: 1411106
注冊: 2011-09-21
專業(yè): 膠體與界面化學(xué)

★
dbb627: 金幣+1, 歡迎交流 2012-04-14 10:59:37

不好意思，沒看清楚，你還有微分方程問題，你可以試試分散點(diǎn)求斜率的方法試試，這樣應(yīng)該能得到多個已知量，然后求解。但是不知道你的各個參數(shù)的物理意義，不知道可不可以那么做。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

v1iv1ia1n

5樓2012-04-14 10:15:59

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

dingd

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

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
cenwanglai: 金幣+5, 計(jì)算強(qiáng)帖+1, 非常感謝熱心應(yīng)助，歡迎常來交流探討~ 2012-04-14 21:17:21
v1iv1ia1n: 金幣+30, ★★★★★最佳答案, 謝謝大俠，你太牛了，我好好研究研究�。。� 2012-04-14 23:32:28

應(yīng)該屬于微分方程擬合了，Cp、Cm的初值都知道嗎？下面是1stOpt求解代碼：

CODE:

Variable t,Ctot;

ODEFunction Ctot' = (ku*1.5)-(kep*Cp)-(kem*Cm);

            Cp' = (ku*1.5) -(km*Cp)-(kep*Cp);

            Cm' = (km*Cp)-(kem*Cm);

Data;//0.5h

1        0.326

2        0.289

3        0.443

Data;// 1h

1        0.362

2        0.38

3        0.417

Data;// 2h

1        0.43

2        0.401

3

Data;//4h

1        0.292

2        0.368

3        0.538

Data;//8h

1        0.507

2        0.492

3        0.513

Data;//12h

1        0.275

2        0.183

3        0.202

Data;//24h

1        0.111

2        0.122

3        0.094

結(jié)果（有多組）：
均方差(RMSE): 7.21477021187873E-12
殘差平方和(SSE): 7.28740728942977E-22
相關(guān)系數(shù)(R): 1
相關(guān)系數(shù)之平方(R^2): 1
決定系數(shù)(DC): 1

參數(shù)       最佳估算
-------------------- -------------
ku       0.0276775305197549
kep       -0.0456535422528521
kem       0.00133632074107808
km       0.514129491172468
cp初值(數(shù)據(jù)文件 - 1) -13.5280340784502
cm初值(數(shù)據(jù)文件 - 1) -306.857892086343
cp初值(數(shù)據(jù)文件 - 2) -1.26798177439112
cm初值(數(shù)據(jù)文件 - 2) -16.1124930625247
cp初值(數(shù)據(jù)文件 - 3) -3.11756100776893
cm初值(數(shù)據(jù)文件 - 3) -30.9936055993614
cp初值(數(shù)據(jù)文件 - 4) -6.62174164284473
cm初值(數(shù)據(jù)文件 - 4) -204.492835490761
cp初值(數(shù)據(jù)文件 - 5) -2.4771057659164
cm初值(數(shù)據(jù)文件 - 5) -24.1373425222873
cp初值(數(shù)據(jù)文件 - 6) -7.81801022205547
cm初值(數(shù)據(jù)文件 - 6) -111.07846956739
cp初值(數(shù)據(jù)文件 - 7) 2.86865604413812
cm初值(數(shù)據(jù)文件 - 7) 101.125291957442

====== 輸出結(jié)果 =====

文件：數(shù)據(jù)文件 - 1
No t 目標(biāo)ctot 計(jì)算ctot
1 2 0.289 0.288999999981285
2 3 0.443 0.443000000008992

文件：數(shù)據(jù)文件 - 2
No t 目標(biāo)ctot 計(jì)算ctot
1 2 0.38 0.37999999999826
2 3 0.417 0.417000000000818

文件：數(shù)據(jù)文件 - 3
No t 目標(biāo)ctot 計(jì)算ctot
1 2 0.401 0.400999999995707
2 3 0.417 0.417000000002048

文件：數(shù)據(jù)文件 - 4
No t 目標(biāo)ctot 計(jì)算ctot
1 2 0.368 0.367999999990831
2 3 0.538 0.538000000004393

文件：數(shù)據(jù)文件 - 5
No t 目標(biāo)ctot 計(jì)算ctot
1 2 0.492 0.491999999996592
2 3 0.513 0.513000000001621

文件：數(shù)據(jù)文件 - 6
No t 目標(biāo)ctot 計(jì)算ctot
1 2 0.183 0.182999999989207
2 3 0.202 0.202000000005179

文件：數(shù)據(jù)文件 - 7
No t 目標(biāo)ctot 計(jì)算ctot
1 2 0.122 0.12200000000272
2 3 0.094 0.0939999999982719

贊一下

回復(fù)此樓

6樓2012-04-14 20:29:52

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

v1iv1ia1n

至尊木蟲 (知名作家)

應(yīng)助: 0 (幼兒園)
金幣: 30721.4
散金: 18
紅花: 14
帖子: 5617
在線: 557.7小時(shí)
蟲號: 1193844
注冊: 2011-01-21
專業(yè): 化學(xué)環(huán)境污染與健康

送鮮花一朵

引用回帖:

5樓: Originally posted by 逃之夭夭啦 at 2012-04-14 10:15:59:
不好意思，沒看清楚，你還有微分方程問題，你可以試試分散點(diǎn)求斜率的方法試試，這樣應(yīng)該能得到多個已知量，然后求解。但是不知道你的各個參數(shù)的物理意義，不知道可不可以那么做。

非常謝謝~！~~

回復(fù)此樓

7樓2012-04-14 23:33:44

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

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

[考研] 08工科求調(diào)劑286 +4	tgs_001 2026-03-28	4/200	2026-03-30 21:03 by dophin1985
[考研] 學(xué)碩274求調(diào)劑 +10	Li李魚 2026-03-26	10/500	2026-03-30 20:54 by dophin1985
[考研] 材料科學(xué)與工程求調(diào)劑 +7	深V宿舍吧 2026-03-29	7/350	2026-03-30 20:43 by dophin1985
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +5	12.27 2026-03-27	5/250	2026-03-30 19:31 by 源_2020
[考研] 309求調(diào)劑 +15	誰不是少年 2026-03-29	15/750	2026-03-30 19:28 by 2026材料調(diào)劑
[考研] 311求調(diào)劑 +7	藍(lán)月亮亮 2026-03-30	7/350	2026-03-30 19:24 by 無際的草原
[考研] 349求調(diào)劑 +3	田洪有 2026-03-30	3/150	2026-03-30 19:11 by 迷糊CCPs
[考研] 287求調(diào)劑 +14	land xuxu 2026-03-26	14/700	2026-03-30 18:38 by 544594351
[考研] 材料科學(xué)與工程調(diào)劑 +8	深V宿舍吧 2026-03-30	8/400	2026-03-30 16:19 by 晶體之美
[考研] 求調(diào)劑 +7	青春裁為三截 2026-03-29	7/350	2026-03-30 13:14 by laoshidan
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5	hanamiko 2026-03-29	5/250	2026-03-29 16:46 by 學(xué)員8dgXkO
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8	mikasa的圍巾 2026-03-28	8/400	2026-03-29 12:48 by 無際的草原
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 085602 化工專碩 338分求調(diào)劑 +12	路癡小琪 2026-03-27	12/600	2026-03-28 15:41 by L135790
[考研] 材料求調(diào)劑一志愿哈工大324 +7	閆旭東 2026-03-28	9/450	2026-03-28 08:51 by Xu de nuo
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。�！ +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[考研] 279 分求調(diào)劑 +4	睡個好覺_16 2026-03-24	4/200	2026-03-27 15:05 by 醉在風(fēng)里
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3	zzffylgg 2026-03-24	3/150	2026-03-27 10:42 by shangxh
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518