| 查看: 505 | 回复: 5 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
zhaoshans铁杆木虫 (小有名气)
|
[交流]
【求助】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 ] |
回复此楼» 猜你喜欢
过年走亲戚时感受到了所开私家车的鄙视链
已经有10人回复
-
今年春晚有几个节目很不错，点赞！
已经有12人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有13人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有12人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复
sitonyl
铁杆木虫 (文坛精英)
- 应助: 33 (小学生)
- 金币: 12935.9
- 红花: 3
- 帖子: 10448
- 在线: 338.8小时
- 虫号: 722261
- 注册: 2009-03-14
- 性别: GG
- 专业: 导航、制导与传感技术
4楼2009-05-21 12:05:14
sitonyl
铁杆木虫 (文坛精英)
- 应助: 33 (小学生)
- 金币: 12935.9
- 红花: 3
- 帖子: 10448
- 在线: 338.8小时
- 虫号: 722261
- 注册: 2009-03-14
- 性别: GG
- 专业: 导航、制导与传感技术
2楼2009-05-20 22:14:32
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
zhaoshans
铁杆木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5289.6
- 红花: 1
- 帖子: 260
- 在线: 39小时
- 虫号: 466956
- 注册: 2007-11-24
- 性别: MM
- 专业: 环境工程
5楼2009-05-21 14:46:09