calcdf.m

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

function pp=CalCdf(x, mui, smgi, c3i, c4i, c5i, c6i)

% calculate probability of a specific value
x1=(x-mui)/smgi;
ei=2.7183;
e1=ei^(-x1^2.0/2);
s1=1.0/sqrt(2*pi);

% cdf
b0=normcdf(x1,0,1);
% pdf 
b1=s1*e1;
% 1th 
b2=-x1*s1*e1;
% 2th
b3=(x1^2-1)*s1*e1;
% 3th
b4=(3*x1-x1^3)*s1*e1;
% 4th
b5=(3-6*x1^2+x1^4)*s1*e1;
% 5th
b6=(-15*x1+10*x1^3-x1^5)*s1*e1;
% 6th
b7=(-15+45*x1^2-15*x1^4+x1^6)*s1*e1;


pp=b0+c3i*b3/6+c4i*b4/24+c5i*b5/120+c6i*b6/720;

return