数值分析教教案16.docVIP

3.4 常微分方程初值问题数值解的MATLAB实现 在MATLAB中，有多个求解常微分方程数值解的函数命令，在此介绍几个常用的函数。 3.4.1 求常微分方程初值问题数值解的函数 ode23和ode45是求解常微分方程数值解最常用的两个函数，命令中的“ode”是英文常微分方程“Ordinary Differential Equation”的缩写，它们都采用龙格-库塔公式进行数值求解，23和45分别表示使用是2/3阶和4/5阶龙格-库塔公式。它们的调用格式基本相同。在此仅以ode23为例来说明函数的用法。函数ode23的调用格式为： [x,y]=ode23(Fun,Tspan,y0,options) 1. 该命令适用于一阶常微分方程组，，如遇到高阶常微分方程，必须先把它们变成一阶常微分方程组，即状态方程，方可使用。 2. 输入参数“Fun”为定义微分方程组的M-函数文件名，可以在文件名前加写@，或用英文格式单引号界定文件名。 3. 在编辑调试窗口中编写一阶常微分方程组的M-函数文件时，每个微分方程的格式必须都与一致，即等号左端为待求函数的一阶导数，右端函数的变量严格以“先自变量、后函数”的固定顺序输入，的下标表示微分方程的序数。 4. 输入参数“Tspan”规定了常微分方程的自变量取值范围，它以矩阵的形式输入，表示自变量。 5. 输入参数“”表示初始条件向量，。微分方程组中的方程个数