降维处理:PCA和LDA

news/2025/1/11 10:54:33/

1,主成分分析

主成分分析(Principal Component Analysis,PCA)是一种无监督问题,是降维中最常用的一种手段,基于方差(方差最大的方向)提取最有价值的信息再进一步分类。降维后数据特征的意义发生变化。

2,向量的表示及基变换

2.1 内积:
在这里插入图片描述设向量B的模为1,则A与B的内积值等于A向B所在直线投影的矢量长度:
在这里插入图片描述
在这里插入图片描述

2.2 向量的线性组合

向量可以表示为(3,2)实际上表示线性组合:
在这里插入图片描述
其中(1,0)和(0,1)叫做二维空间中的一组基
在这里插入图片描述
2.3 基变换

基是正交的(即内积为0,或直观说相互垂直,线性无关),更好的表达坐标轴上的数据。
在这里插入图片描述
变换:数据与一个基做内积运算,结果作为第一个新的坐标分量,然后与第二个基做内积运算,结果作为第二个新坐标的分量,比如:
在这里插入图片描述
推广到多维:

在这里插入图片描述两个矩阵相乘的意义是将右边矩阵a中的每一列列向量ai变换到左边矩阵p中每一行行向量pi为基所表示的空间中去(尽可能去除线性相关性)。

2.4 协方差矩阵

方向:选择的这个方向(或者说基)应该尽量保留最多的原始信息,即一种直观的看法是:希望投影后的投影值尽可能分散。

方差(单个特征发散的程度):
在这里插入图片描述
寻找一个一维基,使得所有数据变换为这个基上的坐标表示后,方差值最大(尽可能分散)

协方差(假设均值为0时(以0为中心化),两个特征之间的关系):

在这里插入图片描述
如果单纯只选择方差最大的方向,后续方向应该会和方差最大的方向接近重合。而为了让两个字段尽可能表示更多的原始信息,我们是不希望它们之间存在(线性)相关性的,可以用两个字段的协方差表示其相关性,当协方差为0时,表示两个字段完全独立。为了让协方差为0,选择第二个基时只能在与第一个基正交的方向上选择,然后在保证第二个基在该方向上的方差尽可能大。因此最终选择的两个方向一定是正交的。

3,优化目标

将一组N维向量降为K维(K大于0,小于N),目标是选择K个单位正交基,使原始数据变换到这组基上后,各字段两两间协方差为0,字段的方差则尽可能大。

协方差矩阵:

在这里插入图片描述m代表数据的个数,x为数据本身,矩阵对角线上的两个元素分别是两个字段各自的方差,而其它元素是a和b的协方差

协方差矩阵对角化:即除对角线外的其它元素化为0(协方差为0),并且在对角线上将元素按大小从上到下排列(方差尽可能大)
在这里插入图片描述
一个n行n列的实对称矩阵一定可以找到n个单位正交特征向量:
在这里插入图片描述
实对称阵可进行对角化:
在这里插入图片描述
根据特征值的从大到小,将特征向量从上到下排列,则用前K行组成的矩阵乘以原始数据矩阵X,就得到了我们需要的降维后的数据矩阵Y

4,PCA实例

(1)原始数据共5个数据,每个数据有2个特征,并进行0中心化
(2)计算协方差矩阵,对角线为特征的方差,非对角线为特征的协方差(特征间的关系):
(3)进行矩阵分解求特征值和特征向量,找出对应特征值最大的特征向量,并对该向量进行归一化处理(基)
(4)将基向量与原始数据作内积得到降维后的结果
在这里插入图片描述5,线性判别分析

线性判别分析(Linear Discriminant Analysis,LDA)是由Ronald A. Fisher在1936年提出了线性判别方法,用于数据预处理中的降维,完成分类任务,LDA关心的是能够最大化类间区分度的坐标轴成分,即将特征空间(数据集中的多维样本)投影(找到更合适分类的空间k)到一个维度更小的 k 维子空间中,同时保持区分类别的信息。所以与PCA不同,LDA更关心分类而不是方差。

