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

《电力系统规划方案实验程序发电系统裕度表生成及可靠性指标计算编程.doc》由会员分享，可在线阅读，更多相关《电力系统规划方案实验程序发电系统裕度表生成及可靠性指标计算编程.doc（14页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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