四阶Runge－Kutta法

rushuidetian

有没有求解微分方程组的四阶Runge－Kutta法 Matlab程序 比如Du/Dt=Mu,其中u为69*1阶向量,M为69*69阶矩阵

--------------------------------
请注意你的标题

[ 本帖最后由 sigma665 于 2008-5-27 17:35 编辑 ]

ch_j1985

原帖由 rushuidetian 于 2008-5-27 16:21 发表

                               
登录/注册后可看大图

有没有求解微分方程组的四阶Runge－Kutta法 Matlab程序 比如Du/Dt=Mu,其中u为69*1阶向量,M为69*69阶矩阵

可以参考此帖：
http://www.chinavib.com/forum/viewthread.php?tid=15685&extra=page%3D1%26amp%3Bfilter%3Ddigest

gh688

ode45就是4阶Runge－Kutta法，最简单的办法：
ode45(@myfun,tspan,x0,opts,A);
里面的这几个参数自己定义一下
function xdot=myfun(t,x,A)
xdot=A*x

rushuidetian

我调用库函数 [T,V]=ODE45(@fun,[0,0.2],V0,[ ],M0,M1);
其中  function f=fun(t,V,M0,M1)
          %其中V为列向量M0，M1为矩阵
      f=(M0+M1)*V;
运行结果显示
??? Attempt to execute SCRIPT ode45 as a function.

Error in ==> C:\Documents and Settings\hp\桌面\ode45.m
On line 47  ==> [T,V]=ODE45(@fun,[0,0.2],V0,[],M0,M1);
请问这是怎么回事呢?

[ 本帖最后由 sigma665 于 2008-5-28 13:52 编辑 ]

sigma665

C:\Documents and Settings\hp\桌面\ode45.m

路径怎么会在桌面上？

咕噜噜

路径倒不是很大关系，有有时候勉强也可以，但是你那个ode45写错了，小写的，不是大写的
或者新版软件用大写的？:lol

rushuidetian

路径改了,ode也改成小写了，可是程序运行结果显示和上面还一样。我不明白这很明显是调用库函数，为什么运行显示说要将ode45作为函数处理?不好意思，我才刚开始接触Matlab，请多多指教！

sigma665

which ode45
看看有几个

rushuidetian

Matlab里面应该只有一个ode45   就是四阶Runge-kutta法

sigma665

你看了没

如果有多余的，保留自带的
删除你自己的

rushuidetian

谢谢各位的热心帮忙，我再仔细看看！