| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1994 人关注过本帖
标题:ode45解3阶微分方程,2阶我会,3阶就晕了.
只看楼主 加入收藏
lvemealone
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2007-5-30
收藏
 问题点数:0 回复次数:2 
ode45解3阶微分方程,2阶我会,3阶就晕了.
方程如下:x'''+tx''-xx'+1=0
我的m文件:
function xdot=vdpol(t,x)
xdot(3)=-t*x(3)+x(1)*x(2)-1;
xdot(2)=x(3);
xdot(1)=x(2);
x(1)=x;
窗口命令:
>> t0=0;tf=20;
>> x0=[0;0];
>> [t,x]=ode45('vdpol',[t0,tf],x0);
错误提示:
??? Attempted to access x(3); index out of bounds because numel(x)=2.

Error in ==> vdpol at 2
xdot(3)=-t*x(3)+x(1)*x(2)-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, odeFcn, ...

搜索更多相关主题的帖子: 微分方程 xdot args funfun 
2007-05-30 21:02
abingchem
Rank: 6Rank: 6
等 级:贵宾
威 望:24
帖 子:716
专家分:0
注 册:2004-12-30
收藏
得分:0 
function xdot=vdpol(t,x)
xdot=zeros(3,1);
xdot(3)=-t.*x(3)+x(1)*x(2)-1;
xdot(2)=x(3);
xdot(1)=x(2);
?x0=[0;0;0];
?[t,x]=ode45('vdpol',[t0,tf],x0);

5.3版本下运行

此人已消失
2007-06-20 09:23
seidel
Rank: 4
等 级:贵宾
威 望:11
帖 子:217
专家分:20
注 册:2007-1-31
收藏
得分:0 
边界条件少啊!“index out of bounds because numel(x)=2.”

究竟离别只是一场退出还是美好记忆的延续,    最美好的爱究竟是成全还是守候?
2007-06-22 10:29
快速回复:ode45解3阶微分方程,2阶我会,3阶就晕了.
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.028336 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved