谭升
非主流人工智能科学家 我和外面那些妖艳的货不一样

【线性代数】3-4:方程组的完整解( $Ax=b$ )

方程组的完整解

$Ax=b$

之前我们已经研究了 $Ax=0$的相关内容,值得说一下的是,列空间和nullspace是有些区别的,列空间指的是b所在的空间,而nullspace是x所在的空间,这个要区别一下,这些所有空间都是针对矩阵的。

特解(Particular Solution)

搞不懂Particular Solution和Sceptical Solution有啥区别的可以仔细看看了,之前我也没发现其有什么根本不同,Particular Solution是把所有的free variables设置为0,来一个完整的例子
$$
\begin{bmatrix}
1&3&0&2\\
0&0&1&4\\
1&3&1&6
\end{bmatrix}
\begin{bmatrix}
x_1\\x_2\\x_3\\x_4
\end{bmatrix}
=\begin{bmatrix}
1\\6\\7
\end{bmatrix}
$$

Augmented Matrix:

$$
\begin{bmatrix}
1&3&0&2&1\\
0&0&1&4&6\\
1&3&1&6&7
\end{bmatrix}
=\begin{bmatrix}
A&b
\end{bmatrix}
$$

经过消元,对原始矩阵方程消元和对Augment Matrix消元得到结果相同,这里就只写augment matrix:

$$
\begin{bmatrix}
1&3&0&2&1\\
0&0&1&4&6\\
0&0&0&0&0
\end{bmatrix}
=\begin{bmatrix}
R&d
\end{bmatrix}
$$
通过定位pivot可以确定free columns是第二列和第四列,如果我们把$x_2,x_4$设置为零,那么就得到了x为
$$
\begin{bmatrix}
x_1\\x_2\\x_3\\x_4
\end{bmatrix}
=\begin{bmatrix}
1\\0\\6\\0
\end{bmatrix}
=x_{particular}
$$
可以通过回代验证$Rx=d$或者$Ax=b$这里的free variables都是0,那么就说这个解是particular的,而且可以发现pivot对应的x位置与d有关系,按顺序就是d对应的元素,这不是巧合,观察矩阵可以得到相应的结论。


这是我们到目前位置讨论的方程组的两种解,一个当b=0的时候,一个当b不等于0的时候,如果我们把这两个方程组左右相加,那么就得到
$$
Ax_p=b\\
Ax_n=0\\
Ax_p+Ax_n=b+0\\
A(x_p+x_n)=b
$$
那么完整的解(式子中的数字来自上面的例子,nullspace我没有写出来,大家可以自行验证):

为什么完整解是上看的式子呢,可以看下一节的详细介绍。

完整解(The Complete Solution)

前面我们确定了完整的解就是$x_p+x_n$,那么我们到底有多少个解呢?

序号 m&r n&r Matrix Shape $Ax=b$ Solution
1 $r=m$ $r=n$ Square and Invertible $Ax=b$ 1
2 $r=m$ $r<n$ Short and Wide $Ax=b$ $\infty$
3 $r<m$ $r=n$ Tall and Thin $Ax=b$ 0 or 1
4 $r<m$ $r<n$ Not full rank $Ax=b$ 0 or $\infty$

这四种情况,简单讲解下:


第一种:是方阵,我们最标准的方程组,rank与m,n相等,那么就一个解。这种情况下,消元后的R是单位矩阵:
$$
R=\begin{bmatrix}I\end{bmatrix}
$$


第二种:如果方程组的个数小于未知数的个数,而rank与行相同,$rank=m<n$ ,消元后的R是如下矩阵:
$$
R=\begin{bmatrix}I&F\end{bmatrix}
$$


第三种:如果方程组的个数大于未知数的个数,而rank与列数相同,$rank=n<m$ ,消元后的R是如下矩阵:
$$
R=\begin{bmatrix}I\\0\end{bmatrix}
$$
如果b与R下面全是0的行对应的行不是0,那么就没有解


第四种:如果有效的方程组的个数小于未知数的个数,与情况二相似,但是有效方程数还小于远方程数,$rank<m$ , $rank<n$ ,消元后的R是如下矩阵:

$$
R=\begin{bmatrix}I&F\\0&0\end{bmatrix}
$$
如果b与R下面全是0的行对应的行不是0,那么就没有解

自此,$Ax=b$ 的所有解以及对应的情况都已经搞定了,但是为啥complete solution是$x_p+x_n$呢?这就是个悬案了,后面一定给出答案,现在就暂时记住就好了(回归课堂教育,先记住,就没有然后了)。

Conclusion

这篇也比较凌乱,不像写算法那种博客很流畅,解释,代码,总结这种套路在这不太实用,因为知识错综复杂,本文上面的四种解的情况才是重点,包括什么时候没有解,什么时候无数个解,如果按照表和矩阵R的形式来看,就会豁然开朗,后面继续。。。

Share

You may also like...

说点什么

avatar
  Subscribe  
提醒

试录了几节线性代数视频课程,欢迎大家支持,
点击试看
点击购买更多内容

由于博客移至wordpress,部分公式和代码显示不正常,博主正在努力修改,如发现公式显示错误,请及时在文章下留言,感谢您的帮助,尽请原谅!