| 查看: 2220 | 回复: 2 | |||
[交流]
【求助】用Matlab画常微分方程的分支图
|
|
我画分支图常常先差分，然后运行很长时间取最后的几百步就可以了，例如以下 clear; T=1000;dt=0.05;a1=5.0;a2=0.1;b2=2.0;d1=0.4;d2=0.01;e=1.0;r=0; X=zeros(1,T);Y=zeros(1,T);Z=zeros(1,T); X(1)=0.4;Y(1)=0.6;Z(1)=5; hold on for b1=2:0.02:9 for t=2:T X(t)=X(t-1)+(X(t-1)*(1-X(t-1))-a1*X(t-1)*Y(t-1)/(1+b1*X(t-1)))*dt; Y(t)=Y(t-1)+(a1*X(t-1)*Y(t-1)/(1+b1*X(t-1))-a2*Y(t-1)^(e)*Z(t-1)/(1+b2*Y(t-1)^(e)+r*Z(t-1))-d1*Y(t-1))*dt; Z(t)=Z(t-1)+(a2*Y(t-1)^(e)*Z(t-1)/(1+b2*Y(t-1)^(e)+r*Z(t-1))-d2*Z(t-1))*dt; end plot(b1,Z(T-500:T),'k.','MarkerSize',4) end 但这样做的结果是周期和混沌无法分辨，有人用极值命令做，好像可以解决这方面的问题。 请教一下高手，帮我把上面程序‘plot(b1,Z(T-500:T),'k.','MarkerSize',4) ’这部分用极大值命令代替一下，万分感谢！ |
回复此楼» 猜你喜欢
怎么查啊
已经有6人回复
-
基金系统什么内容也没有
已经有8人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有10人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然面上复盘~欢迎讨论
已经有12人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab数值求解二阶常微分方程
已经有9人回复
- 如何用 数值解法 求 二阶微分方程的本征值
已经有18人回复
- matlab求解常微分方程
已经有6人回复
- matlab处理常微分方程作图问题
已经有7人回复
- 求解两个简单非线性常微分方程~~
已经有4人回复
- 非线性方程组的解法
已经有5人回复
- 一阶常微分方程组，matlab进行数值求解
已经有3人回复
- 微分方程组的稳定性分析，急！
已经有10人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab 线性常微分方程求解
已经有5人回复
- matlab-常微分方程参数估计
已经有12人回复
- matlab 求解隐式常微分方程时用ode15i 函数出现的问题求助
已经有3人回复
- 数值方法求解偏微分方程的问题
已经有4人回复
- 常微分方程(ODEs)的MATLAB数值解法
已经有119人回复
- 【求助】matlab解常微分方程，暂停和继续
已经有5人回复
2楼2010-11-14 02:09:12
3楼2010-11-21 14:25:57
10