33节点潮流编程

程序编写Institute of Electrical and Electronics Engineers (IEEE) 美国电气和电子工程师协会。IEEE-33节点标准系统潮流程序的设计采用牛顿拉夫逊法潮流分布的计算原理，运用电力系统分析关于系统网络计算的相关知识，在MATLAB里编程实现潮流计算的验证。PQ节点：这类节点的有功功率P和无功功率Q是给定的，节点电压和相位（V，δ）是待求量。通常变电所都是这一类型的节点。

clc;
clear all;
close all;
B=[120.4930+0.2511j010;230.3660+0.1864j010;
   340.3811+0.1941j010;450.8190+0.7070j010;
   560.1872+0.6188j010;670.7114+0.2351j010;
   781.0300+0.7400j010;891.0400+0.7400j010;
   9100.1966+0.0650j010;10110.3744+0.1238j010;
   11121.4680+1.1550j010;12130.5416+0.7129j010;
   13140.5910+0.5260j010;14150.7463+0.5450j010;
   15161.2890+1.7210j010;16170.3720+0.5740j010;
   1180.1640+0.1565j010;18191.5042+1.3554j010;
   19200.4095+0.4784j010;20210.7089+0.9373j010;
   2220.4512+0.3083j010;22230.8980+0.7091j010;
   23240.8960+0.7011j010;5250.2030+0.1034j010;
   25260.2842+0.1447j010;26271.0590+0.9377j010;
   27280.8042+0.7006j010;28290.5075+0.2585j010;
   29300.9744+0.9630j010;30310.3105+0.3619j010; 
   31320.3410+0.5312j010;1330.0922+0.0407j010;];
X=[10;20;30;40;50;60;70;80;90;100;110;120;130;140;150;160;170;
  180;190;200;210;220;230;240;250;260;270;280;290;300;310;320;330;];
U=ones（331）;
U=12.66*U;
m=zeros（331）;
S=[9040;12080;6030;6020;200100;200100; 6020; 6020;4530; 6035; 6035;12080;6010; 6020; 6020;9040; 9040; 9040;9040;
9040; 9050; 420200; 420200; 6025;6025;6020; 12070; 200600;15070;210100;6040;10060;];
Y=zeros（33）;
for i=1:33
    if X（i2）~=0;
        p=X（i1）;
        Y（pp）=1./X（i2）;
    end
end
for i=1:32
 if B（i6）==0
    p=B（i1）;q=B（i2）;
    else
     p=B（i2）;q=B（i1）;
 end
Y（pq）=Y（pq）-1./B（i3）*B（i5）;
Y（qp）=Y（pq）;
Y