G11(i,li(i))=-H(li(i),lj(i)); G11(i,lj(i))=H(li(i),lj(i)); G12(i,li(i))=0.01*2*branch(i,14)/bus(li(i),8)-N(li(i),lj(i)); %%化为标幺值 G12(i,lj(i))=N(li(i),lj(i)); G21(i,li(i))=-M(li(i),lj(i)); G21(i,lj(i))=M(li(i),lj(i)); G22(i,li(i))=0.01*2*branch(i,15)/bus(li(i),8)-H(li(i),lj(i)); %%化为标幺值 G22(i,lj(i))=H(li(i),lj(i));
时间: 2023-10-28 12:05:54 浏览: 28
这部分代码是在构建支路潮流的灵敏度矩阵中的各个分块矩阵,具体包括:
- `G11` 矩阵表示支路有功对角度的倒数,根据潮流方程和导纳矩阵的定义,可得支路有功对角度的导数为 `-Re(Y(i,j))`,因此 `G11(i,li(i))=-H(li(i),lj(i))` 和 `G11(i,lj(i))=H(li(i),lj(i))`,其中 `H` 是节点导纳矩阵的实部;
- `G12` 矩阵表示支路有功对电压的导数,根据潮流方程和导纳矩阵的定义,可得支路有功对电压的导数为 `Re(Y(i,j))*V(i)`,其中 `V(i)` 表示节点电压幅值,因此 `G12(i,li(i))=0.01*2*branch(i,14)/bus(li(i),8)-N(li(i),lj(i))` 和 `G12(i,lj(i))=N(li(i),lj(i))`,其中 `N` 是节点导纳矩阵的虚部,`branch(i,14)` 表示第 `i` 条支路的电阻,`bus(li(i),8)` 表示第 `li(i)` 个节点的基准电压;
- `G21` 矩阵表示支路无功对角度的导数,根据潮流方程和导纳矩阵的定义,可得支路无功对角度的导数为 `-Im(Y(i,j))*V(i)`,因此 `G21(i,li(i))=-M(li(i),lj(i))` 和 `G21(i,lj(i))=M(li(i),lj(i))`,其中 `M` 是节点导纳矩阵的虚部;
- `G22` 矩阵表示支路无功对电压的导数,根据潮流方程和导纳矩阵的定义,可得支路无功对电压的导数为 `-Im(Y(i,j))*V(i)^2`,因此 `G22(i,li(i))=0.01*2*branch(i,15)/bus(li(i),8)-H(li(i),lj(i))` 和 `G22(i,lj(i))=H(li(i),lj(i))`,其中 `branch(i,15)` 表示第 `i` 条支路的电抗。