数值分析课件 第5章6-7节.pptVIP

5.6 误差分析 5.6.2 迭代改善法 5.7 矩阵的正交三角化及应用 5.7.2 平面旋转矩阵 5.7.3 矩阵的QR分解 5.7.4 求解超定方程组 5.7.1 初等反射阵 定义9 设向量 且 ， 为初等反射阵(或称为豪斯霍尔德（Householder）变换). 如果记 , 则 称矩阵 定理25 设有初等反射阵 , (1) 是对称矩阵，即 (2) 是正交矩阵，即 (3) 设 为对称矩阵，那么 也是对 称矩阵.  证明 只证 的正交性，其他都可通过验证得到. 其中 则： 是一个初等反射阵. 初等反射阵的几何意义. 考虑以 为法向量且过原点 的超平面 . 设任意向量 ， 则 , 其中 . 设向量 , 则显然 于是 图5-1 对于 ， 从而对任意向量 ， 其中 为 关于平面S的镜面反射(见图5-1). 总有 定理26 设 为两个不相等的 维向量， 则存在一个初等反射阵 ， 证明 令 ， 而且 使 则得到一个初等反射阵 因为 所以 是使 成立的惟一长度等于1的向量(不计符号). 定理27 设 ， 则存在初等反射阵 , 使 ， (约化定理) 其中 证明 记 , 设 , 取 ， 于是由定理22存在 变换 其中 , 使 记  其中 则有 于是 显然 如果 和 异号，那么计算 时有可能出现两相 近数相减的情况，有效数字可能损失. 取 和 有相同的符号，即取 在计算 时，为了避免上溢或下溢，将 规范化 则有 使 , 其中 算法6 设 ， 及 , 本算法计算 的分量冲掉 的分量. 使 关于 的计算， 设 , 为 的第 列向量， 其中 则 因此计算 就是要计算 于是计算 只需计算两向量的数量积和两向量的加法. 计算 只需作 次乘法运算. 设 , 则变换 是平面上向量的一个旋转变换，其中 为正交矩阵. 其中 中变换： 而 称为 中平面 的旋转变换(或称为吉文斯（Givens） 变换)， 称为平面旋转矩阵. 显然， 具有如下性质： (1) 与单位阵 只是在 位置 元素不一样，其他相同. (2) 为正交矩阵 (3) (左乘)只需计算第 行与第 行元素， 即对  有 其中 (4) (右乘)只需计算第 列与第 列元素 利用平面旋转变换，可使向量 中的指定元素变为零. 其中 证明 取 . 由  , 定理28 设 , 其中 不全为零， 则可选择平面旋转阵 (约化定理) 使 利用矩阵乘法，显然有 由 的取法得 设 且为非零矩阵，则存在初等反射阵 设 使 如果 , 取 ， (1) 第1步约化： 即这一步不需要约化， 不妨设 ， 于是可选取初等反射阵 , 于是 使 其中 (2) 第 步约化： 设已完成对 上述第1步到第 步的约化，再进行第 步约化. 即存在初等反射阵 使 其中 其中 为 阶上三角阵，