武汉科技大学计算机科学与技术学院数值计算基础实验指导书
数值计算基础
实验指导书
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、使用VC语言编写出相应的程序并调试验证通过
六、实验报告要求
1、统一使用《武汉科技大学实验报告》本书写,实验报告的内容要求有:实验目的、实验内容实验实验二 插值方法
一、实验目的
掌握拉格郎日插值法与牛顿插值法构造插值多项式。
二、实验内容
分别写出拉格郎日插值法与牛顿插值法的算法,编写程序上机调试出结果,要求所编程序适用于任何一组插值节点,即能解决这一类问题,而不是某一个问题。实验中以下列数据验证程序的正确性。
已知下列函数表
xi0.561600.562800.564010.56521yi 0.827410.826590.825770.82495求x=0.5635时的函数值。
三、实验仪器设备与材料
主流微型计算机
四、实验原理
已知n个插值节点的函数值,则可由拉格郎日插值公式与牛顿插值公式构造出插值多项式,从而由该插值多项式求出所要求点的函数值。拉格郎日插值公式与牛顿插值公式如下:
1、Lagrange插值公式
2、Newton插值公式
五、实验步骤
1、理解并掌握拉格郎日插值法与牛顿插值法的公式;
2、画出拉格郎日插值法与牛顿插值法算法的流程图;
3、使用VC编写出相应的程序并调试验证通过。
六、实验报告要求
1、统一使用《武汉科技大学实验报告》本书写,实验报告的内容要求有:实验目的、实验内容实验实验三 数值积分
一、实验目的
掌握复化梯形法与龙贝格法计算定积分。
二、实验内容
分别写出变步长梯形法与Romberge法计算定积分的算法,编写程序上机调试出结果,要求所编程序适用于任何类型的定积分,即能解决这一类问题,而不是某一个问题。实验中以下列数据验证程序的正确性。
求。
三、实验仪器设备与材料
主流微型计算机
四、实验原理
通过变步长梯形法与龙贝格法,我们只要知道已知n个求积节点的函数值,则可由相应的公式求出该函数的积分值,从而不需要求该函数的原函数。变步长梯形法与龙贝格法公式如下:
1、变步长梯形法
用来控制精度
2、龙贝格法
梯形法则二阶公式四阶公式六阶公式八阶公式
用来控制精度
五、实验步骤
1、理解并掌握变步长梯形法与龙贝格法的公式;
2、画出变步长梯形法与龙贝格法的流程图
3、使用VC语言编写出相应的程序并调试验证通过
六、实验报告要求
1、统一使用《武汉科技大学实验报告》本书写,实验报告的内容要求有:实验目的、实验内容实验积分中,被积函数在x=0点函数值为1,对该点在程序设计中应注意对其的定义。
八、思考题
使用复化梯形法与复化Simpson法来计算该问题有何缺点?
实验四 常微分方程的数值解
一、实验目的
掌握改进欧拉法与四阶龙格-库塔求解一阶常微分方程的初值问题。
二、实验内容
分别写出改进欧拉法与四阶龙格-库塔求解的算法,编写程序上机调试出结果,要求所编程序适用于任何一阶常微分方程的数值解问题,即能解决这一类问题,而不是某一个问题。实验中以下列数据验证程序的正确性。
求步长h=0.25。
三、实验仪器设备与材料
主流微型计算机
四、实验原理
常微分方程的数值解主要采用“步进式”,即求解过程顺着节点排列次序一步一步向前推进,在单步法中改进欧拉法和四阶龙格-库塔法公式如下:
1、改进欧拉法
2、四阶龙格-库塔法
五、实验步骤
1、理解并掌握改进欧拉法与四阶龙格-库塔法的公式;
2、画出改进欧拉法与四阶龙格-库塔法的流