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

登錄注冊

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

24小時熱門版塊排行榜

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

返回列表

yupeifeng

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 960.3
散金: 34
紅花: 1
帖子: 220
在線: 78.4小時
蟲號: 1421266
注冊: 2011-09-28
專業(yè): 無機材料化學(xué)

[求助] matlab Rung-kutta算法

由于做論文，做動力學(xué)這一塊，有動力學(xué)模型，有幾組濃度隨時間變化的曲線數(shù)據(jù)，來求動力學(xué)模型里的參數(shù)，進而得到動力學(xué)方程。
數(shù)據(jù)是：(c,t)=（0,0.75）（2,0.7）（3,0.65）（4,0.64）（5,0.62）
動力學(xué)模型是：-dc/dt=k(c-c‘)，c'指平衡濃度也就是t=5時的值0.62，采用rung-kutta算法在matlab中求出k進而得到此方程。希望得到高手的指點，不勝感激！我的郵箱是yan.feng21@163.com

回復(fù)此樓

» 本帖附件資源列表

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

2012-06-18 21:28:46, 134.98 K

» 猜你喜歡

【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有267人回復(fù)
北工大宋曉艷教授團隊招收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ù)
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有9人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)

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

【分享】《MATLAB 7.0基礎(chǔ)教程》【轉(zhuǎn)載】已經(jīng)有109人回復(fù)
【請教】 2D Gaussian Fitting based on Matlab 已經(jīng)有4人回復(fù)
matlab論文分享（包括神經(jīng)網(wǎng)絡(luò)，遺傳算法，粒子群算法）已經(jīng)有442人回復(fù)
【分享】《matlab解密》教程PDF版【轉(zhuǎn)載】已經(jīng)有705人回復(fù)
求粒子群優(yōu)化算法進行變量優(yōu)選的matlab程序已經(jīng)有20人回復(fù)
MATLAB經(jīng)典算法程序已經(jīng)有531人回復(fù)
matlab的函數(shù)調(diào)用問題已經(jīng)有3人回復(fù)
求回歸方程參數(shù)估計的最優(yōu)化算法matlab代碼已經(jīng)有9人回復(fù)
數(shù)學(xué)建模matlab算法大全已經(jīng)有720人回復(fù)
關(guān)于MATLAB和OpenCV對相同算法的運算速度對比問題已經(jīng)有20人回復(fù)
【求助】matlab 算法求教已經(jīng)有4人回復(fù)
【求助】多目標(biāo)pso算法的MATLAB程序已經(jīng)有3人回復(fù)
【求助】matlab 遺傳算法工具箱繪圖命令是什么？已經(jīng)有6人回復(fù)

憧憬未來

1樓 2012-06-18 21:29:16

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

leguminosae

2樓2012-06-18 21:57:24

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

ymh900423

木蟲 (正式寫手)

應(yīng)助: 13 (小學(xué)生)
金幣: 4331.5
紅花: 1
帖子: 733
在線: 293.2小時
蟲號: 1359669
注冊: 2011-08-03
性別: GG
專業(yè): 有機合成

資料不錯資料不錯資料不錯資料不錯

贊一下

回復(fù)此樓

勝利就在前方

3樓2012-06-18 22:12:57

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

dingd

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

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

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

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

1stOpt求解：

CODE:

Variable t,c;

ODEFunction c'=-k*(c-0.62);

Data;

0,0.75

2,0.7

3,0.65

4,0.64

5,0.62

參數(shù) 最佳估算
-------------------- -------------
k 0.414791664003293

贊一下

回復(fù)此樓

4樓2012-06-19 08:44:43

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

yupeifeng

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 960.3
散金: 34
紅花: 1
帖子: 220
在線: 78.4小時
蟲號: 1421266
注冊: 2011-09-28
專業(yè): 無機材料化學(xué)

引用回帖:

4樓: Originally posted by dingd at 2012-06-19 08:44:43
1stOpt求解：

Variable t,c;
ODEFunction c'=-k*(c-0.62);
Data;
0,0.75
2,0.7
3,0.65
4,0.64
5,0.62

參數(shù) 最佳估算
-------------------- -------------
k 0.414791664003293...

最好能用Rung-kutta算法來求，謝謝啦！

贊一下

回復(fù)此樓

憧憬未來

5樓2012-06-19 09:34:18

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

z蘿莉

銀蟲 (小有名氣)

應(yīng)助: 12 (小學(xué)生)
金幣: 825.7
紅花: 4
帖子: 72
在線: 125.8小時
蟲號: 1615342
注冊: 2012-02-13
性別: MM
專業(yè): 無機非金屬基復(fù)合材料

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
yupeifeng: 金幣+2, ★有幫助 2012-06-19 11:16:23

