电力系统规划实验程序

1、实验三：发电系统裕度表生成及可靠性指标计算clear;%形成第一台发电机组停运表%C仁input(请输入第一台发电机的容量:C仁)；q1=input(请输入第一台发电机的强迫停运率：q仁)；t=input(请输入发电机的平均修复时间:t=);DX=input(请输入所设置的计算步长：DX=);Na=C1/DX;AF=zeros(Na+1,6);ui=1/t;fori=0:Na戸+1;AF(j,1)=i;endfori=0:Na戸+1;AF(j,2)=AF(j,1)*DX;endfori=0:Na戸+1;ifAF(j,2)=C1AF(j,5)=q1;AF(j,4)=q1*ui;elseAF(j,

2、5)=0;AF(j,4)=q1*ui;endendendfori=0:Nasum=0;forj=i:Naj=j+1;sum=sum+AF(j,5);endAF(i+1,3)=sum;endfori=0:Naifi=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);fori=0:NbBF(i+1,1)=i;endfori=0:NbB

3、F(i+1,2)=BF(i+1,1)*DX;endfori=0:NbifBF(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;endendfori=0:Nbsum=0;forj=i:Nbsum=sum+BF(j+1,5);endBF(i+1,3)=sum;end%第一台和第二台发电机并联后的停运表%Nd=Na+Nb;DF=zeros(Nd+1,5);fori=0:NdDF(i+1,1)=i;DF(i+1,2)=DF(i+1,1)*DX;endfork=0:Ndfori=0:kifk=NaDF(k+1,5

4、)=AF(i+1,5)*BF(k+1-i,5)+DF(k+1,5);elsesum=0;fori=k-Nb:Nasum=AF(i+1,5)*BF(k+1-i,5)+sum;endDF(k+1,5)=sum;endendendfork=0:Ndfori=0:Naifk0DF(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;fori=k-Nb:Nasum2=AF(i+1,5)*BF(k+1-i,4)+AF(i+1,6)*BF(k+

5、1-i,3)+sum2;endDF(k+1,4)=sum2;endendendfork=0:Ndsum1=0;forj=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);fori=0:NcCF(i+1,1)=i;endfori=0:NcCF(i+1,2)=CF(i+1,1)*DX;endfori=0:NcifCF(i+1,2)=C3CF(i+1,5)=q3;CF(

6、i+1,4)=q3*ui;elseCF(i+1,5)=0;CF(i+1,4)=q3*ui;endendendfori=0:Ncsum=0;forj=i:Ncsum=sum+CF(j+1,5);endCF(i+1,3)=sum;endfori=0:Ncifi=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);fori=0:NeEF(i+1,1)=i;EF(i+1,2)=EF(i+1,1)*DX;endfork=0:Nefori=0:kifk=NcE

7、F(k+1,5)=CF(i+1,5)*DF(k+1-i,5)+EF(k+1,5);elseifk=Ndsum=0;fori=0:Ncsum=CF(i+1,5)*DF(k+1-i,5)+sum;endEF(k+1,5)=sum;elsesum=0;fori=k-Nd:Ncsum=CF(i+1,5)*DF(k+1-i,5)+sum;endEF(k+1,5)=sum;endendendendfork=0:Nefori=0:Ncifk0EF(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

8、,6)+EF(k+1,4);endelseifkNS1K=NS1+1;forK=NS1+1:NSf(K,T+2)=f(K,T+2)+1;f(K,I)=1;endNS=NS1;I=I+1;elseK=NS;forK=NS:NS1-1f(K,I)=-1;endNS=NS1;I=I+1;endendfori=0:Nlforj=1:Tf(i+1,T+1)=f(i+1,T+1)+f(i+1,j);endendfori=0:NlPOWER(i+1,4)=f(i+1,T+2);endfori=0:Nlifi=NlPOWER(i+1,6)=POWER(i+1,4)-POWER(i+2,4);elsePOWER

9、(i+1,6)=POWER(i+1,4);endend%形成发电系统裕度表%Nf=Ne+Nl;FF=zeros(Nf+1,5);fork=5:Nffori=0:kifk=NlFF(k+1,5)=POWER(i+1,5)*EF(k+1-i,5)+FF(k+1,5);elseifk=Nesum=0;fori=0:Nlsum=POWER(i+1,5)*EF(k+1-i,5)+sum;endFF(k+1,5)=sum;elsesum=0;fori=k-Ne:Nlsum=POWER(i+1,5)*EF(k+1-i,5)+sum;endFF(k+1,5)=sum;endendendfork=0:Nffor

10、i=0:Nlifk0FF(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);endelseifk=Nesum1=0;fori=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;fori=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;endendendfori=0:Nfforj=i:NfFF(i+1,3)=FF(i+1,3)+FF(j+1,5);endendfori=0:NfFF(i+1,1)=i;endfori=0:NfFF(i+1,2)=(Ne-i