主成分分析(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)对应于两个最大的特征值。

相关推荐
CUC-MenG3 小时前
2025牛客多校第六场 D.漂亮矩阵 K.最大gcd C.栈 L.最小括号串 个人题解
c语言·c++·算法·矩阵
C++、Java和Python的菜鸟13 小时前
第二章 矩阵
人工智能·算法·矩阵
绿炮火15 小时前
【MATLAB】(八)矩阵
人工智能·matlab·矩阵
蜉蝣之翼❉1 天前
图像处理拉普拉斯算子
图像处理·线性代数·矩阵
郝学胜-神的一滴2 天前
能表示旋转的矩阵是一个流形吗?
线性代数·矩阵·图形渲染
ykjhr_3d2 天前
华锐矩阵世界平台与海外客户洽谈合作
线性代数·矩阵
weisian1513 天前
力扣经典算法篇-42-矩阵置零(辅助数组标记法,使用两个标记变量)
算法·leetcode·矩阵
设计师小聂!3 天前
力扣热题100--------240.搜索二维矩阵
java·算法·leetcode·矩阵
盛寒3 天前
二次型 线性代数
线性代数
Ghost-Face4 天前
Ctorch开发日志——矩阵乘法优化及数学原理
矩阵·ctorch