E-mail： jinxz865@nenu.edu.cn 授课: 60 学分：3.

Presentation on theme: "E-mail： jinxz865@nenu.edu.cn 授课: 60 学分：3."— Presentation transcript:

1 E-mail： jinxz865@nenu.edu.cn
授课: 60 学分：3

2 为什么要开设这个课呢？ 在数学发展中，理论和计算是紧密联系的。现代计算机的出现为大规模的数值计算创造了条件，集中而系统的研究适用于计算机的数值方法变得十分迫切和必要。数值计算方法正是在大量的数值计算实践和理论分析工作的基础上发展起来的，它不仅仅是一些数值方法的简单积累，而且揭示了包含在多种多样的数值方法之间的相同的结构和统一的原理。数值算法是进行科学计算必不可缺少的起码常识；更为重要的是通过对它们的讨论，能够使人们掌握设计数值算法的基本方法和一般原理，为在计算机上解决科学计算问题打下基础。因此，计算方法已经成为理工科大学生的必修课程。

3 1. 认识建立算法和对每个算法进行理论分析是基本 任务，主动适应“公式多”的特点； 2. 注重各章建立算法的问题的提法，搞清问题的基
如何进行学习？ 1. 认识建立算法和对每个算法进行理论分析是基本 任务，主动适应“公式多”的特点； 2. 注重各章建立算法的问题的提法，搞清问题的基 本提法，逐步深入； 3. 理解每个算法建立的数学背景，数学原理和基本 线索，对最基本的算法要非常熟悉； 4. 认真进行数值计算的训练，学习各章算法完全是 为用于实际计算，必须真会算。

4 科学素质：拓宽对21世纪科学的了解； 加深对数学思想的理解； 培养用数学思考世界的习惯 数学能力：数学知识的运用能力； 对专业中问题建立数学求解方法与 实际计算能力 应用问题中数学创造性能力 计算知识：常用算法的数学理论； 在“误差、存贮、速度”之下的实 际计算方法； 对结果的数值分析方法

5 几点要求 • 记好课堂笔记 • 保证课堂纪律 • 按时上课，不迟到早退 • 按时完成作业 test

6 数值分析讲述的基本内容 如何把数学模型归结为数值问题 如何制定快速的算法 如何估计一个给定算法的精度 分析误差在计算过程中的积累和传播
如何构造精度更高的算法 如何使算法较少的占用存储量 如何分析算法的优缺点

7 本课程的基本要求 掌握数值方法的基本原理 掌握常用的科学与工程计算的基本方法 能用所学方法在计算机上算出正确结果

8 第一章计算方法与误差 本章内容 §1 引言 §2 误差的来源及分类 §3 误差的度量 §4 误差的传播 §5 减少运算误差的原则 小结

9 第一章计算方法与误差 要求掌握的内容 概念 包括有效数字、绝对误差、绝对误差限、相对误差、相对误差限等 误差
截断误差、舍入误差的详细内容，误差种类等 分析运算误差的方法和减少运算误差的若干原则 test

10 1.1 引言 数值分析又称计算方法, 它是研究各种数 学问题的数值解法及其理论的一门学科。 数值分析的任务 实际问题 数学模型 数值计算方法
数值分析又称计算方法, 它是研究各种数 学问题的数值解法及其理论的一门学科。 数值分析的任务 实际问题 数学模型 数值计算方法 数值结果 上机计算 程序设计 根据数学模型提出求解的数值计算方法直到编出程序上机算出结果，这一过程便是数值分析研究的对象

11 数值方法解题的一般过程 1. 对于要解决的问题建立数学模型 2. 研究用于求解该数学问题近似解的算法和过程
3. 按照2进行计算，得到计算结果 建立数 学模型 转化为 数值公式 进行计算

12 计算作为工程技术研究方法 第一章计算方法与误差
数值计算以及计算机模拟（包括当前流行的虚拟现实的方法），已经是在工程技术研究和经济、社会科学中广泛应用的方法，带来巨大的经济效益 天气预报与亿次计算机 波音777的无纸设计与有限元 CT、核磁共振 计算流体力学与爆炸工程 能源问题与大型计算 计算作为工程技术研究方法

13 第一章计算方法与误差 计算方法课程主要讨论如何构造求数学模型近似解的算法，讨论算法的数学原理、误差和复杂性，配合程序设计进行计算试验并分析试验结果。 与纯数学的理论方法不同，用数值计算方法所求出的结果一般不是解的精确值或者准确的解析表达式，而是所求真解的某些近似值或近似曲线。

14 例如方程 x2=2sinx，在区间(1,2)内有唯一根, 但找不出求根的解析式, 只能用数值计算方法求其近似解。有些数学问题虽有理论上的准确的公式解, 但不一定实用, 例如行列式解法的Cramer法则原则上可用来求解线性方程组,用这种方法解一个n元方程组,要算n+1个n阶行列式的值,总共需要n!(n-1)(n+1)次乘法,当n=20时,其乘除法运算次数约需1021次方,即使用每秒千亿次的计算机也得需要上百年,而用高斯（Guass）消去法约需2660次乘除法运算,并且愈大,相差就愈大。可见研究和选择好的算法是非常重要的。

