[理学]计算法方法 第一章 预备知识_深色.pptVIP

第一章 预备知识 聊城大学计算机学院 孙 剑 本章主要内容： 1.1 计算方法引论 1.2 误差 1.3 算法的数值稳定性 1.4 向量和矩阵的范数 1.2 科学计算中误差的来源 1.2.1 浮点数及其运算特点 1.2.2 误差的来源与分类 1.2.3 误差的基本估计方法 1.2.1 浮点数及其运算特点 纯数学中a + b = a 10位十进制计算机中：1010+4=1010 1. 规格化浮点数： x =±0.d1d2…dt×βn 1.2.1 浮点数及其运算特点 十进制情形，β为10，d1,d2,…,dt为0到9中任一数字。 数0在计算机中尾数为0，阶码任意。 一台计算机能表示的浮点数的全体，记作F。实数x在计算机中用F中最接近x的一个浮点数表示。 例1-1 t = 4, β= 10，即4位十进制计算机中： 1.2.1 浮点数及其运算特点 2.浮点数的运算过程: (1) 规格化两个加数 (2)两数相加前先对阶，统一为较大的阶。 (3)结果规格化 例1-2 设t=4, β=10, 求 x + y 解： x =0.3127×10 –6, y=0.4153×10 –4 则x + y = 0.3127×10 –6+ 0.4153×10 –4 1.2.1 浮点数及其运算特点 例1-3 设在5位十进制计算机上，t=5, β=10, x = 0.37569×10 4, y = 0.96331×10 –5 例1-4 10位十进制计算机中计算：1010+4=？ 1010+4 = 0.100 000 000 0×1011 + 0.400 000 000 0×101 = 0.100 000 000 0×1011 + 0.000 000 000 0×1011 = 0.100 000 000 0×1011 =1010 1.2.2 误差的来源与分类 用计算机解决科学计算问题经历以下过程： 实际问题 建立数学模型 构造数值计算方法 程序设计 上机计算结果 分析结果 。 1.2.2 误差的来源与分类 误差的来源主要有四类： 1.模型误差---客观量的准确值与数学模型的准确解的差。 2.观测误差---由观测数据而产生的误差。 3.截断误差(方法误差)---数学模型的准确解与利用近似计算方法得到的解之差。 4.舍入误差---由于计算机字长有限，将数据进行舍入而产生。 1.2.2 误差的来源与分类 1.2.3 误差的基本估计方法 1.2.3.1 绝对误差和绝对误差限 1.2.3.2 相对误差和相对误差限 1.2.3.3 有效数字 1.2.3.4 算术运算的误差 1.2.3.1 绝对误差和绝对误差限 定义1-1：设某准确值x*近似值为x ， 则称ε(x) = x* – x 为x的绝对误差。 在同一量x*的不同近似值中，|ε(x)|越小，近似值x 的精确度越高。 准确值x*一般无法得到，通常用|ε(x)|适当小的上界ξ来度量绝对误差，称ξ为近似值x的绝对误差限（简称误差限）。 绝对误差限ξ不唯一！ 1.2.3.1 绝对误差和绝对误差限 例 试估计误差限。 解： |ε(x)|=|x*–x|=0.000333….0.0005 x的绝对误差限为0.0005 常用记法： x* = x±ξ 表示 x-ξ≤ x*≤x+ξ ？某商品标注重量为 27±0.5g，实际重量是多少？ ！实际重量在26.5g到27.5g之间。 1.2.3.2 相对误差和相对误差限 定义1-2：设准确值x*的近似值为x,（x*≠0） 则称 为x的相对误差。 在不同近似值中，|εr (x)|越小，x的精确度越高 1.2.3.2 相对误差和相对误差限 常用计算公式（ξ=|x*-x| ，绝对误差限） （为什么？） 当ξ已知时，|εr (x)|= 当η已知时，有|ε(x)|=|εr (x)| |x|≤η|x| 1.2.3.2 相对误差和相对误差限 例1-7：设 x、y的近似值哪一个精度高些? 解：x*的近似值x=1, x的绝对误差限ε (x)=0.5, 相对误差限η(x)= ε (x) / x= 0.5/1=0.5 y*的近似值y=10000, y的绝对误差限ε (y)=5,