数值分析1详解.ppt

*/18 数值计算的背景和内容 关于计算误差讨论 浮点数与有效数字 算术运算的误差估计 《数值分析》1 ? ? ? ? 具划时代意义的一篇论文：《高阶矩阵的数值求逆》 (冯·诺依曼 1947年) 数值计算 —— 以计算机为工具, 求解数学问题的方法. 1958年, 前苏联载人飞船 1969年, 美国Apollo 登月 1994年, 美国GPS运行 冯·诺依曼 图灵 现代数值计算特点: 应用范围广、问题规模大。对计算问题和方法研究，几乎涉及所有的数学分支。 中国北斗系统 本课程要学习和研究的内容: 求解非线性方程组问题的数值计算方法; 求解大规模线性方程组问题的数值计算方法; 解决数据插值问题的数值计算方法; 解决数据拟合问题的数值计算方法; 求各类积分问题的数值计算方法; 求解微分方程问题的数值方法简介. 掌握方法的数学思想、了解各类方法的应用 通信卫星信号覆盖率 将地球设为球体, R为半径, H为卫星高度, 单颗卫星信号覆盖面 三颗同步卫星信号在北极和南极有信号盲区 北极信号盲区模拟图 数学模型 数值方法 数值模拟 结果分析 引例: 圆内接正多边形边长计算Pi方法 评价算法的主要指标: 速度和精度 数值计算中常用的迭代技术: 初始数据获取、计算格式构造、迭代序列的收敛性分析、计算复杂性分析，…… n L error 192 3.1414524 1.4e-004 384 3.1415576 3.5e-005 3.1415926 4.6e-010 模型误差: 建立数学模型时所引起的误差； 误差分类： 舍入误差:计算机表示的数的位数有限，通常用四舍五入的办法取近似值，由此引起的误差. 截断误差:求解数学模型时，用简单代替复杂,或者用有限过程代替无限过程所引起的误差 观测误差:测量工具的限制或在数据的获取时随机因素所引起的物理量的误差； 假设某一数据的准确值为 x*,其近似值 为 x，则称 而称 为 x 的相对误差 误差的有关概念 e(x)= x - x* 为 x 的绝对误差 如果存在一个适当小的正数ε ，使得 则称ε为绝对误差限。 称εr为相对误差限。 如果存在一个适当小的正数εr ，使得 十进制浮点数表示 一台微机价格:￥3999.00, 浮点数表示:0.3999×104 地球半径: 6378137m, (6.378137e+006) 浮点数表示: 0.6378137×107 光速: 2+008 浮点数表示: 0.299792458×109 尾数部 阶码部 其中,正负号占2个位,尾数占52个位,阶码占10个位.对应十进制数字长15,阶码308 尾数部 阶码部 二进制数1.b1b2×2m ( – 4≤ m ≤3 )分布实验 二进制浮点数表示(IEEE754双精度) 有效数字概念: 取? 的有限位数如下(? ≈3.1415926) 取 x1 = 3，误差限不超过0.5; 取 x2 = 3.14,误差限不超过0.005 ; 若近似值 x 的绝对误差限是某一位上的半个 单位，该位到 x 的第一位非零数字一共有 n 位，则称近似值 x 有 n 位有效数字. 取 x3 = 3.1416,误差限不超过0.00005 ; 一个有n 位有效数字的数 绝对误差限满足： 相对误差限满足： 例1 已知 的十进制浮点数第一位是5，要使近似值的相对误差限小于0.1%，问浮点数的有效数字的位数至少应该为多少？ 解: a1=5,利用不等式 所以,浮点数的有效数字位数至少应取3位。 取n≥3,有 |er(x)|≤10-3