 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2289 | 回復(fù): 12 | |||||
qinzhong6138至尊木蟲 (著名寫手)
灰�；页１康男”康�
|
[求助]
請教一下各位大神怎么求解 Runge–Kutta method的微分方程的各個系數(shù) 已有1人參與
|
||||
最近在做一點研究，但是被幾個微分方程難住了，已經(jīng)頭疼好幾天了 �，F(xiàn)在發(fā)帖子求助一下各位牛人，也請伸出援手，幫一下我啊。我的目的是用matlab求解得到k1--k8的值，一共是6個微分方程的8個系數(shù)，如圖所示。請大家指點一下��！我問別人要了一段代碼，但是這段代碼還不夠完整，現(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)) （不會寫了） % ------------------------------------------------------------------ 111.jpg 222.jpg 333.jpg 444.jpg |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : Product_sampling_during_transient_continuou.pdf
2014-07-12 22:08:43, 2.11 M
» 猜你喜歡
2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有168人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有5人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有30人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有56人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求教一個不可壓流動投影法的3階Runge-Kutta顯式時間推進的問題
已經(jīng)有1人回復(fù)
- 請問Runge-Kutta算法的步長
已經(jīng)有3人回復(fù)
- 【資源】Runge-Kutta解ODE的一些程序及應(yīng)用（from MathWorks）
已經(jīng)有3人回復(fù)
- 【求助】如何用Runge-Kutta迭代求解二階常微分方程組【已解決】
已經(jīng)有9人回復(fù)
月只藍
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(qū)
|
@dingd，不好意思，樓主給出的方程中有錯誤，10樓的程序就不用看了，麻煩跑一下下面的程序吧： Parameters K1,K2,K3,K4,K5,K6,K7,K8； Variable t,T,D,M,G,W,F; ODEFunction T’= -K7*T*M+K8*D^2-K1*T*W+K2*D*F； D’= 2*K7*T*M-2*K8*D^2-K1*T*W+K2*D*F-K3*D*W+K4*M*F； M’=K7*T*M+K8*D^2+K3*D*W-K4*M*F-K5*M*W+K6*G*F； W’=-K1*T*W+K2*D*F-K3*D*W+K4*M*F-K5*M*W+K6*G*F； F’=K1*T*W-K2*D*F+K3*D*W+K4*M*F+K5*M*W+K6*G*F； G’= K5*M*W-K6*G*F； data； 0 99.04 0.81 0 0 0.15 0 1 68.34 18.69 0.96 3.9633 12.01 11.86 2 53.58 23.71 2.05 6.8178 20.66 20.51 3 46.46 26.96 3.02 7.7748 23.56 23.41 4 40.76 28.83 4.23 8.6394 26.18 26.03 5 38.11 29.21 4.52 9.2928 28.16 28.01 6 35.7 29.37 4.75 9.9594 30.18 30.03 8 32.51 30.29 4.80 10.692 32.4 32.25 10 31.04 30.78 4.99 10.9527 33.19 33.04 12 30.81 31.24 5.11 10.8372 32.84 32.69 24 28.84 29.08 6.09 11.8767 35.99 35.84 |
» 本帖已獲得的紅花（最新10朵）
qinzhong613811樓2014-07-14 10:18:56
月只藍
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(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
感謝參與，應(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 |
2樓2014-07-12 23:13:57
qinzhong6138
至尊木蟲 (著名寫手)
灰常灰常笨的小笨蛋
- 應(yīng)助: 9 (幼兒園)
- 金幣: 22749.4
- 散金: 3540
- 紅花: 14
- 帖子: 2194
- 在線: 2261.1小時
- 蟲號: 585846
- 注冊: 2008-08-08
- 專業(yè): 功能與智能高分子
 先謝謝版主了，我不懂matlab，所以對matlab代碼也不是很懂。了解到1stOpt對解微分方程也很好，ps:dingd是這方面的行家的，多請教了。昨晚花了幾個小時，大致了解一下這個軟件，因此，自己編寫了一些1stOpt格式的程序和數(shù)據(jù)，不知道您對這個了解嗎？以后我會經(jīng)常遇到類似的非線性擬合的數(shù)據(jù)，如果1stOpt真的很方便的話，我打算買一套，支持國產(chǎn)，于己于人都方便，不過這個軟件價格有些肉疼啊，年輕老師，收入低。版主，您方便的話，能麻煩您用matlb和1stOpt都能幫我算一下嗎，數(shù)據(jù)馬上上傳。@月只藍，@ dingd |
3樓2014-07-13 11:41:48
qinzhong6138
至尊木蟲 (著名寫手)
灰常灰常笨的小笨蛋
- 應(yīng)助: 9 (幼兒園)
- 金幣: 22749.4
- 散金: 3540
- 紅花: 14
- 帖子: 2194
- 在線: 2261.1小時
- 蟲號: 585846
- 注冊: 2008-08-08
- 專業(yè): 功能與智能高分子
|
不懂怎么傳附件，  ，只好用網(wǎng)盤上傳了，謝謝月版主，D專家有空的話，也麻煩幫助一下啊。@月只藍，@ dingd鏈接: http://pan.baidu.com/s/1dD4Bf7V 密碼: 6mcn |
4樓2014-07-13 11:47:26
qinzhong6138
至尊木蟲 (著名寫手)
灰常灰常笨的小笨蛋
- 應(yīng)助: 9 (幼兒園)
- 金幣: 22749.4
- 散金: 3540
- 紅花: 14
- 帖子: 2194
- 在線: 2261.1小時
- 蟲號: 585846
- 注冊: 2008-08-08
- 專業(yè): 功能與智能高分子
5樓2014-07-13 11:50:57
qinzhong6138
至尊木蟲 (著名寫手)
灰�；页１康男”康�
- 應(yīng)助: 9 (幼兒園)
- 金幣: 22749.4
- 散金: 3540
- 紅花: 14
- 帖子: 2194
- 在線: 2261.1小時
- 蟲號: 585846
- 注冊: 2008-08-08
- 專業(yè): 功能與智能高分子
6樓2014-07-13 11:54:19
qinzhong6138
至尊木蟲 (著名寫手)
灰常灰常笨的小笨蛋
- 應(yīng)助: 9 (幼兒園)
- 金幣: 22749.4
- 散金: 3540
- 紅花: 14
- 帖子: 2194
- 在線: 2261.1小時
- 蟲號: 585846
- 注冊: 2008-08-08
- 專業(yè): 功能與智能高分子
qinzhong6138: 回帖置頂 2014-07-13 12:09:04
|
數(shù)據(jù)在這樓，前面的方程沒寫對，重新編寫和上傳，對不起了！ 還是找不到論壇的附件上傳  鏈接: http://pan.baidu.com/s/1jGsVAcQ 密碼: s14d |
7樓2014-07-13 12:08:56
月只藍
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(qū)
8樓2014-07-13 13:39:17
qinzhong6138
至尊木蟲 (著名寫手)
灰�；页１康男”康�
- 應(yīng)助: 9 (幼兒園)
- 金幣: 22749.4
- 散金: 3540
- 紅花: 14
- 帖子: 2194
- 在線: 2261.1小時
- 蟲號: 585846
- 注冊: 2008-08-08
- 專業(yè): 功能與智能高分子
9樓2014-07-13 14:28:13
月只藍
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計算模擬區(qū)
|
@dingd 麻煩幫忙跑一下： Parameters K1,K2,K3,K4,K5,K6,K7,K8； Variable t,T,D,M,G,W,F; ODEFunction T’= -K7*T*M+K8*D^2-K1*T*W+K2*D*F； D’= 2*K7*T*M-2*K8*D^2-K1*T*W+K2*D*F-K3*D*W+K4*M*F； M’=K7*T*M+K8*D^2+K3*D*H-K4*M*F-K5*M*W+K6*G*F； W’=-K1*T*W+K2*D*F-K3*D*W+K4*M*F-K5*M*W+K6*G*F； F’=K1*T*W-K2*D*F+K3*D*W+K4*M*F+K5*M*W+K6*G*F； G’= K5*M*W-K6*G*F； data； 0 99.04 0.81 0 0 0.15 0 1 68.34 18.69 0.96 3.9633 12.01 11.86 2 53.58 23.71 2.05 6.8178 20.66 20.51 3 46.46 26.96 3.02 7.7748 23.56 23.41 4 40.76 28.83 4.23 8.6394 26.18 26.03 5 38.11 29.21 4.52 9.2928 28.16 28.01 6 35.7 29.37 4.75 9.9594 30.18 30.03 8 32.51 30.29 4.80 10.692 32.4 32.25 10 31.04 30.78 4.99 10.9527 33.19 33.04 12 30.81 31.24 5.11 10.8372 32.84 32.69 24 28.84 29.08 6.09 11.8767 35.99 35.84 |
10樓2014-07-14 10:10:59
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 qinzhong6138 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 11/550 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +4 | 蔡大美女 2026-03-13 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|