有些东西太简单了就不放进来浪费时间了
如果要把一组基化为正交基,那最常用的方法叫做格拉姆-施密特正交化方法(Gram-Schmidt process) ,就是计算量可能有一点大
下面的(α2,β1)意为α2和β1的点积
(α1,α2,α3)是原本的一组不正交的向量
⎩⎨⎧β1=α1β2=α2−(β1,β1)(α2,β1)β1β3=α3−(β1,β1)(α3,β1)β1−(β2,β2)(α3,β2)β2
得到的(β1,β2,β3)就是一组正交向量,两者作为基是等价的
更多的向量以此类推就好
另外,(β1,β1)算的是β1的平方!!!不是它的长度!!!
这有一张相关的图,来自维基百科。其中的v在V2上投影,构造出V3上的正交基β

还有一个概念叫做正交矩阵,一个满足ATA=AAT=E的矩阵A就是正交矩阵,这个矩阵对应的线性变换只旋转而不缩放,于是向量的前后长度是不会变的
可以发现AT=A−1
A−1和At也都是正交矩阵,而且∣A∣=±1,如果B也是同阶的正交矩阵,那AB也是正交矩阵
特征向量就是经过一个线性变换后指向方向没有改变的向量,特征值就是经过这个线性变换后向量长度缩放的值,即
Av=λv
其中v为特征向量,λ为特征值,一组特征向量和特征值也可以表示一个线性变换
一个n阶的方阵会有n个对应的特征值,尽管可能重复
由定义可知(A−λ)v=0,即∣A−λE∣=0
所以只要解∣A−λE∣=0这个方程就好(解∣λE−A∣=0也是一样的
得出的几个λ就是特征值,每个特征值对应的解的数量叫做它的代数重数,对应的线性无关的向量的数量叫做几何重数
几何重数亦可称为是这些向量张成的特征子空间的维数
一个矩阵的特征值正负数量分别被称为它的正惯性指数和负惯性指数
既然已经算出了λ的值,接下来就代入原式(A−λ)v=0算就好了(当然也有可能是(λ−A)v=0
专家研究发现i=1∑naii=i=1∑nλi,也就是主对角线上面所有元素的和等于所有特征值的和(重复的特征值也要重复算!)然后这个和还被起了个新名字叫迹,A的迹记作tr(A),看起来是和秩属于同一批产物
以及∣A∣=i=1∏nλi,用人话讲就是特征值的积等于矩阵行列式的值
如果说λ是n阶方阵A的特征值
那Am的特征值就是λm,kA的特征值就是kλ,而特征向量不变
如果A还是可逆的,那A−1的特征值是λ1,而λA就是A∗的特征值,而特征向量还是不变
转置不改变特征值,但是特征向量不确定
==n阶方阵A可相似对角化的充要条件是A的每个特征值的代数重数等于它的几何重数==
那啥是相似对角化?
相似对角化顾名思义就是把矩阵化为一个相似的对角阵,此时由于对角阵的特性,算其幂次时会相当方便
不过可能是由于我们的水平太拉了,老师只教了怎么把实对称矩阵化为对角阵
对一个实对称矩阵而言,必定与一个对角阵相似,所以都是有结果的
==对于实对称矩阵,不同特征值算出来的特征向量是正交的==
目标很简单,如果有个3阶实对称矩阵A,把它相似对角化就是求一个正交矩阵Q实现
QTAQ=Q−1AQ=Λ=λ1λ2λ3
Q又由三个单位向量组成,即Q=(η1,η2,η3),λ 和 η 的位置必须要相对应,不然不对
接下来就开始算吧
首先用上面的方法算出特征值和特征向量,然后看这几个特征向量是否正交,不正交就统统正交化
最后单位化(长度变为1)就得到(η1,η2,η3)了