08125077赵娅丽(实验五).docVIP

实验五 线性方程组直接解法 姓名 赵娅丽 学 班级 信息0811班 一、实验目的 掌握线性方程组直接解法的基本思路和迭代步骤 培养编程与上机调试能力 二、实验内容 ㈠、分别用列主元消元法与不选主元消元法求解，分析对结果的影响。 （1）， （2）， ㈡、对矩阵做LU分解 三、解答过程 A.建立m文件Gauss.m Gauss.m程序编写如下 function[x,det,flag]=Gauss(A,b) %% 线性方程组的列主元Gauss消去法 % A为方程组的系数矩阵 %b为方程组的右端项 %x为方程组的解 %det为系数矩阵的行列式的值 % falg=为指标向量，flag=failure表示计算失败，flag=‘OK’表示计算成功 [n,m]=size(A);nb=length(b); %当方程组行与列的维数不相等时，停止计算，并输出出错信息 if n~=m error(The row and columns of matrix A must be equal!) return; end % 当方程组与右端项的维数不匹配时，停止计算，并输出出错信息 if m~=nb error(The columns of A must be equal the lengh of b!); return; end % 开始计算