数值线性代数第二版徐树方高立张平文上机习题第一章实验报告.pdfVIP

上机习题 1.先用你所熟悉的的计算机语言将不选主元和列主元 Gauss 消去法编写成通用的子程序； 然 后用你编写的程序求解 84 阶方程组；最后将你的计算结果与方程的精确解进行比较，并就 此谈谈你对 Gauss 消去法的看法。 Sol： （1）先用 matlab 将不选主元和列主元 Gauss 消去法编写成通用的子程序，得到 L ,U , P ： 不选主元 Gauss 消去法： L,U GaussLA( A) 得到 L ,U 满足 A LU 列主元 Gauss 消去法： L ,U , P GaussCol(A) 得到 L ,U , P 满足 PA LU （2 ）用前代法解 Ly b or Pb ，得 y 用回代法解 Ux y ，得 x 求解程序为 x Gauss A,b, L,U , P （P 可缺省，缺省时默认为单位矩阵） （3 ）计算脚本为