实验名称：列主元消去法解方程组1 引言我们知道，高斯消去法是一个古老的解线性方程组的方法。
而在用高斯消去法解Ax=b时，其中设A为非奇异矩阵，可能出现的情况，这时必须进行带行交换的高斯消去法。
但在实际计算中即使但其绝对值很小时，用作除数，会导致中间结果矩阵元素数量级严重增长和舍入误差的扩散，使得最后的结果不可靠。
因此，小主元可能导致计算的失败，我们应该避免采用绝对值很小的主元素。
为此，我们在高斯消去法的每一步应该在系数矩阵或消元后的低阶矩阵中选取绝对值最大的元素作为主元素，保持乘数，以便减少计算过程中舍入误差对计算解的影响。
一种方式是完全主元消去法，这种消去法是在每次选主元时，选择为主元素。
这种方法是解低阶稠密矩阵方程组的有效方法，但这种方法在选取主元时要花费一定的计算机时间。
实际计算中我们常采用部分选主元的的消去法。
列主元消去法即在每次选主元时，仅依次按列选取绝对值最大的元素作为主元素，且仅交换两行，再进行消元计算。
2 实验目的和要求运用matlab编写一个.m文件，要求用列主元消去法求解方程组（实现PA=LU）：要求输出以下内容：（1）计算解x；（2） L,U；（3）整形数组IP（i）（i=1,2，…，n-1）（记录主行信息）3 算法原理与流程图（1）算法原理设有线性方程组Ax=b，其中设A为非奇异矩阵。
方程组的增广矩阵为第1步（k=1）：首先在A的第一列中选取绝对值最大的元素，作为第一步的主元素:，然后交换（A，b）的第1行与第i1行元素，再进行消元计算。
设列主元素消去法已经完成第1步到第k-1步的按列选主元，交换两行，消元计算得到与原方程组等价的方程组第k步计算如下：对于k=1，2，…，n-1（1）按列选主元：即确定ik使（2）如果，则A为非奇异矩阵，停止计算。
（3）如果ik≠k，则交换[A，b]第ik行与第k行元素。
（4）消元计算消元乘数满足:（5）回代求解计算解在常数项b(n)内得到。
高斯列主元消去法2.3高斯列主元消去法解线性方程组一：问题的提出我们都知道，高斯列主元素消去法是计算机上常用来求解线性方程组的一种直接的方法。
就是在不考虑舍入误差的情况下，经过有限步的四则运算可以得到线性方程组的准确解的一类方法。
实际运算的时候因为只能有限小数去计算，因此只能得到近似值。
在实际运算的时候，我们很多时候也常用高斯消去法。
但是高斯消去法在计算机中运算的时候常会碰到两个问题。
1．一旦遇到一些主元等于0，消元过程便无法进行下去。
2．在长期使用中还发现，即使消元过程能进行下去，但是当一些主元的绝对值很小时，求解出的结果与真实结果相差甚远。
为了避免高斯消去法消元过程中出现的上述两个问题，一般采用所谓的选择主元法。
其中又可以分为列选主元和全面选主元两种方法。
目前计算机上常用的按列选主元的方法。
因此我在这里做的也是列选主元高斯消去法。
二、算法的基本思想大家知道，如果一个线性方程组的系数矩阵是上三角矩阵时，即这种方程组我们称之为上三角方程组，它是很容易求解的。
我们只要把方程组的最下面的一个方程求解出来，在把求得的解带入倒数第二个方程，求出第二个解，依次往上回代求解。
然而，现实中大多数线性方程组都不是上面所说的上三角方程组，所以我们有可以把不是上三角的方程通过一定的算法化成上三角方程组，由此我们可以很方便地求出方程组的解。
高斯消元法的目的就是把一般线性方程组简化成上三角方程组。
于是高斯消元法的基本思想是：通过逐次消元将所给的线性方程组化为上三角形方程组，继而通过回代过程求解线性方程组。
三、算法的描述1、设有n元线性方程组如下：＝2、第一步：如果a11!=0,令li1= ai1/a11, I= 2,3,……,n用（-li1）乘第一个方程加到第i个方程上,得同解方程组：a(1)11a(1)12...a(1)1nx1b(1)1a(1)21a(1)22...a(1)2nx2b(1)2 .......=.a(1)n-11 a(1)n-12 . .a(1)n-1nxn-1b(1)n-1a(1)n1a(1)n2. . . a(1)nnxnb(1)n简记为：A(2)x=b(2)其中a(2)ij = a(1)ij – li1 a(1)1j ,I ,j = 2,3,..,nb(2)I = b(1)I – li1 b(1)1 ,I = 2,3,...,n第二步：如果a(2)22!=0,令li2= a（2）i2/a（2）22, I= 3,……,n依据同样的原理，对矩阵进行化间（省略），依次下去，直到完成！最后，得到上三角方程组：a(1)11a(1)12...a(1)1nx1b(1)1a(1)22 ...a(1)2nx2b(1)2 .......=.. .a(n-1)n-1nxn-1b(n-1)n-1. . . a(n)nnxnb(n)n简记为：A(n)x=b(n)最后从方程组的最后一个方程进行回代求解为：n = b(n) / a(n)nni = ( b(k)k - a(k)kjxj ) / a(k)kk以上为高斯消去法的基本过程。
目录摘要 (1)一．用列主元消去法解方程组 (2)1．问题的提出 (2)2．问题的分析 (2)3．问题的解决 (3)二．编写一个列主元消去法求逆矩阵的程序 (4)1．问题的提出 (4)2．问题的分析 (4)3．问题的解决 (5)Ax (5)三．用LU分解法解方程组b1．问题的提出 (5)2．问题的分析 (5)3．问题的解决 (6)四．用改进平方根法解方程组 (7)1．问题的提出 (7)2．问题的分析 (7)3．问题的解决 (8)五．用追赶法解方程组 (9)1．问题的提出 (9)2．问题的分析 (9)3．问题的解决 (10)六．分别用雅可比迭代法与高斯-赛德尔迭代法解方程组 (11)1．问题的提出 (11)2．问题的分析 (11)3．问题的解决 (12)参考文献 (14)个人体会 (15)附录：程序代码 (16)摘要在科技研究和工程技术所提出的计算问题中，经常会遇到线性方程组的求解问题，这里主要是有关线性方程组的直接解法。
解线性方程组的直接法是用有限次运算求出线性方程组Ax=b 的解的方法。
线性方程组的直接法主要有Gauss消元法及其变形、LU(如Doolittle、Crout方法等)分解法和一些求解特殊线性方程组的方法（如追赶法、LDLT法等）。
这里主要有列主元消元法,LU分解法,改进的平方根法，追赶法和雅可比迭代，高斯—塞德尔迭代的构造过程及相应的程序。
线性方程的解法在数值计算中占有极重要的地位，因此，线性方程组的求解是数值分析课程中最基本的内容之一。
关键词：列主元消元法；LU分解；改进平方根法；追赶法；雅可比迭代；高斯—塞德尔迭代一．用列主元消去法解方程组：1．问题的提出：（1）⎪⎪⎩⎪⎪⎨⎧=-++--=+--=+-+=++4323231243432143214321421x x x x x x x x x x x x x x x （2）⎪⎪⎩⎪⎪⎨⎧=++--=++-=-+--=-+-434220332282432132143214321x x x x x x x x x x x x x x x2．问题的分析：列主元消去算法主要分为两个过程：消去过程和回代过程1． 消去过程对1,,2,1-=n k(1)选主元 找k i ∈}{,,n k ，⋯使)()(max k ik ni k ki k a ak ≤≤= (2)若0)(=k a k ik 则停止计算（detA=0）(3)若k i k ≠ 则换行()()k i k E E ↔ （4）消元对i =1,,1++n k)()(k kkk ik a a ik l =对1,,1++=n k j )()()1(k kjik k ij k ija l a a -=+ 2.回代过程（1）若0)(=n nn a 则停止计算（detA=0） (2) )()(1,n nnn n n a a n x +=（3）对1,,1 -=n i)(1)()(1,n iini j jn ij n n i a x a a i x ∑=+=+-3．问题的解决：（1）解：对于()b A |)1(=()b A |=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛------43141321************ 第1步选列主元为,3)1(31=a 31=i ，作变换()()31E E ↔，然后计算667.03221==l ， 333.03131==l ，333.03141-==-l再作变换()()(),414143131321212,,E E l E E E l E E E l E →-→-→-得到())2()2(|b A=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛------3533333.0667.2667.10333.2333.0333.10333.0333.0667.102113 第2步，对)2(A 选列主元为667.135)2(22==a ，22=i ，计算8.05432==l ， 142=l ， 再做变换32323)(E E l E →-，42424)(E E l E →-，得到()⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-----=05/1333030015/3915/9003/13/13/502113)3()3(b A消去过程结束，回代计算得到解10214321===-=x x x x所以原方程组的解为TX )1,0,2,1(-=。
线性方程组的消元法线性方程组的消元法是解决线性方程组的常用方法之一，通过逐步消去未知数的系数，将方程组转化为更简单的形式，从而求得方程组的解。
本文将详细介绍线性方程组的消元法及其应用。
1. 消元法简介消元法是一种通过逐步消除未知数的系数，将线性方程组转化为更简单形式的方法。
它的基本思想是通过不断的代入与消去操作，将方程组转化为三角形式或最简形式，从而求得方程组的解。
2. 线性方程组的一般形式线性方程组的一般形式可以表示为：a₁₁x₁ + a₁₂x₂ + ... + a₁ₙxₙ = b₁a₂₁x₁ + a₂₂x₂ + ... + a₂ₙxₙ = b₂...aₙ₁x₁ + aₙ₂x₂ + ... + aₙₙxₙ = bₙ其中，a₁₁、a₁₂、...、aₙₙ为未知数的系数，b₁、b₂、...、bₙ为常数项。
3. 消元法的步骤（1）选取主元：根据方程组的特点，选择一项作为主元，并将其系数置为1，并且使其所在的其他行对应的列的系数皆为0，这样可以简化计算过程并减少误差。
（2）代入消元：选择一个非主元进行代入，将其代入主元所在的其他方程中，从而消去该未知数。
（3）重复步骤（1）和（2），直至将所有的非主元都消去为止。
（4）最后得到一个三角形形式的线性方程组，可以通过回代法求解该方程组的解。
4. 消元法的应用消元法广泛应用于各个领域，特别是在科学和工程领域中具有重要作用。
以下是几个应用实例：（1）经济学中的输入产出模型：通过消元法可以分析不同产业之间的投入产出关系，从而得出经济模型的解释。
（2）物理学中的电路分析：通过消元法可以简化复杂的电路方程组，从而计算出电路中各个节点的电压和电流。
（3）化学反应平衡问题：通过消元法可以解决化学反应平衡过程中的复杂线性方程组，从而得到反应物和生成物的浓度。
5. 总结消元法是一种解决线性方程组的有效方法，通过逐步消除未知数的系数，将方程组转化为更简单的形式，从而求得方程组的解。
