矩阵QR分解途径的研究.pdfVIP

内江师范学院学报 第22卷第4期 ·18· OF TEACHERSCOLLEGE No．4V01．22 JOURNALNEIJIANG 矩阵Q尺分解途径的研究 刘 秀 梅 (连云港师范高等专科学校数学系， 江苏 连云港222006) 摘 及对矩阵进行列初等变换等方法进行． 关键词：矩阵；QR分解；初等变换；标准正交化 中图分类号：0151。21文献标识码：A 文章编号：167l—1785(2007)04一0018一03 将矩阵A分解为一个正交矩阵Q(或列正交矩阵Q) 若记d—lD．1(U7X)，则y—x一彬． 和一个上三角矩阵R的乘积，称为矩阵A的正交三角分 算法如下(计算y一日X，r一1) 解，简称Q尺分解．下面给出4种求矩阵Q尺分解的方法， ● (1)盯一sing(z，)(∑z；)÷， 以加强对QR分解思想及方法的深刻理解． (2)“】=z1+盯，甜，一z，(i一2，3，…，挖)， 1 利用Householder矩阵变换Ⅲ ￡，=(“1，甜2，…，“。)了’， (3)lD一盯“l， 将矩阵A的列向量依次实施Householder矩阵变换， (4)口一Jp一1(U7X)， 简记H，使之化为以具有1个非零元，2个非零元，…，n个非 (5)y—X一础． 零元作为列向量的上三角矩阵R．即若有H。…日。日。A 在A为n×1矩阵情况下，R=y，Q一日． 一欠，则Q一舅1臀2…H。一1． 1．2 若A为n阶矩阵 1．1 若A为n×1矩阵 设A一(X。，X：，…，X。)为”阶矩阵，方法如下： 取以为x=(z1，z2，…，z。)7，∞×1阶矩阵，y=(z1， z2，…，z，一1，一口，O，…，O)7∈印，其中 (1)d，一sign(z。)(∑z；)专(r一1，2，…，)， 盯一sign(z，)(∑z}){， (2)“，=z，+田，甜。=O(i一1，2，…，r一1)， … “，一z，(i—r+1，…，咒)，U=(“1，“2，…，“。)r， (当时z，；o，规定Sign(O)一1)． 若考虑r=】，y有一个非零元，令 (3)肛一盯，“，， (4)口=lD_1(U7X，)， 盯=sign(铂)(∑z；)专， (5)y，=X，一dU， t一1 则 (6)R一(yl，y2，…，y。)． I7=x—Y 例l 设矩阵 A—f：一：一H 则 【2 —4 5J 利用Householder变换作矩阵A