常微分方程(odes)的matlab数值解法.pdfVIP

Matlab 中常微分方程数值解法讲解 1.ODE 解算器简介3 2.微分方程转换5 3. / 8 刚性非刚性问题 4. 隐式微分方程(IDE) 10 5.微分代数 程 15 方 (DAE) 6.延迟微分方程(DDE) 18 7.边值问题(BVP) 20 2 1.ODE 解算器简介 先来认识下常微分方程(ODE)初值问题解算器(solver) [T,Y,TE, YE,IE] = odesolver(odefun,tspan,y0,options) sxint = deval(sol,xint) 解算器(odesolver) 解算器 问题类型 精确度 说明 非 刚 性 ode45 中等 4-5 阶龙格库塔，对以所有问题的首先解算器 (nonstiff) 基于Bogacki-Shampine 2-3 阶Runge-Kutta 公式，有时对轻度的刚度方程，它 ode23 非刚性 低 可以比ode45 有更好的效果，在相 的精度是，要比ode45 更小的步长 变阶次Adams-Bashforth-Moutlon 算法，此算法使用前几次节点上的值来计算 ode113 非刚性 低到高 当前节点处的解，因此在相 jingdu 下，比ode45 和ode23 更快些，适用用 于高阶或者需要大量计算的问题，但不适合不连续的系统 刚性系统的变阶次多步解法，此算法使用最新的数值差分公式，如果使用 ode15s 刚性(stiff) 低到中 ode45 计算比较慢，可以使用它试试 刚性系统固定阶次的单步解法，由于它是单步方法，故有时要比ode15s 速度 ode23s 刚性 低 快些，对具体问题，到底哪个好些，大家可以 时尝试下 ode23t 中等刚性 低 ode23tb 刚性 低 相关优化选项 参数 ode45 ode23 ode113 ode15s ode23s ode23t ode23tb relto abstol √ √ √ √ √ √ √ normcontrol outputfcn outputsel √ √ √ √ √ √ √ refine stats nonnegative √ √ √ √* × √* √* events √ √ √ √ √ √ √ maxstep √ √ √ √ √ √