admain.m

可以计算电力系统概率潮流

obsline=csvread('obsline.csv');
nobl=length(obsline(:,1));

for i=1:nobl
    nfrom=obsline(i,1);
    nto=obsline(i,2);
    ni=locLN(nfrom, nto, line);
    mui=PLineX(ni,1);
    smgi=sqrt(PLineX(ni,2));
    xi=-2600;
    for j=1:19
        PCL(i,j)=xi+(j-1)*100;
        PPL(i,j)=CalCdf1(PCL(i,j), mui, smgi, PLineC(ni, 3), PLineC(ni,4), PLineC(ni,5), PLineC(ni,6), PLineC(ni, 7), PLineC(ni,8), PLineC(ni,9));
        PDL(i,j)=CalPdf1(PCL(i,j), mui, smgi, PLineC(ni, 3), PLineC(ni,4), PLineC(ni,5), PLineC(ni,6), PLineC(ni, 7), PLineC(ni,8), PLineC(ni,9));
    end    
    figure(i);
    plot(PCL(i,:), PPL(i,:))
    figure(2*i);
    plot(PCL(i,:), PDL(i,:))
end   

PCL1=PCL(1,2:12);
for i=1:11
    PL(1,i)=PPL(1,i+1)-PPL(1,i);
end
figure(2)
plot(PCL1(1,:),PL(1,:));