电动汽车有序充放电

大规模电动汽车作为一种可移动的储能设备，通过V2G技术并入电网，实现对电网的良好互动，成为一种重要的研究内容。

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