5.1 原理
LDA是“有监督”的,它计算的是另一类特定的方向,投影到维度更低的空间中,使得投影后的点,会形成按类别区分,一簇一簇的情况,相同类别的点,将会在投影后的空间中更接近方法
在这里插入图片描述原始数据:
在这里插入图片描述
找到合适的空间进行投影:
在这里插入图片描述
变换数据(投影后的数据):
在这里插入图片描述
从上述来看,LDA分类的一个目标是使得不同类别之间的距离越远越好,同一类别之中的距离越近越好。

每类样例的均值:
在这里插入图片描述

投影后的均值:

在这里插入图片描述
投影后的两类样本中心点尽量分离:
在这里插入图片描述
但是只是最大化J(w),并不能满足要求,如下图,X1的方向可以最大化J(w),但是却分的不好
在这里插入图片描述
散列值描述了样本点的密集程度,值越大,越分散,反之,越集中。同类之间应该越密集些:
在这里插入图片描述
故目标函数为:

在这里插入图片描述
散列值公式展开:
在这里插入图片描述
散列矩阵:
在这里插入图片描述
类内散布矩阵Sw = S1+S2:

在这里插入图片描述
目标函数的分子展开:
在这里插入图片描述??称作类间散布矩阵,故最终目标函数:
在这里插入图片描述
分母进行归一化,如果分子、分母是都可以取任意值的,那就会使得有无穷解,故将分母限制为长度为1

拉格朗日乘子法:
在这里插入图片描述
两边都乘以Sw的逆:
在这里插入图片描述


http://www.ppmy.cn/news/611809.html

相关文章

微信小程序input使用

这个input 自己也是花费了2个多小时,熟悉,下面将一一讲解经常使用的里面的属性 value 这个value 就是input中的内容, 如果我们想获取input中输入的内容,value是必须使用的 使用可以在input中 定义 value{{phone}} 然后在data中初始化就行了,如果写死的input里面就是死的数据…

h265player开发

h265player开发 https://github.com/goldvideo/h265player 简介 随着视频编码技术的发展,相比H.264, H.265同等画质体积仅为一半、带宽占用省一半、画质更细腻等诸多优势。 但Web浏览器还不支持H.265的解码播放,因此基于Web Assembly(封装FFmpeg)、JS…

Android Jetpack Compose - 深入了解 AlertDialog

在开发 Android 应用时,对话框是一个常见的 UI 元素,用于向用户显示信息或获取用户的反馈。在 Jetpack Compose 中,我们可以使用 AlertDialog 组件来创建对话框。在这篇博客中,我们将深入了解如何使用 Jetpack Compose 的 AlertDi…

寒假作业pta3

某地老鼠成灾,现悬赏抓老鼠,每抓到一只奖励10元,于是开始跟老鼠斗智斗勇:每天在墙角可选择以下三个操作:放置一个带有一块奶酪的捕鼠夹(T),或者放置一块奶酪(C),或者什么也不放(X)。捕鼠夹可重复…

微信小程序点击图片切换图片

实现的效果图 实现的方式有2中 先说下思路 1 使用wx:if 写2个图片 点击的显示一张即可, 2 使用数据传递当拿到图片的时候更换图片即可. 看下代码 2中方式:xml <view class"intro"> <image catchtapopen_tap wx:if"{{isShow}}" src"/im…

h265webplayer

h265webplayer https://github.com/ksvc/h265webplayer h265webplayer是金山云的Web端H.265视频播放器&#xff0c;该播放器Web SDK让您可以在支持WebAssembly的浏览器上播放MP4格式的点播视频&#xff0c;FLV http-flv协议的直播视频。 支持的功能 1、mp4格式的点播&#…

java后台的微信小程序支付的解决方案

java的后台主要就写两个controller 一个是 ----------------(1)调用支付统一下单API 生成 prepay_id 二个是 ----------------(2)微信支付完成的回调controller 第一个controller需要传入的参数有如下&#xff0c;如图&#xff1a; 然后&#xff0c;发送到微信端&#xff0c;把…

微信小程序bindtap 与 catchtap 是使用

如果写小程序对二者不理解的,那看到这边博客,将很快帮助到您, 个人总结的一句话:,bindtap点击事件在同一个view中会向上冒泡,而catchtap 不会向上冒泡 下面会有一个demo给出解释, 说他们使用的时候先说下微信小程序的事件分类冒泡事件与非冒泡事件 官网上这样规定的 事件分…