数值计算基实验指导书孙谦.docVIP

HYPERLINK mailto:wust_cs@163.com wust_cs@163.com 123456 数值计算基础 实验指导书 2012年 目录 TOC \o 1-3 \h \z \u HYPERLINK \l _Toc182904057 实验一 直接法解线性方程组的 PAGEREF _Toc182904057 \h 1 HYPERLINK \l _Toc182904058 实验二 插值方法 PAGEREF _Toc182904058 \h 4 HYPERLINK \l _Toc182904059 实验三 数值积分 PAGEREF _Toc182904059 \h 6 HYPERLINK \l _Toc182904060 实验四 常微分方程的数值解 PAGEREF _Toc182904060 \h 8 HYPERLINK \l _Toc182904061 实验五 迭代法解线性方程组与非线性方程 PAGEREF _Toc182904061 \h 10 PAGE 11 实验一 直接法解线性方程组 一、实验目的 掌握列选主元消去法与追赶法解线性方程组。 二、实验内容 分别写出Guass列选主元消去法与追赶法的算法，编写程序上机调试出结果，要求所编程序适用于任何一解线性方程组问题，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。 1、用Guass列选主元消去法求解方程组 SKIPIF 1 0 2、用追赶法求解方程组 SKIPIF 1 0 三、实验仪器设备与材料 主流微型计算机 四、实验原理 1、Guass列选主元消去法 对于AX =B 1）、消元过程：将（A|B）进行变换为 SKIPIF 1 0 ，其中 SKIPIF 1 0 是上三角矩阵。即： SKIPIF 1 0 k从1到n-1 列选主元 选取第k列中绝对值最大元素 SKIPIF 1 0 作为主元。 换行 SKIPIF 1 0 归一化 SKIPIF 1 0 消元 SKIPIF 1 0 2）、回代过程：由 SKIPIF 1 0 解出 SKIPIF 1 0 。 SKIPIF 1 0 2、追赶法 线性方程组为： SKIPIF 1 0 做LU分解为： SKIPIF 1 0 分解公式： SKIPIF 1 0 则 SKIPIF 1 0 回代公式： SKIPIF 1 0 SKIPIF 1 0 五、实验步骤 1、理解并掌握列选主元消去法与追赶法； 2、画出列选主元消去法与追赶法的流程图 3、使用VC语言编写出相应的程序并调试验证通过 六、实验报告要求 1、统一使用《武汉科技大学实验报告》本书写，实验报告的内容要求有：实验目的、实验内容、程序流程图、源程序、运行结果及实验小结六个部分。 2、源程序需打印后粘贴在实验报告册内； 3、运行结果以屏幕截图形式保存并打印后粘贴在实验报告册内。 七、实验注意事项 注意如何定义数据结构以保存矩阵和解以降低算法的复杂性。 八、思考题 若使用全主元消去法，在编程中应如何记录保存对于未知数的调换。 实验二 插值方法 一、实验目的 掌握拉格郎日插值法与牛顿插值法构造插值多项式。 二、实验内容 分别写出拉格郎日插值法与牛顿插值法的算法，编写程序上机调试出结果，要求所编程序适用于任何一组插值节点，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。 已知下列函数表 xi 0.56160 0.56280 0.56401 0.56521 yi 0.82741 0.82659 0.82577 0.82495 求x=0.5635时的函数值。 三、实验仪器设备与材料 主流微型计算机 四、实验原理 已知n个插值节点的函数值，则可由拉格郎日插值公式与牛顿插值公式构造出插值多项式，从而由该插值多项式求出所要求点的函数值。拉格郎日插值公式与牛顿插值公式如下： 1、Lagrange插值公式 SKIPIF 1 0 SKIPIF 1 0 2、Newton插值公式 SKIPIF 1 0 五、实验步骤 1、理解并掌握拉格郎日插值法与牛顿插值法的公式； 2、画出拉格郎日插值法与牛顿插值法算法的流程图； 3、使用VC编写出相应的程序并调试验证通过。 六、实验报告要求 1、统一使用《武汉科技大学实验报告》本书写，实验报告的内容要求有：实验目的、实验内容、程序流程图、源