参考文献  

含风光水的虚拟电厂与配电公司协调调度模型——董文略（复现场景削减部分）

主要内容

代码主要做的是风电、光伏以及电价场景不确定性模拟，首先由一组确定性的方案，生成50种光伏场景，为了避免大规模光伏场景造成的计算困难问题，采用基于概率距离快速削减算法的场景削减法，将场景削减至5个，运行后直接给出削减后的场景以及生成的场景，并给出相应的概率，可移植以及可应用性非常强！

部分程序

%风电出力预测均值E
W=[5.8,6.7,5.8,5.1,6.3,5,6.2,6,4.1,6,7,6.8,6.5,6.9,5,5.6,6,5.8,6.2,4.7,3.3,4.4,5.6,5];
%取标准差为风电出力预测值E的5%-20%,这里x=E*10%
x=W*0.1;%风电出力的标准差x
Ws=[];
%生成一个风电场景，E+x*randn(1,24)，其中randn(1,24)为生成随机数的标准正态分布
    %产生一个随机分布的指定均值和方差的矩阵：将randn产生的结果乘以标准差，然后加上期望均值即可。
    %例如，产生均值为0.6，方差为0.1的一个5*5的随机数方式如下：
      %   x = .6 + sqrt(0.1) * randn(5)
for i=1:50
s=W+x.*randn(1,24);
Ws=[Ws;s];
end

%光伏出力预测均值E
S=[0,0,0,0,0,1,2.5,4,5,5.5,5.8,5.7,5.5,5.3,5.1,5,3.8,2.5,1.2,0,0,0,0,0];
%取标准差为光伏出力预测值E的5%-20%,这里x=E*20%
y=S*0.2;%标准差
Ss=[];
%生成一个风电场景，E+x*randn(1,24)，其中randn(1,24)为生成随机数的标准正态分布
for i=1:400
s=S+y.*randn(1,24);
Ss=[Ss;s];
end

运行结果