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

相关推荐
我是大聪明.13 小时前
Attention机制的数学本质:从Softmax到FlashAttention的演进
线性代数·矩阵
我是大聪明.2 天前
大模型Tokenizer原理:BPE、WordPiece与子词编码的核心机制深度解析
人工智能·线性代数·算法·机器学习·矩阵
xin_nai2 天前
LeetCode热题100(Java)(6)矩阵
java·leetcode·矩阵
萌新小码农‍3 天前
人工智能线性代数基础
人工智能·线性代数·机器学习
生信研究猿3 天前
#P4538.第2题-基于混淆矩阵,推导分类模型的核心评估指标
线性代数·矩阵
小白小宋3 天前
【PUSCH第三期】5G NR QC-LDPC编码深度解析:从协议校验矩阵构造到MATLAB完整实现
5g·matlab·矩阵
啦啦啦_99993 天前
1. 线性回归之 向量&矩阵
算法·矩阵·线性回归
star learning white3 天前
线性代数3
人工智能·线性代数·机器学习
爱吃巧克力的程序媛3 天前
计算机图形学---如何理解模型矩阵、视图矩阵、投影矩阵
数码相机·线性代数·矩阵
做cv的小昊4 天前
【TJU】研究生应用统计学课程笔记(5)——第二章 参数估计(2.3 C-R不等式)
c语言·笔记·线性代数·机器学习·数学建模·r语言·概率论