2014-2015-2数值计算基础实验指导书.docVIP

数值计算基础 实验指导书 2015年 目录 实验一 直接法解线性方程组的 1 实验二 插值方法 4 实验三 数值积分 6 实验四 常微分方程的数值解 8 实验五 迭代法解线性方程组与非线性方程（课外） 10 实验一 直接法解线性方程组 一、实验目的 掌握Guass列选主元消去法、三角分解法解线性方程组。 二、实验内容 分别写出Guass列选主元消去法、三角分解法的算法，编写程序上机调试出结果，要求所编程序适用于任何线性方程组问题，即能解决这一类问题，而不是某一个问题。 实验中以下列数据验证程序的正确性。 1、Guass列选主元消去法 2、Doolittle三角分解法 三、实验仪器设备与材料 主流微型计算机 四、实验原理 1、Guass列选主元消去法 对于AX =B 1）、消元过程：将（A|B）进行变换为，其中是上三角矩阵。即： k从1到n-1 列选主元 选取第k列中绝对值最大元素作为主元。 换行 归一化 消元 2）、回代过程：由解出。 2、三角分解法（Doolittle分解） 将A分解为如下形式 由矩阵乘法原理 a、计算U的第一行，再计算L的第一列 b、设已求出U的1至r-1行，L的1至r-1列。先计算U的第r行，再计算L的第r列。 a）计算U的r行 b）计算L的r列 五、实验步骤 1、理解并掌握列选主元消去法、三角分解法、平方根法与追赶法； 2、画出列选主元消去法、三角分解法、平方根法的流程图 3、使用C语言编写出相应的程序并调试验证通过 六、实验报告要求 1、统一使用《武汉科技大学实验报告》本书写，实验报告的内容要求有：实验目的、实验内容实验实验二 插值方法 一、实验目的 掌握拉格郎日插值法与牛顿插值法构造插值多项式。 二、实验内容 分别写出拉格郎日插值法与牛顿插值法的算法，编写程序上机调试出结果，要求所编程序适用于任何一组插值节点，即能解决这一类问题，而不是某一个问题。 实验中以下列数据验证程序的正确性。 已知下列函数表 xi 0.56160 0.56280 0.56401 0.56521 yi 0.82741 0.82659 0.82577 0.82495 求x=0.5635时的函数值。 三、实验仪器设备与材料 主流微型计算机 四、实验原理 已知n个插值节点的函数值，则可由拉格郎日插值公式与牛顿插值公式构造出插值多项式，从而由该插值多项式求出所要求点的函数值。拉格郎日插值公式与牛顿插值公式如下： 1、Lagrange插值公式 2、Newton插值公式 五、实验步骤 1、理解并掌握拉格郎日插值法与牛顿插值法的公式； 2、画出拉格郎日插值法与牛顿插值法算法的流程图； 3、使用C语言编写出相应的程序并调试验证通过。 六、实验报告要求 1、统一使用《武汉科技大学实验报告》本书写，实验报告的内容要求有：实验目的、实验内容实验实验三 数值积分 一、实验目的 掌握复化梯形法与龙贝格法计算定积分。 二、实验内容 分别写出变步长梯形法与Romberge法计算定积分的算法，编写程序上机调试出结果，要求所编程序适用于任何类型的定积分，即能解决这一类问题，而不是某一个问题。 实验中以下列数据验证程序的正确性。 求。 三、实验仪器设备与材料 主流微型计算机 四、实验原理 通过变步长梯形法与龙贝格法，我们只要知道已知n个求积节点的函数值，则可由相应的公式求出该函数的积分值，从而不需要求该函数的原函数。变步长梯形法与龙贝格法公式如下： 1、变步长梯形法 用来控制精度 2、龙贝格法 梯形法则 二阶公式 四阶公式 六阶公式 八阶公式 用来控制精度 五、实验步骤 1、理解并掌握变步长梯形法与龙贝格法的公式； 2、画出变步长梯形法与龙贝格法的流程图 3、使用C语言编写出相应的程序并调试验证通过 六、实验报告要求 1、统一使用《武汉科技大学实验报告》本书写，实验报告的内容要求有：实验目的、实验内容实验积分中，被积函数在x=0点函数值为1，对该点在程序设计中应注意对其的定义。 八、思考题 使用复化梯形法与复化Simpson法来计算该问题有何缺点？ 实验四 常微分方程的数值解 一、实验目的 掌握改进欧拉法与四阶龙格-库塔求解一阶常微分方程的初值问题。 二、实验内容 分别写出改进欧拉法与四阶龙格-库塔求解的算法，编写程序上机调试出结果，要求所编程序适用于任何一阶常微分方程的数值解问题，即能解决这一类问题，而不是某一个问题。 实验中以下列数据验证程序的正确性。 求步长h=0.25。 三、实验仪器设备与材料 主流微型计算机 四、实验原理 常微分方程的数值解主要采用“步进式”，即求解过程顺着节