matlab中ode的用法,关于matlab 的ode45用法

3.6.2

龙格-

库塔方法

改进的欧拉法比欧拉法精度高的原因在于，它在确定平均斜率时，多取了一个点的斜

率值。这样，如果我们在[Xi,X(i+1)]上多取几个点的斜率值，然后对它们作线性组合得到平均

斜率，则有可能构造出精度更高的计算方法。这就是龙格-库塔法的基本思想。龙格-库塔

法可看作是欧拉法思想的提高，属于精度较高的单步法。

龙格-库塔法是求解常微分方程初值问题的最重要的方法之一。MATLAB中提供了几

个采用龙格-库塔法来求解常微分方程的函数，即ode23，ode45，ode113

，ode23s

，ode15s

等，其中最常用的函数是

ode23(

二三阶龙格-库塔函数)和ode45(

四五阶龙格-库塔函数)，

下面分别对它们进行介绍。

1

．二三阶龙格-

库塔函数(ode23)

函数

ode23

的调用格式如下：

(1)

[T,Y]=ODE23('F',TSPAN,Y0)

输入参数中的'F'

是一个字符串，表示微分方程的形

式，也可以是

f

(x

,

y

)的M

文件。TSPAN=[T0

TFINAL]表示积分区间，Y0表示初始条件。

函数

ode23