实验指导一.docVIP

PAGE PAGE 1实验指导一实验一 牛顿插值法1.1 实验目的掌握牛顿插值法的基本思路和步骤；培养编程与上机调试能力。1.2 算法描述1.2.1 牛顿插值法基本思路 给定插值点序列（。构造牛顿插值多项式。输入要计算的函数点并计算的值，利用牛顿插值公式，当增加一个节点时，只需在后面多计算一项，而前面的计算仍有用；另一方面的各项系数恰好又是各阶差商，而各阶差商可用差商公式来计算。1.2.2 牛顿插值法计算步骤输入值及（；要计算的函数点。对给定的由 计算的值。3．输出。1.3 实验内容给定构造牛顿插值函数计算。实验二 最小二乘法2.1 实验目的掌握最小二乘法的基本思路和拟合步骤；培养编程与上机调试能力。2.2 算法描述1.2.1最小二乘法基本思路已知数据对，求多项式使得为最小，这就是一个最小二乘问题。最小二乘法计算步骤 用线性函数为例，拟合给定数据。算法描述：步骤1：输入值，及。步骤2：建立法方程组。步骤3：解法方程组。步骤4：输出。2.3 实验内容已知一组数据如下，求它的线性拟合曲线。1234544.5688.521311实验四 高斯消去法4.1 实验目的掌握高斯消去法的基本思路和迭代步骤；培养编程与上机调试能力。4.2 算法描述4.2.1 高斯消去法基本思路。设有方程组，设是可逆矩阵。高斯消去法的基本思想就是僵局真的初等行变换作用于方程组的增广矩阵，将其中的变换成一个上三角矩阵，然后求解这个三角形方程组。列主元高斯消去法计算步骤 将方程组用增广矩阵表示。步骤1：消元过程，对选主元，找使得如果，则矩阵奇异，程序结束；否则执行（3）。如果，则交换第行与第行对应元素位置，，。消元，对，计算对，计算步骤 2：回代过程：若则矩阵奇异，程序结束；否则执行（2）。对，计算实验内容解方程组实验四 解线性方程组的迭代法5.1实验目的 = 1 \* GB3 ① 掌握解线性方程组的雅可比迭代和高斯-塞德尔迭代算法； = 2 \* GB3 ② 培养编程与上机调试能力.5.2算法步骤5.2.1迭代法的基本思想根据方程组设计出一个迭代公式，然后将任意选取的一初始向量代入迭代公式，求出，再以代入同一迭代公式，求出，如此反复进行，得到向量序列.当收敛时，其极限即为原方程组的解.5.2.2雅可比（Jacobi）迭代法解方程组设方程组的系数矩阵对角线元素，为最大迭代次数，为容许误差. 雅可比（Jacobi）迭代法解方程组算法步骤如下： = 1 \* GB3 ① 取初始向量，令. = 2 \* GB3 ② 对，计算 . = 3 \* GB3 ③ 如果，则输出，结束；否则执行 = 4 \* GB3 ④ = 4 \* GB3 ④ 如果，则不收敛，终止程序；否则，转 = 2 \* GB3 ②5.2.3高斯-塞德尔(Gauss-Seidel)迭代法在雅可比（Jacobi）迭代法中，如果当新的分量求出后，马上用它来代替旧的分量，则可能会更快地接近方程组的准确解.基于这种设想构造的迭代公式称为高斯-塞德尔(Gauss-Seidel)迭代法. 算法可相应地从雅可比（Jacobi）迭代法改造得到.5.3实验题目及参考结果：题目 应用雅可比迭代和高斯-塞德尔迭代算法解线性方程组参考结果： 雅可比迭代法迭代次数20次，结果如下：5.4实验要求： （1）选择一种计算机语言设计出雅可比（Jacobi）迭代法的算法程序，并且选择不同的迭代次数，观察输出结果； （2）利用Matlab求方程组的解步骤如下： 调用格式： %得到线性方程组的解向量Matlab6.1环境下操作如下：; %输入系数矩阵; %输入常数项 %方程组求解5.5思考： = 1 \* GB3 ① 判别迭代法收敛的充分必要条件及充分条件是什么？ = 2 \* GB3 ② 雅可比（Jacobi）迭代法和高斯-塞德尔(Gauss-Seidel)迭代法收敛性的各种判别条件是什么？实验五 非线性方程求根6.1实验目的： = 1 \* GB3 ① 掌握二分法、牛顿迭代法等常用的非线性方程迭代算法； = 2 \* GB3 ② 培养编程与上机调试能力.6.2算法步骤：6.2.1二分法计算的二分法如下： = 1 \* GB3 ① 输入求根取间和误差控制量，定义函数. 如果 ，转 = 2 \* GB3 ②；否则退出选用其它求根方法 = 2 \* GB3 ② 当时，计算中点以及的值； 分情况处理 ：停止计算，，转 = 4 \* GB3 ④ ：修正区间 ：修正区间 = 3 \*