本文共 2479 字，大约阅读时间需要 8 分钟。

节点排序：1至m为PQ节点，m+1至n-1为PV节点，n为平衡节点。

---

清空变量；

清空控制标志；初始化数据矩阵Y，设置为零矩阵；

对于i从1到n：

如果X(i,2)不等于0：p = X(i,1);Y(p,p) = 1./X(i,2);

对于i从1到nl：

如果B1(i,6)等于0：p = B1(i,1);q = B1(i,2);否则：p = B1(i,2);q = B1(i,1);Y(p,q) = Y(p,q) - 1/(B1(i,3)*B1(i,5));Y(q,p) = Y(p,q);Y(q,q) = Y(q,q) + 1/(B1(i,3)*B1(i,5)^2) + B1(i,4)/2;Y(p,p) = Y(p,p) + 1/B1(i,3) + B1(i,4)/2;

输出导纳矩阵Y；

---

G = 实数部分(Y)；

B = 虚数部分(Y)；

对于i从1到n：

e(i) = 实数部分(B2(i,3));f(i) = 虚数部分(B2(i,3));V(i) = B2(i,4)；

对于i从1到n：

S(i) = B2(i,1) - B2(i,2)；

Ps = 实数部分(S)；

Qs = 虚数部分(S)；

初始化数值和矩阵：

num = 0；control = 1；J = 零矩阵，大小为(2n-2, 2n-2)；b = 零向量，大小为(2n-2,1)；

while control ≠ 0：

对于i从1到n-1：计算C和D：C = G(i,j)*e(j) - B(i,j)*f(j)；D = G(i,j)*f(j) + B(i,j)*e(j)；对于j从1到n：C = C + G(i,j)*e(j) - B(i,j)*f(j)；D = D + G(i,j)*f(j) + B(i,j)*e(j)；P(i) = e(i)*C + f(i)*D；Q(i) = f(i)*C - e(i)*D；V2 = e(i)^2 + f(i)^2；如果i ≤ m：DP = Ps(i) - P(i)；DQ = Qs(i) - Q(i)；对于j从1到n-1：如果i == j：X1 = -C - G(i,i)*e(i) - B(i,i)*f(i)；X2 = -D + B(i,i)*e(i) - G(i,i)*f(i)；X3 = D + B(i,i)*e(i) - G(i,i)*f(i)；X4 = -C + G(i,i)*e(i) + B(i,i)*f(i)；否则：X1 = -G(i,j)*e(i) - B(i,j)*f(i)；X2 = B(i,j)*e(i) - G(i,j)*f(i)；X3 = X2；X4 = -X1；p = 2i-1；q = 2j-1；u = p+1；v = q+1；J(p,q) = X1；J(p,v) = X2；J(u,q) = X3；J(u,v) = X4；b(p,1) = DP；b(u,1) = DQ；如果i > m：DV = V(i)^2 - V2；对于j从1到n-1：如果i == j：X1 = -C - G(i,i)*e(i) - B(i,i)*f(i)；X2 = -D + B(i,i)e(i) - G(i,i)f(i)；

X5 = -2e(i)；

X6 = -2f(i)；否则：X1 = -G(i,j)*e(i) - B(i,j)*f(i)；X2 = B(i,j)*e(i) - G(i,j)*f(i)；X5 = 0；X6 = 0；p = 2i-1；q = 2j-1；u = p+1；v = q+1；J(p,q) = X1；J(p,v) = X2；J(u,q) = X5；J(u,v) = X6；b(p,1) = DP；b(u,1) = DV；

J = -J；  x = J \ b；  对于i从1到n-1：      e(i) = e(i) + x(2i-1)；      f(i) = f(i) + x(2i)；  ABS = 绝对值(b)；  MAX = ABS的最大值；  如果MAX ≤ pr：      control = 0；  否则：      num = num + 1；

对于i从1到num：

Dy(num,i) = sqrt(e(i)^2 + f(i)^2)；

输出迭代次数；

对于i从1到n：

V(i) = sqrt(e(i)^2 + f(i)^2)；O(i) = atan(f(i)/e(i)) * 180/π；

输出各节点的电压大小和电压相角；

对于p从1到n：

C(p) = 0；对于q从1到n：C(p) = C(p) + Y(p,q)*E(q)；S(p) = E(p) * conj(C(p))；

输出各节点的功率；

对于i从1到nl：

如果B1(i,6)等于0：p = B1(i,1)；q = B1(i,2)；否则：p = B1(i,2)；q = B1(i,1)；Si(p,q) = E(p) * (conj(E(p)) * conj(B1(i,4)) + (conj(E(p)/B1(i,5)) - conj(E(q))) * conj(1/(B1(i,3)*B1(i,5))))；

对于i从1到nl：

如果B1(i,6)等于0：p = B1(i,1)；q = B1(i,2)；否则：p = B1(i,2)；q = B1(i,1)；Sj(q,p) = E(q) * (conj(E(q)) * conj(B1(i,4)) + (conj(E(q)/B1(i,5)) - conj(E(p))) * conj(1/(B1(i,3)*B1(i,5))))；

输出各条支路的首端功率和末端功率；

输出各条支路的功率损耗；

输出各节点的功率；

输出各条支路的功率损耗；

输出各节点的功率；

生成迭代次数与电压迭代曲线的图形；