24小时热门版块排行榜

>论坛更新日志 (3914)
>虫友互识 (1226)
>文献求助 (321)
>导师招生 (209)
>休闲灌水 (191)
>考博 (164)
>教师之家 (103)
>硕博家园 (95)
>博后之家 (89)
>论文投稿 (87)
>找工作 (78)
>基金申请 (74)
>论文道贺祈福 (68)
>考研 (65)
>公派出国 (56)
>绿色求助(高悬赏) (46)

返回列表

whqs8426212

铜虫 (正式写手)

应助: 5 (幼儿园)
金币: 13.2
帖子: 782
在线: 685.7小时
虫号: 2572504

[交流] 自己写的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)

回复此楼