《发电系统裕度表生成及可靠性指标计算.doc》由会员分享，可在线阅读，更多相关《发电系统裕度表生成及可靠性指标计算.doc（13页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、【精品文档】如有侵权，请联系网站删除，仅供学习与交流发电系统裕度表生成及可靠性指标计算.精品文档.实验一 发电机组停运表生成一、实验目的1、熟悉发电机组停运表的生成原理；2、掌握用计算机编程形成发电机组停运表的方法。二、实验原始数据及内容1、实验原始数据：某发电系统有A、B、C 三台发电机组，其容量分别为30MW、40MW 和50MW，强迫停运率分别为0.04、0.06 和0.08，平均修复时间为38.0208333 天。2、实验内容：（1）编制形成发电机组停运的程序；（2）形成实验数据给出的三台发电机组停运表。三、 实验程序形成框图四、 实验程序结果1、 输入显示2、 结果显示3、 总结果显
2、示五、 程序代码清单实验一与实验三的程序编写在一个程序中，程序代码在实验三中。六、 心得体会在编写第一个程序时c语言和matlab差距不太大，所要的数据也不多。七、 参考资料1、电力系统规划基础实验二 负荷停运表生成一、实验目的1、熟悉负荷停运表的生成原理；2、掌握用计算机编程形成负荷停运表的方法。二、实验原始数据及实验内容1、实验原始数据某系统最大负荷为100MW，负荷曲线如图1 所示。图1 某系统的日负荷曲线2、实验内容（1）编制形成负荷停运表的程序；（2）形成图1 所示的负荷停运表。三、实验程序形成框图以下两个框图：第一个是整个程序的形成框图，第二个是负荷频率表程序的形成框图。（2）负荷
3、频率表的程序形成框图如下：四、实验程序结果1、输入系统的相关信息：2、负荷频率表的形成矩阵结果如下：3、负荷停运表的形成矩阵结果如下：五、程序代码清单clear;PM=input(请输入系统的日负荷曲线对应的最大负荷PM: PM =);L=input(请输入系统的日负荷曲线对应的负荷L:L=);DX=input(请输入步长DX:DX=);T=length(L);%根据日负荷曲线确定周期T%n=PM/DX+1;FHTYB=zeros(n,5);%定义负荷停运表矩阵初值%M=zeros(n,1);for i=1:nFHTYB(i,1)=i-1;end%使负荷停运表矩阵第一列为序号%for i=1:
4、nFHTYB(i,2)= FHTYB(i,1)*DX;end%负荷停运表矩阵第二列为负荷功率%for i=1:nfor j=1:Tif L(j)=FHTYB(i,2);M(i)=M(i)+1;endendendfor i=1:nif M(i)=0FHTYB(i,5)=M(i)/T;endend%负荷停运表矩阵第五列为确切概率%for i=1:nfor j=i:nFHTYB(i,3)=FHTYB(i,3)+FHTYB(j,5);endend%负荷停运表矩阵第三列为累积概率%fhplb=zeros(n,T+2);NS=L(T)/DX+1;%扫描从零时刻开始，负荷水平取最后一个小时的负荷%for I
5、=1:TNS1=L(I)/DX+1;if NSNS1K=NS1+1;for K=NS1+1:NSfhplb(K,I)=1;fhplb(K,T+2)=fhplb(K,T+2)+1;%只在负荷减少时记录为1，则横向相加就是累积频率%endNS=NS1;I=I+1;elseK=NS;for K=NS:NS1-1fhplb(K,I)=-1;endNS=NS1;I=I+1;endendfor i=1:nfor j=1:Tfhplb(i,T+1)=fhplb(i,T+1)+fhplb(i,j);endendfor i=1:nFHTYB(i,4)=fhplb(i,T+2);end%负荷停运表矩阵第四列列为累
6、积频率% %disp(负荷频率表);fhplbdisp(负荷停运表);FHTYB%使所求的结果表现出负荷频率表矩阵，负荷停运表矩阵%六、参考资料1、电力系统规划基础2、电力系统的MATLAB/SIMULINK 仿真与应用实验三 发电系统裕度表生成及可靠性指标计算一、实验目的1、熟悉发电系统可靠性计算原理；2、掌握用计算机编程形成发电系统裕度表的方法；3、掌握根据发电系统裕度表计算发电系统可靠性指标的方法二、实验原始数据及实验内容1、实验原始数据某发电系统有3 台发电机组，其容量分别为30MW，40MW 和50MW，强迫停运率分别为0.05，0.09 和0.07，平均修复时间均为38.02083
7、33 天。该系统最大负荷为100MW，负荷曲线如图1 所示。2、实验内容（1）编制形成发电系统裕度表及计算发电系统可靠性指标的程序；（2）形成实验数据所给出的发电系统裕度表；（3）计算实验数据所给出的发电系统可靠性指标。三、实验程序形成框图四、实验程序结果五、 程序代码清单#include main()int P1,P2,P3; float q1,q2,q3; float F1,F2,F3; float MTTR1,MTTR2,MTTR3; float u1,u2,u3; float p1,p2,p3; double N110=0,M110=0,O110=0,K110=0,N210=0,M21
8、0=0,O210=0,K210=0,N310=0,M310=0,O310=0,K310=0; double N420=0,M420=0,O420=0,K420=0,N520=0,M520=0,O520=0,K520=0; int i,j; double N620=0,M620=0,O620=0,K620=0; double N730=0,M730=0,O730=0; double LOLP,EENS=0,LOLF; /机组1停运表 printf(请输入机组1容量); scanf(%d,&P1); printf(请输入机组1强迫停运率); scanf(%f,&q1); printf(请输入机组1
9、平均修复时间); scanf(%f,&MTTR1); u1=1/MTTR1; F1=q1*u1; p1=1-q1; printf(it停运容量t累积概率t累积频率t确切概率n); printf(0t0 t1 t0 t%fn,p1); printf(1t%d t%f t%f t%fn,P1,q1,F1,q1); N10=1; for(i=1;i=(P1/10);i+) N1i=q1; M10=0; for(i=1;i=(P1/10);i+) M1i=F1; O10=p1; for(i=1;i(P1/10);i+) O1i=0; O1(P1/10)=q1; K10=0-F1; for(i=1;i(
10、P1/10);i+) K1i=0; K1(P1/10)=F1;/机组2停运表 printf(请输入机组2容量); scanf(%d,&P2); printf(请输入机组2强迫停运率); scanf(%f,&q2); printf(请输入机组2平均修复时间); scanf(%f,&MTTR2); u2=1/MTTR2; F2=q2*u2; p2=1-q2; printf(it停运容量t累积概率t累积频率t确切概率n); printf(0t0 t1 t0 t%fn,p2); printf(1t%d t%f t%f t%fn,P2,q2,F2,q2); N20=1; for(i=1;i=(P2/10
11、);i+) N2i=q2; M20=0; for(i=1;i=(P2/10);i+) M2i=F2; O20=p2; for(i=1;i(P2/10);i+) O2i=0; O2(P2/10)=q2; K20=0-F2; for(i=1;i(P2/10);i+) K2i=0; K2(P2/10)=F2; /机组3停运表 printf(请输入机组3容量); scanf(%d,&P3); printf(请输入机组3强迫停运率); scanf(%f,&q3); printf(请输入机组3平均修复时间); scanf(%f,&MTTR3); u3=1/MTTR3; F3=q3*u3; p3=1-q3;
12、 printf(it停运容量t累积概率t累积频率t确切概率n); printf(0t0 t1 t0 t%fn,p3); printf(1t%d t%f t%f t%fn,P3,q3,F3,q3); N30=1; for(i=1;i=(P3/10);i+) N3i=q3; M30=0; for(i=1;i=(P3/10);i+) M3i=F3; O30=p3; for(i=1;i(P3/10);i+) O3i=0; O3(P3/10)=q3; K30=0-F3; for(i=1;i(P3/10);i+) K3i=0; K3(P3/10)=F3; /1、2机组并联N40=O10*N20+O11*N
13、20+O12*N20+O13*N20;N41=O10*N21+O11*N20+O12*N20+O13*N20;N42=O10*N23+O11*N21+O12*N20+O13*N20;M40=O10*M20+O11*M20+O12*M20+O13*M20+K10*N20+K11*N20+K12*N20+K13*N20;M41=O10*M21+O11*M20+O12*M20+O13*M20+K10*N21+K11*N20+K12*N20+K13*N20;M42=O10*M22+O11*M21+O12*M20+O13*M20+K10*N22+K11*N21+K12*N20+K13*N20;O40=O
14、10*O20+O11*0+O12*0+O13*0;O41=O10*O21+O11*O20+O12*0+O13*0;O42=O10*O21+O11*O21+O12*O20+O13*0;K40=O10*K20+K10*O20;K41=O10*K21+O11*K20+K10*O21+K10*O21;K42=O10*K22+O11*K21+O12*K20+K10*O22+K11*O21+K12*O20;for(i=3;i=(P1/10)+(P2/10);i+)N4i=O10*N2i-0+O11*N2i-1+O12*N2i-2+O13*N2i-3; M4i=O10*M2i-0+O11*M2i-1+O12
15、*M2i-2+O13*M2i-3+K10*N2i-0+K11*N2i-1+K12*N2i-2+K13*N2i-3; O4i=O10*O2i-0+O11*O2i-1+O12*O2i-2+O13*O2i-3; K4i=O10*K2i-0+O11*K2i-1+O12*K2i-2+O13*K2i-3+K10*O2i-0+K11*O2i-1+K12*O2i-2+K13*O2i-3;printf(机组1与机组2并联的停运表n);printf(it停运容量t累积概率t累积频率t确切概率n);for(i=0;i=(P1/10)+(P2/10);i+)j=10*i; printf(%dt%d t%1.10ft%
16、1.10ft%1.10fn,i,j,N4i,M4i,O4i);/全部并联N50=O30*N40+O31*N40+O32*N40+O33*N40+O34*N40+O35*N40;N51=O30*N41+O31*N40+O32*N40+O33*N40+O34*N40+O35*N40;N52=O30*N42+O31*N41+O32*N40+O33*N40+O34*N40+O35*N40;N53=O30*N43+O31*N42+O32*N41+O33*N40+O34*N40+O35*N40;N54=O30*N44+O31*N43+O32*N42+O33*N41+O34*N40+O35*N40;M50=
17、O30*M40+O31*M40+O32*M40+O33*M40+O34*M40+O35*M40+K30*N40+K31*N40+K32*N40+K33*N40+K34*N40+K35*N40;M51=O30*M41+O31*M40+O32*M40+O33*M40+O34*M40+O35*M40+K30*N41+K31*N40+K32*N40+K33*N40+K34*N40+K35*N40;M52=O30*M42+O31*M41+O32*M40+O33*M40+O34*M40+O35*M40+K30*N42+K31*N41+K32*N40+K33*N40+K34*N40+K35*N40;M53=
18、O30*M43+O31*M42+O32*M41+O33*M40+O34*M40+O35*M40+K30*N43+K31*N42+K32*N41+K33*N40+K34*N40+K35*N40;M54=O30*M44+O31*M43+O32*M42+O33*M41+O34*M40+O35*M40+K30*N44+K31*N43+K32*N42+K33*N41+K34*N40+K35*N40;O50=O30*O40;O51=O30*O41+O31*O40;O52=O30*O42+O31*O41+O32*O40;O53=O30*O43+O31*O42+O32*O41+O33*O40;O54=O30*
19、O44+O31*O43+O32*O42+O33*O41+O34*O40;for(i=5;i=(P1/10)+(P2/10)+(P3/10);i+)N5i=O30*N4i-0+O31*N4i-1+O32*N4i-2+O33*N4i-3+O34*N4i-4+O35*N4i-5; M5i=O30*M4i-0+O31*M4i-1+O32*M4i-2+O33*M4i-3+O34*M4i-4+O35*M4i-5+K30*N4i-0+K31*N4i-1+K32*N4i-2+K33*N4i-3+K34*N4i-4+K35*N4i-5; O5i=O30*O4i-0+O31*O4i-1+O32*O4i-2+O33
20、*O4i-3+O34*O4i-4+O35*O4i-5;printf(所有机组并联的停运表n);printf(it停运容量t累积概率t累积频率t确切概率n);for(i=0;i=(P1/10)+(P2/10)+(P3/10);i+)j=10*i; printf(%dt%d t%1.10ft%1.10ft%1.10fn,i,j,N5i,M5i,O5i);/负荷停运表printf(负荷停运表n);N60=1;M60=0;O60=0;N61=1;M61=0;O61=0;N62=1;M62=0;O62=0;N63=1;M63=0;O63=0;N64=1;M64=0;O64=0;N65=1;M65=0;O
21、65=0.0416667;N66=0.958333;M66=1;O66=0.208333;N67=0.75;M67=1;O67=0.125;N68=0.625;M68=1;O68=0.208333;N69=0.416667;M69=4;O69=0.25;N610=0.166667;M610=2;O610=0.166667;K65=-1;K66=0;K67=0;K68=-3;K69=2;K610=2;printf(it停运容量t累积概率t累积频率t确切概率n);for(i=0;i11;i+)j=i*10; printf(%dt%d t%f t%f t%fn,i,j,N6i,M6i,O6i);/发
22、电系统裕度表N70=O65+O66+O67+O68+O69+O610;N71=O65+O66+O67+O68+O69+O610;N72=O65+O66+O67+O68+O69+O610;N73=O65+O66+O67+O68+O69+O610;N74=O65+O66+O67+O68+O69+O610;N75=O65+O66+O67+O68+O69+O610;N76=O65*N51+O66+O67+O68+O69+O610;N77=O65*N52+O66*N51+O67+O68+O69+O610;N78=O65*N53+O66*N52+O67*N51+O68+O69+O610;N79=O65*N
23、54+O66*N53+O67*N52+O68*N51+O69+O610;M70=M71=M72=M73=M74=M75=0;M76=O65*M51+K65*N51+K66*+K67+K68+K69+K610;M77=O65*M52+O66*M51+K65*N52+K66*N51+K67+K68+K69+K610;M78=O65*M53+O66*M52+O67*M51+K65*N53+K66*N52+K67*N51+K68+K69+K610;M79=O65*M54+O66*M53+O67*M52+O68*M51+K65*N54+K66*N53+K67*N52+K68*N51+K69+K610;O
24、70=O71=O72=O73=O74=0;O75=O65*O50;O76=O65*O51+O66*O50;O77=O65*O52+O66*O51+O67*O50;O78=O65*O53+O66*O52+O67*O51+O68*O50;O79=O65*O54+O66*O53+O67*O52+O68*O51+O69*O50;for(i=10;i=(P1/10)+(P2/10)+(P3/10)+10);i+)N7i=O65*N5i-5+O66*N5i-6+O67*N5i-7+O68*N5i-8+O69*N5i-9+O610*N5i-10; M7i=O65*M5i-5+O66*M5i-6+O67*M5
25、i-7+O68*M5i-8+O69*M5i-9+O610*M5i-10+K65*N5i-5+K66*N5i-6+K67*N5i-7+K68*N5i-8+K69*N5i-9+K610*N5i-10; O7i=O65*O5i-5+O66*O5i-6+O67*O5i-7+O68*O5i-8+O69*O5i-9+O610*O5i-10;printf(发电系统的裕度表n);printf(it裕度 t累积概率t累积频率t确切概率n);for(i=0;i=(P1/10)+(P2/10)+(P3/10)+10);i+)j=(P1/10)+(P2/10)+(P3/10)*10-10*i; printf(%dt%
26、d t%1.10ft%1.10ft%1.10fn,i,j,N7i,M7i,O7i);for(i=(P1/10+P2/10+P3/10+1);i=(P1/10+P2/10+P3/10+10);i+)EENS+=N7i;EENS*=10*24;printf(电力不足概率LOLP=%1.10fn电量不足期望值EENS=%1.10fn停电频率LOLF=%1.10fn,N7P1/10+P2/10+P3/10+1,EENS,M7P1/10+P2/10+P3/10+1);六、 心得体会c语言的编写的所有的数据需要定义，且需要定义数组边界，没有matlab编写的随意，有些功能在c语言中很难实现。七、 参考资料1、电力系统规划基础