24小时热门版块排行榜

>论坛更新日志 (3960)
>文献求助 (408)
>虫友互识 (345)
>导师招生 (323)
>博后之家 (166)
>考博 (155)
>硕博家园 (154)
>休闲灌水 (141)
>论文投稿 (118)
>基金申请 (90)
>招聘信息布告栏 (55)
>考研 (54)
>教师之家 (52)
>公派出国 (43)
>绿色求助(高悬赏) (38)
>找工作 (37)

返回列表

liangbao

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 978.7
帖子: 264
在线: 164.1小时
虫号: 868553

[交流] 【求助】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 ]

回复此楼

» 猜你喜欢

到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有7人回复
申请2026年博士已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有6人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2010-11-25 09:25:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

★
liangbao(金币+2):谢谢参与
liangbao(金币+10): 2010-11-25 20:38:17

每次循环把t，y的值覆盖掉了。
循环里把y赋给一个矩阵的一行即可。

赞一下(1人)

回复此楼

2楼2010-11-25 20:01:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

liangbao

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 978.7
帖子: 264
在线: 164.1小时
虫号: 868553

引用回帖:

Originally posted by xiegangmai at 2010-11-25 20:01:10:
每次循环把t，y的值覆盖掉了。
循环里把y赋给一个矩阵的一行即可。

恩，谢谢
我这样改，运行通过了
您看看，没有什么对不？
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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