资源描述：
《电力系统规划实验程序---发电系统裕度表生成及可靠性指标计算编程.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
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;fori=0:Naj=i+1;AF(j,1)=i;endfori
2、=0:Naj=i+1;AF(j,2)=AF(j,1)*DX;endfori=0:Naj=i+1;ifAF(j,2)<=0AF(j,5)=1-q1;AF(j,4)=0;elseifAF(j,2)>=C1AF(j,5)=q1;AF(j,4)=q1*ui;elseAF(j,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)=A
3、F(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:NbBF(i+1,2)=BF(i+1,1)*DX;endfori=0:NbifBF(i+1,2)<=0BF(i+1,5)=1-q2;BF(i+1,4)=
4、0;elseifBF(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(
5、i+1,1)*DX;endfork=0:Ndfori=0:kifk<=NaDF(k+1,5)=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:Naifk<=Naifk+1-i>0DF(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);
6、endelsesum2=0;fori=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;endendendfork=0:Ndsum1=0;forj=k:Ndsum1=DF(j+1,5)+sum1;endDF(k+1,3)=sum1;end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%形成第三台发电机组停运表%%%%%%%%%%%%%%%%%%%%%%%%%%C3=input('请输入第三台发电机的容量:C3=');q3=inpu
7、t('请输入第三台发电机的强迫停运率: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)<=0CF(i+1,5)=1-q3;CF(i+1,4)=0;elseifCF(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;endendendfori=0:Ncsum=0;fo
