计及风电、负荷不确定性的概率潮流matlab程序,蒙特卡洛模拟半不变量级数展开(Gram-Charlie)

电力系统随机潮流概率潮流计算MATLAB程序包含 蒙特卡洛模拟法、半不变量法＋级数展开（Gram-Charlie，Cornish-Fisher）； 考虑风电不确定性（webull分布）负荷不确定性（正态分布），以IEEE34节点为例，计算节点电压、支路潮流概率密度、累计概率并绘制曲线。 有注释，附带参考文献，直观清晰。

(风俗webull分布服从尺度系数为11.8,形状参数4.15,切入切出额定风速分别为4,25,15)

可分别做某节点电压幅值（某支路有功无功）概率密度或累计分布曲线，也可以将所有节点电压（所有支路有功无功）概率密度曲线或累计分布曲线。蒙特卡洛与半步变量法做对比图。

该程序每一部分都含有子程序，相应的计算过程及解释都包含在程序中，以及参考论文对应的页码。

%威尔布分布的风速随机生成器
vci=4;  vr=15;  vco=25;             %切入、额定、切出风速vci=4;  vr=15;  vco=25; 
PnW=0.2;  lanta=0.95;             %额定功率200kw, 功率因数
m=10000;
k1=PnW/(vr-vci);
k2=-k1*vci;
FD_c = 11.8003; FD_k = 4.1416;       %c为尺度系数、k为形状系数，风速均值标准差可通过[mu,sigma]=wblstat(13.2003,4.1416)确定，约120kw(FD_c = 11.3003; FD_k = 15.1416;  )
v_FD = wblrnd(FD_c,FD_k,m,1);        %输入为尺度参数、形状参数韦伯函数生成随机变量
p_FD = zeros(1,m);
q_FD = zeros(1,m);
for i=1:m                                                             %实际发出功率很小
    if v_FD(i)<vci||v_FD(i)>vco
        p_FD(i)=0;

        ....
end

风电输出功率概率密度曲线以及累计分布曲线

节点电压概率密度曲线

无功概率密度累计分布曲线

 有功概率密度曲线累计分布曲线

风电接入点概率密度曲线