数值计算基础实指导书.docVIP

数值计算基础 实验指导书 2012年 目录 实验一 直接法解线性方程组的 1 实验二 插值方法 4 实验三 数值积分 6 实验四 常微分方程的数值解 8 实验五 迭代法解线性方程组与非线性方程 10 实验一 直接法解线性方程组 一、实验目的 掌握列选主元消去法与追赶法解线性方程组。 二、实验内容 分别写出Guass列选主元消去法与追赶法的算法，编写程序上机调试出结果，要求所编程序适用于任何一解线性方程组问题，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。 1、用Guass列选主元消去法求解方程组 2、用追赶法求解方程组 三、实验仪器设备与材料 主流微型计算机 四、实验原理 1、Guass列选主元消去法 对于AX =B 1）、消元过程：将（A|B）进行变换为，其中是上三角矩阵。即： k从1到n-1 列选主元 选取第k列中绝对值最大元素作为主元。 换行 归一化 消元 2）、回代过程：由解出。 2、追赶法 线性方程组为： 做LU分解为： 分解公式： 则 回代公式： 五、实验步骤 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。 三、实验仪器设备与材料 主流微型计算机 四、实验原理 常微分方程的数值解主要采用“步进式”，即求解过程顺着节点排列次序一步一步向前推进，在单步法中改进欧拉法和四阶龙格-库塔法公式如下： 1、改进欧拉法 2、四阶龙格-库塔法 五、实验步骤