calc1.m

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

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

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

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;   
    pci(1,7)=-u7/(sgm^7)+21*u5/(sgm^5)-105*u3/sgm^3;
    pci(1,8)=u8/(sgm^8)-28*u6/sgm^6+210*u4/sgm^4-315;
    pci(1,9)=-u9/(sgm^9)+36*u7/sgm^7-378*u5/sgm^5+1260*u3/sgm^3;
end

return