rushuidetian 发表于 2008-5-27 16:21

四阶Runge－Kutta法

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

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

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

ch_j1985 发表于 2008-5-27 16:34

原帖由 rushuidetian 于 2008-5-27 16:21 发表 http://www.chinavib.com/forum/images/common/back.gif
有没有求解微分方程组的四阶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 发表于 2008-5-27 16:52

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

rushuidetian 发表于 2008-5-28 11:11

我调用库函数 =ODE45(@fun,,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==> =ODE45(@fun,,V0,[],M0,M1);
请问这是怎么回事呢?

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

sigma665 发表于 2008-5-28 13:51

回复 4楼 的帖子

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

路径怎么会在桌面上？

咕噜噜 发表于 2008-5-28 14:53

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

rushuidetian 发表于 2008-5-28 15:18

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

sigma665 发表于 2008-5-28 15:54

which ode45
看看有几个

rushuidetian 发表于 2008-5-28 16:15

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

sigma665 发表于 2008-5-28 18:05

回复 9楼 的帖子

你看了没

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

rushuidetian 发表于 2008-5-28 18:30

谢谢各位的热心帮忙，我再仔细看看！
页: [1]
查看完整版本: 四阶Runge－Kutta法