线性代数_矩阵零空间的维度与奇异值的关系
1. 矩阵的零空间
矩阵A的零空间就Ax=0的解的集合。
零空间的求法:对矩阵A进行消元求得主变量和自由变量;给自由变量赋值得到特解;对特解进行线性组合得到零空间。
假设矩阵如下:
对矩阵A进行高斯消元得到上三角矩阵U,继续化简得到最简矩阵R:
由于方程Ax=0的右侧是零向量,所以只对矩阵A进行消元不会影响解,因此不需要增广矩阵,所以有:
从上面的高斯消元的结果可以看出,矩阵A的秩为2,其中第1,3列为主元列,2,4列为自由列,对应于方程主来说,形式转变如下:
从上式可以看出,x2,x4是自由变量,我们可以随意赋值,
或者
可以分别得到两个特解(几个自由变量就有几个特解):
,
然后我们将两组特解进行线性组合就得到了矩阵A的零空间:
上面我们从数值解的角度描述了矩阵零空间的求法,下面从公式角度分析:
上面我们经过消元( 行变换,不改变行空间和零空间,只改变列空间)得到了最简形式R。我们将R经过列变换得到如下矩阵:
我们可以对方程式作如下变形:
我们之所以进行上述变换,是为了有更好的表示形式(不进行列变换也行,但是要记住哪一列是单位矩阵I中的,哪一列是自由变量矩阵F中的):
这样我们代入方程式可以得到零空间矩阵:
所以,可得:
从上面的推导可以看出,得到的零空间矩阵的每一列就是我们前面的特解(注意要变换顺序!交换第2,3行,结果便和前面相同)。因此,我们可以从通过消元法得到最简式R,然后就可以直接得到零空间矩阵,则 零空间就是零空间矩阵各列向量的线性组合,而不需要像前面那样先给x2,x4赋值,然后回代到方程中得到两个特解,从而得到矩阵的零空间。
下面再举一例:
由于R本来就具有很好的形式,就不用进行列变换了:
于是通过解方程得到零空间矩阵:
注:最简矩阵R和零空间矩阵x在MATLAB中可以分别用命令rref(A),null(A,'r')得到
2. 矩阵的零空间和奇异值
什么叫做零空间的维度?是未知量x的维度(或者说未知数的个数n)吗?NO.
- 零空间就是齐次线性方程组Ax=0的全部解,
- 零空间的基就是基础解系,
- 零空间的维数是n-r(A), n是未知元的个数, r是A的秩.
- 性质:矩阵的秩等于矩阵的非零奇异值的个数(注:矩阵的秩不一定等于非零特征值的个数,例如非零的幂零矩阵的所有特征值都为0,但它的秩不为0)
- 所以:可以采用SVD分解,求出奇异值,奇异值中有几个0,那么零空间的维数就是几。
矩阵A
奇异值分解如下
S对角线的元素就是奇异值的3个值,零空间则是指Ax=0时,x的取值,有几个基础解系就代表这个零空间有几个维度。
如上图,秩为2,因此基础解析个数为 3 - 2 = 1,因此零空间中只有一个基础解系,也代表零空间只有一个维度。
参考:
奇异值和零空间_AutoGalaxy的博客-CSDN博客_由奇异值矩阵求零空间
【线性代数】矩阵的零空间_nineheaded_bird的博客-CSDN博客_线性代数零空间
矩阵的奇异值分解 - 知乎
线性代数_矩阵零空间的维度与奇异值的关系相关推荐
- python矩阵运算与线形代数_[译] 线性代数:矩阵基本运算
线性代数:矩阵基本运算 在本文中,我们将介绍矩阵的大部分基本运算,依次是矩阵的加减法.矩阵的标量乘法.矩阵与矩阵的乘法.求转置矩阵,以及深入了解矩阵的行列式运算.本文将不会涉及逆矩阵.矩阵的秩等概念, ...
- 矩阵乘法 算法训练 试题_线性代数入门——矩阵乘法的定义及其意义
系列简介:这个系列文章讲解线性代数的基础内容,注重学习方法的培养.线性代数课程的一个重要特点(也是难点)是概念众多,而且各概念间有着千丝万缕的联系,对于初学者不易理解的问题我们会不惜笔墨加以解释.在内 ...
- matlab求零空间,【线性代数】矩阵的零空间
矩阵A的零空间就Ax=0的解的集合. 零空间的求法:对矩阵A进行消元求得主变量和自由变量:给自由变量赋值得到特解:对特解进行线性组合得到零空间. 假设矩阵如下: 对矩阵A进行高斯消元得到上三角矩阵U, ...
- 【线性代数】矩阵的零空间
矩阵A的零空间就Ax=0的解的集合. 零空间的求法:对矩阵A进行消元求得主变量和自由变量:给自由变量赋值得到特解:对特解进行线性组合得到零空间. 假设矩阵如下: 对矩阵A进行高斯消元得到上三角矩阵U, ...
- 【线性代数】矩阵及其特性
[线性代数]矩阵及其特性 写在前面 只拉伸不旋转的方向 特征值和特征向量 相似和对角化 正交,对称矩阵 拉伸最大的方向 二次型理论 从曲线而来 合同矩阵 正定二次型 参考资料 写在前面 本文是笔者用于 ...
- 线性代数学习之向量空间,维度,和四大子空间
空间: 在上一次https://www.cnblogs.com/webor2006/p/14306046.html学习了诸多在线性代数中非常核心的概念(线性组合.线性相关.线性无关.生成空间,空间的基 ...
- Codeforces 947E Perpetual Subtraction (线性代数、矩阵对角化、DP)
Codeforces 947E Perpetual Subtraction (线性代数.矩阵对角化.DP) 手动博客搬家: 本文发表于20181212 09:37:21, 原地址https://blo ...
- 线性代数之 矩阵的特征值,特征向量,特征分解
线性代数之 矩阵的特征值,特征向量和特征分解 前言 特征值和特征向量 求矩阵特征值 矩阵的特征分解 补充:实对称矩阵 后记 前言 矩阵的特征分解是比较基础的知识了,但是应用却十分广泛,比如主成分分析. ...
- Python中使用元组对ndarray矩阵的某个维度进行选取和调序的操作
大家应该都知道,numpy库中的ndarray可以用所谓的三帽号规则进行切片操作(详情可见我之前写的博文,链接 https://blog.csdn.net/wenhao_ir/article/deta ...
最新文章
- 网络模型:(0)背景
- 腾讯会议用户突破1亿,发布企业版最高支持2000人同时参会
- Activity传递对象的方法
- 程序人生之回顾大学前两年----第一篇
- IDC报告:谁是桌面虚拟化的王者
- python 线程中出现执行错乱_多处理会导致Python崩溃,并在调用fork()时在另一个线程中出现错误...
- 【maven】pom常用配置
- Python3 使用 pymysql 连接 MySQL 建表时出现 Warning3719 UTF8 警告
- xposed框架在分机安装失败_pokemon go 和 虚拟gps
- linux无线电工具grax,开源软件无线电GNU Radio
- 史话上:量子物理学的前世今生
- Mac批量恢复废纸篓文件
- 鹏城实验室“鲲鹏”人才引育计划(2022-2025年)
- UPC Haywire(模拟退火 || 随机数法)
- 计算机从一级到四级要学多久,从一级到四级 全国计算机等级考试全攻略
- 仙童半导体和“八叛逆”所缔造的硅谷模式
- 图像化界面开发之QT入门
- filezilla检查新版本mac_FileZilla for Mac(ftp客户端工具)中文版
- matlab关系运算
- Ubuntu16.04深度学习环境之TensorFlow1.4CPU/GPU安装实测