calc8.m

来自「可以计算电力系统概率潮流」· M 代码 · 共 17 行

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

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

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;
end

return