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