现代数值计算PDF电子书下载

第1章　科学计算与Matlab 1

1.1　科学计算的意义 1

1.2　误差基础知识 2

1.2.1　误差的来源 2

1.2.2　误差度量 2

1.2.3　有效数字 3

1.2.4 计算机的浮点数系 3

1.2.5　一个实例 4

1.2.6　数值计算中应注意的几个问题 4

1.3　Matlab软件 7

1.3.1　简介 7

1.3.2　向量和矩阵的基本运算 8

1.3.3　流程控制 15

1.3.4　脚本文件和函数文件 18

1.3.5　帮助系统 22

1.3.6 图功能 26

1.3.7　数据操作 30

习题一 33

数值实验一 33

第2章　线性方程组的直接解法 35

2.1　高斯消去法 35

2.2　矩阵的三角分解 39

2.2.1　LU分解和LDU分解 39

2.2.2　乔列斯基分解 42

2.2.3　追赶法 44

2.2.4　分块三角分解 46

2.3　QR分解和奇异值分解 47

2.3.1　正交矩阵 47

2.3.2　QR分解 50

2.3.3　奇异值分解 52

习题二 53

数值实验二 54

第3章　多项式插值与样条插值 56

3.1　多项式插值 56

3.1.1 多项式插值问题的定义 56

3.1.2　插值多项式的存在唯一性 57

3.1.3　插值基函数 57

3.2　拉格朗日插值 58

3.2.1　拉格朗日插值基函数 58

3.2.2　拉格朗日插值多项式 58

3.2.3　插值余项 60

3.3 牛顿插值 62

3.3.1　差商 62

3.3.2 牛顿插值公式及其余项 64

3.3.3　差分与等距节点的插值公式 65

3.4　埃尔米特插值 66

3.4.1　两点三次埃尔米特插值 66

3.4.2　埃尔米特插值多项式的余项 68

3.4.3　n+1个点2n+1次埃尔米特插值多项式H2n+1（x）及其余项R2n+1（x） 68

3.5　三次样条插值 68

3.5.1　样条插值概念的产生 68

3.5.2　三次样条函数 70

习题三 78

数值实验三 79

第4章　函数逼近 81

4.1　内积与正交多项式 81

4.1.1　权函数和内积 81

4.1.2　正交函数系 82

4.1.3　勒让德多项式 82

4.1.4　切比雪夫多项式 84

4.1.5　其他正交多项式 85

4.2　最佳一致逼近与切比雪夫展开 86

4.2.1　最佳一致逼近多项式 86

4.2.2　线性最佳逼近多项式的求法 87

4.2.3　切比雪夫展开与近似最佳逼近多项式 88

4.3　最佳平方逼近 90

4.3.1　预备知识 90

4.3.2　最佳平方逼近 90

4.4　曲线拟合的最小二乘法 94

4.4.1　最小二乘法 94

4.4.2　利用正交多项式作最小二乘拟合 97

4.4.3　非线性最小二乘问题 99

4.4.4　矛盾方程组 102

4.5　周期函数逼近与快速傅里叶变换 103

4.5.1　周期函数的最佳平方逼近 103

4.5.2　快速傅里叶变换（FFT） 105

习题四 107

数值实验四 108

第5章　数值积分与数值微分 109

5.1　引言 109

5.2　几个常用积分公式及其复合公式 110

5.2.1　几个常用积分公式 110

5.2.2　代数精度 111

5.2.3　积分公式的复合 112

5.3　变步长方法与外推加速技术 118

5.3.1　变步长梯形法 118

5.3.2　外推加速技术与龙贝格求积方法 119

5.4 牛顿-科茨公式 121

5.5　高斯公式 123

5.5.1　高斯公式的定义及性质 123

5.5.2　常用高斯型公式 126

5.5.3　高斯型公式的应用 132

5.6　多重积分的计算 135

5.6.1　二重积分的计算 135

