线性方程组直接解法课程.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 % 开始计算，并赋初值 flag= OK;det=1;x=zeros(n,1); for k=1:n-1 %选主元 max1=0; for i=k:n if abs(A(i,k))max1 max1=abs(A(i,k));r=i; end end if max1le-10 flag=failure;return; end %交换两行 if rk for j=k:n z=A(k,j);A(r,j)=A(r,j);A(r,j)=z; end z=b(k);b(k)=b(r);b(r)=z;det=-det; end %消元过程 for i=k+1:n m=A(i,k)/A(k,k); for j=k+1:n A(i,j)=A(i,j)-m*A(k,j); end b(i)=b(i)-m*b(k); end det=det*A(k,k); end det=det*A(n,n) %回代过程 if abs(A(n,n))le-10 flag=failure; return; end for k=n:-1;1 for j=k+1:n b(k)=b(k)-A(k,j)*x(j); end x(k)=b(k)/A(k,k); end x(k)=b(k)/A(k,k); end B.建立m文件Gauss_1.m Gauss_1.m程序编写如下 function[x,det,flag]=Gauss_1(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 % 开始计算，并赋初值 flag= OK;det=1;x=zeros(n,1); for k=1:n-1 %消元过程 for i=k+1:n m=A(i,k)/A(k,k); for j=k+1:n A(i,j)=A(i,j)-m*A(k,j); end b(i)=b(i)-m*b(k); end det=det*A(k,k);