基于EFC、ELCC定义下的风电可信容量计算matlab程序.docxVIP

% 4.1计算风电等值容量（EFC定义下）%相关可靠性计算方法参见百度文库“含风电的发电系统可靠性评估（MC法）matlab程序”clearload processload resultwindload result0% 计算EFC定义下风电等值容量P0=0;%设定初值Pe下限为0R0=EENS2;P=;%设定变量P、R存储每次迭代后Pe和对应的可靠性指标R=;% 计算初值P(1)的可靠性指标,P(1)对应等值容量上限，和风电容量相同P(1)=Generator.Wind.Pr*N;Generator.Norm=Generator.Norm;P(1),0;%第12行【15 0】是EFC定义下的理想机组I=0;sumDNS=0;DNS=zeros(100000,1);K=rand(100000,12);%1-11常规 12理想机组for k1=1:100000 Pout=zeros(12,1);%得到一次抽样常规机组状态 for i=1:12 if K(k1,i)Generator.Norm(i,2) Pout(i)=Generator.Norm(i,1); else Pout(i)=0; end end sumPout=sum(Pout);%发电系统出力 % 统计系统状态参数 % 持续峰荷185MW RBTSload=185; DNS(k1)=max(0;RBTSload-sumPout); if DNS(k1)0 I=I+1; %I用来记录发生却负荷的次数 sumDNS=sumDNS+DNS(k1); endend% 计算指标EENS=sumDNS/100000*8760;R(1)=EENS;% 以下进行可信容量的搜索（弦截法）g=1;deta=10; while abs(EENS1-R(g)=deta P(g+1)=(EENS1-R0)*P(g)/(R(g)-R0); Generator.Norm=Generator.Norm(1:11,:);P(g+1),0; % 重复进行含理想机组的可靠性评估I=0;sumDNS=0;DNS=zeros(100000,1);K=rand(100000,12);%1-11常规 12理想机组for k1=1:100000 Pout=zeros(12,1);%得到一次抽样常规机组状态 for i=1:12 if K(k1,i)Generator.Norm(i,2) Pout(i)=Generator.Norm(i,1); else Pout(i)=0; end end sumPout=sum(Pout);%发电系统出力 % 统计系统状态参数 % 持续峰荷185MW RBTSload=185; DNS(k1)=max(0;RBTSload-sumPout); if DNS(k1)0 I=I+1; %I用来记录发生却负荷的次数 sumDNS=sumDNS+DNS(k1); endendEENS=sumDNS/100000*8760;R(g+1)=EENS;g=g+1;endPeq=P(g);rate_credict=Peq/15;disp(风电等值容量);disp(Peq);disp(风电容量可信度);disp(rate_credict);resultEFC.Peq=Peq;resultEFC.rate_credict=rate_credict;save(result_WindCredict.mat,resultEFC)% 4.2计算风电可信容量（ELCC定义下）clearload result0load processload resultwindD=;R=;% 计算200MW负荷对应的可靠性指标R1I=0;sumDNS=0;DNS=zeros(200000,1);K=rand(200000,12);%1-11常规 12风电pwind=zeros(200000,1);for k=1:200000 Pout=zeros(12,1);%得到一次抽样常规机组状态 for i=1:11 if K(k,i)Generator.Norm(i,2) Pout(i)=Generator.Norm(i,1); else Pout(i)=0; end end %得到一次抽样风电场的出力 t=0; for i=1:M t=t+StateFORWeibull6(i,2); if K(k,12)-t=0 continue else Pout(12)=StateFORWeibull6(i,1); pwind(k)= Pout(12); break end end sumPout=sum(Pout);%发电系统出力 % 统计系统状态参数 % 持续峰荷185MW D(1)=200; RBTSload=D(1); DNS(k)=max(0;RBTSload-sumPout); if DNS(k)0 I=I+1; %I用来记录发生却负荷的次数 sumDNS=sumDNS+DNS(k); endendEENS=sumDNS/200000*8760;R(1)=EENS;% 计算ELCC定义下的可信容量g=1;deta=10; while abs(R(g)-EENS2)deta D(g+1)=(EENS2-EENS1)/(R(g)-EENS1)*(D(g)-185)+185; I=0; %I用来记录发生却负荷的次数 sumDNS=0; DNS=zeros(200000,1); K=rand(200000,12);%1-11常规 12风电 pwind=zeros(200000,1); for k=1:200000 Pout=zeros(12,1);%得到一次抽样常规机组状态 for i=1:11 if K(k,i)Generator.Norm(i,2) Pout(i)=Generator.Norm(i,1); else Pout(i)=0; end end %得到一次抽样风电场的出力 t=0; for i=1:M t=t+StateFORWeibull6(i,2); if K(k,12)-t=0 continue else Pout(12)=StateFORWeibull6(i,1); pwind(k)= Pout(12); break end end sumPout=sum(Pout);%发电系统出力 % 统计系统状态参数 % 持续峰荷185MW RBTSload=D(g+1); DNS(k)=max(0;RBTSload-sumPout); if DNS(k)0 I=I+1; %I用来记录发生却负荷的次数 sumDNS=sumDNS+DNS(k); end end EENS=sumDNS/200000*8760; R(g+1)=EENS; g=g+1; endDeq=D(g)-185;rate_credict=Deq/15;disp(风电等