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

返回列表

當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

liuwul000

鐵桿木蟲(chóng) (著名寫(xiě)手)

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

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

我畫(huà)分支圖常常先差分，然后運(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ú)法分辨，有人用極值命令做，好像可以解決這方面的問(wèn)題。
請(qǐng)教一下高手，幫我把上面程序‘plot(b1,Z(T-500:T),'k.','MarkerSize',4) ’這部分用極大值命令代替一下，萬(wàn)分感謝！

回復(fù)此樓

» 猜你喜歡

271材料工程求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專(zhuān)碩調(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)劑材料與化工專(zhuān)碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)

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

Matlab數(shù)值求解二階常微分方程已經(jīng)有9人回復(fù)
如何用數(shù)值解法求二階微分方程的本征值已經(jīng)有18人回復(fù)
matlab求解常微分方程已經(jīng)有6人回復(fù)
matlab處理常微分方程作圖問(wèn)題已經(jīng)有7人回復(fù)
求解兩個(gè)簡(jiǎn)單非線(xiàn)性常微分方程~~ 已經(jīng)有4人回復(fù)
非線(xiàn)性方程組的解法已經(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 線(xiàn)性常微分方程求解已經(jīng)有5人回復(fù)
matlab-常微分方程參數(shù)估計(jì) 已經(jīng)有12人回復(fù)
matlab 求解隱式常微分方程時(shí)用ode15i 函數(shù)出現(xiàn)的問(wèn)題求助已經(jīng)有3人回復(fù)
數(shù)值方法求解偏微分方程的問(wèn)題已經(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的回帖

wangshaoli2967

金蟲(chóng) (小有名氣)

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

你好！你的問(wèn)題解決了嗎？

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

我這有個(gè)用迭代畫(huà)的，看對(duì)你有幫助沒(méi)有？
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的回帖

查看全部 3 個(gè)回答

onesupeng

金蟲(chóng) (職業(yè)作家)

數(shù)學(xué)EPI: 17
應(yīng)助: 256 (大學(xué)生)
貴賓: 1.36
金幣: 2336.2
帖子: 4583
在線(xiàn): 1303.8小時(shí)
蟲(chóng)號(hào): 394701

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

贊一下

回復(fù)此樓

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

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

查看全部 3 個(gè)回答

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

[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4	守候夕陽(yáng)CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 085410人工智能專(zhuān)碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +5	小材化本科 2026-03-18	5/250	2026-03-18 21:53 by 無(wú)懈可擊111
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +8	SUICHILD 2026-03-12	8/400	2026-03-18 20:55 by winsuccess
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	9/450	2026-03-18 15:31 by vgtyfty
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +7	rare12345 2026-03-18	7/350	2026-03-18 14:31 by laoshidan
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 材料專(zhuān)碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考研] 321求調(diào)劑（食品/專(zhuān)碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321