首页 > 其他分享 >线性代数基础-矩阵奇异值分解-02

线性代数基础-矩阵奇异值分解-02

时间:2023-12-30 16:55:40浏览次数:38  
标签:02 翻折 矩阵 旋转 线性代数 奇异 正交 向量

目录

1. 引入

奇异值分解,singular value deconposition是6种矩阵分解方式中,综合性最强应用最广泛的分解技术,是PCA(主成分分析)的基础

六种矩阵分解技术:

只有矩阵为方阵(m=n)时,才有特征值;但对任何一个矩阵,都能求奇异值,SVD对所有矩阵均适用。

Ax=λx
矩阵乘向量(例如:Ax)是对向量做线性变换,如旋转,拉伸,翻折等。旋转与翻折可以通过正交矩阵(Orthogonal matrices)实现。
二维直角坐标系下旋转与翻折可通过如下两个正交矩阵描述:
旋转:

旋转矩阵
代表逆时针旋转
,此式在坐标变换或二维张量推导证明中常用。
对向量(1,0)与(0,1)进行旋转变换

翻折:

翻折矩阵 代表以
为对称轴,进行翻折。
对向量(1,0)与(0,1)进行翻折变换:

正交矩阵:
每列或每行只有与自身相乘时才为1,否则为0;


重点:对于正交矩阵,我们能瞬间求逆(只需转置),这非常漂亮(通常,矩阵求逆非常耗时)

谱聚类就用到这个知识点:
X1...Xn属于A1..Ak中的某一类,通过A1的维度来观察所有的样本,总共有|A1|个样本是属于A1这一类,对于属于这一类的值标上 1/sqrt(|A|) 这样的A1的向量描述W1就有 相乘=1

2. 几何的角度理解SVD

矩阵所代表的线性变换,可分解为先旋转(或翻折),再缩伸,最后再旋转(或翻折)的形式。
这三个过程至少具有一个,如对角阵只有缩拉。用矩阵语言描述如下:从右到左依次作用,先旋转
,再伸缩
,最后再旋转

详细说明:
1)对一个在
空间中的单位球,
作用为转换一下坐标基,球未变形;只是球旋转了一下而已。

2)
作用为将 n维球变为r维球(只有r个奇异值),再进行伸缩变换,此时球变形了。
3)
作用:再进行基的转变。

矩阵作用举例:

作用在向量
上:


作用在
上:

从这里一下子看不出。
若将原始坐标逆时针旋转45°,再将M作用于旋转后的矩阵,可得到如下结果:

-->

旋转过程实际上是寻找矩阵M特征向量.也就是寻找新的正交基,

其实也就是

这其实就是逆时针旋转45°
对于一组正交基V,经过矩阵M变化后,得到正交基U,要描述这一过程,就会得到奇异值分解。奇异值就是不同基变化的伸缩量。
这个结论太重要了

3. 空间的角度理解

设矩阵
秩为r(r个单位正交基)从矩阵
行空间(row space)中找一个向量
,通过
投影到A的列空间(column space)
方向,
v1 --> u1
其中, 放因子,此处与特征值类似,称为奇异值.

SVD妙处在于,在行空间中找的r个向量是正交的,即 彼此正交,且列向量空间的r个向量也是正交的,即 彼此正交




U为
V为:
对角阵

4 如何求解SVD


也就是如何求解:

U V都是正交矩阵:


为对称阵,可以进行分解:

因此:


特征向量归一化后组成的矩阵



特征向量归一化后组成的矩阵

奇异值的平方
或者特征值

奇异值为正数

例:


5. SVD的应用

将代表图像的矩阵进行SVD,只保留前k个较大的奇异值(
),从而保留了图像的主要信息。这样操作将大大减小图像存储的内存需求。这也可用于图像去噪,因为小的奇异值一般对应着噪声。用数学语言描述如下:前k个奇异值近似A效果是最好的

例如:若原始图像大小为mn,则压缩后大小为:mk+n*k+k。

原始图像:
取前20个奇异值就能保存其主要信息,则压缩后大小为:

压缩比为:7.67,压缩了近8倍


分别取前2,4,8,16,32,64,128个奇异值(当k=32时,已保留大部分信息)

标签:02,翻折,矩阵,旋转,线性代数,奇异,正交,向量
From: https://www.cnblogs.com/cavalier-chen/p/17936372.html

