第4章_微分运动和雅可比矩阵
- 格式:ppt
- 大小:2.89 MB
- 文档页数:70


雅可比矩阵(Jacobi方法)Jacobi 方法Jacobi方法是求对称矩阵的全部特征值以及相应的特征向量的一种方法,它是基于以下两个结论1) 任何实对称矩阵A可以通过正交相似变换成对角型,即存在正交矩阵Q,使得Q T AQ = diag(λ1 ,λ2,…,λn) (3.1)其中λi(i=1,2,…,n)是A的特征值,Q中各列为相应的特征向量。
2) 在正交相似变换下,矩阵元素的平方和不变。
即设A=(aij )n×n,Q交矩阵,记B=Q T AQ=(bij )n×n, 则Jacobi方法的基本思想是通过一次正交变换,将A中的一对非零的非对角化成零并且使得非对角元素的平方和减小。
反复进行上述过程,使变换后的矩阵的非对角元素的平方和趋于零,从而使该矩阵近似为对角矩阵,得到全部特征值和特征向量。
1 矩阵的旋转变换设A为n阶实对称矩阵,考虑矩阵易见 Vij(φ)是正交矩阵, 记注意到B=VijA的第i,j行元素以及的第i,j列元素为可得≠0,取φ使得则有如果aij对A(1)重复上述的过程,可得A(2) ,这样继续下去, 得到一个矩阵序列{A(k) }。
可以证明,虽然这种变换不一定能使矩阵中非对角元素零元素的个数单调增加,但可以保证非对角元素的平方和递减,我们以A与A(1)为例进行讨论。
设由式(3.4)可得这表明,在上述旋转变换下,非对角元素的平方和严格单调递减,因而由(3.2)可知,对角元素的平方和单调增加。
2. Jacobi方法通过一系列旋转变换将A变成A(k+1) ,求得A的全部特征值与特征向量的方法称为Jacobi方法。
计算过程如下1)令k=0, A(k) =A2) 求整数i,j, 使得3) 计算旋转矩阵4) 计算A(k+1)5) 计算6) 若E(A(k+1))<ε, 则为特征值,Q T = (V(0) V(1)…V(k+1))T的各列为相应的特征向量;否则,k+1=>k返回2,重复上述过程。
雅可比(Jacobian)矩阵2008-12-05 11:29在向量微积分中,雅可比矩阵是一阶偏导数以一定方式排列成的矩阵,其行列式称为雅可比行列式。
还有,在代数几何中,代数曲线的雅可比量表示雅可比簇:伴随该曲线的一个群簇,曲线可以嵌入其中。
它们全部都以数学家卡尔·雅可比命名;雅可比矩阵雅可比矩阵的重要性在于它体现了一个可微方程与给出点的最优线性逼近。
因此,雅可比矩阵类似于多元函数的导数。
假设F:Rn→Rm 是一个从欧式n维空间转换到欧式m维空间的函数。
这个函数由m个实函数组成: y1(x1,...,xn), ..., ym(x1,...,xn). 这些函数的偏导数(如果存在)可以组成一个m行n列的矩阵,这就是所谓的雅可比矩阵:此矩阵表示为:,或者这个矩阵的第i行是由梯度函数的转置y i(i=1,...,m)表示的如果p是Rn中的一点,F在p点可微分,那么在这一点的导数由J F(p)给出(这是求该点导数最简便的方法)。
在此情况下,由F(p)描述的线性算子即接近点p的F的最优线性逼近,x逼近与p例子由球坐标系到直角坐标系的转化由F函数给出:R × [0,π] × [0,2π] → R3此坐标变换的雅可比矩阵是R4的f函数:其雅可比矩阵为:此例子说明雅可比矩阵不一定为方矩阵。
在动态系统中考虑形为x' = F(x)的动态系统,F : R n→ R n。
如果F(x0) = 0,那么x0是一个驻点。
系统接近驻点时的表现通常可以从JF(x0)的特征值来决定。
雅可比行列式如果m = n,那么F是从n维空间到n维空间的函数,且它的雅可比矩阵是一个方块矩阵。
于是我们可以取它的行列式,称为雅可比行列式。
在某个给定点的雅可比行列式提供了F在接近该点时的表现的重要信息。
例如,如果连续可微函数F在p点的雅可比行列式不是零,那么它在该点具有反函数。
这称为反函数定理。
更进一步,如果p点的雅可比行列式是正数,则F在p 点的取向不变;如果是负数,则F的取向相反。
雅可比矩阵(Jacobi方法)Jacobi 方法Jacobi方法是求对称矩阵的全部特征值以及相应的特征向量的一种方法,它是基于以下两个结论1) 任何实对称矩阵A可以通过正交相似变换成对角型,即存在正交矩阵Q,使得Q T AQ = diag(λ1 ,λ2,…,λn) (3.1)其中λi(i=1,2,…,n)是A的特征值,Q中各列为相应的特征向量。
2) 在正交相似变换下,矩阵元素的平方和不变。
即设A=(aij )n×n,Q交矩阵,记B=Q T AQ=(bij )n×n, 则Jacobi方法的基本思想是通过一次正交变换,将A中的一对非零的非对角化成零并且使得非对角元素的平方和减小。
反复进行上述过程,使变换后的矩阵的非对角元素的平方和趋于零,从而使该矩阵近似为对角矩阵,得到全部特征值和特征向量。
1 矩阵的旋转变换设A为n阶实对称矩阵,考虑矩阵易见 Vij(φ)是正交矩阵, 记注意到B=VijA的第i,j行元素以及的第i,j列元素为可得≠0,取φ使得则有如果aij对A(1)重复上述的过程,可得A(2) ,这样继续下去, 得到一个矩阵序列{A(k) }。
可以证明,虽然这种变换不一定能使矩阵中非对角元素零元素的个数单调增加,但可以保证非对角元素的平方和递减,我们以A与A(1)为例进行讨论。
设由式(3.4)可得这表明,在上述旋转变换下,非对角元素的平方和严格单调递减,因而由(3.2)可知,对角元素的平方和单调增加。
2. Jacobi方法通过一系列旋转变换将A变成A(k+1) ,求得A的全部特征值与特征向量的方法称为Jacobi方法。
计算过程如下1)令k=0, A(k) =A2) 求整数i,j, 使得3) 计算旋转矩阵4) 计算A(k+1)5) 计算6) 若E(A(k+1))<ε, 则为特征值,Q T = (V(0) V(1)…V(k+1))T的各列为相应的特征向量;否则,k+1=>k返回2,重复上述过程。