计算方法方程组b.pptVIP

2.2 Tridiagonal System 2.2　追赶法 将高斯消去法 直接应用于 三对角方程 三对角方程组 /* Tridiagonal Linear System */ Step 1: 追——即消元 Step 2: 赶——即回代： 一旦bi-ui-1ai=0算法中断 2.2 Tridiagonal System Hey, what does diagonally dominant mean??? It means that the diagonal entries of the matrix are very LARGE. Well, how large is LARGE? They satisfy the following inequality: 证明：只需证明公式中分母bi-ui-1ai≠0,i=2,3,…,n-1。设ki=bi-ui-1ai,根据右式和对角占优条件，有 定理 若 A 为对角占优 /* diagonally dominant */ 的三对角阵，既满足 　　　，则追赶法可解以 A 为系数矩阵的方程组。 2.2 Tridiagonal System 例 求解三对角方程 i ai bi ci di bi- uiai ui yi xi 1 0 2 -1 1 -－ -0.5 0.5 1 2 -1 2 -1 0 1.5 -.667 .333 1 3 -1 2 -1 0 1.333 0.75 0.25 1 4 -1 2 0 1 1.25 - 1 1 解： 2.2 Tridiagonal System 2.3 Matrix Factorization 如果线性代数方程组能化成以下形式 启示 2.3　三角分解法 /* Matrix Factorization */ 上三角阵Upper-triangular 下三角阵Lower-triangular Could we factorize A into triangular matrix directly? Sure! for example Doolittle Factorization, Crout Factorization etc. ? 高斯消元法的矩阵形式 /* Matrix Form of G.E. */： Step 1: 记 L1 = ，则 Step n ? 1: 其中 Lk = 2.3 Matrix Factorization – Matrix Form of G.E. 记为 L 单位下三角阵 /* unitary lower-triangular matrix */ 记 U = A 的 LU 分解 /* LU factorization */ 2.3 Matrix Factorization – Matrix Form of G.E. 2.3 Matrix Factorization – Doolittle ? 道立特分解法 /* Doolittle Factorization */： —— LU 分解的紧凑格式 /* compact form */ 通过比较法直接导出L 和 U 的计算公式。 思路 根据矩阵相乘法则，有 lii=1, i=1,2,…,n 2.3 Matrix Factorization – Doolittle u11=a11, u12=a12, …, u1n=a1n l21=a21/u11 l31=a31/u11 …… …… ln1=an1/u11 u22=a22－l21u12, … , u2n=a2n－l21u1n l32=(a32－l31u12)/u22 …… ln2=(an2－ln1u12)/u22 uij lji 2.3 Matrix Factorization – Doolittle Algorithm: Doolittle Factorization Step 1: u1j = a1j ; lj1= aj1/u11; ( j = 1, …, n ) Step 2: compute uij and lji for i = 2, …, n?1; uij lji 2.3 Matrix Factorization – Choleski ? 平方根法 /* Choleski’s Method */： ——对称 /* symmetric */