| 查看: 631 | 回复: 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 ] |
回复此楼» 猜你喜欢
确定了，国自然21号放榜
已经有4人回复
-
静等基金结果
已经有16人回复
-
filecode
已经有6人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
-
国自然结果
已经有8人回复
-
FileCode能看出啥？
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有22人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
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
sitonyl
铁杆木虫 (文坛精英)
- 应助: 33 (小学生)
- 金币: 12935.9
- 红花: 3
- 帖子: 10448
- 在线: 338.8小时
- 虫号: 722261
- 注册: 2009-03-14
- 性别: GG
- 专业: 导航、制导与传感技术
4楼2009-05-21 12:05:14
zhaoshans
铁杆木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5289.6
- 红花: 1
- 帖子: 260
- 在线: 39小时
- 虫号: 466956
- 注册: 2007-11-24
- 性别: MM
- 专业: 环境工程
5楼2009-05-21 14:46:09
sxhys
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1096.9
- 红花: 1
- 帖子: 184
- 在线: 22.8小时
- 虫号: 531677
- 注册: 2008-03-23
- 性别: GG
- 专业: 通信理论与系统
6楼2009-05-30 10:07:06
