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