calc7.m

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

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

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

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

return