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

% 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:100000Pout=zeros(12,1);%得到一次抽样常规机组状态for i=1:12if K(k1,i)Generator.Norm(i,2)Pout(i)=Generator.Norm(i,1);elsePout(i)=0;endendsumPout=sum(Pout);%发电系统出力% 统计系统状态参数% 持续峰荷 185MWRBTSload=185;DNS(k1)=max(0;RBTSload-sumPout);if DNS(k1)0I=I+1; %I 用来记录发生却负荷的次数sumDNS=sumDNS+DNS(k1);endend% 计算指标EENS=sumDNS/100000*8760;R(1)=EENS;% 以下进行可信容量的搜索（弦截法）g=1;deta=10; while abs(EENS1-R(g)=detaP(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:100000Pout=zeros(12,1);%得到一次抽样常规机组状态for i=1:12if K(k1,i)Generator.Norm(i,2)Pout(i)=Generator.Norm(i,1);elsePout(i)=0;endendsumPout=sum(Pout);%发电系统出力% 统计系统状态参数% 持续峰荷 185MWRBTSload=185;DNS(k1)=max(0;RBTSload-sumPout);if DNS(k1)0I=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:200000Pout=zeros(12,1);%得到一次抽样常规机组状态for i=1:11if K(k,i)Generator.Norm(i,2)Pout(i)=Generator.Norm(i,1);elsePout(i)=0;endend%得到一次抽样风电场的出力t=0;for i=1:Mt=t+StateFORWeibull6(i,2);if K(k,12)-t=0continueelsePout(12)=StateFORWeibull6(i,1);pwind(k)= Pout(12);breakendendsumPout=sum(Pout);%发电系统出力% 统计系统状态参数% 持续峰荷 185MWD(1)=200;RBTSload=D(1);DNS(k)=max(0;RBTSload-sumPout);if DNS(k)0I=I+1; %I 用来记录发生却负荷的次数sumDNS=sumDNS+DNS(k);endendEENS=sumDNS/200000*8760;R(1)=EENS;% 计算 ELCC 定义下的可信容量g=1;deta=10; while abs(R(g)-EENS2)detaD(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:200000Pout=zeros(12,1);%得到一次抽样常规机组状态for i=1:11if K(k,i)Generator.Norm(i,2)Pout(i)=Generator.Norm(i,1);elsePout(i)=0;endend%得到一次抽样风电场的出力t=0;for i=1:Mt=t+StateFORWeibull6(i,2);if K(k,12)-t=0continueelsePout(12)=StateFORWeibull6(i,1);pwind(k)= Pout(12);breakendendsumPout=sum(Pout);%发电系统出力% 统计系统状态参数% 持续峰荷 185MWRBTSload=D(g+1);DNS(k)=max(0;RBTSload-sumPout);if DNS(k)0I=I+1; %I 用来记录发生却负荷的次数sumDNS=sumDNS+DNS(k);endendEENS=sumDNS/200000*8760;R(g+1)=EENS;g=g+1;endDeq=D(g)-185;rate_credict=Deq/15;disp(风电等值容量)