| 查看: 998 | 回复: 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 ] |
回复此楼» 猜你喜欢
271材料工程求调剂
已经有5人回复
-
281求调剂（0805）
已经有16人回复
-
304求调剂
已经有6人回复
-
材料工程专硕调剂
已经有6人回复
-
一志愿天大材料与化工（085600）总分338
已经有4人回复
-
085700资源与环境308求调剂
已经有3人回复
-
求材料调剂
已经有8人回复
-
294求调剂材料与化工专硕
已经有5人回复
-
一志愿华中科技大学，080502，354分求调剂
已经有4人回复
-
一志愿吉林大学材料学硕321求调剂
已经有6人回复
» 抢金币啦！回帖就可以得到:
-
祈福2026国基面上项目
+5/350
-
化学化工学院 招收化工、化学、材料等相关方向研究生（学硕、专硕都有调剂名额）
+1/89
-
天津理工大学国家杰青王铁课题组招收2026年博士研究生
+1/85
-
西北大学化学与材料科学学院博士招生（还有两个名额，3月30日截至）
+1/68
-
青岛高校 研究生招生
+1/55
-
山东师范大学有机化学专业胡忠燕老师课题组招收2026届硕士研究生以及调剂生
+1/40
-
华南师范大学（211）博士招生- 电子、自动化、机械、生物学、物理相关专业
+2/38
-
西北工业大学光子学研究团队诚聘博士后！
+2/34
-
澳大利亚西澳大学（UWA）张金强课题组招收全奖博士生（化工新能源方向）
+1/32
-
26年博士自荐
+1/17
-
【预硕】【预博】青岛大学生物医药创新转化
+1/16
-
Call for paper: [Fermentation – IF 3.3]
+2/14
-
新加坡国立大学药学系化学生物学课题组招PhD
+1/13
-
南京医科大学-膜蛋白结构与功能课题组-招收2026级博士研究生
+1/11
-
新加坡国立大学药学系化学生物学课题组招PhD
+1/8
-
课题组招收环境及相关专业调剂硕士研究生（欢迎优秀学生加入）
+1/6
-
东北林业大学（211，双一流）宋老师课题组招收高分子、化学、复合材料等方向研究生
+1/4
-
南方科技大学生命科学学院基础免疫与微生物学系招收博士后
+1/3
-
巴塞罗那自治大学（UAB）-巴塞罗那材料研究所（ICMAB）招收博士研究生
+1/2
-
福建师范大学环境微生物技术课题组接受2026级硕士调剂
+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