5.6.2　蒙特卡罗模拟求积法简介 138

5.7　数值微分 141

5.7.1　基于拉格朗日插值多项式的求导方法 141

5.7.2　基于样条函数的求导方法 144

习题五 147

数值实验五 149

第6章　线性方程组的迭代解法 151

6.1　范数和条件数 151

6.1.1　向量范数和矩阵范数 151

6.1.2　扰动分析和条件数 153

6.2　基本迭代法 155

6.2.1　雅可比迭代法 156

6.2.2　高斯-赛德尔迭代法 157

6.2.3　超松弛（SOR）迭代法 158

6.2.4　迭代的收敛性分析和误差估计 160

6.3　不定常迭代法 165

6.3.1　最速下降法 165

6.3.2　共轭梯度法 168

6.3.3　广义极小残量法 172

6.3.4　预处理技术 176

习题六 178

数值实验六 179

第7章　非线性方程求根 180

7.1　非线性方程求根的基本问题 180

7.2　二分法 182

7.3　不动点迭代方法 183

7.4　迭代加速 186

7.5 牛顿法 188

7.6　割线法 194

7.7　非线性方程组简介 196

7.8　非线性最小二乘问题 199

7.9　大范围求解方法 201

习题七 204

数值实验七 205

第8章　矩阵特征值与特征向量的计算 206

8.1　前言 206

8.2　幂方法 208

8.2.1　乘幂法 208

8.2.2　反幂法 212

8.2.3　结合原点平移的反幂法 213

8.3　QR方法 214

习题八 216

数值实验八 217

第9章　常微分方程初边值问题数值解 218

9.1　欧拉公式及其改进 218

9.1.1　欧拉公式 218

9.1.2　数值积分与多步法 220

9.1.3　预估校正格式 223

9.2　龙格-库塔公式 225

9.3　收敛性与稳定性 230

9.3.1　单步法的收敛性 230

9.3.2　单步法的稳定性 232

9.4　微分方程组和刚性问题 233

9.5　有限差分法 236

习题九 239

数值实验九 239

参考文献 241

索引 242

表3-1　一般插值数据表 56

表3-2　函数y=ln x数据表 59

表3-3　差商表 63

表3-4　函数y=ln x的差商表 65

表3-5 两点埃尔米特插值数据 66

表3-6　给定数据表 76

表3-7　例3.5.3的差商表 76

表4-1　切比雪夫展开多项式系数 89

表4-2　给定的数据表 95

表4-3　二次最小二乘拟合数据 96

表4-4　多项式拟合数据表 98

表4-5　未知函数数据表 99

表4-6　函数y=a sin bx拟合数据 100

表5-1　函数值表 109

表5-2　函数值表 115

表5-3　三种不同方法的计算结果 116

表5-4　区间折半法的计算结果 119

表5-5　龙贝格方法计算结果 120

表5-6　科茨系数 122

表5-7　勒让德高斯点及高斯系数 127

表5-8　部分拉盖尔高斯点及高斯系数 129

表5-9　埃尔米特高斯点及高斯系数 131

表5-10　函数f（x）数据表 143

表7-1　二分法计算结果 183

表7-2　x3-x-1=0的四个不同的迭代方法 184

表7-3 方程xlnx=1的不同迭代法 186

表7-4 加速方法计算效果 187

表7-5 艾特肯加速效果 188

表7-6 ？的不同近似值 191

表7-7 牛顿法和割线法的计算效果 195

表7-8　方程组的不同迭代效果比较 198

表8-1　乘幂法的计算结果 210

表8-2　改进乘幂法的计算结果 212

表8-3　反幂法的计算结果 213

表8-4　带原点平移的反幂法的计算结果 214

表9-1　欧拉法的计算结果 220

表9-2　梯形公式的计算结果 222

表9-3　改进欧拉法的计算结果 224

表9-4　龙格-库塔公式阶数和次数的关系 228

