| 查看: 491 | 回复: 1 | |||
whqs8426212铜虫 (正式写手)
|
[求助]
matlab用ode236元解微分方程组，运行长时间busy ，求大家指点
|
|
主界面键入： [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) 子函数见二楼 |
回复此楼» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有172人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一下各位大神怎么求解 Runge–Kutta method的微分方程的各个系数
已经有12人回复
- 呼叫版主，在线紧急求助，关于matlab中微分方程组参数拟合得问题！
已经有12人回复
- 用matlab求解一个隐三角函数方程
已经有10人回复
- 请问这个微分方程怎么解，右边积分函数中被积函数是不是1
已经有5人回复
- 非线性二阶微分方程组求解，matlab
已经有12人回复
- Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！
已经有4人回复
- MATLAB解含三角函数方程，求大家指点！
已经有4人回复
- 利用matlab求一阶变系数微分方程组的数值解。
已经有7人回复
- matlab如何求解含有bessel方程中的参数
已经有5人回复
- 求助，matlab求解二元二阶的常微分方程组
已经有9人回复
- Matlab: 用bvp解二阶微分方程（边界值）遇到困难，请高手指点
已经有5人回复
- matlab解一阶常微分方程组
已经有4人回复
- 请教用matlab求解一个非线性偏微分方程组的数值解
已经有8人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- 如何理解变系数微分方程
已经有6人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab解微分方程组
已经有15人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
whqs8426212
铜虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 13.2
- 散金: 375
- 红花: 7
- 帖子: 782
- 在线: 685.7小时
- 虫号: 2572504
- 注册: 2013-07-30
- 专业: 数论
|
****************************************************子函数***************************************************** 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 **************************************************************************************************************** |
2楼2014-08-04 00:14:47
