calc.m

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

function pci=CalC(x1, x2, x3, x4, x5, x6)

% calculate Gram Charles index 
pci=zeros(1,6);

u1=0;
u2=x2;
u3=x3;
u4=x4+3*(u2^2);
u5=x5+10*u2*u3;
u6=x6+15*u2*u4+10*(u3^2)-30*(u2^3);

sgm=sqrt(x2);

if (sgm~=0)
    pci(1,1)=0;
    pci(1,2)=0;
    pci(1,3)=-u3/(sgm^3);   
    pci(1,4)=u4/(sgm^4)-3;
    pci(1,5)=-u5/(sgm^5)+10*u3/(sgm^3);
    pci(1,6)=u6/(sgm^6)-15*u4/(sgm^4)+30;   
end

return