试题清华第五版数值分析第1章ppt课件
数值分析第一章ppt
s 某商品标注重量为 27±0.5kg, 实际重量是多少？
}
1.2.2 相对误差和相对误差限
x*的相对误差
r
x x x
在不同近似值中，|εr (x)|越小，x*的精确度越高。
r(x)
| ( x) | |x|
x x x
——x*的相对误差限
常用计算公式： r ( x)
( x)
x*Βιβλιοθήκη x x* x* ,}
（2）相对误差：
r ( x1
x2 )
( x1 x2 )
x1 x2
( x1 ) ( x2 )
x1 x2
x1 x1 x2
( x1 )
x1
x2 x1 x2
(x2 )
x2
x1
x1
x2
r
(
x1
)
x1
x2
x2
r
(
x2
)
当x1≈x2时, x1 – x2 ≈0,所以相近两数之差的相对误差将很大 。
}
1.2 误差的基本估计方法
= 1.2.1 绝对误差和绝对误差限 = 1.2.2 相对误差和相对误差限 = 1.2.3 有效数字 = 1.2.4 算术运算的误差
}
1.2.1 绝对误差和绝对误差限
设某准确值x近似值为x*。 x*的绝对误差 ε(x)=x–x*
在同一量的不同近似值中，|ε(x)|越小，x*的精确度越高。
sin x x x3 x5 x7 , x 3! 5! 7!
用近似计算公式 sin x x
截断误差 sin x x x3 x5 x7 cos x3
3! 5! 7!
3!
sin x x 1 x 3 6
数值分析第一章PPT
1.1.2 计算数学与科学计算 现代科学的三个组成部分: 科学理论, 科学实验, 科学计算 科学计算 的核心内容是以现代化的计算机及数学软件 （Matlab, Mathematica, Maple, MathCAD etc. ）为工具，以数学 模型为基础进行模拟研究。
一些边缘学科的相继出现：
计算数学,计算物理学,计算力学,计算化学,计算生物学, 计算地质学,计算经济学,等等
取 0 e
1
x2
dx S4 ,
S4
R4
/* Remainder */
1 1 1 1 由留下部分 称为截断误差 /* Truncation Error */ 4! 9 5! 11 /* included terms */ 1 1 这里 R4 引起.005 0 由截去部分 4! 9 /* excluded terms */ 1 1 1 S4 1 1 0 .333 0 .1 0 .024 0 .743 引起 3 10 42 | 舍入误差 /* Roundoff Error */ | 0.0005 2 0.001
数值分析
第1章
数值分析与科学计算引论
§1.1 数值分析的对象、作用与特点
1.1.1 什么是数值分析 数值分析是计算数学的主要部分,计算数学是数学 科学的一个分支,它研究用计算机求解各种数学问题的 数值计算方法及其理论与软件实现.这门课程又称为(数 值)计算方法、科学与工程计算等。
•
在电子计算机成为数值计算的主要工具的今天， 需要研究适合计算机使用的数值计算方法。使用计 算机解决科学计算问题时大致要经历如下几个过程:
造成这种情况的是不稳定的算法 /* unstable algorithm */ 我们有责任改变。
高等数学第五版第一章ppt1(4)
1 此时对 M = , ∃δ > 0, 使得当 ∴ ∀ε > 0,
ε
1 ∴ 当x → x 0时, 为无穷小. f ( x)
0 < x − x0 < δ时,有 f ( x ) > M = ,即 1 < ε. ε f ( x)
1
17
无穷小与无穷大
定理4 在同一过程中,无穷大的倒数为无穷小; 定理4 在同一过程中,无穷大的倒数为无穷小; 恒不为零的无穷小的倒数为无穷大. 恒不为零的无穷小的倒数为无穷大.
积仍为无穷大; 积仍为无穷大 用无零值有界变量去除无穷大仍为无穷大. ∗ 用无零值有界变量去除无穷大仍为无穷大 例 求 lim ( x + 1 +
x → +∞
x)
解 lim ( x + 1 +
x → +∞
x ) = +∞
19
无穷小与无穷大
四、小结
无穷小的概念; 无穷小的概念 无穷小与函数极限的关系; 无穷小与函数极限的关系 无穷小的运算; 无穷小的运算; 无穷大的概念; 无穷大的概念 无穷小与无穷大的关系. 无穷小与无穷大的关系
⇐
于是
设 f ( x ) = A + α( x ),
| f ( x ) − A |=| α( x ) |
其中A 其中 是常数 , α( x )是当x → x0时的无穷小,
∀ε > 0, ∃δ > 0,当0 <| x − x0 |< δ , 恒有
| α( x ) |< ε
lim 即 | f ( x ) − A |< ε. ∴ x → x f ( x ) = A.
<M⋅
数值分析(清华大学出版社)第1,2章
2.
x 的相对误差是
x x d x er ( x ) d ln x x x
它是对数函数的微分。
设 u = xy , 则 lnu=lnx+lny , 因而 dlnu = dlnx + dlny
e r ( u ) e r ( x ) e r ( y ) r ( u ) r ( x ) r ( y )
即 m- n = - 2, m=1, n = 3, 所以 x = 3.14 作为 近似值 时, 就有3 位有效数字。
16
四、 相对误差限与有效数字的关系
定理1
设近似值
x 0.a1a2 an 10m
有n 位有效数字， a1 0 。则其相对误差限为 1 n 1 r (x ) 10 2a1 x 0.a1a2 an 10m 故 证明
a1 10
m 1
| x | (a1 1) 10
m 1
r (x )
x x x
0 .5 10m n 1 10 n1 a1 10m 1 2a1
17
定理2 设近似值 x 0.a1a2 an 10 的相对误差限
m
1 10 n1 ,则它至少有n 位有效数字。 不大于 2( a 1) 1
25
对多元函数 y f ( x1 , x2 , , xn ), 自变量的近似值为 x1 , x2 ,, xn , y 的近似值为 y f ( x1 , x2 , , xn ),
y 的运算误差为 函数值
e ( y ) e[ f ( x1 , x2 , , xn )] df ( x1 , x2 , , xn )
课件-数值分析(第五版)1-3章
x x
f ( x) f ( x* ) f ( x)
x x
xf ( x) f ( x)
C p 10 即认为是病态
f ( x) x n
9 第1章 数值分析与科学计算引论
研究对象 作用特点
数值计算 误差
误差分析 避免危害
数值计算 算法设计
数学软件
2. 算法的数值稳定性 定义3 一个算法如果输入数据有误差，而在计算过程中舍入误 差不增长，则称此算法是数值稳定的，否则称此算法为不稳定 的。 例1.1：P.９ I n e
x 0.003
y 1
2017/3/12
1000
1.00314 , y * 1.003
6 第1章 数值分析与科学计算引论
研究对象 作用特点
数值计算 误差
误差分析 避免危害
数值计算 算法设计
数学软件
注: 有效位数与小数点后有多少位无关； m相同情况下，有效位数越多，误差限越小； 相对误差及相对误差限是无量纲的，绝对误差及误差限是有量纲的。
数值计算 算法设计
数学软件
1.1 数值分析的对象、作用与特点
1 研究对象
用计算机求解数学问题的数值计算方法、理论及软件实现
实际问题 数学模型 数值计算方法 程序设计（数学软件） 上机计算求出结果
应用数学
计算数学即数值分析
数值分析（计算方法） 插值与函数逼近（2、3）数值微分与数值积分（4） 的研究对象
第一章习题
1, 5,７,12,1４
谢
谢 !
2017/3/12
14 第1章 数值分析与科学计算引论
第2章 插值法
引言
拉格朗日（Lagrange）插值 均差与牛顿（Newton）插值 埃尔米特（Hermite）插值 分段低次插值 三次样条插值
数值分析第五版
数值分析第五版第一章绪论1．设0x >,x 的相对误差为δ，求ln x 的误差。
解：近似值*x 的相对误差为*****r e x x e x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈ 进而有(ln *)x εδ≈2．设x 的相对误差为2%，求nx 的相对误差。
解：设()n f x x =，则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=, 1||n p x nx C n n-?∴== 又((*))(*)r p r x n C x εε≈?且(*)r e x 为2((*))0.02n r x n ε∴≈5计算球体积要使相对误差限为1，问度量半径R 时允许的相对误差限是多少？解：球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===(*)(*)3(*)r p r r V C R R εεε∴≈=又(*)1r V ε=故度量半径R 时允许的相对误差限为1(*)10.333r R ε=≈ 9．正方形的边长大约为了100cm ，应怎样测量才能使其面积误差不超过21cm ？解：正方形的面积函数为2()A x x =(*)2*(*)A A x εε∴=.当*100x =时，若(*)1A ε≤,则21(*)102x ε-≤故测量中边长误差限不超过0.005cm 时，才能使其面积误差不超过21cm第二章插值法1．当1,1,2x =-时，()0,3,4f x =-,求()f x 的二次插值多项式。
解：0120121200102021101201220211,1,2,()0,()3,()4;()()1()(1)(2)()()2()()1()(1)(2)()()6()()1()(1)(1)()()3x x x f x f x f x x x x x l x x x x x x x x x x x l x x x x x x x x x x xl x x x x x x x ==-===-=--==-+-----==------==-+--则二次拉格朗日插值多项式为220()()k k k L x y l x ==∑0223()4()14(1)(2)(1)(1)23537623l x l x x x x x x x =-+=---+-+=+- 2．给出()ln f x x =的数值表 X 0.4 0.5 0.6 0.7 0.8 lnx-0.916291-0.693147-0.510826-0.356675-0.223144用线性插值及二次插值计算ln 0.54的近似值。
数值分析课件第1章
解：
(s ) l (d ) d (l )
110 (0.1) 80 (0.2) 27( m 2 )
r
(
s
)
(s)
s
(s)
ld
27 0.31% 8800
2、函数误差 当自变量有误差时计算函数值也产生误差，可以利用
函数的泰勒展开式进行估计。
工科研究生公共课程数学系列
（f (x)） f (x)（x）.
例1-4 设x0,x的相对误差为，求lnx的误差。
解：
lnx*
-lnx
1 x*
(x*
-
x), 即有
e(lnx) er(x)
进而有(ln(x)) 。
工科研究生公共课程数学系列
机动 上页 下页 首页 结束
1.3 误差定性分析与避免误差危害
一、几种定性分析误差的方法 1、概率分析法：考虑到误差分布的随机性，用概率统计的
二、数值分析的特点
• 面向计算机，要根据计算机的特点提供切实可行的有效算 法。
• 有可靠的理论分析，能任意逼近并达到精度要求，对近似 算法要保证收敛性和数值稳定性，还要对误差进行分析。 这些都是建立在数学理论的基础上，因此不应片面的将数 值分析理解为各种数值方法的简单罗列和堆积。
• 要有好的计算复杂性，时间复杂性好是指节省时间，空间 复杂性好是指节省存上实现。
• 要有数值实验，即任何一个算法除了从理论上要满足上述 三点外，还要通过数值实验证明是行之有效的。
工科研究生公共课程数学系列
机动 上页 下页 首页 结束
三、数值分析的学习方法 初学可能仍会觉得公式多，理论分析复杂。给出如下的 几点学习方法。
• 认识建立算法和对每个算法进行理论分析是基本任务，主 动适应公式多和讲究理论分析的特点。
数值分析原理课件第一章
第一章 绪 论本章以误差为主线，介绍了计算方法课程的特点，并概略描述了与算法相关的基本概念，如收敛性、稳定性，其次给出了误差的度量方法以及误差的传播规律，最后，结合数值实验指出了算法设计时应注意的问题．§1.1 引 言计算方法以科学与工程等领域所建立的数学模型为求解对象，目的是在有限的时间段内利用有限的计算工具计算出模型的有效解答。
由于科学与工程问题的多样性和复杂性，所建立的数学模型也是各种各样的、复杂的. 复杂性表现在如下几个方面：求解系统的规模很大，多种因素之间的非线性耦合，海量的数据处理等等，这样就使得在其它课程中学到的分析求解方法因计算量庞大而不能得到计算结果，且更多的复杂数学模型没有分析求解方法． 这门课程则是针对从各种各样的数学模型中抽象出或转化出的典型问题，介绍有效的串行求解算法，它们包括(1)非线性方程的近似求解方法；(2)线性代数方程组的求解方法；(3)函数的插值近似和数据的拟合近似；(4)积分和微分的近似计算方法；(5)常微分方程初值问题的数值解法；(6)优化问题的近似解法；等等从如上内容可以看出，计算方法的显著特点之一是“近似”． 之所以要进行近似计算，这与我们使用的工具、追求的目标、以及参与计算的数据来源等因素有关．计算机只能处理有限数据，只能区分、存储有限信息，而实数包含有无穷多个数据，这样，当把原始数据、中间数据、以及最终计算结果用机器数表示时就不可避免的引入了误差，称之为舍入误差．我们需要在有限的时间段内得到运算结果，就需要将无穷的计算过程截断，从而产生截断误差． 如的计算是无穷过程，当用作为的 +++=!21!111e !1!21!111n e n ++++= e 近似时，则需要进行有限过程的计算，但产生了截断误差．e e n - 当用计算机计算时，因为舍入误差的存在，我们也只能得到的近似值，也就是n e n e *e 说最终用近似，该近似值既包含有舍入误差，也包含有截断误差．*e e 当参与计算的原始数据是从仪器中观测得来时，也不可避免得有观测误差．由于这些误差的大量存在，我们得到的只能是近似结果，进而对这些结果的“可靠性”进行分析就是必须的，它成为计算方法的第二个显著特点． 可靠性分析包括原问题的适定性和算法的收敛性、稳定性．所谓适定性问题是指解存在、惟一，且解对原始数据具有连续依赖性的问题． 对于非适定问题的求解，通常需要作特殊的预处理，然后才能做数值计算． 在这里，如无特殊说明，都是对适定的问题进行求解．对于给定的算法，若有限步内得不到精确解，则需研究其收敛性． 收敛性是研究当允许计算时间越来越长时，是否能够得到越来越可靠的结果，也就是研究截断误差是否能够趋于零．对于给定的算法，稳定性分析是指随着计算过程的逐步向前推进，研究观测误差、舍入对于同一类模型问题的求解算法可能不止一种，常希望从中选出高效可靠的求解算法． 如我国南宋时期著名的数学家秦九韶就提出求n 次多项式值0111a x a x a x a n n nn ++++-- 的如下快速算法；n a s =； k n a t -=t sx s +=),,2,1(n k =它通过n 次乘法和n 次加法就计算出了任意n 次多项式的值． 再如幂函数可以通过如下64x 快速算法计算出其值；x s =；循环6次s s s ⋅=如上算法仅用了6次乘法运算，就得到运算结果．算法最终需要在计算机上运行相应程序，才能得到结果，这样就要关注算法的时间复杂度(计算机运行程序所需时间的度量)、空间复杂度(程序、数据对存储空间需求的度量)和逻辑复杂度(关联程序的开发周期、可维护性以及可扩展性)． 事实上，每一种算法都有自己的局限性和优点，仅仅理论分析是很不够的，大量的实际计算也非常重要，结合理论分析以及相当的数值算例结果才有可能选择出适合自己关心问题的有效求解算法． 也正因如此，只有理论分析结合实际计算才能真正把握准算法．§1.2 误差的度量与传播一、误差的度量误差的度量方式有绝对误差、相对误差和有效数字．定义1.1 用作为量的近似，则称为近似值的绝对误差．*x x )(:**x e x x =-*x 由于量x 的真值通常未知，所以绝对误差不能依据定义求得，但根据测量工具或计算情况，可以估计出绝对误差绝对值的一个较小上界，即有ε (1.1)ε≤-=x x x e **)(称正数为近似值的绝对误差限，简称误差． 这样得到不等式ε*x εε+≤≤-**x x x 工程中常用ε±=*x x 表示近似值的精度或真值x 所在的范围．*x 误差是有量纲的，所以仅误差数值的大小不足以刻划近似的准确程度． 如量 (1.2)m m cm s μ50001230000005.023.15.0123±=±=±=为此，我们需要引入相对误差定义1.2 用作为量的近似，称为近似值的相对误差． 当0*≠x x )(:**x e xxx r =-*x 是x 的较好近似时，也可以用如下公式计算相对误差*x (1.3)***)(xx x x e r -= 显然，相对误差是一个无量纲量，它不随使用单位变化． 如式(1.2)中的量s 的近似，无论使用何种单位，它的相对误差都是同一个值．同样地，因为量x 的真值未知，我们需要引入近似值的相对误差限，它是相*x )(*x r ε对误差绝对值的较小上界． 结合式(1.1)和(1.3)，相对误差限可通过绝对误差限除以近似*x 值的绝对值得到，即(1.4)***)()(xx x r εε=为给出近似数的一种表示法，使之既能表示其大小，又能体现其精确程度，需引入有效数字以及有效数的概念．定义1.3 设量的近似值有如下标准形式x *x p n ma a a a x 21*.010⨯±= (1.5)()pm p n m n m m a a a a ----⨯++⨯++⨯+⨯±101010102211 ＝其中且，m 为近似值的量级． 如果使不等式}9,,1,0{}{1 ⊂=pi i a 01≠a (1.6)n m x x -⨯≤-1021*成立的最大整数为n ，则称近似值具有n 位有效数字，它们分别是、、… 和 ． *x 1a 2a n a 特别地，如果有，即最后一位数字也是有效数字，则称是有效数．p n =*x 从定义可以看出，近似数是有效数的充分必要条件是末位数字所在位置的单位一半是绝对误差限． 利用该定义也可以证明，对真值进行“四舍五入”得到的是有效数． 对于有效数，有效数字的位数等于从第一位非零数字开始算起，该近似数具有的位数． 注意，不能给有效数的末位之后随意添加零，否则就改变了它的精度．例1.1 设量，其近似值，，． 试回答这三个近π=x 141.3*1=x 142.3*2=x 722*3=x 似值分别有几位有效数字，它们是有效数吗？解 这三个近似值的量级，因为有1=m 312*110211021005.000059.0--⨯=⨯=≤=- x x 413*2102110210005.00004.0--⨯=⨯=≤=- x x571428571428.3*3=x 312*310211021005.0001.0--⨯=⨯=≤=- x x 所以和都有3位有效数字，但不是有效数． 具有4位有效数字，是有效数．*1x *3x *2x 二、误差的传播这里仅介绍初值误差传播，即假设自变量带有误差，函数值的计算不引入新的误差． 对于函数有近似值，利用在点处),,,(21n x x x f y =),,,(**2*1*n x x x f y =),,,(**2*1n x x x 的泰勒公式(Taylor Formula)，可以得到 )(),,,()(*1**2*1**i i ni n i x x x x xf y y y e -≈-=∑=(1.7))(),,,(*1**2*1i ni n i x e x x xf ∑== 其中，是的近似值，是的绝对误差． 式(1.7)表明函ii x f f ∂∂=:*i x i x )(*i x e *i x ),,2,1(n i =数值的绝对误差近似等于自变量绝对误差的线性组合，组合系数为相应的偏导数值． 从式(1.7)也可以推得如下函数值的相对误差传播近似计算公式 (1.8))(),,,()(***1**2*1*i r i ni ni r x e yx x x x f y e ∑=≈对于一元函数，从式(1.7)和(1.8)可得到如下初值误差传播近似计算公式)(x f y = (1.9))()()(***x e x f y e '≈ (1.10))()()(*****x e yx x f y e r r '≈式(1.9)表明，当导数值的绝对值很大时，即使自变量的绝对误差比较小，函数值的绝对误差也可能很大．例1.2 试建立函数的绝对误差(限)、相对误差n n x x x x x x f y +++== 2121),,,(的近似传播公式，以及时的相对误差限传播公式．{}ni i x 1*0=> 解 由公式(1.7)和(1.8)可分别推得和的绝对误差、相对误差传播公式如下(1.11)∑∑==≈ni i ini nix e x e x x xf y e 1**1**2*1*)()(),,,()(＝ (1.12)∑∑==≈ni i r i i r i ni ni r x e yx x e y x x x x f y e 1******1**2*1*)()(),,,()(＝ 进而有∑∑∑===≤≤≈ni i n i ini ix x e xe y e 1*1*1**)()()()(ε于是有和的绝对误差限近似传播公式 ∑=≈ni ixy 1**)()(εε当时，由式(1.3)推得相对误差限的近似传播公式{}ni i x 1*=>)(max )(max )(max )()()(*11***11***11****1**i r ni ni i ir n i ni i i r n i ni i r i ni ir x yx x yx x x y x yxy εεεεεε≤≤=≤≤=≤≤====≤=≈∑∑∑∑ 例1.3使用足够长且最小刻度为1mm 的尺子，量得某桌面长的近似值3.1304*=a mm ，宽的近似值mm (数据的最后一位均为估计值)． 试求桌子面积近似值的绝8.704*=b 对误差限和相对误差限．解 长和宽的近似值的最后一位都是估计位，尺子的最小刻度是毫米，故有误差限 mm ，mm 5.0)(*=a ε5.0)(*=b ε面积，由式(1.7)得到近似值的绝对误差近似为ab S =***b a S = )()()(*****b e a a e b S e +≈进而有绝对误差限 mm 255.10045.03.13045.08.704)()()(*****=⨯+⨯=+≈b a a b S εεε相对误差限 %11.00011.08.7043.130455.1004)()(***=≈⨯=≈S S S r εε§1.3 数值实验与算法性能比较本节通过几个简单算例说明解决同一个问题可以有不同的算法，但算法的性能并不完全相同，他们各自有自己的适用范围，并进而指出算法设计时应该注意的事项． 算例1.1 表达式，在计算过程中保留7位有效数字，研究对不同)1(1111+=+-x x x x 的x ，两种计算公式的计算精度的差异．说明1：Matlab 软件采用IEEE 规定的双精度浮点系统，即64位浮点系统，其中尾数占52位，阶码占10位，尾数以及阶码的符号各占1位． 机器数的相对误差限(机器精度)eps=2－52≈2.220446×10－16，能够表示的数的绝对值在区间(2.2250739×10－308，1.797693×10308)内，该区间内的数能够近似表达，但有舍入误差，能够保留至少15位有效数字． 其原理可参阅参考文献[2, 4]．分析算法1： 和算法2： 的误差时，精确解用双精111)(1+-=x x x y )1(1)(2+=x x x y 度的计算结果代替． 我们选取点集中的点作为x ，比较两种方法误差的差异．301}{=i i π 从图1.1可以看出，当x 不是很大时，两种算法的精度相当，但当x 很大时算法2的精度明显高于算法1． 这是因为，当x 很大时，和是相近数，用算法1进行计算时出x 111+x 现相近数相减，相同的有效数字相减后变成零，于是有效数字位数急剧减少，自然相对误差增大． 这一事实也可以从误差传播公式(1.12)分析出． 鉴于此，算法设计时，应该避免相近数相减．在图1.2中我们给出了当x 接近时，两种算法的精度比较，其中变量x 依次取为1-． 从图中可以看出两种方法的相对误差基本上都为，因而二者的精度相当．{}3011=--i iπ710-图1.1 算例1.1中两种算法的相对误差图（）+∞→x图1.2 算例1.1中两种算法的精度比较)1(-→x 算例1.2 试用不同位数的浮点数系统求解如下线性方程组⎩⎨⎧=+=+2321200001.02121x x x x 说明2：浮点数系统中的加减法在运算时，首先按较大的阶对齐，其次对尾数实施相应的加减法运算，最后规范化存入计算机．算法1 首先用第一个方程乘以适当的系数加至第二个方程，使得第二个方程的的系1x 数为零，这时可解出；其次将带入第一个方程，进而求得(在第三章中称该方法为高2x 2x 1x 斯消元法)． 当用4位和7位尾数的浮点运算实现该算法，分别记之为算法1a 和算法1b ． 算法 2 首先交换两个方程的位置，其次按算法1计算未知数 (第三章中称其为选主元的高斯消元法)． 当用4位和7位尾数的浮点运算实现该算法，分别记之为算法2a 和算法2b ．方程组的精确解为， ，用不同的算法计算出的...25000187.01=x ...49999874.02=x 结果见表1.1．表1.1 对算例1.2用不同算法的计算结果比较算例1.2*1x )(*1x r ε*2x )(*2x r ε算法1a 0.00000.10×1010.50000.25×10－7算法2a 0.25000.75×10－70.50000.25×10－7算法1b 0.26000000.40×10－10.49999870.10×10－6算法2b0.25000200.50×10－80.50000000.25×10－7对于算例1.2，表中的数据表明，当用4位尾数计算时，算法1给出错误的结果，算法2则给出解很好的近似． 这是因为在实现算法1时，需要给第一个方程乘以加00001.0/2－至第二个方程，从而削去第二个方程中的系数，但在计算的系数时需做如下运算1x 2x(1.13)661610000003.0104.0103.0104.03200001.02⨯⨯⨯⨯=+⨯＋＝－＋－－对上式用4位尾数进行计算，其结果为． 因为舍入误差，给相对较大的数加以6104.0⨯－相对较小的数时，出现大数“吃掉”小数的现象． 计算右端项时，需做如下运算(1.14)661610000002.0102.0102.0102.02100001.02⨯⨯⨯⨯=+⨯＋＝－＋－－同样出现了大数吃小数现象，其结果为． 这样，得到的变形方程组6102.0⨯－⎩⎨⎧⨯-=⨯-⨯=⨯+⨯62612114102.0104.0101.0102.0101.0x x x 中没有原方程组中第二个方程的信息，因而其解远偏离于原方程组的解． 该算法中之所以出现较大数的原因是因为运算，因而算法设计中尽可能避免用绝对值较大的数00001.0/2－除以绝对值较小的数． 其实当分子的量级远远大于分母的量级时，除法运算还会导致溢出，计算机终止运行．虽从单纯的一步计算来看，大数吃掉小数，只是精度有所损失，但多次的大数吃小数，累计起来可能带来巨大的误差，甚至导致错误． 例如在算法1a 中出现了两次大数吃小数现象，带来严重的后果． 因而尽可能避免大数吃小数的出现在算法设计中也是非常必要的． 当用较多的尾数位数进行计算，舍入误差减小，算法1和2的结果都有所改善，算法1的改进幅度更大些．算例1.3 计算积分有递推公式，已知⎰+=1055dx x x I n ),2,1(511 =-=-n I nI n n ． 采用IEEE 双精度浮点数，分别用如下两种算法计算的近似值．56ln 0=I 30I算法1 取的近似值为，按递推公式计算0I 6793950.18232155*0=I *1*51--=n n I nI *30I 算法2 因为，取的近似值为)139(5156)139(611039103939+⨯=<<=+⨯⎰⎰dx x I dx x 39I ，按递推公式计算3333330.004583332001240121*39≈⎪⎭⎫ ⎝⎛+=I ⎪⎭⎫ ⎝⎛-=-**1151n n I n I *30I 算法1和算法2 的计算结果见表1.2． 误差绝对值的对数图见图1.3．表1.2 算例1.3的计算结果算法1算法2n *nI n n I I -*n *nI nn I I -*18.8392e-002 1.9429e-01639 4.5833e-0033.9959e-0042 5.8039e-0029.8532e-016384.2115e-0037.9919e-0053 4.3139e-002 4.9197e-01537 4.4209e-003 1.5984e-0054 3.4306e-002 2.4605e-01436 4.5212e-003 3.1967e-0065 2.8468e-002 1.2304e-01335 4.6513e-003 6.3935e-0076 2.4325e-002 6.1520e-01334 4.7840e-003 1.2787e-007………33 4.9255e-003 2.5574e-00825 1.1740e+001 1.1734e+00132 5.0755e-003 5.1148e-00926-5.8664e+001 5.8670e+00131 5.2349e-003 1.0230e-00927 2.9336e+002 2.9335e+002 305.4046e-003 2.0459e-01028-1.4667e+003 1.4668e+003 297.3338e+0037.3338e+003 30-3.6669e+004 3.6669e+004图1.3 算例1.3用不同算法计算结果的误差绝对值的对数图从表1.2中的计算结果可以看出，算法1随着计算过程的推进，绝对误差几乎不断地以5的倍数增长，即有0*02*221*1*555I I I I I I I I n n n n n n n -≈≈-≈-≈----- 成立． 对于逐步向前推进的算法，若随着过程的进行，相对误差在不断增长，导致产生不可靠的结果，这种算法称之为数值不稳定的算法． 对于算法1绝对误差按5的幂次增长，但真值的绝对值却在不断变小且小于1，相对误差增长的速度快于5的幂次，导致产生错误的结果，因而算法1数值不稳定，不能使用． 而算法2随着计算过程的推进，绝对误差几乎不断地缩小为上一步的1/5，即有m m n m n n n n n n n I I I I I I I I 5/5/5/*22*21*1*++++++-≈≈-≈-≈- 成立． 绝对误差不断变小，真值的绝对值随着过程向前推进却在变大，这样相对误差也越来越小，这样的方法称之为数值稳定的算法． 算法1和算法2的误差对数示意图见图1.3． 这个算例告诉我们应该选用数值稳定的算法．知识结构图⎪⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎪⎨⎧⎪⎩⎪⎨⎧⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧⎩⎨⎧⎪⎩⎪⎨⎧⎩⎨⎧算法设计要点数值方法的稳定性数值方法的收敛性算法多元函数一元函数传播有效数字相对误差（限）绝对误差（限）度量截断误差舍入误差误差的产生误差误差与算法习题一1 已知有效数，，． 试给出各个近似值的绝对误105.3*1-=x 4*210125.0⨯=x 010.0*3=x 差限和相对误差限，并指出它们各有几位有效数字．2 证明当近似值是x 的较好近似时，计算相对误差的计算公式和相差一个*x x x x -***xxx -和同阶的无穷小量．2*⎪⎪⎭⎫⎝⎛-x x x 3 设x 的近似值具有如式(1.5)的表示形式，试证明*x 1)若具有n 位有效数字，则相对误差；*x n r a x e -⨯≤11*1021)(2)若相对误差，则至少具有n 位有效数字．n r a x e -⨯+≤11*10)1(21)(*x 4 试建立二元算术运算的绝对误差限传播近似计算公式．5 试建立如下表达式的相对误差限近似传播公式，并针对第1题中数据，求下列各近似值的相对误差限．1) ； 2) ； 3) *3*2*1*1x x x y +=3*2*2x y =*3*2*3/x x y =6若例题1.3中使用的尺子长度是80mm ，最小刻度为1mm ，量得某桌面长的近似值mm ，宽的近似值mm ． 试估计桌子长度、宽度的绝对误差限，并3.1304*=a 8.704*=b 求用该近似数据计算出的桌子面积的绝对误差限和相对误差限．7 改变如下计算公式，使其计算结果更为精确．1) 且0,cos 1≠-x xx1<<x 2)1,1ln )1ln()1(ln 1>>--++=⎰+N N N N N xdx N N3) 1,133>>-+x x x 8 (数值试验)试通过分析和数值试验两种手段，比较如下三种计算近似值算法的可靠性．1-e 算法1 ；∑=--≈mn nn e 01!)1( 算法2 ；101!1-=-⎪⎭⎫ ⎝⎛≈∑m n n e算法3 ；101)!(1-=-⎪⎪⎭⎫ ⎝⎛-≈∑m n n m e9 (数值试验)设某应用问题归结为如下递推计算公式 ，，72.280=y 251-=-n n y y,2,1=n 在计算时取为具有5位有效数字的有效数． 试分析近似计算公式的2*c **1*5c y y n n -=-绝对误差传播以及相对误差传播情况，并通过数值实验验证 (准确值可以用IEEE 双精度浮点运算结果代替)，该算法可靠可用吗？。
《数值分析》第1章
b
上两式作用得到：
4T ( h) − T ( 2h) = 3 I + O (h4 )
忽略高阶项得， I ≈ T (h) + (T (h) − T (2h)) . 公式的精度为 O (h4 ) .
1 3
此
其中 c1 , c2 ,L与 h 无关，则有，
19
20
§3 误差来源与误差分析的重要性
误差来源（或分类）
(1) 模型误差：建立数学模型时忽略一些次要 因素而引起的与真实情况的误差.
(2) 测量误差：数学模型中的一些已知参数， 由于受到测量工具或其它主观因素的影 响所带来的误差.
21
(3) 截断误差：数学模型常难以求解，往往要 用近似、易于求解的问题代替，这种简化 引起的误差.
P ( x ) = a0 x n + L + an −1 x + an 已知，对输入
的x，要计算P(x)的值，采取方法
u0 = 0 ⎧ t 1 = 1, ⎪ ⎨ t k = xt k − 1 , k = 2 , L , n ⎪u = u k = 1, L , n k −1 + a n− k tk , ⎩ k
29 30
例 15. 为使 20 的相对误差小于 0.1% ，要取几 位有效数字.
例 16. 用 3. 1416 表示π 的近似值，求其相对误 差？
解：因为 a1 = 3, n = 5 ，所以
er ( x ) ≤
1 1 × 10−5 + 1 = × 10−4 2× 3 6
解： 由 er ≤ 只需
1 × 10− n + 1 且 a1 = 4 ， 为使 er ≤ 0.1% ， 2a1
数值分析：第一章绪论PPT课件
x
*
是指对每一个 1 i
n
都有lim k
xi( k )
x
* i
可以。理解为 | |
x
(
k
)
x*
||
0
定义1.2.3
若存在常数
C1、C2
>
0
使得，
C1 || x ||B || x ||A C2 || x ||B
则称 || ·||A 和|| ·||B 等价。
可以理解为对任何
向量范数都成立。
数值分析课程中所讲述的各种数值方 法在科学与工程计算、信息科学、管理 科学、生命科学等交叉学科中有着广泛 的应用
第3页/共44页
应用问题举例
第4页/共44页
1、一个两千年前的例子
今有上禾三秉，中禾二秉，下禾一秉， 实三十九斗；
上禾二秉，中禾三秉，下禾一秉， 实三十四斗；
上禾一秉，中禾二秉，下禾三秉， 实二十六斗。 问上、中、下禾实一秉各几何？ 答曰：上禾一秉九斗四分斗之一。中禾 一秉四斗四分斗之一。下禾一秉二斗四 分斗之三。-------《九章算术》
定理1.2.1 Rn 上一切范数都等价。
第27页/共44页
二. 矩阵范数
定义1.2.4
Rmn空间的矩阵范数 || ·|| 对任意A, B R满mn足： (1) || A || 0 ; || A || 0 A 0 （正定性)
(2) || A || | | || A || 对任意 C (齐次性) (3) || A B || || A|| || B || （三角不等式)
1 1
(1
I1*
)
0.63
212056
第24页/共44页
我们仅仅是幸运吗?
数值分析第一章PPT课件
= f ’( )(x* x)
x* 与 x 非常接近时，可认为 f ’( ) f ’(x*) ，则有：
|e*(y)| | f ’(x*)|·|e*(x)|
即：x*产生的误差经过 f 作用后被放大/缩小了| f ’(x*)| 倍。故称| f ’(x*)|为放大因子 /* amplification factor */ 或 绝对条件数 /* absolute condition number */.
r* (x ) ln x * r* (y )
11 0n1lnx*0.1% 2a1
n4
.
10
1.3 避免误差危害的若干原则
算法的数值稳定性
用一个算法进行计算，如果初始数据误差在计算中 传播使计算结果的误差增长很快，这个算法就是数值不 稳定的.
.
11
1.3 避免误差危害的若干原则
病态问题与条件数
Cp
x f (x) f (x)
x nxn1 xn
n,
它表示相对误差可能放大 n倍.
如 n10，有 f(1 ) 1 ,f(1 .0)2 1 .2，4 若取 x 1, x*1.02, 自变量相对误差为 2% ，函数值相对误差为 24%， 这时问题可以认为是病态的.
一般情况下,条件数
Cp
10就认为是病态，
εr*21 a11 0n10.0 0% 1
已知 a1 = 3，则从以上不等式可解得 n > 6 log6，即
n 6，应取 * = 3.14159。
.
8
1.2 数值计算的误差
问题：对于y = f (x)，若用x* 取代x，将对y 产生什么影响？
分析：e*(y) = f (x*) f (x)
e*(x) = x* x
数值分析第5版第一章课件李庆扬著.
2
在 m相同的情况下， 越n 大则 10越m小n1，故有效位数越 多，绝对误差限越小.
x x * 1 10 . mn1 2
（2.2）
24
定理1 设近似数 x *表示为
x* 10m (a1 a2 101 al 10(l1) ), (2.1) 其中 ai (i 1,, l) 是0到9中的一个数字，a1 0, m为整数.
2
介于x, x * 之间，
取绝对值得 f (x) f (x*) f (x*) (x*) f ( ) 2 (x*).
2
30
假定 f (x*) 与 f (x*)的比值不太大，可忽略 (x*)
的高阶项，于是可得计算函数的误差限
( f (x*)) f (x*) (x*).
31
当 f 为多元函数，如计算 A f (x1, xn ) 时. 如果 x1,, xn 的近似值为 x1*,, xn*，则 A 的近似值为
需要注意的是误差限的大小并不能完全表示近似值的好坏.
13
例如，有两个量 x 10 1， y 1000 5, 则
x* 10,
* x
1;
y* 1000,
虽然
* y
比
*大
x
4
倍，但
* y
/
y*
5
/1000
0.5%
比
* y
5.
* x
/
x*
1 / 10
10%
要小得多，这说明 y *近似 y的程度比 x *近似 x的程度好.
Rn (x)
f (x) Pn (x)
f (n) ( ) xn1,
(n 1)!
在0与x之间.
有了计算公式后，在用计算机做数值计算时，还要受 计算机字长的限制，原始数据在计算机上表示会产生误差， 计算过程又可能产生新的误差，这种误差称为舍入误差.
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
4
装备一个铸造车间，需要熔炼设备、 造型及 制芯设 备、砂 处理设 备、铸 件清洗 设备以 及各种 运输机 械，通 风除尘 设备等 。只有 设备配 套，才 能形成 生产能 力。
误差来源
• 模型误差 • 方法误差 • 观测误差 • 舍入误差
2020/9/30
© Wuhan University Confidential
上述几种误差都会影响计算结果的精确性,因
而了解和研究这些误差对数值计算是有帮助
的.但是研究前俩种误差对计算结果的影响
往往不是计算工作者所能独立完成的,所以
我们一般只研究截断误差和舍入误差对计算
结果的影响.这俩种误差在数值计算中会产
生什么样的影响?这是我们这门课要重视的
问题.
We should think much
f (x)dx,
d f (x), ......
a
dx
近似解
2020/9/30
计算机
© Wuhan University Confidential
计算方法
方法 可行
3
装备一个铸造车间，需要熔炼设备、 造型及 制芯设 备、砂 处理设 备、铸 件清洗 设备以 及各种 运输机 械，通 风除尘 设备等 。只有 设备配 套，才 能形成 生产能 力。
• 方法可行性分析包含以下内容： 1.计算速度。 例 如,求解一个20 阶线性方
程组,用消元法需3000 次乘法运 算;而用克莱姆法则要进行20 10 7 . 9 ×次运算,如用每秒1 亿次乘法运算 的计算机要30 万年。
2.存储量。 大型问题有必要考虑。
3.精度。
2020/9/30
© Wuhan University Confidential
2020/9/30
© Wuhan University Confidential
6
装备一个铸造车间，需要熔炼设备、 造型及 制芯设 备、砂 处理设 备、铸 件清洗 设备以 及各种 运输机 械，通 风除尘 设备等 。只有 设备配 套，才 能形成 生产能 力。
来源二: 观测误差
• 测量工具精度与测量方法限制
of the influence of
these four errors!
2020/9/30
© Wuhan University Confidential
10
装备一个铸造车间，需要熔炼设备、 造型及 制芯设 备、砂 处理设 备、铸 件清洗 设备以 及各种 运输机 械，通 风除尘 设备等 。只有 设备配 套，才 能形成 生产能 力。
5
装备一个铸造车间，需要熔炼设备、 造型及 制芯设 备、砂 处理设 备、铸 件清洗 设备以 及各种 运输机 械，通 风除尘 设备等 。只有 设备配 套，才 能形成 生产能 力。
来源一 : 模型误差
• 模型误差:在建立数学模型过程中,不可能将 所有因素均考虑,必然要进行必要的简化,这 就带来了与实际问题的误差。
2020/9/30
© Wuhan University Confidential
2
装备一个铸造车间，需要熔炼设备、 造型及 制芯设 备、砂 处理设 备、铸 件清洗 设备以 及各种 运输机 械，通 风除尘 设备等 。只有 设备配 套，才 能形成 生产能 力。
•
提出问题
x,
ax,
lnx,
Ax b,
b
误差的表示方法
• 绝对误差绝对误差限 • 相对误差相对误差限 • 有效数字
2020/9/30
© Wuhan University Confidential
11
装备一个铸造车间，需要熔炼设备、 造型及 制芯设 备、砂 处理设 备、铸 件清洗 设备以 及各种 运输机 械，通 风除尘 设备等 。只有 设备配 套，才 能形成 生产能 力。
1
装备一个铸造车间，需要熔炼设备、 造型及 制芯设 备、砂 处理设 备、铸 件清洗 设备以 及各种 运输机 械，通 风除尘 设备等 。只有 设备配 套，才 能形成 生产能 力。
提问：数值分析是做什么用的？
这门课程的主要内容是研究使用计算机 求解各种数学问题的数值方法，(要求 方法能在计算机上实行，计算机只能 做加减乘除和逻辑运算)对求得解的精 度进行评估，以及如何在计算机上实 现求解等。
绝对误差（限）
定义：准确值 x, 近似值 x*, 则e*= x*- x 称为近似 值x* 的绝对误差，简称误差。
注意：误差可正可负，有量纲（有效数字后面举例）。
定义： |x*-x|≤ε*，ε*绝对误差限。
例: e-2.718=0.0002818…≤0.0003
2020/9/30
© Wuhan University Confidential
装备一个铸造车间，需要熔炼设备、 造型及 制芯设 备、砂 处理设 备、铸 件清洗 设备以 及各种 运输机 械，通 风除尘 设备等 。只有 设备配 套，才 能形成 生产能 力。
第一章 基本知识（一）
一 绪论 二 误差的来源 三 误差危害的防止
2020/9/30
© Wuhan University Confidential
Axb xk1Bkxg
ex
1xx2
x3
2! 3!
若将前若干项的部分和作为函数值的近似公式, 由于以后各项都舍弃了,自然产生uhan University Confidential
8
来源四: 装备一个铸造车间，需要熔炼设备、造型及制芯设备、砂处理设备、铸件清洗设备以及各种运输机械，通风除尘设备等。只有设备配套，才能形成生产能力。
2020/9/30
© Wuhan University Confidential
7
装备一个铸造车间，需要熔炼设备、 造型及 制芯设 备、砂 处理设 备、铸 件清洗 设备以 及各种 运输机 械，通 风除尘 设备等 。只有 设备配 套，才 能形成 生产能 力。
来源三:方法误差（截断误差)
• 计算方法本身的原因
舍入误差
• 计算机长有限
3.14159
• 注意:少量运算的舍入误差一般是微不足道的,但是 在计算 机上完成千白万次运算后
误差的积累很惊人.
2020/9/30
© Wuhan University Confidential
9
装备一个铸造车间，需要熔炼设备、 造型及 制芯设 备、砂 处理设 备、铸 件清洗 设备以 及各种 运输机 械，通 风除尘 设备等 。只有 设备配 套，才 能形成 生产能 力。