function [V]=disflowtry()% 算例為 IEEE33 ，只計算初始網絡的節點電壓，程序不完備，目的是便于閱讀BranchData = [    0   1   0.0922  0.0470;    1   2   0.4930  0.2511;    2   3   0.3660  0.1864;    3   4   0.3811  0.1941;    4   5   0.8190  0.7070;    5   6   0.1872  0.6188;    6   7   0.7114  0.2351;    7   8   1.0300  0.7400;    8   9   1.0440  0.7400;    9   10  0.1966  0.0650;    10  11  0.3744  0.1238;    11  12  1.4680  1.1550;    12  13  0.5416  0.7129;    13  14  0.5910  0.5260;    14  15  0.7463  0.5450;    15  16  1.2890  1.7210;    16  17  0.7320  0.5740;    1   18  0.1640  0.1565;    18  19  1.5042  1.3554;    19  20  0.4095  0.4784;    20  21  0.7089  0.9373;    2   22  0.4512  0.3083;    22  23  0.8980  0.7091;    23  24  0.8960  0.7011;    5   25  0.2030  0.1034;    25  26  0.2842  0.1447;    26  27  1.0590  0.9337;    27  28  0.8042  0.7006;    28  29  0.5075  0.2585;    29  30  0.9744  0.9630;    30  31  0.3105  0.3619;    31  32  0.3410  0.5302;];                                     % 支路，阻抗NodeData = [    1   100.00  60.00;    2   90.00   40.00;    3   120.00  80.00;    4   60.00   30.00;    5   60.00   20.00;    6   200.00  100.00;    7   200.00  100.00;    8   60.00   20.00;    9   60.00   20.00;    10  45.00   30.00;    11  60.00   35.00;    12  60.00   35.00;    13  120.00  80.00;    14  60.00   10.00;    15  60.00   20.00;    16  60.00   20.00;    17  90.00   40.00;    18  90.00   40.00;    19  90.00   40.00;    20  90.00   40.00;    21  90.00   40.00;    22  90.00   50.00;    23  420.00  200.00;    24  420.00  200.00;    25  60.00   25.00;    26  60.00   25.00;    27  60.00   20.00;    28  120.00  70.00;    29  200.00  600.00;    30  150.00  70.00;    31  210.00  100.00;    32  60.00   40.00;];                                      % 節點，負荷UB = 12.66;                             % 電壓基準 kVSB = 10;                                % 功率基準 MVAZB = UB^2/SB;                           % 阻抗基準 ohmBranchData(:,[3,4]) = BranchData(:,[3,4]) / ZB;     % 阻抗標幺化NodeData(:,[2,3]) = NodeData(:,[2,3]) / SB / 1000;% 功率標幺化BranchData(:,[1,2]) = BranchData(:,[1,2]) + 1;        % 節點從1開始NN = 33;                                % 節點數A0 = zeros(NN);for n = 1:32    A0(BranchData(n,1),BranchData(n,2)) = 1;end                                     % 形成 A0A0T = A0';                              % 形成 A0 的轉置S = [0;-NodeData(:,2) - i*NodeData(:,3)];       % 形成 SZL = [0;BranchData(:,3) + i*BranchData(:,4)];       % 形成 ZLV = ones(NN,1);V(1) = 1.05;                            % 各個節點電壓賦初值IL(NN,1) = -conj(S(NN) / V(NN));          % 最末支路電流賦初值Delta = 1;                              % 收斂判據賦初值TempV = V;                              % 賦初值，用于記憶上次迭代結果while Delta > 1e-8    IN = conj(S ./ V);                    % 節點注入電流        for n = 1:NN-1        IL(NN-n) = A0(NN-n,NN-n+1:end) * IL(NN-n+1:end) - IN(NN-n);    end                                 % 電流回代過程        for n = 2:NN        V(n) = A0T(n,1:n-1) * V(1:n-1) - ZL(n) * IL(n);    end                                 % 電壓前推過程        Delta = max(abs(V-TempV));          % 更新收斂判據    TempV = V;                          % 記憶迭代結果end