caltable.m

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

% calculate probability of a specific value
i=0;
for x1=0:0.1:1
i=i+1;    
ei=2.7183;
e1=ei^(-x1^2.0/2);
s1=1.0/sqrt(2*pi);

% cdf
b(i,1)=normcdf(x1,0,1);
% pdf 
b(i,2)=s1*e1;
% 1th 
b(i,3)=-x1*s1*e1;
% 2th
b(i,4)=(x1^2-1)*s1*e1;
% 3th
b(i,5)=(3*x1-x1^3)*s1*e1;
% 4th
b(i,6)=(3-6*x1^2+x1^4)*s1*e1;
% 5th
b(i,7)=(-15*x1+10*x1^3-x1^5)*s1*e1;
% 6th
b(i,8)=(-15+45*x1^2-15*x1^4+x1^6)*s1*e1;
% 7th
b(i,9)=(-x1^7+21*x1^5-105*x1^3+105*x1)*s1*e1;
% 8th
b(i,10)=(x1^8-28*x1^6+210*x1^4-420*x1^2+105)*s1*e1;
% 9th
b(i,11)=(-x1^9+36*x1^7-378*x1^5+1260*x1^3-945*x1)*s1*e1;

end