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

实验三：发电系统裕度表生成及可靠性指标计算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;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;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: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(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;endfor 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+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,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);endfor 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: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)= (Ne-i)*DX;enddisp(第一台发电机组停运表AF);AFdisp(第二台发电机组停运表BF);BFdisp(第一台和第二台发电机组并联后停运表DF);DFdisp(第三台发电机组停运表CF);CFdisp(形成系统停运表EF);EFdisp(形成负荷停运表POWER);POWERdisp(形成发电系统裕度表FF);FF