1. A X = 0 AX=0 AX=0
一个图可以由节点和边组成,假设我们有一个节点notes :n=4,边edges:m=5的有向图,表示如下
- 通过以上电路图可以得到关联矩阵(incident matrix),我们定义边,开始端用-1表示,结束端用1表示,用行表示边,用列表示节点。可得如下
比如edge1是由从节点1出发,到节点2 ,用行向量表示可得
[ − 1 1 0 0 ] ⇒ e d g e 1 (1) \begin{bmatrix}-1&1&0&0\end{bmatrix}\Rightarrow edge1\tag{1} [−1100]⇒edge1(1) - 所以矩阵A根据上图可得如下:
A = [ − 1 1 0 0 0 − 1 1 0 − 1 0 1 0 − 1 0 0 1 0 0 − 1 1 ] (2) A=\begin{bmatrix}-1& 1& 0& 0\\\\ 0 &-1& 1& 0\\\\ -1& 0& 1& 0\\\\ -1& 0& 0& 1\\\\ 0 & 0&-1& 1\end{bmatrix}\tag{2} A= −10−1−101−10000110−100011 (2)-
我们可以看到由边1,2,3组成一个子图,那么看看由1,2,3组成的矩阵表示如下:
[ − 1 1 0 0 0 − 1 1 0 − 1 0 1 0 ] ⇒ r o w 1 + r o w 2 = r o w 3 (3) \begin{bmatrix}-1& 1& 0& 0\\\\ 0 &-1& 1& 0\\\\ -1& 0& 1& 0 \end{bmatrix}\Rightarrow row_1+row_2=row_3\tag{3} −10−11−10011000 ⇒row1+row2=row3(3)
也就是说,如果在矩阵中3行能达到等式关系,可以等效与在图上形成一个子图。
- 那么矩阵A的零空间代表什么呢?我们将AX=0 用列的形式表示如下:
A X = [ − 1 1 0 0 0 − 1 1 0 − 1 0 1 0 − 1 0 0 1 0 0 − 1 1 ] [ x 1 x 2 x 3 x 4 ] = [ x 2 − x 1 x 3 − x 2 x 3 − x 1 x 4 − x 3 ] (4) AX=\begin{bmatrix}-1& 1& 0& 0\\\\ 0 &-1& 1& 0\\\\ -1& 0& 1& 0\\\\ -1& 0& 0& 1\\\\ 0 & 0&-1& 1\end{bmatrix}\begin{bmatrix}x_1\\\\x_2\\\\x_3\\\\x_4\end{bmatrix}=\begin{bmatrix}x_2-x_1\\\\x_3-x_2\\\\x_3-x_1\\\\x_4-x_3\end{bmatrix}\tag{4} AX= −10−1−101−10000110−100011 x1x2x3x4 = x2−x1x3−x2x3−x1x4−x3 (4)
A X = [ − 1 0 − 1 − 1 0 ] x 1 + [ 1 − 1 0 0 0 ] x 2 + [ 0 1 1 0 − 1 ] x 3 + [ 0 0 0 1 1 ] x 4 (5) AX=\begin{bmatrix}-1\\\\0\\\\-1\\\\-1\\\\0\end{bmatrix}x_1+\begin{bmatrix}1\\\\-1\\\\0\\\\0\\\\0\end{bmatrix}x_2+\begin{bmatrix}0\\\\1\\\\1\\\\0\\\\-1\end{bmatrix}x_3+\begin{bmatrix}0\\\\0\\\\0\\\\1\\\\1\end{bmatrix}x_4\tag{5} AX= −10−1−10 x1+ 1−1000 x2+ 0110−1 x3+ 00011 x4(5)
- 那么矩阵A的零空间代表什么呢?我们将AX=0 用列的形式表示如下:
-
- AX=0的零空间X表示的是列向量的线性组合。X的物理意义表示的每个节点的电势值[potential at the node ],那么 x 2 − x 1 x_2-x_1 x2−x1表示的是edge1的电势差,那什么时候可以保证电势差为0呢,那只有每个结点的电势相同即可,所以可以得到解如下,所以 dim N ( A ) = 1 , R a n k ( A ) = 4 − 1 = 3 N(A)=1,Rank(A)=4-1=3 N(A)=1,Rank(A)=4−1=3
X = c [ 1 1 1 1 ] (6) X=c\begin{bmatrix}1\\\\1\\\\1\\\\1\end{bmatrix}\tag{6} X=c 1111 (6)
2. A T Y = 0 A^TY=0 ATY=0
- 矩阵表示形式如下, dim N ( A T ) = m − r ( A T ) = 5 − 3 = 2 N(A^T)=m-r(A^T)=5-3=2 N(AT)=m−r(AT)=5−3=2,物理上表示 y 1 … y 5 y_1\dots y_5 y1…y5为每个边edge上的电流大小:
A T Y = [ − 1 0 − 1 − 1 0 1 − 1 0 0 0 0 1 1 0 − 1 0 0 0 1 1 ] [ y 1 y 2 y 3 y 4 y 5 ] = [ − y 1 − y 3 − y 4 y 1 − y 2 y 2 + y 3 − y 5 y 4 + y 5 ] = [ 0 0 0 0 ] (7) A^TY=\begin{bmatrix}-1& 0& -1& -1&0\\\\ 1& -1& 0& 0&0\\\\ 0& 1& 1& 0&-1\\\\ 0& 0& 0& 1&1 \end{bmatrix}\begin{bmatrix}y_1\\\\y_2\\\\y_3\\\\y_4\\\\y_5\end{bmatrix}=\begin{bmatrix}-y_1-y_3-y_4\\\\y_1-y_2\\\\y_2+y_3-y_5\\\\y_4+y_5\end{bmatrix}=\begin{bmatrix}0\\\\0\\\\0\\\\0\end{bmatrix}\tag{7} ATY= −11000−110−1010−100100−11 y1y2y3y4y5 = −y1−y3−y4y1−y2y2+y3−y5y4+y5 = 0000 (7) - 函数物理意义:
− y 1 − y 3 − y 4 = 0 ⇒ y 1 + y 3 = − y 4 (8) -y_1-y_3-y_4=0\Rightarrow y_1+y_3=-y_4\tag{8} −y1−y3−y4=0⇒y1+y3=−y4(8)
表示的对于节点1来说,流入的电流为 − y 4 -y_4 −y4,流出的电流为 y 1 + y 3 y_1+y_3 y1+y3,所以可以得到基尔霍夫电流定律,流入节点电流和等于流出节点的电流和
- N( A T A^T AT)=m-r=5-3=2,我们看由1,2,3组成的节点环1,如果想成为闭环,可以得到如下:
y 1 + y 2 − y 3 = 0 (9) y_1+y_2-y_3=0\tag{9} y1+y2−y3=0(9) - A T Y = 0 A^TY=0 ATY=0的特解中当 y 1 = 1 , y 2 = 1 , y 3 = − 1 , y 4 = 0 , y 5 = 0 y_1=1,y_2=1,y_3=-1,y_4=0,y_5=0 y1=1,y2=1,y3=−1,y4=0,y5=0,表示 L o o p 123 Loop_{123} Loop123
- A T Y = 0 A^TY=0 ATY=0的特解中当 y 1 = 0 , y 2 = 0 , y 3 = 1 , y 4 = − 1 , y 5 = 1 y_1=0,y_2=0,y_3=1,y_4=-1,y_5=1 y1=0,y2=0,y3=1,y4=−1,y5=1,表示 L o o p 345 Loop_{345} Loop345
- A T Y = 0 A^TY=0 ATY=0的特解中当 y 1 = 1 , y 2 = 1 , y 3 = 0 , y 4 = − 1 , y 5 = 1 y_1=1,y_2=1,y_3=0,y_4=-1,y_5=1 y1=1,y2=1,y3=0,y4=−1,y5=1,表示 L o o p 1234 Loop_{1234} Loop1234
- 线性无关组表示不形成环时的边组合。这样可以看出来矩阵和图的对应关系。很直观。
- 矩阵的相关就是表示图形中的环,矩阵的不相关就是表示不成环,形成树。
- 矩阵的秩和图之间的关系
- D i m { N ( A T ) } Dim \{N(A^T)\} Dim{N(AT)}=#loops;
- m=#edges;
- r=n-1,节点的数量-1
D i m { N ( A T ) } = m − r (10) Dim \{N(A^T)\}= m-r\tag{10} Dim{N(AT)}=m−r(10)
综上所述,可以得到如下
# n o d e s − # e d g e s + # L o o p s = 1 (11) \#{nodes}-\#{edges}+\#{Loops}=1\tag{11} #nodes−#edges+#Loops=1(11)
3. A X = 0 AX=0 AX=0和 A T Y = 0 A^TY=0 ATY=0的关系
-
在 A X = 0 AX=0 AX=0 的方程中,我们定义每个边两端的电势差E
E = A X (12) E=AX\tag{12} E=AX(12) -
在 A T Y = 0 A^TY=0 ATY=0的方程中,我们定义每个边中的电流为y,通过欧姆定理,定义电势差E,和常数C,电流Y之间的关系如下:
Y = C E (13) Y=CE\tag{13} Y=CE(13) -
电流方程如下:
A T Y = 0 (14) A^TY=0\tag{14} ATY=0(14) -
将方程12,13 代入到方程14中可得如下:
A T C A X = 0 (15) A^TCAX=0\tag{15} ATCAX=0(15) -
我们定义外部电流从原来的0变为f,整理公式15可得平衡方程,并且 A T C A A^TCA ATCA对称:
A T C A X = f (16) A^TCAX=f\tag{16} ATCAX=f(16)