clc;
clear;

%%?抽樣
load（‘Parameters.mat‘）;
interval=zeros（196）;
EVNoCount=0;??%可提供充電的總數目
EVstarttime=zeros（1SampleNo）;
EVssoc=zeros（1SampleNo）;
EVendtime=zeros（1SampleNo）;


startt=12.25;
for?t=1:96?%96個時段?從中午12點開始計算，也就是從第50個時段開始計算?12.25為第一個時間段
????starttime=startt+0.25*（t-1）;
?%%對電動汽車的到達時間及車輛數進行進行抽樣
????if?starttime>24
????????starttime=starttime-24;
????end
????o=3.41;u=17.47;
????if?starttime<=（u-12）
???????liangshu=@（x）（1./（o.*（（2.*pi）.^0.5）））.*exp（-（（x+24-u）.^2）./（2.*（o.^2）））;%每個刷新時段內的入網量數的概率密度函數
????else??
???????liangshu=@（x）（1./（o.*（（2.*pi）.^0.5）））.*exp（-（（x-u）.^2）./（2.*（o.^2）））;%每個刷新時段內的入網量數的概率密度函數
????end
????s1=quadl（liangshustarttime-0.25starttime）;%0.25代表刷新時段為0.25小時
????intervalNo=round（s1*SampleNo）;?%該段時間的接入的車輛數目

????interval（t）=intervalNo;
???
????%%對初始SOC進行抽樣
?????reasonnumber=0;
????for?i=1:intervalNo???
??????????????
????MU=2.98;
????SIGMA=1.14;
??????R?=?lognrnd（MUSIGMA）;
?????socs=?0.9-?15*R/（100*Ci）;

????while?socs<0.1?||?socs>=0.3
?????????R?=?lognrnd（MUSIGMA）;
?????????socs=?0.9-?15*R/（100*Ci）;
????end
????socs=fix（socs*100）/100;
????endtimeone=normrnd（8.923.24）;
????while?endtimeone>=12?||?endtimeone<=2
????????endtimeone=normrnd（8.923.24）;
????end
????endtimeone=ceil（（endtimeone*60）/15）+48;
????if?（starttime>12）||（starttime????????maxJ=endtimeone-t;
???????if?（0.25*maxJ*Pi>（（0.9-socs）*Ci））
???????reasonnumber=reasonnumber+1;
???????EVstarttime（1EVNoCount+reasonnumber）=t;
???????EVssoc（1EVNoCount+reasonnumber）=socs;
???????EVendtime（1EVNoCount+reasonnumber）=endtimeone;
???????end???
????end
????end
?????EVNoCount=EVNoCount+reasonnumber;
end

RealEVstarttime=zeros（1EVNoCount）;
RealEVssoc=zeros（1EVNoCount）;
RealEVendtime=zeros（1EVNoCount）;

for?i=1:EVNoCount
????RealEVstarttime（1i）=?EVstarttime（1i）;
???RealEVssoc（1i）=EVssoc（1i）;
???RealEVendtime（1i）=EVendtime（1i）;
end

clear?EVstarttime?EVssoc?EVendtime;

%%全局優化

dt=0.25;%時間間隔

Snj=binvar（EVNoCount96）;?%控制變量，0-1型

Transformer=ones（196）;

Aj=ones（961）;

re=0.9*ones（EVNoCount1）;

Cend=ones（EVNoCount1）;

Price=[?0.687?0.687?0.687?0.687?0.687?0.869?0.869?0.869?0.869?0.687?0.687?0.687?0.365?0.365?0.365?0.365?0.365?0.365?0.365?0.365??0.869?0.869?0.869?0.869];
pricekuo=zeros（196）;
for?i=1:95
????pricekuo（i）=Price（floor（i/4）+1）;
????if（floor（（i-1）/4）???????pricekuo（i）=Price（floor（i/4））;????
????end
end
pricekuo（96）=Price（24）;
%f=Pni*dt*sum（Snj*pricekuo‘）;


f=max（fh+Pni*sum（Snj））-min（fh+Pni*sum（Snj））;
F=[];
F=set（（Pni*sum（Snj））<=?（Brl*Transformer-fh））;%構建約束條件

F=F+set（（Pi*dt*Snj*Aj+Ci*（RealEVssoc‘））>=（Ci*re））+set（（Pi*dt*Snj*Aj+Ci*（RealEVssoc‘））<=（Ci*Cend））;%多個約束條件

for?i?=?1:EVNoCount
????for?k=1:96
????????if?（k?RealEVendtime（i））
????????F?=?[FSnj（ik）<1];
????????end
????end
end

options?=?sdpsettings（‘solver‘‘cplex‘‘verbose‘1）;%‘solver‘?參數指定程序用cplex求解器（如果已經安裝，否則會報錯），如果不指定?‘solver’?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3928??2014-08-24?08:41??minimum?peak-valley\Chouyang150.m

?????文件???????3357??2014-03-18?14:22??minimum?peak-valley\Chouyang250.m

?????文件????????528??2014-03-18?14:28??minimum?peak-valley\clone0.log

?????文件???????1748??2014-03-18?14:28??minimum?peak-valley\clone1.log

?????文件????????528??2014-03-18?14:28??minimum?peak-valley\clone2.log

?????文件????????528??2014-03-18?14:28??minimum?peak-valley\clone3.log

?????文件????????987??2014-07-27?19:46??minimum?peak-valley\Parameters.m

?????文件????????988??2014-03-18?14:21??minimum?peak-valley\parameters.mat

?????文件????????953??2014-03-18?11:11??minimum?peak-valley\Parameters250.m

?????文件????????989??2014-03-18?14:22??minimum?peak-valley\parameters250.mat

?????文件???????1160??2014-07-27?19:54??minimum?peak-valley\result.m

?????文件????????823??2014-08-24?14:01??minimum?peak-valley\wx300.mat

?????文件????????823??2014-07-27?19:59??minimum?peak-valley\wx500.mat

?????目錄??????????0??2014-11-08?19:36??minimum?peak-valley

-----------?---------??----------?-----??----

????????????????17340????????????????????14