Appilied energy论文复现:含多类型充电桩的电动汽车充电站优化配置方法程序代码!

本程序参考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)。为了处理由多类型充电桩间相互影响导致的条件场景约束,同时降低电动汽车充电站优化配置模型的复杂度,提出了针对条件场景约束的二步等价方法,并应用二阶锥松弛技术将混合整数非线性规划问题转化为混合整数二阶锥规划问题。

文中结果:

程序结果:

部分程序:

clc;clear all;
tic;
T = 96;%时段数为1小时
nb = 31;%节点数,根节点为33
nl = 30;%支路数
P1=1e-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]';
Q1=1e-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;
%电动汽车
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)];%测试数据,下为正式数据,运行速度太慢
%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 i=1:nbif sty_jd(i)==1pload(i,1:T)=P1(i).*pc_jm_w;%春 工作日pload(i,T+1:2*T)=P1(i).*pc_jm_wd;%春 周末pload(i,2*T+1:3*T)=P1(i).*px_jm_w;%夏 工作日pload(i,3*T+1:4*T)=P1(i).*px_jm_wd;%夏 周末pload(i,4*T+1:5*T)=P1(i).*pq_jm_w;%秋 工作日pload(i,5*T+1:6*T)=P1(i).*pq_jm_wd;%秋 周末pload(i,6*T+1:7*T)=P1(i).*pd_jm_w;%冬 工作日pload(i,7*T+1:8*T)=P1(i).*pd_jm_wd;%冬 周末elseif sty_jd(i)==2pload(i,1:T)=P1(i).*pc_sc_w;%春 工作日pload(i,T+1:2*T)=P1(i).*pc_sc_wd;%春 周末pload(i,2*T+1:3*T)=P1(i).*px_sc_w;%夏 工作日pload(i,3*T+1:4*T)=P1(i).*px_sc_wd;%夏 周末pload(i,4*T+1:5*T)=P1(i).*pq_sc_w;%秋 工作日pload(i,5*T+1:6*T)=P1(i).*pq_sc_wd;%秋 周末pload(i,6*T+1:7*T)=P1(i).*pd_sc_w;%冬 工作日pload(i,7*T+1:8*T)=P1(i).*pd_sc_wd;%冬 周末elsepload(i,1:T)=P1(i).*pc_bg_w;%春 工作日pload(i,T+1:2*T)=P1(i).*pc_bg_wd;%春 周末pload(i,2*T+1:3*T)=P1(i).*px_bg_w;%夏 工作日pload(i,3*T+1:4*T)=P1(i).*px_bg_wd;%夏 周末pload(i,4*T+1:5*T)=P1(i).*pq_bg_w;%秋 工作日pload(i,5*T+1:6*T)=P1(i).*pq_bg_wd;%秋 周末pload(i,6*T+1:7*T)=P1(i).*pd_bg_w;%冬 工作日pload(i,7*T+1:8*T)=P1(i).*pd_bg_wd;%冬 周末end
end
qload=repmat(Q1,1,8*T);
num_w=[];num_wd=[];
for i=1: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_w=sum(num_w);
sum_num_wd=sum(num_wd);
max_num=max(sum_num_w,sum_num_wd);
% k=1;
% for t=1:T
%     k1=1;
%     for i=1:nb
%         if stay_time_w(i,t)~=0%计算每个电动汽车的充电时长
%             for y=1:stay_time_w(i,t)
%                 st(k1,t)=t*0.25;
%             end
%             k1=k1+1;
%         end
%     end
% end
k=1;k1=1;
for t=1:Tfor i=1:nbif num_w(i,t)~=0for jj=1:num_w(i,t)soc=rand;st(k)=round(1+95*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)timeless=round(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_end=t+round(prl*(1-soc)/cdz(1,3));%超级快速充电结束时间t_end2=min(t_end+timeless,T);%不能越限超过时间尺度timeup(k1,:)=[t,zeros(1,T)];timeup(k1,t_end+2:t_end2+1)=1;%化为时序上可以节约充电桩数量 k1=k1+1;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);
%             endk=k+1;endendend
end

