数值代数qr方法实验报告.docVIP

一、引言 很多情况下在工程中抽象出来的数学方程组是超定的，没有精确解，这样就需要找一个在某种意义下最接近精确解的解。设A是m*n的实矩阵，,这样就等价与,为方便求解，需要是上三角矩阵，这样引入分解就比较求解方便。 二、豪斯霍尔德变换 豪斯霍尔德变换（Householder transformation）又称初等反射（Elementary reflection），最初由A.C Aitken在1932年提出，Alston Scott Householder在1958年指出了这一变换在数值线性代数上的意义。这一变换将一个向量变换为由一个超平面反射的镜像，是一种线性变换。其变换矩阵被称作豪斯霍尔德矩阵，在一般内积空间中的类比被称作豪斯霍尔德算子，超平面的法向量被称作豪斯霍尔德向量。 三、理论依据 householder变换：任意的一个向量，经过一个正交变换 后总可以变为一个与之范数相等的另一个向量。 如上图中所示，记，，上述H既为所要求的householder变换。 具体操作时将需要变化的矩阵的每一列当做一个向量，第一列变为除了第一个元素都为0的向量…第i列变为除了前i个元素都为0的向量...为了保证在每次变化时不改变已经变好的0元素，第i次只变化每列第i个元素到第n个元素，每个变化矩阵的形式是[I 0;0 H]。 算法如下：qr_Houserhoilder.m funct