admain1.m

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

obsline=csvread('obsline.csv');

nfrom=obsline(1,1);
nto=obsline(1,2);
ni=locLN(nfrom, nto, line);
mui=PLineX(ni,1);
smgi=sqrt(PLineX(ni,2));
xi=-160;
for j=1:400
    PCL(j,1)=xi+(j-1);
    PCL(j,2)=CalCdf1(PCL(j,1), mui, smgi, PLineC(ni, 3), PLineC(ni,4), PLineC(ni,5), PLineC(ni,6), PLineC(ni, 7), PLineC(ni,8), PLineC(ni,9));
    PCL(j,3)=CalPdf1(PCL(j,1), mui, smgi, PLineC(ni, 3), PLineC(ni,4), PLineC(ni,5), PLineC(ni,6), PLineC(ni, 7), PLineC(ni,8), PLineC(ni,9));
end    
csvwrite('Npc3dfline107_108.csv',PCL);