| 查看: 1743 | 回复: 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) |
回复此楼» 猜你喜欢
如果工作了想读博，可以边工作边读全日制嘛？
已经有6人回复
-
青B发送上会通知了吗
已经有9人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有6人回复
-
论文撤稿了
已经有8人回复
-
博士申请
已经有7人回复
-
化学专业申博
已经有4人回复
-
评审有感
已经有27人回复
-
河北省自然科学基金
已经有9人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
-
26/27申博自荐
已经有9人回复
» 抢金币啦！回帖就可以得到:
-
广西锰资源高值化利用重点实验室 2026年博士人才招聘通告
+1/969
-
祝大家学业有成，工作顺利，生意兴隆
+1/945
-
中科院深圳先进技术研究院集成电路先进封装博士后招聘
+1/283
-
坐标深圳，愿遇良缘
+1/213
-
山东第一医科大学泰山学者团队2026年招收副教授、博士后（不限专业）
+1/81
-
真诚寻找另一半 觉得合适可以相互了解
+1/81
-
南京农业大学工学院博士招生 1个名额
+1/74
-
河南省医学科学院王宁利院士科研团队2026年博士、博后、硕士招聘
+1/36
-
海南医科大学肖华教授(国家级青年人才)招收博士生1名---报名从速
+4/32
-
【王宁利院士团队招聘】眼科学博士3人
+1/32
-
安徽农业大学招聘师资博士后以及人才引进教师
+1/28
-
紧急招收2026年秋季入学博士生1名（湘潭大学 固体废弃物低碳利用湖南省工程研究中心）
+1/15
-
紧急招收2026年秋季入学博士生1名（湘潭大学 固体废弃物低碳利用湖南省工程研究中心）
+1/13
-
哪位大神能使用lammps软件的reaxff力场计算煤和塑料热解的分子动力学，有偿
+1/11
-
西安交通大学陕西省科技创新团队招聘2027年免试保送生和考研生，2-4个名额
+1/11
-
紧急招收2026年秋季入学博士生1名（河北工大/北京科技大学联合 增材制造/生物材料）
+1/6
-
杭州师范大学博士后招聘
+1/6
-
北航国新院季梦奇副教授招聘博士(智能感知、立体视觉等)
+1/4
-
重庆新桥医院临床研究助理招聘
+1/4
-
博士去绵阳华丰科技怎么样，待遇如何，这个企业靠谱吗？
+3/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): 谢谢参与