calu6.m

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

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

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

afa(1,1)=x1;
afa(1,2)=x2+x1^2;
afa(1,3)=x3+3*x1*x2+x1^3;
afa(1,4)=x4+3*x2^2+4*x1*x3+6*x1^2*x2+x1^4;
afa(1,5)=x5+5*x1*x4+10*x2*x3+10*x1^2*x3+15*x1*x2^2+10*x1^3*x2+x1^5;
afa(1,6)=x6+10*x3^2+15*x2*x4+6*x1*x5+15*x2^3+15*x1^2*x4+60*x1*x2*x3+20*x1^3*x3+45*x1^2*x2^2+15*x1^4*x2+x1^6;

pui=zeros(1,6);
pui(1,1)=0;
pui(1,2)=afa(1,2)-afa(1,1)^2;
pui(1,3)=afa(1,3)-3*afa(1,1)*afa(1,2)+2*afa(1,1)^3;
pui(1,4)=afa(1,4)-4*afa(1,1)*afa(1,3)+6*afa(1,1)^2*afa(1,2)-3*afa(1,1)^4;
pui(1,5)=afa(1,5)-5*afa(1,1)*afa(1,4)+10*afa(1,1)^2*afa(1,3)-10*afa(1,1)^3*afa(1,2)+4*afa(1,1)^5;
pui(1,6)=afa(1,6)-6*afa(1,1)*afa(1,5)+15*afa(1,1)^2*afa(1,4)-20*afa(1,1)^3*afa(1,3)+15*afa(1,1)^4*afa(1,2)-5*afa(1,1)^6;

return