《数值计算方法及其程序实现》由编著者多年以来承担的暨南大学物理系硕士研究生必修课"数值计算方法"的讲授内容汇集而成,其内容包括七个部分:绪论、误差和数据处理、线性方程组的数值解法、非线性方程(组)的数值解法、数值积分与微分、常微分方程(组)的数值解法、偏微分方程的数值解法。这些内容通过例题分多个步骤予以展现。首先简要介绍数
值计算的基本方法和理论,再给出实现数值计算的逻辑流程构建,进而在Fortran和Matlab环境下编制计算程序,并分别VisualFortran6.0及Matlab6.5环境下运行,最终获得数值计算结果及其图示,同时提供了Fortran和Matlab两种计算机语言编写的相关程序。本书可作为数值计算方法课程的教材或参考书,也可作为计算物理及其相关学科的基础参考书。
计算物理与理论物理、实验物理相互依存、相互补充,是当今物理学不可缺少的三大板块之一。计算物理是基于物理原理、结合数值计算方法、利用计算机进行计算模拟,从而解释物理现象、描述复杂的物理过程、预测可能的物理结果的学科。计算物理的一个重要特点是利用计算机进行计算模拟(实验),得到物理过程中需要了解的中间信息记录,给出物理过程的详细结果描述。因此,数值计算方法及其程序实现是不可缺少的基础。
计算物理经近一个世纪的发展已趋于成型,不仅在物理及其相关学科中得到了广泛的应用,而且在丁程项目中也应用广泛。随着计算机的发展,计算物理在自然科学研究和社会经济发展中起着越来越重要的作用,因此迫切需要培养计算物理专业的人才。这些人才将以数学建模和数值计算为方法,以计算机高级语言(Fortran、c、Matlab等)为工具,以编制计算模拟程序或运用通用计算模拟程序为手段,进行相关应用基础研究或大量数据处理的工作,为相关学科的研究、工程技术应用和社会服务等提供直接应用或可供参考的计算结果。因此,这些人才的培养,少不了数值计算方法及其程序实现基础课程的学习。
本书由编著者多年以来承担的暨南大学物理系硕士研究生必修课“数值计算方法”的讲授内容汇集而成,其内容包括七个部分:绪论、误差和数据处理、线性方程组的数值解法、非线性方程(组)的数值解法、数值积分与微分、常微分方程(组)的数值解法、偏微分方程的数值解法。这些内容通过例题分多个步骤予以展现。首先简要介绍数值计算的基本方法和理论,再给出实现数值计算的逻辑流程构建,进而在Forrar和Mat—lab环境下编制计算程序,并分别在VisualFoltran6.0及Matlab6.5环境下运行,最终获得数值计算结果及其图示,同时提供了Fortran和Matlab两种计算机语言编写的相关程序。本书可作为数值计算方法课程的教材或参考书,也可作为计算物理及其相关学科的基础参考书。
在本书的编写完成之际,我心怀感激:感谢四川大学、南京大学、中国科学技术大学、美国麻省理工学院给予的教育和成长的空间,使我有机会在人生的进程中自由欣赏无限的风景,不断深入认知和体会生命的意义;感谢西北核技术研究所提供的工作空间,让我有机会从1987年就开始了解和从事与计算物理相关的研究工作,不断增长与历练相关能力和才干,使我热爱计算物理的程度与日俱增;感谢暨南大学提供的讲台,让我有机会将多年积累和掌握的计算物理相关基础知识与更多的人分享;感谢暨南大学出版社总编辑史小军教授对本书出版的大力支持。
1 绪论
1.1 数值计算方法及其技巧
1.2 计算物理简介
1.3 计算物理研究问题的方法和步骤
1.4 举例说明计算物理中数值问题解法
习题
2 误差和数据处理
2.1 测量数据的误差和分布
2.1.1 测量数据的误差
2.1.2 等精度测量数据的误差
2.1.3 非等精度测量数据的误差
2.1.4 测量数据的分布
2.1.5 应用实例
2.2 插值法
2.2.1 拉格朗日插值(Lagrange)
2.2.2 分段插值
2.2.3 二元函数插值
2.2.4 插值法在Matlab中的实现
2.3 最小二乘拟合
2.3.1 最小二乘原理
2.3.2 线性最小二乘拟合
2.3.3 直线最小二乘拟合
2.3.4 多项式最小二乘拟合
2.3.5 非线性函数最小二乘拟合
习题
3 线性方程组的数值解法
3.1 引言
3.2 直接解法
3.2.1 高斯(Gauss)消去法
3.2.2 高斯—约当(Gauss—Jordan)消去法
3.2.3 追赶法
3.3 迭代解法
3.3.1 雅可比(Jacobi)迭代法
3.3.2高斯一塞德尔(Gauss—Seidel)迭代法
习题
4非线性方程(组)的数值解法
4.1引言
4.2二分法
4.2.1确定有根区间
4.2.2 二分法
4.3 迭代法
4.3.1 不动点迭代法
4.3.2 牛顿(Newton)迭代法
习题
5 数值积分与微分
5.1 引言
5.2 等距节点求积公式
5.2.1 矩形求积公式
5.2.2 梯形求积公式
5.2.3 辛普森(Simpson)求积公式
5.2.4 牛顿—柯特斯(Newton—Cotes)求积公式
5.3 求积公式拓展
5.3.1 龙贝格(Romberg)求积公式
5.3.2 数值多重积分
5.4 数值微分
5.4.1 两点公式
5.4.2 三点公式
习题
6 常微分方程(组)的数值解法
6.1 常微分方程的离散化方法
6.2 一阶方程初值问题的数值解法
6.2.1 欧拉(Euler)方法和改进的欧拉方法
6.2.2 龙格—库塔(Runge—Kutta)方法
6.2.3 阿达姆斯(Adams)方法
6.3 一阶方程组和高阶方程的数值解法
6.3.1 一阶方程组的数值解法
6.3.2 高阶方程的数值解法
6.4 常微分方程边值问题的数值解法
6.4.1 化为初值问题的方法
6.4.2 边值问题的差分方法
习题
7 偏微分方程的数值解法
7.1 偏微分方程的离散化方法
7.1.1 偏微分方程的分类
7.1.2 偏导数的差分表示
7.2 拉普拉斯(Laplace)方程的差分解法
7.2.1 拉普拉斯方程的差分格式
7.2.2 特殊边界的处理
7.3 热传导方程的差分解法
7.3.1 显式、隐式差分格式
7.3.2 显隐交替差分格式
7.4 波动方程的差分解法
7.4.1 显式、隐式差分格式
7.4.2 显隐交替差分格式
习题
参考文献