WT_PV_Load_Scenario 针对电力系统中风电和光伏出力以及负荷不确定性问题

针对电力系统中风电和光伏出力以及负荷不确定性问题，该程序通过拉丁超立方抽样和样本削减将其转化为场景分析问题。

function syst
clear all
clc
%风、光、负荷场景生成
[s1p1]=lhsreduce（wtlhs（1006）5）;
[s2p2]=lhsreduce（pvlhs（1006）5）;
[s3p3]=lhsreduce（llhs（30）4）;
disp（‘风电出力各场景及概率‘）
s1 % 风电各场景
p1 % 风电各场景概率

disp（‘光伏出力各场景及概率‘）
s2 %光伏各场景
p2 %光伏各场景概率

disp（‘负荷各场景及概率‘）
s3 %负荷各场景
p3 %负荷各场景概率


% 场景的削减
function [xp]=lhsreduce（x1x2）
[n1p1]=size（x1）;
%样本概率初始化
for i=1:n1
    p（i）=1/n1;
end
%样本削减
while （n1~=x2）
    %计算每个样本的欧氏距离
    for i=1:n1
        l（ii）=inf;
        for j=1:n1
            if（i~=j）
                l（ij）=0;
                for k=1:p1
                    l（ij）=l（ij）+（x1（ik）-x1（jk））^2;
                end
                l（ij）=sqrt（l（ij））;
            end
        end
    end
    %计算样本的相应值
    for i=1:n1
        [m（i）n（i）]=min（l（i:））;
        d（i）=p（i）*m（i）;
    end
    %找出样本最小的
    [sI]=min（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2233  2015-12-12 21:45  WT_PV_Load_Scenario\WT_PV_Load_Scenario.m

     目录          0  2015-12-12 21:42  WT_PV_Load_Scenario

----------- ---------  ---------- -----  ----

                 2233                    2