欢迎感兴趣的小伙伴关注下方公众号获取完整版代码,小编会不定期更新高质量的学习资料、文章和程序代码,为您的科研加油助力!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/216139.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Python实现软件中英文对照表功能

一、通过字典实现 import locale# 定义中英文对照的字典 translation_dict {account: {en: Account, zh: 账号},password: {en: Password, zh: 密码},logon: {en: Logon, zh: 登录}# 其他需要翻译的内容... }def get_translation_dict():"""根据系统语言选择对…

adb命令学习记录

1、 adb ( android debug bridge)安卓调试桥&#xff0c;用于完成电脑和手机之间的通信控制。 xcode来完成对于ios设备的操控&#xff0c;前提是有个mac电脑。 安卓系统是基于linux内核来进行开发的。 2、adb的安装: 本身 adb是 android SDK 其中自带的工具&#xff0c;用于完…

排坑指南之STM32串口接收队列定时异常导致接收失败

背景: 公司的项目,今天讲的这部分功能主要是和IC卡读取板进行串口通讯,然后将读取回来的IC卡保存在本地。我在调试的过程中发现了一个问题,上电刚开始的阶段,程序是好用的,能读取回来IC卡卡号,然后运行一段时间之后,就读取不回来卡号了,刷卡没有响应。 摘要: 讲述STM…

C语言--每日练习题--Day38

第一题 1. 下列代码的运行结果&#xff08;&#xff09; short i 65537; int j i 1; printf("i%d,j%d\n", i, j); A&#xff1a;i 65537&#xff0c;j 65538 B&#xff1a;i 1&#xff0c;j 2 C&#xff1a;i -1&#xff0c;j 0 D&#xff1a;i 1&#xff…

【基础知识】大数据概述

关键词—分布式 化整为零&#xff0c;再化零为整 大数据的定义 传统数据库处理起来困难的数据集。 发展历程 中国开源生态图谱2023 参考内容 中国开源生态图谱 2023.pdf 技术组件说明 数据集成 sqoop、dataX、flume 数据存储 hdfs、kafka 数据处理 mapreduce、hive…

30天入门Python(基础篇)——第7天:Python中运算符详解(算数运算符、赋值运算符、比较运算符、位运算符、逻辑运算符、成员运算符、身份运算符)

文章目录 专栏导读什么是算数运算符Python中的算数运算符Python中的算数运算符——+、-、*、/,和数学中一样关于【/】小知识点Python中的算数运算符—— % 取余数Python中的算数运算符—— ** 幂Python中的算数运算符—— // 取整除(向下取,也叫地板除)Python中的赋值运算…

出入职场,小白为何频繁换工作?

频繁换工作的原因 出入职场&#xff0c;有些朋友频繁更换工作&#xff0c;其原因由多方面组成。有的人是因为入职时心理预期过高&#xff0c;工作之后存在心理落差&#xff0c;有的人是因为找不到自己的职业定位&#xff0c;在工作时做得不顺手&#xff0c;还有的人是因为在职…

树莓派5安装opencv

1 建立虚拟环境 参考网站 https://www.piwheels.org/faq.html#venv 虚拟环境建立过程&#xff1a; To create a virtual environment: $ sudo apt install virtualenv python3-virtualenv -y $ virtualenv -p /usr/bin/python3 testpip sudo apt install virtualenv pytho…

C++学习笔记(十一)------has_a和use_a关系

文章目录 前言 一、has_a关系 1.1 has_a概念 1.2 has_a中构造和析构的顺序 1.3 has_a对象的内存情况 二、use_a关系&#xff08;友元关系&#xff09; 1.友元函数&#xff1a; 2.友元类 3 使用多文件编程的方式重新编辑上述代码 总结 前言 随着技术的革新&#xff0c;出现各种各…

如何设定一个N层BiLSTM的Layer,LSTM初始化

