版塊導(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í)熱門版塊排行榜

返回列表

liangbao

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 978.7
帖子: 264
在線: 164.1小時(shí)
蟲號(hào): 868553

[交流] 【求助】ode45數(shù)值求解常微分方程總是錯(cuò)，找不出原因

程序如下：S1=0.4;
S2=0.8;
y=zeros(42,3);
y(1,1)=6e10;
y(2,1)=4e10;
y(3,1)=0;
for i=0:5;
for j=1:7
      if(j==6)|(j==7)
      [t,y]=ode45('tumor',[i*7+j  i*7+j+1],[y(i*7+j,1) y(i*7+j,2) y(i*7+j,3)]);
      else
         a=y(i*7+j,1);
         y(i*7+j,1)=S1*a;
         b=y(i*7+j,2);
         y(i*7+j,2)=S2*b;
         y(i*7+j,3)=(1-S1)*a+(1-S2)*b;
      [t,y]=ode45('tumor',[i*7+j i*7+j+1],[y(i*7+j,1) y(i*7+j,2) y(i*7+j,3)]);
      end
end
end
plot(t,y(1,

,'*',t,y(2,

,'.',t,y(3,

,'+');
錯(cuò)誤為：Index exceeds matrix dimensions.

Error in ==> number1 at 12
[t,y]=ode45('tumor',[i*7+j i*7+j+1],[y(i*7+j,1) y(i*7+j,2) y(i*7+j,3)]);
我估計(jì)是向量傳遞出錯(cuò)了
發(fā)現(xiàn)ode45返回的是t和y(:,1) y(:,2) y(:,3)的列向量，
如果想得到向量的最后一個(gè)值應(yīng)該怎么辦？

[ Last edited by liangbao on 2010-11-25 at 09:51 ]

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有9人回復(fù)
293求調(diào)劑已經(jīng)有14人回復(fù)
0703化學(xué)336分求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有9人回復(fù)
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請(qǐng) 已經(jīng)有3人回復(fù)

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

查看全部散金貼

1樓 2010-11-25 09:25:06

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

xiegangmai

版主 (職業(yè)作家)

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
帖子: 4046
在線: 1330.2小時(shí)
蟲號(hào): 653607

★
liangbao(金幣+2):謝謝參與
liangbao(金幣+10): 2010-11-25 20:38:17

每次循環(huán)把t，y的值覆蓋掉了。
循環(huán)里把y賦給一個(gè)矩陣的一行即可。

贊一下(1人)

回復(fù)此樓

2樓2010-11-25 20:01:10

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

liangbao

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 978.7
帖子: 264
在線: 164.1小時(shí)
蟲號(hào): 868553

引用回帖:

Originally posted by xiegangmai at 2010-11-25 20:01:10:
每次循環(huán)把t，y的值覆蓋掉了。
循環(huán)里把y賦給一個(gè)矩陣的一行即可。

恩，謝謝
我這樣改，運(yùn)行通過(guò)了
您看看，沒(méi)有什么對(duì)不？
S1=0.4;
S2=0.8;
y=zeros(100,42);
y(1,1)=6e10;
y(1,2)=4e10;
y(1,3)=0;
for i=0:5;
for j=1:7
if(j==6)|(j==7)
tspan=i*7+j:0.05:i*7+j+1;
[t,p]=ode45('tumor',tspan,[y(i*7+j,1) y(i*7+j,2) y(i*7+j,3)]);
else
      y(i*7+j,1)=S1*y(i*7+j,1);
      y(i*7+j,2)=S2*y(i*7+j,2);
      y(i*7+j,3)=((1-S1)/S1)*y(i*7+j,1)+((1-S2)/S2)*y(i*7+j,2);
tspan=i*7+j:0.05:i*7+j+1;
[t,p]=ode45('tumor',tspan,[y(i*7+j,1) y(i*7+j,2) y(i*7+j,3)]);
end
plot(t,p(:,1),'*',t,p(:,2),'.',t,p(:,3));
hold on;
y(i*7+j+1,1)=p(21,1);
y(i*7+j+1,2)=p(21,2);
y(i*7+j+1,3)=p(21,3);
end
end

贊一下

回復(fù)此樓

3樓2010-11-25 20:38:04

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

xiegangmai

版主 (職業(yè)作家)

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
帖子: 4046
在線: 1330.2小時(shí)
蟲號(hào): 653607

★
liangbao(金幣+2):謝謝參與

沒(méi)'tumor'函數(shù)，不好判斷運(yùn)行結(jié)果正確不正確。

贊一下(1人)

回復(fù)此樓

4樓2010-11-26 08:34:13

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

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

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 268求調(diào)劑 +3	簡(jiǎn)單點(diǎn)0 2026-03-17	3/150	2026-03-17 23:15 by wangjihu
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 341求調(diào)劑 +5	搗蛋豬豬 2026-03-11	7/350	2026-03-17 19:09 by 搗蛋豬豬
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287 +3	小李0854 2026-03-16	3/150	2026-03-17 13:40 by 熱情沙漠
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 311求調(diào)劑 +5	26研0 2026-03-15	5/250	2026-03-16 16:21 by a不易
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +5	ZY-05 2026-03-13	6/300	2026-03-14 10:51 by Jy?
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5	瑪卡巴卡啊哈 2026-03-11	5/250	2026-03-13 21:40 by JourneyLucky
[考研] 315求調(diào)劑 +9	小羊小羊_ 2026-03-11	10/500	2026-03-13 21:13 by SXNU李老師
[考研] 求調(diào)劑 +7	18880831720 2026-03-11	7/350	2026-03-13 16:10 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過(guò) 本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz

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

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

liangbao

[交流] 【求助】ode45數(shù)值求解常微分方程總是錯(cuò)，找不出原因

» 猜你喜歡

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

xiegangmai

liangbao

xiegangmai

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

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

查看全部散金貼