Sample TextSample Text
我是這么算的，先把公式化成-1*d(c-c')/dt=k(c-c');d(c-c')/(c-c')=-1*k*dt;
ln(c-c')=-k*t;c=exp(-k*t)+c';然后利用matlab非線性最小二乘法進行擬合，方法如下先建立一個名為jack的函數(shù)，內(nèi)容
function c=jack(k,t)
c=exp(-k*t)+0.62;
然后在matlab主窗口中輸入t=[0,2,3,4,5];
c=[0.75,0.7,0.65,0.64,0.62];
[xx,res]=lsqcurvefit(@jack,[0],t,c)
結(jié)果為xx=1.2359
res=0.7571
我把擬合的函數(shù)與原數(shù)據(jù)進行對比發(fā)現(xiàn)t=0的點不太符合，其他的符合的非常好，誤差雖然是0.7571但是主要是t=0點處的誤差。下圖是對比圖其中紅色是擬合出來的圖。

贊一下

回復(fù)此樓

做好自己該做的事

6樓2012-06-19 10:00:56

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

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

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

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

function Untitled10
t0=[0 2 3 4 5];
c0=[0.75 0.7 0.65 0.64 0.62];
[a,b]=ode45(@kfcn,t0,0.75)
polytool(c0,b)
corrcoef(c0,b)
end

function dc=kfcn(t,c)
dc=-(c-0.62);
end

計算得到斜率為0.68左右，相關(guān)系數(shù)0.8879

贊一下

回復(fù)此樓

showmethemoney

7樓2012-06-19 11:50:09

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

flingbaby

鐵桿木蟲 (著名寫手)

小蟲子

應(yīng)助: 4 (幼兒園)
金幣: 6595.2
散金: 10
紅花: 3
帖子: 1209
在線: 404.1小時
蟲號: 904212
注冊: 2009-11-16
性別: GG
專業(yè): 聚合物共混與復(fù)合材料

這個問題解決了嗎？好像前面那些解決方法有些不妥啊，不是RungKutta法啊

贊一下

回復(fù)此樓

讓我早點畢業(yè)吧��！

8樓2014-09-28 16:37:20

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

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

返回列表

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

[考研] 抱歉 +4	田洪有 2026-03-30	4/200	2026-03-30 21:26 by mumin1990
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-30	7/350	2026-03-30 20:25 by 12022246319
[有機交流] 考研調(diào)劑 +8	watb 2026-03-26	8/400	2026-03-30 18:40 by 544594351
[考研] 297求調(diào)劑 +17	田洪有 2026-03-26	18/900	2026-03-30 18:32 by nothing投稿中
[考研] 08工科，295，接受跨專業(yè)調(diào)劑 +3	lmnlzy 2026-03-30	3/150	2026-03-30 17:49 by wangjy2002
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +6	來財aa 2026-03-25	6/300	2026-03-30 09:48 by 青海小西牛
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥求調(diào)劑 +5	篤行致遠. 2026-03-28	5/250	2026-03-30 06:56 by ilovexiaobin
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8	柳溪邊 2026-03-26	8/400	2026-03-29 20:42 by 唐沐兒
[考研] 312，生物學(xué)求調(diào)劑 +3	小譯同學(xué)abc 2026-03-28	3/150	2026-03-28 15:32 by 落睿可思
[考研] 復(fù)試調(diào)劑 +3	raojunqi0129 2026-03-28	3/150	2026-03-28 15:27 by 落�？伤�
[考研] 283求調(diào)劑 +7	A child 2026-03-28	7/350	2026-03-28 12:05 by zllcz
[考研] 285求調(diào)劑 +4	AZMK 2026-03-27	7/350	2026-03-27 20:59 by AZMK
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考博] 26申博 +3	加油沖��！ 2026-03-26	3/150	2026-03-27 15:38 by cls512
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麥當(dāng)當(dāng)
[考研] 348求調(diào)劑 +4	小懶蟲不懶了 2026-03-27	5/250	2026-03-27 12:47 by 果果媽咪
[論文投稿] Journal of Mechanical Science and Technology +3	Russ_ss 2026-03-25	5/250	2026-03-27 10:49 by 陸小果畫大餅
[考研] 0703化學(xué)338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚的貓
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4	gdlf9999 2026-03-24	4/200	2026-03-25 23:01 by boxking200

24小時熱門版塊排行榜

[求助] matlab Rung-kutta算法

» 本帖附件資源列表

» 猜你喜歡

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

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

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

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

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