matlab ode45函数 常微分方程中参数为矩阵不是固定值怎么调整代码

chen10217 2017-12-11 12:43:08
新手入门Matlab 实践练习完成一个存在竞争的两个种群增长的模型:

m_1为耐缺氧的细胞种群A的数量
m_2为不耐缺氧细胞种群B的的数量
r_1为A种群的固定增长率,
r_2为B种群的固定增长率
N_x为距离球心x处种群A,B的最大容量(目前设为定值)
σ_1 (O)为B种群相对于A的相对竞争参数,
σ_2 (O)为A种群相对于B的相对竞争参数,它们是含氧量O(s)的函数

我使用ode45解常微分方程,关键问题是O是求解函数后得到的1X3501的矩阵= = 实在改了半天不知道怎么调整。。
代码:

其中的fun函数:

现在出来的图实在是#¥%*&。。。:


直接设为常数画出来的图:


求大神指点!!1
...全文
295 回复 点赞 打赏 收藏 举报
写回复
回复
切换为时间正序
当前发帖距今超过3年,不再开放新的回复
发表回复
相关推荐
发帖
其他开发语言
加入

3386

社区成员

其他开发语言 其他开发语言
申请成为版主
帖子事件
创建了帖子
2017-12-11 12:43
社区公告
暂无公告