电动汽车有序充放电.rar

电动汽车的充放电模型，考虑到电动汽车的日行驶规律，建立充电模型和放电模型，得到日负荷曲线，画出参与V2G的电动汽车的总功率需求曲线，程序可运行，带数据

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