关于矩阵类中利用雅克比法求出的特征值和求矩阵特征向量量的对应关系问题.

热门搜索:
您现在的位置是: &
> 雅可比法计算特征向量程序雅可比过关法求实对称矩阵特征值与特征向量
雅可比法计算特征向量程序雅可比过关法求实对称矩阵特征值与特征向量
资 源 简 介
雅可比法计算特征向量程序雅可比过关法求实对称矩阵特征值与特征向量
VIP 专区(每个包含40-100个资源包)
您 可 能 感 兴 趣 的
相 关 代 码
相 关 资 源
该 用 户 还 上 传
月度VIP会员
月度VIP会员
月度VIP会员
月度VIP会员
月度VIP会员
月度VIP会员
月度VIP会员
月度VIP会员
月度VIP会员
月度VIP会员请问如何用雅克比法求解矩阵特征值和特征向量
xiadd414確
雅可比方法的基本思想是通过一系列的由平面旋转矩阵构成的正交变换将实对称矩阵逐步化为对角阵,从而得到 的全部特征值及其相应的特征向量.首先引进 中的平面旋转变换.变换(7)记为 ,其中 (8)则称 为 中 平面内的一个平面旋转变换,称为 平面内的平面旋转矩阵.容易证明 具有如下简单性质:① 为正交矩阵.② 的主对角线元素中除第 个与第 个元素为 外,其它元素均为1;非对角线元素中除第 行第 列元素为 ,第 行第 列元素为 外,其它元素均为零.③ 只改变 的第 行与第 行元素,只改变 的第 列与第 列元素,所以 只改变 的第 行、第 行、第 列、第 列元素.设 为 阶实对称矩阵,为一对非对角线元素.令则 为实对称矩阵,且 与 有相同的特征值.通过直接计算知(9)当取 满足关系式(10)时,,且(11)由于在正交相似变换下,矩阵元素的平方和不变,所以若用 表示矩阵 的对角线元素平方和,用 表示 的非对角线元素平方和,则由(11)式得(12) 这说明用 对 作正交相似变换化为 后,的对角线元素平方和比 的对角线元素平方和增加了 ,的非对角线元素平方和比 的非对角线元素平方和减少了 ,且将事先选定的非对角线元素消去了(即 ).因此,只要我们逐次地用这种变换,就可以使得矩阵 的非对角线元素平方和趋于零,也即使得矩阵 逐步化为对角阵.这里需要说明一点:并不是对矩阵 的每一对非对角线非零元素进行一次这样的变换就能得到对角阵.因为在用变换消去 的时候,只有第 行、第 行、第 列、第 列元素在变化,如果 或 为零,经变换后又往往不是零了.雅可比方法就是逐步对矩阵 进行正交相似变换,消去非对角线上的非零元素,直到将 的非对角线元素化为接近于零为止,从而求得 的全部特征值,把逐次的正交相似变换矩阵乘起来,便是所要求的特征向量.雅可比方法的计算步骤归纳如下:第一步 在矩阵 的非对角线元素中选取一个非零元素 .一般说来,取绝对值最大的非对角线元素;第二步 由公式 求出 ,从而得平面旋转矩阵 ;第三步 ,的元素由公式(9)计算.第四步 以 代替 ,重复第一、二、三步求出 及 ,继续重复这一过程,直到 的非对角线元素全化为充分小(即小于允许误差)时为止.第五步 的对角线元素为 的全部特征值的近似值,的第j列为对应于特征值 ( 为 的对角线上第j个元素)的特征向量.
为您推荐:
其他类似问题
用雅可比法(Jacobi)计算对称矩阵的特征值和特征向量 ' 参数: n - Integer型变量,对称矩阵的阶数。 ' dblA - Double型二维数组,体积为n x n。我也许说的不清楚,其实我想弄明白如何用雅克比法求解实对称矩阵的特征值,具体的步骤是什么?谢谢啊...
扫描下载二维码用fortran语言编写: 用雅克比方法求矩阵特征值和特征向量_百度知道
用fortran语言编写: 用雅克比方法求矩阵特征值和特征向量
用雅克比方法求矩阵特征值和特征向量
我有更好的答案
SUBROUTINE JCB(N,A,S,EPS)
THIS SUBROUTINE COMPUTS EIGENVALUES AND standard EIGENVECTORS OF A
DIMENSION A(N,N),S(N,N)
DO 30 I=1,N
DO 30 J=1,I
IF(I-J) 20,10,20
DO 40 I=2,N
DO 40 J=1,I1
G=G+2.*A(I,J)*A(I,J)
S1=SQRT(G)
S2=EPS/FLOAT(N)*S1
S3=S3/FLOAT(N)
DO 130 IQ=2,N
DO 130 IP=1,IQ1
IF(ABS(A(IP,IQ)).LT.S3) GOTO 130
V1=A(IP,IP)
V2=A(IP,IQ)
V3=A(IQ,IQ)
U=0.5*(V1-V3)
其他类似问题
为您推荐:
您可能关注的推广
矩阵特征值的相关知识
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁这是个机器人猖狂的时代,请输一下验证码,证明咱是正常人~Jacobi矩阵 特征值 雅克比网上论文出下一个正特征值的雅克比矩阵,这个是啥意思?难道还分正负吗?
吾乃萝莉控247
Jacobi矩阵一般是坐标变换时用的——这个坐标系的坐标依次对那个坐标系的坐标求偏导数即得Jacobi矩阵.积分时从直角坐标变到球坐标、极坐标肯定要在积分号里乘上Jacobi矩阵的行列式.若其特征值全为正,则此Jacobi矩阵正定,即它的任何二次型都大于零.并且如果是对称Jacobi矩阵,还可以实相合于一个单位阵.是比较强的条件.矩阵A的特征值是这个关于x的方程的所有根det(A-xI)=0
为您推荐:
其他类似问题
扫描下载二维码

我要回帖

更多关于 求矩阵特征向量 的文章

 

随机推荐