fitness.zip_微电网·_微电网优化_微电网系统_风光发电互补

function LOLP=XDC(x) f=x(1); g=x(2); xd=x(3); %风速建模 a=0.5;%修正系数 h=19;%实际塔架高度 h0=10;%设测风塔高度为10m V0=[3.5 3.6 4 4.7 4.3 3.7 3.1 2.9 3.1 3.4 3.7 3.5]; Vmean=(h/h0)^a*V0;%实际风速 BC=100/743; v=linspace(0,100,744); n=0; for j=1:12 n=n+1; h=1.62/(Vmean(1,j)^2)*v.*exp(-((0.9/Vmean(1,j))^2*v.^2)); p1=polyfit(v,h,20); %20是拟合阶数； h1=polyval(p1,v); %计算拟合的h值; for m=1:744 H(n,m)=h1(1,m);%拟合值矩阵 end end for i=1:12 for m=1:744 if H(i,m)<0 H(i,m)=0; end end end H;%每天每时速率分布 %每天风力发电机发电量 for i=1:12 for m=1:744 t(i,m)=24*H(i,m); end end Pw=linspace(0,0,744); for j=1:744 pw=ex5_1(v(1,j)); Pw(1,j)=pw; end for n=1:12 Ew1=0; for i=1:744 Ew1=Ew1+Pw(1,i)*t(n,i); end for j=1:31 Ew(n,j)=Ew1; end T=linspace(1,31,31); end Ew;%单台风力发电机日发电量 %一年中单台风力发电机发电量 Ef=1:365;%单台风力发电机年发电量 for i=1:31 Ef(1,i)=Ew(1,i); end for i=1:28 Ef(1,i+31)=Ew(2,i); end for i=1:31 Ef(1,i+59)=Ew(3,i); end for i=1:30 Ef(1,i+89)=Ew(4,i); end for i=1:31 Ef(1,i+120)=Ew(5,i); end for i=1:30 Ef(1,i+150)=Ew(6,i); end for i=1:31 Ef(1,i+181)=Ew(7,i); end for i=1:31 Ef(1,i+212)=Ew(8,i); end for i=1:30 Ef(1,i+242)=Ew(9,i); end for i=1:31 Ef(1,i+273)=Ew(10,i); end for i=1:30 Ef(1,i+304)=Ew(11,i); end for i=1:31 Ef(1,i+334)=Ew(12,i); end Ef;%单台风力发电机年发电量 %太阳能辐射量计算 Efficiency=0.14;%太阳能电池板效率 Ecc=3.6;%能量转换系数 Q0=[304.67 368.93 535.68 643.68 707.54 681.48 635.04 575.86 505.44 422.96 306.72 264.49];%月平均太阳能辐射 Dr=[28 25 27 25 23 28 20 21 22 25 26 30];%当月晴天数 DAY=[31 28 31 30 31 30 31 31 30 31 30 31];%每月天数； PSH=Efficiency*(Ecc^(-1))*Q0./Dr; %每月每日的PSH值 %天气模型 n=0; TQMX=linspace(0,0,31); for i=1:12; n=n+1; B=zeros(1,DAY(1,i)); idx=randperm(DAY(1,i)); B(idx(1:Dr(1,i)))=1; for m=1:DAY(1,i) TQMX(n,m)=B(1,m); %12行31列，每一行代表该月天气情况； end end TQMX;%天气模型 %光伏电池日发电量 Np=1;%假设太阳能电池数量为1； S=1.94;%单块有效发电面积; Ep1=Np*S*PSH*1;%单块晴天的日发电量，雨天的为0； for i=1:12 Ep(i,1:31)=Ep1(1,i).*TQMX(i,1:31);%每月每日光伏发电量； T=linspace(1,31,31); end Ep;%每月每日单块太阳能发电量； %一年中光伏电池板发电量 Ep;%12*31的矩阵；每月每日单块太阳能发电量； Eg=1:365;%光伏电池板年发电量矩阵 for i=1:31 Eg(1,i)=Ep(1,i); end for i=1:28 Eg(1,i+31)=Ep(2,i); end for i=1:31 Eg(1,i+59)=Ep(3,i); end for i=1:30 Eg(1,i+89)=Ep(4,i); end for i=1:31 Eg(1,i+120)=Ep(5,i); end for i=1:30 Eg(1,i+150)=Ep(6,i); end for i=1:31 Eg(1,i+181)=Ep(7,i); end for i=1:31 Eg(1,i+212)=Ep(8,i); end for i=1:30 Eg(1,i+242)=Ep(9,i); end for i=1:31 Eg(1,i+273)=Ep(10,i); end for i=1:30 Eg(1,i+304)=Ep(11,i); end for i=1:31 Eg(1,i+334)=Ep(12,i); end %失负荷概率 Ql=25*24*60*60;%当日的负荷需求量 Qb=xd*480*0.75;%蓄电池可用容量； Eb=zeros(1,356);%每日的蓄电池剩余电量 Wt=zeros(1,356);%能量不平衡量 Wti=zeros(1,356);%充电后的能量不平衡量 Slolp=zeros(1,356);%甩负荷情况的布尔量； Sloep=zeros(1,356);%电能浪费量的布尔量 %蓄电池建模； for i=1:365 Wt(1,i)=Eb(1,i)+Ef(1,i)*f+Eg(1,i)*g-Ql; Wti(1,i)=Eb(1,i)+(Ef(1,i)*f+Eg(1,i)*g-Ql)*0.85; if Wt(1,i)<=0 Eb(1,i+1)=0; else if Wt(1,i)>0&&Wt(1,i)<=Eb(1,i) Eb(1,i+1)=Wt(1,i); else if Wt(1,i)>Eb(1,i)&&Wt(1,i)<Qb Eb(1,i+1)=Wti(1,i); else if Wt(1,i)>=Qb Eb(1,i+1)=Qb; end end end end end %失负荷概率 for i=1:365 if Wt(1,i)>=0 Slolp(1,i)=0; else Slolp(1,i)=1; end end L=0; for i=1:365 L=L+((0-Slolp(1,i))*Wt(1,i)); end Qli=365*25*24*60*60; LOLP=L/Qli %清洁能源浪费概率 for i=1:365 if Wti(1,i)<=Qb Sloep(1,i)=0; else Sloep(1,i)=1; end end Le=0; for i=1:365 Le=Le+(Sloep(1,i)*(Wti(1,i)-Qb)); end Qli2=365*25*24*60*60*0.85; LOEP= Le/Qli2 %适应函数 function F = fitness(x) F = (150000*x(1)+10000*x(2)+1000*x(3))/1000000000; function Pw=ex5_1(v) %风速修正； a=0.5;z=19;z0=10; V=(z/z0)^a*v; %空气密度修正； midu0=1.293; midu=midu0*exp(-0.297*19/3048);%密度修正，且正在使用； %功率计算； Cp=0.593;Rfj=6; Vci=4.8;Vco=25;Vr=14; if 4.8<V<14 Pw=0.5*pi*Rfj^2*Cp*midu*10^-3*(V^3); elseif 14<V<25 Pw=35; else Pw=0; end