雅克比矩阵知识介绍

1、雅可比矩阵Jacobi方法)Jacobi方法Jacobi方法是求对称矩阵的全部特征值以及相应的特征向量的一种方法，它是基于以下两个结论1) 任何实对称矩阵A可以通过正交相似变换成对角型，即存在正交矩阵Q,使得QtAQ=diag(九，入，入)(3.1)12n其中九(i=1,2,n)是A的特征值，Q中各列为相应的特征向量。i2) 在正交相似变换下，矩阵元素的平方和不变。即设A=(aij)x,Q交矩阵,记B=QiAQ=(b)x,则八“ijnxnnn三Ej二2陽Jacobi方法的基本思想是通过一次正交变换，将A中的一对非零的非对角化成零并且使得非对角元素的平方和减小。反复进行上述过程,使变换后的矩阵的

2、非对角元素的平方和趋于零，从而使该矩阵近似为对角矩阵，得到全部特征值和特征向量。1矩阵的旋转变换设A为n阶实对称矩阵，考虑矩阵记易见V.3)是正交矩阵,注意到B=VA的第i,j行元素以及的第i,j列元素为ij可得诂=丐icos<b+幻y<b十丐于sinl<b|a)*?*Ajj-殆jsin妙ha力cos4>一a巧sinl妙如果a和，取使得c'<-7则有,这样继续下去，得到一个矩阵序列A（k）。可得A(2)对A（1）重复上述的过程，可以证明，虽然这种变换不一定能使矩阵中非对角元素零元素的个数单调增加但可以保证非对角元素的平方和递减，我们以A与A（i）为例进行讨

