matlab代码：基于概率距离快速削减法的风光场景生成与削减方法

采用蒙特卡洛进行场景生成，并再次进行场景缩减。

clear;clc;
%风电出力预测均值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;
Ws=[];
%生成一个风电场景，E+x*randn(1,24)，其中randn(1,24)为生成随机数的标准正态分布
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*10%
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