数值分析
（ Numerical Analysis）
武汉大学数学与统计学院信息与计算科学系数值分析课程建设小组
教材 (Text Book)
数值计算方法 郑慧娆等 编著
（武汉大学出版社）
辅导教材 (Tutorial Text Book)
数值计算方法学习指导书 邹秀芬等 编著
（武汉大学出版社）
参考书目 (Reference)
Numerical Analysis:Mathematics of Scientific
Computing (Third Edition)
数值分析 （英文版 第 3版 ）
David Kincaid & Ward Cheney
（ 机械工业出版社 )
Numerical Analysis (Seventh Edition)
数值分析 （第七版 影印版）
Richard L,Burden & J,Douglas Faires
（ 高等教育出版社）
网络资源
1、武汉大学,数值分析,教学专题网站
http://jpkc.whu.edu.cn/jpkcsite/szfx
2、国外,数值分析,课程网站
http://math.fullerton.edu/mathews/numerical.html
基础知识和工具
微积分
线性代数
常微分方程
VC程序设计语言
Matlab数学软件学时
理论教学,95学时（ 5学分 )
实践教学,36学时 ( 2学分，由实验教师单独给成绩,具体的方法见数值分析实践教学大纲 ）
考试方法
1,期终闭卷考试占 70%;
2,平时成绩占 20％，包括作业和课堂回答问题；
3,创新成绩占 10％，根据课堂内容所进行的创新活动，如科技小论文、心得体会、对课程改革的建议等，以读书报告的形式提交两次,
数值分析能够做什么？
Introduction
研究使用计算机求解各种科学与工程计算问题的数值方法（近似方法），对求得的解的精度进行评估，以及如何在计算机上实现求解等。
数值分析课程中所讲述的各种数值方法在科学与工程计算、信息科学、管理科学、生命科学等交叉学科中有着广泛的应用应用问题举例
2632
3432
3923
zyx
zyx
zyx
今有上禾三秉，中禾二秉，下禾一秉，
实三十九斗；
上禾二秉，中禾三秉，下禾一秉，
实三十四斗；
上禾一秉，中禾二秉，下禾三秉，
实二十六斗。
问上、中、下禾实一秉各几何？
答曰：上禾一秉九斗四分斗之一。中禾一秉四斗四分斗之一。下禾一秉二斗四分斗之三。 -------,九章算术,
1、一个两千年前的例子
nnnn
n
n
aaa
aaa
aaa
21
22221
11211
bxA?
nn
b
b
b
x
x
x
2
1
2
1
本课程第二章的内容：
线性方程组的数值方法！
2、天体力学中的 Kepler方程
x是行星运动的轨道，它是时间 t 的函数
s i n 0,0 1x x t
本课程第三章的内容：
非线性方程的数值解法全球定位系统：
在地球的任何一个位置，至少可以同时收到 4颗以上卫星发射的信号
3、全球定位系统（ Global Positioning System,GPS)
0
2
4
6
8
0
5
10
0
2
4
6
8
图 7,8
H
e
i
g
h
t
S6
S3
S4
S2
S1
R
S5
N - S p o s i t i o n s
表示地球上一个接收点 R的当前位置，卫星 Si的位置为
，则得到下列非线性方程组
(,,,)x y z t
(,,,)i i i ix y z t
2 2 2
1 1 1 1
2 2 2
2 2 2 2
2 2 2
3 3 3 3
2 2 2
4 4 4 4
2 2 2
5 5 5 5
2 2 2
6 6 6 6
( ) ( ) ( ) ( t - t ) 0
( ) ( ) ( ) ( t - t ) 0
( ) ( ) ( ) ( t - t ) 0
( ) ( ) ( ) ( t - t ) 0
( ) ( ) ( ) ( t - t ) 0
( ) ( ) ( ) ( t - t )
x x y y z z c
x x y y z z c
x x y y z z c
x x y y z z c
x x y y z z c
x x y y z z
0c
本课程第三章的内容：
非线性方程组的数值方法
1 1 2
2 1 2
12
(,,) 0
(,,) 0
(,,) 0
n
n
nn
f x x x
f x x x
f x x x
( ) 0Fx?记为其中,,
nnF D R R
12(,,,)
T
nx x x x?
4、已经测得在某处海洋不同深度处的水温如下：
深度（ M） 466 741 950 1422 1634
水温（ oC） 7.04 4.28 3.40 2.54 2.13
根据这些数据，希望合理地估计出其它深度（如 500米，
600米,1000米 … ）处的水温本课程第四章的内容：插值法
5、用比较简单的函数代替复杂的函数误差为最小，即距离为最小
（在不同的度量意义下）
本课程第五章的内容：函数逼近
6、人口预测下面给出的是中国 1900
年到 2000年的人口数，
我们的目标是预测未来的人口数（数据量较大时）
1950 55196
1960 66207
1970 82992
1980 98705
1990 114333
2000 126743
432231 ttty
30/)1 97 9( ts
432231 sssy
本课程第六章的内容：曲线拟合
7、铝制波纹瓦的长度问题建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的,
假若要求波纹瓦长 4英尺,每个波纹的高度 (从中心线 )为 1英寸,且每个波纹以近似 2π 英寸为一个周期,求制做一块波纹瓦所需铝板的长度 L.
这个问题就是要求由函数 f(x)=sin x给定的曲线从 x=0到 x=48英寸间的 弧长 L.
由微积分学我们知道,所求的弧长可表示为,
dxxdxxfL 480 2480 2' )( co s1))((1
上述积分称为第二类椭圆积分,它不能用普通方法来计算,
本课程第七章的内容：数值积分
8、生物化学反应的例子
A,B,C是三种蛋白质，其反应如下：
1
2
3
a
a
a
AB
B B C B
B C A C
我们通过建模可以得到如下方程组
'1 1 1 3 2 3y a y a y y
1 (0 ) 1y?
'22 1 1 3 2 3 2 2y a y a y y a y2 (0) 0y?
'3y? 2
22ay
3 (0) 0y?
A：
B：
C：
本课程第八章的内容：
常微分方程的数值方法
xxG?T
1ex?T
G,Google Matrix,
“the world’s largest matrix computation”,
4,300,000,000
x,PageRank vector
“The $25,000,000,000 Eigenvector”
9,Google搜索引擎
London,England,Millennium ('Wobbly') Bridge (1998-2002,
Norman Foster and Partners and Arup Associates)
… the natural modes and frequencies of a structure are the
solution of an eigenvalue problem that is quadratic when
damping effects are included in the model,(F,Tisseur,K,
Meerbergen,The quadratic Eigenvalue Problem,SiREV 43,2000,
pp.235-286)
本课程第九章的内容：
矩阵特征值问题的数值方法用计算机解决实际问题的步骤建立数学模型选择数值方法编写程序上机计算结果数值分析的特点
1、方法是近似的；
2、与计算机不能分离：上机实习
（掌握一门语言,C语言或 Fortran语言，
会用一种数学软件,Matlab或 Mathematica，
Maple）
在我们今后的讨论中,误差 将不可回避，
上机实习是需要大家创造条件完成的
1.2 误差 /* Error */
§ 1 误差的背景介绍 /* Introduction */
1,来源与分类 /* Source & Classification */
从实际问题中抽象出数学模型
—— 模型误差 /* Modeling Error */
通过测量得到模型中参数的值
—— 观测误差 /* Measurement Error */
求近似解 —— 方法误差 (截断误差 Truncation Error）
机器字长有限 —— 舍入误差 /* Roundoff Error */
§ 1.2.4 误差与有效数字
(Error and Significant Digits)
绝对误差 /* absolute error */
**e x x 其中 x
*为精确值,x为 x*的近似值。
10 006074302,.dxe x
例如：
**xx ε工程上常记为
|| *e *ε的上限记为,称为 绝对误差限 /* accuracy */，
相对误差 ( relative error ) **
*r
ee
x?
*
*
||r
εε
x
x 的 相对误差上限 定义为
有效数字 （ significant digits ）
用科学计数法，记 (其中 ）
若 （即 的截取按四舍五入规则），则称 为有 n 位有效数字，精确到 。
120 1 0 mnx,a a a
01?a
nm.xx 1050|| * na
x nm?10
1 4 1 5.3*;8 9 7 9 3 21 4 1 5 9 2 6 5 3 5.3
例：
问,有几位有效数字？请证明你的结论。*?
1
3 1 4
0 31 41 5 10
a nd * 0 5 10 0 5 10
*
π *,,
| π π |,,
证明：
有 4 位有效数字，精确到小数点后第 3 位。
（ Numerical Analysis）
武汉大学数学与统计学院信息与计算科学系数值分析课程建设小组
教材 (Text Book)
数值计算方法 郑慧娆等 编著
（武汉大学出版社）
辅导教材 (Tutorial Text Book)
数值计算方法学习指导书 邹秀芬等 编著
（武汉大学出版社）
参考书目 (Reference)
Numerical Analysis:Mathematics of Scientific
Computing (Third Edition)
数值分析 （英文版 第 3版 ）
David Kincaid & Ward Cheney
（ 机械工业出版社 )
Numerical Analysis (Seventh Edition)
数值分析 （第七版 影印版）
Richard L,Burden & J,Douglas Faires
（ 高等教育出版社）
网络资源
1、武汉大学,数值分析,教学专题网站
http://jpkc.whu.edu.cn/jpkcsite/szfx
2、国外,数值分析,课程网站
http://math.fullerton.edu/mathews/numerical.html
基础知识和工具
微积分
线性代数
常微分方程
VC程序设计语言
Matlab数学软件学时
理论教学,95学时（ 5学分 )
实践教学,36学时 ( 2学分，由实验教师单独给成绩,具体的方法见数值分析实践教学大纲 ）
考试方法
1,期终闭卷考试占 70%;
2,平时成绩占 20％，包括作业和课堂回答问题；
3,创新成绩占 10％，根据课堂内容所进行的创新活动，如科技小论文、心得体会、对课程改革的建议等，以读书报告的形式提交两次,
数值分析能够做什么？
Introduction
研究使用计算机求解各种科学与工程计算问题的数值方法（近似方法），对求得的解的精度进行评估，以及如何在计算机上实现求解等。
数值分析课程中所讲述的各种数值方法在科学与工程计算、信息科学、管理科学、生命科学等交叉学科中有着广泛的应用应用问题举例
2632
3432
3923
zyx
zyx
zyx
今有上禾三秉，中禾二秉，下禾一秉，
实三十九斗；
上禾二秉，中禾三秉，下禾一秉，
实三十四斗；
上禾一秉，中禾二秉，下禾三秉，
实二十六斗。
问上、中、下禾实一秉各几何？
答曰：上禾一秉九斗四分斗之一。中禾一秉四斗四分斗之一。下禾一秉二斗四分斗之三。 -------,九章算术,
1、一个两千年前的例子
nnnn
n
n
aaa
aaa
aaa
21
22221
11211
bxA?
nn
b
b
b
x
x
x
2
1
2
1
本课程第二章的内容：
线性方程组的数值方法！
2、天体力学中的 Kepler方程
x是行星运动的轨道，它是时间 t 的函数
s i n 0,0 1x x t
本课程第三章的内容：
非线性方程的数值解法全球定位系统：
在地球的任何一个位置，至少可以同时收到 4颗以上卫星发射的信号
3、全球定位系统（ Global Positioning System,GPS)
0
2
4
6
8
0
5
10
0
2
4
6
8
图 7,8
H
e
i
g
h
t
S6
S3
S4
S2
S1
R
S5
N - S p o s i t i o n s
表示地球上一个接收点 R的当前位置，卫星 Si的位置为
，则得到下列非线性方程组
(,,,)x y z t
(,,,)i i i ix y z t
2 2 2
1 1 1 1
2 2 2
2 2 2 2
2 2 2
3 3 3 3
2 2 2
4 4 4 4
2 2 2
5 5 5 5
2 2 2
6 6 6 6
( ) ( ) ( ) ( t - t ) 0
( ) ( ) ( ) ( t - t ) 0
( ) ( ) ( ) ( t - t ) 0
( ) ( ) ( ) ( t - t ) 0
( ) ( ) ( ) ( t - t ) 0
( ) ( ) ( ) ( t - t )
x x y y z z c
x x y y z z c
x x y y z z c
x x y y z z c
x x y y z z c
x x y y z z
0c
本课程第三章的内容：
非线性方程组的数值方法
1 1 2
2 1 2
12
(,,) 0
(,,) 0
(,,) 0
n
n
nn
f x x x
f x x x
f x x x
( ) 0Fx?记为其中,,
nnF D R R
12(,,,)
T
nx x x x?
4、已经测得在某处海洋不同深度处的水温如下：
深度（ M） 466 741 950 1422 1634
水温（ oC） 7.04 4.28 3.40 2.54 2.13
根据这些数据，希望合理地估计出其它深度（如 500米，
600米,1000米 … ）处的水温本课程第四章的内容：插值法
5、用比较简单的函数代替复杂的函数误差为最小，即距离为最小
（在不同的度量意义下）
本课程第五章的内容：函数逼近
6、人口预测下面给出的是中国 1900
年到 2000年的人口数，
我们的目标是预测未来的人口数（数据量较大时）
1950 55196
1960 66207
1970 82992
1980 98705
1990 114333
2000 126743
432231 ttty
30/)1 97 9( ts
432231 sssy
本课程第六章的内容：曲线拟合
7、铝制波纹瓦的长度问题建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的,
假若要求波纹瓦长 4英尺,每个波纹的高度 (从中心线 )为 1英寸,且每个波纹以近似 2π 英寸为一个周期,求制做一块波纹瓦所需铝板的长度 L.
这个问题就是要求由函数 f(x)=sin x给定的曲线从 x=0到 x=48英寸间的 弧长 L.
由微积分学我们知道,所求的弧长可表示为,
dxxdxxfL 480 2480 2' )( co s1))((1
上述积分称为第二类椭圆积分,它不能用普通方法来计算,
本课程第七章的内容：数值积分
8、生物化学反应的例子
A,B,C是三种蛋白质，其反应如下：
1
2
3
a
a
a
AB
B B C B
B C A C
我们通过建模可以得到如下方程组
'1 1 1 3 2 3y a y a y y
1 (0 ) 1y?
'22 1 1 3 2 3 2 2y a y a y y a y2 (0) 0y?
'3y? 2
22ay
3 (0) 0y?
A：
B：
C：
本课程第八章的内容：
常微分方程的数值方法
xxG?T
1ex?T
G,Google Matrix,
“the world’s largest matrix computation”,
4,300,000,000
x,PageRank vector
“The $25,000,000,000 Eigenvector”
9,Google搜索引擎
London,England,Millennium ('Wobbly') Bridge (1998-2002,
Norman Foster and Partners and Arup Associates)
… the natural modes and frequencies of a structure are the
solution of an eigenvalue problem that is quadratic when
damping effects are included in the model,(F,Tisseur,K,
Meerbergen,The quadratic Eigenvalue Problem,SiREV 43,2000,
pp.235-286)
本课程第九章的内容：
矩阵特征值问题的数值方法用计算机解决实际问题的步骤建立数学模型选择数值方法编写程序上机计算结果数值分析的特点
1、方法是近似的；
2、与计算机不能分离：上机实习
（掌握一门语言,C语言或 Fortran语言，
会用一种数学软件,Matlab或 Mathematica，
Maple）
在我们今后的讨论中,误差 将不可回避，
上机实习是需要大家创造条件完成的
1.2 误差 /* Error */
§ 1 误差的背景介绍 /* Introduction */
1,来源与分类 /* Source & Classification */
从实际问题中抽象出数学模型
—— 模型误差 /* Modeling Error */
通过测量得到模型中参数的值
—— 观测误差 /* Measurement Error */
求近似解 —— 方法误差 (截断误差 Truncation Error）
机器字长有限 —— 舍入误差 /* Roundoff Error */
§ 1.2.4 误差与有效数字
(Error and Significant Digits)
绝对误差 /* absolute error */
**e x x 其中 x
*为精确值,x为 x*的近似值。
10 006074302,.dxe x
例如：
**xx ε工程上常记为
|| *e *ε的上限记为,称为 绝对误差限 /* accuracy */，
相对误差 ( relative error ) **
*r
ee
x?
*
*
||r
εε
x
x 的 相对误差上限 定义为
有效数字 （ significant digits ）
用科学计数法，记 (其中 ）
若 （即 的截取按四舍五入规则），则称 为有 n 位有效数字，精确到 。
120 1 0 mnx,a a a
01?a
nm.xx 1050|| * na
x nm?10
1 4 1 5.3*;8 9 7 9 3 21 4 1 5 9 2 6 5 3 5.3
例：
问,有几位有效数字？请证明你的结论。*?
1
3 1 4
0 31 41 5 10
a nd * 0 5 10 0 5 10
*
π *,,
| π π |,,
证明：
有 4 位有效数字，精确到小数点后第 3 位。