3、论。E二E.设由式(3.4)可得应）二Z+2£（智乜針）喘二砖尸=0二刀希+2工+咗）二E-喝U）klki,jk,1fj这表明,在上述旋转变换下,非对角元素的平方和严格单调递减,因而由（3.2）可知，对角元素的平方和单调增加。2.Jacobi方法通过一系列旋转变换将A变成A(k+i),求得A的全部特征值与特征向量的方法称为Jacobi方法。计算过程如下1)令k=0,A(k)=Amax1W2,页刍?础二2) 求整数i,j,使得3) 计算旋转矩阵4)计算A(k+1)(fr+l:%(fr+1:%=古如1-a*u=a°=0硏1j2(fr)©%+d2S丄OJ3初+2dcj住

4、十1)i2gd%十/3r>j閃切-2血切3.鼻聞-<7码7十»+i;=-谢+竝呂JEUf=E右餐5) 计算6) 若E(A(k+i)<g则为特征值，Qt=（V（0）VV（k+1）T的特征值和特征向量。一般地，Jacobi法不能在有限步内将A化成对角阵，但有下面的定理。定理3设A为n阶使对称矩阵，对A用Jacobi法得到序列A（k）,其中A（o）=A,则Ji®疋防）二0证明由Jacobi法计算过程故有>1M閃）（3.5）有计算A的公式可以得到另一方面，-a_c十21斟力(fr-J-J(a+2)+1|7XI/1QJJ3J?aac/Is+1于是有,代入式(3

5、.5)得)W(l厅因为一爲6所以TimE(A>Wk>科1-/FC4)盘too雅可比矩阵以m个n元函数u=u(x,x,x)(i=1,2,m)的偏导数(j=1,2,n)ii12n为元素的矩阵血13肚18諷928巧饥%dx2如果把原来的函数组看作由点x=(x,x,，x)到点u=(u,u,，u)的一个12n12m变换T,则在偏导数都连续的前提之下，u随x的变化由相应的微分方程组来描述。这是一个关于微分的线性方程组，其系数矩阵便是雅可比矩阵J),因而可写成矩阵形式这隐含着(J)具有微分系数的某些性质，类似于一元函数的导数。而在m=n=1的情形,它又恰好是一个一元函数的导数；所以它也是一个一元

6、函数的导数到m个n元函数的一种推广。因此，(J)作为微分系数或导数的推广，有时也被当作变换T的“导数”看待并记为T卜(x)=(J)。变换T的进一步的数量描述需要雅可比行列式。定义任给一个n维向量X,其范数凶是一个满足下列三个条件的实数:(1) 对于任意向量X,11X0,且凶=06X=0;(2) 对于任意实数九及任意向量X,甌Xl=|九I凶；(3) 对于任意向量X和Y,IIX+Y怪凶+MI；对于这样的,叫雅克比矩阵定义。雅克比矩阵证明关于这个的一般性证明稍微复杂点,现在就给你证明为什么二维的dx(u,v)dy(u,v)=Jdudv成立证明：对于曲面x=x(u,v),y=y(u,v),取它的微元，

7、即小曲边四边形ABCD,其中A(u,v),B(u+Au,v),C(u+Au,v+Av),D(u,v+Av),那么这个曲边四边形ABCD可以近似看成是微小向量B(u+Au,v)-A(u,v)和D(u,v+Av)-A(u,v)张成的。利用中值定理可知：(u+u,v)-(u,v)=Mdu(u,v+v)-(u,v)=Ndv这里的M,N是偏导数的形式,不好打出,你可以自己算出来,很简单的当变化量很小时，我们把(u+Au,v)-(u,v)近似看成dx(u,v),(u,v+Av)-(u,v)看成dy(u,v),所以，dx(u,v)dy(u,v)=M*Ndudv而其中的M*N刚好就是二维Jacobi行列式的展

8、开形式。由此问题得证。雅可比矩阵在向量微积分中，雅可比矩阵是一阶偏导数以一定方式排列成的矩阵，其行列式成为雅可比行列式。还有，在代数几何中，代数曲线的雅可比量表示雅可比簇：伴随该曲线的一个群簇，曲线可以嵌入其中。它们全部都以数学家雅可比命名；英文雅可比量"Jacobian"可以发音为ja'kobin或者屯kobin。雅可比矩阵的重要性在于它体现了一个可微方程与给出点的最优线性逼近。因此，雅可比矩阵类似于多元函数的导数。雅可比矩阵定义：雅可比矩阵定义为向量对向量的微分矩阵，定义式如下：见所附jpg图片。例：MATLAB中jacobian是用来计算Jacobi矩阵的函数

9、。symsrlfx=r*cos(l)*cos(f);y=r*cos(l)*sin(f);z=r*sin(l);J=jacobian(x;y;z,rlf)结果：J=cos(l)*cos(f),-r*sin(l)*cos(f),-r*cos(l)*sin(f)cos(l)*sin(f),-r*sin(l)*sin(f),r*cos(l)*cos(f)sin(l),r*cos(l),0Hessian矩阵就是一个多元实函数的二阶导数，设f=f(x1,x2.xn)二阶导数(d八2f/d(xi)d(xj)构成矩阵，在优化分析中常用到。Jacobi矩阵就是一个多元矢量函数的一阶导数，如仁(f1(x1,x2.

10、xn),.,fm(x1,x2.xn),相应矩阵元素为d(fi)/d(xj)。在稳定点附近的稳定性分析常用到它。本质上说，以上两者是相关的Jacobi可以看作是一个多元实函数的梯度(一阶导数)的导数。在向量微积分中，雅可比矩阵是一阶偏导数以一定方式排列成的矩阵，其行列式称为雅可比行列式。还有，在代数几何中，代数曲线的雅可比量表示雅可比簇：伴随该曲线的一个代数群，曲线可以嵌入其中。它们全部都以数学家卡尔雅可比命名；英文雅可比量"Jacobian"可以发音为jakobin或者屯kobin。雅可比矩阵雅可比矩阵的重要性在于它体现了一个可微方程与给出点的最优线性逼近。因此，雅可比矩阵

11、类似于多元函数的导数。假设F：RfR是一个从欧式n维空间转换到欧式m维空间的函数。这个函数由m个实函数组成:y1(x1,.,xn),.,ym(x1,.,xn).这些函数的偏导数(如果存在)可以组成一个m行n列的矩阵，这就是所谓的雅可比矩阵：此矩阵表示为：这个矩阵的第i行是由梯度函数的转置y(i=l,.,m)表示的如果p是R中的一点，F在p点可微分，那么在这一点的导数由J(p)给出(这是nF求该点导数最简便的方法)。在此情况下，由(p)描述的线性算子即接近点p的FF的最优线性逼近，x逼近与p例子由球坐标系到直角坐标系的转化由F函数给出：RX0,nX0,2n-Rs=rsin9sin©&#

12、163;3=rcosff此坐标变换的雅可比矩阵是R4的f函数:其雅可比矩阵为:;0002t)1252354511&11141MsMs血力如亦_-.!1/Th2rX/It.JF此例子说明雅可比矩阵不一定为方矩阵。在动力系统中考虑形为X'=F(x)的动力系统，F:RnfRn。如果F(x)=0,那么x是一个00驻点。系统接近驻点时的表现通常可以从J(x)的特征值来决定。F0雅可比行列式如果m=n，那么F是从n维空间到n维空间的函数，且它的雅可比矩阵是一个方块矩阵。于是我们可以取它的行列式，称为雅可比行列式。在某个给定点的雅可比行列式提供了F在接近该点时的表现的重要信息。例如，如果连续

13、可微函数F在p点的雅可比行列式不是零，那么它在该点附近具有反函数。这称为反函数定理。更进一步，如果p点的雅可比行列式是正数，贝吓在p点的取向不变；如果是负数，则F的取向相反。而从雅可比行列式的绝对值，就可以知道函数F在p点的缩放因子；这就是为什么它出现在换元积分法中。例子设有函数F：R3fR3，其分量为:如=5t2则它的雅可比行列式为：从中我们可以看到，当x和x同号时，F的取向相反；该函数处处具有反函数,12除了在x=0和灭=0时以外。12海森矩阵Jacobi阵和Hessian阵可谓是用途广泛。不仅是在优化问题中常用，在各种多元问题中一般都常遇到。比如在求解非线性方程组时，两者会经常被使用。关

14、于它们的定义和用法，推荐参看李庆扬等写的非线性方程组的数值解法。楼主可参看R.A.Horn的矩阵分析卷一雅可比矩阵雅可比矩阵定义：雅可比矩阵定义为向量对向量的微分矩阵，定义式如下:见所附jpg图片。雅克比矩阵在有限单元法中指的是全局坐标对局部坐标的偏导数，我们在由节点位移求解节点应变的时候会碰到形函数对整体坐标的偏导问题，即求解B矩阵，由于形函数对整体坐标的偏导比较难求，我们可以先求形函数对局部坐标的偏导数，然后利用雅克比矩阵转化为形函数对整体坐标的偏导数。雅克比矩阵其实质和泛函求导数中泛函中的变量对自变量的导数是一样的。只不过一个是矩阵，一个数数值而已。在空间问题8节点的线性单元中，雅克比矩

15、阵是3*3得矩阵。我们在计算过程中还经常用到雅克比行列式的值，它用来判断单元是否畸形，一般雅克比行列式为正，则说明单元形态较好，反之单元形态不好。雅克比方法是解决线性方程组的一种迭代方法，当现行方程组的阶数较高时，用直接法解方程组时可能误差较大，就要采用迭代方法。不过雅克比迭代不是一种很高效的迭代方法，高斯-塞德尔迭代方法较其效率要高。雅克比矩阵必然是n*n的矩阵，因为局部坐标和全局坐标之间的变量的数量永远是相同的在一般的应用过程中，局部坐标和整体坐标线性无关且数目相等，所以雅克比矩阵是方阵，并且行列式不恒为零。但是，从纯数学的角度讲，如果允许广义坐标之间线性相关，那么雅克比矩阵可能不是方阵，

16、即使是方阵，行列式也可能恒为零。讨论Jacobi矩阵不能仅限于有限元。应用数学上是这样定义的：设有n个变元的m个函数yi=fi(x1,x2,xn)(i=1,2,m)，A二D(y1,.,ym)/D(x1,.xn)称为上式的Jacobi矩阵。m和n是可以不同的。其实局部坐标和全局坐标的变量的数量是可以不同的，例如三角形单元的面积坐标和四面体单元的体积坐标，因为面积坐标和体积坐标都不是完全独立的，其分量的和都是1。在有限元中，通过变量替换，使得Jacobi矩阵成为方阵，是为了求其逆矩阵的需要。只是在有限元中是这样用的而已。不过Jacobi矩阵J的行列式|J|何时为常数？我一直弄不清楚。书上说二维情形

17、下的矩形和平行四边形单元的|J|是常数，三维情形下的正六面体或平行六面体单元的|J|不是常数，请问从理论上怎么理解呢？最好不是数学推导的结果。回楼上的兄弟：雅克比矩阵的行列式|J|为常数，就是里头的每个元素都为常数，以二维等参单元的自然坐标和平面坐标为例，J(1,1)=Ni'Z(Z，n)Xi，只要插值函数Ni是Z的一次函数，则不管任何点(Z，n)上的J(1,1)都是常数。所以插值函数只要是单维线性的，雅克比矩阵的行列式就是常数。ljz07022006-7-2311:32雅可比矩阵的作用是什么?如题,请教雅可比矩阵和弹塑性矩阵的区别,以及在材料的本构关系中为什么有的时候用雅可比矩阵,有时

18、用弹塑性矩阵?谢谢不吝赐教!camus2006-7-2717:22在等参单元中J反映的是整体坐标与局部坐标的转换关系。；ZeQ-odCn弹塑性矩阵反映的是材料的性质与弹模和泊松比有关;在b计算刚度/b矩阵时用到了BTDBJ行列式的局部坐标下的积分b|H2w'K而弹塑性矩阵是在b计算应力/b时用到De怎么理解海森矩阵和雅可比矩阵首先类比一下一维。Jacobian相当于一阶导数，Hessian相当于二阶导数。一维函数的导数的motivation是很明显的。二阶导数的零点就是一阶导数的极值点。对于很多应用，我们不仅关心一阶导数的零点(也就是函数的极值点)，也关心一阶导数的极值点，比如信号处理

19、中，信号的一阶导数的极值点反映信号变化的最剧烈程度。极值点寻求在编程时不方便，不如找二阶导数的零点。Jacobian对于标量函数f:Rn->R1，实际是个向量，这个向量实际上就是函数的梯度gradient。gradient根据Cauchy-Swartz公式，指向的是在某处方向导数取极大值的方向。在二维图像处理中，可用gradient来检测灰度值的边缘。对于向量场F:Rn->Rm,Jacobian的每一行实际都是一个梯度。且有F(X)=F(P)+J(P)(X-P)+O(|X-P|)这个式子的每一行都是一个分量的局部线性化。考虑一个二维的数字图像线性变换(Homography，imag

20、ewarping),以有限差分代替微分，可作类似分析。H:像素(x,y)->像素(u,v)u=u(x,y)v=v(x,y)则其Jacobian为u'(x)u'(y)v'(x)v'(y)反映了局部图像的变形程度。最理想的情况u'(x)=1,v'(y)=1,u'(y)=0,v'(x)=0.说明图像维持原状。由于dudv=|det(Jacobian(x,y)|dxdy(此式的有效性可参考换元法)注：有的书上称det(Jacobian(x,y)为Jacobian.说明面积微元改变的程度由|det(Jacobian(x,y)|决定当|

21、det(Jacobian(x,y)|=1时，说明面积不变，当|det(Jacobian(x,y)|<1时，说明面积压缩，出现了像素丢失现象。当|det(Jacobian(x,y)|>1时，说明面积扩张，需要进行像素插值。另外，由Jacobian矩阵的特征值或奇异值，可作类似说明。可参考Wielandt-Hoffman定理Hessian矩阵定义在标量函数上，对于矢量函数，则成为一个rank3的张量。ls的解释很好Jacobian和Hessian就好比单变量标量函数情况下的一阶导数和二阶导数能很好的概括函数的极值点和单调性Hessian矩阵有一个特例是FisherInformation

22、Matrix,也叫InformationMatrix,是对对数似然函数的Hessian矩阵求期望得到的，衡量了分布中信息量的大小,在统计中非常有用.雅克比矩阵(Jacobiamatrix)以m个n元函数的偏导数为元素的矩阵如果把原来的函数组看作由点到点的一个变换T，则在偏导数都连续的前提之下，u随x的变化由相应的微分方程组来描述。这是一个关于微分的线性方程组，其系数矩阵便是雅克比矩阵(J),因而可写成矩阵形式正交矩阵正交矩阵是实数特殊化的酉矩阵，因此总是正规矩阵。尽管我们在这里只考虑实数矩阵，这个定义可用于其元素来自任何域的矩阵。正交矩阵毕竟是从内积自然引出的，对于复数的矩阵这导致了归一要求。

23、定义定义1如果：AA'=E(E为单位矩阵，A'表示“矩阵A的转置矩阵”。)或A，A=E，则n阶实矩阵A称为正交矩阵，若A为正交阵，则满足以下条件：1) A是正交矩阵2) AA，=E(E为单位矩阵)3) A，是正交矩阵4) A的各行是单位向量且两两正交5) A的各列是单位向量且两两正交6) (Ax,Ay)=(x,y)x,yeR正交矩阵通常用字母Q表示。举例：A=r11r12r13;r21r22r23;r31r32r33则有：rlL2+r12八2+r13八2二r2广2+r22八2+r23八2二r3广2+r32八2+r33八2=1r11*r12+r21*r22+r31*r32=0等性

24、质正交方阵是欧氏空间中标准正交基到标准正交基的过渡矩阵。在矩阵论中，实数正交矩阵是方块矩阵Q,它的转置矩阵是它的逆矩阵：，如果正交矩阵的行列式为+1，则我们称之为特殊正交矩阵:概述要看出与内积的联系，考虑在n维实数内积空间中的关于正交基写出的向量v。v的长度的平方是vv。如果矩阵形式为Qv的线性变换保持了向量长度，则所以有限维线性等距同构，比如旋转、反射和它们的组合，都产生正交矩阵。反过来也成立：正交矩阵蕴涵了正交变换。但是，线性代数包括了在既不是有限维的也不是同样维度的空间之间的正交变换，它们没有等价的正交矩阵。有多种原由使正交矩阵对理论和实践是重要的。nxn正交矩阵形成了一个群，即指示为O

25、(n)的正交群，它和它的子群广泛的用在数学和物理科学中。例如，分子的点群是O(3)的子群。因为浮点版本的正交矩阵有有利的性质，它们是字数值线性代数中很多算法比如QR分解的关键，通过适当的规范化，离散余弦变换(用于MP3压缩)可用正交矩阵表示。例子下面是一些小正交矩阵的例子和可能的解释。恒等变换。旋转16.26°。针对x轴反射。旋转反演(rotoinversion):轴(0,-3/5,4/5)，角度90°。置换坐标轴。基本构造低维度最简单的正交矩阵是1x1矩阵1和-1,它们可分别解释为恒等和实数线针对原点的反射。如下形式的2x2矩阵它的正交性要求满足三个方程在考虑第一个方程时

26、，不丢失一般性而设P=cos0,q=sinG；因此要么t=-q,u=p要么t=q,u=-p。我们可以解释第一种情况为旋转0(0=0是单位矩阵)，第二个解释为针对在角0/2的直线的反射。旋转反射在45°的反射对换x和y；它是置换矩阵，在每列和每行带有一个单一的1(其他都是0):单位矩阵也是置换矩阵。反射是它自己的逆，这蕴涵了反射矩阵是对称的（等于它的转置矩阵）也是正交的。两个旋转矩阵的积是一个旋转矩阵，两个反射矩阵的积也是旋转矩阵。更高维度不管维度，总是可能把正交矩阵按纯旋转与否来分类，但是对于3x3矩阵和更高维度矩阵要比反射复杂多了。例如，和表示通过原点的反演和关于z轴的旋转反演（逆

27、时针旋转90°后针对x-y平面反射，或逆时针旋转270°后对原点反演）。旋转也变得更加复杂；它们不再由一个角来刻画，并可能影响多于一个平面子空间。尽管经常以一个轴和角来描述3x3旋转矩阵，在这个维度旋转轴的存在是偶然的性质而不适用于其他维度。但是，我们有了一般适用的基本建造板块如置换、反射、和旋转。基本变换最基本的置换是换位（transposition），通过交换单位矩阵的两行得到。任何nxn置换矩阵都可以构造为最多n-1次换位的积。构造自非零向量v的Householder反射为这里的分子是对称矩阵，而分母是v的平方量的一个数。这是在垂直于v的超平面上的反射（取负平行于v任

28、何向量分量）。如果v是单位向量，则Q=I-2vv就足够了。Householder反射典型的用于同时置零一列的较低部分。任何nxn正交矩阵都可以构造为最多n次这种反射的积。Givens旋转作用于由两个坐标轴所生成的二维（平面）子空间上，按选定角度旋转。它典型的用来置零一个单一的次对角线元素（subdiagonalentry）。任何nxn的旋转矩阵都可以构造为最多n（n-1）/2次这种旋转的积。在3x3矩阵的情况下，三个这种旋转就足够了；并且通过固定这个序列，我们可以用经常叫做欧拉角的三个角来（尽管不唯一）描述所有3x3旋转矩阵。雅可比旋转有同Givens旋转一样的形式，但是被用做相似变换，选择来

29、置零2x2子矩阵的两个远离对角元素（off-diagonalentry）。性质矩阵性质实数方块矩阵是正交的，当且仅当它的列形成了带有普通欧几里得点积的欧几里得空间R的正交规范基，它为真当且仅当它的行形成R的正交基。假设带有正交（非正交规范）列的矩阵叫正交矩阵可能是诱人的，但是这种矩阵没有特殊价值而没有特殊名字；他们只是MM=D，D是对角矩阵。任何正交矩阵的行列式是+1或-1。这可从关于行列式的如下基本事实得出:反过来不是真的；有+1行列式不保证正交性，即使带有正交列，可由下列反例证实。对于置换矩阵，行列式是+1还是-1匹配置换是偶还是奇的标志，行列式是行的交替函数。比行列式限制更强的是正交矩阵

30、总可以是在复数上可对角化来展示特征值的完全的集合，它们全都必须有（复数）绝对值1。群性质正交矩阵的逆是正交的，两个正交矩阵的积是正交的。事实上，所有nxn正交矩阵的集合满足群的所有公理。它是n(n-1)/2维的紧致李群，叫做正交群并指示为O(n)。行列式为+1的正交矩阵形成了路径连通的子群指标为2的O(n)正规子群，叫做旋转的特殊正交群S0(n)。商群0(n)/S0(n)同构于0(1),带有依据行列式选择+1或-1的投影映射。带有行列式-1的正交矩阵不包括单位矩阵，所以不形成子群而只是陪集；它也是(分离的)连通的。所以每个正交群被分为两个部分；因为投影映射分裂，0(n)是S0(n)与0(1)的

31、半直积。用实用术语说，一个相当的陈述是任何正交矩阵可以通过采用一个旋转矩阵并可能取负它的一列来生成，如我们在2x2矩阵中看到的。如果n是奇数，则半直积实际上是直积，任何正交矩阵可以通过采用一个旋转矩阵并可能取负它的所有列来生成。现在考虑(n+1)x(n+1)右底元素等于1的正交矩阵。最后一列(和最后一行)的余下元素必须是零，而任何两个这种矩阵的积有同样的形式。余下的矩阵是nxn正交矩阵；因此0(n)是0(n+1)(和所有更高维群)的子群。因为Householder正交矩阵形式的基本反射可把任何正交矩阵简约成这种约束形式，一系列的这种反射可以把任何正交矩阵变回单位矩阵；因此正交群是反射群。最后一

32、列可以被固定为任何单位向量，并且每种选择给出不同的0(n)在0(n+1)中的复本；以这种方式0(n+1)是在单位球S与纤维0(n)上的丛。类似的，S0(n)是S0(n+1)的子群；任何特定正交矩阵可以使用类似过程通过Givens平面旋转来生成。丛结构持续：S0(n)S0(n+1)-S。一个单一旋转可以在最后一列的第一行生成一个零，而n-1次旋转序列将置零nxn旋转矩阵的除了最后一列的最后一行的所有元素。因为平面是固定的，每次旋转只有一个自由度，就是它的角度。通过归纳，S0(n)因此有自由度，0(n)也是。置换矩阵简单一些；它们不形成李群，只是一个有限群，n!次对称群Sn。通过同类的讨论，Sn是

33、Sn+1的子群。偶置换生成行列式+1的置换矩阵的子群，n!/2次交错群。规范形式更广泛的说，任何正交矩阵的效果分离到在正交二维空间上的独立动作。就是说，如果Q是狭义正交的，则你可以找到(旋转)改变基的一个正交矩阵P，把Q带回到分块对角形式：(n偶数)，(n奇数)。这里的矩阵R1,.,Rk是2x2旋转矩阵，而余下的元素是零。作为例外，一个旋转块可以是对角的，±1。因此如果需要的话取负一列，并注意2x2反射可对角化为+1和-1，任何正交矩阵可变为如下形式,矩阵R1,，Rk给出位于复平面中单位圆上的特征值的共轭对；所以这个分解复合确定所有带有绝对值1的特征值。如果n是奇数，至少有一个实数特

34、征值+1或-1；对于3x3旋转，关联着+1的特征向量是旋转轴。数值线性代数利益数值分析自然的利用了正交矩阵的很多数值线性代数的性质。例如，经常需要计算空间的正交基，或基的正交变更；二者都采用了正交矩阵的形式。有行列式±1和所有模为1的特征值是对数值稳定性非常有利的。一个蕴涵是条件数为1(这是极小的)，所以在乘以正交矩阵的时候错误不放大。很多算法为此使用正交矩阵如Householder反射和Givens旋转。有帮助的不只是正交矩阵是可逆的，还有它的逆矩阵本质上是免花费的，只需要对换索引(下标)。置换是很多算法成功的根本，包括有局部定支点(partialpivoting)的运算繁重的高斯消去法(这里的置换用来定支点)。但是它们很少明显作为矩阵出现；它们的特殊形式允许更有限的表示，比如n个索引的列表。同样的，使用Householder和Givens矩阵的算法典型的使用特殊方法的乘法和存储。例如，Givens旋转只影响它所乘的矩阵的两行，替代完全的n次的矩阵乘法为更有效的n次运算。在使用这些反射和旋转向矩阵介入零的时候，腾出的空间足够存储充足的数据来重生成这个变换。分解一些重要的矩阵分解(Golub&VanLoan