潮流上机课程设计报告华电.docVIP

《电力系统潮流上机》课程设计报告 院系 班级： 学号： 学生姓名： 指导教师： 设计周数 成绩： 日期：年月日一、课程设计的目的与要求 培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识 二、设计正文（详细内容见附录） 1.手算： 要求应用牛顿-拉夫逊法或P-Q分解法手算求解，要求精度为0.001MW。 节点1为平衡节点，电压，节点2为PQ节点，负荷功率，节点3是PV节点，，两条支路分别为，，对地支路。 2.计算机计算：编写潮流计算程序，要求如下： 2.1据给定的潮流计算任务书整理潮流计算的基础数据：节点的分类，线路模型，等值变压器模型，电压等级的归算，标幺值的计算； 2.2基础数据的计算机存储：节点数据，支路数据（包括变压器）； 2.3用牛顿-拉夫逊法计算； 2.4根据所选潮流计算方法画流程图，划分出功能模块，有数据输入模块，导纳阵形成模块，解线性方程组模块，计算不平衡功率模块，形成雅可比矩阵模块，解修正方程模块，计算线路潮流，网损，PV节点无功功率和平衡节点功率，数据输出模块； 2.5据上述模块编制程序并上机调试程序，得出潮流计算结果； 2.6源程序及其程序中的符号说明集、程序流图 简单系统如下图所示，支路数据如下： ，，， ，， 节点数据如下： ，， ，， 1)节点导纳阵 #include stdio.h #include math.h #include fstream.h #include LF.h //form node conductance matrix intMakeY( intnB, intnL, Line* sL, double** YG, double** YB ) { inti,j,l; double r,x,d1,g,b,t; for(i=0;inB;i++) for(j=0;jnB;j++) { YG[i][j]=0.0; YB[i][j]=0.0;} for(i=0;inL;i++) {r=sL[i].R; x=sL[i].X; g=r/(r*r+x*x); b=-x/(r*r+x*x); switch(sL[i].Type) {case 1://Line break; case 2://Transformer g*=1/sL[i].K; b*=1/sL[i].K; break;} YG[sL[i].NumI][sL[i].NumI]+=g; YG[sL[i].NumJ][sL[i].NumJ]+=g; YG[sL[i].NumI][sL[i].NumJ]-=g; YG[sL[i].NumJ][sL[i].NumI]-=g; YB[sL[i].NumI][sL[i].NumI]+=b+sL[i].B; YB[sL[i].NumJ][sL[i].NumJ]+=b+sL[i].B; YB[sL[i].NumI][sL[i].NumJ]-=b; YB[sL[i].NumJ][sL[i].NumI]-=b; } printf(实部：\n); for(i=0;inB;i++) {for(j=0;jnB;j++) printf(%lf\t,YG[i][j]); printf(\n);} printf(虚部：\n); for(i=0;inB;i++) {for(j=0;jnB;j++) printf(%lf\t,YB[i][j]); printf(\n);} /* Check the Y matrix */ ofstreamfout(out.txt); fout --------------Y Matrix-------------------- endl; for(i=0;inB;i++) {for(j=0;jnB;j++) fout YG[i][j] +j YB[i][j] \t; foutendl;} fout.close(); return 0;} 2）计算功率不平衡量 #include stdio.h #include math.h #include fstream.h #include LF.h //form delta p and delta q intCalDeltaPQ( intnpv, intnpq, Bus* bus, double** YG, double** YB, int* p_Jtobus, double* deltaf ) { intk,i,j; for(k=0;knpv+npq*2;k++) { i=p_Jtobus[k]; if(knpv) {