如何设定一个N层BiLSTM的Layer&#xff0c;LSTM初始化 前言干货代码解释 结语 前言 目前人工智能已经融入到我们的生活&#xff0c;大数据模型也是层出不穷。那我们就学习一些简单的模型设置。 干货 代码 # This is the original BiLSTM layer setup, def build_bilstm(inp…

【力扣100】54.螺旋矩阵

添加链接描述 class Solution:def spiralOrder(self, matrix: List[List[int]]) -> List[int]:if not matrix or not matrix[0]:return list()rows, columns len(matrix), len(matrix[0])order list()left, right, top, bottom 0, columns - 1, 0, rows - 1while left &…

通俗易懂,什么是.NET Core以及.NET Core能做什么

我们都知道.NET Core是一个可以用来构建现代、可伸缩和高性能的跨平台软件应用程序的通用开发框架。可用于为Windows、Linux和MacOS构建软件应用程序。与其他软件框架不同&#xff0c;.NET Core是最通用的框架&#xff0c;可用于构建各种软件&#xff0c;包括Web应用程序、移动…

多线程(Thread)

一、实现多线程 多线程是为了同步完成多项任务&#xff0c;提高资源使用率提高系统使用率。 1、继承Thread类 void run()在线程开启后&#xff0c;此方法将被调用执行&#xff0c;run()封装了被线程执行的代码void start()启动线程&#xff0c;Java虚拟机会调用run方法()、即…

搭建商城系统的构架如何选择?

近期有很多网友在csdn、gitee、知乎的评论区留言&#xff0c;搭建商城系统是选择单体架构还是微服务架构&#xff0c;这里先说结论&#xff0c;如果是纯电商的话&#xff0c;商城系统的架构建议选择单体架构。我们分析下微服务和单体架构的优劣势&#xff0c;就知道了。 一、什…

网易有道强力开源中英双语语音克隆

项目地址&#xff08;基于PromptTTS&#xff09;&#xff1a; https://github.com/netease-youdao/EmotiVoice EmotiVoice Docker镜像 尝试EmotiVoice最简单的方法是运行docker镜像。你需要一台带有NVidia GPU的机器。先按照Linux和Windows WSL2平台的说明安装NVidia容器工具…

ECCV 2024 European Conference on Computer Vision

目录 DatesPeople 官网&#xff1a; https://eccv2024.ecva.net/ 链接 Dates Thursday, February 29, 2024 authors Abstract registration deadlineThursday, March 7, 2024 authors Submission deadlineThursday, March 14, 2024 authors Supplementary materials deadlineM…

17.(vue3.x+vite)组件间通信方式之作用域插槽

前端技术社区总目录(订阅之前请先查看该博客) 示例效果 作用域插槽 父组件中的插槽内容是无法访问到子组件中的数据的,而作用域插槽就是解决获取子组件数据。 父组件代码 <template><div><div>父组件--Hello World!</div><Child>

c++ 快速排序

快速排序 #include <iostream> #include <vector> using namespace std;int partition(vector<int> &nums ,int left ,int right) {int pivot nums[left];int ileft; int jright;while (i<j){while (i<j && nums[j] >pivot) j--;nums[…

Numpy:数据转换成image并保存

Matplotlib库 若没有matplotlib库&#xff0c;则使用一下命令进行安装: pip install matplotlibconda install matplotlib要将 float64 类型的 NumPy 数据转换为图像并保存&#xff0c;可以使用 matplotlib 库。以下是使用 matplotlib 的示例代码&#xff1a; import numpy …

扁平疣的影响?劲松中西医结合医院结合临床经验分析

提到扁平疣&#xff0c;相信都不陌生&#xff0c;要知道它通常不会引起严重的危害&#xff0c;但仍然会给人们带来一些不适和困扰。除了影响美观之外&#xff0c;扁平疣当然还存在其他一些负面影响。为了让更多人知道扁平疣对健康的影响&#xff0c;今日特邀劲松中西医结合医院…