相关文章

  • 矩阵行列式
    定义与形式给定一个大小为\(n\timesn\)的矩阵\(A\),则行列式\[\det(A)=|A|=\sum_{p}(-1)^{\pi(p)}\prodA_{i,p_i}\]其中的\(p\)是一个\(1\simn\)的排列,\(\pi(p)\)为排列\(p\)的逆序对数。一些性质对于矩阵\(A\)来说:以主对角线为对称轴翻转之后的行列式和......
  • 2023 我的年度总结
    2023年,对于我来说是一个崭新的一年,今年我从专科毕业,专科的生活也让我受益良多,在专科时,我主要的学习方向是Java,学习的内容是数据结构和算法Java基础Mysql基础JavaWebSSMSpringBootSpringCloud和SpringCloudAlibabaRedisMysql调优的部分知识在我大二上学期学习......
  • 2024朋友圈跨年文案准备好啦~
    今年会遇到心软的财神爷吗?2023重在参与,2024重在搞钱。龙年愿望:大吉大利,有钱有你。祝我们发财被爱,一路好运常在。新年新气象,祝大家财源滚滚来,一路旺到底。今年目标,比去年健康,比去年快乐,比去年有钱。咔!这一年杀青了。希望2024,全糖去冰。人生何所求,暴富和自由。新年新对象,年年没谈上......
  • 每日总结2023年12月30日
    今天在配置Vue的路由,给大家分享一下经验和遇到的问题。首先我并不是提前安装好的router,是后面手动安装的。安装指令如下:[email protected]安装完毕之后配置路径新建router文件夹在src目录下,在该目录下新建一个index.js文件在一切配置完成之后我发现我的路由并没有起......
  • 2023.12.30做题纪要
    SAM模板评价:逆天纸糊串,学不会一点。#include<bits/stdc++.h>constintMAXN=3e6+100;intN;charch[MAXN];longlonganswer;classSuffix_Automaton{private:inttot,last,root;intchild[MAXN][26],link[MAXN],length[MAXN];longlongcnt......
  • 2023-12-30 训练总结
    返回C组做题,然后发现自己挂分了。T1寻找道路[NOIP2014提高组]寻找道路题目背景NOIP2014提高组D2T2题目描述在有向图\(G\)中,每条边的长度均为\(1\),现给定起点和终点,请你在图中找一条从起点到终点的路径,该路径满足以下条件:路径上的所有点的出边所指向的点都直接......
  • 02 AXI4总线axi-lite-slave
    软件版本:vitis2021.1(vivado2021.1)操作系统:WIN1064bit硬件平台:适用XILINXA7/K7/Z7/ZU/KU系列FPGA登录"米联客"FPGA社区-www.uisrc.com视频课程、答疑解惑!1概述使用XILINX的软件工具VIVADO以及XILINX的7代以上的FPGA或者SOC掌握AXI-4总线协议,并且可以灵活使用AXI-4总线技......
  • [软件测试] 02 白盒 逻辑覆盖测试 习题
    提示:篇幅较长,可以使用Ctrl+F,在页面中快速查询关键词(或者你要找的题目)并跳转到指定的位置。关键词:白盒测试,逻辑覆盖测试单选题以下不属于逻辑覆盖测试的是(C)A.语句覆盖B.判定覆盖C.基本路径覆盖D.条件覆盖根据不同的测试要求,逻辑覆盖测试可以分为语句覆盖、判定覆......
  • 【2023.12.30】PVE的PCIE直通改VGPU授权
    之前使用直通有个坏处,就是其他的CT和虚拟机用不了GPU,只能使用核显在这里参考的链接是https://gitlab.com/polloloco/vgpu-proxmoxaptupdateaptdist-upgradeaptinstall-ygitbuild-essentialdkmspve-headersmdevctlgitclonehttps://gitlab.com/polloloco/vgpu-prox......
  • [软件测试] 02 白盒 逻辑覆盖测试 习题
    提示:篇幅较长,可以使用Ctrl+F,在页面中快速查询关键词(或者你要找的题目)并跳转到指定的位置。关键词:白盒测试,逻辑覆盖测试单选题以下不属于逻辑覆盖测试的是(C)A.语句覆盖B.判定覆盖C.基本路径覆盖D.条件覆盖根据不同的测试要求,逻辑覆盖测试可以分为语句覆盖、判定覆......