表9-5　标准四阶龙格-库塔方法的计算结果 229

表9-6　龙格-库塔方法的计算结果 235

表9-7　不同N情况下有限差分解的误差 238

图1-1　河渠的图形 4

图1-2　命令plot演示 26

图1-3　命令plot3演示 27

图1-4　命令surf演示 28

图1-5　命令surf和contour用于二元函数 28

图1-6　图像的标注 29

图2-1　Givens变换 48

图2-2　Householder变换 50

图3-1　多项式插值的几何意义 57

图3-2　龙格函数及其10阶拉格朗日插值多项式与分段线性多项式 69

图4-1　勒让德多项式的图形，n=0,1,2,3 84

图4-2 切比雪夫多项式的图形，n=0,1,2,3,4 85

图4-3 最佳一次逼近多项式的几何意义 88

图4-4 例4.4.4数据点的分布情况 98

图4-5 非线性最小二乘拟合 101

图5-1 三种不同方法的收敛速度的比较 116

图5-2 被积函数的光滑性与高斯公式收敛速度的关系 133

图5-3 节点数与高斯公式的误差的关系 134

图5-4 二重积分求积区域 135

图5-5 误差RN随节点N的变化规律 140

图7-1 函数f(x)=x sin x/x2+1的图像 180

图7-2　非线性方程求根的不同困难程度 181

图7-3　不同的序列收敛速度 182

图7-4　不动点迭代的局部收敛性 185

图7-5 牛顿法的几何意义 189

图7-6 函数f（x）=x2+sin 10x-1的图像 193

图7-7　割线法的几何意义 195

图7-8　三个方程同伦算法的解曲线 202

图9-1　欧拉法的几何意义 219

函数mysort.m 18

函数mysort2.m 19

函数mysort3.m 20

函数zhouchang.m 22

函数fib.m 22

函数fib.m（带注释） 24

函数tridiagsolver.m 45

函数lagrange.m 59

函数chashang.m 63

函数nlfit.m 100

函数nlfitb.m 101

函数fmid.m 113

函数ftrapz.m 113

函数fsimpson.m 114

函数coeflege.m 128

函数gauss_lege.m 128

函数gauss_laguerre.m 130

函数gauss_her.m 131

函数herval.m 132

函数jacobi.m 156

函数gs.m 157

函数sor.m 158

函数cg.m 170

函数gmres.m 174

函数newton.m 192

函数homo.m 202

函数eigIPower.m 211

函数odeEuler.m 219

函数odeIEuler.m 224

算法2.1.1 高斯消去法 36

算法2.2.4　杜利脱尔算法 41

算法2.2.5　克洛脱算法 41

算法2.2.7　乔列斯基算法 43

算法2.2.9　追赶法 44

算法5.3.1　区间折半法 118

算法6.2.1　雅可比迭代算法 156

算法6.2.2　高斯-赛德尔迭代算法 157

算法6.2.3　SOR.迭代算法 158

算法6.3.2　最速下降法 167

算法6.3.6　共轭梯度法 169

算法6.3.10　广义极小残量法 173

算法6.3.13　左预处理共轭梯度法 177

算法7.2.1　二分法 183

算法7.5.5 牛顿下山法 192

算法7.7.1　非线性方程组的GS方法 196

算法7.7.4　拟牛顿法 198

算法7.7.6　改进拟牛顿法 199

算法7.8.1　高斯-牛顿法 200

算法7.8.2　LM算法 201

算法8.2.1　原始乘幂法 209

算法8.2.3　改进乘幂法 211

算法8.2.5　反幂法 212

算法8.2.7　带原点平移的反幂法 213

算法8.3.1　QR方法 214

算法9.1.1　欧拉方法 219

算法9.1.3　梯形方法 221

算法9.1.7　改进欧拉方法 224

算法9.2.1 标准四阶四段龙格-库塔方法 229

算法9.5.1　有限差分法 237