%%%%%%%%%%%%%%%%%%%%%%%%輸入電動汽車相關原始數據%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%以下為各種類型電動汽車數量，
%家用電動汽車數量
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）;?????%隨機抽取起始充電時刻
????????