数值实验期实习题目2011.docVIP

数值实验2011暑期实习题目 ********************************************************************实验一：迭代法求方程根的收敛性 实验目的：体会在非线性方程求根的迭代法中，迭代函数和初始迭代值的选取对迭代收敛性的影响。 实验内容：考虑一个简单的代数方程 SKIPIF 1 0 ，针对上述方程，可以构造多种迭代法，如下列几种迭代格式：① SKIPIF 1 0 ；② SKIPIF 1 0 ；③ SKIPIF 1 0 ；④ SKIPIF 1 0 。 实验要求： （1）取定某个初始值 SKIPIF 1 0 ，按如上四种迭代格式进行计算，它们的收敛性如何？重复选取不同的初始值，反复实验，并自行设计一种比较形象的记录方式（如利用Matlab的图形功能等），分析四种迭代法的收敛性与初值选取的关系。 （2）选取第④种迭代格式（Newton迭代法），取不同的初始值进行迭代，结果如何？并分析迭代法对不同的初值是否有差异。 （3）对上述四种迭代格式，编制Steffensen迭代程序，选取不同的初始值，输出迭代次数和方程的根，并与（1）、（2）中的结果进行比较。 ********************************************************************实验二：多项式插值的振荡现象，即插值的龙格现象 实验目的：体会在用多项式插值方法逼近一个函数时，使用的节点越多，插值多项式的次数就越高，随着插值多项式次数的增加，插值函数是否也更加逼近给定的函数。 实验内容： 考虑龙格函数 SKIPIF 1 0 在区间 SKIPIF 1 0 的一个等距划分，分点取为 SKIPIF 1 0 则拉格朗日插值多项式为 SKIPIF 1 0 其中的 SKIPIF 1 0 是n次拉格朗日插值基函数。 实验要求： （l）选择不断增大的分点数目 SKIPIF 1 0 时，画出原函数 SKIPIF 1 0 及插值多项式函数 SKIPIF 1 0 在 SKIPIF 1 0 上的图像，比较并分析实验结果。 （2）选择其他的函数，例如定义在区间[-5，5]上的函数 SKIPIF 1 0 重复上述（1）的过程，比较并分析实验结果。 （3）区间 SKIPIF 1 0 上切比雪夫点的定义为 SKIPIF 1 0 以 SKIPIF 1 0 为插值节点，构造上述各函数的拉格朗日插值多项式，比较并分析实验结果。 （4）利用Newton插值方法重复上述实验过程，比较并分析实验结果。 ******************************************************************** 实验三： 用Gauss消去法求解线性方程组 实验目的：体会在Gauss消去法中,主元的选择对数值算法的影响以及方程组的系数矩阵的条件数对解的影响。 实验内容：考虑线性方程组 SKIPIF 1 0 编制求解线性方程组的Gauss列主元消去法。 实验要求： 取三对角矩阵 SKIPIF 1 0 ， SKIPIF 1 0 ，则方程组的解为 SKIPIF 1 0 。取 SKIPIF 1 0 ，计算矩阵的条件数，并在消去过程分别按自然顺序、按模最小的元素和按模最大的元素作为主元，比较并分析实验结果。 注：若方程组阶数比较大时，输出结果要求用文件存储。 ******************************************************************** 实验四：三次样条插值方法的实现 实验目的：体会三次样条插值函数的构造方法及其应用，并比较利用该方法得到的分段多项式函数和被逼近函数的近似程度（用图形的形式比较）。 实验内容： 考虑函数 SKIPIF 1 0 在区间 SKIPIF 1 0 上，分成 10等分。利用三次样条插值的三弯矩方法（ SKIPIF 1 0 方法），编制用于第一种和第二种边界条件的程序。设已知数据为第一种边界条件： SKIPIF 1 0 第二种边界条件： SKIPIF 1 0 实验要求： 分别用所编程序求解，输出各插值节点的弯矩值 SKIPIF 1 0 和插值中点的样条函数值，并