CP.12图、网络、关联矩阵

1.图

图经常被应用于各种领域，用结点和边来描述变量之间的关系

上图是具有4个节点5条边的图，可以用矩阵来表示，这门课里面把矩阵称作关联矩阵，其他地方称作邻接矩阵，总之是个矩阵用来描述边与节点的关系。
$$A=\left[\begin{array}{cccc}-1& 1& 0& 0\\ 0& -1& 1& 0\\ -1& 0& 1& 0\\ -1& 0& 0& 1\\ 0& 0& -1& 1\end{array}\right]$$
关联矩阵的每一行代表路径，每一列代表节点。例如，第1行$\left[\begin{array}{cccc}-1& 1& 0& 0\end{array}\right]$是指路径1从节点1流出，流入节点2，这里将流出表达为“-1”，流入表达为“1”。这是一种有向图，无向图的矩阵表达只有相连或者不相连。

1.1零空间

上图中可以看作是电路节点与电压电流的关系，以此为基础考察矩阵的零空间，也就是$AX=0$
$$AX=\left[\begin{array}{cccc}-1& 1& 0& 0\\ 0& -1& 1& 0\\ -1& 0& 1& 0\\ -1& 0& 0& 1\\ 0& 0& -1& 1\end{array}\right]\left[\begin{array}{c}{x}_1\\ x_2\\ x_3\\ x_4\end{array}\right]=0$$
现在问题转换成了求解一组节点电势，能够满足：
$$AX=\left[\begin{array}{c}{x}_2-x_1\\ x_3-x_2\\ x_3-x_1\\ x_4-x_1\\ x_4-x_3\end{array}\right]=0$$
也就是要求解能让各点电势平衡的X。通过MATLAB的rank(A)命令，可以求去矩阵的秩为3，因此有一个自由列，$dimN(A)=1$。观察可知第4列是自由列，对应$x_4$是自由变量，赋值$x_4=1$可以解得$$X=\left[\begin{array}{c}1\\ 1\\ 1\\ 1\end{array}\right]$$
因此矩阵$A$的零空间是$cX$，说明当各点的电势相等时，各支路上不会有电流。

1.2左零空间

从前面的知识可以知道，左零空间的表达是$A^{T}y=0$，矩阵$A^T$有5列，$rank(A^T)=3$，因此左零空间的维度是2。因此$A^T$中有存在2种线性相关的列，$A^T$的列就是图中的路径，路径相关代表形成了环路，如节点1，2，3。
$$A^{T}y=\left[\begin{array}{ccccc}-1& 0& -1& -1& 0\\ 1& -1& 0& 0& 0\\ 0& 1& 1& 0& -1\\ 0& 0& 0& 1& 1\end{array}\right]\left[\begin{array}{c}{y}_1\\ y_2\\ y_3\\ y_4\end{array}\right]=0$$
用$y_i$代表支路上的电流，那么方程表征的就是基尔霍夫电流定律了。