计算方法课件_易大义主编
计算方法第一章数值计算方法.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
计算方法_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. 一元函数情形
计算方法ppt
再按消元公式计算 ai(j(3) i,j=3,…,n)。
然后每一步类似的都在右下角方阵中的第一列中选
主元,再经行对调, 将主元素换到右下脚方阵中左上
角的位置。再按下一步计算 ai(jk() i,j=k,…,n)
直至将方程组化成上三角方程组,再进行回代就可 求得解。
21
全主元消去法
a11 a12 a1n b1
+
a (1) 1n
xn
+
...
+
a(2) 2n
xn
=
b(1) 1
=
b(2) 2
a(2) 32
x2
+
...
+
a(2) 3n
xn
=
b(2) 3
......
a(2) n2
xn
+ ... +
a(2) nn
xn
=
b(2) n
A(2) x = b(2)
其中 a(2) ij
=
a (1) ij
− li1a1(1j)
右端列向量视为方程组A的第n+1列,直接对矩阵A
(指现在的n行,n+1列的增广矩阵)进行行初等变换,
将其变换为上三角形矩阵,从而回代求解得到方程组
的解。
15
Gauss消去法(续)
Remark2:可以统计出,如果A为n阶方阵,则Gauss顺
序消去法消去过程所需的乘除运算次数为
∑n−1 (n − k + (n − k)(n − k +1)) = n3
(3.1)
2
引言(续)
其矩阵表示形式为: AX = b
a11 a12 a1n
x1
数值计算方法的义内容
3. 避免大数吃小数
例7:用单精度计算 精确解为
算法1:利用求根公式
的根。
在计则:1 = 0.0000000001 1010,取单精度时就成为: 109+1=0.10000000 1010+0.00000000 1010=0.10000000 1010
·算法2:先解出 再利用
注:求和时从小到大相加,可使和的误差减小。
输出 x1, x2
结束
二、算法的优劣
计算量小
例:用行列式解法求解线性方程组: n阶方程组,要计算n + 1个n阶行列式的值, 总共需要做n! (n - 1) (n + 1) 次乘法运算。
n=20 需要运 算多少次?
存贮量少 逻辑结构简单
n=100?
§ 3 数值计算中的误差
一、 误差的背景介绍
1. 来源与分类
现算法代的研科究和学应用研正是究本的课程三的主大题 !支柱
计算数学
理论研科究学实验科学计算
21世纪信息社会的两个主要特征: “计算机无处不在” “数学无处不在”
21世纪信息社会对科技人才的要求: - - 会“用数学”解决实际问题 - - 会用计算机进行科学计算
建立数学模型
选取计算方法
编写上机程序
数值计算方法
内容和方法
§ 1 数值计算方法的意义、内容与方法
20 世纪最伟大的科学技术发明---计算机 计算机是对人脑的模拟,它强化了人的思维智能; 计算机的发展和应用,已不仅仅是一种科学技术 现象,而且成了一种政治、军事、经济和社会现象;
软件的核心就是算法。 算法犹如乐谱, 软件犹如CD盘片, 而硬件如同CD唱机。
数值代数:方程求根、线性方程组求解、 特征值的求解;
计算方法第二章ppt
当方程组的系数矩阵为非奇异 矩阵(即满秩矩阵)时,高斯 消元法可求得唯一解。
列主元高斯消元法
列主元高斯消元法的 基本思想
在高斯消元法的基础上,每次选取列 中绝对值最大的元素作为主元进行消 元,以避免出现小主元导致的误差放 大问题。
列主元高斯消元法的 步骤
首先选取第一列中绝对值最大的元素 作为主元,通过行交换将其移到第一 行第一列位置,然后进行高斯消元。 在后续的消元过程中,每次均选取当 前列中绝对值最大的元素作为主元进 行消元。
100%
数值解法
通过计算机求解常微分方程的近 似解的方法,主要包括欧拉方法 和龙格-库塔方法等。
80%
离散化与步长
将连续的时间或空间域离散化, 取离散点上的函数值作为近似解 ,步长是相邻离散点间的距离。
欧拉方法
显式欧拉法
一种简单的数值解法,通过前 一步的函数值及其导数来推算 下一步的函数值。
隐式欧拉法
通过求解一个非线性方程来得 到下一步的函数值,具有较高 的精度和稳定性。
改进欧拉法
结合显式欧拉法和隐式欧拉法 的优点,提高算法的精度和效 率。
龙格-库塔方法
龙格-库塔法基本思想
自适应步长龙格-库塔法
通过多步计算并利用泰勒级数展开式, 得到更高精度的近似解。
根据误差估计自动调整步长,实现精 度和计算效率的动态平衡。
标准四阶龙格-库塔法
一种常用的高精度数值解法,具有局 部截断误差为$O(h^5)$的优点。
常微分方程数值解法误差分析
局部截断误差
数值解法在单步计算中所产生的误差,可以通过泰勒级数展开式进行估计。
全局误差
数值解法在整个计算过程中所产生的累积误差,与算法稳定性、步长选择等因素有关。
计算方法课件
(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课件
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.2 误差与有效数字
x 定义1 设 x 是某量的准确值,* 是 x 的一个 近似值,则称 e* = x* − x为近似值的误差或绝 对误差。 * 的绝对值的上界,即满足 x* − x ≤ ε * 的 ε *, e 称为近似值 x* 的误差限。 误差与精确值的比值称为相对误差。即 * er = ( x* − x) / x ,如果 ( x* − x) / x ≤ ε r*,则 ε r 称 为相对误差限。 实际使用中以 er* = ( x* − x) / x*为相对误差。
*
ε r (s ) ≈
ε (s )
*
s
*
27 = * * ≈ = 0.31% ld 8800
ε (s )
*
1.3 误差定性分析与避免误差危害
前面讨论的误差限的方法是最坏情况 对于千万次的运算可用概率统计的方法 20世纪60年代后提出
向后误差分析法 区间分析法
目前尚无有效方法对误差做出定量分析 本节讨论:
数值分析
Numerical Analysis
主讲教师: 主讲教师: 郭策安
Instructor: GUO CEAN E-mail: guocean@
教材
(Text Book)
TUP & Springer Press
李庆扬、王能超、 李庆扬、王能超、易大义 编
数值分析
参考书目 (Reference)
In = e (x e
I0 = e
−1
−1
n x 1 0
− n ∫ x n −1e x dx) = 1 − nI n −1 (n = 1, 2,L)0来自1∫1
0
计算方法课件 1.4
计算方法第一章绪论§1.4 差分方程§1.4 差分方程,,(),i i i t a ih a b x x t 设为给定函数,记dx a,b C :i iIx x 1i i Ex x 差分方程:微分方程的离散形式1i i ix x x 在节点x i 处以h 为步长的一阶向前差分:一阶向后差分:1i i i x x x ::向前差分算子向后差分算子恒等算子位移算子I :E : 算子均是线性算子且可两两交换. 这些算子的和的幂可按二项式展开.,,I ,E a t 0t 1t 2...t i-1t it i+1...t N-1t N b0x x 12x ix 1N x -Nx -1i x 1i x +2差分算子的关系及运算11,.mm mm i i i i x x x x1 i i i i i i x x x Ex Ix E I x ,算子之间的关系: ,E I 1=. I E 高阶差分:00100 1=11=1=,,!.!!mmm j j m m ji i i m i j j j m m m m j m j m m ji i i i j m j j m m x E I x E x x j j m m x I E x E x x j j m m j j m j其中 差分的计算:3k 阶差分方程形如10 ;,,,,,,,kdn n n n n n k F n x x x x x x C且均显含于上式中的方程称为k 阶差分方程。
n n k x x , 110;,,,,,,,.dn n n k n n n k G n x x x x x x C 本课程主要涉及如下k 阶线性差分方程012 (1.19)(),,,,,k jn jn j n a n xb 其中系数00(),,()().j n k a n b C a n a n 且 n x .011k x x x ,,,,若给定其k 个初始值则由方程即可求出其解序列可等价地化为如下k 阶差分方程若则称方程为齐次的,否则称之为非齐次的。
计算方法课件第一章 绪论
在我们今后的讨论中,误差将不可回避,
上机实习是需要大家创造条件完成的
物联网工程学院
第一章 绪论
用计算机解决实际问题的步骤
建立数学模型 选择数值方法
编写程序 上机计算结果
结果分析
物联网工程学院
第一章 绪论
小结:用计算机求数学问题的数值解不是简单地构造算法,它涉
及多方面的理论,例如,算法的收敛性和稳定性等。除理论分析外,
物联网工程学院
第一章 绪论
课程学习结束后你具备的能力
1. 对具体的数值计算问题,你会选择合适的 算法,并通过计算机计算出正确结果; 2. 对给定的算法会从理论上分析其优劣性; 3. 会根据原理构造解决较简单数值计算问题 的算法。
构造性
方法的构造,解的存在唯一性的证明
递推性
复杂计算过程转化成简单的计算过程的多次重复 (适合计算机计算)
近似替代
在误差允许的范围内,无限次的计算用 有限次计算替代
物联网工程学院
ห้องสมุดไป่ตู้
第一章 绪论
特点:
1、方法是近似的;
2、与计算机不能分离:上机实习
(掌握一门语言:C语言或Fortran语言, 会用一种数学软件:Matlab 或 Mathematica ,Maple)
研究对象由数学模型提出求解的
数值计算方法并编程计算出结果, 然后进行误差分析。
内容
非线性方程求解 (Ch. 2) 线性方程组求解 (Ch.3) 曲线拟合、插值法 (Ch.5) 数值积分 (Ch.6) 常微分方程组 (Ch.7)
物联网工程学院
第一章 绪论
方法:
离散化
计算离散点上的近似值
《计算方法(1)》课程教学大纲
计算方法课程教学大纲(Calculation Method)一、课程概况课程代码:0821018学分:3学时:48(其中:讲授学时32 ,实验学时16 ,上机学时0)先修课程:数学分析,高等代数等适用专业:小学教育(理)专业建议教材:《计算方法》,易大义,浙江大学出版社,2017.5课程归口:理学院课程的性质与任务:本课程是小学教育(理)专业的一门重要基础课。
通过本课程的学习,使学生系统地获得计算方法的基本知识、必要的基础理论;提高学生的数学视野、数学思维能力、逻辑推理能力;提高学生的数学素养,为学生学习后续相关课程及终身学习奠定必要的数学基础。
二、课程目标目标1.能够获得课程基本概念与性质。
目标2. 能够掌握本课程要求的计算方法。
目标3. 能够具有一定的抽象概括、逻辑推理等能力。
目标4. 能够具有一定的运算能力。
目标5. 能够具有一定的数学思维与分析能力。
本课程支撑专业人才培养方案中毕业要求3-1、毕业要求3-2,毕业要求6-2对应关系如表所示。
三、课程内容及要求(一)数值计算的基本概念1.教学内容(1)能够了解数值计算的研究对象和内容(2)能够了解数值算法的基本概念(3)能够了解误差的基本理论(4)能够了解数值算法设计的若干原则2.基本要求(1)重点与难点:误差的计算。
(2)教学方法:启发式互动讲授结合多媒体辅助;适当课堂练习;及时了解学生的作业状况并对共同的问题作及时解答;安排好课后答疑。
3.思政内容注重理论联系实际,尊重客观规律,树立社会主义核心价值观,增强专业素养,强调理论对实践的指导意义。
(二)非线性方程的迭代法1.教学内容(1)能够了解二分法(2)能够掌握Picard迭代法(3)能够掌握牛顿型迭代法2.基本要求(1)重点与难点:Picard迭代法、牛顿型迭代法及其实现。
(2)教学方法:启发式互动讲授结合多媒体辅助;适当课堂练习;及时了解学生的作业状况并对共同的问题作及时解答;安排好课后答疑。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
则 p(x) = q(x),
以下我们分别作画
p(x) 与 q(x)的图。
18/47
计算方法——09计11、61
p(x)
q(x)
上例说明,即使数学上的恒等公式,用计 算机来算,结果也是不一样的。
计算方法——09计11、61
19/47
§1.2 误差的种类及其来源
一. 误差来源
实 际 问 题 ( 数 学 模 型
高等数学、线性代数、程序设计语言
计算方法——09计11、61
5/47
教材
《计算方法》,
易大义等,
浙江大学出版社, 2002年第2版
计算方法——09计11、61
6/47
主要参考书
1.《数值分析引论》,
易大义 陈道琦,
浙江大学出版社, 1998年第1版
计算方法——09计11、61
7/47
主要参考书
(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 * 因
总成绩:100%。
计算方法——09计11、61
4/47
课程基本信息
总学时 56
其中上课40学时;上机16学时。
主要讲授如何解决各种工程技术问题,用数学语 言描述问题,即建立数学模型,将之转化为一个 数学问题,寻求合适的近似计算方法,编程计算, 充分发挥计算机的记忆和快速运算功能,寻求最 佳方案。 本课程先修课程为:
插值法 数值逼近 本 课 程 的 内 容 数据拟合的最小二乘法 数值积分和数值微分* 线性方程组的求解 数值代数 非线性方程组的求解 矩阵特征值* 常微分方程的数值方法
计算方法——09计11、61
11/47
第一章
数值计算中的误差
3 学时
本章内容
§1.1 引言 §1.2 误差的种类及其来源 §1.3 绝对误差和相对误差 §1.4 有效数字及其与误差的关系 §1.5 误差的传播与估计 §1.6 选用算法应遵循的原则 小结 作业与实验
值班军官对连长: 根据营长的命令, 明晚8点哈雷彗
星将在操场上空出现。如果下雨的话, 就让士兵穿 着野战服列队前往礼堂, 这一罕见的现象将在那里 出现。
计算方法——09计11、61
24/47
§1.2 误差的种类及其来源
连长对排长:
根据营长的命令, 明晚 8 点, 非凡的哈 雷彗星将身穿野战服在礼堂中出现。如果操场上 下雨, 营长将下达另一个命令, 这种命令每隔 76 年 才会出现一次。
有限过程代替无限过程的误差(无穷级数
右端是截断误差。
x x si n x ( x ) ...... 3! 5!
计算方法——09计11、61
22/47
3
5
§1.2 误差的种类及其来源
4.
舍入误差 /* Roundoff Error */
计算机字长有限,一般实数不能精确存储,于 是产生舍入误差。 例如:在 10 位十进制数限制下:
25/47
班长对士兵:
§1.3 绝对误差和相对误差
一.绝对误差 /* absolute error */
设
称
x
——准确值,x * ——近似值。
*
*
e( x ) x x 为 x 的绝对误差(简称误差) | e( x ) | 为 x * 的绝对误差限。
二.相对误差 /* relative error */
注: 四舍五入规则修正为四舍五以上入, 五时若前一
位是偶数则5舍去, 奇数则进一, 以减少积累误差。
计算方法——09计11、61
28/47
§1.4 有效数字及其与误差的关系
例4:x 32.93,x * 32.89, = 1 x x * 0.04 0.05 101 2 1 x x * 102 3 2 故 x * 有 3 位有效数字,分别是 3,,。 2 8 由于 x * 中的数字 9 不是有效数字, 故 x * 不是有效数。
计算方法——09计11、61
13/47
本章要求
1. 熟悉计算方法在解决实际问题中所处的地位, 熟悉计算方法是以计算机为工具求近似解的数 值方法; 2. 熟悉绝对误差(限),相对误差(限)及有 效数字概念; 3. 熟悉公式;
4. 熟悉选用算法应遵循的原则。
计算方法——09计11、61
14/47
计算方法——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%)。
2.《数值分析基础教程》,
李庆杨,
高等教育出版社,
2001年第1版
计算方法——09计11、61
8/47
主要参考书
3.《数值方法和MATLAB实现与应用》, (美) Gerald Recktenwald 著 伍卫国 万群 张辉 等译, 机械工业出版社,
2004年第1版
其他各类有关 “数值分析” 和 “计算方法” 的书
有n位有效数字时, 0.5 10m n1 1 r 10 n1 a1 10m 2a1 x 反之,由 1 x r (a1 1) 10 10 n1 2(a1 1)
设 x (0.a1a 2 a n a p )10m (a1 0, p ) 1 若 x * x 10m n ( P10 1.4.2) 2 则说 x *具有n位有效数字,分别是 1 , a 2 , , a n a 若n=p,则称 x *为有效数。
计算方法
计算机学院
陈克建
诚信声明
本课件中大量采用网络及其他渠道搜集的相关文 字信息和图片信息,这些信息因数量较大,无法 一一列明出处,本人在此郑重声明:这些相关资 料的版权归原作者所有,本文引用仅仅用于教学 目的。如有不妥,请与本人联系,联系方式: xznuckj@。 在此对相关资料的作者所付出的辛勤劳动表示衷 心的感谢,并对作者表示诚挚敬意! 本人郑重承诺:尊重知识,尊重劳动,尊重版权, 学术诚信。
计算方法——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
排长对班长:
明晚 8 点, 营长将带着哈雷彗星在礼 堂中出现, 这是每隔 76 年才有的事。如果下雨的 话,营长将命令彗星穿上野战服到操场上去。 在明晚 8 点下雨的时候, 著名的 76 岁 哈雷将军将在营长的陪同下身着野战服, 开着他那 “彗星”牌汽车, 经过操场前往礼堂。
计算方法——09计11、61
§1.1 引言
解决科学技术和工程问题的步骤:
实际问题 数学问题 提供计算方法
程序设计
的数学模 型简化成一系列算术运算和逻辑运算, 以便 在计算机上求解, 并对算法的稳定性、收敛 性和误差进行分析。
计算方法——09计11、61
15/47
§1.1 引言
舍入误差很小,本课程将研究它在运算过程中 是否能有效控制。
计算方法——09计11、61
23/47
§1.2 误差的种类及其来源
据说, 美军1910年的一次部队的命令传递是这样的:
营长对值班军官: 明晚大约 8点钟左右, 哈雷彗星将
可能在这个地区看到, 这种彗星每隔 76年才能看见 一次。命令所有士兵着野战服在操场上集合, 我将 向他们解释这一罕见的现象。如果下雨的话, 就在 礼堂集合, 我为他们放一部有关彗星的影片。
e( x ) 称 er ( x ) 为 x * 的相对误差。 x * e( x ) 实用中,常用 表示 x 的相对误差。 x 称 er ( x ) 为 x * 的相对误差限。
计算方法——09计11、61
26/47
§1.4 有效数字及其与误差的关系
一.有效数字 /* significant digits */
例1,例2 的结果的根源 建 立 算 法 截方 断法 误误 差差 ) ( 上 机 计 算 舍 入 误 差
20/47
结
果
模 型 误 差
观 测 误 差
初 值 误 差
)
计算方法——09计11、61
§1.2 误差的种类及其来源
二. 误差分类
1.
模型误差(描述误差)/* Modeling Error */ 题之差。
计算方法——09计11、61
9/47
《计算方法》课程体系
第一章
数值计算中的误差
第二章
第三章 第四章 第五章 第六章 第七章
插值法
曲线拟合的最小二乘法 数值积分 非线性方程的数值解法 方程组的数值解法 常微分方程数值解法