主成分分析(PCA)例题——给定协方差矩阵

向量 x x x的相关矩阵为

R x = [ 0.3 0.1 0.1 0.1 0.3 − 0.1 0.1 − 0.1 0.3 ] {\bm R}_x = \begin{bmatrix} 0.3 & 0.1 & 0.1 \\ 0.1 & 0.3 & -0.1 \\ 0.1 & -0.1 & 0.3 \end{bmatrix} Rx= 0.30.10.10.10.3−0.10.1−0.10.3

计算输入向量的 KL 变换。


解答

R x {\bm R}_x Rx的特征值为 λ 0 = 0.1 \lambda_0 = 0.1 λ0=0.1, λ 1 = λ 2 = 0.4 \lambda_1 = \lambda_2 = 0.4 λ1=λ2=0.4。

既然 R x {\bm R}_x Rx是对称的,可以构建正交特征向量。

在这个例子中,有

u 0 = 1 3 [ 1 − 1 − 1 ] , u 1 = 1 6 [ 2 1 1 ] , u 2 = 1 2 [ 0 1 − 1 ] {\bm u}_0 = \frac{1}{\sqrt{3}} \begin{bmatrix} 1 \\ -1 \\ -1 \end{bmatrix}, \quad {\bm u}_1 = \frac{1}{\sqrt{6}} \begin{bmatrix} 2 \\ 1 \\ 1 \end{bmatrix}, \quad {\bm u}_2 = \frac{1}{\sqrt{2}} \begin{bmatrix} 0 \\ 1 \\ -1 \end{bmatrix} u0=3 1 1−1−1 ,u1=6 1 211 ,u2=2 1 01−1

则 KL 变换为

y ( 0 ) y ( 1 ) y ( 2 ) \] = \[ 2 / 6 1 / 6 1 / 6 0 1 / 2 − 1 / 2 1 / 3 − 1 / 3 − 1 / 3 \] \[ x ( 0 ) x ( 1 ) x ( 2 ) \] \\begin{bmatrix} y(0) \\\\ y(1) \\\\ y(2) \\end{bmatrix} = \\begin{bmatrix} 2/\\sqrt{6} \& 1/\\sqrt{6} \& 1/\\sqrt{6} \\\\ 0 \& 1/\\sqrt{2} \& -1/\\sqrt{2} \\\\ 1/\\sqrt{3} \& -1/\\sqrt{3} \& -1/\\sqrt{3} \\end{bmatrix} \\begin{bmatrix} x(0) \\\\ x(1) \\\\ x(2) \\end{bmatrix} y(0)y(1)y(2) = 2/6 01/3 1/6 1/2 −1/3 1/6 −1/2 −1/3 x(0)x(1)x(2) 其中 y ( 0 ) , y ( 1 ) y(0), y(1) y(0),y(1)对应于两个最大的特征值。

相关推荐
撩妹小狗6 天前
科赫雪花--Python--数学原理--turtle绘图
python·线性代数·几何学
张3蜂6 天前
线性代数与AI的关系
人工智能·线性代数
汉克老师6 天前
GESP2023年12月认证C++二级( 第三部分编程题(2) 小杨的H字矩阵)
c++·算法·矩阵·循环结构·gesp二级·gesp2级
AI科技星6 天前
物理世界的几何建构:论统一场论的本体论革命与概念生成
人工智能·opencv·线性代数·算法·矩阵
没有bug.的程序员6 天前
订单系统重构史诗:从单体巨兽到微服务矩阵的演进、数据一致性内核与分布式事务
java·微服务·矩阵·重构·分布式事务·数据一致性·订单系统
super_lzb6 天前
【线性代数】矩阵第一讲:矩阵与矩阵的运算
线性代数·矩阵·考研数学·矩阵的计算
newbiai6 天前
TikTok矩阵账号引流怎么解决效率低成本高?
python·线性代数·矩阵
逆境不可逃6 天前
【从零入门23种设计模式08】结构型之组合模式(含电商业务场景)
线性代数·算法·设计模式·职场和发展·矩阵·组合模式
菜鸡儿齐7 天前
leetcode-搜索二维矩阵
算法·leetcode·矩阵
炽烈小老头7 天前
【每天学习一点算法 2026/02/24】矩阵置零
学习·算法·矩阵