7 微分方程求解.pptxVIP

第七章常微分方程的差分法本章内容§7.1 引言§7.2 欧拉法§7.3 龙格—库塔法§7.4 亚当斯法§7.5 收敛性与稳定性小结本章要求1. 熟悉各种Euler格式2.掌握龙格—库塔法3.了解亚当斯法3.了解差分格式的收敛性和稳定性§7.1 引言本节内容一. 问题提出二. 两类求定解问题三. 数值解法含义四. 差商公式返回章节目录§7.1 引言一. 问题提出 有一个或多个导数及其函数的方程式称为微分方程，在工程中常遇到求解微分方程的问题。§7.1 引言 很多微分方程的解不能用初等函数来表示，有时即使能够用解析式表示其解，但计算量太大而不实用(表达式过于复杂)。 需要用数值方法来求解，一般只要求得到若干个点上的近似值或者解的简单的近似表达式(精度要求满足即可)。§7.1 引言二. 两类求定解问题实际中求解常微分方程的所谓定解问题有两类:初值问题和边值问题定解指已知因变量和/或其导数在某些点上是已知的(约束条件)。1. 边值问题 约束条件为已知，在自变量的任一非初值上，已知函数值和/或其导数值，如 常常可以将边解问题转化为初值问题求解。§7.1 引言2. 初值问题§7.1 引言单步法： 利用前一个单步的信息(一个点)，在y=f(x) 上找下一点yi， 有欧拉法，龙格－库格法等。多步法： 利用前面两步或两步以上的信息，在y=f(x) 上找下一点yi ， 有改进欧拉法，亚当斯法等。初值问题的常见解法§7.1 引言§7.1 引言§7.1 引言李普希兹条件: 一条曲线上任意两点连线的斜率的绝对值都小于某一个数L§7.1 引言三. 数值解法含义所谓数值解法, 就是设法将常微分方程离散化, 建立差分方程, 给出解在一些离散点上的近似值。微分方程的数值解: 设方程问题的解y(x)的存在区间是[a,b], 令a= x0 x1… xn =b, 其中hk=xk+1-xk, 如是等距节点h=(b-a)/n, h称为步长。由于y(x)的解析表达式不容易得到或根本无法得到, 我们用数值方法求得y(x)在每个节点xk上y(xk)的近似值, 用yk表示, 即yk≈y(xk), 这样y0, y1, ..., yn称为微分方程的数值解。§3.1 引言四. 数值求导的差商公式§7.1 引言函数y(x)在某一点x=a的导数若取差商作为导数的近似值一阶差商§7.2 欧拉方法本节内容一. 欧拉格式二. 梯形格式三. 改进的欧拉格式四. 欧拉法的精度返回章节目录§7.2 欧拉方法一. 欧拉（Euler）格式近似关系式差分格式§7.2 欧拉方法左矩形复习：y梯形公式梯形公式y=?(x)左矩形公式oaxb右矩形公式右矩形公式中矩形公式中矩形公式左矩形公式§7.2 欧拉方法/47向前差商近似导数记为§7.2 欧拉方法? 欧拉格式：/* Euler’s Method */向后差商近似导数?+y(x)yhf(x,y(x))1011=+=-yyhf(x,y)(i0,...,n1)?隐式欧拉格式 /* implicit Euler method */+++111iiii§7.2 欧拉方法由于未知数 yi+1 同时出现在等式的两边，不能直接得到，故称为隐式 /* implicit */ 欧拉格式，而前者称为显式 /* explicit */ 欧拉格式。一般先用显式计算一个初值，再迭代求解。§3.2 欧拉方法当i=0，即x0=0时， y1=0+0.2(1-0×0)=0.2000当i=1 ，即x1=0.2时， y2=0.2+0.2(1-0.2×0.2)=0.3920当i=2 ，即x2=0.4时， y3=0.392+0.2(1-0.4×0.392)=0.56064当i=3 ，即x3=0.6时， y4=0.56064+0.2(1-0.6×0.56064)=0.6933632当i=4 ，即x4=0.8时， y4=0.6933632+0.2(10.8×0.6933632)=0.782425088列表计算如下：§7.2 欧拉方法i012345xi00.20.40.60.81yi0.00000.2000.39200.56060.69340.7824/47§7.2 欧拉方法hxiyiy(xi)y(xi)-yih=0.20.000.400.801.201.602.000.000000.376310.542280.527090.466320.406820.000000.344830.487800.491800.449440.400000.00000-0.03148-0.05448-0.03529-0.01689-0.00682h=0.10.000.400.801.201.602.000.000000.360850.513710.509610.458720.404190.0000