本文共 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))))；输出各条支路的首端功率和末端功率；
输出各条支路的功率损耗；
输出各节点的功率；
输出各条支路的功率损耗；
输出各节点的功率；
生成迭代次数与电压迭代曲线的图形；
转载地址：http://htgfk.baihongyu.com/