学习视频来源:麻省理工公开课_线性代数导论 讲师:Gilbert Strang

http://open.163.com/special/opencourse/daishu.html


Lecture 7 求解Ax=0:主变量、特殊解

我们将从定义转向算法,求解Ax=0的算法是怎样的?这节的主要内容是零空间。

已知矩阵A=(1222246836810),解方程组{x1+2x2+2x3+2x4=02x1+4x4+6x3+8x4=03x1+6x2+8x3+10x4=0。在对矩阵进行消元的过程中,零空间不会改变(方程组的解不变),改变的是列空间。
麻省理工大学线性代数导论笔记 - Lecture 7 求解Ax=0:主变量、特殊解
- 首先经过消元得到一个echelon form 阶梯形式的矩阵U,本例中只有两个主元。现在我们得到了矩阵里最重要的数字:矩阵的rank 秩(主元的数量)

  • 要想求出Ux=0的解,我们进行下一步:找出pivot variables 主变量
    麻省理工大学线性代数导论笔记 - Lecture 7 求解Ax=0:主变量、特殊解
  • 本例中我们有两列pivot columns 主列,和剩下的两列free columns 自由列。这些自由列表示,可以自由或任意分配数值给未知数x2x4,即列2和列4的乘数是任意的。然后我们只需求解主变量x1x3

  • 假设给定x2=1x4=0,那么方程组
    Ux=0 <=> {x1+2x2+2x3+2x4=02x3+4x4=0变为{x1+2+2x3+0=02x3+0=0
    解得{x1=2x3=0

    假设给定x2=0x4=1,那么方程组变为{x1+0+2x3+2=02x3+4=0
    解得{x1=2x3=2

  • 现在我们有了两个解向量(2100)(2021)。这两个向量被称为special solutions 特殊解。special 特殊之处在于给自由变量分配的special numbers 特殊值(0、1)。通过取两个特殊解的所有线性组合,即构造出了整个零空间。那么特殊解的数量有多少呢?每个自由变量对应一个特殊解。而秩(r)表示主变量的个数,列数减去秩(nr)就是自由变量的个数。

    • 现在U是阶梯型矩阵,我们来将其进一步简化为reduced row echelon form 简化行阶梯形式R,向上消元,将主元上下全变为 0,并简化为 1。
      麻省理工大学线性代数导论笔记 - Lecture 7 求解Ax=0:主变量、特殊解
  • R以最简形式包含了所有信息:主元,主行,主列,以及一个单位阵(位于主行和主列交汇处)。现在我们求解Rx=0 <=> {x1+2x22x4=0x3+2x4=0(从Ax=0Ux=0Rx=0,消元不改变解)。
  • 当我们给自由变量分配特殊值并回代后,主列构成的单位阵 I=(1001),自由列构成的部分 F=(2202),恰能构成两个特殊解(符号相反):(2100)(2021)

    • 抽象一点来证明上述的结论:假设有简化行阶梯形式矩阵R=(IF00),主列在前,自由列在后,下面是一些零行。Ir×r矩阵,有r个主元,nr个自由列。
  • 要满足Rx=0,首先要找到一个null space matrix 零空间矩阵(N(A),使得RN=0N的各列由特殊解组成。那么易看出N=(FI)I放在自由变量部分,F放在主变量部分。

  • 再给出另一个例子把这个算法过一遍。

    已知矩阵A=(1232462682810),消元至U=(123022000000),可得r=2,2列主列,1列自由列。继续消元至R=(101011000000),可得I=(1001)F=(11)。则A的零空间N(A)=c(111)c为常数)。回代Ax=0验证无误。

相关文章:

  • 2022-12-23
  • 2021-04-12
  • 2022-12-23
  • 2022-12-23
  • 2021-12-29
  • 2021-09-11
  • 2022-12-23
  • 2022-02-09
猜你喜欢
  • 2021-10-16
  • 2021-11-17
  • 2022-12-23
  • 2021-06-29
  • 2021-06-24
  • 2021-12-10
  • 2021-09-18
相关资源
相似解决方案