%% 3.计算含风电场的发电系统可靠性指标（非序贯MC）
clc
clear
loadresult_WindFarmOutput
%文件“result_WindFarmOutput.mat”构成了风电场出力的状态模型【风力状态状态概率】相关状态计算查看百度文库“风电场出力模型matlab程序”
% 3.1 求出常规机组的出力模型，按类构成多状态模型
% RBTS发电系统中共有6类常规机组，%11台常规机组数据
% %2台5MW水电机组%% %1台10MW热电机组%% %4台20MW水电机组%% %1台20MW 热电机组%% %1台40MW水电机组%% %2台40MW热电机组%
Generator.Norm=[5 0.01
5 0.01
10 0.02
20 0.015
20 0.015
20 0.015
20 0.015
20 0.025
40 0.02
40 0.03
40 0.03];
save('process.mat');
% 3.2MC抽样机组确定机组状态
% 3.2.1计算含风电场的RBTS可靠性
% 共有7类机组，常规机组状态在StateNorm【出力概率】元胞数组中，风电状态在StateFORWeibull6【出力概率】
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);
%得到一次抽样常规机组状态
fori=1:11
if K(k,i)>Generator.Norm(i,2)
Pout(i)=Generator.Norm(i,1);
else
Pout(i)=0;
end
end
%得到一次抽样风电场的出力
t=0;
fori=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=185;
DNS(k)=max([0;RBTSload-sumPout]);
if DNS(k)>0
I=I+1; %I用来记录发生却负荷的次数sumDNS=sumDNS+DNS(k);
end
end
LOLP1=I/200000;
EENS1=sumDNS/200000*8760;
LOLE1=I/200000*8760;
disp('LOLP1=');disp(LOLP1);
disp('EENS1=');disp(EENS1);
disp('LOLE1=');disp(LOLE1);
save('resultwind.mat','LOLP1','EENS1','LOLE1'); %% 3.2.2计算常规RBTS系统可靠性
clc
clear
load process
I=0;
sumDNS=0;
DNS=zeros(100000,1);
K=rand(100000,11);%1-11常规7风电
for k=1:100000
Pout=zeros(11,1);
%得到一次抽样常规机组状态
fori=1:11
if K(k,i)>Generator.Norm(i,2)
Pout(i)=Generator.Norm(i,1);
else
Pout(i)=0;
end
end
sumPout=sum(Pout);%发电系统出力
% 统计系统状态参数
% 持续峰荷185MW
RBTSload=185;
DNS(k)=max([0;RBTSload-sumPout]);
if DNS(k)>0
I=I+1; %I用来记录发生却负荷的次数sumDNS=sumDNS+DNS(k);
end
end
% 计算LOLP EENS LOLE
LOLP2=I/100000;
EENS2=sumDNS/100000*8760;
LOLE2=I/100000*8760;
disp('LOLP2=');disp(LOLP2);
disp('EENS2=');disp(EENS2);
disp('LOLE2=');disp(LOLE2);
save('result0.mat','LOLP2','EENS2','LOLE2');。