 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1001 | 回復(fù): 3 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
【求助】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ù)此樓» 猜你喜歡
南京大學(xué)化學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑!
已經(jīng)有5人回復(fù)
-
283求調(diào)劑（080500）
已經(jīng)有7人回復(fù)
-
311求調(diào)劑
已經(jīng)有6人回復(fù)
-
275求調(diào)劑
已經(jīng)有20人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有9人回復(fù)
-
294分080500材料科學(xué)與工程求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿北京工業(yè)大學(xué)，324分求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
已經(jīng)有11人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所博士后招聘
+1/280
-
山東建筑大學(xué) （市政與環(huán)境工程學(xué)院）碩士研究生招生啦 靜電紡絲纖維與鹽水分離應(yīng)用
+1/86
-
江南大學(xué)招收具有3年海外博后經(jīng)歷的青年教師
+1/80
-
重慶交通大學(xué)院士團(tuán)隊(duì)招收材料、能源動(dòng)力專業(yè)碩士研究生
+2/48
-
南方油�！Y源與環(huán)境專碩招生，化學(xué)環(huán)境生物類專業(yè)均可，就業(yè)好，有專業(yè)壁壘
+1/47
-
大連工業(yè)大學(xué) 高分子材料 接收考研調(diào)劑
+1/43
-
考研調(diào)劑!
+1/42
-
2026延安大學(xué)設(shè)施蔬菜環(huán)境調(diào)控與分子機(jī)制研究團(tuán)隊(duì)團(tuán)隊(duì)招收相關(guān)專業(yè)調(diào)劑研究生
+1/39
-
招收2026級(jí)環(huán)境科學(xué)方向（專業(yè)代碼083001）碩士研究生1-2名
+1/37
-
濟(jì)南大學(xué)水利與環(huán)境學(xué)院趙艷俠教授簡(jiǎn)介
+1/27
-
食品學(xué)碩考研調(diào)劑
+1/11
-
中科院和北京工商大學(xué)招收2026博士（化學(xué)背景）
+1/11
-
26年申博自薦-計(jì)算機(jī)視覺
+1/9
-
長(zhǎng)江師范學(xué)院化學(xué)工程方向接收大量調(diào)劑考生
+1/9
-
蘇州初創(chuàng)公司誠(chéng)聘有機(jī)合成化學(xué)家
+1/7
-
22408一志愿211，329分，本科績(jī)優(yōu)，計(jì)算機(jī)科班
+1/6
-
香港中文大學(xué)（深圳）生物智能制造方向招生-高分子功能材料、智能/仿生材料、生物材料
+1/5
-
武漢紡織大學(xué)沈愛國(guó)教授團(tuán)隊(duì)招收化學(xué)、生物、材料專業(yè)研究生
+1/4
-
三峽大學(xué)材化學(xué)院磷硅精細(xì)化工課題組招收化學(xué)類、材料類、材料與化工類碩士
+1/3
-
武漢理工大學(xué)船海與能源動(dòng)力工程學(xué)院劉鵬飛院士團(tuán)隊(duì)博士后招聘啟事
+1/2
|
恩，謝謝 我這樣改，運(yùn)行通過了 您看看，沒有什么對(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
2樓2010-11-25 20:01:10
xiegangmai
版主 (職業(yè)作家)
- 仿真EPI: 3
- 應(yīng)助: 157 (高中生)
- 貴賓: 5.735
- 金幣: 40609
- 帖子: 4046
- 在線: 1330.2小時(shí)
- 蟲號(hào): 653607
4樓2010-11-26 08:34:13
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +8 | 景隨風(fēng) 2026-03-29 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
|
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +4 | 12.27 2026-03-27 | 4/200 |
|
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[考研] 085600，材料與化工321分求調(diào)劑 +9 | 大饞小子 2026-03-28 | 9/450 |
|
|
|
[考研]
 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑
+10
|
y7czhao 2026-03-26 | 10/500 |
|
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么！ 2026-03-26 | 7/350 |
|
|
|
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
|
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|