版塊導(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)劑公告

返回列表

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

qinzhong6138

至尊木蟲(chóng) (著名寫(xiě)手)

灰�；页１康男”康�

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

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

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

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

）和文獻(xiàn)上面的公式貼上來(lái)，原始的文獻(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ì)寫(xiě)了）
% ------------------------------------------------------------------

111.jpg

222.jpg

333.jpg

444.jpg

回復(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 : Product_sampling_during_transient_continuou.pdf

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

» 猜你喜歡

計(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)有4人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有79人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有29人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有26人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
UJN物理學(xué)專業(yè)調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名已經(jīng)有0人回復(fù)

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

求教一個(gè)不可壓流動(dòng)投影法的3階Runge-Kutta顯式時(shí)間推進(jìn)的問(wèn)題已經(jīng)有1人回復(fù)
請(qǐng)問(wèn)Runge-Kutta算法的步長(zhǎng) 已經(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

至尊木蟲(chóng) (著名寫(xiě)手)

灰�；页１康男”康�

應(yīng)助: 9 (幼兒園)
金幣: 22749.4
散金: 3540
紅花: 14
帖子: 2194
在線: 2261.1小時(shí)
蟲(chóng)號(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ā)上來(lái)。
如果只是要程序代碼的話，參見(jiàn)：http://www.gaoyang168.com/bbs/viewthread.php?tid=7575773&authorid=1122189

不懂怎么傳附件，

，只好用網(wǎng)盤(pán)上傳了，謝謝月版主，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的回帖

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

月只藍(lán)

主管區(qū)長(zhǎng) (職業(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í)
蟲(chóng)號(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ā)上來(lái)。
如果只是要程序代碼的話，參見(jiàn)：http://www.gaoyang168.com/bbs/viewthread.php?tid=7575773&authorid=1122189

贊一下

回復(fù)此樓

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

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

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

qinzhong6138

至尊木蟲(chóng) (著名寫(xiě)手)

灰�；页１康男”康�

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

引用回帖:

先謝謝版主了，我不懂matlab，所以對(duì)matlab代碼也不是很懂。了解到1stOpt對(duì)解微分方程也很好，ps:dingd是這方面的行家的，多請(qǐng)教了。昨晚花了幾個(gè)小時(shí)，大致了解一下這個(gè)軟件，因此，自己編寫(xiě)了一些1stOpt格式的程序和數(shù)據(jù)，不知道您對(duì)這個(gè)了解嗎？以后我會(huì)經(jīng)常遇到類似的非線性擬合的數(shù)據(jù)，如果1stOpt真的很方便的話，我打算買(mǎi)一套，支持國(guó)產(chǎn)，于己于人都方便，不過(guò)這個(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

至尊木蟲(chóng) (著名寫(xiě)手)

灰�；页１康男”康�

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

帖子的數(shù)據(jù)出錯(cuò)了，重新上傳，不好意思

鏈接: http://pan.baidu.com/s/1jGJucmi 密碼: pkqv

贊一下

回復(fù)此樓

加油

5樓2014-07-13 11:50:57

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

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

普通表情龍兔虎貓

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

[考研] 289求調(diào)劑 +23	新時(shí)代材料 2026-03-27	26/1300	2026-04-02 10:29 by oooqiao
[考研] 324分 085600材料與化工 +20	呆鵝oor 2026-03-27	20/1000	2026-04-02 10:13 by oooqiao
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +14	傾____蕭 2026-03-27	15/750	2026-04-02 09:21 by olim
[考研] 274求調(diào)劑 +8	薛定諤的虎。 2026-04-01	8/400	2026-04-02 09:21 by 晴空210210
[考研] 求生物學(xué)調(diào)劑 +8	15172915737 2026-04-01	8/400	2026-04-02 06:49 by ilovexiaobin
[考研] 310分求調(diào)劑 +4	成功上岸wang 2026-04-01	4/200	2026-04-01 20:35 by liu823948201
[考研] 材料科學(xué)與工程339求調(diào)劑 +11	hyz0119 2026-03-31	12/600	2026-04-01 18:40 by 偉大河北
[考研] 353求調(diào)劑 +4	拉鉤不許變 2026-04-01	4/200	2026-04-01 18:10 by 記事本2026
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥求調(diào)劑 +6	篤行致遠(yuǎn). 2026-03-28	6/300	2026-04-01 09:28 by oooqiao
[考研] 求調(diào)劑：085600材料與化工，考材科基，總分319 +17	678lucky 2026-03-31	21/1050	2026-04-01 01:40 by 1018329917
[考研] 335求調(diào)劑 +3	321* 2026-03-31	4/200	2026-04-01 00:00 by 321*
[考研] 求收留 +8	1943443204 2026-03-28	8/400	2026-03-31 15:00 by -迷了路啊路
[考研] 287求調(diào)劑 +17	land xuxu 2026-03-26	17/850	2026-03-31 11:16 by Zzxxxs
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +6	大饞小子 2026-03-29	6/300	2026-03-31 10:03 by 氯化亞硝酰
[考研] 085602化工求調(diào)劑（331分） +8	111@127 2026-03-30	8/400	2026-03-30 21:23 by 研究僧導(dǎo)導(dǎo)
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6	零八# 2026-03-28	6/300	2026-03-29 21:20 by nanaliuyun
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3	@taotao 2026-03-28	3/150	2026-03-28 10:26 by JourneyLucky
[考研] 0856調(diào)劑 +5	求求讓我有書(shū)讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 一志愿吉大071010，316分求調(diào)劑 +3	xgbiknn 2026-03-27	3/150	2026-03-27 10:36 by guoweigw