数值分析全套c课件.pptVIP

科学计算的背景 关于计算误差讨论 浮点数与有效数字 算术运算的误差估计 《数值分析》1 ? ? ? ? von Neumann and Goldstine：“高阶矩阵的数值求逆” (1947年) 数值分析——研究用计算机求解 1958年, 前苏联载人飞船 1969年, 美国Apollo 登月 1994年, 美国GPS运行 2/16 数学问题的数值计算方法及其理论 方程组求解、方程求根、数据插值、 数据拟合、数值积分、微分方程求解 von Neumann 引例: 圆内接正多边形边长计算Pi方法(P.42177) 评价算法的主要指标: 速度和精度 求未知数据的迭代计算技术: 初始猜测数据、迭代计算格式、迭代序列的收敛性分析、计算复杂性分析，…… 3/16 n L error 192 3.1414524 1.4e-004 384 3.1415576 3.5e-005 3.1415926 4.6e-010 通信卫星覆盖地球面积 数学模型 实际问题 获取数据 数值方法、程序 数据结果 将地球考虑成一 个球体, 设R为地 球半径,h为卫星 高度,D为覆盖面 在切痕平面上的 投影(积分区域) 4/16 参考P.190 模型误差: 建立数学模型时所引起的误差； 误差分类： 舍入误差:计算机表示的数的位数有限，通常用四舍五入的办法取近似值，由此引起的误差. 截断误差:求解数学模型时，用简单代替复杂,或者用有限过程代替无限过程所引起的误差 观测误差:测量工具的限制或在数据的获取时随机因素所引起的物理量的误差； 5/16 假设某一数据的准确值为 x*,其近似值 为 x，则称 而称 为 x 的相对误差 误差的有关概念 e(x)= x - x* 为 x 的绝对误差 6/16 如果存在一个适当小的正数ε ，使得 则称ε为绝对误差限。 称εr为相对误差限。 如果存在一个适当小的正数εr ，使得 7/16 十进制浮点数表示 一台微机价格:￥3999.00, 浮点数表示:0.3999×104 地球半径: 6378137m, (6.378137e+006) 浮点数表示: 0.6378137×107 光速: 2+008 浮点数表示: 0.299792458×109 尾数部 阶码部 8/16 其中,正负号占2个位,尾数占52个位,阶码占10个位.对应十进制数字长15,阶码308 尾数部 阶码部 二进制数1.b1b2×2m ( – 4≤ m ≤3 )分布实验 二进制浮点数表示(IEEE754双精度) 9/16 有效数字概念: 取? 的有限位数如下(? ≈3.1415926) 取 x1 = 3，误差限不超过0.5; 取 x2 = 3.14,误差限不超过0.005 ; 若近似值 x 的绝对误差限是某一位上的半个 单位，该位到 x 的第一位非零数字一共有 n 位，则称近似值 x 有 n 位有效数字. 取 x3 = 3.1416,误差限不超过0.00005 ; 10/16 一个有n 位有效数字的数 绝对误差限满足： 相对误差限满足： 11/16 例1 已知 的十进制浮点数第一位是5，要使近似值的相对误差限小于0.1%，问浮点数的有效数字的位数至少应该为多少？ 解: a1=5,利用不等式 所以,浮点数的有效数字位数至少应取3位。 取n≥3,有 |er(x)|≤10-3 12/16 13/16 例2.圆面积计算的误差估计 圆面积计算公式: 全微分近似: 取 r = 50 cm, 则有 cm ≈2×1%=2% ≈150 cm2, 反问题:利用 估计 1.一元函数 y=f(x)误差分析( 准确值 y*=f(x*) ) 由Taylor 公式 同理: 所以 14/16 反问题:估计 2.多元函数 z = f(x1,x2,···,xn)误差分析 (1) (3) (2) 15/16 数据误差对算术运算影响 例3.二次方程 x2 – 16 x + 1 = 0, 取 求 使具有4位有效数 解:直接计算 x1≈8 – 7.937 = 0.063 修改算法 4位有效数 16/16 计算出的x1 具有两位有效数 参考文献 [1]李庆扬 关治 白峰杉, 数值计算原理(清华) [2]蔡大用 白峰杉, 现代科学计算 [3]蔡大用, 数值分析与实验学习指导 [4]孙志忠,计算方法典型例题分析 [5]车刚明等, 数值分析典