电动汽车充电负荷预测

%%%%%%%%%%%%%%%%%%%%%%%%输入电动汽车相关原始数据%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%以下为各种类型电动汽车数量，
%家用电动汽车数量
Nh=958;%输入总数量 
a1=0.6;a2=0.3;a3=0.1;%设置一充二充三充比例；
Nh1=Nh*a1;Nh2=Nh*a2;Nh3=Nh*a3;%计算各数
Nh1=round（Nh1）;Nh2=round（Nh2）;Nh3=round（Nh3）;
%%%以下为各种类型电动汽车电池容量，暂设为唯一确定量；从界面输入
Eh=16;   %BYD F3DM
%%%以下为各种类型电动汽车快、中、慢充电功率，从界面输入
Pch=3.5; %家用慢充，充电功率，220V*16A=3.5kW
Pcm1=7;   %商场、停车场，L-1模式，220V*32A=7kW
Pcm2=12;  %商场、停车场，L-2模式，380V*32A=12kW
Pcm3=24;  %商场、停车场，L-3模式，380V*63A=24kW
%%%以下为家用电动汽车充电随机分布，期望值需要从界面输入
%一天一充
h1=1140;h2=0.6; %充电起始时刻；%起始荷电状态
%一天二充
h3=540;h4=0.6;  %第1次充电起始时刻；%第1次起始荷电状态
h5=1140;h6=0.6; %第2次充电起始时刻；%第2次起始荷电状态
%一天三充
h7=540;h8=0.6;   %第1次充电起始时刻；%第1次起始荷电状态
h9=840;h10=0.6;  %第2次充电起始时刻；%第2次起始荷电状态
h11=1140;h12=0.6;%第3次充电起始时刻；%第3次起始荷电状态
%%标准差，设为已知；
%一天一充
hh1=120;hh2=0.1;
%一天二充
hh3=60;hh4=0.1;hh5=120;hh6=0.1;
%一天三充
hh7=60;hh8=0.1;hh9=120;hh10=0.1;hh11=120;hh12=0.1;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%蒙特卡洛仿真参数输入%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
M=100;%monte carol 仿真次数
Bh=zeros（M1440）;%创建存放M次仿真的电动汽车充电功率原始矩阵
T=（1:1440）;%时间，以分钟为单位，横坐标
r=0; %monte carol仿真循环时，M矩阵每行下移用
%%%%%%%%%%%%%%%%%%%%%开始进行M次蒙特卡洛仿真%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:1:M
    Ph=zeros（11440）;%创建用于存放电动汽车充电功率的一维数组，初始值设为0
    r=r+1;
%%%%%%%%以下为家用电动汽车充电功率计算程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%一天一充
 for j=1:1:Nh1                
        Ts=normrnd（h1hh1）;     %随机抽取起始充电时刻
        Ts=round（Ts）;          %对随机抽取的时间取整    
        Csoc=normrnd（h2hh2）;   %随机抽取起始充电荷电状态
        Tc=（1-Csoc）*Eh*60/（0.9*Pch）;%计算充电时长，0.9为充电效率
        Tc=round（Tc）;            %对充电时长取整
        k=Ts+Tc;
       if k>1440
           k=k-1440;
           Ph（1:k）=Ph（1:k）+Pch;
           Ph（Ts:1440）=Ph（Ts:1440）+Pch;
       else 
           Ph（Ts:k）=Ph（Ts:k）+Pch;    %充电时段加上充电功率
       end
           Ph（1:1440）=Ph（1:1440）;
 end
       for t=1:1:1440        
           Bh（rt）=Ph（t）;%把每次运行Ph的monte carol仿真结果存入B矩阵的行
       end
%%一天二充
 for j=1:1:Nh2                
        %第1次充，停车场
        Ts=normrnd（h3hh3）;     %随机抽取起始充电时刻
        Ts=round（Ts）;           %对随机抽取的时间取整    
        Csoc=normrnd（h4hh4）;   %随机抽取起始充电荷电状态
        %Tc=120;                %充电时长限定在2小时
        Tc=（1-Csoc）*Eh*60*0.8/（0.9*Pcm1）;%计算充电时长，Pcs为单位充电桩
        Tc=round（Tc）;           %对充电时长取整
        k=Ts+Tc;
        Ph（Ts:k）=Ph（Ts:k）+Pcm1;    %充电时段加上充电功率
        %第2次充，家用
        Ts=normrnd（h5hh5）;     %随机抽取起始充电时刻