plotpmf.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=-200;
    for j=1:38
        PCL(i,j)=xi+(j-1)*10;
        PPL(i,j)=CalCdf(PCL(i,j), mui, smgi, PLineC(ni, 3), PLineC(ni,4), PLineC(ni,5), PLineC(ni,6));
    end    
    % figure(i);
    % plot(PCL(i,:), PPL(i,:));
end   

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