24小时热门版块排行榜    

查看: 493  |  回复: 0

hoddescle

木虫 (小有名气)

老和部队

[求助] matlab计算动力学方程求代码指点

我做毕业设计,是关于陶瓷烧结的。我想试试固相烧结动力学的问题
粗略的方程是
2 surface=1 interface+x hole
我想用ode指令求解surface和interface关于时间的变化情况。
当反应达到平衡,烧结结束
比如surfer的初值是2,其他的都是0
ds/dt=-k1*s^2
di/dt=k2*i*h^x
dh/dt=k3*i*h^x
为简化起见x看做定值
但是在设定的时间范围内每次都会出现要么s一直减小到负数,而其他的都会增加的很多很多。我想求助解这个方程组的matlab代码,然后能画出surface和interface关于时间的变化情况。
关键问题是怎么让电脑明白s减小之后变成了i和h
i和h减小后增加到了s上

[ Last edited by hoddescle on 2012-2-18 at 21:01 ]
回复此楼
忠贤随明主,奸佞觅昏君
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 hoddescle 的主题更新
信息提示
请填处理意见