1、实验三:发电系统裕度表生成及可靠性指标计算clear;%形成第一台发电机组停运表%C1=input(请输入第一台发电机容量:C1=);q1=input(请输入第一台发电机逼迫停运率:q1=);t=input(请输入发电机平均修复时间:t=);DX=input(请输入所设立计算步长:DX=);Na=C1/DX;AF=zeros(Na+1,6);ui=1/t;for i=0:Naj=i+1;AF(j,1)=i;endfor i=0:Naj=i+1;AF(j,2)= AF(j,1)*DX;endfor i=0:Naj=i+1;if AF(j,2)=C1AF(j,5)=q1;AF(j,4)=q1*ui
2、;elseAF(j,5)=0;AF(j,4)=q1*ui;endendendfor i=0:Nasum=0;for j=i:Naj=j+1;sum=sum+AF(j,5);endAF(i+1,3)=sum;endfor i=0:Naif i=NaAF(i+1,6)=AF(i+1,4)-AF(i+2,4);elseAF(i+1,6)=AF(i+1,4);endend%形成第二台发电机组停运表%C2=input(请输入第二台发电机容量:C2=);q2=input(请输入第二台发电机逼迫停运率:q2=);Nb=C2/DX;BF=zeros(Nb+1,5);for i=0:NbBF(i+1,1)=i;
3、endfor i=0:NbBF(i+1,2)= BF(i+1,1)*DX;endfor i=0:Nbif BF(i+1,2)=C2BF(i+1,5)=q2;BF(i+1,4)=q2*ui;elseBF(i+1,5)=0;BF(i+1,4)=q2*ui;endendfor i=0:Nbsum=0;for j=i:Nbsum=sum+BF(j+1,5);endBF(i+1,3)=sum;end%第一台和第二台发电机并联后停运表%Nd=Na+Nb;DF=zeros(Nd+1,5);for i=0:NdDF(i+1,1)=i;DF(i+1,2)=DF(i+1,1)*DX;endfor k=0:Ndfo
4、r i=0:kif k=NaDF(k+1,5)=AF(i+1,5)*BF(k+1-i,5)+DF(k+1,5);elsesum=0;for i=k-Nb:Nasum=AF(i+1,5)*BF(k+1-i,5)+sum;endDF(k+1,5)=sum;endendendfor k=0:Ndfor i=0:Naif k0DF(k+1,4)=AF(i+1,5)*BF(k+1-i,4)+AF(i+1,6)*BF(k+1-i,3)+DF(k+1,4);elseDF(k+1,4)=AF(i+1,6)+DF(k+1,4);endelsesum2=0;for i=k-Nb:Nasum2=AF(i+1,5)*
5、BF(k+1-i,4)+AF(i+1,6)*BF(k+1-i,3)+sum2;endDF(k+1,4)=sum2;endendendfor k=0:Ndsum1=0;for j=k:Ndsum1=DF(j+1,5)+sum1;endDF(k+1,3)=sum1;end%形成第三台发电机组停运表%C3=input(请输入第三台发电机容量:C3=);q3=input(请输入第三台发电机逼迫停运率:q3=);Nc=C3/DX;CF=zeros(Nc+1,6);for i=0:NcCF(i+1,1)=i;endfor i=0:NcCF(i+1,2)= CF(i+1,1)*DX;endfor i=0:N
6、cif CF(i+1,2)=C3CF(i+1,5)=q3;CF(i+1,4)=q3*ui;elseCF(i+1,5)=0;CF(i+1,4)=q3*ui;endendendfor i=0:Ncsum=0;for j=i:Ncsum=sum+CF(j+1,5);endCF(i+1,3)=sum;endfor i=0:Ncif i=NcCF(i+1,6)=CF(i+1,4)-CF(i+2,4);elseCF(i+1,6)=CF(i+1,4);endend%形成发电系统停运表%Ne=Nd+Nc;EF=zeros(Ne+1,5);for i=0:NeEF(i+1,1)=i;EF(i+1,2)=EF(i
7、+1,1)*DX;endfor k=0:Nefor i=0:kif k=NcEF(k+1,5)=CF(i+1,5)*DF(k+1-i,5)+EF(k+1,5);else if k=Ndsum=0;for i=0:Ncsum=CF(i+1,5)*DF(k+1-i,5)+sum;endEF(k+1,5)=sum;elsesum=0;for i=k-Nd:Ncsum=CF(i+1,5)*DF(k+1-i,5)+sum;endEF(k+1,5)=sum;endendendendfor k=0:Nefor i=0:Ncif k0EF(k+1,4)=CF(i+1,5)*DF(k+1-i,4)+CF(i+1
8、,6)*DF(k+1-i,3)+EF(k+1,4);elseEF(k+1,4)=CF(i+1,6)+EF(k+1,4);endelseif kNS1K=NS1+1;for K=NS1+1:NSf(K,T+2)=f(K,T+2)+1;f(K,I)=1;endNS=NS1;I=I+1;elseK=NS;for K=NS:NS1-1f(K,I)=-1;endNS=NS1;I=I+1;endendfor i=0:Nlfor j=1:Tf(i+1,T+1)=f(i+1,T+1)+f(i+1,j);endendfor i=0:NlPOWER(i+1,4)=f(i+1,T+2);endfor i=0:Nli
9、f i=NlPOWER(i+1,6)=POWER(i+1,4)-POWER(i+2,4);elsePOWER(i+1,6)=POWER(i+1,4);endend%形成发电系统裕度表%Nf=Ne+Nl;FF=zeros(Nf+1,5);for k=5:Nffor i=0:kif k=NlFF(k+1,5)= POWER (i+1,5)*EF(k+1-i,5)+FF(k+1,5);elseif k=Nesum=0;for i=0:Nlsum= POWER (i+1,5)*EF(k+1-i,5)+sum;endFF(k+1,5)=sum;elsesum=0;for i=k-Ne:Nlsum= PO
10、WER (i+1,5)*EF(k+1-i,5)+sum;endFF(k+1,5)=sum;endendendfor k=0:Nffor i=0:Nlif k0FF(k+1,4)=POWER(i+1,5)*EF(k+1-i,4)+POWER(i+1,6)*EF(k+1-i,3)+FF(k+1,4);elseFF(k+1,4)= POWER (i+1,6)+FF(k+1,4);endelseif k=Nesum1=0;for i=0:Nlsum1=POWER(i+1,5)*EF(k+1-i,4)+POWER(i+1,6)*EF(k+1-i,3)+sum1;endFF(k+1,4)=sum1;els
11、esum2=0;for i=k-Ne:Nlsum2=POWER(i+1,5)*EF(k+1-i,4)+POWER(i+1,6)*EF(k+1-i,3)+sum2;endFF(k+1,4)=sum2;endendendfor i=0:Nffor j=i:NfFF(i+1,3)=FF(i+1,3)+FF(j+1,5);endendfor i=0:NfFF(i+1,1)=i;endfor i=0:NfFF(i+1,2)= (Ne-i)*DX;enddisp(第一台发电机组停运表AF);AFdisp(第二台发电机组停运表BF);BFdisp(第一台和第二台发电机组并联后停运表DF);DFdisp(第三台发电机组停运表CF);CFdisp(形成系统停运表EF);EFdisp(形成负荷停运表POWER);POWERdisp(形成发电系统裕度表FF);FF