計算方法实验报告习题2(浙大版).docVIP

计算方法实验报告 实验名称： 实验2 列主元素消去法解方程组 引言 很小，用作除数会导致中间结果矩阵元素数量级严重增长和舍入误差的扩散，使结果不可靠；采用选主元素的三角分解法可以避免此类问题。高斯消去法的消去过程，实质上是将A分解为两个三角矩阵的乘积A=LU，并求解Ly=b的过程。回带过程就是求解上三角方程组Ux=y。所以在实际的运算中，矩阵L和U可以直接计算出，而不需要任何中间步骤，从而在计算过程中将高斯消去法的步骤进行了进一步的简略，大大提高了运算速度，这就是三角分解法。采用选主元的方式与列主元高斯消去法一样，也是为了避免除数过小，从而保证了计算的精确度。 实验目的和要求算法原理与流程图 原理 经过k-1步分解后，可变成如下形式： 第k步分解，为了避免用绝对值很小的数作除数，引进量 ，于是有=。如果，则将矩阵的第t行与第k行元素互换，将（i,j）位置的新元素仍记为或，然后再做第k步分解，这时 流程图 程序代码及注释 %运用列主元素消去法求解方程组并实现PA=LU function [x,L,U,P,IP]=lzylu(A,b) %定义函数 [m,n]=size(A); %计算系数矩阵A的行列数 if m~=n error(系数矩阵不是方阵); return; end if det(A)==0 %计算矩阵A的行列式，若为零则A是奇异