MIT线性代数：7.求解Ax=0：主变量、特解

1.零空间（Ax=0）

上面为矩阵A，然后我们对A进行消元，因为消元是行变换不会改变Ax=0的解，零空间也不会改变，会改变的是列空间

最后得到消元结果（上三角矩阵）

   

这个U又可以说是阶梯形式的矩阵(echelon form)，第一列和第三列为主元列，而其余的列为自由列，这时定义了一个rank秩，它是矩阵中主元个个数，本例rank=2，现在我们变成了Ux=0，解和零空间不变。

自由列的意思是，自由列上的变量可以自由的分配任意数值，因此我们第二列和第四列的变量x2和x4可以自由赋值，然后求解x1和x3即可。

我们先取x2=1，x4=0，利用回代可以求出x1和x3：

于是我们又发现乘以一个常数解仍然成立

同理我们再找一个特解（自由变量为我们赋予特定的值求出来的解叫特解）

于是Ax=0所有的解就出来了。零空间所包含的正好是特解的线性组合。我们发现特解的数目就是自由列的数目，如果一个矩阵为m*n，那么就有n个变量，去除主变量的个数r，那么自由列的个数为n-r。

总结算法：1.消元完得到主元的数量r，剩下n-r个自由变量。2.令这些自由变量为0和1得出特解。所有的特解构成了零空间的基，特解的线性组合即构成了整个零空间。

2.简化行阶梯（R(reduced)）

我们把上述的U进行化简：

最简的形式包含了所有信息：

1）主行（行一，行二）；

2）主列（列一，列三），自由列；

3）一个单位阵，主元上下均为0，而且主元为1，单位阵位于主列和主行的交汇处。以上是一个2×2的单位阵；

4）一个全为0的行，全为0的行总表示，该行的原行是其他行的线性组合；

5）从Ax=0变为Ux=0再变为Rx=0的解均相同。

我们可以进行列交换把主元列和自由列重新排列：

让单位矩阵在左侧，其余的在右侧，单位阵记作I，自由阵记作F，变成下图：

于是可以容易得出零空间矩阵N：

验证该零空间矩阵：

    

所以零空间为：

注意：如果在形成R=［ I，F］的时候进行了列交换，本例进行了第二列和第三列的交换，那么就要在N中交换第二行和第三行得出最后的解。

本节课最后一个例子：