| 查看: 999 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】ode45数值求解常微分方程总是错，找不出原因
|
|||
|
程序如下：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)]); 我估计是向量传递出错了 发现ode45返回的是t和y(:,1) y(:,2) y(:,3)的列向量， 如果想得到向量的最后一个值应该怎么办？ [ Last edited by liangbao on 2010-11-25 at 09:51 ] |
回复此楼» 猜你喜欢
材料专硕英一数二306
已经有5人回复
-
085600材料与化工调剂 324分
已经有9人回复
-
0703化学调剂
已经有12人回复
-
0703化学 305求调剂
已经有4人回复
-
0703化学调剂，求各位老师收留
已经有10人回复
-
271材料工程求调剂
已经有5人回复
-
281求调剂（0805）
已经有16人回复
-
304求调剂
已经有6人回复
-
材料工程专硕调剂
已经有6人回复
-
一志愿天大材料与化工（085600）总分338
已经有4人回复
» 抢金币啦！回帖就可以得到:
-
青岛高校2026招生硕士研究生
+5/840
-
深圳技术大学集成电路与光电芯片学院郝俊杰课题组诚聘博士后、研究助理、访问学生
+2/120
-
【技术服务】研究所横向技术开发及服务
+1/99
-
散金币，求好运，祝面上顺利！
+2/92
-
化学化工学院 招收化工、化学、材料等相关方向研究生（学硕、专硕都有调剂名额）
+1/88
-
愿遇见你，遇见不止擦肩而过的缘分
+1/75
-
中国药科大学211-课题组招收2026博士/硕士
+1/55
-
PBAT溶解做成乳液
+1/40
-
上海交通大学沈道智副教授招收2026年联培博士/硕士研究生（微纳器件与制造方向）
+1/33
-
天津工业大学航空航天博士招生
+1/32
-
华南师范大学（211）申请审核制博士招生- 电子、自动化、机械、生物学、物理相关专业
+1/25
-
【武汉高校】【新能源材料与器件课题组】 招收优秀硕士调剂生！！！
+5/15
-
南开大学罗景山教授课题组招收2026年入学博士生（光/电催化方向）
+2/14
-
新加坡国立大学药学系化学生物学课题组招PhD
+1/8
-
聊城大学化院青椒招收化学相关专业调剂学生
+1/7
-
125603工业工程与管理调剂
+1/2
-
国家杰青低维材料与器件力学团队2026年招收博士研究生
+1/2
-
339求调剂
+1/1
-
福建师范大学环境微生物技术课题组接受2026级硕士调剂
+1/1
-
Top-88悉尼科技大学数据科学/AI 招收2027年入学 校奖 博士生1到2名(国际和本地学生）
+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
