特征值与特征向量
Eigen Values & Eigen Vectors
Part III:如何求解特征向量与特征值
The Key Equation
对于一般矩阵A,如何找到他的特征值与特征向量?
Step I: Find λ first!
首先,我们有方程:
但这里有两个未知数,因此我们把上面的方程改写一下:
这个齐次方程的解就是矩阵(A-I)的零空间,抛开平凡解全0向量不说。要想让矩阵的零空间存在非零向量,则矩阵的A必为奇异矩阵,即不可逆矩阵。同时,结合之前学到的行列式的概念,若一个矩阵是奇异矩阵,则矩阵的行列式为0。这样一来,我们就不用考虑未知数x,也就是特征向量,先求未知数,也就是特征值。如下:
这个方程是一个非常重要的方程(Key equation),叫特征值方程(Eigen-value equation)或者叫特征方程(characteristic equation)。
Step II: Substitue λ and solve equations!
求解完特征值方程后会得到n个(可能会有相同的),把这些代入到经过改写后的方程组中,求解齐次方程组,或者说是求解零空间,得到相应的特征向量。
举个例子来说明上述求解过程:
对于上面提到过的置换矩阵A,第一步,我们先求det(A-λI)=0:
得到:
由此得到两个特征根:
把第一个特征值λ=1代入改写后的方程:
其中(A-I)为:
对该矩阵消元得到矩阵U,找到主元列和自由列,其中与自由列位置所对应的向量x中的元素就是自由变量。设自由变量为一个任意数,求出相应的特解。
消元矩阵U的秩为1,存在一个特解。第一列为主元列,第二列为自由列。因此,对应的自由变量为x2。令x2=1, 则x1=1得到一个特解,也就是我们要找的特征向量x:
相应的代入第二个特征值λ=-1:
令自由变量x2=1,则x1=-1:
这里我们顺便基于矩阵的行列式和迹来验证一下两个特征根:
特例 I:多个重复/相同的特征值会引发特征向量的短缺
令矩阵A为:
先求特征值:
这会得到两个重复的特征值:
带入其中一个特征值,并求特征向量:
令自由变量x1=1,得到:
现在问题来了,如果你代入第二个特征向量(还是3),并重复上述操作,你是无法找到一个与第一个特征向量线性无关的另一个特征向量了。最多你只是令自由变量x1为一个其他数,然后得到一个特征向量,但最终你得到的这个特征向量与之前得到的那个是线性相关的。
这个例子说明了一个矩阵的特征向量的线性相关性即取决于矩阵的维数也取决于是否存在重复的特征值。对于一个 2x2 的方阵,如果它有两个不同的特征值,那么就能得到两个线性无关的特征向量。但如果是两个重复的特征值,就只能得到一个特征向量或两个线性相关的特征向量。
特例 II:特征根为复数的情况
已知旋转矩阵Q,求Q的特征向量和特征值:
这个旋转矩阵Q的作用是把任意向量旋转90°,这里我们会遇到问题。根据特征向量的意义,当矩阵Q作用于特征向量x后,向量Qx的方向应当与x相同。但这个要求明显与旋转矩阵Q的作用向左。换句话说,我们需要找出一个向量,他旋转90°仍在原处。
其次,根据特征值的性质,我们同样也会发现类似的"Bug"。首先,这是一个2x2矩阵,因此有两个特征值。其次,我们还知道矩阵的迹等于两个特征值的和,矩阵的行列式等于两个特征值的乘积:
从第一个等式上说必定是一正一负,而如果这样的话他们的乘积必定是负数,这和第二个等式相违背。
现在,我们按照一般的方法去求解:
最终得到两个复根:
这两个复数的特征值正好满足和为0,且乘积为1。
带入第一个特征根i:
令自由变量x2=1,得到x1=1/-i=i:
带入第二个特征根-i:
令自由变量x2=1,得到x1=1/i=-i:
(全文完)
作者 --- 松下J27
参考文献(鸣谢):
1,Introduction to Linear Algebra,Fifth Edition - Gilbert Strang
2,麻省理工Gilbert Strang教授线代大师-线性代数(全)_哔哩哔哩_bilibili
(配图与本文无关) 文章来源:https://www.toymoban.com/news/detail-839702.html
版权声明:文中的部分图片,文字或者其他素材,可能来自很多不同的网站和说明,在此没法一一列出,如有侵权,请告知,立即删除。欢迎大家转载,但是,如果有人引用或者COPY我的文章,必须在你的文章中注明你所使用的图片或者文字来自于我的文章,否则,侵权必究。 ----松下J27文章来源地址https://www.toymoban.com/news/detail-839702.html
到了这里,关于线性代数 --- 特征值与特征向量(下)的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!