1.《计算方法》-误差

《计算方法》教案

（第一章 误差）

选用教材：普通高等教育“十一五”国家级规划教材

《计算方法引论》（第三版）

徐箤薇 孙绳武 编著

主讲老师：刘鸣放

2010年3月于河南大学

一． 基本内容提要

1. 误差的来源

2. 浮点数、误差、误差限和有效数字

3. 相对误差和相对误差限

4. 误差的传播

5. 在近似计算中需要注意的一些问题

二． 教学目的和要求

1. 熟练掌握绝对误差、绝对误差限、相对误差、相对误差限和有效数字的概念及其相互关系；

2. 了解误差的来源以及误差传播的情况，掌握在基本算术运算中误差传播后对运算结果误差限的计算方法和函数求值中的误差估计；

3. 理解并掌握几种减少误差避免错误结果应采取的措施，了解选用数值稳定的算法的重要性。

三． 教学重点

1. 绝对误差、绝对误差限、相对误差、相对误差限和有效数字的概念及其相互关系，误差传播，减少误差避免错误结果应采取的措施。

四． 教学难点

1. 误差传播； 2. 数值稳定算法的选用。

五． 课程类型

新知识理论课；

六． 教学方法

结合课堂提问，以讲授为主。

七．教学过程如下：

Introduction

1.《计算方法》课程介绍

计算方法是用数值的方法研究研究科学与工程中的计算问题；它的内容主要包括：近似值的计算和误差估计两个方面；主要工具：计算机；地位：这门课已成为工科各专业，特别是计算机科学与技术、土木工程、机械、数学等专业的必修基础课。

2. 发展状况

几十年来，计算方法效率的提高是与计算机速度的提高几乎同步地、同比例地前进的。这里简述一下国家重点基础研究计划项目（简称973项目）“大规模科学计算研究”（1999-2004）的主要内容，可以帮助同学们了解我国科学计算界所关心的问题。此项目由石钟慈院士等人为首组织，集中了我国计算数学、计算物理、计算力学、计算机、以及材料、环境能源等领域60多名专家，跨学科，跨部门通力合作研究以下几个方面的主要内容：

（1） 复杂流体的高精度计算，含天气预报数值模拟研究；

（2） 新材料的物理性质机理多尺度计算研究，含超导、超硬度合金等问题的计算研究；

（3） 地质油藏模拟与波动问题及其反问题计算研究；

（4） 基础计算方法的理论创新与发展；

（5） 大规模计算软件系统的基础理论和实施。

计算规模是解决百万级节点问题（即求解百万至千万个未知数的方程组）。项目的目标是在几个重大科技难题的计算研究中取得突破性进展，并在国际科学计算的学科前沿取得重要的一席之地。

3. 我国计算科学主要专家简介

冯康（1920-1993）院士，早年学习物理，后去苏联研究拓扑学和函数论，有深厚的数学和物理功底。1957年受命组建国家计算中心，为推动我国计算科学进入国际前列，培养一批批优秀人才，做出了不朽的贡献。冯院士一生中作出两项国际公认的重大创造，早在20世纪60年代，独立提出变分差分格式（即有限元），用于计算水坝很成功，1964年在国际上首次证明有限元的收敛性（中文发表），比M.Zlamal早四年，接着10年文化大革命，冯康的工作几乎被淹没，到70年代末才被国际公认为有限元理论的开拓者。1984年又开创了Hamilton系统的辛几何算法。

周毓麟（1923- ）院士，著名的拓扑学、偏微分方程与计算科学家。1945年毕业于大同大学数学系，后在北京大学任教，1953年留学莫斯科大学研究偏微分方程，是我国最早用先验估计和拓扑方法研究偏微分方程的专家。1957年回国在北大任教，1960年调国防科工委从事核武器与数值模拟的研究，为“两弹一星”的成功做出了默默的贡献，这20年的工作鲜为人知，仍公开发表百余篇论文，他对差

分法首次建立离散Sobolev空间的嵌入理论，并用于偏微分方程的研究，建立了全新体系。

石钟慈（1933- ）院士，早年在浙江大学和复旦大学学习基础数学，后留学苏联学习计算数学，在中国科学院成为冯康的得力助手、合作者和接班人。与冯康合作研究弹性组合结构的数学理论，获得国家自然科学奖。20世纪80年代初，是我国最早赴德国洪堡基金会研究的学者之一。在对四阶板问题的非协调元的研究中，首次提出了非协调元收敛的F-E-M检验准则，成为以后一系列研究的基石，也培养了一批批优秀的人才。今年来对瀑布式网格法研究中，证明了最佳收敛性，并引发了一系列工作。1993年冯康院士去世后，石钟慈院士成为我国计算数学的领头人，为继续推动我国科学计算进入世界强国做出了贡献。

我国计算数学界的院士还有林群和催俊芝。

4. 河南省计算数学发展状况介绍

5. 河南大学计算数学发展状况介绍

6. 课堂纪律要求，作业收交安排和答疑办法。

第一章 误 差

1.1误差的来源

用数学作为工具解决实际问题过程：

实际问题→数学模型→数值计算方法

→程序设计→上机计算结果。

