當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

qinzhong6138

灰�；页１康男”康�

應(yīng)助: 9 (幼兒園)
金幣: 22749.4
散金: 3540
紅花: 14
帖子: 2194
在線: 2261.1小時(shí)
蟲號(hào): 585846
注冊(cè): 2008-08-08
專業(yè): 功能與智能高分子

[求助] 請(qǐng)教一下各位大神怎么求解 Runge–Kutta method的微分方程的各個(gè)系數(shù) 已有1人參與

最近在做一點(diǎn)研究，但是被幾個(gè)微分方程難住了，已經(jīng)頭疼好幾天了

。現(xiàn)在發(fā)帖子求助一下各位牛人，也請(qǐng)伸出援手，幫一下我啊。
我的目的是用matlab求解得到k1--k8的值，一共是6個(gè)微分方程的8個(gè)系數(shù)，如圖所示。請(qǐng)大家指點(diǎn)一下�。∥覇杽e人要了一段代碼，但是這段代碼還不夠完整，現(xiàn)將其貼上來（本人不會(huì)寫

）和文獻(xiàn)上面的公式貼上來，原始的文獻(xiàn)也在附件。我在這里先謝謝大家了啊�。。�

600金幣表示一下心意，不夠的話，可以加啊

function KineticsEst1_int
clear all
clc
global CAm CBm CCm CDm CEm
t = [ ];
CAm = []’;
CBm = []';
CCm = []';
CDm=[]';
CEm=[]';
beta0 = [ ];
tspan = [ ];
CA0 = 0.11;
CB0 = 0;
CC0 = 0;
CD0= 0.03;
CE0= 0.106;
C0=[0.31 0 0 0.02 0.10];
[beta,resnorm,resid,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@OptObjFunc,beta0,[0 0 0 0 0 0],[],[],tspan,[CA0 CB0 CC0 CD0 CE0])
ci = nlparci(beta,resid,jacobian)
[t4plot CA4plot] = ode45(@KineticsEqs,[tspan(1) tspan(end)],[CA0 CB0 CC0 CD0 CE0],[],beta);
plot(tspan,CAm,'bo',t4plot,CA4plot(:,1),'k-')
hold on
plot(tspan,CBm,'bo',t4plot,CA4plot(:,2),'k-')
hold on
plot(tspan,CCm,'bo',t4plot,CA4plot(:,3),'k-')
legend('Exp','Model')
xlabel('time, h')
ylabel('concentration, mol/L')
legend('Exp','Model')
fprintf('Estimated Parameters:\n')
fprintf('\tk1 = %.20f ¡À %.20f\n', beta(1),ci(1,2)-beta(1))
fprintf('\tk2 = %.20f ¡À %.20f\n', beta(2),ci(2,2)-beta(2))
fprintf('\tk3 = %.20f ¡À %.20f\n',beta(3),ci(3,2)-beta(3))
fprintf('\tk4 = %.20f ¡À %.20f\n',beta(4),ci(4,2)-beta(4))
fprintf('\tk5 = %.20f ¡À %.20f\n',beta(5),ci(5,2)-beta(5))
fprintf('\tk6 = %.20f ¡À %.20f\n',beta(6),ci(6,2)-beta(6))
（不會(huì)寫了）
% ------------------------------------------------------------------

111.jpg

222.jpg

333.jpg

444.jpg

回復(fù)此樓

» 本帖附件資源列表

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

2014-07-12 22:08:43, 2.11 M

» 猜你喜歡

中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有236人回復(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é)工程國際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有1人回復(fù)

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

求教一個(gè)不可壓流動(dòng)投影法的3階Runge-Kutta顯式時(shí)間推進(jìn)的問題已經(jīng)有1人回復(fù)
請(qǐng)問Runge-Kutta算法的步長已經(jīng)有3人回復(fù)
【資源】Runge-Kutta解ODE的一些程序及應(yīng)用（from MathWorks）已經(jīng)有3人回復(fù)
【求助】如何用Runge-Kutta迭代求解二階常微分方程組【已解決】已經(jīng)有9人回復(fù)

加油

1樓 2014-07-12 22:05:49

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

qinzhong6138

至尊木蟲 (著名寫手)

灰常灰常笨的小笨蛋

應(yīng)助: 9 (幼兒園)
金幣: 22749.4
散金: 3540
紅花: 14
帖子: 2194
在線: 2261.1小時(shí)
蟲號(hào): 585846
注冊(cè): 2008-08-08
專業(yè): 功能與智能高分子

引用回帖:

8樓: Originally posted by 月只藍(lán) at 2014-07-13 13:39:17
1stopt軟件要高版本的才能做常微分方程的擬合，而且這款軟件的擬合效果常比普通編寫的MATLAB程序強(qiáng)。你可以把完整的1stopt的程序?qū)懞�，讓dingd專家跑一下。...

我已經(jīng)寫了，在七樓的網(wǎng)盤鏈接中，不知道版主有空幫我跑下matlab嗎？

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

加油

9樓2014-07-13 14:28:13

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

