数值实验与计算机模拟PDF电子书下载

第1篇 引论 1

1.1 数值实验及其步骤 1

1.2 数值实验的重要性 1

1.3 误差的来源与误差的基本概念 3

1.3.1 误差的来源与类型 3

1.3.2 绝对误差与绝对误差限 4

1.3.3 相对误差与相对误差限 5

1.3.4 有效数字 5

1.4 数值计算中需要注意的问题 7

1.4.1 避免两个相近的数相减 7

1.4.2 防止大数“吃掉”小数 8

1.4.3 注意简化计算步骤，减少运算次数 9

1.5 程序设计中应该注意的问题 10

1.5.1 对程序设计的一些建议 10

1.5.2 程序设计中的一些经验 11

第2篇 数值实验 17

2.1 Fortran语言程序的调试运行过程 17

2.1.1 Fortran语言简介 17

2.1.2 Fortran程序在V—F环境下的调试运行过程 17

2.1.3 调试运行例程序 18

2.1.4 实验与思考 19

2.2 实验数据处理Ⅰ与插值法 20

2.2.1 实验目的 20

2.2.2 算法描述 20

2.2.3 实验内容 20

2.2.4 参考程序 20

2.3 实验数据处理Ⅱ与曲线拟合 25

2.3.1 实验目的 25

2.3.2 算法描述 25

2.3.3 实验内容 26

2.3.4 参考程序 26

2.4 解线性方程组的迭代法 30

2.4.1 实验目的 30

2.4.2 算法描述 30

2.4.3 实验内容 31

2.4.4 参考程序 31

2.5 解线性方程组的消去法 33

2.5.1 实验目的 33

2.5.2 算法描述 34

2.5.3 实验内容 34

2.5.4 参考程序 34

2.6 非线性方程求根 38

2.6.1 实验目的 38

2.6.2 算法步骤 38

2.6.3 实验内容与要求 38

2.6.4 参考程序 39

2.7 矩阵特征值的计算 43

2.7.1 实验目的 43

2.7.2 幂法算法描述 43

2.7.3 实验内容 44

2.7.4 参考程序 44

2.8 常微分方程（组）的数值解法 49

2.8.1 实验目的 49

2.8.2 算法描述 50

2.8.3 实验内容 50

2.8.4 参考程序 50

2.8.5 思考练习 56

2.9 随机数的产生与检验 56

2.9.1 实验目的 56

2.9.2 算法描述 56

2.9.3 实验内容 58

2.9.4 参考程序 58

2.9.5 思考练习 59

2.10 蒙特卡罗方法的应用 60

2.10.1 实验目的 60

2.10.2 算法描述 60

2.10.3 实验内容 64

2.10.4 参考程序 64

第3篇 研究专题 70

3.1 无规网格点方法求解LJ团簇的基态几何结构 70

3.1.1 Lennard-Jones（LJ）团簇简介 70

3.1.2 LJ团簇的势函数描述 70

3.1.3 无规格点方法、模型构造与优化 71

3.1.4 最速下降法 72

3.1.5 无规格点方法求解LJ团簇的基态几何结构 72

3.1.6 程序使用说明 72

3.1.7 源程序 73

3.2 蒙特卡罗模拟退火方法求解LJ团簇的基态几何结构 81

3.2.1 蒙特卡罗方法 81

3.2.2 Metropolis蒙特卡罗模拟退火方法 81

3.2.3 蒙特卡罗模拟退火方法求解LJ团簇的基态几何结构 81

3.2.4 程序使用说明 82

3.2.5 源程序 83

3.3 中子星星体结构的计算 88

3.3.1 星体结构微分方程（组） 88

3.3.2 广义相对论效应 89

3.3.3 物态方程及微分方程（组）的化简 89

3.3.4 中子星结构特性计算 90

3.3.5 参考程序 90

3.4 逾渗模型的研究 98

3.4.1 逾渗理论的描述 99

3.4.2 二维正方逾渗模型的临界点研究 101

3.4.3 三维正方逾渗模型的临界点研究 108

3.5 有限扩散凝聚模型（DLA） 116

3.5.1 DLA模型简介 117

3.5.2 DLA生长过程的M-C模拟 118

3.5.3 参考程序 122

3.6 低能离子注入植物种子的深度-浓度分布的模拟 129

3.6.1 有关离子注入的LSS理论 129

3.6.2 近似假设 130

3.6.3 参考程序 131

3.7 中子输运问题的M-C模拟 136

3.7.1 中子输送问题 136

3.7.2 中子穿透平板模型 136

3.7.3 直接法与加权法模拟 137

3.7.4 参考程序 139

3.8 地下水渗流问题中最紧凑存储的有限单元方法 145

3.8.1 地下水渗流运动的有限元方程 145

3.8.2 按结点集成有限元方程组及一维最紧凑存储的方法 147

3.8.3 非线性问题的处理 150

3.8.4 参考程序 150

附录1 FORTRAN语句表 154

附录2 FORTRAN函数表 156

参考文献 159