计算方法课件_易大义主编
合集下载
计算方法第一章数值计算方法.ppt
x1
a22b1
a12b2 D
S4 输出计算的结果 x1, x2
x2
a11b2
a21b1 D
开始
输入
a11, a12 , a21, a22 ,b1,b2
D=a11a22-a12a21
Yes D=0
No
x1 (b1a22 b2a12 ) / D x2 (b2a11 b1a21) / D
输出无解信息
…
…
第一章计算方法与误差
本章内容
§1 引言 §2 误差的来源及分类 §3 误差的度量 §4 误差的传播 §5 减少运算误差的原则
小结
第一章计算方法与误差
要求掌握的内容
概念 包括有效数字、绝对误差、绝对误差限、 相对误差、相对误差限等
误差 截断误差、舍入误差的详细内容，误差种 类等
分析运算误差的方法和减少运算误差的若 干原则
常用的两种复杂性有：计算时间复杂性和空间复杂性。
二、算法的优劣
➢ 计算量小 例：用行列式解法求解线性方程组:
n阶方程组，要计算n + 1个n阶行列式的值，
总共需要做n! (n - 1) (n + 1) 次乘法运算。
n=20 需要运 算多少次？
n=100?
计算量大小是衡量算法优劣的一项重要标准。
在估计计算量时，我们将区分主次抓住计算过程中费时较多的 环节。比如，由于加减操作的机器时间比乘除少得多，对和式
例：求解二元一次联立方程组
aa1211xx11
a12 x2 a22 x2
b1 b2
用行列式解法：首先判别
D a11a22 a21a12
是否为零，存在两种可能：
（1）如果 D 0，则令计算机计算
计算方法易大义主编 ppt课件
29
§1.4 有效数字及其与误差的关系
例3： 设x*＝0.0270是 某 数x经 “ 四 舍 五 入 ” 所 得 ， 则 误 差e(x*)不 超 过x*末 位 的 的 半 个 单 位 ，：即 x x * 1 104 2 又 x* 101(0.270) 故 该 不 等 式 又 可 写 为 x x * 1 1013 2 由 有 效 数 字 定 义 可, x知* 有3位 有 效 数 字 ,分 别 是2，7，0。
解决科学技术和工程问题的步骤:
实际问题 数学问题 提供计算方法 程序设计 上解的数学模 型简化成一系列算术运算和逻辑运算, 以便 在计算机上求解, 并对算法的稳定性、收敛 性和误差进行分析。
计算方法——09计11、61
17
§1.1 引言
简单地说，就是研究如何用计算机有效地 解决一个数学问题。
值班军官对连长: 根据营长的命令, 明晚8点哈雷彗 星将在操场上空出现。如果下雨的话, 就让士兵穿 着野战服列队前往礼堂, 这一罕见的现象将在那里 出现。
计算方法——09计11、61
26
§1.2 误差的种类及其来源
连长对排长: 根据营长的命令, 明晚 8 点, 非凡的哈 雷彗星将身穿野战服在礼堂中出现。如果操场上 下雨, 营长将下达另一个命令, 这种命令每隔 76 年 才会出现一次。
若x*x 110mn 2
(P101.4.2)
则 说x*具 有 n位 有 效 数 字 ，a分 1,a2别 ,,是 an
若n＝p， 则 称 x*为 有 效 数 。
➢ 一定要从规格化后的数来判断其位数
➢ 有效位数与第一个非 0 项后的数字个数是不 一致的。 四舍五入所得到的数是一致的。
计算方法——09计11、61
计算方法
1 3 1 2 ≈ n + n − n. 3 3
►高斯消元法(Gauss)： 运算量(乘除)
取 n = 20,
Gauss: 3060次；
Cramer: ( 20+ 1) ≈ 5. × 1019 ≈ 16200年(108 次 / 秒) ! 1
2011-2-25 5
本课程的任务： 本课程的任务：
● 建立各种数学问题的数值计算算法的方法和理论。 通俗地讲，就是为各种实际问题提供有效的数值近似 解方法。 ● 提供在计算机上实际可行的、理论可靠的、计算复杂 性好的各种常用算法。
绝对误差限不超过该近似数末位的半个单位。
2011-2-25
13
相对误差
先看两种产品的不合格率： 8/1000=0.8% ； 12/2000=0.6% （相对误差的值） 定义：设 x* 为精确值，x 为近似值，相对误差为
e x * −x er = = . x* x*
当绝对误差 e 较小时，相对误差可写为
一、防止相近的两数相减
(会耗失许多有效数字,可以用数学公式化简后再做).
例1: 各有五位有效数字的23.034与22.993相减. 23.034-22.993=0.041 0.041只有两位有效数字,有效数字的耗失,说明准确度减小,因此, 在计算时需要加工计算公式,以免这种情况发生. 例2:当较大时,计算 x +1 − x
2011-2-25
11
误差的概念
(绝对)误差 设 x* 为精确值（或准确值），x 为 x* 的近似值，称 e = x*－ x 为近似值x的(绝对)误差。 (绝对)误差限 （ ε ） 如果精确值 x* 与近似值 x 的误差的绝对值不超过某个正 数 ε，即 | e |= | x*－ x |≤ ε。 于是，精确值也可表示为 x* = x ± ε， 或
►高斯消元法(Gauss)： 运算量(乘除)
取 n = 20,
Gauss: 3060次；
Cramer: ( 20+ 1) ≈ 5. × 1019 ≈ 16200年(108 次 / 秒) ! 1
2011-2-25 5
本课程的任务： 本课程的任务：
● 建立各种数学问题的数值计算算法的方法和理论。 通俗地讲，就是为各种实际问题提供有效的数值近似 解方法。 ● 提供在计算机上实际可行的、理论可靠的、计算复杂 性好的各种常用算法。
绝对误差限不超过该近似数末位的半个单位。
2011-2-25
13
相对误差
先看两种产品的不合格率： 8/1000=0.8% ； 12/2000=0.6% （相对误差的值） 定义：设 x* 为精确值，x 为近似值，相对误差为
e x * −x er = = . x* x*
当绝对误差 e 较小时，相对误差可写为
一、防止相近的两数相减
(会耗失许多有效数字,可以用数学公式化简后再做).
例1: 各有五位有效数字的23.034与22.993相减. 23.034-22.993=0.041 0.041只有两位有效数字,有效数字的耗失,说明准确度减小,因此, 在计算时需要加工计算公式,以免这种情况发生. 例2:当较大时,计算 x +1 − x
2011-2-25
11
误差的概念
(绝对)误差 设 x* 为精确值（或准确值），x 为 x* 的近似值，称 e = x*－ x 为近似值x的(绝对)误差。 (绝对)误差限 （ ε ） 如果精确值 x* 与近似值 x 的误差的绝对值不超过某个正 数 ε，即 | e |= | x*－ x |≤ ε。 于是，精确值也可表示为 x* = x ± ε， 或
计算方法_6.ppt
随
时
孕
间育
呈
S长期形变化饱和。
例如：反应过程中的产物浓期度、相变过期程
中新相的含量、植物的生长量、市场某产
品的销售量等与时间的关系都为S形曲线。
S形生长反映了一个动力学过程中“孕 育期”、“生长期”和“饱和期”的发展 规律，是材料研究中经常涉及到的模型。
第六章 非线性拟合问题
硬度，HV 再结晶体积分数，%
▪ 取对数分离变量，如： y axb , y aebx 取对数，分别得：ln y = ln a + b ln x、 ln y = ln a + b x 再通过变量替换，转换为线性模型
v=a+bu
第六章 非线性拟合问题
“非线性” 问题
但是在许多情况下，拟合模型过于复杂。直接采用 最小二乘法，“法方程组” (p.90, 4.56式) 可能是不可解 的，或者是没有相应的解析表达式的。
算模型参数：s、EF
300
200
rapidly solidified FeSi2 100 950°C 50MPa 30’ HUP
800°C 20hr annealed
00
100 200 300 400 500
Temperature, °C
问题： Fermi 积分中包含模型参数 Fermi 积分无解析解 不能给出最小二乘“法方程组”
xi b
x
2 i
g
yi
xi a
x
2 i
b
x
3 i
g
xi yi
x
2 i
a
x
3 i
b
x
4 i
g
x
2 i
yi
已证明这个关于 a、b、g 的线性方程组有唯一解。
计算方法 易大义 陈道琦
主要内容
插值法 数值逼近 数值积分 数值微分 解线性方程组 非线性方程（组）数值解法 矩阵特征值与特征向量的计算 常微分方程数值解法
基础知识
高等数学（一元微积分） 线性代数（多项式、矩阵、线性方程组）
第一章 数值计算引论
§1 数值分析研究对象
实际问题 程序设计
则乘法次数仅为n.
2.防止大数“吃掉”小数
当|a|远远大于|b|时，尽量避免a+b 。例如，假设 计算机 只能存放10位尾数的十进制数，则
108 0.04 109 0.1 109 0.00000000004 108
3.尽量避免相近数相减
例如，当x很大时，应
，
x1 x
ex / y 1 ex x ey
y
y2
er x / y er x er y 3. f x, y x y
ex y ex ey
er x
y
x
x
y
er x
x
y
y
er y
例1.3：测得某桌面的长a的近似值a*=120cm,宽b的
a
b
b * e(a*) a * e(b*)
| e(s*) || b* || e(a*) | | a* || e(b*) | 60 0.2 120 0.1 24cm2
相对误差限为
e(s*) 24
| er (s*) ||
s*
|
0.33%
120 60
§5 数值计算中应该注意的几个问题
1
101n
，则
x* 至少有n位有效数字。
§4 求函数值的误差估计
1. 一元函数情形
插值法 数值逼近 数值积分 数值微分 解线性方程组 非线性方程（组）数值解法 矩阵特征值与特征向量的计算 常微分方程数值解法
基础知识
高等数学（一元微积分） 线性代数（多项式、矩阵、线性方程组）
第一章 数值计算引论
§1 数值分析研究对象
实际问题 程序设计
则乘法次数仅为n.
2.防止大数“吃掉”小数
当|a|远远大于|b|时，尽量避免a+b 。例如，假设 计算机 只能存放10位尾数的十进制数，则
108 0.04 109 0.1 109 0.00000000004 108
3.尽量避免相近数相减
例如，当x很大时，应
，
x1 x
ex / y 1 ex x ey
y
y2
er x / y er x er y 3. f x, y x y
ex y ex ey
er x
y
x
x
y
er x
x
y
y
er y
例1.3：测得某桌面的长a的近似值a*=120cm,宽b的
a
b
b * e(a*) a * e(b*)
| e(s*) || b* || e(a*) | | a* || e(b*) | 60 0.2 120 0.1 24cm2
相对误差限为
e(s*) 24
| er (s*) ||
s*
|
0.33%
120 60
§5 数值计算中应该注意的几个问题
1
101n
，则
x* 至少有n位有效数字。
§4 求函数值的误差估计
1. 一元函数情形
计算方法第二章ppt
当方程组的系数矩阵为非奇异 矩阵（即满秩矩阵）时，高斯 消元法可求得唯一解。
列主元高斯消元法
列主元高斯消元法的 基本思想
在高斯消元法的基础上，每次选取列 中绝对值最大的元素作为主元进行消 元，以避免出现小主元导致的误差放 大问题。
列主元高斯消元法的 步骤
首先选取第一列中绝对值最大的元素 作为主元，通过行交换将其移到第一 行第一列位置，然后进行高斯消元。 在后续的消元过程中，每次均选取当 前列中绝对值最大的元素作为主元进 行消元。
100%
数值解法
通过计算机求解常微分方程的近 似解的方法，主要包括欧拉方法 和龙格-库塔方法等。
80%
离散化与步长
将连续的时间或空间域离散化， 取离散点上的函数值作为近似解 ，步长是相邻离散点间的距离。
欧拉方法
显式欧拉法
一种简单的数值解法，通过前 一步的函数值及其导数来推算 下一步的函数值。
隐式欧拉法
通过求解一个非线性方程来得 到下一步的函数值，具有较高 的精度和稳定性。
改进欧拉法
结合显式欧拉法和隐式欧拉法 的优点，提高算法的精度和效 率。
龙格-库塔方法
龙格-库塔法基本思想
自适应步长龙格-库塔法
通过多步计算并利用泰勒级数展开式， 得到更高精度的近似解。
根据误差估计自动调整步长，实现精 度和计算效率的动态平衡。
标准四阶龙格-库塔法
一种常用的高精度数值解法，具有局 部截断误差为$O(h^5)$的优点。
常微分方程数值解法误差分析
局部截断误差
数值解法在单步计算中所产生的误差，可以通过泰勒级数展开式进行估计。
全局误差
数值解法在整个计算过程中所产生的累积误差，与算法稳定性、步长选择等因素有关。
计算方法课件
x g(x)
（2） 构造迭代公式
xk1 g(xk ), (k 0,1,2,)
在有根区间【a,b】上取一点 x0 (初始近似根)作为方程的近似值，代 入上面公式右端，求得 x1 g(x0 ) ，在把 x1 作为预测值，得到
x2 g(x1), 如此反复进行下去，得到一个近似根的序列
| g(x) | q 1
则迭代公式 xk1 g(xk ) 对于任意初值 x0 [a,b] 均收敛于方程
x g(x) 在区间 [a,b] 上的惟一根 x*, 且有如下误差估计式
q
qk
| x * xk | 1 q | xk xk 1 | 1 q | x1 x0 |
二. 误差
1. 误差的来源 模型误差：由于数学建模过程中往往要忽略一些次要的因 素而产生的，它是实际问题的客观量与其理论数学模型 精确解之间的差。（不讨论）
观测误差：在数学模型中的数据如果是从观测得到的，由 此产生的误差叫做观测误差。（不讨论）
截断误差（方法误差）：数学模型的精确解与近似算法的 解之间的差成为该算法的截断误差。
断 f (x) 的符号，逐步将有根区间缩小，直至有根区间足够的小，便求 出满足精度要求的近似根。
a x*
b
a
x*
b
■
bk
ak
bk 1
ak1 2
ba 2k
令
xk
1 2
(ak
bk ),
则
lim
k
xk
x*
误差估计式： | x * xk | (bk ak ) / 2 (b a) / 2(k1)
（3） 若 | x1 x0 | , 则停止计算；否则 x0 x1, 转（2）。
（2） 构造迭代公式
xk1 g(xk ), (k 0,1,2,)
在有根区间【a,b】上取一点 x0 (初始近似根)作为方程的近似值，代 入上面公式右端，求得 x1 g(x0 ) ，在把 x1 作为预测值，得到
x2 g(x1), 如此反复进行下去，得到一个近似根的序列
| g(x) | q 1
则迭代公式 xk1 g(xk ) 对于任意初值 x0 [a,b] 均收敛于方程
x g(x) 在区间 [a,b] 上的惟一根 x*, 且有如下误差估计式
q
qk
| x * xk | 1 q | xk xk 1 | 1 q | x1 x0 |
二. 误差
1. 误差的来源 模型误差：由于数学建模过程中往往要忽略一些次要的因 素而产生的，它是实际问题的客观量与其理论数学模型 精确解之间的差。（不讨论）
观测误差：在数学模型中的数据如果是从观测得到的，由 此产生的误差叫做观测误差。（不讨论）
截断误差（方法误差）：数学模型的精确解与近似算法的 解之间的差成为该算法的截断误差。
断 f (x) 的符号，逐步将有根区间缩小，直至有根区间足够的小，便求 出满足精度要求的近似根。
a x*
b
a
x*
b
■
bk
ak
bk 1
ak1 2
ba 2k
令
xk
1 2
(ak
bk ),
则
lim
k
xk
x*
误差估计式： | x * xk | (bk ak ) / 2 (b a) / 2(k1)
（3） 若 | x1 x0 | , 则停止计算；否则 x0 x1, 转（2）。
计算方法上课用PPT课件
1. 方程组阶数n很大，例如n=20,计算机运算速度 1亿次/秒,用不好的方法,大约需算30多万年; 好方法不到一分钟。另外，有计算结果可靠性 问题。
2. 特征值定义 A x x ( x 0 ) A xx0(AI)x0 | AI|0
14
3. f ( x) 形式复杂时求根和求积分很困难。
4.线性微分方程易解， 如
“计算方法"研究对象与特点
“计算方法"是计算数学的一个分支，它研究用计算机求解数学问题的数
值计算方法及其软件实现.计算数学几乎与数学科学的一切
分支有联系，它利用数学领域的成果发展了新的更有效的
算法及其理论，反过来很多数学分支都需要探讨和研究适
用于计算机的数值方法.因此，"计算方法"内容十分广泛.但
实际问题 程序设计
数学问题 上机计算
提供计算方法 结果分析
12
基本的数学问题:
1.大型线性代数方程组Ax=b求解;
2.矩阵A的特征值和特征向量计算;
3.非线性方程 f ( x ) 0 求解(求根);
4.积分 b a
f
( x)dx计算;
5.常微分方程初值问题求解;
6.其它。
13
求精确解(值)一般非常困难。例如：
17
截断误差 在求解过程中，往往以近似替代， 化繁为简，这样产生的误差称为截断误差。
舍入误差 在计算机上运算时受机器字长的 限制，一般必须进行舍入，此时产生的误 差称为舍入误差。
18
3. 截断误差，如
sin xxx3 x5 ....,.. 3! 5!
7
数值计算方法或数值分析主要是研究如何 运用计算机去获得数学问题的数值解的理 论和方法.对那些在经典数学中,用解析方法 在理论上已作出解的存在,但要求出他的解 析解又十分困难,甚至是不可能的这类数学 问题,数值解法就显得不可缺少,同时有十分 有效.
2. 特征值定义 A x x ( x 0 ) A xx0(AI)x0 | AI|0
14
3. f ( x) 形式复杂时求根和求积分很困难。
4.线性微分方程易解， 如
“计算方法"研究对象与特点
“计算方法"是计算数学的一个分支，它研究用计算机求解数学问题的数
值计算方法及其软件实现.计算数学几乎与数学科学的一切
分支有联系，它利用数学领域的成果发展了新的更有效的
算法及其理论，反过来很多数学分支都需要探讨和研究适
用于计算机的数值方法.因此，"计算方法"内容十分广泛.但
实际问题 程序设计
数学问题 上机计算
提供计算方法 结果分析
12
基本的数学问题:
1.大型线性代数方程组Ax=b求解;
2.矩阵A的特征值和特征向量计算;
3.非线性方程 f ( x ) 0 求解(求根);
4.积分 b a
f
( x)dx计算;
5.常微分方程初值问题求解;
6.其它。
13
求精确解(值)一般非常困难。例如：
17
截断误差 在求解过程中，往往以近似替代， 化繁为简，这样产生的误差称为截断误差。
舍入误差 在计算机上运算时受机器字长的 限制，一般必须进行舍入，此时产生的误 差称为舍入误差。
18
3. 截断误差，如
sin xxx3 x5 ....,.. 3! 5!
7
数值计算方法或数值分析主要是研究如何 运用计算机去获得数学问题的数值解的理 论和方法.对那些在经典数学中,用解析方法 在理论上已作出解的存在,但要求出他的解 析解又十分困难,甚至是不可能的这类数学 问题,数值解法就显得不可缺少,同时有十分 有效.
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
高等数学、线性代数、程序设计语言
计算方法——09计11、61
5/47
教材
《计算方法》,
易大义等,
浙江大学出版社, 2002年第2版
计算方法——09计11、61
6/47
主要参考书
1.《数值分析引论》,
易大义 陈道琦,
浙江大学出版社, 1998年第1版
计算方法——09计11、61
7/47
主要参考书
简化，抽象问题后建立的数学模型与实际问
2.
观测误差 /* Measurement Error */ 流、温度等）
观测和实验得到的参量（物理量为电压、电
计算方法——09计11、61
21/47
§1.2 误差的种类及其来源
3.
截断误差（方法误差）/*Truncation Error*/ 求和，只能取前面有限项求和来近似代 替）。这种计算方法本身出现的误差，所 以也称为方法误差。如 3 5 x x si n x x ......, 3! 5!
有限过程代替无限过程的误差（无穷级数
右端是截断误差。
x x si n x ( x ) ...... 3! 5!
计算方法——09计11、61
22/47
3
5
§1.2 误差的种类及其来源
4.
舍入误差 /* Roundoff Error */
计算机字长有限，一般实数不能精确存储，于 是产生舍入误差。 例如：在 10 位十进制数限制下：
e( x ) 称 er ( x ) 为 x * 的相对误差。 x * e( x ) 实用中，常用 表示 x 的相对误差。 x 称 er ( x ) 为 x * 的相对误差限。
计算方法——09计11、61
26/47
§1.4 有效数字及其与误差的关系
一．有效数字 /* significant digits */
1、1 3 0.3333333333 (本应1 3 0.3333333333 3 ） 2、(1.000002 ) 2 1.000004 0
2 (本应（ 1.000002 ） 1.000004
1.0000040000 04 1.000004 0.0000000000 04 4 1012）
计算方法——09计11、61
2/47
联系方式
课程网站
http://202.195.66.5/start/StudyNA/
(校内访问)
E－mail地址
xznuckj@
办公室
泉山9＃701室
计算方法——09计11、61
3/47
本课程成绩的组成
平时成绩（占15％）：包括出勤、课堂提 问、讨论情况等。 实验成绩（占25％）：包括出勤、实验报 告（预习报告）等。 期末成绩（占60％）。
一定要从规格化后的数来判断其位数
有效位数与第一个非 0 项后的数字个数是不
一致的。 四舍五入所得到的数是一致的。
计算方法——09计11、61
27/47
§1.4 有效数字及其与误差的关系
例3：设 x * ＝ 0.0270是某数 x经“四舍五入”所得，
则误差 e ( x *) 不超过 x *末位的的半个单位，即： 1 x x * 10 4 2 又 x * 101 (0.270) 故该不等式又可写为 1 x x * 101 3 2 由有效数字定义可知, x * 有3位有效数字, 分别是2， 7， 0。
2.《数值分析基础教程》,
李庆杨,
高等教育出版社,
2001年第1版
Байду номын сангаас
计算方法——09计11、61
8/47
主要参考书
3.《数值方法和MATLAB实现与应用》, （美） Gerald Recktenwald 著 伍卫国 万群 张辉 等译, 机械工业出版社,
2004年第1版
其他各类有关 “数值分析” 和 “计算方法” 的书
排长对班长:
明晚 8 点, 营长将带着哈雷彗星在礼 堂中出现, 这是每隔 76 年才有的事。如果下雨的 话,营长将命令彗星穿上野战服到操场上去。 在明晚 8 点下雨的时候, 著名的 76 岁 哈雷将军将在营长的陪同下身着野战服, 开着他那 “彗星”牌汽车, 经过操场前往礼堂。
计算方法——09计11、61
数值计算方法
数理学院数学系
诚信声明
本课件中大量采用网络及其他渠道搜集的相关文 字信息和图片信息，这些信息因数量较大，无法 一一列明出处，本人在此郑重声明：这些相关资 料的版权归原作者所有，本文引用仅仅用于教学 目的。如有不妥，请与本人联系，联系方式： xznuckj@。 在此对相关资料的作者所付出的辛勤劳动表示衷 心的感谢，并对作者表示诚挚敬意！ 本人郑重承诺：尊重知识，尊重劳动，尊重版权， 学术诚信。
舍入误差很小，本课程将研究它在运算过程中 是否能有效控制。
计算方法——09计11、61
23/47
§1.2 误差的种类及其来源
据说, 美军1910年的一次部队的命令传递是这样的:
营长对值班军官: 明晚大约 8点钟左右, 哈雷彗星将
可能在这个地区看到, 这种彗星每隔 76年才能看见 一次。命令所有士兵着野战服在操场上集合, 我将 向他们解释这一罕见的现象。如果下雨的话, 就在 礼堂集合, 我为他们放一部有关彗星的影片。
设 x (0.a1a 2 a n a p )10m (a1 0, p ) 1 若 x * x 10m n ( P10 1.4.2) 2 则说 x *具有n位有效数字，分别是 a1 , a 2 , , a n 若n＝p，则称 x *为有效数。
插值法 数值逼近 本 课 程 的 内 容 数据拟合的最小二乘法 数值积分和数值微分* 线性方程组的求解 数值代数 非线性方程组的求解 矩阵特征值* 常微分方程的数值方法
计算方法——09计11、61
11/47
第一章
数值计算中的误差
3 学时
本章内容
§1.1 引言 §1.2 误差的种类及其来源 §1.3 绝对误差和相对误差 §1.4 有效数字及其与误差的关系 §1.5 误差的传播与估计 §1.6 选用算法应遵循的原则 小结 作业与实验
25/47
班长对士兵:
§1.3 绝对误差和相对误差
一．绝对误差 /* absolute error */
设
称
x
——准确值，x * ——近似值。
*
*
e( x ) x x 为 x 的绝对误差(简称误差) | e( x ) | 为 x * 的绝对误差限。
二．相对误差 /* relative error */
计算方法——09计11、61
13/47
本章要求
1. 熟悉计算方法在解决实际问题中所处的地位, 熟悉计算方法是以计算机为工具求近似解的数 值方法； 2. 熟悉绝对误差（限），相对误差（限）及有 效数字概念； 3. 熟悉公式；
4. 熟悉选用算法应遵循的原则。
计算方法——09计11、61
14/47
值班军官对连长: 根据营长的命令, 明晚8点哈雷彗
星将在操场上空出现。如果下雨的话, 就让士兵穿 着野战服列队前往礼堂, 这一罕见的现象将在那里 出现。
计算方法——09计11、61
24/47
§1.2 误差的种类及其来源
连长对排长:
根据营长的命令, 明晚 8 点, 非凡的哈 雷彗星将身穿野战服在礼堂中出现。如果操场上 下雨, 营长将下达另一个命令, 这种命令每隔 76 年 才会出现一次。
例1，例2 的结果的根源 建 立 算 法 截方 断法 误误 差差 ) ( 上 机 计 算 结
果
模 型 误 差
观 测 误 差
初 值 误 差
舍 入 误 差
20/47
)
计算方法——09计11、61
§1.2 误差的种类及其来源
二. 误差分类
1.
模型误差（描述误差）/* Modeling Error */ 题之差。
（P11 1.4.4 ）
1 反之，若 er ( x ) 10 ( n1) 2(a1 1) 则x 至少有n位有效数字。
计算方法——09计11、61
30/47
§1.4 有效数字及其与误差的关系
证： 因a1 10m x * (a1 1) 10m，故当x *
有n位有效数字时， 0.5 10m n1 1 n1 r 10 a1 10m 2a1 x 反之，由 1 x r (a1 1) 10 10 n1 2(a1 1)
§1.1 引言
解决科学技术和工程问题的步骤:
实际问题 数学问题 提供计算方法
程序设计
上机计算
结果分析
什么是数值计算方法: 将所预求解的数学模 型简化成一系列算术运算和逻辑运算, 以便 在计算机上求解, 并对算法的稳定性、收敛 性和误差进行分析。
计算方法——09计11、61
15/47
§1.1 引言
计算方法——09计11、61
29/47
§1.4 有效数字及其与误差的关系
二．有效数位与误差的关系 e ( x *)
1. 有效数位n越多，则绝对误差e ( x ) 越小 （由定义1.4.2 ） 2. 定理：若近似数x 具有n位有效数字，则 1 ( n 1 ) er ( x ) 10 2a1
计算方法——09计11、61
17/47
计算方法——09计11、61
5/47
教材
《计算方法》,
易大义等,
浙江大学出版社, 2002年第2版
计算方法——09计11、61
6/47
主要参考书
1.《数值分析引论》,
易大义 陈道琦,
浙江大学出版社, 1998年第1版
计算方法——09计11、61
7/47
主要参考书
简化，抽象问题后建立的数学模型与实际问
2.
观测误差 /* Measurement Error */ 流、温度等）
观测和实验得到的参量（物理量为电压、电
计算方法——09计11、61
21/47
§1.2 误差的种类及其来源
3.
截断误差（方法误差）/*Truncation Error*/ 求和，只能取前面有限项求和来近似代 替）。这种计算方法本身出现的误差，所 以也称为方法误差。如 3 5 x x si n x x ......, 3! 5!
有限过程代替无限过程的误差（无穷级数
右端是截断误差。
x x si n x ( x ) ...... 3! 5!
计算方法——09计11、61
22/47
3
5
§1.2 误差的种类及其来源
4.
舍入误差 /* Roundoff Error */
计算机字长有限，一般实数不能精确存储，于 是产生舍入误差。 例如：在 10 位十进制数限制下：
e( x ) 称 er ( x ) 为 x * 的相对误差。 x * e( x ) 实用中，常用 表示 x 的相对误差。 x 称 er ( x ) 为 x * 的相对误差限。
计算方法——09计11、61
26/47
§1.4 有效数字及其与误差的关系
一．有效数字 /* significant digits */
1、1 3 0.3333333333 (本应1 3 0.3333333333 3 ） 2、(1.000002 ) 2 1.000004 0
2 (本应（ 1.000002 ） 1.000004
1.0000040000 04 1.000004 0.0000000000 04 4 1012）
计算方法——09计11、61
2/47
联系方式
课程网站
http://202.195.66.5/start/StudyNA/
(校内访问)
E－mail地址
xznuckj@
办公室
泉山9＃701室
计算方法——09计11、61
3/47
本课程成绩的组成
平时成绩（占15％）：包括出勤、课堂提 问、讨论情况等。 实验成绩（占25％）：包括出勤、实验报 告（预习报告）等。 期末成绩（占60％）。
一定要从规格化后的数来判断其位数
有效位数与第一个非 0 项后的数字个数是不
一致的。 四舍五入所得到的数是一致的。
计算方法——09计11、61
27/47
§1.4 有效数字及其与误差的关系
例3：设 x * ＝ 0.0270是某数 x经“四舍五入”所得，
则误差 e ( x *) 不超过 x *末位的的半个单位，即： 1 x x * 10 4 2 又 x * 101 (0.270) 故该不等式又可写为 1 x x * 101 3 2 由有效数字定义可知, x * 有3位有效数字, 分别是2， 7， 0。
2.《数值分析基础教程》,
李庆杨,
高等教育出版社,
2001年第1版
Байду номын сангаас
计算方法——09计11、61
8/47
主要参考书
3.《数值方法和MATLAB实现与应用》, （美） Gerald Recktenwald 著 伍卫国 万群 张辉 等译, 机械工业出版社,
2004年第1版
其他各类有关 “数值分析” 和 “计算方法” 的书
排长对班长:
明晚 8 点, 营长将带着哈雷彗星在礼 堂中出现, 这是每隔 76 年才有的事。如果下雨的 话,营长将命令彗星穿上野战服到操场上去。 在明晚 8 点下雨的时候, 著名的 76 岁 哈雷将军将在营长的陪同下身着野战服, 开着他那 “彗星”牌汽车, 经过操场前往礼堂。
计算方法——09计11、61
数值计算方法
数理学院数学系
诚信声明
本课件中大量采用网络及其他渠道搜集的相关文 字信息和图片信息，这些信息因数量较大，无法 一一列明出处，本人在此郑重声明：这些相关资 料的版权归原作者所有，本文引用仅仅用于教学 目的。如有不妥，请与本人联系，联系方式： xznuckj@。 在此对相关资料的作者所付出的辛勤劳动表示衷 心的感谢，并对作者表示诚挚敬意！ 本人郑重承诺：尊重知识，尊重劳动，尊重版权， 学术诚信。
舍入误差很小，本课程将研究它在运算过程中 是否能有效控制。
计算方法——09计11、61
23/47
§1.2 误差的种类及其来源
据说, 美军1910年的一次部队的命令传递是这样的:
营长对值班军官: 明晚大约 8点钟左右, 哈雷彗星将
可能在这个地区看到, 这种彗星每隔 76年才能看见 一次。命令所有士兵着野战服在操场上集合, 我将 向他们解释这一罕见的现象。如果下雨的话, 就在 礼堂集合, 我为他们放一部有关彗星的影片。
设 x (0.a1a 2 a n a p )10m (a1 0, p ) 1 若 x * x 10m n ( P10 1.4.2) 2 则说 x *具有n位有效数字，分别是 a1 , a 2 , , a n 若n＝p，则称 x *为有效数。
插值法 数值逼近 本 课 程 的 内 容 数据拟合的最小二乘法 数值积分和数值微分* 线性方程组的求解 数值代数 非线性方程组的求解 矩阵特征值* 常微分方程的数值方法
计算方法——09计11、61
11/47
第一章
数值计算中的误差
3 学时
本章内容
§1.1 引言 §1.2 误差的种类及其来源 §1.3 绝对误差和相对误差 §1.4 有效数字及其与误差的关系 §1.5 误差的传播与估计 §1.6 选用算法应遵循的原则 小结 作业与实验
25/47
班长对士兵:
§1.3 绝对误差和相对误差
一．绝对误差 /* absolute error */
设
称
x
——准确值，x * ——近似值。
*
*
e( x ) x x 为 x 的绝对误差(简称误差) | e( x ) | 为 x * 的绝对误差限。
二．相对误差 /* relative error */
计算方法——09计11、61
13/47
本章要求
1. 熟悉计算方法在解决实际问题中所处的地位, 熟悉计算方法是以计算机为工具求近似解的数 值方法； 2. 熟悉绝对误差（限），相对误差（限）及有 效数字概念； 3. 熟悉公式；
4. 熟悉选用算法应遵循的原则。
计算方法——09计11、61
14/47
值班军官对连长: 根据营长的命令, 明晚8点哈雷彗
星将在操场上空出现。如果下雨的话, 就让士兵穿 着野战服列队前往礼堂, 这一罕见的现象将在那里 出现。
计算方法——09计11、61
24/47
§1.2 误差的种类及其来源
连长对排长:
根据营长的命令, 明晚 8 点, 非凡的哈 雷彗星将身穿野战服在礼堂中出现。如果操场上 下雨, 营长将下达另一个命令, 这种命令每隔 76 年 才会出现一次。
例1，例2 的结果的根源 建 立 算 法 截方 断法 误误 差差 ) ( 上 机 计 算 结
果
模 型 误 差
观 测 误 差
初 值 误 差
舍 入 误 差
20/47
)
计算方法——09计11、61
§1.2 误差的种类及其来源
二. 误差分类
1.
模型误差（描述误差）/* Modeling Error */ 题之差。
（P11 1.4.4 ）
1 反之，若 er ( x ) 10 ( n1) 2(a1 1) 则x 至少有n位有效数字。
计算方法——09计11、61
30/47
§1.4 有效数字及其与误差的关系
证： 因a1 10m x * (a1 1) 10m，故当x *
有n位有效数字时， 0.5 10m n1 1 n1 r 10 a1 10m 2a1 x 反之，由 1 x r (a1 1) 10 10 n1 2(a1 1)
§1.1 引言
解决科学技术和工程问题的步骤:
实际问题 数学问题 提供计算方法
程序设计
上机计算
结果分析
什么是数值计算方法: 将所预求解的数学模 型简化成一系列算术运算和逻辑运算, 以便 在计算机上求解, 并对算法的稳定性、收敛 性和误差进行分析。
计算方法——09计11、61
15/47
§1.1 引言
计算方法——09计11、61
29/47
§1.4 有效数字及其与误差的关系
二．有效数位与误差的关系 e ( x *)
1. 有效数位n越多，则绝对误差e ( x ) 越小 （由定义1.4.2 ） 2. 定理：若近似数x 具有n位有效数字，则 1 ( n 1 ) er ( x ) 10 2a1
计算方法——09计11、61
17/47