2迭代改善法教学课件.pptVIP

5.6.2 迭代改善法 5.7 矩阵的正交三角化及应用 5.7.2 平面旋转矩阵 5.7.3 矩阵的QR分解 5.7.4 求解超定方程组 (2) 由设及定理29，存在初等反射阵 记 , 即 其中 为正交矩阵， 为上三角阵. 证明 (1) 由定理29可得. 使 则上式为 唯一性. 设有 其中 为正交阵， 为非奇异上三角阵，且 具有正对角元素， 由假设及对称正定矩阵 的Cholesky分解的唯一性， 则 . 也可以用平面旋转变换来约化矩阵. 则 从而可得 定理31 (1) 存在正交矩阵 设 为非奇异矩阵， (用吉文斯变换计算矩阵的QR分解) 则 使 (2) 有QR分解 其中 为正交阵， 为非奇异上三角阵，且当 对角元素 都为正时，分解是唯一的. 证明 由设有 使 . 若 ， (1) 第1步约化： 则可选择吉文斯变换 使 可简记为 , 其中 (2) 第 步约化： 设上述过程已完成第1步至第 步， 于是有 由设有 使 ， 若 ， 则可选择吉文斯变换 ， 使 其中 (3) 继续上述约化过程，最后则有 其中 为正交阵(为一系列平面旋转阵的乘积). 记 , 则 有QR分解 其中 为正交矩阵， 为非奇异上三角阵. 用 变换实现 的正交三角约化需要 次乘法， 而用吉文斯变换约化计算约需要 次 次开方运算， 乘法， 次开方运算. 这说明，对一般矩阵 用吉文斯变换实现 的 正交三角约化比用 变换实现 的正交三角约化，计算量 要大一倍. 但是，如果 为三对角阵或上海森伯格阵，则需要约 化为零的元素比较少，这时用吉文斯变换实现 的正交三 角约化就显得简单、经济. 设 ，其中 . 当 时称为超定方 线性最小二乘问题： 寻求 使 如果使上式成立的 存在，称 为超定方程组最小 二乘解. 程组， 一般地，没有通常意义下的解. 对超定方程组， 记残差向量 ， 于是 设 ，其中 为非奇异矩阵，且为病态 方程组(但不过分病态). 本节研究的问题是，当求得方程组的近似解 后， 首先用选主元三角分解法实现分解计算 其中 为置换阵， 为单位下三角阵， 为上三角阵， 且求得计算解 . 如何对其精度进行改善. 计算剩余向量 （6.14） 现利用 的剩余向量来提高 的精度. 求解 , 得到的解记为 . （6.15） 然后改善 如果(6.14)，(6.15)及解 的计算没有误差，则 说明 就是 的精确解. 但在实际计算中，由于有舍入误差， 只是方程组的 近似解，重复(6.14)，(6.15)的过程，就产生一近似解序列 ，有时可能得到比较好的近似. 算法5 设 ，其中 为非 奇异矩阵，且 为病态方程组(但不过分病态)，用选 主元分解法得到 及计算解 . 本算法用迭代改善法提高近似解 的精度. (迭代改善法) 设计算机字长为 ，用数组 保存 元素，数组 保存三角矩阵 及 ，用 记录行交换信息， 存贮 及 ,  保存 或 . 1. 用选主元三角分解实行分解计算 且求计 算解 (用单精度) 2. 对于 (1) 计算 (用原始 及双精度计算) (2) 求解 , 即 (用单精度) (3) 如果 则输出 , 停机