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

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

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

v1iv1ia1n

至尊木蟲(chóng) (知名作家)

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

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

求助：
   不太會(huì)用matlab，請(qǐng)路過(guò)的大俠拔刀相助��！

   想求ku、kep、kem和km這幾個(gè)參數(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)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : data.xls

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

» 猜你喜歡

中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有114人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(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é)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有2人回復(fù)

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

matlab擬合方程參數(shù)時(shí)初值的選擇已經(jīng)有15人回復(fù)
matlab擬合求值已經(jīng)有17人回復(fù)
在matlab中擬合曲面的表達(dá)式已經(jīng)有29人回復(fù)
matlab擬合擬均相動(dòng)力學(xué)參數(shù) 已經(jīng)有6人回復(fù)
需要用matlab擬合數(shù)據(jù)的外行人急需幫助已經(jīng)有10人回復(fù)
matlab非線性參數(shù)擬合問(wèn)題已經(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)過(guò)一個(gè)點(diǎn) 已經(jīng)有8人回復(fù)
【求助】使用MATLAB怎么實(shí)現(xiàn)擬合力場(chǎng)參數(shù)的程序？【已解決】已經(jīng)有15人回復(fù)

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

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

回帖置頂 ( 共有1個(gè) )

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
cenwanglai: 金幣+5, 計(jì)算強(qiáng)帖+1, 非常感謝熱心應(yīng)助，歡迎常來(lái)交流探討~ 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

至尊木蟲(chóng) (知名作家)

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

能提供寫(xiě)思路也行

贊一下

回復(fù)此樓

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

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

逃之夭夭啦

銅蟲(chóng) (初入文壇)

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

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

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

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

贊一下

回復(fù)此樓

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

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

逃之夭夭啦

銅蟲(chóng) (初入文壇)

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

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

贊一下

回復(fù)此樓

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

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

逃之夭夭啦

銅蟲(chóng) (初入文壇)

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

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

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

贊一下

回復(fù)此樓

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

v1iv1ia1n

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

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

v1iv1ia1n

至尊木蟲(chóng) (知名作家)

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

送鮮花一朵

引用回帖:

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

非常謝謝~！~~

回復(fù)此樓

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

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主學(xué)員ifMwFv 的主題更新

返回列表

普通表情龍兔虎貓

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

[考研] 310求調(diào)劑 +10	爭(zhēng)取九點(diǎn)睡 2026-03-30	10/500	2026-03-30 16:45 by ztnimte
[有機(jī)交流] 考研調(diào)劑 +6	watb 2026-03-26	6/300	2026-03-30 16:32 by herarysara
[考研] 0703一志愿9，初試成績(jī)：338，四六級(jí)已過(guò)，有科研經(jīng)歷，求調(diào)劑！ +6	Zuhui0306 2026-03-25	6/300	2026-03-30 16:31 by nothing投稿中
[考研] 11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè) +3	天賦帶到THU 2026-03-29	4/200	2026-03-30 16:24 by nihaoar
[考研] 322求調(diào)劑 +9	宋明欣 2026-03-27	9/450	2026-03-30 16:15 by wangjy2002
[考研] 298求調(diào)劑 +3	什么是胖頭魚(yú) 2026-03-30	5/250	2026-03-30 14:41 by 青海小西牛
[考研] 一志愿南航 335分 \| 0856材料化工 \| GPA 4.07 \| 有科研經(jīng)歷 +8	cccchenso 2026-03-23	8/400	2026-03-30 14:14 by w蟲(chóng)蟲(chóng)123
[基金申請(qǐng)] 面上5B能上會(huì)嗎？ +7	redcom 2026-03-29	7/350	2026-03-30 12:44 by 愛(ài)搞研究的小孩
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +7	nhj_ 2026-03-25	7/350	2026-03-30 12:44 by fangnagu
[考研] 生物技術(shù)與工程 +7	1294608413 2026-03-25	8/400	2026-03-30 11:36 by 唐沐兒
[考研] 297求調(diào)劑 +13	田洪有 2026-03-26	13/650	2026-03-30 11:33 by 替代品000
[考研] 求調(diào)劑 +10	張zz111 2026-03-27	11/550	2026-03-30 09:17 by 無(wú)際的草原
[考研] 考研調(diào)劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學(xué)員8dgXkO
[考研] 329求調(diào)劑 +10	鈕恩雪 2026-03-25	10/500	2026-03-29 13:32 by peike
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +7	@taotao 2026-03-29	7/350	2026-03-29 12:03 by longlotian
[考研] 343求調(diào)劑 +6	愛(ài)羈絆 2026-03-29	6/300	2026-03-29 12:00 by 無(wú)際的草原
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 311（085601）求調(diào)劑 +4	liziyeyeye 2026-03-28	4/200	2026-03-28 18:50 by 535743368
[考研] 0856，材料與化工321分求調(diào)劑 +12	大饞小子 2026-03-27	13/650	2026-03-28 10:56 by self2008
[考研] 315調(diào)劑 +4	0860求調(diào)劑 2026-03-26	5/250	2026-03-27 11:23 by wangjy2002