| 查看: 2328 | 回复: 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函数不支持这样吗？我要怎么解决。欢迎大家探讨，谢谢大家指点！！ |
回复此楼» 猜你喜欢
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
要不要辞职读博？
已经有3人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
磺酰氟产物，毕不了业了！
已经有6人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！！对于调用龙格库塔法如何对结果判断反馈
已经有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