| 查看: 2325 | 回复: 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函数不支持这样吗？我要怎么解决。欢迎大家探讨，谢谢大家指点！！ |
回复此楼» 猜你喜欢
求个博导看看
已经有12人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有5人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有7人回复
-
需要合成515-64-0，50g，能接单的留言
已经有4人回复
-
自荐读博
已经有4人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有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