| 2 | 1/1 | 返回列表 |
| 查看: 488 | 回復: 1 | ||
whqs8426212銅蟲 (正式寫手)
|
[求助]
matlab用ode236元解微分方程組,運行長時間busy ,求大家指點
|
|
主界面鍵入: [ts ys]=ode23(@KineticsEqs,0:60:3600,[0.002 0 0.002 0 0 0]) 然后就沒反應了,一直busy。。。。。。這是什么情況啊, 隔了10分鐘ctrl+C,報錯:Operation terminated by user during ode23 (line 419) 子函數(shù)見二樓 |
銅蟲 (正式寫手)
|
****************************************************子函數(shù)***************************************************** function dCdt = kinetiCsEqs(t,C) % ODE模型方程 k=[76 0.01 (1.0e+7)*2.7 (1.0e+5)*1.58 1.0e+10 3.2*(1.0e+8) 1.2*(1.0e+6) 3.1*(1.0e+5) (1.0e+7) 5*(1.0e+7) 4.2*(1.0e+9) 8.3*(1.0e+5) (1.0e+10) (1.0e+10) 9.7*(1.0e+7)]; %參數(shù) dCAdt= -k(1)*C(1)*C(3)+k(2)*C(2)*C(3)-k(6)*C(1)*C(4)-k(7)*C(1)*C(5)+k(8)*C(2)*C(5)-k(9)*C(1)*C(6)+k(10)*C(2)*C(6); dCBdt= k(1)*C(1)*C(3)+k(2)*C(2)*C(3)-k(6)*C(1)*C(4)-k(7)*C(1)*C(5)+k(8)*C(2)*C(5)-k(9)*C(1)*C(6)+k(10)*C(2)*C(6); dCCdt= -k(1)*C(1)*C(3)-k(2)*C(2)*C(3)-k(3)*C(4)*C(3)+k(7)*C(5)*C(1)+k(9)*C(1)*C(6)+k(11)*C(4)*C(4)+k(12)*C(5)*C(5)-k(15)*C(5)*C(6); dCDdt=k(1)*C(1)*C(3)-k(3)*C(4)*C(3)-k(6)*C(1)*C(4)-k(11)*C(4)*C(4)-k(13)*C(4)*C(5)-k(14)*C(4)*C(6); dCEdt= k(2)*C(2)*C(3)-k(4)*C(5)+k(5)*0.01*C(6)-k(7)*C(1)*C(5)-k(8)*C(2)*C(5)-k(12)*C(5)*C(5)-k(13)*C(4)*C(5)- k(15)*C(5)*C(6); dCFdt= k(4)*C(5)-k(5)*0.001*C(6)-k(9)*C(1)*C(6)-k(10)*C(2)*C(6)-k(14)*C(6)*(1.0e-11)-k(15)*C(5)*C(6); %列出6個反應方程組 dCdt=[ dCAdt;dCBdt;dCCdt;dCDEdt; dCEdt;dCFdt]; %輸出方程組 end **************************************************************************************************************** |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求材料,環(huán)境專業(yè)調劑 +3 | 18567500178 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿吉大化學322求調劑 +3 | 17501029541 2026-03-23 | 4/200 |
|
|
[考研] 384求調劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
[考研] 323求調劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[考研] 石河子大學(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 08工科 320總分 求調劑 +11 | 梨花珞晚風 2026-03-17 | 11/550 |
|
|
[考研] 尋找調劑 +4 | 倔強芒? 2026-03-21 | 4/200 |
|
|
[考研] 求調劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 資源與環(huán)境 調劑申請(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境總分308求調劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 求調劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 一志愿天津大學化學工藝專業(yè)(081702)315分求調劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
[考研] 304求調劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
[考研] 一志愿蘇州大學材料求調劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 材料學碩297已過四六級求調劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 085410人工智能專碩317求調劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 考研調劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|