| 查看: 1661 | 回复: 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人回复
-
国家基金申请书模板内插入图片不可调整大小?
已经有6人回复
-
招博士
已经有6人回复
-
限项规定
已经有8人回复
-
国家级人才课题组招收2026年入学博士
已经有5人回复
-
Fe3O4@SiO2合成
已经有6人回复
-
青年基金C终止
已经有4人回复
-
青椒八年已不青，大家都被折磨成啥样了？
已经有7人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
哈尔滨工程大学青岛创新发展基地招聘青年教师
+1/474
-
海南大学海洋清洁能源创新团队（光电催化方向）招收2026年博士研究生（第二批）
+5/305
-
nanomaterials特刊：纳米材料在电化学能源领域的应用（SCI 中科院三区，IF4.5）征稿
+5/210
-
【2026/2027 哈工大计算机类博士招生】
+1/78
-
广东工业大学自动化学院鲁仁全教授团队刘勇华老师招收2026年博士研究生（申请制）
+1/73
-
南京林业大学特聘教授团队招聘博后和2026博士研究生
+1/72
-
自驾淄博回四川遂宁过年，寻找小伙伴一起
+1/56
-
上海市“光探测材料与器件”工程技术研究中心（上海应用技术大学）招聘优秀研究人员
+1/46
-
教育部重点实验室和清华大学某国家重点实验室，联合培养硕生、博生，并长期招博士后
+1/33
-
26申博自荐-电磁屏蔽材料方向
+3/32
-
青岛大学化学院泰山学者冯佳课题组申请考核制博士招生
+1/13
-
博士后, 博士招生 美国大纽约地区 NJIT
+1/11
-
山东大学集成电路学院王凌云研究员招收2026年硕士生及联合培养硕士生
+1/6
-
中山大学柔性电子学院黄维院士团队诚招柔性可穿戴电子方向博士生（2026年9月入学）
+1/6
-
华南理工大学宋波教授招聘材料和化学方向博士后（长期有效）
+1/6
-
华南理工大学宋波教授招收2026年博士生（二氧化碳转化方向优先）
+1/3
-
福建师范大学柔性电子学院 院士团队招2026级博士 光电器件、发光传感忆阻器
+1/3
-
广东以色列理工学院博士/硕士招生-通过稀疏观测用数据驱动方法预测湍流
+1/2
-
山东大学集成电路学院太赫兹团队博士招生(2026年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): 谢谢参与