| 查看: 2076 | 回复: 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) ’这部分用极大值命令代替一下,万分感谢! |
» 猜你喜欢
生物学308求调剂(一志愿华东师大)
已经有6人回复
调剂求收留
已经有7人回复
272分材料子求调剂
已经有36人回复
275求调剂
已经有8人回复
一志愿211,化学学硕,310分,本科重点双非,求调剂
已经有20人回复
070300化学学硕311分求调剂
已经有19人回复
材料与化工调剂
已经有13人回复
材料与化工调剂
已经有33人回复
复试调剂
已经有7人回复
一志愿哈工大 085600 277 12材科基求调剂
已经有17人回复
» 本主题相关价值贴推荐,对您同样有帮助:
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人回复
» 抢金币啦!回帖就可以得到:
长江大学化工学院2026年硕士研究生调剂公告(理、工均收)
+1/90
江西科技师范大学物理化学课题组招生,07,08都可,还有名额
+1/89
华北理工大学+矿业工程+招收调剂13名【调剂系统已开】
+1/89
大湾区大学李红庚与中国科学技术大学王柳教授联合招收博士后
+1/87
大湾区大学李红庚与中国科学技术大学王柳教授联合招收博士后
+1/86
坐标上海和长三角 诚征女友
+1/74
接收调剂:组内氛围轻松活跃,欢迎加入
+1/44
化学、化学工程、生物与医药,微生物学、生物化学与分子生物学专业调剂
+1/43
调剂系统即将关闭-广西大学农学院硕士研究生调剂(报考专业为09开头)
+2/40
2026-大连大学~六盘水师范学院联合培养硕士研究生,招收化学、环境等相关专业
+1/18
太原科技大学光学工程学科接受调剂
+1/11
深圳理工大学博士后招聘-先进拉曼技术课题组
+1/10
中山大学课题组招实习生有机合成经验优先
+1/7
Top-88悉尼科技大学数据科学/AI 招收2027年入学 校奖 博士生1到2名(国际和本地学生)
+1/5
广东省五邑大学MOFs催化功能材料课题组招收化学/材料背景调剂生
+1/5
烟台大学精准材料高等研究院26年材料方向研究生招生
+1/5
福建师范大学环境微生物技术课题组接受2026级硕士调剂
+1/5
招收化学/材料调剂同学
+1/4
中山大学课题组招实习生本硕皆可,有细胞培养,分子生物,动物实验经验优先
+1/2
安徽工程大学化学与环境工程学院2026级硕士研究生还有部分调剂名额-资源与环境专硕
+1/1
2楼2010-11-14 02:09:12
3楼2010-11-21 14:25:57













回复此楼