| 查看: 1798 | 回复: 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) |
回复此楼» 猜你喜欢
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
-
filecode
已经有4人回复
-
求各位大神看下
已经有27人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有11人回复
-
一个有机合成实验室都需要哪些设备？
已经有11人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有30人回复
-
2027年申博
已经有4人回复
-
好消息？这个有何含义？？？
已经有10人回复
-
【2027博士申请】纳米药物递送方向
已经有4人回复
» 抢金币啦！回帖就可以得到:
-
中科院招聘博士后及科研助理
+2/244
-
山东征女友，坐标济南
+1/168
-
招聘机械领域专利代理人
+1/100
-
3-羟基丙腈（CAS：109-78-4）
+1/90
-
丁烯二醇（CAS: 497-06-3）
+1/89
-
三(三甲基硅基)甘油醚（CAS：6787-10-6 ）
+1/89
-
十二烷基溴化吡啶（CAS：104-73-4）
+1/89
-
三氟甲基亚磺酸锂
+1/88
-
2-氯-1,3,2-二氧磷杂环戊烷（CAS：822-39-9）
+1/88
-
具身智能机器人公司招聘硬件工程师岗位 预算100万内
+1/81
-
化工学院邱惠斌教授联合同仁医院招聘博士后
+1/65
-
中国科学院西双版纳热带植物园人事处招聘
+1/47
-
中科院东莞材料所依托中科大招收2027届推免生
+1/33
-
上海交通大学诚聘有机化学方向博士生/科研助理
+1/33
-
双一流高校招收博士（现代技术管理/环境技术经济管理）及硕士研究生（长期有效）
+1/16
-
希望大家都高中！！
+1/10
-
北理工集成电路杰青团队 | 诚招科助理
+1/4
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/3
-
北京师范大学（珠海校区）赵亚然副教授招收2027年电催化推免硕士/博士研究生
+1/3
-
国自然发榜在即：此刻的你，是什么心情？
+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): 谢谢参与