当前位置: 首页 > java >正文

PCA例题

已知样本集合的协方差矩阵为

C x = 1 10 [ 3 1 1 1 3 − 1 1 − 1 3 ] {\bm C}_x = \frac{1}{10} \begin{bmatrix} 3 & 1 & 1 \\ 1 & 3 & -1 \\ 1 & -1 & 3 \end{bmatrix} Cx=101 311131113

使用PCA方法将样本向量降到二维 。


求解
计算 C x {\bm C}_x Cx的特征值得到:

λ 1 = 0.1 , λ 2 = λ 3 = 0.4 \lambda_1 = 0.1, \quad \lambda_2 = \lambda_3 = 0.4 λ1=0.1,λ2=λ3=0.4

其对应的特征向量分别为:

u 1 = 1 3 [ 1 − 1 − 1 ] , u 2 = 1 6 [ 2 1 1 ] , u 3 = 1 2 [ 0 1 − 1 ] {\bm u}_1 = \frac{1}{\sqrt{3}} \begin{bmatrix} 1 \\ -1 \\ -1 \end{bmatrix}, \quad {\bm u}_2 = \frac{1}{\sqrt{6}} \begin{bmatrix} 2 \\ 1 \\ 1 \end{bmatrix}, \quad {\bm u}_3 = \frac{1}{\sqrt{2}} \begin{bmatrix} 0 \\ 1 \\ -1 \end{bmatrix} u1=3 1 111 ,u2=6 1 211 ,u3=2 1 011

取最大的两个特征值 λ 2 , λ 3 \lambda_2, \lambda_3 λ2,λ3对应的特征向量构成压缩空间,形成新的向量:

y ~ = [ u 2 u 3 ] ⊤ x = [ 2 / 6 1 / 6 1 / 6 0 1 / 2 − 1 / 2 ] x {\tilde {\bm y}} = \begin{bmatrix} {\bm u}_2 & {\bm u}_3 \end{bmatrix}^\top {\bm x} = \begin{bmatrix} 2/\sqrt{6} & 1/\sqrt{6} & 1/\sqrt{6} \\ 0 & 1/\sqrt{2} & -1/\sqrt{2} \end{bmatrix} {\bm x} y~=[u2u3]x=[2/6 01/6 1/2 1/6 1/2 ]x


K-L 变换为:

y = W ⊤ x = [ u 1 u 2 u 3 ] ⊤ x = [ 1 / 3 − 1 / 3 − 1 / 3 2 / 6 1 / 6 1 / 6 0 1 / 2 − 1 / 2 ] x {\bm y} = {\bm W}^\top {\bm x} = \begin{bmatrix} {\bm u}_1 & {\bm u}_2 & {\bm u}_3 \end{bmatrix}^\top {\bm x} = \begin{bmatrix} 1/\sqrt{3} & -1/\sqrt{3} & -1/\sqrt{3} \\ 2/\sqrt{6} & 1/\sqrt{6} & 1/\sqrt{6} \\ 0 & 1/\sqrt{2} & -1/\sqrt{2} \end{bmatrix} {\bm x} y=Wx=[u1u2u3]x= 1/3 2/6 01/3 1/6 1/2 1/3 1/6 1/2 x

http://www.xdnf.cn/news/8325.html

相关文章:

  • UG903 (v2024.2)的中文版
  • 楼宇自控系统赋予建筑智慧动能,加速智慧低碳建筑建设进程
  • 代码随想录算法训练营第四十八天
  • 昆仑芯超节点创新设计:1U 4 卡高密算力,无缝适配各类机房环境
  • Linux之Ext系列文件系统(含动静态库)
  • ansible剧本和角色的使用,部署lnmp
  • 搭建自己的语音对话系统:开源 S2S 流水线深度解析与实战
  • 李宏毅《深度学习》:Self-attention 自注意力机制
  • redis 进行缓存实战-18
  • 第J2周:ResNet50V2 算法实战与解析
  • Python爬虫(35)Python爬虫高阶:基于Docker集群的动态页面自动化采集系统实战
  • 内网渗透——红日靶场四
  • 从逻辑视角学习信息论:概念框架与实践指南
  • 127. 单词接龙
  • WDS 无线桥接
  • 交安安全员:交通工程安全领域的关键角色
  • 无人机桥梁检测如何通过数据存储、边缘AI、无线通讯等技术路线,提升检测效率
  • Seata分布式事物案例及详解
  • R语言开始绘图--柱状图
  • 业务场景中使用 SQL 实现快速数据更新与插入
  • MyBatis-Plus 中 QueryWrapper 的 Limit 实现
  • ceph osd 磁盘分区对齐
  • TCP与UDP区别及应用场景详解
  • 力扣HOT100之图论:200. 岛屿数量
  • 【LangChain大模型应用与多智能体开发 ① 初识LangChain 】
  • 用户意图识别模块
  • 跟Gemini学做PPT:字号选择
  • MyBatisPlus使用教程
  • Ubuntu 上进行树莓派交叉编译
  • hadoop 无法存储数据到hbase里面 已经解决