复矩阵的Givens变换及其QR分解.doc

 ｇｙ ［文章编号］１６７１－９７２７（２０１１）０６－０６９３－０４ 复矩阵的 Ｇｉｖｅｎｓ变换及其ＱＲ 分解 杜 鹃 冯思臣 （成都理工大学 管理科学学院，成都 ６１００５９） ［摘要］实矩阵有成熟的三角分解算法，复矩阵尚无好的三角分解算法。 为解决复矩阵的三角 分解与ＱＲ 分解问题，采用科学类比，重新拓展定义，演绎计算的方法，给出复 Ｇｉｖｅｎｓ矩阵的定 义，推导出了复 Ｇｉｖｅｎｓ矩阵是酉矩阵，得到了用有限个复 Ｇｉｖｅｎｓ变换将一个ｎ 维复向 量 旋 转 到任何一个给定方向的方法，证明了任何一个非奇异复矩阵能够通过有限次复 Ｇｉｖｅｎｓ变 换， 分解为一个酉矩阵与一个复非奇异上三解矩阵的乘积，利 用 复 Ｇｉｖｅｎｓ变 换 解决了复矩阵的 ＱＲ 分解问题。 ［关键词］复 Ｇｉｖｅｎｓ变换；Ｇｉｖｅｎｓ矩阵；ＱＲ 分解 ［分类号］Ｏ１５１．２１ ［文献标志码］Ａ １ 预备知识 随着现代化科学技术的迅速发展，矩 阵 的 分 解在控制理论、信息论、系统识别和信息处理、优 化理 论、最 小 二 乘 问 题 中 都 是 十 分 重 要 的 工 具。 参考文献［１～４］涉 及 到 矩 阵 的 ＱＲ 分 解，并 且 目 前的大多数文 献 中 也 只 对 实 矩 阵 利 用 Ｇｉｖｅｎｓ矩 阵变换、Ｈｏｕｓｅｈｏｌｄｅｒ矩阵变换、Ｄｏｏｌｉｔｔｌｅ分解得 到ＱＲ 分解公式。文献［５～８］中给出了一些更好 的算法途径，而对复矩阵的 Ｇｉｖｅｎｓ矩阵变换及其 ＱＲ 分 解，尚 没 有 具 体 方 法。 但 这 一 问 题 在 工 程 技术应 用 中 是 非 常 有 实 用 价 值 的。 本 文 就 从 复 Ｇｉｖｅｎｓ矩阵变换入 手，给 出 复 矩 阵 的 ＱＲ 分 解 方 法。 定义１ 设实数ｃ 与ｓ 满 足ｃ２ ＋ｓ２ ＝１，称 矩 阵 熿 燄  １ ｃｅｊθ１ ｓｅｊθ２ （ｉ） １ Ｕｉｋ ＝  １ ｃｅｊθ４ （ｋ） －ｓｅｊθ３ １  燀 燅 ［ｋ］ 为复数形式的 Ｇｉｖｅｎｓ矩阵（初等旋转矩阵），其中 ｃ＝ｃｏｓθ＞θ＞θ 为旋转角，θ１ ＋θ４ ＝θ２ ＋ θ 。可以得到 Ｕ ＝ｅｊ（θ ＋θ ）当 ｎπ 时 θ４ ＝ －θ１ ＋２ １ ４ ３ ｉｋ ｄｅｔＵｉｋ ＝１。 定义２ 如果实（复）非奇异矩阵 Ａ 能够化成 正交（酉）矩 阵 Ｑ 与 实 （复）非奇异上三角矩阵 Ｒ 的乘积，即Ａ＝ＱＲ，则称为Ａ 的ＱＲ 分解。 ［收稿日期］２０１０－０５－２８ ［基金项目］国家自然科学基金资助项目 ；四川省教育厅自然科学重点资助项目 （０８ＺＡ１１４） ［作者简介］杜 鹃 （１９６１－），女 ，副 教 授 ，主 要 研 究 方 向 ：数 值 代 数 ，Ｅ ｍａｉｌ：ｄ４０７８ １２６．ｃｏｍ。 性质１ ＵＨＵ ＝Ｉ，即Ｕ 是酉矩阵，且ＵＨ 也 则有 ，ηｋ 。 ηｉ＝ 槡 ＞ ＝０ ｉｋ ｉｋ ｉｋ ｉｋ 是复 Ｇｉｖｅｎｓ矩阵。 性质２ 当θ４＝－θ１，θ３＝－θ２ 时，则 证 可由性质３推得。 主要结论 ２ １ 熿 燄 定理１ 设ｘ＝（ξ１，ξ２，…，ξｎ）Ｔ ≠０，ｘ∈Ｃｎ ，则 １ 存在有限个复数形式的 Ｇｉｖｅｎｓ矩阵的乘积，记为 ｃｅｊθ１ ｓｅｊθ２ Ｕ，使Ｕｘ＝‖ｘ‖ ｅ ２ １。 １ 证 若ξ１≠０，构造复 Ｇｉｖｅｎｓ矩阵Ｕ１２， Ｕｉｋ ＝  令ｃ＝ ，ｓ＝ ， １ ｃｅ－ｊθ１ 槡 槡 θ２ １＝ －ａｒ ξ１， ２＝－ａｒ ξ２，由性质４可得Ｕ１２ｘ＝ θ ｇ θ ｇ １ （槡 ，０，ξ３，…，ξｎ ）Ｔ ，再 对 Ｕ１２ｘ 构 造  复 Ｇｉｖｅｎｓ矩阵Ｕ１３， １燅 燀 槡 ， 以下讨论都假设θ４＝－θ１，θ３＝－θ２。 槡 对于不全为零的复数α 和β，可选取 ξ３ ， ｓ＝ ，ｓ＝ ， ｃ＝ 槡 槡 槡 θ２＝－ａｒｇξ３， θ１＝－ａｒｇ槡 θ１＝ －ａｒｇα，θ２ ＝ －ａｒｇβ，则αｃｅｊθ１ ＋βｓｅｊθ２ ＞ αｓｅ－ｊθ２ ＋βｃｅ－ｊθ１ ＝０。 证 αｃｅｊθ１ ＋βｓｅｊθ２ Ｕ１３（Ｕ１２ｘ）＝ （槡 ，０，０， ， ξ４ …，ξｎ）Ｔ ，如此继续下去，最后一次对Ｕ１ｎ－１ …Ｕ１２ｘ 构造Ｕ１ｎ ， ２ － ｅｊθ１ ＋β ｅｊθ２ ， ＝α ｃ＝ 槡 槡 ２ ２ 槡 ｎ αα珔 ββ珋 ＝ ＋ ｓ＝ 槡 槡 槡 ２ ， ＞ θ１ ＝ －ａｒｇ 槡 ＝ 槡 ＝ 槡 θ ＝－ａｒｇξｎ ， Ｕ１ｎ（Ｕ１，ｎ－１…Ｕ１２ｘ）＝ 同理 －αｓｅ－ｊθ２ ＋βｃｅ－ｊθ１ （槡 ，０，０，…， ０）Ｔ ，令Ｕ＝Ｕ１ｎＵ１，ｎ－１…Ｕ１２，则Ｕｘ＝ ‖ｘ‖２ｅ１，如 ｅ－ｊθ２ ＋β ｅ－ｊθ１ ＝－