function?Eiva=OptimalCharging（k）
%k為每天待充電的電動汽車的數量
%該程序表示電動汽車兩天充電一次的情況

t1=datenum（‘01:00‘）;
t2=datenum（‘24:00‘）;
tx=linspace（t1t224）;
k=2140;?????????????????????????????%車輛的數量

%為變量賦初值
n=5000;
Fai=zeros（241）;
Psy=zeros（241）;
x=zeros（91）;
sd=zeros（n1）;
ED=zeros（n1）;???
tc=zeros（n1）;
N=zeros（n1）;
sumN=0;

%求取充電時間長度的概率函數
for?j=1:n
????sd（j）=2*lognrnd（3.460.51）;?????%以日行駛里程sd（x）為概率分布取n個樣本點
????ED（j）=0.15*sd（j）;???????????????%第i輛EV充電電量需求為每公里耗電量乘以日行駛里程
????tc（j）=ED（j）/（3.6*0.88）;?????????%充電持續時刻為充電電量需求除以充電功率，充電功率可近似為常數
end
for?i=1:24
????for?j=1:n
?????????if?tc（j）????????????N（i）=0;
????????else?
????????????N（i）=1;
????????end
????????sumN=sumN+N（i）;
????end????
?????Fai（i）=sumN/n;?????????????????%?充電時長不小于i的概率
?????sumN=0;
end
Fai;

%配電網負

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????144??2020-07-27?16:28??Load.txt
?????文件????????2949??2020-07-27?16:28??Process.m