| 查看: 2187 | 回复: 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) ’这部分用极大值命令代替一下，万分感谢！ |
回复此楼» 猜你喜欢
听说今天filecode变了
已经有33人回复
-
基金中了
已经有6人回复
-
filecode
已经有14人回复
-
固定端突然变了，今天
已经有7人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
» 抢金币啦！回帖就可以得到:
-
山东征女友，坐标济南
+1/163
-
1,2-双(三甲基硅氧基)乙烷（CAS：7381-30-8 ）
+1/87
-
3,4-二(三甲基硅氧基)-1-丁烯（CAS：89543-84-0）
+1/87
-
2-氯-1,3,2-二氧磷杂环戊烷（CAS：822-39-9）
+1/86
-
1,2-双(三甲基硅氧基)苯（CAS：5075-52-5）
+1/85
-
三(三甲基硅烷)磷酸酯（CAS：10497-05-9）
+1/85
-
招聘机械领域专利代理人
+1/81
-
具身智能机器人公司招聘机械结构岗位（工程师、专家、负责人均可）薪酬最高150万
+1/80
-
西安交大马伟/王亚洲课题组诚聘科研助理/2027入学博士-有机生物电子学方向
+5/55
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/44
-
2027 Chemical Engineering Ph.D 招生
+1/40
-
同济大学环境学院肖倩课题组2026年补充招收推免生
+1/31
-
北理工集成电路杰青团队 | 诚招科助理
+1/19
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/18
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
+1/12
-
《上海高薪》VS《中科院》
+1/11
-
希望大家都高中！！
+1/7
-
北理工集成电路杰青团队 | 诚招科助理
+1/5
-
国家纳米科学中心段鹏飞课题组招收2027年入学博士研究生
+1/4
-
高弯曲强度玻璃制品的特殊加工工艺
+1/4
2楼2010-11-14 02:09:12
3楼2010-11-21 14:25:57
