ieee33.rar

ieee33.rar

clear all;
Bus=[100   ;
2100 60;
39040;
412080;
56030;
66020;
7 200 100 ;
8200  100 ;
9 60 20;
1060 20;
114530;
1260 35;
1360 35;
1412080;
1560 10;
1660 20;
1760 20;
1890 40;
1990 40;
2090 40;
2190 40;
2290 40; 
2390 40; 
24420 200;
25420 200 ;
2660 25;
2760 25;
2860 20;
29120 70;
30200600 ;
31150 70;
32210 100 ;
336040;
];
Branch=[1 1 2 0.09220.0407;
  2 2 3 0.49300.2511;
  3 3 4 0.36600.1864;
4 4 5 0.38110.1941;
5 5 6 0.81900.7070;
6 6 7 0.18720.6188;
7 7 8 0.71440.2351;
8 8 9 1.03000.7400;
9 9 101.04400.7400;
1010110.19660.065;
1111120.37440.1238;
1212131.46801.1550;
1313140.54160.7129;
1414150.59100.5260;
1515160.74630.5450;
1616171.28901.7210;
1717180.73200.5740;
182 190.16400.1565;%这不1到19的吗
1919201.50421.3554;
2020210.40950.4784;
2121220.70890.9373;
223 230.45120.3083;%2到23？
2323240.89800.7091;
2424250.89600.7011;
256 260.20300.1034;%5到26？
2626270.28420.1447;
2727281.05900.9337;
2828290.80420.7006;
2929300.50750.2585;
3030310.97440.9630;
3131320.31050.3619;
3232330.34100.5302;
 ];
[busnumrow]=size（Bus）;%返回33行3列
[branchnumrow]=size（Branch）;
soubus=Branch（:2）;%返回支路矩阵第二列所有元素
mobus=Branch（:3）;
Vbus=ones（busnum1）;%产生一个33行1列的全1阵
Vbus（:1）=12.66;%Vbus阵中所有元素全变为12.66，即假定所有母线电压为12.66
Vbus1=Vbus;
Ploss=zeros（busnum1）;%产生一个33行1列的全0阵
Qloss=zeros（busnum1）;
e=1;
k=0;
Branch1=Branch;
n=1;
while ~isempty（Branch1）
m=1;
[srow]=size（Branch1）;
while s>0
t=find（Branch1（:2）== Branch1（s3））;%寻找所有行第二列等于33的位置，返回‘行’
if isempty（t）
T1（n:）= Branch1（s:）;%把32行的所有元素赋给第n行，第一次n=1
n=n+1;
else
T2（m:）= Branch1（s:）;
m=m+1;
end
s=s-1;%313029，...0
end
Branch1=T2;
T2=[];
end
while e>1.0e-05
P=zeros（busnum1）;%产生一个33行1列的全0阵
Q=zeros（busnum1）;
for s=1:branchnum %s=1:32
i=T1（s2）;
j=T1（s3）;
R=T1（s4）;
X=T1（s5）;
Pload=Bus（j2）;
Qload=Bus（j3）;
II=（（Pload+P（j））^2+（Qload+Q（j））^2）/（Vbus（j）^2*1000）;
Ploss（ij）=II*R;
Qloss（ij）=II*X;
P（ij）=Pload+Ploss（ij）+P（j）;
Q（ij）=Qload+Ploss（ij）+Q（j）;
P（i）=P（i）+P（ij）;
Q（i）=Q（i）+Q（ij）;
end
for s=branchnum:-1:1%s=323130...1
i=T1（s2）;
j=T1（s3）;
R=T1（s4）;
X=T1（s5）;
Vbus（j）=（Vbus（i）-（P（ij）*R+Q（ij）*X）/（Vbus（i）*1000））^2+（（P（ij）*X-Q（ij）*R）/（Vbus（i）*1000））^2;
Vbus（j）=sqrt（Vbus（j））;
end
e=max（abs（Vbus1-Vbus））;
Vbus1=Vbus;
k=k+1;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      27788  2018-10-27 23:44  ieee33.slx

     文件       2619  2018-10-29 19:46  mieee33.m

----------- ---------  ---------- -----  ----

                30407                    2