(完整版)统计专业数值分析上机实习.docx

上机实验 李雪娇 章社生 武汉理工大学统计系 二零一二年二月 前言 数值分析 (numerical analysis) 是研究分析用计算机求解数学计算问题的数值计算方 法及其理论的学科， 是数学的一个分支， 它以数字计算机求解数学问题的理论和方法为研究 对象。 运用数值分析解决问题的过程： 实际问题→数学模型→数值计算方法→程序设计→上机 计算求出结果。数值分析这门学科有如下特点： 1. 面向计算机； 2. 有可靠的理论分析；  3. 要有好的计算复杂性； 4. 要有数值实验； 5. 要对算法进行误差分析。其主要内容有插值法， 函数逼近，曲线拟合，数值积分，数值微分，解线性方程组的直接方法，解线性方程组的迭代法，非线性方程求根，常微分方程的数值解法。 数值计算离不开好的计算方法， 也离不开好的编程能力。 人类社会已经迈入电子计算机时代，熟练地运用计算机进行科学计算，是广大科技工作者的一项基本技能。本书基于MATLAB软件和 C 语言给出各种算法的程序代码，在理论知识和实际应用间建立起桥梁。本 书以一个学生的角度编写， 贴近大学生的思维方式， 能够激励学生独立思考问题和创新的意 识；把基础理论的研究、方法构造的思路与利用 MATLAB和 C 语言有机地结合起来。虽然现 在很多软件具有数值计算的能力， 但是使用者不能根据自己的实际情况来进行操作， 存在着 固有的局限性。 本书详细地给出各种算法的代码 （除特别说明外， 每节中的例题指令是独立 完整的，因此读者可以很容易地在自己机器上实践） ，读者可以根据自己的使用情况加以修 改，留给了读者巨大的自我创新空间。 本书利用算法编程分析、 算法语言实现与说明文字配合的方式， 旨在培养读者的科学精 神和创新思维习惯， 重视培养读者的收集处理信息的能力、 获取知识的能力和分析解决问题 的能力。 对于那些经过大学基本数学教程的读者来说，通过本书， MATLAB精良完善的计算指令， 自然易读的程序将使他们感悟“教程”数学的基础地位和局限性，看到从“理想化”简单算例通向科学研究和工程设计实际问题的一条途径。 本书参加编写的还有李雪娇，阮向奎，肖先国，张维，张亮，江斌，朱海清，陈晓静，曹亮成，高丽莎，李营，于建华，郑洋，陈攀，马超，郑晓亮，劳展杰。 限于编著者水平，书中难免有怒当或失误之处，敬请广大读者朋友批评指正。 目录 第一章 秦九韶算法 5 第二章 插值方法 11 2.1 泰勒插值 11 2.2 拉格朗日插值 14 2.3 埃尔米特插值 17 2.4 分段插值 20 2.5 牛顿插值公式 25 2.6 埃特金（ Aitken ）算法 32 2.7 三次埃乐米特插值公式 39 第三章 数值积分 46 3.1 数值积分各公式的导出 46 3.2 各积分公式的复化形式 47 3.3 龙贝格积分 55 3.4 两点高斯公式及其复化公式 59 3.5 三点导数公式 61 第四章 常微分方程的差分方法 67 4.1 一步欧拉方法 67 4.2 二步欧拉方法 69 4.3 隐式欧拉方法及其改进方法 70 4.4 龙格 - 库塔方法 73 4.5 四阶亚当姆斯方法 95 第五章 方程求根的方法 101 5.1 二分法 101 5.2 简单迭代法 103 5.3 牛顿法 105 5.4 弦截法 107 第六章 线性方程组求根的方法 109 6.1 雅克比迭代法 109 6.2 高斯 -- 塞德尔迭代法 110 6.3 高斯消元法 112 6.4 追赶法 114 6.5 矩阵行范数 117 6.6 矩阵的列范数 120 6.7 矩阵的 2 范数 122 6.8 矩阵的特征根根 123 6.9 主元消元法 125 6.10 高斯消元法 3 138 3 章引论； 章． 1 ）牛顿插值公式 -江斌； 2）三次埃尔米特插值 -郑洋； 3）三次样条函数； 4）直 线拟合 -郑晓亮； 5）多项式拟合 -郑晓亮； 6）多元回归； 2.4 埃特金算法 -肖先国 2 章． 1）龙贝格积分 - 江斌 , 陈晓静； 2 ）随机抽样积分法； 2 ）二维重积分计算； 3 ）三维重积分计算； 4 ）三点导数公式 - 高丽莎； 5 ）五点导数公式 - 高丽莎； 6 ）中点加速法求导数 - 高丽莎 3 章． 1）二阶亚当姆斯方法 - 张亮； 2 ）三阶亚当姆斯方法 - 张亮； 3 ）四阶亚当姆斯方 法- 张亮； 4 ）预报校正亚当姆斯方法 - 张亮； 5 ）高阶方程差分法； 6 ）方程组差分法 - 于建华； 7 ）二阶线性边值问题差分法； 4 章． F(x,y)=0 计算 5 章． 1）矩阵的行范数 -朱海清； 2）矩阵的列范数 -朱海清； 3）矩阵的 2-范数 -朱海清； 4） 最大特征根计算 -曹亮成； 5）；最大