| 查看: 1685 | 回复: 12 | |||
[交流]
自己写的ode23解微分方程，运行长时间busy
|
|||
|
****************************************************子函数***************************************************** 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)]; %参数 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 **************************************************************************************************************** 在matlab主界面键入： [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) |
回复此楼» 猜你喜欢
【博士招生】太原理工大学2026化工博士
已经有4人回复
-
什么是人一生最重要的？
已经有5人回复
-
280求调剂
已经有3人回复
-
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？
已经有8人回复
-
面上可以超过30页吧？
已经有11人回复
-
版面费该交吗
已经有15人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有18人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有10人回复
» 抢金币啦！回帖就可以得到:
-
单曲循环久了很多都会变
+1/399
-
侯旭课题组(化学、化工、环境、能源相关方向)欢迎你的加入
+1/182
-
87 年东北人定居苏州，诚寻携手余生的你
+1/164
-
【最后机会】深圳大学26级土木工程博士急招
+5/160
-
2026时光标本=已经注册了 5608 天，合计 16 年
+1/138
-
【最后机会】深圳大学2026级土木工程博士急招
+5/120
-
澳门理工大学人工智能智慧康养方向26 年9月入学博士招生有奖学金
+1/35
-
澳门理工大学人工智能智慧康养方向26 年9月入学博士招生 奖学金
+1/31
-
南通大学生物医药方向国家级人才团队招收“申请-考核”制博士研究生20260222
+2/30
-
211大学【2026学博】补招
+1/28
-
深圳理工大学刘鑫课题组高薪招聘博士后
+1/14
-
中国科学院金属研究所何杰团队招收2026年9月入学博士生（具有铝合金研究经验者优先）
+1/12
-
诚邀有志之士加入江苏大学环境学院合成生物学课题组！(长期接收学硕专硕、硕博连读)
+1/11
-
中科院合肥物质研究院招聘水系锌电博士后（长期有效）诚聘海内外优秀博士后
+1/7
-
【2026年】【材料力学方向】中科院力学所招1名硕士
+1/7
-
【全奖博士招生】美国科罗拉多大学科罗拉多斯普林斯分校----固体与结构力学
+1/7
-
蛋白结构质谱实验室202609入学博士生（1名，202603月底报名）
+1/6
-
集美大学轮机工程学院2026年博士研究生招生
+1/2
-
新加坡南洋理工大学- 光电/ 智能传感/ 脑机接口方向 博士后
+1/2
-
中国科学技术大学 工程科学学院 国家级人才团队 诚聘博士后
+1/1
3楼2014-08-03 21:43:53
7楼2014-08-03 22:11:42
9楼2014-08-03 22:19:54
13楼2014-08-08 17:26:44
简单回复
pzx10042楼
2014-08-03 21:37
回复
whqs8426212(金币+1): 谢谢参与
祝福 [ 发自小木虫客户端 ]
2014-08-03 21:58
回复
whqs8426212(金币+1): 谢谢参与
chyzm5楼
2014-08-03 21:58
回复
whqs8426212(金币+1): 谢谢参与
tzynew6楼
2014-08-03 22:03
回复
whqs8426212(金币+1): 谢谢参与
SCI-ICS8楼
2014-08-03 22:12
回复
whqs8426212(金币+1): 谢谢参与
祝福
shangshui12310楼
2014-08-03 22:19
回复
whqs8426212(金币+1): 谢谢参与
84317229011楼
2014-08-03 22:25
回复
whqs8426212(金币+1): 谢谢参与
232898093312楼
2014-08-03 22:30
回复
whqs8426212(金币+1): 谢谢参与