APP下载

矩阵剖析与PCA实例

2020-08-19万子峰蔡泽雄帅晨珊陈梦琪

黑龙江科学 2020年16期
关键词:对角线协方差特征向量

万子峰,蔡泽雄,帅晨珊,陈梦琪

(南昌航空大学 a.信息工程学院; b.国际教育学院,南昌 330001)

1 矩阵变换

f(0,1)=(3,5,9)=3(1,0,0)+5(0,1,0)+9(0,0,1)

f(1,0)=(1,2,7)=1(1,0,0)+2(0,1,0)+7(0,0,1)

简而言之,B在f(A)上的投影为T(映射)。故有定义“矩阵表达的就是对线性运动的一种描述[1]”。

显然,旧基能线性表达新基且矩阵T为过渡矩阵。

同时过渡矩阵T为满秩矩阵,是非奇异矩阵,由“非奇异矩阵行列式不为0”的性质可推T可逆得证。

图1 向量变换立体图Fig.1 Stereogram of vector transformation

注意:原向量只能对应新向量在新坐标系中的所占权重,因为基之间并不能保证相互正交且长度为单位长度。

2 相似

教科书上对“相似”的定义:存在可逆矩阵p,使B=p-1Ap成立,那么说明A相似于B[2]。证明这个定义需要之前证明过的两个推理:

3.若第1条推理成立,则存在旧基坐标x与新基坐标y的转换关系y=T-1x,设有关系图表S,

图2 关系图表SFig.2 Relationship Chart S

通过图2,有:

①λ=T-1u(推理2)。

②λ=M2η(第2节的左乘推理)。

联立①②,得到③M2η=T-1u。

④由η=T-1j,两侧左乘T,得到Tη=j。

⑤u=M1j(左乘推理)。

联立③④,得到⑥M2T-1j=T-1u。

联立⑤⑥,得到⑦M2T-1j=T-1M1j。

⑦化简,得到结果M2=T-1M1T。

证毕。

所以相似矩阵也是一种对同一个线性运动在不同基的矩阵描述。纵使对矩阵进行了各种线性变化后,只要保证结果矩阵与原矩阵维持相似状态,则运动的本质就不会发生改变。

3 特征值与特征向量

AAT=U∑VT·V∑UT=U∑2U-1奇异值分解(SVD)定义:一个矩阵B可以分解成B=U∑VT,其中U的列向量两两正交且模为1,V列向量两两正交且模为1,即UT=U-1,VT=V-1,∑矩阵的主对角线为奇异值,其他元素为0。特征值是相对于方阵来说的,而奇异值分解则解除了这一限制,它适用于任意的矩阵。特征值与奇异值分解的关系:对任意矩阵A作svd,再对AAT作svd,AAT=U′∑V′T。又有U′=U,U-1=UT=V′T∴U′=V′=U,∑′=∑2。因此,对AAT作SVD,U′为其特征向量(为A中的U)。同理,对ATA作SVD,U′为其特征向量(为A中的V),∑′为特征值对角矩阵。

4 协方差矩阵

通常处理的数据都是多维的,“协方差”就是这样用来度量2个变量的关系,如果维数多于2,就用协方差矩阵。协方差定义:

举例,假设有1组数据,用矩阵表示为:

根据协方差定义:

cov(x,y)=E[(1-E(x))(2-E(y))+(35-E(x))(25-E(y))+…]

cov(x,z)=E[(1-E(x))(3-E(z))+(35-E(x))(55-E(z))+…]

……

通过数据初始化,使E(x)=E(y)=E(z)=0,则X的协方差可表示为:

我们发现,对角线的元素就是各个变量的方差,即:

5 PCA过程

猜你喜欢

对角线协方差特征向量
二年制职教本科线性代数课程的几何化教学设计——以特征值和特征向量为例
克罗内克积的特征向量
一类三阶矩阵特征向量的特殊求法
高效秩-μ更新自动协方差矩阵自适应演化策略
基于子集重采样的高维资产组合的构建
用于检验散斑协方差矩阵估计性能的白化度评价方法
二维随机变量边缘分布函数的教学探索
边、角、对角线与平行四边形的关系
看四边形对角线的“气质”
数学题