| 查看: 1656 | 回复: 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) |
回复此楼» 猜你喜欢
请教限项目规定
已经有5人回复
-
拟解决的关键科学问题还要不要写
已经有8人回复
-
最失望的一年
已经有16人回复
-
存款400万可以在学校里躺平吗
已经有33人回复
-
求助一下有机合成大神
已经有3人回复
-
求推荐英文EI期刊
已经有5人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
疑惑？
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
南京大学智能驱动与感知材料实验室诚招硕士/科研助理
+1/271
-
东北石油大学三亚海洋油气研究院|地学硕士|地质资源与地质工程、地质学、地质工程等
+1/189
-
海南大学专任教师招聘：海洋工程与海洋智能装备相关方向【长期有效】
+2/132
-
中南林科大-自然资源与生态环境管理-博士招生（导师本人发布）
+1/88
-
大连理工大学智能系统实验室优秀硕博研究生招生
+1/81
-
中国地质大学（武汉）地质学、地质资源与地质工程、资源与环境方向招生，请尽快联系！
+1/70
-
211双一流北京工业大学招计算机、AI、自动控制、电子信息等方向博士生（长期有效）
+1/40
-
上海交通大学 碳纳米管、阵列、器件应用博士后招聘
+1/34
-
华中科技大学周英教授招博后
+1/34
-
北欧岗位制博士招生（锂离子电池正极材料和粘结剂）
+1/33
-
南京林业大学木质纤维功能材料国际联合创新中心招收2026级博士生（申请-考核制）
+2/20
-
2026申博自荐 本硕双一流学科 纳米药物递送方向 一篇一区TOP 两个国家奖学金
+1/18
-
云南大学智能光电子课题组招收申请-考核博士研究生1名（光电、光学）
+1/16
-
2025版《中国药典》方法测定二甲基亚砜含量偏高
+1/6
-
南方科技大学田雷蕾课题组招收2026年博士生
+1/5
-
2026年南方医科大学基础医学院李琳课题组招收“申请-审核”生物信息学博士研究生
+1/4
-
意大利米兰理工大学Lorenza Petrini教授招收玛丽居里MEDALLOY博士
+1/3
-
生殖医学与子代健康全国重点实验室华鹏课题组招收博士后及研究生(长期有效)
+1/2
-
重庆大学诚招2026年生物材料方向博士生
+1/1
-
【博士后/科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+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): 谢谢参与