版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

liuwul000

鐵桿木蟲 (著名寫手)

應(yīng)助: 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
但這樣做的結(jié)果是周期和混沌無法分辨，有人用極值命令做，好像可以解決這方面的問題。
請教一下高手，幫我把上面程序‘plot(b1,Z(T-500:T),'k.','MarkerSize',4) ’這部分用極大值命令代替一下，萬分感謝！

回復(fù)此樓

» 猜你喜歡

271材料工程求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

Matlab數(shù)值求解二階常微分方程已經(jīng)有9人回復(fù)
如何用數(shù)值解法求二階微分方程的本征值已經(jīng)有18人回復(fù)
matlab求解常微分方程已經(jīng)有6人回復(fù)
matlab處理常微分方程作圖問題已經(jīng)有7人回復(fù)
求解兩個簡單非線性常微分方程~~ 已經(jīng)有4人回復(fù)
非線性方程組的解法已經(jīng)有5人回復(fù)
一階常微分方程組，matlab進(jìn)行數(shù)值求解已經(jīng)有3人回復(fù)
微分方程組的穩(wěn)定性分析，急！已經(jīng)有10人回復(fù)
高手進(jìn)來用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復(fù)
matlab 線性常微分方程求解已經(jīng)有5人回復(fù)
matlab-常微分方程參數(shù)估計已經(jīng)有12人回復(fù)
matlab 求解隱式常微分方程時用ode15i 函數(shù)出現(xiàn)的問題求助已經(jīng)有3人回復(fù)
數(shù)值方法求解偏微分方程的問題已經(jīng)有4人回復(fù)
常微分方程(ODEs)的MATLAB數(shù)值解法已經(jīng)有119人回復(fù)
【求助】matlab解常微分方程，暫停和繼續(xù) 已經(jīng)有5人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2010-11-13 21:28:42

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

onesupeng

金蟲 (職業(yè)作家)

你這個精度不是太好，換一個積分格式

贊一下

回復(fù)此樓

2樓2010-11-14 02:09:12

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangshaoli2967

金蟲 (小有名氣)

應(yīng)助: 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

贊一下(1人)

回復(fù)此樓

3樓2010-11-21 14:25:57

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 liuwul000 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿武理材料305分求調(diào)劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無際的草原
[考研] 085600材料與化工 +5	安全上岸！ 2026-03-16	5/250	2026-03-18 15:33 by cmz0325
[教師之家] 焦慮 +8	水冰月月野兔 2026-03-13	12/600	2026-03-18 15:27 by 咪嗚喵嗚
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 299求調(diào)劑 +5	△小透明* 2026-03-17	5/250	2026-03-18 11:49 by 盡舜堯1
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[基金申請] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

liuwul000

[交流] 【求助】用Matlab畫常微分方程的分支圖

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

onesupeng

wangshaoli2967

你好！你的問題解決了嗎？

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

你好！你的問題解決了嗎？