IEEE标准14节点PQ分解法MATLAB程序

% 程序功能:14节点潮流计算PQ分解法
clc;


% NODE INFO 
% NODE TYPE    PQ:1     PV:2    PH:3   所有发电机节点均被视为PV节点，节点1为平衡节点
NODE=[
% TYPE           V        ANGEL      PG         QG        PL       QL     SHUNT C（节点并联电容）
        3    1.0600         0         0         0         0         0         0
        2    1.0450         0   40.0000         0   21.7000   12.7000         0
        2    1.0100         0         0         0   94.2000   19.0000         0
        1    1.0000         0         0         0   47.8000   -3.9000         0
        1    1.0000         0         0         0    7.6000    1.6000         0
        2    1.0700         0         0         0   11.2000    7.5000         0
        1    1.0000         0         0         0         0         0         0
        2    1.0900         0         0         0         0         0         0
        1    1.0000         0         0         0   29.5000   16.6000    0.1900
        1    1.0000         0         0         0    9.0000    5.8000         0
        1    1.0000         0         0         0    3.5000    1.8000         0
        1    1.0000         0         0         0    6.1000    1.6000         0
        1    1.0000         0         0         0   13.5000    5.8000         0
        1    1.0000         0         0         0   14.9000    5.0000         0
 ] 



% 功率为有名值，基准功率为100MVA
% 线路对地电纳参数为B，不是B/2
%  BRANCH DATA  IEEE 14 NODE
%     i    j               R            X            B             K（变比）
BRANCH=[
      2     1           .01938       .05917       .0264*2        1.000;
      3     2           .04699       .19797       .0219*2        1.000;
      4     2           .05811       .17632       .0187*2        1.000;
      5     1           .05403       .22304       .0246*2        1.000;
      5     2           .05695       .17388       .0170*2        1.000;
      4     3           .06701       .17103       .0173*2        1.000;
      5     4           .01335       .04211       .0064*2        1.000;
      6     5           .00000       .25202       .0000          0.932;    %变压器支路
      7     4           .00000       .20912       .0000          0.978;    %变压器支路
      8     7           .00000       .17615       .0000          1.000;
      9     4           .00000       .55618       .0000          0.969;    %变压器支路
      9     7           .00000       .11001       .0000          1.000;
     10     9           .03181       .08450       .0000          1.000;
     11     6           .09498       .19890       .0000          1.000;
     12     6           .12291       .15581       .0000          1.000;
     13     6           .06615       .13027       .0000          1.000;
     14     9           .12711       .27038       .0000          1.000;
     11    10           .08205       .19207       .0000          1.000;
     13    12           .22092       .19988       .0000          1.000;
     14    13           .17093       .34802       .0000          1.000;