%%?私家車充電負荷計算
close?all;?%關閉打開了的所有圖形窗口
clc;?%清屏命令
clear;?%清除工作空間中所有變量
tic;?%啟動秒表計時
Pn=[7?28];?%恒功率充電，私家車常規充電功率為7kw，快速充電功率為28kw
dm=175;?%私家車最大日行駛里程
Cn=35;?%鋰電池充電容量
k=0.9;?%充電效率
N=5000;?%循環次數
M1=10000;?%私家車數量
soc1（1M1）=0;?%給矩陣預配置內存空間
while?1
???for?n=1:N
????????P1=zeros（148）;?%私家車
????????d=lognrnd（3.20.88M11）;?%抽取日行駛里程
????????soc1=1-d./dm;?%初始soc百分數
????????t2=round（（1-soc1）.*Cn*2/Pn（11）/k）;?%充電時間長度-取整假設私家車一天只充1次電
????????n1=find（t2>11）;?%返回矩陣中大于6的元素所在位置
????????t2（n1（1:end））=11;
????????for?m=1:M1
????????????????while?1
?????????????????????t1=2*normrnd（17.63.411）;?%開始充電時刻，正態分布
?????????????????????t1=round（t1）;?%取整；
?????????????????????if?t1>48||t1<=0
????????????????????????continue;
?????????????????????end
?????????????????????break;
???????????