Matlab关于微分方程解法.pdfVIP

在下面的初值问题中，有两个未知函数：x1(t)和x2(t),并用以下式子表达其微... 页码，1/11 Matlab 于微分方程的解法 MATLAB使用龙格 库塔 芬尔格（ ）方法来解ODE 问题。在有限点内计算求解。而 - - Runge-Kutta-Fehlberg 这些点的间距有解的本身来决定。当解比较平滑时，区间内 用的点数少一些，在解变化很快时，区间内应 用较多的点。 l 为了得到更多的有关何时使用哪种解法和算法的信息，推荐使用helpdesk 。所有求解方程通用的语法或句法在 命令集中头两行给出。时间间隔将以向量t=[t0,tt]给出。 a i 命令ode23可以求解(2,3)阶的常微分方程组，函数ode45使用(4,5)阶的龙格-库塔-芬尔格方法。注意，在这种情 r 况下 的微分不 的转置。 x x x 在命令集中solver将被诸如ode45 函数所取代。 T 命令集 龙格 库塔 芬尔格方法 - - 计算 或由字符串 给定的 的值，部分解已在向量 中给出。在向量 [time,x]=solver(str,t,x0) ODE str ODE time time F m 中给出部分解，包含的 时间值。还有部分解在矩阵x 中给出，x 的列向量是每 o 个方程在这些值下的解。对于标量问题，方程的解将在向量 中给出。这些解 D c x . 在时间区间t(1)到t(2)上计算得到。其初始值 x0 即x(t(1)).此方程组有str指定的 d 文件中函数表示出。这个函数需要两个参数：标量 和向量 应该返回向量 M P r t x, 即 的导数 。因为对标量 来说， 和 都是标量。在 文件中输入 x ( x ) ODE x x M odefile a 可得到更多信息。同时可以用命令numj ac来计算Jacobi 函数。 z 此方程的求解过程同上，结构 包含用户给 的命令。参见 和表 ，可得