QR方法的探讨(07-5)依拉木江苏甫尔20070109058.docVIP

　　　　　　　　　　　　　　　　　　　　　编号 学士学位论文 QR方法的探讨 学生姓名： 依拉木江·苏甫尔 学 号： 20070109058 系 部： 数学系 专 业： 信息与计算科学 年 级： 07-5班 指导教师： 阿米娜·沙比尔 完成日期： 2012 年 5月 日 中文摘要 本文主要讨论矩阵的QR分解的三种方法,即Schmidt正交化、Givens变换、Householder变换等方法,具体构造正交矩阵Q。首先提出理论基础及算法，求出数值例题的QR分解。然后编写MATLAB程序，把如上数值例题在MATLAB环境下实现，得到的数值QR分解。最后通过比较，验证本文提出的三种方法的应用性可靠性。 关键词：QR分解；Schmidt正交化；Givens变换；Householder变换  目　录 中文摘要 1 引言 2 1 QR分解的概念 2 2 Givens方法 7 3豪斯霍尔德方法 （镜像变换） 15 3.1 Householder 矩阵 和Householder变换 15 3.2 QR算法 18 4 QR 分解的应用 20 总结 22 参考文献 23 致谢 24 引言 矩阵A的QR分解与基于QR分解的QR方法是求矩阵A的特征值的一种有效方法。对阶非奇异矩阵A的个列向量，使用Schmidt正交化程序可以得到A的QR分解；使用Givens变换和 Householder变换也可以直接将正交矩阵Q具体构造出来。 本文主要讨论得到A的QR分解的三种方法，即Schmidt正交化，使用Givens方法和 Householder方法等。任何实(复)非奇异阶矩阵可以分解成正交(西)矩阵和实(复)非奇异上三角矩阵的乘积. 1 QR分解的概念 定理1.1 如果实（复）非奇异矩阵能化成正交（酉）矩阵与实（复）非奇异上三角矩阵的乘积， 即 （1.1） 则称式（1.1）是的分解。 引理1.1 任何实的非奇异阶矩阵可以分解成正交矩阵和上三角矩阵的乘积，且除去相差一个对角线元素之绝对值全等于1的对角矩阵因子外，分解式（1.1）是唯一的。 证明：首先证明存在性 设的各列向量依次，由于非奇异，所以 线性无关，将它们按照施密特正交化法正交化，得到个标准正交的向量，且 这里都是常数，且由正交过程知。写成矩阵形式有,即 其中 是上三角矩阵。显然，可逆，而且也是上三角矩阵；由于的各列标准正交，所以为正交矩阵。从而有。 现在要证明唯一性,设有两种形式如（1.1）的分解式： （1.2） 其中和都是正交矩阵，和都是非奇异上三角矩阵，由式（1.2）得，式中仍为实非奇异上三角矩阵,于是 （1.3） 设 代入式（1.3）并与单位矩阵相比较，得 从而有 ，即 ， 这表明不仅是正交矩阵，而且还是对角线元素的绝对值全为1的对角矩阵，再由式（1.2）不难得 ， 显然，当规定上三角矩阵和对角线上元素为正实数时，则,从而分解唯一。 设 写出的分解。 解：用施密特（Schmidt）方法，令 单位化后得 所以 定理1.2 设为复矩阵（），且个列向量线性无关，则具有分解 （1.4） 其中是复矩阵，且满足 ，是阶复非奇异上三角矩阵，且除去相差一个对角元素的模全为1的对角矩阵因子外，分解式（1.4）唯一的。 施密特（Schmidt）方法的MATLAB程序如下： function c=schmidt_sum(A,y1,n) c=zeros(length(A),1); for i=1:n-1 xx=y1(:,i) *A(:,n); yy=y1(:,i)*y1(:,i); c=c+(xx/yy)*y1(:,i) ; end y1(:,n)=A(:,n)-c; c=y1(:,n); % end function A=[1 1 0;1 -1 1;0 0 2]; m=length(A); y1=zeros(size(A)); y=y1; y1(:,1)=A(:,1); for i=2:m y1(:,i)=schmidt_sum(A,y1,i); en