Gauss列主元消去法.docVIP

贵州师范大学数学与计算机科学学院学生实验报告 课程名称： 数值分析 班级： 实验日期： 年 月 日 学 号： 姓名： 指导教师： 实验成绩： 一、实验名称 二、实验目的及要求 三、实验环境四、实验内容 取初始向量,迭代控制条件为 请绘制出迭代次数与松弛因子关系的函数曲线,给出最佳松弛因子.SOR迭代的收敛速度是否一定比Gauss-Seidel迭代快? 2. 编制列主元 Gauss 消去法函数 子程序 ,并用于解 要求输出方程组的解和消元后的增广矩阵. 注：题2必须写实验报告 五、算法描述及实验步骤 . 输入 ,,. 输出 方程组的解或失败信息. 步1 对执行步2→步4 . 步2 调选列主元模块 . 步3 若，则“消去法失败”，结束 . 步4 对执行步5→步6 . 步5 对执行 . 步6 . 步7 . 步8 对执行 . 步9 输出 . 选列主元模块： 功能 选列主元 . 输入 . 输出 . 步1 ； . 步2 对执行若则； . 步3 若，则交换和，；交换和 . 步4 返回主模块 . 六、调试过程及实验结果 A [12,-3,3;-18,3,-1;1,1,1]; b [15;-15;6]; x Gauss1 A,b Ab -18.0000 3.0000 -1.0000 -15.0000 0 1.1667 0.9444 5.1667 0 0 3.1429 9.4286 index 1 x 1.0000 2.0000 3.0000 七、总结在Gauss消去法中有着突出的作用，第k步消元时，要用作除数，如果 0消元会失败，即使主元0，但很小时，舍入误差也会使计算结果面目全非，避免这种缺陷的基本方法就是选主元。 通过选主元，就可避免绝对值小的数作除数，从而避免舍入误差的恶性增长，使得Gauss列主元消去法是解中小规模的线性方程组和某些大型稀疏线性方程组的有效方法。 八、附录（源程序清单）function [x,index] Gauss1 A,b [n,m] size A ;x zeros n,1 ; index 1 for k 1:n-1 a_max 0; for i k:n if abs A i,k a_max a_max abs A i,k ;r i; end end if a_max 1e-10 index 0;return; end if r k for j k:n z A k,j ;A k,j A r,j ;A r,j z; end z b k ;b k b r ;b r z; end for i k+1:n m A i,k /A k,k ; for j k:n A i,j A i,j -m*A k,j ; end b i b i -m*b k ; end end if abs A n,n 0 index 0;return; end Ab [A,b] x n b n /A n,n ; for i n-1:-1:1 for j i+1:n b i b i -A i,j *x j ; end x i b i /A i,i ; end