| 查看: 2434 | 回复: 1 | |||
[求助]
关于利用ode45函数求解微分方程组
|
|
我想利用ode45函数求解某二阶微分方程组。编程如下： ！！！M函数： function dx=sl(t,x) xx=[1,2,3,,4,5,6,7,8,9]; yy=[0,-1,2,5,-1,4,7,10,-3]; dx=[x(2) 2*x(4)+x(1)+xx(t/0.02) x(4) -2*x(2)+x(3)+3*yy(t/0.02)]; ！！！Matlab主程序 x0=[0;0;0;0]; options=odeset('reltol',1e-8); n=9; A=cell(1,n); for i=1:9 tic [t,A{1,i}]=ode45(@sl,[0,i*0.01,i*0.02],x0,options); toc end 通过运行发现报错： ??? Attempted to access xx(0); index must be a positive integer or logical. Error in ==> sl at 4 dx=[x(2) 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, ... 通过分析发现，其说明对xx矩阵要调用xx（0），这个不符合逻辑，因而出错。我本意是想通过t/0.02而调用xx中的xx(t/0.02)元素。难道ode45函数不支持这样吗？我要怎么解决。欢迎大家探讨，谢谢大家指点！！ |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！！对于调用龙格库塔法如何对结果判断反馈
已经有6人回复
- 高数关于不显含x和y的微分方程求解？
已经有8人回复
- 哪位大神帮忙看一下MATLAB求解微分方程的问题？
已经有3人回复
- 微分方程求解（通解中含hypergeom函数）
已经有17人回复
- 求助Matlab求解含已知变量的微分方程组（变量随时间变化的离散值已知）
已经有12人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- matlab解微分方程组
已经有15人回复
- matlab 计算过程中工作空间变量保存及读取 求解决内存不足的方法
已经有5人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】急求各位大侠，关于用ode45解微分方程
已经有8人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】常微分方程组求解中系数与某变量值关联的问题
已经有12人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
2楼2014-12-10 21:06:37
200