列主元高斯消去法和列主元三角分解法解线性方程

计算方法实验报告1【课题名称】用列主元高斯消去法和列主元三角分解法解线性方程【目的和意义】高斯消去法是一个古老的求解线性方程组的方法，但由它改进得到的选主元的高斯消去法则是目前计算机上常用的解低阶稠密矩阵方程组的有效方法。用高斯消去法解线性方程组的基本思想时用矩阵行的初等变换将系数矩阵A约化为具有简单形式的矩阵（上三角矩阵、单位矩阵等），而三角形方程组则可以直接回带求解用高斯消去法解线性方程组（其中A∈Rn×n）的计算量为：乘除法运算步骤为，加减运算步骤为。相比之下，传统的克莱姆法则则较为繁琐，如求解20阶线性方程组，克莱姆法则大约要次乘法，而用高斯消去法只需要3060次乘除法。在高斯消去法运算的过程中，如果出现abs(A(i,i))等于零或过小的情况，则会导致矩阵元素数量级严重增长和舍入误差的扩散，使得最后的计算结果不可靠，所以目前计算机上常用的解低阶稠密矩阵方程的快速有效的方法时列主元高斯消去法，从而使计算结果更加精确。2、列主元三角分解法高斯消去法的消去过程，实质上是将A分解为两个三角矩阵的乘积A=LU，并求解Ly=b的过程。回带过程就是求解上三角方程组Ux=y。所以在实际的运算中，矩阵L和U可以直接计算出，而不需要任何中间步骤，从而在计算过程中将高斯消去法的步骤进行了进一步的简略，大大提高了运算速度，这就是三角分解法采用选主元的方式与列主元高斯消去法一样，也是为了避免除数过小，从而保证了计算的精确度【计算公式】列主元高斯消去法设有线性方程组Ax=b，其中设A为非奇异矩阵。方程组的增广矩阵为第1步（k=1）：首先在A的第一列中选取绝对值最大的元素，作为第一步的主元素:然后交换（A，b）的第1行与第l行元素，再进行消元计算。设列主元素消去法已经完成第1步到第k-1步的按列选主元，交换两行，消元计算得到与原方程组等价的方程组A(k)x=b(k)第k步计算如下：对于k=1，2，…，n-1（1）按列选主元：即确定t使（2）如果t≠k，则交换[A，b]第t行与第k行元素。（3）消元计算消元乘数mik满足:（4）回代求解列主元三角分解法对方程组的增广矩阵经过k-1步分解后，可变成如下形式：clearx;disp('AX=b的解x是')x=b;【调用函数解题】【列主元三角分解法Matlab主程序】clear;augm=input('Pleaseinputtheaugumentalmatric:');[n,t]=size(augm);fori=1:n%d控制列三角主元的层数b=augm(i,:);p=i;%4到10行是选列主元并交换forw=i+1:nifb(1,1)<=augm(w,i)b=augm(w,:);p=w;endendc=augm(i,:);augm(i,:)=b;augm(p,:)=c;ifp~=i%只有发生了换行才将这种效果输出augm,endforj=i:t%首先变换与augm（i,i）同一行的元素，其列指标从i到ts=0;fork=1:i-1s=s+augm(i,k)*augm(k,j);endaugm(i,j)=(augm(i,j)-s);endforI=i+1:n%再变换与augm（i,i）同一列的元素，其行指标I从i+1到n,列数为is=0;fork=1:i-1%下三角部分列数与第i层一致s=s+augm(I,k)*augm(k,i);endaugm(I,i)=(augm(I,i)-s)/augm(i,i);endaugm,endx(n)=augm(n,t)/augm(n,n);%回代fori=n-1:-1:1s=0;forj=n:-1:i+1s=s+x(j)*augm(i,j);endx(i)=(augm(i,t)-s)/augm(i,i),end【调用函数解题】输入[Ab]可以清晰地看到augm（A）的变化过程：结果输出：【列主元三角分解法程序流程图】【编程疑难】这是第一次用matlab编程，对matlab的语句还不是非常熟悉，因此在编程过程中，出现了许多错误提示。并且此次编程的两种方法对矩阵的运算也比较复杂。问题主要集中在循环控制中，循环次数多了一次或者缺少了