24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

liuwul000

铁杆木虫 (著名写手)

应助: 2 (幼儿园)
金币: 6534.9
帖子: 1704
在线: 70.2小时
虫号: 312430

[交流] 【求助】用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) ’这部分用极大值命令代替一下，万分感谢！

回复此楼

» 猜你喜欢

求环氧树脂研发1名已经有10人回复
280求调剂已经有5人回复
什么是人一生最重要的？已经有10人回复
面上可以超过30页吧？已经有13人回复
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？已经有10人回复
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有13人回复
版面费该交吗已经有17人回复
【博士招生】太原理工大学2026化工博士已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

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楼 2010-11-13 21:28:42

已阅关注TA 给TA发消息送TA红花 TA的回帖

wangshaoli2967

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 15.7
帖子: 212
在线: 114.1小时
虫号: 997547

你好！你的问题解决了吗？

★ ★ ★ ★ ★
wuguocheng(金币+5):多谢帮忙 2010-11-23 17:11:48

我这有个用迭代画的，看对你有帮助没有？
close all;
clear all;
c=0;
y=0.0;

hold on
while c < 4,
      for i=1:100;
         y = y.^2 -c; %converge the iteration
      end
      for i=1:20

         y = y.^2 - c;
         plot(c,y,'.'); % plot the converged points
      end
      c=c+0.01;
end