第六章 线性方程组的数值解.pptVIP

§1 引言 §2. Gauss 消去法 二.一般线性方程组的Gauss 消去法计算过程 第k-1次消元得到 k次消元 回代过程算法 定理1 (高斯消去法)设Ax=b，其中A∈Rn×n。如果约化的主元素a kk≠0(k=1,2…,n)，则可通过高斯消去法(不进行交换两行的初等交换)将方程组Ax=b约化为三角形矩阵方程组且消元和求解公式为 四、 Gauss消去法乘法计算量 (1)消元计算：在第k步(k=1,2….n-1) 1、计算乘数：需作n-k次除法运算； 2、消元：需作有(n-k)2次乘法运算 3、计算b (k):需作(n-k)次乘法运算； 6.3 选主元高斯消去法 第k步 从 的第k列 ， ， 中选取绝对值最大项，记录所在行，即 若 交换第k行与 行的所有对应元素，再进行顺序消元。 具体如下: 第1步(k=1)：首先在A中选主元素，即选择i1, j1使 再交换(A,b)的第1行 与第i1行元素，交换A的第1列与第j1列元 素，将ai1 j1调到(1,1)位置(交换后增广阵为简 单起见仍记为(A,b))；然后,进行消元计算。 矩阵的三角分解 §4 三角分解法 二.Doolittle分解法 ： 2.一般计算公式 三.LU 分解求解线性方程组 四.例题 定理3：设A为非奇异阵，则必存在置换矩阵P，使得 其中L为单位下三角阵，U为非奇异的上三角阵。 §5 求解三对角方程组的追赶法 §5 求解三对角方程组的追赶法 在实际问题中如样条插值及常微分方程边值问题的数值解中，会遇到求解三对角线形线组： 追赶法公式实际上就是把高斯消元法用到求解三对角线方程组上去的结果，这时由于A特别简单，因此使得求解的计算公式非常简单，而且计算量仅有5n-4次乘除法。3n-3次加减法，工作量小，电算时，需要3个一维数组存储A的系数，两个一维数组保存中间结果。 例6 用追赶法解方程组 练习: 用追赶法解线组 §6对称正定阵的平方根法 设有方程组Ax=b，其中，A∈Rn×n。若A满足下述条件，则称A为对称正定矩阵。 回顾：对称正定阵A的几个重要性质 （1）A?1 亦对称正定，且 aii 0 （2）A 的顺序主子阵 Ak 亦对称正定 （3）A 的特征值 ?i 0 (i=1,2,…,n) （4）A 的全部顺序主子式 det ( Ak ) 0(i=1,2,…,n) 三 对称正定阵的LDLT分解 设A为对称正定矩阵，则由定理知，A有惟一的三角分解 于是，对角阵D还可以分解为 定理6(对称正定阵的三角分解) 设A为n阶对称正定矩阵，则有三角分解： ①A=LDLT，其中L为单位下三角阵，D为对角阵，或 ②A=LLT，其中，L为下三角阵且当限定L的对角元素为正时，这种分解是惟一的，这种矩阵分解称为(Cholesky)分解。 由矩阵乘法，则有L的第1列元素 同理，可确定L的第j列元素lij(i=j,…,n)。 由分解公式有 所以 例7:用平方根法求以下方程组的解. Cholesky分解法要用到开方运算，为避免开方运算，可将A分解为A=LDLT(其中L为单位下三角矩阵)，再分别解方程组LY=b及DLTX=Y或，这种方法称为改进平方根法. §7.向量和矩阵的范数 为了研究线性方程组近似解的误差估计和迭代法的收敛性，我们需要对Rn(n维向量空间)中的向量或Rnxn中矩阵的“大小”引入一种度量，——向量和矩阵的范数。 定义7.2 设x=(x1, x2,…, xn)T∈Rn，定义Rn上3种常用的向量范数 定义3 (向量序列的极限) 设{x(k)}为向量序列，记x(k)=(x1(k), x2(k),…, xn(k))T∈Rn及x*=(x1*,…, xn*)T ∈Rn。如果n个数列极限存在且 则称{x(k)}收敛于x*，记为 定理7 设{x(k)}是Rn中一向量序列，且x*∈Rn，则 证只就υ=∞证明。 显然有 1.定义7.4 (矩阵的范数) 2.相容范数 3.算子范数 三 例题 §8 线性代数方程组的迭代解法 例题11：用高斯-塞德尔迭代公式(GS)解方程组 10 病态方程组和迭代法的改善 迭代结果分析：并不是每一个迭代公式所构造的迭代序列都收敛。于是，计算方法的目的就是要寻求一个使得构造序列收敛的方法。因此产生了各种各样的迭代方法，根据迭代矩阵的不同构造方法，形成了不同的迭代方法。这里介绍两种迭代方法：雅可比迭代方法、高斯-塞德尔迭