从上述过程看，影响计算精度的误差可分为两类：一类是“过失误差”，人为造成的，可以避免；另一类是“非过失误差”，无法避免。按来源的不同，分为下面几种：

1. 模型误差：数学模型与实际问题之间出现的不可避免的误差。

将实际问题转化为数学问题，即建立数学模型时，对被描述的实际问题进行了抽象和简化，忽略了一些次要因素，数学模型只是客观世界的一种近似描述，之间存在一定的差别。属于“非过失误差”。

例如，用

212sgt

描述自由落体规律，就是一个数学模型，此模型建立时忽略了空气阻力等因素。如果用()st表示真正的运动规律，则模型误差为

21()2stgt。

2. 观测误差：建立模型和数值计算过程中，通常用到一些观测数据，由于仪器设备精度的限制，观测值和实际值之间的误差称为观测误差，也叫数据误差。

3. 截断误差：在计算中常遇到只有通过无限过程才能得到的结

果，但在实际计算中只能用有限过程来计算，于是产生了有限过程代替无限过程的误差，这种误差称为截断误差，也叫方法误差。

本课程主要研究该误差。

例如，指数函数在0x点有Taylor展式

212!!nxxxexn (||1)x

实际计算时只能取前面的有限项（例如n项）

10!knkxk

截断误差 100!!kknkkxxkk 。

4. 舍入误差：数值计算过程中遇到的数据可能位数很多，也可能是无穷小数，但在计算时只能对有限位数进行运算，往往要进行四舍五入，这样产生的误差称为舍入误差。

例如，12,,3等，在计算时只能取有限位数进行运算，要进行四舍五入。

少量舍入误差微不足道，但在计算机上完成百千万次运算后，舍入误差的积累有时是很大的。

前两种误差是客观存在的，后面两种是计算方法所引起的，本课程所研究的内容只涉及后两种误差。

1.2 浮点数、误差、误差限和有限数字

一． 浮点数

1. 浮点数：（1）浮点数属于有理数集的某特定子集，为该特定子集中数的数字表示；（2）浮点数可以用来近似表示任意一个实数，这种表示方法类似于基数为10的科学计数法；（3）如何一个浮点数x均可表示为

120.,JJtx LJU

其中，叫做这个数x的基（10进制中10，二进制中2），J是阶，取整数，是尾数，由t位小数构成，t又称精度，01(1,2,)iit。

所谓浮点数就是小数点在逻辑上是不固定的，而定点数只能表示小数点固定的数值，具用浮点数或定点数表示某哪一种数要看用户赋予了这个数的意义是什么。

浮点计算是指浮点数参与的运算，这种运算通常伴随着无法精确表示而进行的近似或舍入。

2. 浮点数的规格化（normalize）

浮点数在计算机中的表示是基于科学计数法（Scientific

Notation）.例如，32676用科学计数法可写成：43.267610，3.2676称为尾数（Mantissa \ mæn'tisə\, 或者叫Significand）.4称为指数（Exponent）,10为基数（Radix）.浮点数在计算机中的表示于此类似，只不过基数是2而不是10.例如:

021717.0100.1710,

类似地

052217(10001)2(0.10001)2，

可以看出，每个浮点数的表示都不唯一，这样给计算机处理数据增加了复杂性。为了解决这个问题，规定尾数部分的最高位必须是1，也就是说尾数必须以0.1开头，只对指数做相应的调整，这称为正规化，也叫规格化。

注：二进制：除2取余！例如，100用二进制表示（%表示求余数，/表示整数除法，忽略余数）：

100%2=0 \二进制最后一位\， 100/2=50;

50%2=0 \二进制倒数第二位\, 50/2=25;

25%2=1 \二进制倒数第三位\， 25/2=12;(忽略余数)

12%2=0 \二进制倒数第四位\， 12/2=6;

6%2=0 \二进制倒数第五位\， 6/2=3;

3%2=1 \二进制倒数第六位\， 3/2=1;

1%2=1 \二进制倒数第七位\， 结束。

100用二进制表示为

2561100100121212。

二． 误差、误差限和有效数字

1.误差：数x的一个近似值x与准确值x的差称为误差。用e来表

示，即

exx，

误差可证可负，近似值大于准确值，误差为正，叫“强近似”；近似值小于准确值，误差为负，叫“弱近似”。

2.误差限：误差绝对值的上限。用表示，即

||||exx

xxx

因此，可以用误差限表示近似值x的精确度

xx.

例. 用有毫米刻度的尺子测量桌子长度，1235xmm，是实际长度的一个近似值，由米尺的精度知，误差不会超过半个毫米，则有

1|||1235|2xxx

即

1234.51235.5x，

写成

(12350.5)xmm.

对于一个近似值，除了用误差表示精确程度外，还希望这个近似值本身就能表示出它的准确程度。于是引入有效数字的概念。

引例，31.732050808x，下面通过四舍五入取近似值：

取3位： 1.73x， ()0.005x；

取5位： 1.7321x， ()0.00005x，

这种近似值取法的特点是它们的误差都不超过末位数字的半个单位，即

241|31.73|10,21|31.7321|10.2

3.有效数字：如果近似值x的误差限是某位上的半个单位时，称近似值“准确”到这一位，且该位直到x的第一位非零数字一共有n位，则称x有n位有效数字。