clc;
clear;
finish=zeros（10001）;
distance=zeros（10001）;
disch=zeros（10001）;
ch=zeros（10001）;
T=zeros（10001）;
de=zeros（100024）;
tot=zeros（124）;

for?i=1:1:1000
????S=20+10*rand（1）;%容量
????P=rand（1）+2;%充電功率
????finish（i1）=17.6+3.4*randn（1）;%最后出行結束時刻
????if?finish（i1）>24
????????finish（i1）=finish（i1）-24;
????end
????while（finish（i1）<0）
????????finish（i1）=17.6+3.4*randn（1）;
????end
????finish（i1）=ceil（finish（i1））;
????distance（i1）=lognrnd（3.20.88）;
????while（distance（i1）>160/1.61）
????????distance（i1）=lognrnd（3.20.88）;%日行駛里程
????end
????%%放電
????temp=9*rand（1）;
????if?（temp<=6）
????????disch（i1）=10+temp;
????else
????????disch（i1）=20-6+temp;
????end
????disch（i1）=ceil（disch（i1））;%放電開始時刻
????ch（i1）=22+7*rand（1）;
????ch（i1）=ceil（ch（i1））;
????ch（i1）=ceil（ch（i1））;%充電開始時刻
????while（finish（i1）>disch（i1）&&finish（i1）<24）%放電在最后出行結束后開始
????????temp=9*rand（1）;
????????if?（temp<=6）
????????????disch（i1）=10+temp;
????????else
?????????????disch（i1）=20-6+temp;
????????end
????????disch（i1）=ceil（disch（i1））;
????end
????S2=S-distance（i1）/100*15*1.61;
????S2=round（S2）;
????if?S2>0&&finish（i1）<24
????????td=ceil（S2/P）;
????else
????????td=0;
????end
????if?td<24-disch（i1）&&td>0
????????de（idisch（i1）:disch（i1）+td-1）=-1*P;
????elseif?td>=24-disch（i1）&&td>0
????????td=24-disch（i1）;
????????de（idisch（i1）:disch（i1）+td-1）=-1*P;
????end
????S3=S2-td*P;
????tc=ceil（（S-S3）/P）;
????if?tc+ch（i1）-1>24&&ch（i1）<=24
????????de（ich（i1）:24）=P;
????????de（i1:ch（i1）+tc-1-24）=P;
????elseif?tc+ch（i1）-1<24&&ch（i1）<=24?
????????de（ich（i1）:ch（i1）+tc-1）=P;
????elseif?ch（i1）>24?
????????de（i（ch（i1）-24）:（ch（i1）+tc-1-24））=P;
????end
????tot（1:）=de（i:）+tot（1:）;
end
sigma=std（de）;
figure（1）
x=1:1:24;
loadcurve=plot（xtot（1:）/i‘-rs‘xsigma（1:）‘-s‘）;
set（loadcurve‘linewidth‘2）;
xlim（[024]）;
title（‘V2G功率需求‘）;
xlabel（‘時間/h‘）;
ylabel（‘功率需求/kW‘）;
legend（‘期望‘‘標準差‘）;
figure（2）
x=1:24;
V2G=plot（x1000*tot（1:）/i‘-rs‘x1000*tot（1:）/i+3*sigma（1:）*1000^0.5‘--‘x1000*tot（1:）/i-3*sigma（1:）*1000^0.5‘--‘）;
set（V2G‘linewidth‘2）;
title（‘V2G功率曲線‘）;
xlabel（‘時間/h‘）;
ylabel（‘V2G總功率需求/kW‘）;
xlim（[024]）;
????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????21504??2014-04-26?12:01??電動汽車有序充放電\dailyloadcurve.xls

?????文件???????2319??2015-07-02?21:21??電動汽車有序充放電\V2GPOWER.m

?????目錄??????????0??2015-12-19?16:57??電動汽車有序充放電

-----------?---------??----------?-----??----

????????????????23823????????????????????3