| 查看: 715 | 回复: 11 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
monitor2885至尊木虫 (知名作家)
队长
|
[交流]
【求助】u顺次取值
|
||
|
function test global u t=0:1:10; [t,y]=ode45(@ivpodefun,t,[0 2]); plot(t,y(:,1),'ro-',t,y(:,2),'b^-') function dydt=ivpodefun(t,y,u) global u u=xlsread('data.xls');% u在excel里面是一列数，1；2；3；4；5；6；5；4；3；2；1； dydt=zeros(2,1); dydt(1)=y(2); dydt(2)=u*(1-y(1)^2)*y(2)-y(1); 目的：让时间运行一个值，参数u就取一个值，顺次取值。u的这11个数是从excel写入的。上面这些程序似乎要加几句，谢谢帮忙。 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★ ★ ★ ★ ★ ★
monitor2885(金币+5,VIP+0):y(end,:)中的end和冒号分别指代谁呢？ 10-15 00:10
sunxiao(金币+2,VIP+0):很好，很有耐心 10-15 05:00
monitor2885(金币+5,VIP+0):y(end,:)中的end和冒号分别指代谁呢？ 10-15 00:10
sunxiao(金币+2,VIP+0):很好，很有耐心 10-15 05:00
|
for i=1:length(t0)-1 t=[t0(i), t0(i+1)]; %u取值的不同，将积分区间依次分成了10段 u=u0(i); [t,y]=ode45(@ivpodefun,t,y0); y0=y(end,: ); % 设置下一积分区间的初值，因为y0=[0 2]只是t=0时的初值。在积分区间t=[t0(i+1), t0(i+2)]时，初值应该变为t0(i+1)对应的y(end,: )值。 tt=[tt;t]; % 将本次积分区间点与前面的积分区间点合并 yy=[yy;y]; % 将本次积分区间点对应的函数值与前面的函数值合并 end 有一个缺点就是由于分段区间是闭区间，所以在分段点(1 2 3 4 5 6 7 8 9)上, tt，yy在有二次重复 可以 改成 for i=1:length(t0)-1 t=[t0(i), t0(i+1)]; u=u0(i); [t,y]=ode45(@ivpodefun,t,y0); y0=y(end,: ); tt=[tt;t(1:end-1)]; yy=[yy;y(1:end-1,: )]; end tt=[tt;t(end)]; yy=[yy;y(end,: )]; 具体的区别你可以显示一下看看 |
8楼2009-10-14 18:57:13
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
2楼2009-10-14 12:32:25
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
3楼2009-10-14 12:36:53
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
nono2009(金币+3,VIP+0):谢谢热心帮助。欢迎常来交流。 10-14 13:12
monitor2885(金币+10,VIP+0):的确牛B，还得继续请教。 10-14 14:13
nono2009(金币+3,VIP+0):谢谢热心帮助。欢迎常来交流。 10-14 13:12
monitor2885(金币+10,VIP+0):的确牛B，还得继续请教。 10-14 14:13
|
function test clear clc global u u0=xlsread('data.xls');% u在excel里面是一列数，1；2；3；4；5；6；5；4；3；2；1； t0=0:1:10; y0=[0 2]; tt=[];yy=[]; for i=1:length(t0)-1 t=[t0(i), t0(i+1)]; u=u0(i); [t,y]=ode45(@ivpodefun,t,y0); y0=y(end,: ); tt=[tt;t]; yy=[yy;y]; end figure plot(tt,yy(:,1),'ro-',tt,yy(:,2),'b^-') function dydt=ivpodefun(t,y) global u dydt=[y(2); u*(1-y(1)^2)*y(2)-y(1)]; [ Last edited by change0618 on 2009-10-14 at 13:01 ] |
4楼2009-10-14 12:49:47
