clc;
clear;
finish=zeros(1000,1);
distance=zeros(1000,1);
disch=zeros(1000,1);
ch=zeros(1000,1);
T=zeros(1000,1);
de=zeros(1000,24);
tot=zeros(1,24);
for i=1:1:1000
S=20+10*rand(1);%容量
P=rand(1)+2;%充电功率
finish(i,1)=17.6+3.4*randn(1);%最后出行结束时刻
if finish(i,1)>24
finish(i,1)=finish(i,1)-24;
end
while(finish(i,1)<0)
finish(i,1)=17.6+3.4*randn(1);
end
finish(i,1)=ceil(finish(i,1));
distance(i,1)=lognrnd(3.2,0.88);
while(distance(i,1)>160/1.61)
distance(i,1)=lognrnd(3.2,0.88);%日行驶里程
end
%%放电
temp=9*rand(1);
if (temp<=6)
disch(i,1)=10+temp;
else
disch(i,1)=20-6+temp;
end
disch(i,1)=ceil(disch(i,1));%放电开始时刻
ch(i,1)=22+7*rand(1);
ch(i,1)=ceil(ch(i,1));
ch(i,1)=ceil(ch(i,1));%充电开始时刻
while(finish(i,1)>disch(i,1)&&finish(i,1)<24)%放电在最后出行结束后开始
temp=9*rand(1);
if (temp<=6)
disch(i,1)=10+temp;
else
disch(i,1)=20-6+temp;
end
disch(i,1)=ceil(disch(i,1));
end
S2=S-distance(i,1)/100*15*1.61;
S2=round(S2);
if S2>0&&finish(i,1)<24
td=ceil(S2/P);
else
td=0;
end
if td<24-disch(i,1)&&td>0
de(i,disch(i,1):disch(i,1)+td-1)=-1*P;
elseif td>=24-disch(i,1)&&td>0
td=24-disch(i,1);
de(i,disch(i,1):disch(i,1)+td-1)=-1*P;
end
S3=S2-td*P;
tc=ceil((S-S3)/P);
if tc+ch(i,1)-1>24&&ch(i,1)<=24
de(i,ch(i,1):24)=P;
de(i,1:ch(i,1)+tc-1-24)=P;
elseif tc+ch(i,1)-1<24&&ch(i,1)<=24
de(i,ch(i,1):ch(i,1)+tc-1)=P;
elseif ch(i,1)>24
de(i,(ch(i,1)-24):(ch(i,1)+tc-1-24))=P;
end
tot(1,:)=de(i,:)+tot(1,:);
end
sigma=std(de);
figure(1)
x=1:1:24;
loadcurve=plot(x,tot(1,:)/i,'-rs',x,sigma(1,:),'-s');
set(loadcurve,'linewidth',2);
xlim([0,24]);
title('V2G功率需求');
xlabel('时间/h');
ylabel('功率需求/kW');
legend('期望','标准差');
figure(2)
x=1:24;
V2G=plot(x,1000*tot(1,:)/i,'-rs',x,1000*tot(1,:)/i+3*sigma(1,:)*1000^0.5,'--',x,1000*tot(1,:)/i-3*sigma(1,:)*1000^0.5,'--');
set(V2G,'linewidth',2);
title('V2G功率曲线');
xlabel('时间/h');
ylabel('V2G总功率需求/kW');
xlim([0,24]);
- 1
- 2
- 3
前往页