24小时热门版块排行榜

返回列表

当前主题已经存档。

zhaoshans

铁杆木虫 (小有名气)

应助: 0 (幼儿园)
金币: 5289.6
红花: 1
帖子: 260
在线: 39小时
虫号: 466956
注册: 2007-11-24
性别: MM
专业: 环境工程

[交流] 【求助】ode45问题【已完成】

2个M文件的时间范围和步长一样，在第2个M文件中参数u 是变化的，且u也是21*1的，u 的第一个值是10，从第二个开始都是第1个M文件中y(1) 的解，（通过输入y(:,1)可以得到所有21个解）。在第2个M文件中运行有问题，请高手指点，帮忙修改程序，而且能运行得了。谢谢

%第一个M文件
function dy=ivpode1(t,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=(1-y(1)^2)*y(2)-y(1);

%第二个M文件
function df=ivpode2(t,f)
global y
u=ones(length(t),1)*10;
for i=1:1:20
u(i<=t&t<(i+1))=y(i,1);
end
df=zeros(2,1);
df(1)=f(2);
df(2)=u*(1-f(1)^2)*f(2)-f(1);

t=0:1:20;
>> [t,y]=ode45(@ivpode1,t,[2,0]);
>> [t,f]=ode45(@ivpode2,t,[2,0]);
??? Attempted to access y(7,1); index out of bounds because size(y)=[6,1].

Error in ==> ivpode2 at 5
u(i<=t&t<(i+1))=y(i,1);

Error in ==> funfun\private\odearguments at 110
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...

[ Last edited by nono2009 on 2009-9-25 at 07:02 ]

回复此楼

» 猜你喜欢

1楼 2009-05-20 15:53:26

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

sitonyl

铁杆木虫 (文坛精英)

应助: 33 (小学生)
金币: 12935.9
红花: 3
帖子: 10448
在线: 338.8小时
虫号: 722261
注册: 2009-03-14
性别: GG
专业: 导航、制导与传感技术

★ ★
woshilsh(金币+2,VIP+0):感谢帮助，期待大家点评 5-20 23:25

按照下述顺序调用即可解决问题：

global y
t=0:1:20;
[t,y]=ode45(@ivpode1,t,[2,0]);
[t,f]=ode45(@ivpode2,t,[2,0]);

赞一下(10人)

回复此楼

2楼2009-05-20 22:14:32

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

zhaoshans

铁杆木虫 (小有名气)

应助: 0 (幼儿园)
金币: 5289.6
红花: 1
帖子: 260
在线: 39小时
虫号: 466956
注册: 2007-11-24
性别: MM
专业: 环境工程

还得麻烦你帮忙，谢谢

谢谢关注，但是我按照下面指令输入后，还是运行有问题
global y
t=0:1:20;
[t,y]=ode45(@ivpode1,t,[2,0]);
[t,f]=ode45(@ivpode2,t,[2,0]);

Warning: Failure at t=5.188906e+000. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.421085e-014) at time t.
> In ode45 at 355

赞一下

回复此楼

3楼2009-05-21 11:44:16

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

sitonyl

铁杆木虫 (文坛精英)

应助: 33 (小学生)
金币: 12935.9
红花: 3
帖子: 10448
在线: 338.8小时
虫号: 722261
注册: 2009-03-14
性别: GG
专业: 导航、制导与传感技术

★ ★ ★ ★ ★ ★ ★ ★
zhaoshans(金币+5,VIP+0):很耐心的帮忙 5-21 14:46
woshilsh(金币+3,VIP+0):问题得到楼主认可，嘉奖3金币！ 5-21 15:19
woshilsh(金币+0,VIP+0):兄弟以后常来哈！ 5-21 15:19

引用回帖:

Originally posted by zhaoshans at 2009-5-21 11:44:
谢谢关注，但是我按照下面指令输入后，还是运行有问题
global y
t=0:1:20;
[t,y]=ode45(@ivpode1,t,[2,0]);
[t,f]=ode45(@ivpode2,t,[2,0]);

Warning: Failure at t=5.188906e+000. Unable to meet int ...

那是计算的精度不够，也就是积分计算中选择的步长太大的原因，可以试着更改一下积分步长，如
t=0:0.01:20;
应该就没问题了。

赞一下(21人)

回复此楼