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

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

24小時(shí)熱門版塊排行榜

返回列表

liuwul000

鐵桿木蟲 (著名寫手)

應(yīng)助: 2 (幼兒園)
金幣: 6534.9
帖子: 1704
在線: 70.2小時(shí)
蟲號(hào): 312430

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

我畫分支圖常常先差分，然后運(yùn)行很長(zhǎng)時(shí)間取最后的幾百步就可以了，例如以下
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é)果是周期和混沌無(wú)法分辨，有人用極值命令做，好像可以解決這方面的問題。
請(qǐng)教一下高手，幫我把上面程序‘plot(b1,Z(T-500:T),'k.','MarkerSize',4) ’這部分用極大值命令代替一下，萬(wàn)分感謝！

回復(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ù)
085410人工智能專碩317求調(diào)劑（0854都可以）已經(jīng)有3人回復(fù)

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

Matlab數(shù)值求解二階常微分方程已經(jīng)有9人回復(fù)
如何用數(shù)值解法求二階微分方程的本征值已經(jīng)有18人回復(fù)
matlab求解常微分方程已經(jīng)有6人回復(fù)
matlab處理常微分方程作圖問題已經(jīng)有7人回復(fù)
求解兩個(gè)簡(jiǎn)單非線性常微分方程~~ 已經(jīng)有4人回復(fù)
非線性方程組的解法已經(jīng)有5人回復(fù)
一階常微分方程組，matlab進(jìn)行數(shù)值求解已經(jīng)有3人回復(fù)
微分方程組的穩(wěn)定性分析，急！已經(jīng)有10人回復(fù)
高手進(jìn)來(lái)用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復(fù)
matlab 線性常微分方程求解已經(jīng)有5人回復(fù)
matlab-常微分方程參數(shù)估計(jì) 已經(jīng)有12人回復(fù)
matlab 求解隱式常微分方程時(shí)用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è)作家)

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

贊一下

回復(fù)此樓

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

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

wangshaoli2967

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 15.7
帖子: 212
在線: 114.1小時(shí)
蟲號(hào): 997547

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

★ ★ ★ ★ ★
wuguocheng(金幣+5):多謝幫忙 2010-11-23 17:11:48

我這有個(gè)用迭代畫的，看對(duì)你有幫助沒有？
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ā)表

[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	10/500	2026-03-18 20:41 by Wangjingyue
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6	薛云鵬 2026-03-13	6/300	2026-03-18 14:14 by 脫穎而出
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4	ioodiiij 2026-03-17	4/200	2026-03-18 12:36 by Linda Hu
[考研] 268求調(diào)劑 +7	好運(yùn)連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考博] 26申博 +4	八旬速覽 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 304求調(diào)劑 +5	素年祭語(yǔ) 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月

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

24小時(shí)熱門版塊排行榜

liuwul000

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

» 猜你喜歡

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

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

onesupeng

wangshaoli2967

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

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

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

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

查看全部散金貼

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