六安政务中心网站,门户网站做商城的,我要外包网站,直播平台创建公会本程序参考Applied energy论文《Optimal planning of electric vehicle charging stations comprising multi-types of charging facilities》#xff0c;文中主要对多类型充电桩的电动汽车充电站进行优化配置#xff0c;程序较为简单和基础#xff0c;具有较强的可扩展性和…本程序参考Applied energy论文《Optimal planning of electric vehicle charging stations comprising multi-types of charging facilities》文中主要对多类型充电桩的电动汽车充电站进行优化配置程序较为简单和基础具有较强的可扩展性和适用性注释清晰干货满满下面对文章和程序作简要介绍。 问题
随着越来越多的不同类型不同功率的充电设施接入配电系统即便是在同一地点、同一充电场所也逐渐呈现出多种类型充电设施共存的状态。这些充电设施在满足电动汽车车主多样化的充电需求的同时也改变了电动汽车负荷的时空分布情况使得不同类型充电设施间的相互影响日益显著并逐渐成为电动汽车充电站优化配置过程中不可忽略的因素。
创新点解决办法 本程序选择3种典型的电动汽车充电桩研究了含多类型充电桩的电动汽车充电站优化配置方法。按照充电功率的升序排列所选取的3种电动汽车充电桩依次代表了慢速充电设施Slow charging facilitySCF、快速充电设施Fast charging facilityFCF和超级快速充电设施Ultra-fast charging facilityUCF。为了处理由多类型充电桩间相互影响导致的条件场景约束同时降低电动汽车充电站优化配置模型的复杂度提出了针对条件场景约束的二步等价方法并应用二阶锥松弛技术将混合整数非线性规划问题转化为混合整数二阶锥规划问题。
文中结果 程序结果 部分程序
clc;clear all;
tic;
T 96;%时段数为1小时
nb 31;%节点数,根节点为33
nl 30;%支路数
P11e-5.*[500 200 250 180 180 120 240 240 240 240 240 160 160 400 320 300 300 300 300 800 0 80 200 100 60 40 500 200 80 30 0];
Q11e-5.*[300 90 160 90 60 60 80 80 80 120 120 90 100 75 100 125 100 100 100 320 0 40 100 40 25 15 160 120 30 10 0];
load pc_jm_w;
load pc_jm_wd;
load px_jm_w;
load px_jm_wd;
load pq_jm_w;
load pq_jm_wd;
load pd_jm_w;
load pd_jm_wd;
load pc_sc_w;
load pc_sc_wd;
load px_sc_w;
load px_sc_wd;
load pq_sc_w;
load pq_sc_wd;
load pd_sc_w;
load pd_sc_wd;
load pc_bg_w;
load pc_bg_wd;
load px_bg_w;
load px_bg_wd;
load pq_bg_w;
load pq_bg_wd;
load pd_bg_w;
load pd_bg_wd;
%电动汽车
prl100;%容量
cdz[7 30 60;400 3250 5600;40 325 560;10 10 10];%充电桩数据
cr106.5;%单位容量增容成本
cl80;%单位网损费用
plmax600*1e-4;%线路最大电流
d0.03;%折现率
num_peak[1 2 5 1.*ones(1,28)];%测试数据下为正式数据运行速度太慢
%num_peak[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];%峰值电动汽车停车数量
load arr_jm_w;
load arr_jm_wd;
load arr_sc_w;
load arr_sc_wd;
load arr_bg_w;
load arr_bg_wd;
load stay_jm_w;
load stay_jm_wd;
load stay_sc_w;
load stay_sc_wd;
load stay_bg_w;
load stay_bg_wd;
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];%节点区域类型
%构建负荷矩阵
for i1:nbif sty_jd(i)1pload(i,1:T)P1(i).*pc_jm_w;%春 工作日pload(i,T1:2*T)P1(i).*pc_jm_wd;%春 周末pload(i,2*T1:3*T)P1(i).*px_jm_w;%夏 工作日pload(i,3*T1:4*T)P1(i).*px_jm_wd;%夏 周末pload(i,4*T1:5*T)P1(i).*pq_jm_w;%秋 工作日pload(i,5*T1:6*T)P1(i).*pq_jm_wd;%秋 周末pload(i,6*T1:7*T)P1(i).*pd_jm_w;%冬 工作日pload(i,7*T1:8*T)P1(i).*pd_jm_wd;%冬 周末elseif sty_jd(i)2pload(i,1:T)P1(i).*pc_sc_w;%春 工作日pload(i,T1:2*T)P1(i).*pc_sc_wd;%春 周末pload(i,2*T1:3*T)P1(i).*px_sc_w;%夏 工作日pload(i,3*T1:4*T)P1(i).*px_sc_wd;%夏 周末pload(i,4*T1:5*T)P1(i).*pq_sc_w;%秋 工作日pload(i,5*T1:6*T)P1(i).*pq_sc_wd;%秋 周末pload(i,6*T1:7*T)P1(i).*pd_sc_w;%冬 工作日pload(i,7*T1:8*T)P1(i).*pd_sc_wd;%冬 周末elsepload(i,1:T)P1(i).*pc_bg_w;%春 工作日pload(i,T1:2*T)P1(i).*pc_bg_wd;%春 周末pload(i,2*T1:3*T)P1(i).*px_bg_w;%夏 工作日pload(i,3*T1:4*T)P1(i).*px_bg_wd;%夏 周末pload(i,4*T1:5*T)P1(i).*pq_bg_w;%秋 工作日pload(i,5*T1:6*T)P1(i).*pq_bg_wd;%秋 周末pload(i,6*T1:7*T)P1(i).*pd_bg_w;%冬 工作日pload(i,7*T1:8*T)P1(i).*pd_bg_wd;%冬 周末end
end
qloadrepmat(Q1,1,8*T);
num_w[];num_wd[];
for i1:31%分别计算工作日和周末时序停车数量if sty_jd(i)1num_w(i,:)round(num_peak(i).*arr_jm_w./max(arr_jm_w));num_wd(i,:)round(num_peak(i).*arr_jm_wd./max(arr_jm_wd));%停车时长分布,样本数量不足没法用停车时长概率曲线来计算随机产生停车时长elseif sty_jd(i)2num_w(i,:)round(num_peak(i).*arr_sc_w./max(arr_sc_w));num_wd(i,:)round(num_peak(i).*arr_sc_wd./max(arr_sc_wd));elsenum_w(i,:)round(num_peak(i).*arr_bg_w./max(arr_bg_w));num_wd(i,:)round(num_peak(i).*arr_bg_wd./max(arr_bg_wd));end
end
%建立节点电动汽车矩阵
sum_num_wsum(num_w);
sum_num_wdsum(num_wd);
max_nummax(sum_num_w,sum_num_wd);
% k1;
% for t1:T
% k11;
% for i1:nb
% if stay_time_w(i,t)~0%计算每个电动汽车的充电时长
% for y1:stay_time_w(i,t)
% st(k1,t)t*0.25;
% end
% k1k11;
% end
% end
% end
k1;k11;
for t1:Tfor i1:nbif num_w(i,t)~0for jj1:num_w(i,t)socrand;st(k)round(195*rand)*0.25;%充电时长if st(k)*cdz(1,1)prl*(1-soc)evjd_w(k,:)[t,i,soc,0,st(k),1,sty_jd(i),0];%时间节点soc充电节点充电时长充电桩选择节点区域类型节省时间elseif st(k)*cdz(1,1)prl*(1-soc) prl*(1-soc)st(k)*cdz(1,2)timelessround(prl*(1-soc)/cdz(1,2)-prl*(1-soc)/cdz(1,3));evjd_w(k,:)[t,i,soc,0,st(k),3,sty_jd(i),timeless];%时间节点soc充电节点充电时长充电桩选择节点区域类型 t_endtround(prl*(1-soc)/cdz(1,3));%超级快速充电结束时间t_end2min(t_endtimeless,T);%不能越限超过时间尺度timeup(k1,:)[t,zeros(1,T)];timeup(k1,t_end2:t_end21)1;%化为时序上可以节约充电桩数量 k1k11;elseevjd_w(k,:)[t,i,soc,0,st(k),3,sty_jd(i),0];%时间节点soc充电节点充电时长充电桩选择节点区域类型 end
% if gdch(i)~0
% evjd_w(k,4)gdch(i);
% endkk1;endendend
end
欢迎感兴趣的小伙伴关注下方公众号获取完整版代码小编会不定期更新高质量的学习资料、文章和程序代码为您的科研加油助力
