应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 应用数学与力学 » 【求助】用Matlab画常微分方程的分支图
51/1返回列表
查看: 2078  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

liuwul000

铁杆木虫 (著名写手)

[交流] 【求助】用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楼 2010-11-13 21:28:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wangshaoli2967

金虫 (小有名气)

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

★ ★ ★ ★ ★
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
赞一下(1人) 回复此楼
3楼2010-11-21 14:25:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

onesupeng

金虫 (职业作家)
你这个精度不是太好，换一个积分格式
赞一下 回复此楼
2楼2010-11-14 02:09:12
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人女孩 +5 吼吼， 2026-04-10 6/300 2026-04-11 00:09 by 任文亮
[考研] 一志愿211，化学310分，本科重点双非，求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 材料085601调剂 +25 何润采123 2026-04-10 27/1350 2026-04-10 23:17 by Ftglcn90
[考研] 080500求调剂 +16 黄宇博 2026-04-06 16/800 2026-04-10 22:33 by Ftglcn90
[考研] 0856专硕求调剂 希望是a区院校 +21 好好休息好不好 2026-04-09 24/1200 2026-04-10 16:58 by luoyongfeng
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +20 慕绝cc 2026-04-09 24/1200 2026-04-10 10:22 by xujun0624
[考研] 293求调剂 +5 勇远库爱314 2026-04-08 5/250 2026-04-10 08:46 by vgtyfty
[考研] 314求调剂 +14 weltZeng 2026-04-09 14/700 2026-04-09 23:14 by wolf97
[考研] 291分调剂 +3 上岸小莹加油 2026-04-09 4/200 2026-04-09 17:28 by 1753564080
[考研] 305求调剂 +4 77Qi 2026-04-07 4/200 2026-04-09 17:27 by wp06
[考研] 301求调剂 +18 121. 2026-04-04 18/900 2026-04-07 17:49 by 蓝云思雨
[考研] 一志愿北京化工085600 310分求调剂 +20 0856材料与化工3 2026-04-04 22/1100 2026-04-07 15:14 by 上岸快快
[考研] 一志愿华中农业大学0710（A）初试329分 求调剂 +5 一名26考研生 2026-04-04 5/250 2026-04-07 08:54 by 18828373951
[考研] 机械专硕274求调剂，不挑专业学校 +6 泛泛2333 2026-04-05 8/400 2026-04-06 18:06 by 泛泛2333
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂 +12 万事宜臻 2026-04-04 12/600 2026-04-06 07:46 by 无际的草原
[考研] 考研调剂 +5 美丽的youth_ 2026-04-04 6/300 2026-04-06 06:57 by houyaoxu
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
信息提示
关闭
请填处理意见
关闭