资源描述：
《列主元高斯消去法和列主元三角分解法解线性方程(20210420081514).docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、列主元高斯消去法和列主元三角分解法解线性方程-CAL-FENGHAI.-(YICAI)-CompanyOne1计算方法实验报告1［课题名称］用列主元高斯消去法和列主元三角分解法解线性方程【目的和意义】高斯消去法是一个古老的求解线性方程组的方法，但山它改进得到的选主元的高斯消去法则是LI前计算机上常用的解低阶稠密矩阵方程组的有效方法。用高斯消去法解线性方程组的基本思想时用矩阵行的初等变换将系数矩阵A约化为具有简单形式的矩阵(上三角矩阵、单位矩阵等)，而三角形方程组则可以直接回带求解用高斯消去法解线性方程组Ax=b(其中A
2、WRnXn)的计算量为：乘除法运"(畀―1)(/I—l)n(2n—1)n(n—1)n(n+1)n2理22.u.-txr,“M1)=++F=4-H——加减运算步骤为262233,av_(//—i)n(2/7—1)n(n—1)n(n—1)_(/?—l)n(2.n+5)算步骤为'=2=6。相比之下，传统的克莱姆法则则较为繁琐，如求解20阶线性方程组，克莱姆法则大约要5x10"次乘法，而用高斯消去法只需要3060次乘除法。在高斯消去法运算的过程中，如果出现abs(A(i,i))等于零或过小的惜况，则会导致矩阵元素数量级严重增长
3、和舍入误差的扩散，使得最后的计算结果不可靠，所以LI前计算机上常用的解低阶稠密矩阵方程的快速有效的方法时列主元高斯消去法，从而使计算结果更加精确。2、列主元三角分解法高斯消去法的消去过程，实质上是将A分解为两个三角矩阵的乘积A二LU,并求解Ly=b的过程。回带过程就是求解上三角方程组Ux=yo所以在实际的运算中，矩阵L和U可以直接计算出，而不需要任何中间步骤，从而在计算过程中将高斯消去法的步骤进行了进一步的简略，大大提高了运算速度，这就是三角分解法采用选主元的方式与列主元高斯消去法一样，也是为了避免除数过小，从而保证了
4、计算的精确度【计算公式】1、列主元高斯消去法设有线性方程组Ax=b,其中设A为非奇异矩阵。方程组的增广矩阵为nlnnb第1步（kJ）:首先在A的第一列中选取绝对值最大的元素切，作为第一步的主兀素:=maxan
5、工0然后交换（A,b）的第1行与第I行元素，再进行消元计算。设列主元素消去法已经完成第1步到第k]=第k步计算如下：对于k=l»2,…，n-1（1）按列选主元：即确定t使（2）如果tHk,则交换［A,b］
6、第t行与第k行元素。（3）消元计算%J片加=——,（i=k+1,・••,〃）Clkkaj%+叫,（：j=*+1,…/）bj叫片（j=k+l,・・・M）消元乘数mik满足:（4）回代求解2、列主元三角龙卿•法对方程组的增广矩卩轻=［人如过k-1步分解后，可变成如下形式:M11W12…5"%•••%…％>1W22…"心li2k♦••皿j…"2”■••■•■■■■•■•■■■•■■••91.2…"mUk-.k3•••…"in…♦♦•…%bk■■•■•••■■■■♦■■■■•■■13第k步分解，为了避免用绝对值很小的数5作除
7、数，引进量'女-】ukj一E—%(丿・=£,k+1,…M；*=1,2,…/)m-1<&-1/汝=(血一)/%(，=k+1,£+2,・・•/;&=1,2,•••丿).wr-lSi=7一工.％W=k、k+X・・m)ni-l,于是有％=吐。如果T忠x
8、a」则将矩阵的第t行与第k行元素互换，将(i,j)位置的新元素仍记为'力或5,然后再做第k步分解，这时%=sk(sk即交换前的耳)，lik=6Isk(j=k+l,k+2,・・・,〃)且
9、Z;J<1(i=R+l,R+2,•••,“)，【列主元高斯消去法程序流程图】13输入万程：系
10、数矩P5A，b以恐主元允许最小值ck=l,2,3n-1扌安列迭取主元及其下标[孙切=绅心(小e(A(k:n,k)))并将在ACk:n,k)中的行号转换対在卫中的行号'即q-q+k-1擁行tompl-A(fc,:);A〔k,J=A〔mJ；ACct，:)=templ;tatnp2=b(k/);b(k,b(q,:)=temp2;%列主元法高斯消去法解线性方程Ax=b%判断输入的方程组是否有误【列主元高斯消去法Matlab主程序】functionx=gaussl(A,b,c)if(length(A)"=length(b))di
11、sp('输入方程有误！')return;enddisp(*原方程为AX=b：')%显示方程组AbdispC•)n=length(A);fork=l:n-l%找列主元[p,q]=max(abs(A(k:n,k)));%找出第k列中的最大值，其下标为[p,q]q=q+k-l;%q在A(k:n,k)中的行号转换为在A中的行号13讦abs
