码农知识堂 - 1000bd
  •   Python
  •   PHP
  •   JS/TS
  •   JAVA
  •   C/C++
  •   C#
  •   GO
  •   Kotlin
  •   Swift
  • 主成分分析及matlab实现


    主成分分析及matlab实现

    • 一、主成分分析的基本原理
    • 二、定义
    • 三、主成分分析的计算步骤
      • 1. 将原始数据标准化
      • 2. 建立变量的相关系数阵:
      • 3. 求R的特征根及相应的单位特征向量:
      • 4.写出主成分
      • 5.计算主成分贡献率及累计贡献率
    • 四、 实例
    • 五、主成分分析的应用

    在实际问题研究中,多变量问题是经常会遇到的。变量太多,无疑会增加分析问题的难度与复杂性,而且在许多实际问题中,多个变量之间是具有一定的相关关系的。因此,人们会很自然地想到,能否在相关分析的基础上,用较少的新变量代替原来较多的旧变量,而且使这些较少的新变量尽可能多地保留原来变量所反映的信息?
    事实上,这种想法是可以实现的,主成分分析方法就是综合处理这种问题的一种强有力的工具。主成分分析是把原来多个变量划为少数几个综合指标的一种统计分析方法。从数学角度来看,这是一种降维处理技术。

    例如,某人要做一件上衣要测量很多尺寸,如身长、袖长、胸围、腰围、肩宽、肩厚等十几项指标,但某服装厂要生产一批新型服装绝不可能把尺寸的型号分得过多
    ?而是从多种指标中综合成几个少数的综合指标,做为分类的型号,利用主成分分析将十几项指标综合成3项指标,一项是反映长度的指标,一项是反映胖瘦的指标,一项是反映特体的指标。

    一、主成分分析的基本原理

    假定有n个样本,每个样本共有p个变量,构成一个n×p阶的数据矩阵
    在这里插入图片描述
    当p较大时,在p维空间中考察问题比较麻烦。为了克服这一困难,就需要进行降维处理,即用较少的几个综合指标代替原来较多的变量指标,而且使这些较少的综合指标既能尽量多地反映原来较多变量指标所反映的信息,同时它们之间又是彼此独立的。

    二、定义

    记x1,x2,…,xP为原变量指标,z1,z2,…,zm(m≤p)为新变量指标
    在这里插入图片描述
    系数lij的确定原则:
    ① zi与zj(i≠j;i,j=1,2,…,m)相互无关;
    ② z1是x1,x2,…,xP的一切线性组合中方差最大者,z2是与z1不相关的x1,x2,…,xP的所有线性组合中方差最大者;
    zm是与z1,z2,……,zm-1都不相关的x1,x2,…xP, 的所有线性组合中方差最大者。
    则新变量指标z1,z2,…,zm分别称为原变量指标x1,x2,…,xP的第1,第2,…,第m主成分。
    从以上的分析可以看出,主成分分析的实质就是确定原来变量xj(j=1,2 ,…, p)在诸主成分zi(i=1,2,…,m)上的荷载 lij( i=1,2,…,m; j=1,2 ,…,p)。 从数学上可以证明,它们分别是相关矩阵m个较大的特征值所对应的特征向量。

    三、主成分分析的计算步骤

    设有 n 个样品,每个样品观测 p 个指标,将原始数据写成矩阵
    在这里插入图片描述

    1. 将原始数据标准化

    这里不妨设上边矩阵已标准化了

    2. 建立变量的相关系数阵:

    在这里插入图片描述
    在这里插入图片描述

    3. 求R的特征根及相应的单位特征向量:

    在这里插入图片描述
    在这里插入图片描述

    4.写出主成分

    在这里插入图片描述
    在这里插入图片描述

    5.计算主成分贡献率及累计贡献率

    在这里插入图片描述

    四、 实例

    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

    五、主成分分析的应用

    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

    在这里插入图片描述
    在这里插入图片描述

    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

  • 相关阅读:
    [python]十九、进程、线程和协程
    谈谈数据决策平台搭建的必要性
    swagger在线api文档搭建指南
    More effective C++:条款1.仔细区别pointers和references及条款2:最好使用C++转型操作符
    基于虚拟同步发电机控制的双机并联Simulink仿真模型
    Python花瓣雨
    YOLOv7模型训练
    反转链表-js
    【滤波器】基于matlab实现微波带低通高通带通滤波器设计
    beyond Compare连接 openWrt 和 VsCode
  • 原文地址:https://blog.csdn.net/weixin_56462041/article/details/126468168
  • 最新文章
  • 攻防演习之三天拿下官网站群
    数据安全治理学习——前期安全规划和安全管理体系建设
    企业安全 | 企业内一次钓鱼演练准备过程
    内网渗透测试 | Kerberos协议及其部分攻击手法
    0day的产生 | 不懂代码的"代码审计"
    安装scrcpy-client模块av模块异常,环境问题解决方案
    leetcode hot100【LeetCode 279. 完全平方数】java实现
    OpenWrt下安装Mosquitto
    AnatoMask论文汇总
    【AI日记】24.11.01 LangChain、openai api和github copilot
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1
正则表达式工具 cron表达式工具 密码生成工具

京公网安备 11010502049817号