82号资源-源程序：论文可在知网下载-电动汽车充电负荷时空分布预测-本人博客有解读

clc; clear; close all; load dataTRN.mat; load dataEV.mat; load dataWq.mat; G=graph(W,'upper'); %生成路-网耦合图G，W是每条边的权重值，W为32*32阶矩阵 Charge=zeros(size(Mcar,1),5); %生成一个size(Mcar,1)=1000*5的零矩阵，其中size(Mcar,1)=1000。******Charge代表什么意思*******电量吗？ TP_carnumber=zeros(size(W,1),3); %生成一个size(W,1)=32*3的零矩阵，W是每条边的权重值，W为32*32阶矩阵,TP_carnumber(x,1)、TP_carnumber(x,2)、TP_carnumber(x,3)应该是每条边上/每个交通节点私家车1，公交车2，私家车3的数量 Pcharge=zeros(24,3); %生成一个24*3的零矩阵，******PCharge代表24小时内私家车1，公交车2，私家车3待充电电量？？？？ Pchar_slow=12; %设置慢充充电功率为12KW Pchar_fast=48; %设置慢充充电功率为48KW Pntcharge=zeros(24,32); %生成一个24*32的零矩阵，PntCharge代表24个小时内32个配网节点的供电负荷=电动汽车耗电量 %% MC for icar=1:size(Mcar,1) %EV的编号，从1-1000，共1000台EV，size(Mcar,1)=1000，一个for icar=1:size(Mcar,1) 循环代表把第一台EV的 %1-9个特征变量（车辆编号icar、EV种类Icar_kind、EV初始位置Birthland、 %目的地Destination、初始时刻tBirth1、返程时刻tBirth2、ev容量Cbat、初始soc SOC0、速度Vcar）全部采集一遍 %1000个for icar=1:size(Mcar,1) 循环代表把1000台EV的 %1-9个特征变量（车辆编号icar、EV种类Icar_kind、EV初始位置Birthland、 %目的地Destination、初始时刻tBirth1、返程时刻tBirth2、ev容量Cbat、初始soc SOC0、速度Vcar）全部采集一遍 icar Icar_kind=Mcar(icar,2);%EV种类；%size(Mcar,2)=1000 Birthland=Mcar(icar,3);%EV初始位置 Destination=Mcar(icar,4);%目的地 tBirth1=Mcar(icar,5);%初始时刻 tBirth2=Mcar(icar,6);%返程时刻 Cbat=Mcar(icar,7);%ev容量 SOC0=Mcar(icar,8);%初始soc Vcar=Mcar(icar,9);%速度 TP_carnumber(Birthland,Icar_kind)=TP_carnumber(Birthland,Icar_kind)+1;%第Birthland条边/第Birthland个交通节点下，对应的第Icar_kind种车型加1，此行代码是计数的作用 if Mcar(icar,3)==Mcar(icar,4) %如果EV初始位置Birthland=car(icar,3)==EV目的地Destination=Mcar(icar,4) continue % end % [path,distance]=shortestpath(G,Birthland,Destination); %如果EV初始位置Birthland=car(icar,3) 不等于 EV目的地Destination=Mcar(icar,4) %则输出EV初始位置Birthland=car(icar,3)与EV目的地Destination=Mcar(icar,4)最短路径path与最短距离distance %% C1 私家车 if Icar_kind==1 %考虑环境温度和速度 耗电量 dE=1.5*distance/Vcar+(0.21-1e-3*Vcar+1.531/Vcar)*distance; %由最短行驶距离distance，电动汽车行驶实时车速Vcar=Mcar(icar,9)，计算出EV初始位置到目的地的耗电量 SOC1=SOC0-dE/0.9/Cbat; %由初始soc即SOC0=Mcar(icar,8) 与 EV初始位置到目的地的耗电量 与 充电效率0.9 与 ev容量Cbat=Mcar(icar,7) SOC2=SOC1-dE/0.9/Cbat; %计算出电动汽车剩余电量SOC2 tdest1=tBirth1+distance/Vcar; %到达目的地时间（tdest1）=初始时刻+行驶时间 if SOC2<0 %判断私家车是否充电的条件 % 触发慢充/计及排队时间 Tmmc=Wq(ceil(tdest1)); %ceil为向上取整，*****Wq(ceil(tdest1))代表什么意思，Tmmc代表什么意思？****** T80=tdest1+(0.8-SOC1)*Cbat/Pchar_slow+Tmmc; %充电到0.8EV容量时间（T80）=到达目的地时间+慢充到0.8EV容量时间+Tmmc,*****这个Tmmc是什么意思？****** if T80<tBirth2 %如果充电到0.8EV容量时间＜到达第二个目的地的时间 Charge(icar,2)=1; %Icar_kind=Mcar(icar,2);EV种类；所以Charge(icar,2)赋值为私家车Icar_kind==1 Charge(icar,3)=floor(tdest1); %EV初始位置Mcar(icar,3)赋值为向下取整tdest1，即32个交通节点中的一个 Charge(icar,4)=ceil(T80); %目的地Mcar(icar,4)赋值为向上取整T80，即32个交通节点中的一个 Charge(icar,5)=ceil(Pchar_slow); %初始时刻Mcar(icar,5)赋值为向上取整Pchar_slow Pntcharge(Charge(icar,3):Charge(icar,4),Destination)=Pntcharge(Charge(icar,3):Charge(icar,4),Destination)+Charge(icar,5);%列表切片，Pntcharge(Charge(icar,3):Charge(icar,4),Destination)即为取Charge(icar,3):Charge(icar,4)所在行与Destination所在列的交叉元素 else Charge(icar,2)=1; %Icar_kind=Mcar(icar,2);EV种类；所以Charge(icar,2)赋值为私家车Icar_kind==1 Charge(icar,3)=floor(tdest1); %EV初始位置Mcar(icar,3)赋值为向下取整tdest1 Charge(icar,4)=ceil(tBirth2); %目的地Mcar(icar,4)赋值为向上取整tBirth2，因为此时T80＞tBirth2成立 Charge(icar,5)=ceil(Pchar_slow); %初始时刻Mcar(icar,5)赋值为向上取整Pchar_slow Pntcharge(Charge(icar,3):Charge(icar,4),Destination)=Pntcharge(Charge(icar,3):Charge(icar,4),Destination)+Charge(icar,5);%列表切片，Pntcharge(Charge(icar,3):Charge(icar,4),Destination)即为取Charge(icar,3):Charge(icar,4)所在行与Destination所在列的交叉元素 %这一行应该是有问题的，应该是充电桩快/慢充充电功率*充电时间=电动汽车充电功率 %同理，类似的这总公式都有问题，应该将64行的Charge(icar,5)=ceil(Pchar_slow); %更改为Charge(icar,5)=ceil(dE); end end end %% 出租车 if Icar_kind==2 for iload=2:length(path) % 考虑环境温度和速度 耗电量 dE=1.5*RL(path(iload-1),path(iload))/Vcar+(0.21-1e-3*Vcar+1.531/Vcar)*RL(path(iload-1),path(iload));%由每段最短路径path，电动汽车行驶实时车速Vcar=Mcar(icar,9)，路程模型权重RL，计算出EV每段最短路径path的耗电量 SOC2=SOC0-dE/0.9/Cbat; %由初始soc即SOC0=Mcar(icar,8) 与 EV初始位置到目的地的耗电量 与 充电效率0.9 与 ev容量Cbat=Mcar(icar,7)，计算出电动汽车剩余电量SOC2 t2=tBirth1+RL(path(iload-1),path(iload))/Vcar; %到达目的地时间（t2）=初始时刻+每段最短路径path行驶时间 if SOC2<0.25 %判断出租车是否充电的条件 % 触发充电/计及排队时间 Tmmc=Wq(ceil(tBirth1)); %ceil为向上取整，*****Wq(ceil(tBirth1))代表什么意思，Tmmc代表什么意思？****** T50=tBirth1+(0.5-SOC0)*Cbat/Pchar_fast+Tmmc; %充电到0.5EV容量时间（T50）=到达目的地时间+慢充到0.5EV容量时间+Tmmc,*****这个Tmmc是什么意思？****** Charge(icar,2)=2; %Icar_kind=Mcar(icar,2);EV种类；所以Charge(icar,2)赋值为出租车Icar_kind==2 Charge(icar,3)=floor(t2); %EV初始位置Mcar(icar,3)赋值为向下取整t2，即32个交通节点中的一个 Charge(icar,4)=ceil(T50); %目的地Mcar(icar,4)赋值为向上取整T50，即32个交通节点中的一个 Charge(icar,5)=ceil(Pchar_fast); %初始时刻Mcar(icar,5)赋值为向上取整Pchar_fast Pntcharge(Charge(icar,3):Charge(icar,4),Birthland)=Pntcharge(Charge(icar,3):Charge(icar,4),Birthland)+Charge(icar,5);%列表切片，Pntcharge(Charge(icar,3):Charge(icar,4),Destination)即为取Charge(icar,3