minimum-peak-valley 利用cplex编程

利用cplex编程，在Matlab环境里实现电动汽车有序充电，以最小化峰谷差。

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;

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

     文件       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