第5章-常微分方程.pptVIP

例 用Euler两步法求解初值问题 解： (1)取步长 h=0.1，首先用Euler方法计算初值， (2)用它代入Euler两步法公式，得 1.0954 1.1832 1.2649 1.3416 1.4142 1.4832 1.5492 1.6125 1.6733 1.7321 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 y(xn) xn 1.1000 1.1918 1.2774 1.3582 1.4351 1.5090 1.5803 1.6498 1.7178 1.7848 Euler 1.1000 1.1836 1.2691 1.3429 1.4186 1.4856 1.5542 1.6163 1.6794 1.7378 Euler两步 例 用Euler两步法的预测校正方法求解初值问题 解： (1)取步长 h=0.1，首先用Euler方法计算初值， (2)用它代入Euler两步法公式，得 (3)用它代入梯形公式，得 1.0954 1.1832 1.2649 1.3416 1.4142 1.4832 1.5492 1.6125 1.6733 1.7321 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 y(xn) xn 1.1000 1.1918 1.2774 1.3582 1.4351 1.5090 1.5803 1.6498 1.7178 1.7848 Euler 1.1000 1.1836 1.2691 1.3429 1.4186 1.4856 1.5542 1.6163 1.6794 1.7378 Euler2步 1.0959 1.1838 1.2656 1.3425 1.4153 1.4845 1.5507 1.6142 1.6754 1.7345 预估校正 ans =1/(x-1+2*exp(-x)) ■ MATLAB 解常微分方程初值问题命令 解析解命令:dsolve(eqn1, ...) 解析解: 例 syms x y dsolve(Dy=y-x*y^2,y(0)=1,x) ■ MATLAB 解常微分方程初值问题命令 数值解命令:ode23(‘f, [a,b], y0) 例 f=inline(y-x.*y.^2); [x,y]=ode23(f,[0,2],1) 1.0000 1.0796 1.2485 1.3755 1.4402 1.4209 1.3374 1.2196 1.0923 0.9709 0.8620 0.7870 0 0.08 0.2621 0.4437 0.6437 0.8437 1.0437 1.2437 1.4437 1.6437 1.8437 2.0000 y x 作 业 1．思考题1中的（a）—（f） 2．习题3 * * 马尔萨斯人口模型：假设某特定区域在 t0 时刻的人口p(t0) = p0为已知的，该区域人口的自然增长率为α。人口的增长与人口的总数成正比，所以 t 时刻的人口总数 p(t)满足如下的微分方程： 生活中常常有这样一类问题： 问 题 的 提 出 这些常微分方程有各种各样的解析方法，但解析方法只能用来求解一些特殊类型的问题，实际问题中归结出来的微分方程主要靠数值解法。 解析法：给出精确解析解。只适合少数简单情况。 近似解法：给出解的近似表达式。如级数法,逐步逼近法。 数值方法：给出方程在离散点上的近似解。它适合计算机求解,应用广泛,具有理论应用价值。 常微分方程的解法： 内容分类： 定解问题 初值问题 边值问题 单步法 Euler方法 Taylor方法和Runge-Kutta方法 多步法 Adams方法和一般线性多部法 线性多部法的收敛性与稳定性 一阶常微分方程初值问题的一般形式： 问题: 求函数 满足 其中: f (x, y) 为已知函数, α是已知值. (可能是观察值或实验值) 基本条件： f (x,y)在D上连续； f (x,y)在D上关于变量y满足Lipschitz连续条件： 设 满足解的存在唯一 对求解区域[a,b]做剖分 构造数值解法的基本思想 在区间[xk, xk+1]上对微分方程做积分,则有 常用等步长: , 则有 将微分方程的准确解记为y(x), 称为步长。 的近似解记为 能不能将微分转化为积分？ 因此，建立节点处近似值yn满足的差分公式 称之为Euler公式. 对右边的积分应用左矩形公式，则有 Euler公式的几何意义 特点：简单，精度低. 例 求解初值问题 解： Euler公式的具体形式为 取步长 h=0.1，那么即可计算该微分方程。 具体结果见下页。