calc6.m

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

function pci=CalC1(sgm, u1, u2, u3, u4, u5, u6)

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

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