| 查看: 2188 | 回复: 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
已经有18人回复
听说今天filecode变了
已经有44人回复
基金中了
已经有9人回复
固定端突然变了,今天
已经有8人回复
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人回复
» 抢金币啦!回帖就可以得到:
一边期待,一边忐忑:发榜前的科研人正在想什么?
+5/175
2026年量子算法、密码学与数据分析国际会议(QACDA 2026)
+1/93
3-羟基丙腈(CAS:109-78-4)
+1/88
1,3,2-Dioxaphospholane, 2-(2,2,3,3,3-pentafluoropropoxy)-
+1/87
三氟甲基亚磺酸锂
+1/87
丁烯二醇(CAS: 497-06-3)
+1/86
香港城市大学郑星课题组诚聘博士后 (生态学、环境科学、植物学、大气科学背景)
+1/82
化工学院邱惠斌教授联合同仁医院招聘博士后
+1/65
西安交大马伟/王亚洲课题组诚聘科研助理/2027入学博士-有机生物电子学方向
+5/40
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/34
双一流高校招收博士(现代技术管理/环境技术经济管理)及硕士研究生(长期有效)
+1/14
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/9
福建师范大学海峡柔性电子学院-2027年9月入学博士-考核制-电催、能源
+1/7
有合作举办scopus会议的吗?
+1/6
香港理工大学殷骏课题组诚聘1 名博士后以及2 名博士生
+1/6
POSTDOCTORAL POSITION in Structural Virology at the University of Minnesota
+1/3
北京师范大学(珠海校区)赵亚然副教授招收2027年电催化推免硕士/博士研究生
+1/2
北理工集成电路杰青团队 | 诚招科助理
+1/2
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/2
国自然发榜在即:此刻的你,是什么心情?
+1/1
2楼2010-11-14 02:09:12
3楼2010-11-21 14:25:57










回复此楼