function?[FcostCCSfcsucsNchIEpCVTCVTICDLCDLI]=VorCostCDEV（xybbcsBL）
%%?說明：計算電動汽車規(guī)劃目標函數(shù)，即適應值
%輸入
%x，y為充電站坐標；為外界交互變換參數(shù)。x和y保證是列矩陣。
%b為充電需求點坐標和充電需求，第12列為橫坐標和縱坐標，第3列為充電需求，為固定參數(shù)；
%bcs為配送區(qū)域的集中充電站坐標，第12列為橫、縱坐標。
%BL為圖坐標與實際坐標的比例，為固定參數(shù)。

%輸出
%Fcost為年社會成本
%CCS為建設(shè)運行成本
%fcs為固定建設(shè)成本
%ucs為運行管理成本，取fcs的0.1
%NchI為各充電站規(guī)模，即充電機數(shù)量
%Ep為個充電站服務(wù)的快充車輛數(shù)，即充電需求
%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