www.pudn.com > powerflow.zip > aco2.m, change:2012-06-07,size:2450b


clear; 
clc; 
branch=xlsread('IEEE33.xls');%初始网络结构数据 
n=1; 
times=zeros(32,11); 
timeslost=zeros(1,2); 
for n=1:10000 
    P1=normrnd(100,15); 
    P2=normrnd(90,13.5); 
    P3=normrnd(120,18); 
    P4=normrnd(60,9); 
    P5=normrnd(60,9); 
    P6=normrnd(200,30); 
    P7=normrnd(200,30); 
    P8=normrnd(60,9); 
    P9=normrnd(60,9); 
    P10=normrnd(45,6.75); 
    P11=normrnd(60,9); 
    P12=normrnd(60,9); 
    P13=normrnd(120,18); 
    P14=normrnd(60,9); 
    P15=normrnd(60,9); 
    P16=normrnd(60,9); 
    P17=normrnd(90,13.5); 
    P18=normrnd(90,13.5); 
    P19=normrnd(90,13.5); 
    P20=normrnd(90,13.5); 
    P21=normrnd(90,13.5); 
    P22=normrnd(90,13.5); 
    P23=normrnd(420,63); 
    P24=normrnd(420,63); 
    P25=normrnd(60,9); 
    P26=normrnd(60,9); 
    P27=normrnd(60,9); 
    P28=normrnd(120,18); 
    P29=normrnd(200,30); 
    P30=normrnd(150,22.5); 
    P31=normrnd(210,31.5); 
    P32=normrnd(60,9); 
    Q1=normrnd(60,9); 
    Q2=normrnd(40,6); 
    Q3=normrnd(80,12); 
    Q4=normrnd(30,4.5); 
    Q5=normrnd(20,3); 
    Q6=normrnd(100,15); 
    Q7=normrnd(100,15); 
    Q8=normrnd(20,3); 
    Q9=normrnd(20,3); 
    Q10=normrnd(30,4.5); 
    Q11=normrnd(35,5.25); 
    Q12=normrnd(35,5.25); 
    Q13=normrnd(80,12); 
    Q14=normrnd(10,1.5); 
    Q15=normrnd(20,3); 
    Q16=normrnd(20,3); 
    Q17=normrnd(40,6); 
    Q18=normrnd(40,6); 
    Q19=normrnd(40,6); 
    Q20=normrnd(40,6); 
    Q21=normrnd(40,6); 
    Q22=normrnd(50,7.5); 
    Q23=normrnd(200,30); 
    Q24=normrnd(200,30); 
    Q25=normrnd(25,3.75); 
    Q26=normrnd(25,3.75); 
    Q27=normrnd(20,3); 
    Q28=normrnd(70,10.5); 
    Q29=normrnd(600,90); 
    Q30=normrnd(70,10.5); 
    Q31=normrnd(100,15); 
    Q32=normrnd(40,6); 
    branch(:,6)=[P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 P11 P12 P13 P14 P15 P16 P17 P18 P19 P20 P21 P22 P23 P24 P25 P26 P27 P28 P29 P30 P31 P32]; 
    branch(:,7)=[Q1 Q2 Q3 Q4 Q5 Q6 Q7 Q8 Q9 Q10 Q11 Q12 Q13 Q14 Q15 Q16 Q17 Q18 Q19 Q20 Q21 Q22 Q23 Q24 Q25 Q26 Q27 Q28 Q29 Q30 Q31 Q32]; 
    branch0=branch;   %转移负荷用 
    nbranch=length(branch);%初始网络支路数 
    maxnext=4;             %最多下接支路数 
    powerflow;             %初始网络潮流计算 
end 
prob=zeros(32,11); 
problost=zeros(1,2); 
prob=times/10000 
problost=timeslost/10000 
xlswrite('不含光伏电源的结果.xls',prob,'','B1'); 
xlswrite('不含光伏电源的结果.xls',problost,'','M1') 
S0=S;                  %初始网络潮流 
Plose0=Plose;          %初始网络损耗