线性代数——矩阵的逆


线性代数—矩阵的逆


brycezou@163.com

1、计算二阶矩阵的逆


[acbd]1=1adbc[dcba]

2、计算矩阵的逆(余子式法)


A=101021151

  • 矩阵A的余子式(matrix of minors)

    B=215101110215015111111015012111011002=312505212

  • 矩阵A的代数余子式(matrix of cofactors)

    C=B×+++++=312505212

  • 矩阵A的伴随矩阵(the adjoint of matrix)

    A=CT=352101252

  • 矩阵A的行列式

    |A|=+A11B11A12B12+A13B13=5|A|=A21B21+A22B22A23B23=5|A|=+A31B31A32B32+A33B33=5

  • 矩阵A的逆(the inverse of matrix)

    A1=1|A|A=15352101252

3、计算矩阵的逆(行初等变换法)





        每作一次初等行变换,就相当于给增广矩阵的左侧和右侧分别左乘一个初等变换矩阵,当左侧变为单位矩阵时,右侧就变为原矩阵的逆矩阵。

4、解线性方程组


{3a+2b=76a+6b=6     [3626][ab]=[76]

[ab]=[3626]1[76]=[12]

5、子空间窥探


        已知
a⃗ =[36],  b⃗ =[26],  c⃗ =[76]

        问是否存在系数 x,y 使得
a⃗ x+b⃗ y=c⃗ 

        求解可得
[xy]=[12]

        这可以理解为,在向量a和b所张成的子空间中,向量c的坐标为(1,2)

6、判断矩阵是否可逆


A1=1|A|A

        由上述定义可知,当 |A|=0 时,矩阵 A 的逆不存在。下面以二阶矩阵为例来阐述其背后的原因,设
A=[acbd]

        如果 |A|=adbc=0 ,则有
ac=bd       ab=cd

        因次,下列线性方程组无解
[acbd][xy]=[ef]=>{ax+by=ecx+dy=f=>{y=abx+eby=cdx+fd

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库的文档。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值