跳转到内容

欧氏空间、特征向量和相似对角化?

有些东西太简单了就不放进来浪费时间了

如果要把一组化为正交基,那最常用的方法叫做格拉姆-施密特正交化方法(Gram-Schmidt process) ,就是计算量可能有一点大

下面的(α2,β1)(\alpha_2,\beta_1)意为α2\alpha_2β1\beta_1的点积

(α1,α2,α3)(\alpha_1,\alpha_2,\alpha_3)是原本的一组不正交的向量

{β1=α1β2=α2(α2,β1)(β1,β1)β1β3=α3(α3,β1)(β1,β1)β1(α3,β2)(β2,β2)β2\left\{ \begin{aligned} &\beta_1 = \alpha_1 \\ &\beta_2 = \alpha_2 - \frac{(\alpha_2,\beta_1)}{(\beta_1,\beta_1)}\beta_1 \\ &\beta_3 = \alpha_3 - \frac{(\alpha_3,\beta_1)}{(\beta_1,\beta_1)}\beta_1 - \frac{(\alpha_3,\beta_2)}{(\beta_2,\beta_2)}\beta_2 \end{aligned} \right.

得到的(β1,β2,β3)(\beta_1,\beta_2,\beta_3)就是一组正交向量,两者作为是等价的

更多的向量以此类推就好

另外,(β1,β1)(\beta_1,\beta_1)算的是β1\beta_1平方!!!不是它的长度!!!

这有一张相关的图,来自维基百科。其中的vvV2V^2上投影,构造出V3V^3上的正交基β\beta

格拉姆-施密特正交化的图

还有一个概念叫做正交矩阵,一个满足ATA=AAT=EA^TA=AA^T=E的矩阵AA就是正交矩阵,这个矩阵对应的线性变换只旋转而不缩放,于是向量的前后长度是不会变的

可以发现AT=A1A^T=A^{-1}

A1A^{-1}AtA^t也都是正交矩阵,而且A=±1\vert A\vert=±1,如果BB也是同阶的正交矩阵,那ABAB也是正交矩阵

特征向量就是经过一个线性变换后指向方向没有改变的向量,特征值就是经过这个线性变换后向量长度缩放的值,即

Av=λvAv=\lambda v 其中vv为特征向量,λ\lambda为特征值,一组特征向量和特征值也可以表示一个线性变换

一个nn阶的方阵会有nn个对应的特征值,尽管可能重复

由定义可知(Aλ)v=0(A-\lambda)v=0,即AλE=0\left \vert A-\lambda E \right\vert = 0

所以只要解AλE=0\vert A-\lambda E\vert=0这个方程就好(解λEA=0\vert\lambda E-A\vert=0也是一样的

得出的几个λ\lambda就是特征值,每个特征值对应的解的数量叫做它的代数重数,对应的线性无关的向量的数量叫做几何重数

几何重数亦可称为是这些向量张成的特征子空间维数

一个矩阵的特征值正负数量分别被称为它的正惯性指数负惯性指数

既然已经算出了λ\lambda的值,接下来就代入原式(Aλ)v=0(A-\lambda)v=0算就好了(当然也有可能是(λA)v=0(\lambda-A)v=0

专家研究发现i=1naii=i=1nλi\displaystyle \sum_{i=1}^n a_{ii} = \displaystyle \sum_{i=1}^n \lambda_i,也就是主对角线上面所有元素的和等于所有特征值的和重复的特征值也要重复算!)然后这个和还被起了个新名字叫AA的迹记作tr(A)tr(A),看起来是和秩属于同一批产物

以及A=i=1nλi\vert A \vert =\displaystyle \prod_{i=1}^n \lambda_i,用人话讲就是特征值的积等于矩阵行列式的值

如果说λ\lambdann阶方阵AA的特征值

AmA^m的特征值就是λm\lambda^mkAkA的特征值就是kλk\lambda,而特征向量不变

如果AA还是可逆的,那A1A^{-1}的特征值是1λ\displaystyle\frac{1}{\lambda},而Aλ\displaystyle\frac{A}{\lambda}就是AA^*的特征值,而特征向量还是不变

转置不改变特征值,但是特征向量不确定

==nn阶方阵AA可相似对角化的充要条件是AA每个特征值代数重数等于它的几何重数==

那啥是相似对角化?

相似对角化顾名思义就是把矩阵化为一个相似的对角阵,此时由于对角阵的特性,算其幂次时会相当方便

不过可能是由于我们的水平太拉了,老师只教了怎么把实对称矩阵化为对角阵

对一个实对称矩阵而言,必定与一个对角阵相似,所以都是有结果的

==对于实对称矩阵,不同特征值算出来的特征向量是正交的==

目标很简单,如果有个33阶实对称矩阵AA,把它相似对角化就是求一个正交矩阵QQ实现

QTAQ=Q1AQ=Λ=[λ1λ2λ3]Q^TAQ=Q^{-1}AQ=\Lambda= \left[ \begin{matrix} \lambda_1 & & \\ & \lambda_2 & \\ & & \lambda_3 \end{matrix}\right]

QQ又由三个单位向量组成,即Q=(η1,η2,η3)Q=(\eta_1,\eta_2,\eta_3)λ\lambdaη\eta 的位置必须要相对应,不然不对

接下来就开始算吧

首先用上面的方法算出特征值和特征向量,然后看这几个特征向量是否正交,不正交就统统正交化

最后单位化(长度变为1)就得到(η1,η2,η3)(\eta_1,\eta_2,\eta_3)