 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1002 | 回復: 3 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[交流]
【求助】ode45數(shù)值求解常微分方程總是錯，找不出原因
|
|||
|
程序如下：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,,'+');錯誤為：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)]); 我估計是向量傳遞出錯了 發(fā)現(xiàn)ode45返回的是t和y(:,1) y(:,2) y(:,3)的列向量， 如果想得到向量的最后一個值應該怎么辦？ [ Last edited by liangbao on 2010-11-25 at 09:51 ] |
回復此樓» 猜你喜歡
070300化學354求調劑
已經有7人回復
-
286分調劑
已經有6人回復
-
287求調劑
已經有11人回復
-
求調劑
已經有7人回復
-
一志愿鄭大材料工程290求調劑
已經有4人回復
-
材料科學與工程調劑
已經有4人回復
-
化學308分調劑
已經有8人回復
-
300求調劑，材料科學英一數(shù)二
已經有16人回復
-
085701求調劑初試286分
已經有5人回復
-
289求調劑
已經有14人回復
» 搶金幣啦！回帖就可以得到:
-
中國石油大學（北京） 化學工程專業(yè)招2026級博士
+1/487
-
耐熱氧化耐熱黃變效果好的聚醚多元醇制備的聚氨酯都有哪些潛在應用領域？
+2/136
-
下一站，欽州
+1/78
-
武漢紡織大學全國重點實驗室招收材料、化學、紡織相關研究生
+1/73
-
青島農業(yè)大學海洋學院招調劑碩士研究生
+1/72
-
青島大學固態(tài)電池研究團隊招收材料、化學、物理等方向博士研究生
+2/44
-
大連工業(yè)大學 高分子材料 接收考研調劑
+1/41
-
中山大學地理科學與規(guī)劃學院退化土壤環(huán)境修復課題組招聘博士后及科研助理
+1/40
-
五邑大學環(huán)境與化學工程學院先進有機高分子材料團隊
+1/39
-
南京信息工程大學 化材院接受調劑研究生
+1/37
-
中國石油大學（北京）國家級大人才團隊博士招生1名：化學、材料、石油工程：油田化學
+1/36
-
南通大學國家級人才團隊招收衛(wèi)生檢驗、醫(yī)學檢驗相關方向調劑研究生
+2/36
-
澳大利亞昆士蘭大學化工學院招收全獎博士生2名
+1/22
-
汕頭大學海洋與生命學科李升康教授團隊招收生物類調劑生
+1/18
-
哈爾濱工業(yè)大學-大灣區(qū)大學2026年聯(lián)合培養(yǎng)博士研究生春季招生-報名3月31號截至
+1/12
-
南京林業(yè)大學化學工程學院柏惺峰課題組招收碩士研究生1名
+1/5
-
招收從事大氣復合污染成因領域研究生
+1/3
-
調劑來這邊
+1/2
-
【MP期刊】青年編委招募開啟！加入我們，一起推動材料科學前沿
+1/1
-
內蒙古科技大學鐵電新能源材料與器件團隊簡介
+1/1
4樓2010-11-26 08:34:13
2樓2010-11-25 20:01:10
|
恩，謝謝 我這樣改，運行通過了 您看看，沒有什么對不？ 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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調劑 +11 | 芋泥小鈴鐺 2026-03-28 | 11/550 |
|
|---|---|---|---|---|
|
|
[考研] 求調劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
|
[考研] 085600 材料與化工 329分求調劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|
[考研] 一志愿鄭州大學，080500學碩，總分317分求調劑 +8 | 舉個栗子oi 2026-03-24 | 9/450 |
|
|
|
[考研] 340求調劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調劑，任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
|
[考研] 332求調劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
|
[考研] 復試調劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 292求調劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[有機交流]
高溫高壓反應求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
|
[考研] 一志愿上海理工能源動力（085800）310分求調劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 266求調劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 272求調劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 材料考研求調劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
|
[考研]
 085602 289分求調劑
+8
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
|
[考研] 309求調劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 07化學303求調劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[考研] 考研一志愿蘇州大學初始315（英一）求調劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
|
[考研] 318求調劑 +3 | plum李子 2026-03-23 | 3/150 |
|