2阶实对称矩阵特征值和特征向量的简单求解方法

  1. 2阶实对称矩阵特性

定理:2阶实对称矩阵H的特征值是实数

H=[a,b;b,c]    

    a,b,c是实数,λ 是特征值
A=[a-λ,b;b,c-λ]    

    特征值求解方法为:(a- λ )(c- λ) - b2  = 0

    求解方程得到两个根为:
λ=(a+c)±(a+c)2-4(ac-b2)2  

                   (a+c)2-4ac-b2=a-c2+4b2≥0

    所以,在a、b、c为实数时,特征值也是实数。

       

2、特征向量

根据特征值和特征向量的定义:HX=λX,(H-λE)X = 0;因此方程若有解,则

det(H-λE)=0;


A=[a-λ,b;b,c-λ]    

则有-b/(a-λ) = (c-λ)/b,   线性齐次方程组AX=0有非零解,其中之一解向量 [1,-b/(a-λ)],归一化后得到标准解。

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
LAPACK是一种线性代数库,用于解决各种数值线性代数问题,包括矩阵特征值特征向量求解。下面我们来介绍一下LAPACK库中求解矩阵特征值特征向量的函数。 1. DSYEV函数 DSYEV函数用于求解实对称矩阵特征值特征向量。该函数的原型如下: ``` void dsyev(char jobz, char uplo, int n, double* a, int lda, double* w, double* work, int lwork, int* info); ``` 参数说明: - jobz:指定计算特征值还是特征向量,取值为'N'(只计算特征值)或'V'(同时计算特征值特征向量)。 - uplo:指定矩阵的上三角部分还是下三角部分存储在数组a中,取值为'U'(上三角部分)或'L'(下三角部分)。 - n:矩阵的维数。 - a:存储矩阵的一维数组。 - lda:指定a数组中每个列向量的存储长度(通常为n)。 - w:存储特征值的一维数组。 - work:工作空间数组。 - lwork:指定work数组的长度(通常为3n)。 - info:返回求解结果,取值为0表示成功,其他值表示出错。 2. ZGEEV函数 ZGEEV函数用于求解矩阵特征值特征向量。该函数的原型如下: ``` void zgeev(char jobvl, char jobvr, int n, std::complex<double>* a, int lda, std::complex<double>* w, std::complex<double>* vl, int ldvl, std::complex<double>* vr, int ldvr, std::complex<double>* work, int lwork, double* rwork, int* info); ``` 参数说明: - jobvl:指定是否计算左特征向量,取值为'N'(不计算)或'V'(计算)。 - jobvr:指定是否计算右特征向量,取值为'N'或'V'。 - n:矩阵的维数。 - a:存储矩阵的一维数组。 - lda:指定a数组中每个列向量的存储长度(通常为n)。 - w:存储特征值的一维数组。 - vl:存储左特征向量的一维数组。 - ldvl:指定vl数组中每个列向量的存储长度(通常为n)。 - vr:存储右特征向量的一维数组。 - ldvr:指定vr数组中每个列向量的存储长度(通常为n)。 - work:工作空间数组。 - lwork:指定work数组的长度(通常为2n)。 - rwork:实数数组,长度为2n(用于存储中间计算结果)。 - info:返回求解结果,取值为0表示成功,其他值表示出错。 以上就是LAPACK库中求解矩阵特征值特征向量的函数介绍。需要注意的是,在调用这些函数之前,需要先将矩阵按列存储方式存储在一维数组中,并传入一些参数,如矩阵的维数、存储方式等。具体的参数可以参考LAPACK库的文档。
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值