MATLAB微分方程..docVIP

在下面的初值问题中，有两个未知函数：x1(t)和x2(t),并用以下式子表达其微... Matlab关于微分方程的解法 页码，1/11  MATLAB使用龙格-库塔-芬尔格（Runge-Kutta-Fehlberg）方法来解ODE问题。在有限点内计算求解。而 这些点的间距有解的本身来决定。当解比较平滑时，区间内使用的点数少一些，在解变化很快时，区间内应使 用较多的点。 为了得到更多的有关何时使用哪种解法和算法的信息，推荐使用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 可得到更多信息。同时可以用命令numjac来计算Jacobi函数。 [t,x]=solver(str,t,x0,val) 此方程的求解过程同上，结构val包含用户给solver的命令。参见odeset和表1，可得 到更多信息。 Ode45 Ode23 Ode113 Ode23t Ode23s Ode15s Ode23tb 此方法被推荐为首选方法。 这是一个比ode45低阶的方法。 用于更高阶或大的标量计算。 用于解决难度适中的问题。 用于解决难度较大的微分方程组。对于系统中存在常量矩阵的情况也有用。 与ode23相同，但要求的精度更高。 用于解决难度较大的问题，对于系统中存在常量矩阵的情况也有用。 Set=odeset(set1,vak1,set2,val2,…) 返回结构set,其中包含用于ODE求解方程的设置参数， 的信息参见表1。 Odeget(set,’set1’) 返回结构set中设置set1的值。 有关可用设置  有许多设置对odeset控制的ODE解是有用的，参见表1。例如，如果要在求解过程中画出解的图形，可以 输入：inst=odeset(‘outputfcn’,’odeplot’);.也可使用命令odedemo。 表1 ODE求解方程的设置参数 RelTol 给出求解方程允许的相对误差 AbsTol 给出求解方程允许的绝对误差  /matlab-pde/ODEmatlab.htm  2004-7-27 在下面的初值问题中，有两个未知函数：x1(t)和x2(t),并用以下式子表达其微... Refine 给出与输入点数相乘的因子 页码，2/11 OutputFcn 这是一个带有输入函数名的字符串，该字符串将在求解函数执行的每步被调用：odephas2(画出2D 的平面相位图)。Odephas3(画出3D的平面相位图)，odeplot(画出解的图形)，odeprint(显示中间结 果) OutputSel 是一个整型向量。指出哪些元素应该被传递给函数，特别是传递给OutputFcn Stats 如果参数Stats为on,则将统计并显示出计算过程中资源消耗情况 Jacobian… 如果编写ODE文件代码以便F（t,y,jocobian）返回dF/dy,则将jacobian设置为on Jconstant…如果雅可比数df/dy是常量，则将此参数设置为on Jpattern… 如果编写ODE文件的编码以便函数F（[],[],jpattern）返回带有零的稀疏矩阵并输出非零元素dF/fy,则 需将Jpattern设置为on Vectorized…如果编写ODE文件的编码以便函数F(t,[y1,y2……])返回[F(t,y1) F(t,y2)…],则将此参数设置成on Events… 如果ODE文件中带有参数‘events’，则将此参数设置成on Mass… 如果编写ODE文件编码以实现函数F(t,[],‘mass’)返回M和M(t),应将此参数设置成on MassConstant…如果矩阵M(t)是常量，则将此参数设置成o