2.4直接三角分解法
线性方程组的几种求解方法
甘肃政法学院本科学年论文（设计）题目浅议线性方程组的几种求解方法学号：姓名：指导教师：成绩：__________________完成时间： 2012 年 11 月目录第一章引言 (1)第二章线性方程组的几种解法 (1)2.1 斯消元法 (1)2.1.1 消元过程 (1)2.1.2 回代过程 (2)2.1.3 解的判断 (2)2.2 克莱姆法则 (3)2.3 LU分解法 (4)2.4 追赶法 (6)第三章结束语 (8)致谢 (8)参考文献 (9)摘要:线性方程组是线性代数的核心内容之一，其解法研究是代数学中经典且重要的研究课题.下面将综述几种不同类型的线性方程组的解法，如消元法、克莱姆法则、直接三角形法、、追赶法，并以具体例子介绍不同解法的应用技巧. 在这些解法中，高斯消元法方法，具有表达式清晰，使用范围广的特点.另外，这些方法有利于快速有效地解决线性方程组的求解问题，为解线性方程组提供一个简易平台，促进了理论与实际的结合。
关键词:线性方程组;解法;应用Several methods of solving linear equation groupAbstract: The system of linear equations is one of linear algebra core contents, its solution research is in the algebra the classics also the important research topic. This article summarized several kind of different type system of linear equations solution, like the elimination, the Cramer principle, the generalized inverse matrix law, the direct triangle law, the square root method, pursue the law, and by concrete example introduction different solution application skill. In these solutions, the generalized inverse matrix method, has the expression to be clear, use scope broad characteristic. Moreover, these methods favor effectively solve the system of linear equations solution problem fast, provides a simple platform for the solution system of linear equations, promoted the theory and the actual union.Key word: Linear equations; Solution ; Example第一章 引言线性方程组理论是高等数学中十分重要的内容，而线性方程组的解法是利用线性方程组理论解决问题的关键.下面将介绍线性方程组的消元法、追赶法、直接三角形法等求解方法，为求解线性方程组提供一个平台。
2010版数值分析教学大纲(最新)
《数值分析》课程教学大纲课程代码：090141031课程英文名称：Numerical Analysis课程总学时：64 讲课：64 实验：0 上机：0适用专业：信息与计算科学专业大纲编写（修订）时间：2010.07一、大纲使用说明（一）课程的地位及教学目标《数值分析》是为信息与计算科学专业学生开设的必修课。
在实验方法和理论方法之后，科学计算已成为科学研究的第三种方法。
学习和掌握计算机上常用的数值计算方法已成为现代科学教育的重要内容。
通过本课程的学习，使学生了解和掌握这门课程所涉及的各种常用的数值计算公式、数值方法的构造原理及适用范围，为今后用计算机去有效地解决实际问题打下基础。
通过本课程的学习，学生将达到以下要求：1.掌握数值计算的基本理论和基本方法，提高数学素养；2.具有运用Matlab等工具进行具有一定难度和复杂度的数值解运算的技能，提高应用计算机进行科学与工程计算的能力；3.树立正确的算法设计理念；4.了解数值计算方法的新发展。
（二）知识、能力及技能方面的基本要求1.知识方面的基本要求：掌握算法的基本原理和思想，包括算法的构造、算法处理的技巧、误差分析、收敛性和稳定性等基本理论。
2．基本理论和方法：误差与有效数字定义、函数插值与逼近的方法、积分与微分的数值计算方法、线性方程组的直接解法、线性方程组的迭代法、非线性方程根的求解方法、常微分方程初值问题的数值解法等3.基本能力：使用各种数值方法解决实际计算问题。
不仅要学会“怎样算”，而且必须做到“真会算”，即不仅要知道问题的解是存在的，还必须能求出具体的结果。
具有应用计算机进行科学与工程计算和解决实际问题的能力。
（三）实施说明1．教学方法：课堂讲授中要重点对算法的构造、算法处理技巧和误差分析的讲解；采用启发式教学，培养学生思考问题、分析问题和解决问题的能力；引导和鼓励学生通过实践和自学获取知识，培养学生的自学能力；增加讨论课，调动学生学习的主观能动性；讲课要联系实际并注重培养学生的创新能力。
三角分解法
ai j =
min( i , j ) k =1
∑l
ik
uk j
ai j =
min( i , j )
一般采用列主元 对换, 将 i ,j 对换,对 j = i, i+1, …, n 有 一般采用列主元 ii a ji = ∑ l jk uki + l ji u k= k =1 法增强稳定性. 法增强稳定性.但注意 v i 1 b 也必须做相应的 l ji = ( a ji ∑ l jk uki ) / uii b 行交换. 行交换. k =1
=I
Upper-triangular
Lower-triangular With diagonal entries 1
注: L 为一般下三角阵而 U 为单位上三角阵的分解称为 单位上三角阵的分解称为 Crout 分解. 分解. ~~ 分解, 实际上只要考虑 A* 的 LU 分解,即A* = L U ,则 ~ ~ A= U * L* 即是 A 的 Crout 分解. 分解. =
(
)
n1 Step 6 Set l = 运算量为2 O(n3/6), 比普通 ann ∑ k =1 lnk ; , 比普通LU nn Step 7 Output ( lij for j = 1, …, i and i = 1, …, n );A = LDLT 分解少一半, 次开方. 分解少一半,但有 n 次开方.用
mn1
v A b
( 2) (2)
(1 (1 ( a11) a12) ... a11) n Step n 1: (2 ( v a22) ... a22) n Ln1Ln2 ... L1 A b = ... . . . (n ann)
电力系统分析课程设计报告_4
电力系统分析课程设计报告题目: 电力系统三相对称短路计算专业: 电气工程及其自动化班级:姓名:学号:指导教师:目录电力系统分析........................................................................................................................... - 0 -第一章设计目的与任务 ......................................................................................................... - 2 -1.1设计目的.................................................................................................................... - 2 -1.2设计任务.................................................................................................................... - 2 -第二章基础理论与原理 ......................................................................................................... - 2 -2.1 对称短路计算的基本方法 ....................................................................................... - 2 -2.2 用节点阻抗矩阵的计算方法 ................................................................................... - 4 -2.3 用节点导纳矩阵的计算方法 ................................................................................... - 6 -2.4 用三角分解法求解节点阻抗矩阵 ........................................................................... - 7 -2.5 短路发生在线路上任意处的计算方法 ................................................................... - 8 -第三章程序设计..................................................................................................................... - 9 -3.1 变量说明................................................................................................................... - 9 -3.2 程序流程图............................................................................................................. - 10 -3.2.1主程序流程图 .............................................................................................. - 11 -3.2.2导纳矩阵流程图 .......................................................................................... - 12 -3.2.3三角分解法流程图 ...................................................................................... - 13 -3.3 程序源代码见附录1 ............................................................................................ - 14 -第四章结果分析................................................................................................................... - 14 -第五章收获与建议............................................................................................................... - 15 -参考文献................................................................................................................................. - 17 -附录......................................................................................................................................... - 17 -附录1: 程序源代码..................................................................................................... - 18 - 附录2: 测试系统数据与系统图 ................................................................................... - 23 - 附录3: 测试系统的运行结果- 25 -第一章设计目的与任务1.1设计目的1、加深理解并巩固电力系统发生短路的基本知识。
数值分析-第二章-学习小结
数值分析-第二章-学习小结(总9页) --本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--第2章线性方程组的解法--------学习小结一、本章学习体会本章主要学习的是线性方程组的解法。
而我们则主要学习了高斯消去法、直接三角分解法以及迭代法三种方法。
这三种方法的优缺点以及适用范围各有不同。
高斯消去法中，我们又学习了顺序高斯消去法以及列主元素高斯消去法。
顺序高斯消去法可以得到方程组的精确解，但要求系数矩阵的主对角线元素不为零，而且该方法的数值稳定性没有保证。
但列主元素高斯消去法因为方程顺序的调整，其有较好的数值稳定性。
直接三角分解法中，我们主要学习了Doolitte分解法与Crout分解法。
其思想主要是：令系数矩阵A=UL，其中L为下三角矩阵，U是上三角矩阵，为求AX=b 的解，则引进Ly=b,Ux=y两个方程，以求X得解向量。
这种方法计算量较小，但是条件苛刻，且不具有数值稳定性。
迭代法（逐次逼近法）是从一个初始向量出发，按照一定的计算格式，构造一个向量的无穷序列，其极限才是所求问题的精确解，只经过有限次运算得不到精确解。
该方法要求迭代收敛，而且只经过有限次迭代，减少了运算次数，但是该方法无法得到方程组的精确解。
二、本章知识梳理针对解线性方程组，求解线性方程组的方法可分为两大类：直接法和迭代法，直接法（精确法）：指在没有舍入误差的情况下经过有限次运算就能得到精确解。
迭代法（逐次逼近法）：从一个初始向量出发，按照一定的计算格式，构造一个向量的无穷序列，其极限才是所求问题的精确解，只经过有限次运算得不到精确解。
我们以前用的是克莱姆法则，对于计算机来说，这种方法运算量比较大，因此我们学习了几种减少运算次数的方法，有高斯消去法、直接三角分解法，同时针对病态方程组，也提出了几种不同的解法。
Gauss消去法Gauss消去法由消元和回代两个过程组成，消元过程是指针对方程组的增广矩阵，做有限次初等行变化，使它系数矩阵变为上三角矩阵。
用矩阵的直接三角分解法解方程组
用矩阵的直接三角分解法解方程组矩阵的直接三角分解法（LU分解法）是解线性方程组的一种常用方法。
该方法通过将系数矩阵分解为一个下三角矩阵和一个上三角矩阵的乘积，从而简化方程组的求解过程。
下面我们就来详细分步骤地介绍一下这种方法的求解过程。
第一步，将原线性方程组表示为矩阵形式，即将系数矩阵、未知量矩阵和常数矩阵分别表示为A、X和B。
我们的目标是找到一个下三角矩阵L和一个上三角矩阵U，使得方程组可以表示为LUx = B的形式。
第二步，通过高斯消元法将系数矩阵A化为上三角矩阵U。
具体地，我们将系数矩阵A变换为U的过程可以分解为一系列的初等矩阵变换，例如交换两行、乘以一个非零常数和将某一行加上另一行的若干倍等等。
这些初等矩阵变换可以表示为一个矩阵M的乘积，即A =M1M2...MnU。
从而，我们得到了上三角矩阵U。
第三步，同样通过一系列初等矩阵变换将U转化为下三角矩阵L。
这些初等矩阵变换可以表示为一个矩阵N的乘积，即U = NL1L2...Lm。
从而，我们得到了下三角矩阵L。
第四步，将方程组表示为LUx = B的形式。
具体地，我们将A, X 和B分解为L, U和x的乘积，即A = LU，X = UL，B = Ux。
从而，原方程组可以表示为LUx = B，即L(Ux) = B。
第五步，解方程组L(Ux) = B。
由于L是下三角矩阵，因此可以通过前代法求解得到Ux。
具体地，我们先通过Lw = B求解出向量w，然后再通过Ux = w求解出未知量向量x。
总的来说，矩阵的直接三角分解法（LU分解法）是一种常用的解线性方程组的方法。
它将原方程组表示为LUx = B的形式，然后通过前代法和回代法求解得到未知量向量x。
这种方法具有求解速度快、计算量小的优点，因此在实际应用中得到了广泛的应用。
(整理)线性方程组的直接法
第二章线性方程组的直接法在近代数学数值计算和工程应用中，求解线性方程组是重要的课题。
例如，样条插值中形成的关系式，曲线拟合形成的法方程等，都落实到解一个元线性方程组，尤其是大型方程组的求解，即求线性方程组（2.1）的未知量的数值。
（2.1）其中ai j，bi为常数。
上式可写成矩阵形式Ax = b，即（2.2）其中，为系数矩阵，为解向量，为常数向量。
当detA=D0时，由线性代数中的克莱姆法则，方程组的解存在且惟一，且有为系数矩阵的第列元素以代替的矩阵的行列式的值。
克莱姆法则在建立线性方程组解的理论基础中功不可没，但是在实际计算中，我们难以承受它的计算量。
例如，解一个100阶的线性方程组，乘除法次数约为（101·100！·99），即使以每秒的运算速度，也需要近年的时间。
在石油勘探、天气预报等问题中常常出现成百上千阶的方程组，也就产生了各种形式方程组数值解法的需求。
研究大型方程组的解是目前计算数学中的一个重要方向和课题。
解方程组的方法可归纳为直接解法和迭代解法。
从理论上来说，直接法经过有限次四则运算，假定每一步运算过程中没有舍入误差，那么，最后得到方程组的解就是精确解。
但是，这只是理想化的假定，在计算过程中，完全杜绝舍入误差是不可能的，只能控制和约束由有限位算术运算带来的舍入误差的增长和危害，这样直接法得到的解也不一定是绝对精确的。
迭代法是将方程组的解看作某种极限过程的向量极限的值，像第2章中非线性方程求解一样，计算极限过程是用迭代过程完成的，只不过将迭代式中单变量换成向量而已。
在用迭代算法时，我们不可能将极限过程算到底，只能将迭代进行有限多次，得到满足一定精度要求的方程组的近似解。
在数值计算历史上，直接解法和迭代解法交替生辉。
一种解法的兴旺与计算机的硬件环境和问题规模是密切相关的。
一般说来，对同等规模的线性方程组，直接法对计算机的要求高于迭代法。
对于中等规模的线性方程组，由于直接法的准确性和可靠性高，一般都用直接法求解。
矩阵理论及方法(谢冬秀,雷纪刚,陈桂芝编著)PPT模板
第7章矩阵的广义逆与直积及其应用
7.1矩阵的几种广义逆
01
7.1.1广义逆矩阵的基本概 念
03
7.1.3自反减号逆A<sup></sup><sub>r</sub>
05
7.1.5最小二乘广义逆A<sup></sup><sub>l</sub>
02
7.1.2减号逆
04
7.1.4极小范数广义逆A<sup></sup><sub>m</sub>
01
习题6
06
6.2随机矩 阵与双随 02 机 矩 阵
6 . 5 T o e p l 05 itz矩阵与 Hankel
矩阵
04
6.4广义对 角占优矩阵
6.3M矩
03
阵与 Stieltje
s矩阵
第6章几类特殊矩阵
6.1非负矩阵
6.1.2非负矩 阵谱半径的 界
6.1.1Perron -Frobenius 定理
2.4.3常用的 直接三角分 解法
第2章矩阵的变换与分解
2.5QR分解
2.5.1QR分解的概念
2.5.2QR分解的实际求 法
2.5.3基于QR分解的参 数估计问题
2.5.4矩阵与Hessenberg矩 阵的正交相似问题
04 第3章矩阵范数及其应用
第3章矩阵范数及其应用
3.1向量范数
3.2矩阵范数
06
7.1.6加号逆 A<sup>+</sup>
第7章矩阵的广义逆与直积及其应用
7.2广义逆与线性方程组的解
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
§4 直接三角分解法 一、教学设计1．教学内容：Doolittle 分解法、Crout 分解法，紧凑格式的Doolittle 分解法、部分选主元的Doolittle 分解法。
2．重点难点：紧凑格式的Doolittle 分解法、部分选主元的Doolittle 分解法。
3．教学目标：了解直接三角分解法的基本思想，掌握基本三角分解法及其各种变形。
4．教学方法：讲授与讨论。
二、教学过程在上节中我们用矩阵初等变换来分析Gauss 消去法，得到了重要的矩阵LU 分解定理（定理 3.1，3.2）。
由此我们将得到Gauss 消去法的变形：直接三角分解法。
直接三角分解法的基本想法是，一旦实现了矩阵A 的LU 分解，那么求解方程组b x =A 的问题就等价于求解两个三角形方程组 （1）b y =L ，求y ； （2）y x =U ，求x 。
而这两个三角形方程组的求解是容易的。
下面我们先给出这两个三角形方程组的求解公式；然后研究在LU A =或LU PA =时，U L ,的元素与A 的元素之间的直接关系。
4－0 三角形线性方程组的解法 设⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=nn n n l l l l l l L21222111， 11121222n n nn u u u u u U u ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦则b y =L 为下三角形方程组，它的第i 个方程为),2,1(11,22111n i b y l y l y l y l y l ii ii i i i i i ij j ij ==++++=--=∑假定0≠ii l ，按n y y y ,,,21 的顺序解得：⎪⎪⎩⎪⎪⎨⎧=+-==∑-=),,3,2(/111111n i l b y l y l b y ii i i j j ij i上三角形方程组y x =U 的第i 个方程为),2,1(11,n i y x u x u x u x u in in i i i i ii nij j ij ==+++=++=∑假定0≠ii u ，按121,,,,x x x x n n -的顺序求解得：⎪⎪⎩⎪⎪⎨⎧--=+-==∑+=)1,2,,2,1(/1n n i u y x u x u y x ii i n i j j ij i nnn n4－1 基本的三角分解法设矩阵n n ij a A ⨯=)(的各阶顺序主子式均不为零，由定理3.1，A 有LU 分解如下⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡nn nrn n rn rrr r n rn ra a a a a a a a a a a a a a a a 2121222221111211 LU u u u u u u u u u u l l l l l l nn rn rrn r n r nrn n r r =⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡= 22221112112121211111比较等式两边的第1行上的相应元素，得U 的第1行元素：),,2,1(11n i a u i i == （4.4）比较等式两边的第1列上的相应元素（从第2开始），得L 的第1列元素：),,3,2(1111n i u a l i i ==（4.5） 比较等式两边的第2行上的相应元素（从第2开始），得i i i u u l a 21212+=故得U 的第2行元素：),,3,2(12122n i u l a u i i i =-=比较等式两边的第2列上的相应元素（从第3开始），得2222121i i i a u l u l =+，得L 的第2列元素：),,4,3(2212122n i u u l a l i i i =-=假设已求得U 的第1至1-r 行，L 的第1至1-r 列，比较等式两边的第r 行上的相应元素（从第r 开始），得rnrn n r r r n r n r r r r r r r r r r r r r rrrr r r r r r r r r a u u l u l u l a u u l u l u l a u u l u l u l =++++=++++=++++--+++--++--,11,22111,1,1,11,1,221,11,11,2211即),,1,(11n r r i a u u l riri r k ki rk +==+∑-=故得U 的第r 行元素：),,3,2;,,1,(11n r n r r i u l a u r k ki rk ri ri =+=-=∑-= （4.6）比较等式两边的第r 列上的相应元素（从第1+r 开始），得nrrr nr r r r n r n r n r r rr r r r r r r r r r r r r rr r r r r r r r r r r a u l u l u l u l a u l u l u l u l a u l u l u l u l =++++=++++=++++--++--+++++--+++,11,2211,2,2,11,222,211,2,1,1,11,122,111,1即),,2,1(,11n r r i a u l u l ir rr ir r k kr ik ++==+∑-=故得L 的第r 列元素：)1,,3,2;,,2,1(,11-=++=-=∑-=n r n r r i u u l a l rrr k krik ir ir （4.7）称由（4.4）－（4.7）式所表示的矩阵分解为Doolittle 分解。
它的计算顺序为：U 的第1行→L 的第1列→U 的第2行→L 的第2列→…→U 的第r 行→L 的第r 列→…→U 的第1-n 行→L 的第1-n 列→U 的第n 行。
(r 从1到n 循环即可！)实现了系数矩阵A 的Doolittle 分解后，求解方程组b x =A 的问题就等价于求解两个三角形方程组（1）b y =L ，求y ；（2）y x =U ，求x 。
利用4－0节有关结论，立得Doolittle 分解相应的求解公式为：⎪⎩⎪⎨⎧=-==∑-=),,3,2(1111n r y l b y b y r i i ri r r ⎪⎪⎩⎪⎪⎨⎧--=-==∑+=)1,2,,2,1(/1 n n r u x u y x u y x rr n r i i ri r r nnn n 类似地，可以推导Crout 分解的计算公式：),,2,1;,,1,(11n r n r r i u l a l r k kr ik ir ir =+=-=∑-=)1,,2,1;,,1(11-=+=-=∑-=n r n r i l u l a u rrr k kirk ri ri它的计算顺序为：L 的第1列→U 的第1行→L 的第2列→U 的第2行→…→L 的第r 列→U 的第r 行→…→L 的第1-n 列→U 的第1-n 行→L 的第n 行Crout 分解相应的求解公式为：⎪⎪⎩⎪⎪⎨⎧=-==∑-=),,3,2(/111111n r l y l b y l b y rr r i i ri r r ⎪⎩⎪⎨⎧--=-==∑+=)1,2,,2,1(1 n n r x u y x y x nr i i ri r r n n小结：Doolittle 分解计算：A ＝LU 1．计算U 的第1行元素： ),,2,1(11n i a u i i ==如果011=u ，则计算停止。
2．计算L 的第1列元素：),,3,2(1111n i u a l i i ==3．对于n r ,,3,2 =（1）计算U 的第r 行元素：),,1,(11n r r i u l a u r k ki rk ri ri +=-=∑-=如果0=rr u ，则计算停止。
（2）计算L 的第r 列元素：),,2,1(,11n r r i u u l a l rrr k krik ir ir ++=-=∑-=4．求解计算：y x b y ==U L ,⎪⎩⎪⎨⎧=-==∑-=),,3,2(1111n r y l b y b y r i iri r r⎪⎪⎩⎪⎪⎨⎧--=-==∑+=)1,,2,1(/1n n r u x u y x u y x rrnr i iri r r nn n n例：用Doolittle 法解方程组⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------725101391444321131243301024321x x x x 直接利用公式计算，从中总结 Doolittle 分解的计算规律（图示中还可清楚地看到存贮情况）： Gauss 消去法逐步对A 施加变换，逐步获得L ，U ，而Doolittle 分解公式则集中了前r-1次变换，一次性付诸实施。
U 的第1行与L 的第1列：⎥⎥⎥⎥⎦⎤⎢⎢⎣nn nn a a a a 2222U的第2行与L的第2列：⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡nn nin n i l l l u 2112111⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎣nn nin n in iin i n ia a a l a a u u u u212222211 U的第3行与L 的第3列：⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡nn ni n n n in iii i i n in i n ia a a l l a a a l l a a a l l u u u u l u u u u u 321321333332312223222111131211⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎣nn nin n n in iin in i n ia a a l l a a u u u u u u 321332211 U的第r 行与L 的第r 列（1,,3,2-=n r ）：⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎣---nn ninrr n n in iirn rin r i r n ia a a l l a a u u u u1,1,1,111最后计算U 的第n 行（实际上只有一个元素nn u ）。
注意在计算U 和解y b 看作A 的第1+n 列，它们遵循着相同的规则。
事实上，由y b L LU A ==,，得),(),(1b y A L U -=，据此，我们可得到紧凑格式的Doolittle 分解：即将上述分解过程处理的对象由A 改为增广矩阵),(b A ，用同样的方法处理第1+n 列，从而在获得系数矩阵的Doolittle 分解U L ,的同时，存放在增广矩阵第1+n 列中的元素就是y 。
(注意此时在公式（4.4）（4.6）中计算ri u 时，i 应到1+n 为止，编程计算时，应注意所有y x ,,,U L 的元素都存贮在增广矩阵中)紧凑格式的Doolittle 分解的算法设计： 小结：分解计算：A ＝LU 1．对于n r ,,3,2,1 =（1）计算U 的第r 行元素：)1,,1,(1111++=-=-=←∑∑-=-=n r r i a a a u l a u a r k ki rk ri r k ki rk ri ri ri如果0=rr a ，则计算停止。