电网中直流潮流Distribution Factor的计算，bus间有多条输电线

电网中DC power flow的计算依赖与distribution factor的计算，而DF的计算是由电网结构决定的，电网主要有bus与line组成的，以一个简单的3-bus系统为例：Bus 数量：N
Line 数量：L
line1: bus1-bus2
line2: bus1-bus3
line3: bus2-bus3
incidence matrix A: 表示line的起始信息，行表示line，相关列中1是表示起始bus，-1表示终点bus，L×N-matrix。
$A=\left[\begin{array}{ccc}1& -1& 0\\ 1& 0& -1\\ 0& 1& -1\end{array}\right]$

admittance matrix Y：N×N-matrix 导纳矩阵
$\mathbf{Y}={\mathbf{A}}^T\cdot {\mathbf{Y}}_d\cdot \mathbf{A}$
${\mathbf{Y}}_d$ a L×L-diagonal matrix with the line admittances on the diagonal (i.e., the primitive admittance matrix).
阻抗的实数部分是电阻，虚数部分是电抗。
导纳是电导和电纳的复合形式，忽略了电阻的话，电纳就是电抗的倒数。
susceptance matrix B:电纳矩阵
DC power flow中，忽略了电阻，主要关注的是线路的电抗。
${\mathbf{p}}_L={\mathbf{B}}_d\cdot \mathbf{A}\cdot {\mathbf{\delta }}_N$

${\mathbf{B}}_d$是LxL的对角矩阵，具体形式为：
${\mathbf{B}}_d=\left[\begin{array}{ccc}5& 0& 0\\ 0& 10& 0\\ 0& 0& 5\end{array}\right]$

line L上的flow为：$P_L=B_L\left({\delta }_N-{\delta }_Q\right)=\frac{1}{X_{NQ}}({\delta }_N-{\delta }_Q)$

从bus N中出去的flow总和等于bus B的节点注入功率，
$\begin{array}{rl}{P}_N& =\sum _Q{B}_L\left({\delta }_N-{\delta }_Q\right)\\ {\mathbf{p}}_N& ={\mathbf{A}}^T\cdot {\mathbf{B}}_d\cdot \mathbf{A}\cdot {\mathbf{\delta }}_N\end{array}$

${\mathbf{p}}_L=\left(\left({\mathbf{B}}_d\cdot \mathbf{A}\right)\cdot {\left({\mathbf{A}}^T\cdot {\mathbf{B}}_d\cdot \mathbf{A}\right)}^{-1}\right)\cdot {\mathbf{p}}_N$

${\mathbf{P}\mathbf{T}\mathbf{D}\mathbf{F}}^{L\times N}=\left({\mathbf{B}}_d\cdot \mathbf{A}\right)\cdot {\left({\mathbf{A}}^T\cdot {\mathbf{B}}_d\cdot \mathbf{A}\right)}^{-1}$
需要注意的是，$\left({\mathbf{B}}_d\cdot \mathbf{A}\right)$需要去除reference bus的列, $\left({\mathbf{A}}^T\cdot {\mathbf{B}}_d\cdot \mathbf{A}\right)$是非线性独立的，需要去掉reference bus的行和列。由此得到shift（distribution） factor。
bus 3作为reference，求得：
$PTDF=\left[\begin{array}{ccc}0.2& -0.4& 0\\ 0.8& 0.4& 0\\ 0.2& 0.6& 0\end{array}\right]$

bus之间不止一条输电线的case：
line1: bus1-bus2
line2: bus1-bus3
line3: bus1-bus3
line4: bus2-bus3
${\mathbf{B}}_d=\left[\begin{array}{ccc}5& 0& 0\\ 0& 10& 0\\ 0& 0& 5\end{array}\right]$
$PTDF=\left[\begin{array}{ccc}0.2& -0.4& 0\\ 0.4& 0.2& 0\\ 0.4& 0.2& 0\\ 0.2& 0.6& 0\end{array}\right]$