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

相关推荐
不解风水17 小时前
【教程笔记】KalmanFilter
笔记·学习·算法·矩阵·ekf
CreasyChan1 天前
unity矩阵与变换 - “空间转换的魔术”
unity·矩阵·c#·游戏引擎
Leweslyh1 天前
线性时不变系统传递函数矩阵的状态空间实现理论及其多重性机理研究
线性代数·矩阵
曹文杰15190301121 天前
2025 年大模型背景下应用统计本科 计算机方向 培养方案
python·线性代数·机器学习·学习方法
iAkuya2 天前
(leetcode)力扣100 19螺旋矩阵(方向数组/边界把控)
算法·leetcode·矩阵
闻缺陷则喜何志丹2 天前
【计算几何】仿射变换与齐次矩阵
c++·数学·算法·矩阵·计算几何
闻缺陷则喜何志丹3 天前
【计算几何 线性代数】仿射矩阵的秩及行列式
c++·线性代数·数学·矩阵·计算几何·行列式·仿射矩阵得秩
iAkuya3 天前
(leetcode)力扣100 18矩阵置零(哈希)
leetcode·矩阵·哈希算法
点云侠3 天前
粒子群优化算法求解三维变换矩阵的数学推导
线性代数·算法·矩阵
c#上位机3 天前
halcon计算仿射变换矩阵的逆矩阵
计算机视觉·矩阵·c#