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