N=50000;%輸入電動汽車總數量?
Eh=24;???%電池容量
Pch=3;?%充電功率
%%%%%%%%%%%%%%%%%%%%%%%蒙特卡洛仿真參數輸入%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
M=300;%monte?carol?仿真次數
Bh=zeros（M24）;%創建存放M次仿真的電動汽車充電功率原始矩陣
T=（1:24）;%時間，以小時為單位，橫坐標
r=1;?%monte?carol仿真循環時，M矩陣每行下移用
%%%%%%%%%%%%%%%%%%%%%開始進行M次蒙特卡洛仿真%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for?i=1:1:M
????Ph=zeros（124）;%創建用于存放電動汽車充電功率的一維數組，初始值設為0?
%%%%%%%%以下為家用電動汽車充電功率計算程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????for?j=1:1:N??????????????
????????Ts=normrnd（14.473.41）;?????%隨機抽取起始充電時刻
????????Ts=round（Ts）;??????????%對隨機抽取的時間取整????
????????if?Ts<=0&&Ts>24
???????????;
????????elseif?Ts>0&&Ts<=24
????????Csoc=unifrnd（0.20.3）;???%隨機抽取起始充電荷電狀態
????????Tc=（1-Csoc）*Eh/（0.9*Pch）;