数值计算 第七章.ppt

计算机科学与工程系 计算机科学与工程系 数值计算方法 第七章 常微分方程的数值解法 7.1 引言 常微分方程求解 微分方程的分类 一阶常微分方程初值问题的数值解法 7.1.1 常微分方程求解 微分方程的定义 包含自变量、未知函数以及未知函数的导数或微分的方程叫做微分方程 在求解微分方程时，必须附加某种定解条件 微分方程和定解条件一起组成定解问题 微分方程的解析解法 在数学分析中，求解微分方程采用的方法就是解析法，求得解析解为未知函数y(x) 如果想求微分方程某点的函数值，只需将某点的点值x代入所求的未知函数y(x)，即可求得某点的函数值 7.1.1 常微分方程求解 微分方程的数值解法 在计算方法中求解微分方程采用的方法就是差分法。其数值解为某点函数值y(xi)的近似值yi 即给定一组离散点值x0，x1，x2，…xk，就可求得一组 y(x0)，y(x1)，y(x2)，…，y(xk)，此值为精确值，但实际求得的是近似值y0，y1，y2，…，yk 7.1.2 微分方程的分类 未知变量的个数 在微分方程中只含有一个自变量x的方程称作常微分方程，相当于一元函数的范畴。其应用广泛，求解方法也比较成熟 在微分方程中，含有多个自变量x，y，z，t…的方程称为偏微分方程，它相当于多元函数的范畴，其求解比较困难 7.1.2 微分方程的分类 例如：对一个各向同性的物体，温度u(x, y, z, t)在物体内的分布满足偏微分方程 求解的条件 初值问题：已知函数的导函数和给定的初始条件，按照假设的步长值h，确定一系列的值x0，x1，x2，…，然后求得一系列的函数值y(x0)，y(x1)，y(x2)，…，的近似值y0，y1，y2，…，直到求得所需点值xk处的函数值y(xk)的近似值yk为止 7.1.2 微分方程的分类 边值问题：已知函数的导函数和给定的边界条件，求解满足微分方程式的某点函数值 未知函数的个数 一阶常微分方程：只含有一个未知变量x，未知函数的导数阶数为一阶的，而且只有一个未知函数，简言之，方程中只含有x、y、y’的方程 7.1.2 微分方程的分类 一阶常微分方程组：只含有一个未知变量x，未知函数的导数阶数为一阶的，而未知函数的个数是多个，简言之，方程中含有x，y1，y2，…，yn，y1′，y2′，…，yn′的方程 7.1.2 微分方程的分类 未知函数导数的阶 一阶常微分方程：只含有x，y，y’时，就是最常见的微分方程 高阶常微分方程：除了含有x，y，y’以外，还含有未知函数y的高阶导数y′，y〞，… 7.1.2 微分方程的分类 高阶常微分方程组：凡是由x，y1，y1′， y1〞 ， …，y2，y2′，y2〞，…，等组成的微分方程 7.1.3一阶常微分方程初值问题的数值解法 数值方法的基本思想 在解的存在区间上取n + 1个节点 这里 ，i = 0,1, …, n，称为由xi 到xi+1的步长 在这些节点上采用离散化方法，将上述初值问题化成关于离散变量的相应问题 把这个相应问题的解yn作为y(xn)的近似值 7.1.3一阶常微分方程初值问题的数值解法 定理：常微分方程初值问题 设x0∈[a, b]，f(x,y)对x连续并且关于y满足李普希兹条件：存在常数L，使 对所有x∈[a, b]以及任何实数y1，y2都成立，则初值问题的解y(x)存在并且唯一 7.1.3一阶常微分方程初值问题的数值解法 常微分方程初值问题的数值解法 一步法：一步法是在计算yn+1时，只用到xn+1，xn和yn，即前一步的值。因此，有了初值以后就可以逐步往下计算，其代表是龙格－库塔法 多步法：多步法就是在计算yn+1时，除了用到xn+1，xn和yn以外，还要用到xn-p，yn-p (p=1,2,…,k)，即前面k步的值，其代表是亚当斯法 7.2 欧拉法 欧拉法 两步欧拉公式 改进欧拉法 7.2.1.欧拉法 公式的导出 由于y (x0) = y0已给定，因而可以算出 设x1 –x0= h充分小，则近似地有 记 从而我们可以取 7.2.1.欧拉法 作为y (x1)的近似值 利用y1及f (x1, y1)又可以算出y(x2)的近似值 一般地，在任意点xn+1 = (n + 1)h+x0处y(x)的近似值由下式给出 7.2.1.欧拉法 欧拉法的几何意义 一阶常微分方程 的解 是通过点 的一条曲线，称之为微分方程的积分曲线。积分曲线上每一点 的切线斜率 等于函数 在这点的值 7.2.1.欧拉法 7.2.1.