calu.m

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

function pui=CalU(x1, x2, x3, x4, x5, x6)

% calculate Central Moments
pui=zeros(1,6);

pui(1,1)=0;
pui(1,2)=x2;
pui(1,3)=x3;
pui(1,4)=x4+3*(pui(1,2)^2);
pui(1,5)=x5+10*pui(1,2)*pui(1,3);
pui(1,6)=x6+15*pui(1,2)*pui(1,4)+10*(pui(1,3)^2)-30*(pui(1,2)^3);


return