本书介绍了科学与工程计算中常用的数值计算方法及相关理论。内容包括解线性方程组的直接法和迭代法、插值法、函数最优逼近、数值微积分、非线性方程(组)的迭代解法、矩阵特征值和特征向量的计算、常微分与偏微分方程数值解法等。其中包含了一些在实际中有重要应用的新方法,如求解超定方程组的最小二乘法、求解线性方程组的基于伽辽金原理的迭代法、奇异值分解、广义特征值问题的求解方法等。同时,对数值计算方法的计算效率、稳定性、收敛性、误差估计、适用范围及优缺点也进行了分析和介绍。
样章试读
目录
- 目录
前言
第1章 绪论 1
1.1 数值分析研究的内容与特点 1
1.2 误差 2
1.2.1 误差的来源与分类 2
1.2.2 绝对误差、相对误差与准确数字 2
1.2.3 计算机中数的表示与舍入误差 4
1.2.4 数据误差影响的估计 6
1.3 算法的数值稳定性 8
小结 9
习题 10
第2章 解线性方程组的直接法 12
2.1 高斯消去法 13
2.1.1 高斯消去法 13
2.1.2 高斯消去法中乘除法的运算量 18
2.1.3 高斯消去法顺利进行的条件 18
2.1.4 高斯消去法的算法组织 19
2.1.5 列主元高斯消去法 20
2.2 矩阵的三角分解 22
2.2.1 高斯消去法的矩阵形式 22
2.2.2 矩阵的LU分解 25
2.2.3 平方根法和改进平方根法 30
2.2.4 求解三对角方程组的追赶法 35
2.3 舍入误差对解的影响 37
2.3.1 向量范数与矩阵范数 37
2.3.2 舍入误差对解的影响 44
2.4 正交变换与矩阵的QR分解 49
2.4.1 吉文斯变换与豪斯霍尔德变换 49
2.4.2 矩阵的QR分解 52
2.5 超定方程组 65
2.5.1 线性最小二乘问题 65
2.5.2 最小二乘问题的求解 67
小结 71
习题 72
计算实习 75
第3章 解线性方程组的迭代法 77
3.1 向量序列和矩阵序列的极限 77
3.2 解线性方程组的基本迭代法 78
3.2.1 迭代法的一般格式 78
3.2.2 三种基本迭代法 78
3.3 迭代法的收敛性 83
3.3.1 迭代法的矩阵表示 83
3.3.2 迭代法的收敛性 84
3.4 共轭梯度法 92
3.4.1 求解线性方程组与求解二次函数极小点的等价性 92
3.4.2 共轭梯度法 93
3.5 基于伽辽金原理的迭代法 100
3.5.1 伽辽金原理和克雷洛夫子空间 100
3.5.2 阿诺尔迪过程 101
3.5.3 阿诺尔迪算法 103
3.5.4 广义极小残余算法 106
小结 110
习题 111
计算实习 113
第4章 插值法 115
4.1 多项式插值问题 115
4.2 拉格朗日插值多项式 118
4.3 牛顿插值多项式 120
4.3.1 差商的定义 121
4.3.2 牛顿插值多项式 121
4.3.3 差商的性质 124
4.4 埃尔米特插值多项式 125
4.5 分段低次插值多项式 129
4.5.1 高次插值多项式的缺陷 129
4.5.2 分段低次插值法 130
4.6 三次样条插值函数 132
4.6.1 三次样条插值函数的定义 132
4.6.2 三次样条插值函数的导出 132
4.6.3 三次样条插值函数的收敛性与误差估计 138
小结 138
习题 139
计算实习 141
第5章 函数最优逼近 142
5.1 函数的内积、范数和正交多项式 142
5.1.1 函数的内积和范数 142
5.1.2 正交多项式 144
5.2 最优平方逼近 151
5.2.1 最优平方逼近 151
5.2.2 正规方程组 152
5.3 最优一致逼近 162
5.3.1 最优一致逼近多项式 162
5.3.2 近似最优一致逼近多项式 165
小结 173
习题 174
计算实习 176
第6章 数值积分与数值微分 177
6.1 牛顿-科茨求积公式 177
6.1.1 数值积分的基本思想 177
6.1.2 牛顿-科茨求积公式 178
6.1.3 复化求积公式 180
6.1.4 变步长积分法 183
6.1.5 龙贝格积分法 184
6.2 待定系数法与高斯型求积公式 187
6.2.1 代数精度与待定系数法 187
6.2.2 广义佩亚诺定理 189
6.2.3 高斯型求积公式 191
6.2.4 常用的4种高斯型求积公式 197
6.3 数值积分的稳定性 201
6.4 数值微分 201
6.4.1 插值型数值微分公式 202
6.4.2 待定系数法 204
6.4.3 外推求导法 205
6.4.4 利用三次样条插值函数求导法 208
小结 208
习题 209
计算实习 210
第7章 非线性方程(组)的迭代解法 212
7.1 求解非线性方程的迭代法 212
7.1.1 几种基本迭代法 212
7.1.2 迭代法的收敛性 219
7.1.3 迭代法的收敛速度 224
7.1.4 加速收敛技术 226
7.2 求解非线性代数方程组的迭代法 228
7.2.1 简单迭代法 229
7.2.2 牛顿法 231
7.2.3 弦割法 234
7.2.4 布洛依登法 236
小结 237
习题 238
计算实习 239
第8章 矩阵特征值与特征向量的计算 241
8.1 基本性质 241
8.2 求一般矩阵特征值的计算方法 242
8.2.1 乘幂法及反幂法 242
8.2.2 求矩阵全部特征值与特征向量的QR 方法 245
8.2.3 阿诺尔迪方法 251
8.3 求实对称矩阵特征值的计算方法 253
8.3.1 雅可比方法 253
8.3.2 吉文斯方法 256
8.3.3 兰乔斯方法 258
8.4 奇异值(SVD)的计算 259
8.5 广义特征值问题 261
8.5.1 广义Schur分解 261
8.5.2 对称正定矩阵的广义Schur分解 262
小结 262
习题 263
计算实习 263
第9章 常微分方程数值解法 265
9.1 初值问题常用数值解法的建立与使用 265
9.1.1 基本数值解法的建立与隐式法的求解 265
9.1.2 龙格-库塔法 273
9.1.3 待定系数法、预测—校正公式 278
9.2 数值解中误差的积累、数值方法的收敛性和绝对稳定性 282
9.2.1 数值解中误差的积累和数值方法的收敛性 282
9.2.2 绝对稳定性 286
9.3 一阶微分方程组与高阶方程的数值解法 289
9.3.1 一阶微分方程组 289
9.3.2 高阶常微分方程 291
9.4 边值问题的数值解法 293
9.4.1 有限差分法 293
9.4.2 打靶法 300
小结 302
习题 303
计算实习 304
第10章 偏微分方程的数值解法 305
10.1 椭圆型边值问题 305
10.1.1 差分方程的建立 305
10.1.2 差分解的误差估计与收敛性 307
10.1.3 一般二阶椭圆型方程边值问题 310
10.2 抛物型方程初、边值问题 310
10.2.1 差分方程的建立与求解 311
10.2.2 差分格式的稳定性 313
10.2.3 差分解的误差估计与收敛性 315
10.3 双曲型方程混合问题 316
10.3.1 一阶双曲型方程 316
10.3.2 一阶常系数双曲型方程组 317
10.3.3 二阶双曲型方程 318
10.4 有限元法 320
10.4.1 变分原理 320
10.4.2 伽辽金逼近解 323
10.4.3 单元及形状函数 324
10.4.4 有限元求解步骤 327
小结 329
习题 329
计算实习 332
参考文献 333