15 算法(数值算法):是指有步骤地完成解数值问 题的过程。 数值算法的特点 • 目的性，条件和结论、输入和输出数据均要有明
• 目的性，条件和结论、输入和输出数据均要有明 确的规定与要求。 • 确定性，精确地给出每一步的操作(不一定都是运 算)定义, 不容许有歧义。 • 可执行性，算法中的每个操作都是可执行的 • 有穷性，在有限步内能够结束解题过程 计算机上的算法，按面向求解问题的不同， 分为数值算法和非数值算法。

16 第一章计算方法与误差 1.2 误差的来源及分类 早在中学我们就接触过误差的概念，如在做热力学实验中，从温度计上读出的温度是23.4度，就不是一个精确的值，而是含有误差的近似值。事实上，误差在我们的日常生活中无处不在，无处不有。如量体裁衣，量与裁的结果都不是精确无误的，都含有误差。

17 第一章计算方法与误差 在用数值方法解题过程中可能产生的误差归纳起来有如下几类： 1. 模型误差 2. 观测误差 3. 截断误差
1. 模型误差 2. 观测误差 3. 截断误差 4. 舍入误差

18 1. 模型误差 用数学方法解决一个具体的实际问题，首先要建立数学模型，这就要对实际问题进行抽象、简化，因而数学模型本身总含有误差，这种误差叫做模型误差 数学模型是指那些利用数学语言模拟现实而建立起来的有关量的描述 数学模型的准确解与实际问题的真解不同 实际问题的真解 数学模型的真解 为简化模型忽略次要因素 定理在特定条件下建立与实际条件有别

19 2. 观测误差 在数学模型中通常包含各种各样的参变量，如温度、长度、电压等，这些参数往往是通过观测得到的，因此也带来了误差，这种误差叫观测误差 数学模型中的参数和原始数据，是由观测和试验得到的 由于测量工具的精度、观测方法或客观条件的限制,使数据含有测量误差,这类误差叫做观测误差或数据误差 根据实际情况可以得到误差上下界 数值方法中需要了解观测误差,以便选择合理的数值方法与之适应

20 3. 截断误差 精确公式用近似公式代替时,所产生的误差叫截断误差 例如, 函数f(x)用泰勒(Taylor)多项式
近似代替，则数值方法的截断误差是 （介于0与x之间） 截断误差的大小直接影响计算结果的精度和计算 工作量，是数值计算中必须考虑的一类误差

21 4. 舍入误差 在数值计算中只能对有限位字长的数值进行运算 需要对参数、中间结果、最终结果作有限位字长的处理工作，这种处理工作称作舍入处理 用有限位数字代替精确数，这种误差叫做舍入误差，是数值计算中必须考虑的一类误差

22 第一章计算方法与误差 例如在计算时用3.14159近似代替， 产生的误差R= -3.14159=0.0000026… 就是舍入误差。
上述种种误差都会影响计算结果的准确性，因此需要了解与研究误差，在数值计算中将着重研究截断误差、舍入误差，并对它们的传播与积累作出分析

23 1.3 误差的度量 1.3.1 绝对误差和绝对误差限 定义1.1 设精确值x的近似值为x* ，称差 e(x*) =x-x*
1.3 误差的度量 1.3.1 绝对误差和绝对误差限 定义1.1 设精确值x的近似值为x* ，称差 e(x*) =x-x* 为近似值x*的绝对误差，简称误差。e(x*)又记为e* 当e*>0时，x*称为弱近似值，当e*<0时，x*称为强近似值。 |e*|越小， x*的精度越高。 由于精确值一般是未知的,因而e*不能求出来, 但可以根据测量误差或计算情况设法估计出它的取值范围，即误差绝对值的一个上界或称误差限。

24 则称为近似值的绝对误差限，简称误差限或精度。 来表示近似值的精度或准确值x所在的范围。
1.3 误差的度量 定义1.2 设存在一个正数，使 则称为近似值的绝对误差限，简称误差限或精度。 实际应用中经常使用这个量来衡量误差限, 这就是说, 如果近似数 的误差限为 , 则 表明准确值 x 必落在  上, 常采用下面的写法 来表示近似值的精度或准确值x所在的范围。

25 可见，绝对误差限*不是唯一的，但*越小越好
1.3 误差的度量 例1 设x == … 近似值x* =3.14,它的绝 对误差是 …，有 ‌‌ x-x*= … 0.002=0.210-2 例2 又近似值x* =3.1416，它的绝对误差是 …，有 x-x*= …  =0.810-5 例3 而近似值x* =3.1415，它的绝对误差是 …，有 x-x*= … 0.0001=0.110-3 可见，绝对误差限*不是唯一的，但*越小越好 a-ε a+ε a A

