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

计算方法实验报告1【课题名称】用列主元高斯消去法和列主元三角分解法解线性方程【目的和意义】高斯消去法是一个古老的求解线性方程组的方法，但由它改进得到的选主元的高斯消去法则是目前计算机上常用的解低阶稠密矩阵方程组的有效方法。用高斯消去法解线性方程组的基本思想时用矩阵行的初等变换将系数矩阵A约化为具有简单形式的矩阵（上三角矩阵、单位矩阵等），而三角形方程组则可以直接回带求解用高斯消去法解线性方程组Ax=b（其中ARn×n∈）的计算量为：乘除法运算步骤为，加减运算步骤为。相比之下，传统的克莱姆法则则较为繁琐，如求解20阶线性方程组，克莱姆法则大约要次乘法，而用高斯消去法只需要3060次乘除法。在高斯消去法运算的过程中，如果出现abs(A(i,i))等于零或过小的情况，则会导致矩阵元素数量级严重增长和舍入误差的扩散，使得最后的计算结果不可靠，所以目前计算机上常用的解低阶稠密矩阵方程的快速有效的方法时列主元高斯消去法，从而使计算结果更加精确。2、列主元三角分解法高斯消去法的消去过程，实质上是将A分解为两个三角矩阵的乘积A=LU，并求解Ly=b的过程。回带过程就是求解上三角方程组Ux=y。所以在实际的运算中，矩阵L和U可以直接计算出，而不需要任何中间步骤，从而在计算过程中将高斯消去法的步骤进行了进一步的简略，大大提高了运算速度，这就是三角分解法采用选主元的方式与列主元高斯消去法一样，也是为了避免除数过小，从而保证了计算的精确度【计算公式】1、列主元高斯消去法设有线性方程组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）回代求解2、列主元三角分解法对方程组的增广矩阵经过k-1步分解后，可变成如下形式：{xn←bnannxi←(bi−∑j=i+1naijxj)aii,(i=n−1,n−2,⋯,1)第k步分解，为了避免用绝对值很小的数作除数，引进量，于是有=。如果，则将矩阵的第t行与第k行元素互换，将（i，j）位置的新元素仍记为或，然后再做第k步分解，这时【列主元高斯消去法程序流程图】【列主元高斯消去法Matlab主程序】functionx=gauss1(A,b,c)%列主元法高斯消去法解线性方程Ax=bif(length(A)~=length(b))%判断输入的方程组是否有误disp('输入方程有误！')return;enddisp('原方程为AX=b：')%显示方程组Abdisp('------------------------')n=length(A);fork=1:n-1%找列主元[p,q]=max(abs(A(k:n,k)));%找出第k列中的最大值，其下标为[p,q]q=q+k-1;%q在A(k:n,k)中的行号转换为在A中的行号ifabs(p)ktemp1=A(k,:);%列主元所在行不是当前行，将当前行与列主A(k,:)=A(q,:);元所在行交换(包括b)A(q,:)=temp1;temp2=b(k,:);b(k,:)=b(q,:);b(q,:)=temp2;end%消元fori=k+1:nm(i,k)=A(i,k)/A(k,k);%A(k,k)将A(i,k)消为0所乘系数A(i,k:n)=A(i,k:n)-m(i,k)*A(k,k:n);%第i行消元处理b(i)=b(i)-m(i,k)*b(k);%b消元处理endenddisp('消元后所得到的上三角阵是')A%显示消元后的系数矩阵b(n)=b(n)/A(n,n);%回代求解fori=n-1:-1:1b(i)=(b(i)-sum(A(i,i+1:n)*b(i+1:n)))/A(i,i);endclearx;disp('AX=b的解x是')x=b;【调用函数解题】【列主元三角分解法程序流程图】【列主元三角分解法Matlab主程序】①自己编的程序：functionx=PLU(A,b,eps)%定义函数列主元三角分解法函数if(length(A)~=length(b))%判断输入的方程组是否有误disp('输入方程有误！')return;enddisp('原方程为AX=b：')%显示方程组Abdisp('------------------------')n=length(A);A=[Ab];%将A与b合并，得到增广矩阵forr=1:nifr==1fori=1:n[cd]=max(abs(A(:,1)));%选取最大列向量，并做行交换ifc<=eps%最大值小于e，主元太小，程序结束break;elseendd=d+1-1;p=A(1,:...