Matlab微分方程的解法精选.pdfVIP

在下面的初值问题中，有两个未知函数：x1(t)和x2(t),并用以下式子表达其微... 页码，1/11 Matlab 于微分方程的解法 MATLAB使用龙格 库塔 芬尔格（ ）方法来解ODE 问题。在有限点内计算求解。而 - - Runge-Kutta-Fehlberg 这些点的间距有解的本身来决定。当解比较平滑时，区间内 用的点数少一些，在解变化很快时，区间内应 用较多的点。 为了得到更多的有关何时使用哪种解法和算法的信息，推荐使用helpdesk 。所有求解方程通用的语法或句法在 命令集中头两行给出。时间间隔将以向量t=[t0,tt]给出。 命令ode23可以求解(2,3)阶的常微分方程组，函数ode45使用(4,5)阶的龙格-库塔-芬尔格方法。注意，在这种情 况下 的微分不 的转置。 x x x 在命令集中solver将被诸如ode45 函数所取代。 命令集 龙格 库塔 芬尔格方法 - - 计算 或由字符串 给定的 的值，部分解已在向量 中给出。在向量 [time,x]=solver(str,t,x0) ODE str ODE time time 中给出部分解，包含的 时间值。还有部分解在矩阵x 中给出，x 的列向量是每 个方程在这些值下的解。对于标量问题，方程的解将在向量 中给出。这些解 x 在时间区间t(1)到t(2)上计算得到。其初始值 x0 即x(t(1)).此方程组有str指定的 文件中函数表示出。这个函数需要两个参数：标量 和向量 应该返回向量 M t x, 即 的导数 。因为对标量 来说， 和 都是标量。在 文件中输入 x ( x ) ODE x x M odefile 可得到更多信息。同时可以用命令numj ac来计算Jacobi 函数。 此方程的求解过程同上，结构 包含用户给 的命令。参见 和表 ，可得 [t,x]=solver(str,t,x0,val) val solver odeset 1 到更多信息。 Ode45 此方法被推荐为首选方法。 Ode23 这 一个比ode45低阶的方法。 Ode113 用于更高阶或大的标量计算。 Ode23t 用于解决难度适中的问题。 Ode23s 用于解决难度较大的微分方程组。对于系统中存在常量矩阵的情况也有用。 Ode15s 与ode23相同，但要求的精度更高。 Ode23tb 用于解决难度较大的问题，对于系统中存在常量矩阵的情况也有用。 返回结构 其中包含用于 求解方程的设置参数， 有关可用设置 Set=odeset(set1,vak1,set2,val2, ) set, ODE 的信息参见表1。 返回结构 中设置 的值。 Odeget(set,’set1’) set set1 有许多设置对 控制的 解 有用的，参见表 。例如，如果要在求解过程中画出解的图形，