实用 MATLAB教学资料-27-28常微分方程续.pptVIP

节 点 含 义 设 置 workspace中的变量名 初始网格节点 用于计算网格初值 solinit = bvpinit(x,yinit) 结构数组solint中的solint.x 实际计算网格节点 MATLAB计算时采用的节点 软件根据精度要求自行设置，用户无法修改。 结构数组Sol中的Sol.x 输出网格节点 根据用户需要，对输出数据的疏密进行设置 sx = deval(sol,x,idx) x 三种节点 例7 求二阶常微分方程第一边值问题在[0,4]上的解,要求ε 1×10-5 solinit = bvpinit(x,yinit) options =bvpset(name1,value1,...) sol= bvp4c(odefun,bcfun,solinit,options) sx = deval(sol, x, idx) 实用MATLAB 刚性问题ode solver 求解器Solver 算法 ODE类型 精度 适用场合 ode15s Gear法 刚性 低~中 若因刚性导致ode45失效时，可尝试使用；可求解DAEs(微分代数方程) ode23s 单步2阶Rosenbrock法 刚性 低 解决精度要求低、刚性问题。计算时间比ode15s短 ode23t 梯形Euler公式 中等刚性 低 适用于中等刚性的场合；求解DAEs ode23tb 隐式RK公式 刚性 低 解决精度要求低、刚性问题。精度较低时，计算时间比ode15s短 例4 分别用ode15s和ode45解如下刚性方程 （1）tspan=[0,300] （2）tspan=[0,3000] 5、隐式常微分方程 定解条件：函数和导函数初值 5、隐式常微分方程 ode15i [x,y]=ode15i(odefun,tspan,y0, dydx0,…) 定义隐式方程，三个输入x,y,y’ 函数初值 导函数初值 例5 求隐式常微分方程在[0,12]的解 [x,y]=ode15i(odefun,tspan,y0, dydx0) 6、延迟常微分方程 描述延迟系统： 系统随时间的演化，不仅依赖于系统当前的状态，且依赖于系统过去某一时刻或若干时刻的状态 延迟常数 dde23 适用场合： 解恒定迟滞的延迟微分方程 原理：隐式Runge-Katta算法 dde23 options =ddeset(name1,value1,...) sol=dde23(ddefun,lags,history,tspan,options) lags-延迟常数，行向量 dde23 sol=dde23(ddefun,lags,history,tspan,options) function dydt = ddefun(x,y,Z) 其中z–待解函数的历史状态，矩阵 τj = lags(j) 函数式m文件定义 dde23 sol=dde23(ddefun,lags,history,tspan,options) history-定义历史状态， 待解函数y(t)在t ≤ t0的值， 列向量 options =ddeset(name1,value1,...) 参数设置 RelTol-相对误差 AbsTol-绝对误差 NormControl-范数误差控制 InitialStep-初始步长 MaxStep –最大步长 …… ddesd 适用场合： 处理迟滞是自变量函数的延迟微分方程 ddesd sol=ddesd(ddefun,lags,history,tspan,options) lags:一般延迟,lags=p(t) 用匿名函数或函数式m文件定义 例6 解延迟常微分方程 1、边值问题的描述 在自变量的两端给定边界条件 10.2 边值问题 一、边值问题数值解法 Boundary Value Problems 有限差分法 将边值问题（包括方程及边界条件）离散化，求离散点上函数近似值yi。 2、数值求解方法 Boundary Value Problems bvp4c, bvp5c 功能：解常微分方程边值问题 算法：有限差分法 二、MATLAB功能函数 Boundary Value Problems sol=bvp4c(odefun,bcfun,solinit,options) bvp4c Boundary Value Problems 定义微分方程 定义边界条件 初始猜测 sol=bvp4c(odefun,bcfun,solinit,options) (1)odefun 定义降阶后形成的常微分方程组, 同ode solver 必须返回导函数列向量 sol=bvp4c(odefun,@bcfun,s