简介
本程序参考Applied energy论文《Optimal planning of electric vehicle charging stations comprising multi-types of charging facilities》，文中主要对多类型充电桩的电动汽车充电站进行优化配置，程序较为简单和基础，具有较强的可扩展性和适用性，注释清晰，干货满满，下面对文章和程序作简要介绍。
问题及解决办法
问题背景
随着越来越多的不同类型不同功率的充电设施接入配电系统，即便是在同一地点、同一充电场所，也逐渐呈现出多种类型充电设施共存的状态。这些充电设施在满足电动汽车车主多样化的充电需求的同时，也改变了电动汽车负荷的时空分布情况，使得不同类型充电设施间的相互影响日益显著，并逐渐成为电动汽车充电站优化配置过程中不可忽略的因素。
解决办法
本程序选择3种典型的电动汽车充电桩，研究了含多类型充电桩的电动汽车充电站优化配置方法。按照充电功率的升序排列，所选取的3种电动汽车充电桩依次代表了慢速充电设施（Slow charging facility，SCF）、快速充电设施（Fast charging facility，FCF）和超级快速充电设施（Ultra-fast charging facility，UCF）。为了处理由多类型充电桩间相互影响导致的条件场景约束，同时降低电动汽车充电站优化配置模型的复杂度，提出了针对条件场景约束的二步等价方法，并应用二阶锥松弛技术将混合整数非线性规划问题转化为混合整数二阶锥规划问题。
文中结果
程序结果
部分程序
%电动汽车prl=100;%容量cdz=[7 30 60;400 3250 5600;40 325 560;10 10 10];%充电桩数据cr=106.5;%单位容量增容成本cl=80;%单位网损费用plmax=600*1e-4;%线路最大电流d=0.03;%折现率num_peak=[1 2 5 1.*ones(1,28)];%测试数据，下为正式数据，运行速度太慢=[30 12 15 11 11 7 14 14 14 14 14 10 10 24 19 18 18 18 18 48 0 5 12 6 4 2 30 12 5 2 0]';%峰值电动汽车停车数量sty_jd=[2 2 1 3 3 2 1 1 1 1 1 1 1 2 2 2 2 1 1 2 1 3 2 2 3 3 3 1 2 3 1];%节点区域类型%构建负荷矩阵:T)=P1(i).*pc_jm_w;%春 工作:2*T)=P1(i).*pc_jm_wd;%春 周末:3*T)=P1(i).*px_jm_w;%夏 工作日:4*T)=P1(i).*px_jm_wd;%夏 周末:5*T)=P1(i).*pq_jm_w;%秋 工作日:6*T)=P1(i).*pq_jm_wd;%秋 周末:7*T)=P1(i).*pd_jm_w;%冬 工作日:8*T)=P1(i).*pd_jm_wd;%冬 周末:T)=P1(i).*pc_sc_w;%春 工作日:2*T)=P1(i).*pc_sc_wd;%春 周末:3*T)=P1(i).*px_sc_w;%夏 工作日:4*T)=P1(i).*px_sc_wd;%夏 周末qload=repmat(Q1,1,8*T);num_w=[];num_wd=[];for i=1:31%分别计算工作日和周末时序停车数量==1:)=round(num_peak(i).*arr_jm_w./max(arr_jm_w));:)=round(num_peak(i).*arr_jm_wd./max(arr_jm_wd));%停车时长分布,样本数量不足，没法用停车时长概率曲线来计算，随机产生停车时长==2:)=round(num_peak(i).*arr_sc_w./max(arr_sc_w));:)=round(num_peak(i).*arr_sc_wd./max(arr_sc_wd));:)=round(num_peak(i).*arr_bg_w./max(arr_bg_w));:)=round(num_peak(i).*arr_bg_wd./max(arr_bg_wd));=0%计算每个电动汽车的充电时长=1:stay_time_w(i,t)=t*0.25;=round(1+95*rand)*0.25;%充电时长=prl*(1-soc:)=[t,i,soc,0,st(k),1,sty_jd(i),0];%时间，节点，soc，充电节点，充电时长，充电桩选择，节点区域类型，节省时间=prl*(1-soc) && prl*(1-soc)<=st(k)*cdz(1,2)timeless=round(prl*(1-soc)/cdz(1,2)-prl*(1-soc)/cdz(1,3));:)=[t,i,soc,0,st(k),3,sty_jd(i),timeless];%时间，节点，soc，充电节点，充电时长，充电桩选择，节点区域类型t_end=t+round(prl*(1-soc)/cdz(1,3));%超级快速充电结束时间t_end2=min(t_end+timeless,T);%不能越限超过时间尺度:)=[t,zeros(1,T)];:t_end2+1)=1;%化为时序上可以节约充电桩数量
部分内容源自网络，侵权联系删除！
欢迎感兴趣的小伙伴点击文末阅读原文获取完整版代码，小编会不定期更新高质量的学习资料、文章和程序代码，为您的科研加油助力！
更多创新性预测及优化程序请点击公众号首页：《预测优化》|《综合能源》链接！
