矩阵A的零空间就Ax=0的解的集合。

零空间的求法:对矩阵A进行消元求得主变量和自由变量;给自由变量赋值得到特解;对特解进行线性组合得到零空间。

假设矩阵如下:

对矩阵A进行高斯消元得到上三角矩阵U,继续化简得到最简矩阵R:

由于方程Ax=0的右侧是零向量,所以只对矩阵A进行消元不会影响解,因此不需要增广矩阵,所以有:

从上面的高斯消元的结果可以看出,矩阵A的秩为2,其中第1,3列为主元列,2,4列为自由列,对应于方程主来说,形式转变如下:

从上式可以看出,x2,x4是自由变量,我们可以随意赋值,x2=0,x4=1;x2=1,x4=0可以分别得到两个特解(几个自由变量就有几个特解):

然后我们将两组特解进行线性组合就得到了矩阵A的零空间:

上面我们从数值解的角度描述了矩阵零空间的求法,下面从公式角度分析:

上面我们经过消元(行变换,不改变行空间和零空间,只改变列空间)得到了最简形式R。我们将R经过列变换得到如下矩阵:

我们可以对方程式作如下变形:

我们之所以进行上述变换,是为了有更好的表示形式(不进行列变换也行,但是要记住哪一列是单位矩阵I中的,哪一列是自由变量矩阵F中的):

这样我们代入方程式可以得到零空间矩阵:

从上面的推导可以看出,得到的零空间矩阵的每一列就是我们前面的特解(注意要变换顺序!交换第2,3行,结果便和前面相同)。因此,我们可以从通过消元法得到最简式R,然后就可以直接得到零空间矩阵,则零空间就是零空间矩阵各列向量的线性组合,而不需要像前面那样先给x2,x4赋值,然后回代到方程中得到两个特解,从而得到矩阵的零空间。

下面再举一例:

由于R本来就具有很好的形式,就不用进行列变换了:

于是通过解方程得到零空间矩阵:

注:最简矩阵R和零空间矩阵x在MATLAB中可以分别用命令rref(A),null(A,'r')得到

作者:nineheadedbird

matlab求零空间,【线性代数】矩阵的零空间相关推荐

  1. matlab矩阵保存到表格,将matlab求出的矩阵保存在Excel表格中

    将matlab求出的矩阵保存在Excel表格中 转贴Genial分享程序 function xlswrite(m,header,colnames,filename); % xlswrite Easil ...

  2. matlab可达矩阵 结果,matlab求可达矩阵

    matlab中可达矩阵的计算 2页 1下载券 邻接矩阵的可达矩阵及层... 2页 2下载券 邻接矩阵的一种计算方法 7页 免费 ISM算法(邻接矩阵求可达... 暂无...... JAVA矩阵算法代码 ...

  3. matlab求偏迹,矩阵的偏迹

    矩阵的偏迹运算是量子力学中的一种特殊的运算,它是一种特殊的矩阵迹运算.只不过这种取迹的过程并不是对全空间展开的,而是对某一个子空间. 定义 如果我们有两个希尔伯特空间分别记为HA和HB,它们可以分别用 ...

  4. 一分钟了解 Matlab求两个矩阵的相关程度corr2

    r = corr2(A,B)

  5. 关于matlab求含有nan矩阵的均值

    注意不是mean_a=nanmean(nanmean(a)) 应该是mean_a=nanmean(a(:))

  6. MATLAB求带符号的矩阵的逆

    syms a syms b B=[-sin(a) -sin(b)*cos(a) cos(b)*cos(a);cos(a) -sin(b)*sin(a),cos(b)*sin(a);0 cos(b) s ...

  7. 利用matlab软件求数列的中值,matlab求两个数列的方差

    matlab怎么求矩阵所有元素的期望和方差? 标准差s=std(X(1:end),flag)flag=0,采用1/(N-1)的系数,flag=1,采用1/(N)的系数 matlab用var函数算出的方 ...

  8. matlab参数群体赋值,matlab 带参数的矩阵赋值

    可以借助于符号函数. 代码: % 以前只用matlab求不带参数的矩阵的逆, % 请问能不能用matlab求带参数矩阵的逆, % 参考书上说可以,但是不知道怎么定义参数, % 实际操作显示错误,说参数 ...

  9. matlab求零空间,matlab求矩阵的零空间的一组整数基,该怎样操作?

    匿名用户 1级 2015-09-18 回答 第一部分:矩阵基本知识 一.矩阵的创建 直接输入法 利用Matlab函数创建矩阵 利用文件创建矩阵 二.矩阵的拆分 矩阵元素 矩阵拆分 特殊矩阵 三.矩阵的 ...

最新文章

  1. 陌陌股价过山车背后隐藏了什么?
  2. ubuntu 18.04 配置notebook远程连接的坑
  3. containerd 与安全沙箱的 Kubernetes 初体验
  4. 浅谈equals和==的区别
  5. java Arrays Generic
  6. java外围设计_Java 编程(23 种设计模式)
  7. 配置struts.xml时extends=struts-default会报错,原因和解决
  8. win7中能对窗口的排列方法是_win7系统窗口排列方式怎么修改?修改窗口排列方式方法...
  9. Flutter FocusNode输入框焦点控制概述
  10. Flutter 从配制开发环境再到开发第一个应用
  11. 机器学习资料合计(一)
  12. 基于氢探PowerECU的燃料电池控制系统开发经验
  13. gg 修改器游戏被保护_2.游戏so防封的制作思路
  14. 服务器重启后samba启动不了,提示stop: Unknown instance: smbd start/running, process 76585...
  15. 开发中经常用到的特效效果
  16. 【App自动化测试】(四)Capability配置参数说明
  17. 玩转扫雷(大战即将开始O(∩_∩))
  18. 有三个数abc 要求按大小顺序把它们输出_认识数字、比大小、凑十法… 孩子们做做手工玩玩游戏就OK了...
  19. 全球与中国口腔引导性骨再生(GTR)膜市场深度研究分析报告
  20. Android Studio Chipmunk 现已发布

热门文章

  1. 易语言学习笔记(一)
  2. a span做成按钮样式不选中文字
  3. 罗列几个Android插件化开发框架
  4. [转载] python字符串处理函数汇总
  5. [转载] [Python] np.ones_like(ndarray)和np.zeros_like(ndarray)
  6. [转载] python怎么获取redis中的数据_python操作redis数据库
  7. day006bootstrap的简单学习 + 轮播图
  8. 第四节:EasyUI的一些操作
  9. cookie与session的比较
  10. 14.Xcode8imageview图片圆角不显示的bug