查看全部 13 個(gè)回答

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68128.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
qinzhong6138: 金幣+20, ★★★很有幫助, 先謝謝版主，等待你施以援手！ 2014-07-13 11:45:16
qinzhong6138: 金幣+80, ★★★很有幫助, 謝謝版主的幫助，我從別的地方找到相應(yīng)解決辦法了，該帖子關(guān)了吧，謝謝啦 2014-07-16 14:52:15

你是要擬合常微分方程組的參數(shù)吧。
方便的話，把數(shù)據(jù)也一起發(fā)上來。
如果只是要程序代碼的話，參見：http://www.gaoyang168.com/bbs/viewthread.php?tid=7575773&authorid=1122189

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

2樓2014-07-12 23:13:57

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

qinzhong6138

至尊木蟲 (著名寫手)

灰�；页１康男”康�

應(yīng)助: 9 (幼兒園)
金幣: 22749.4
散金: 3540
紅花: 14
帖子: 2194
在線: 2261.1小時(shí)
蟲號(hào): 585846
注冊(cè): 2008-08-08
專業(yè): 功能與智能高分子

引用回帖:

2樓: Originally posted by 月只藍(lán) at 2014-07-12 23:13:57
你是要擬合常微分方程組的參數(shù)吧。
方便的話，把數(shù)據(jù)也一起發(fā)上來。
如果只是要程序代碼的話，參見：http://www.gaoyang168.com/bbs/viewthread.php?tid=7575773&authorid=1122189

先謝謝版主了，我不懂matlab，所以對(duì)matlab代碼也不是很懂。了解到1stOpt對(duì)解微分方程也很好，ps:dingd是這方面的行家的，多請(qǐng)教了。昨晚花了幾個(gè)小時(shí)，大致了解一下這個(gè)軟件，因此，自己編寫了一些1stOpt格式的程序和數(shù)據(jù)，不知道您對(duì)這個(gè)了解嗎？以后我會(huì)經(jīng)常遇到類似的非線性擬合的數(shù)據(jù)，如果1stOpt真的很方便的話，我打算買一套，支持國產(chǎn)，于己于人都方便，不過這個(gè)軟件價(jià)格有些肉疼啊，年輕老師，收入低。版主，您方便的話，能麻煩您用matlb和1stOpt都能幫我算一下嗎，數(shù)據(jù)馬上上傳。@月只藍(lán)，@ dingd

贊一下

回復(fù)此樓

加油

3樓2014-07-13 11:41:48

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

qinzhong6138

至尊木蟲 (著名寫手)

灰�；页１康男”康�

應(yīng)助: 9 (幼兒園)
金幣: 22749.4
散金: 3540
紅花: 14
帖子: 2194
在線: 2261.1小時(shí)
蟲號(hào): 585846
注冊(cè): 2008-08-08
專業(yè): 功能與智能高分子

引用回帖:

不懂怎么傳附件，

，只好用網(wǎng)盤上傳了，謝謝月版主，D專家有空的話，也麻煩幫助一下啊。@月只藍(lán)，@ dingd
鏈接: http://pan.baidu.com/s/1dD4Bf7V 密碼: 6mcn

贊一下

回復(fù)此樓

加油

4樓2014-07-13 11:47:26

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

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

[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +5	Shern—- 2026-03-28	5/250	2026-03-29 23:32 by 飛行日記西
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +5	傾____蕭 2026-03-27	6/300	2026-03-29 16:35 by 唐沐兒
[考研] 329求調(diào)劑 +10	鈕恩雪 2026-03-25	10/500	2026-03-29 13:32 by peike
[考研] 一志愿211 初試270分求調(diào)劑 +8	谷雨上岸 2026-03-23	9/450	2026-03-29 01:00 by 我是小康
[考研] 11408軟件工程求調(diào)劑 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 21:50 by zhq0425
[考研] 316求調(diào)劑 +7	江辭666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4	小張zxy張 2026-03-26	8/400	2026-03-28 19:25 by lbsjt
[考研] 085701環(huán)境工程求調(diào)劑 +9	多久上課 2026-03-27	9/450	2026-03-28 03:58 by fmesaito
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4	Ripcord上岸 2026-03-27	4/200	2026-03-27 17:00 by zhyzzh
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級(jí)，六級(jí)521，求調(diào)劑 +4	美味可樂雞翅 2026-03-26	4/200	2026-03-27 15:27 by 星空星月
[考研] 359求調(diào)劑 +4	王了個(gè)楠 2026-03-25	4/200	2026-03-27 08:43 by 不吃魚的貓
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4	鋼鐵大炮 2026-03-24	4/200	2026-03-26 16:28 by dick_runner
[考研] 281求調(diào)劑 +6	Koxui 2026-03-24	7/350	2026-03-26 15:37 by 無際的草原
[考研] 求b區(qū)院校調(diào)劑 +4	周56 2026-03-24	5/250	2026-03-25 17:12 by yishunmin
[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 086003食品工程求調(diào)劑 +6	淼淼111 2026-03-24	6/300	2026-03-25 10:29 by 3Strings
[考研] 化工專碩求調(diào)劑 +3	question挽風(fēng) 2026-03-24	3/150	2026-03-24 18:48 by jhhcooi
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106