 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 994 | 回復(fù): 3 | |||
[交流]
【求助】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)有8人回復(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ù)
-
268求調(diào)劑
已經(jīng)有3人回復(fù)
-
311求調(diào)劑
已經(jīng)有10人回復(fù)
-
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開(kāi)始管格式了
已經(jīng)有4人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
山東征女友，坐標(biāo)濟(jì)南
+1/167
-
福建師范大學(xué)招收2026年化學(xué)、材料碩士3-4名
+1/81
-
愿遇見(jiàn)你，遇見(jiàn)不止擦肩而過(guò)的緣分
+1/76
-
華南師范大學(xué)（211）博士招生- 電子、自動(dòng)化、機(jī)械、生物學(xué)、物理相關(guān)專業(yè)
+2/38
-
歡迎交流咨詢
+1/28
-
中山大學(xué)院士團(tuán)隊(duì)微納器件、腦機(jī)接口方向博士招生、博后招聘
+1/21
-
26年博士自薦
+1/17
-
【預(yù)碩】【預(yù)博】青島大學(xué)生物醫(yī)藥創(chuàng)新轉(zhuǎn)化
+1/16
-
北京某研究院生物化學(xué)與分子生物學(xué)相關(guān)專業(yè)申請(qǐng)考核制博士招生（春季）
+1/15
-
2026年西南科技大學(xué)功能涂層課題組簡(jiǎn)介
+1/12
-
【博士招生】武漢科技大學(xué)招收材料、化工、環(huán)境、冶金類2026年“申請(qǐng)-考核”制博士生
+1/8
-
英國(guó)�？巳卮髮W(xué) & 法國(guó)巴黎薩克雷大學(xué)聯(lián)合培養(yǎng)博士
+1/7
-
華中科技大學(xué)管理學(xué)院招聘社會(huì)用工（科研助理）1名
+1/6
-
北京高校副校長(zhǎng)團(tuán)隊(duì)招收機(jī)械類，環(huán)境類學(xué)碩和專碩
+1/6
-
上海第二工業(yè)大學(xué)-朱大海課題組招生（過(guò)線就能調(diào)劑）
+1/5
-
雙一流天津工業(yè)大學(xué)電信學(xué)院李鴻強(qiáng)教授招收2026年申請(qǐng)審核制博士3名
+1/5
-
南方科技大學(xué)生命科學(xué)學(xué)院基礎(chǔ)免疫與微生物學(xué)系招收博士后
+1/4
-
Top-88悉尼科技大學(xué)數(shù)據(jù)科學(xué)/AI 招收2027年入學(xué) 校獎(jiǎng) 博士生1到2名(國(guó)際和本地學(xué)生）
+1/4
-
中國(guó)科學(xué)院地球化學(xué)研究所招科研助理
+1/1
-
北京大學(xué)環(huán)境學(xué)院張子帥研究員招收科研助理
+1/1
xiegangmai
版主 (職業(yè)作家)
- 仿真EPI: 3
- 應(yīng)助: 157 (高中生)
- 貴賓: 5.735
- 金幣: 40609
- 帖子: 4046
- 在線: 1330.2小時(shí)
- 蟲號(hào): 653607
2樓2010-11-25 20:01:10
|
恩，謝謝 我這樣改，運(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 |
3樓2010-11-25 20:38:04
xiegangmai
版主 (職業(yè)作家)
- 仿真EPI: 3
- 應(yīng)助: 157 (高中生)
- 貴賓: 5.735
- 金幣: 40609
- 帖子: 4046
- 在線: 1330.2小時(shí)
- 蟲號(hào): 653607
4樓2010-11-26 08:34:13
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 liangbao 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +7 | yangfz 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +4 | share_joy 2026-03-16 | 4/200 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 312求調(diào)劑 +4 | 陌宸希 2026-03-16 | 5/250 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287 +3 | 小李0854 2026-03-16 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | xc321 2026-03-12 | 6/300 |
|