电力系统规划实验程序---发电系统裕度表生成及可靠性指标计算编程

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)=q

2、1*ui;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

3、,1)=i;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

4、=0:Ndfor 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(

5、i+1,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;end

6、for i=0:Ncif 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

7、+1,2)=EF(i+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

8、,4)+CF(i+1,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);endf

9、or i=0:Nlif 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

10、:Nlsum= POWER (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;elsesum2=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)= (