电动汽车充电站选址定容matlab程序

在一定区域内的电动汽车充电站多目标规划选址定容的matlab程序

function [FcostCCSfcsucsNchIEpCVTCVTICDLCDLI]=VorCostCDEV（xybbcsBL）
%% 说明：计算电动汽车规划目标函数，即适应值
%输入
%x，y为充电站坐标；为外界交互变换参数。x和y保证是列矩阵。
%b为充电需求点坐标和充电需求，第12列为横坐标和纵坐标，第3列为充电需求，为固定参数；
%bcs为配送区域的集中充电站坐标，第12列为横、纵坐标。
%BL为图坐标与实际坐标的比例，为固定参数。

%输出
%Fcost为年社会成本
%CCS为建设运行成本
%fcs为固定建设成本
%ucs为运行管理成本，取fcs的0.1
%NchI为各充电站规模，即充电机数量
%Ep为个充电站服务的快充车辆数，即充电需求
%CVT为用户充电途中耗时成本
%CVTI为各充电站的用户充电途中耗时成本
%CDL为电池年配送成本
%CDLI为各更换站电池配送成本


w=100*1e4;  

q=70*1e4;   

e=1.5*1e4;

beta=25;

lmd=1.2;

v=25;   

Nchmax=12;  
Nchmin=4;

ns=4;

dmax=1200;

Dmin=500;

rr=0.08;

ms=20;

Rz=（rr*（1+rr）.^ms）./（（1+rr）.^ms-1）;

mui=0.6;

dr=1;

m=1.5;


Ep=[];  


vv=VoronoiArea（[xy]3）;


bax=b; 
for k=1:length（x）
    Ai=find（vv（:3）==k）;
    xx=vv（Ai1）;
    yy=vv（Ai2）;
    kk=convhull（xxyy）; 
    in=inpolygon（bax（:1）bax（:2）xx（kk）yy（kk））;
    bax（in4）=k;
end

for i=1:length（x）
    gb=bax（bax（:4）==i:）;
    Ep=[Ep;[sum（gb（:3））i]];
    
    
    dtI=sqrt（（（gb（:1））-x（i））.^2+（（gb（:2））-y（i））.^2）.*BL;
    
    if all（（lmd.*dtI>dmax）==0）==0
        dp=1e5;
    else
        dp=0; 
    end
   
    CVTI（i:）=（1+dp）.*365.*beta.*（sum（gb（:3）.*lmd.*dtI）./（v*1e3））;
end
CVT=sum（CVTI）;


NchI=round（mui.*Ep（:1）./ns）;

fcs=w+q.*NchI+e.*NchI.^2;

Np=1e3; 
NchIdex=find（NchINchmax）;

Dij=dist（[x‘;y‘]）.*BL;
[DaDb]=find（（lmd.*triu（Dij）>=Dmin）~=triu（ones（size（Dij））1））; 

ULI=unique（[Da;Db;NchIdex]）;  

fcs（ULI）=Np.*fcs（ULI）;   


ucs=0.1.*fcs;


CCS=sum（fcs.*Rz+ucs）;


BcsX=bcs（:1）.*ones（length（x）1）;
BcsY=bcs（:2）.*ones（length（x）1）;
rij=sqrt（（BcsX-x）.^2+（BcsY-y）.^2）.*BL./1e3;
CDLI=365.*dr.*m.*Ep（:1）.*lmd.*rij;

CDL=sum（CDLI）;


Fcost=CCS+CVT+CDL;

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

     文件    2675774  2012-03-01 16:59  matlab程序-电动汽车\matlab程序\T3.bmp

     文件       3654  2012-03-21 09:40  matlab程序-电动汽车\matlab程序\VorCostCDEV.asv

     文件       2007  2017-12-31 14:58  matlab程序-电动汽车\matlab程序\VorCostCDEV.m

     文件       8125  2012-03-19 21:57  matlab程序-电动汽车\matlab程序\VorCost_CDEV3.asv

     文件       5332  2017-12-31 15:11  matlab程序-电动汽车\matlab程序\VorCost_CDEV3.m

     文件       5249  2017-12-31 15:27  matlab程序-电动汽车\matlab程序\VorCost_CDEV4.m

     文件       1619  2017-12-31 14:54  matlab程序-电动汽车\matlab程序\VoronoiArea.m

     文件        199  2017-12-31 14:53  matlab程序-电动汽车\matlab程序\VoronoiT.m

     文件        233  2017-12-31 15:12  matlab程序-电动汽车\matlab程序\说明.txt

     目录          0  2017-12-31 14:52  matlab程序-电动汽车\matlab程序

     目录          0  2018-04-22 20:44  matlab程序-电动汽车

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

              2702192                    11