Taylor公式与科学计算.docVIP

Taylor公式与科学计算 ----------导数数值计算 课程名称：工科数学分析 课程总学时：240学时 授课对象：北航工科专业 北京航空航天大学理学院数学系 杨小远 Taylor公式与科学计算 ----------导数数值计算 摘要：本文深入研究了微积分著名Taylor公式的教学内容，我们通过Taylor公式的在导数数值计算的应用，阐述了Taylor在科学计算中的应用，拓宽了学生的知识领域，激发了学生的学习积极性和提高思考问题的能力。主要教学内容如下： 科学计算研究意义 计算机如何实现导数计算以及计算精度分析 导数数值计算存在的问题 Taylor公式解决问题 李查逊外推（Richardson） 进一步的思考 通过实际教学效果看，学生普遍能接受本次教学内容，很好完成了相关的思考题。 科学计算研究的意义 科学计算的兴起是20世纪后半叶最重要的科技进步之一计算与理论及实验相并列已经成为当今世界科学活动的第三种手段为把信息和数据变成知识从而探索科学未知促进技术创新加强国防建设保障国家安全计算将起不可替代的重要作用当前科学计算正在向大规模和高性能发展要达到全物理、全系统、三维、高分辨、高逼真的数值模拟高效的计算方法重要大规模计算提出的世界性难题已形成科学计算的学科前沿求解由实际问题得到的复杂的偏微分方程不仅计算规模大更由于非线性、多尺度、长时间、不适定、多区域、高病态等特点使计算格外困难现有的算法远不能满足需求这正是必须解决的关键问题Taylor公式在科学计算中的作用。使得学生对科学计算有初步认识，拓展学生的知识面。 导数的数值计算与误差分析 我们知道导数的定义如下 （1） 上述过程是无限逼近的过程，但是计算机实现需要的有限次运算，因此给出几种近似的计算公式： 对于任何近似计算都需要对算法的误差进行分析，通过Taylor公式我们有： 类似有 通过（5），（6），（7）可以知道（7）的计算精度更高，我们可以看到阶的运算在分析算法性能是一个重要有利的工具。 在科学计算中阶的运算是分析算法性能的重要工具。 导数数值计算存在的问题 下面的例子给出计算在的导数值的实验结果。 表1 h 逼近值 逼近误差 2 0.3360 -0.012447 1 0.3564 -0.002847 0.2 0.3535 0.000053 0.1 0.3530 0.000553 0.02 0.3550 -0.001447 0.01 0.3500 0.003553 0.002 0.3500 0.003553 0.001 0.3000 0.053553 0.0002 0.3000 -0.146447 从表1我们发现下面两个事实： 从表中看出时候计算效果最佳； h取得比小时，计算的效果越来越差。 这两个现象与我们的设想有着“天壤之别”，按照导数的定义越小应该计算精度越高，导致的原因： 用计算机表示任何数字只能是有限位，计算机实现任何运算都会有舍入误差，在科学计算中必须充分重视舍入误差对计算结果的影响，这也是科学计算重要而十分艰难的重要研究课题。 为此我们引入几个几个概念： 定义1：记为的近似数，为的绝对误差。 为的相对误差。 定义2：假设为整数，，如果，则称有位有效数字。 例如：的近似值具有5位有效数字，的近似有2位有效数字。有效数字是从第一位不等于0的数算起。 由于公式（7）当非常小的时候出现两个相近数的相减，使得有效数字大量减少，致使产生误差。例如有四位有效数字，只有1位有效数字。 因此在计算机计算的过程中应该避免两个相近数字的减法运算。任何一个工程或者科学问题，其数值计算的次数是巨大和海量的，我们必须设计有效的算法控制舍入误差的传播。 解决问题的办法 这里我们将用Taylor公式解决问题。 由Taylor公式，得到 因此有 假设 则有 因此得到 令 依次类推，我们构造序列： 并且有： 上述公式的计算过程如下： (3) （6） (10) （5） （9） （4） （8） 实验结果分析 计算，在的一阶导数值，这里 实验结果如下： 表2 0.0128 696.6346914 623.4601726 625.3455055 625.3334226 0.0064 641.7538023 625.2276722 625.3336144 0.0032 629.35920