第3讲 线性方程组的数值解法_shopo.ppt

按矩阵乘法展开，可逐行求出分解矩阵L的元素，计算公式是对于i=1,2,…,n … … … 移项便得Cholesky分解公式 j=i+1,i+2,…,n 因对称性无需存储 Step1 Step2 Step3 Stepn 的计算过程: 逐 列 计 算 元素 仍然存放在矩阵 的相应位置上 乔累斯基（Cholesky）分解所需要的乘除次数约为 数量级，比LU分解节省近一般的工作量 其分解法的计算过程分为两个步骤： 1. 将对称正定矩阵A分解，有A=LLT, 算出L。 2. 将A= 代入线性方程组Ax=b，有 x=b。因 此求解两个三角形方程组Ly=b及 相应的求解公式 i=1,2,…,n i=n,n-1,…,2,1 例3.12用Cholesky分 解法求解下列方程组 解： 系数矩阵为 Step1 Step2 Step3 求解方程组 求解方程组 3.5 追赶法 常常会遇到三对角线性方程组 其中 ① ② 则 非奇异，且 如果三对角矩阵A 为对角占优矩阵，则可以采用追赶法求解。 三对角矩阵 满足条件 系数A矩阵可作三角分解： 按乘法展开 则可计算 可依次计算 = = 分解为 “追”的过程 方程组求解的计算公式： ①解方程组 ②解方程组 “赶”的过程 a