26 1.3.2 相对误差和相对误差限 只用绝对误差还不能说明数的近似程度,例如甲打字每100个错一个,乙打字每1000个错一个,他们的误差都是错一个,但显然乙要准确些,这就启发我们除了要看绝对误差外,还必须顾及量的本身。 定义1.3 绝对误差与精确值 x 的比值 称为相对误差。 简记为

27 1.3.2 相对误差和相对误差限 定义1.4 设存在一个正数 ，使
相对误差越小,精度就越高,实际计算时,x通常是不知道的,因此可用下列公式计算相对误差 定义1.4 设存在一个正数 ，使 则称 为近似值 的相对误差限。 简记为

28 1.3.2 相对误差和相对误差限 例4. 甲打字每100个错一个，乙打字每1000个 错一个，求其相对误差 解： 根椐定义:甲打字时的相对误差 乙打字时的相对误差

29 1.3.3 有效数字 定义1.5 设x的近似值 其中 是0到9之间的任一个数,但 n是正整数, m是整数,若
有效数字 定义1.5 设x的近似值 其中 是0到9之间的任一个数,但 n是正整数, m是整数,若 则称 为 x 的具有n位有效数字的近似值， 准确到第n位， 是 的有效数字。

30 1.3.3 有效数字 例5. 3.142作为π的近似值时有几位有效数字 解： 3.141592…= 0.3141592…×
有效数字 例 作为π的近似值时有几位有效数字 解： …= …× 3.142 = × m = 1 |π |=| …× × | ＜ × ＜ = × m –n = 1–n = -3 所以 n = 4，具有4位有效数字

31 例6. 当取3.141作为的近似值时 -3.141= … 101 ≤ 101 <0.005=1/2  10-2 m-n=1-n=-2 所以n=3具有3位有效数字 推论 如果近似数x*的误差限是某一位的半个单 位,由该位到x*的第一位非零数字一共有n位， x*就有n位有效数字,也就是说准确到该位。

32 再如3.1416作为的近似值时   = … 101  ≤ 101 ≤ < <0.5  10-4 m-n=1-n=-4 所以 n=5 x*= 有5位有效数字

33 关于有效数字说明 ① 用四舍五入取准确值的前n位x*作为近似值,则 x*必有n位有效数字。如3.142作为 的近似值 有4位有效数字，而3.141为3位有效数字 ② 有效数字相同的两个近似数，绝对误差不一定 相同。例如，设x1*=12345,设x2*=12.345,两者 均有5位有效数字但绝对误差不一样 x- x1* =x  ≤ 0.5= 1/2  100 x- x2* =x ≤0.0005=1/210-3 ③ 把任何数乘以10p(p=0,1,…)不影响有效位数 ④ 准确值具有无穷多位有效数字,如三角形面积 S=1/2ah=0.5ah 因为0.5是真值,没有误差 *=0,因此n,准确值具有无穷位有效数字

34 1.3.4 有效数字与相对误差 定理1.1 若近似数x*=0.x1x2…xn10m具有 n 位 有效数字，则其相对误差
∴ x* ≥x110 m-1 又 ∵ x*具有n位有效数字,则x- x*≤1/210 m - n ∴ 即

35 1.3.4 有效数字与相对误差 ∴有效数字越多，相对误差就越小 例7 取3.14作为的四舍五入的近似值时，求其 相对误差
一般应用中可以取r*=1/2x1 10-(n-1),n越大,r*越小, ∴有效数字越多，相对误差就越小 例7 取3.14作为的四舍五入的近似值时，求其 相对误差 解：3.14=0.314  x1=3 m=1 ∵ 四舍五入的近似值,其各位都是有效数字 ∴ n=3 r*=1/2x1 10-(n-1)=1/2*3 10-2=17%

36 1.3.4 有效数字与相对误差 例8 已知近似数x*有两位有效数字，试求其相 对误差限
解：已知 n=2 代入公式 r*=1/2x1 10-(n-1)得 r*=1/2x1 10-1 x*的第一位有效数字x1没有给出，可进行如下讨论：当 x1=1 r*=1/2x1 10-1=1/2*1 10-1=5% x1=9 r*=1/2x1 10-1=1/2*9 10-1=0.56% 取 x1=1 时相对误差为最大，即 5%

37 由有效数字定义可知, x*具有n位有效数字。证毕
1.3.4 有效数字与相对误差 定理1.2 若近似数x*=  0.x1x2…xn10m的相对误差 则该近似数具有n位有效数字 证: ∵ x* =  0.x1x2…xn10m ∴  x* ≤ (x1+1) 10m-1 由有效数字定义可知, x*具有n位有效数字。证毕

38 例9 已知近似数x*的相对误差限为0.3%，问x*
有几位有效数字？ 解：由 得 ⅰ 当x1=1时,310-3=1/410-(n-1)1210-3=10-(n-1) 上式两边取以10为底的对数得 lg22+lg3+(-3)=-n+1 ∵lg2= lg3=0.4771 2 =-n ∴ n=2.9209 ⅱ 当x1=9时,310-3=1/2010-(n-1) 610-3=10-n lg2+lg3+(-3)=-n ∴ n=2.2219 ∴ x*至少有3位有效数字