06 常微分方程21855.pptVIP

MATLAB数学实验 第六章 常微分方程 陆利正 * 主要内容 6.1 预备知识：常微分方程 6.2 解常微分方程的MATLAB指令 6.3 计算实验：Euler法和刚性方程组 6.4 建模实验: 导弹系统的改进 6.1 预备知识: 常微分方程 1. 微分方程的概念 常微分方程: f (t, y, y ?, y ??, …, y (n)) = 0 微分方程组: 联系一些未知函数的一组微分方程 线性常微分方程 y (n) + a1(t)y (n－1) + … + an－1(t)y ? + an (t)y = b(t) 若ai (t) ( i = 1, …, n) 与 t 无关, 称为常系数的 若 b(t) = 0, 称为齐次的 6.1 预备知识: 常微分方程 2. 初等积分法 分离变量法等 3. 常系数线性微分方程 线性常微分方程的解为一个特解和相应的齐次微分方程通解的叠加 齐次微分方程的解可用特征根法求得 例6.1 求 x ?? + 0.2 x ? + 3.92 x = 0 的通解 特征方程: l2 + 0.2l +3.92 = 0 roots([1 0.2 3.92]) ? + ? i = －0.1 ? 1.9774 i 通解为 x(t) = Ae? t cos( ? t ) + Be? t sin( ? t ) 6.1 预备知识: 常微分方程 4. 初值问题数值解 数值解法 寻求解 y (t) 在一系列离散节点 t0 … tn tf 上的近似值 yk (k = 0, 1, …n) hk = tk+1 － tk 为步长, 通常取为常量 h 6.1 预备知识: 常微分方程 已知一个 n 阶方程 y (n) = f (t, y, y ?, y ??, …, y (n－1)) 设 y1 = y, y2 = y ?, …, yn = y (n－1) 则, 可化为一阶方程组 6.2 解常微分方程MATLAB指令 初值问题求解? 常用格式 [t, y] = ode45(odefun, tspan, y0) odefun：表示 f (t, y) 的函数句柄或Inline函数 t 是标量, y 是标量或向量 tspan： 若为[t0, tf], 表示自变量初值t0和终值tf 若为[t0,t1, …, tn], 表示输出节点列向量 y0： 表示初值向量 t： 表示节点列向量 (t0, t1, …, tn)T y： 数值解矩阵,每一列对应 y 的一个分量 若无输出参数, 则作出图形 6.2 解常微分方程MATLAB指令 初值问题求解? 完整格式 [t, y] = ode45(odefun, tspan, y0, odefun：表示 f (t, y) 的函数句柄或Inline函数 t 是标量, y 是标量或向量 tspan： 若为[t0, tf], 表示自变量初值t0和终值tf 若为[t0,t1, …, tn], 表示输出节点列向量 y0： 表示初值向量 t： 表示节点列向量 (t0, t1, …, tn)T y： 数值解矩阵,每一列对应 y 的一个分量 若无输出参数, 则作出图形 例6.2 解微分方程 y ? = y － 2t/y, y(0) = 1, 0 t 4 odefun = inline(y-2*t/y, t, y); [t, y] = ode45(odefun,[0,4],1) plot(t, y, o-) [t, y] = ode45(odefun, 0:1:4, 1) 例6.3 解微分方程组 例6.4 求解微分方程组 解 引入辅助变量 6.2 解常微分方程MATLAB指令 2. 边值问题解法? sinit = bvpinit(tinit, yinit) 由在粗略节点tinit的预估解yinit生成粗略解网络sinit sol = bvp4c(odefun, bcfun, sinit) bcfun为边值条件函数 sinit是由bvpinit得到的粗略解网络 sol是一个结构, sol.x为求解节点, sol.y是y(t)的数值解 sx = deval(sol, ti) 计算由bvp4c得到的解在ti的值 例6.5 求解边值问题 解 首先改写为标准形式 6.3