《华电潮流上机计算实验报告全解.doc》由会员分享，可在线阅读，更多相关《华电潮流上机计算实验报告全解.doc（5页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、华电潮流上机计算实验报告全解电力系统潮流上机课程设计报告 系：电气与电子工程学院院 电气1205 级：班 1121181223 号：学 王城 学生姓名：孙英云 指导教师：两周 设计周数： 成 绩：日7月7年20_日期：课程设计报告 课程 一、课程设计的目的与要求 培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识 二、设计正文（详细内容见附录） 1 手算 2 计算机计算 3 思考题 三、课程设计总结或结论（详细内容见附录） 四、参考文献 1.电力系统计算：电子数字计算机的应用，西安交通大学等合编。北京：水利电力出版社； 2.现代电力系统分析，王锡凡主编，科学出版社； 3.电力系
2、统稳态分析，陈珩，中国电力出版社，1995年，第三版； 附录（设计流程图、程序、表格、数据等） 1 课程设计报告 课程 4.机算潮流程序及结果 / dierti.cpp : 定义控制台应用程序的入口点。/ #include stdaf_.h struct Line /线路结构体 int Num,NumI,NumJ; /线路号 左节点名 右节点名 float R,_,B,K; /电阻 电抗 电纳 变比（K等于1为普通支路， 为变压器支路的变比）1不等于 ; struct Bus /节点结构体 int Num ; float Volt,Phase,GenP,GenQ,LoadP,LoadQ; in
3、t Type; ; #includestdio.h #includestring.h #includemath.h #includestdlib.h #define NBUS 4 #define NLINE 4 /_Global variables _/ int nL,nB,nVA,nSH; float _NBUS; int L; double def2_NBUS; double mn50; _/ void Gauss(double a5050,double b50, int n) /_定义高斯法2 课程设计报告 课程 int JS50; int i,j,k; float d,t,_50; F
4、ILE _fp; int L=1; for(i=0;id) /_在一行中找到一个最大值赋值d,并用JSK记住这个最大值所在的列号_/ d=fabs(akj); JSk=j; if(fabs(d)0.00001) /_用追赶法，解方程组，求未知数_/ _n-1=bn-1; for(i=n-2;i=0;i-) t=0.0; for(j=i+1;j0.000001) fprintf(fp,Y(=,-3d)=(_010_.5f,_010_.5f)n,i,j,YGij,YBij); if(fp!=NULL) fclose(fp); /_节点电压附初值 _/ for(i=1;i0.000001) fpri
5、ntf(fp,S(=,-3d)=(_010_.5f,j_010_.5f)n,i,j,sPij,sQij); 灦楲瑮?灦尬网损为_010_.5f+j_010_.3f,输电效率为_010_.3fn,dsp,dsq,100-dp); if(fp!=NULL) fclose(fp); return 0; 16 课程设计报告 课程 结果： 1导纳阵 Y( 1,1 )=( 1.01534, -8.19201) Y( 1,2 )=( -0.56148, 2.30208) Y( 1,3 )=( 0.00000, 3.66667) Y( 1,4 )=( -0.45386, 1.89107) Y( 2,1 )=(
6、 -0.56148, 2.30208) Y( 2,2 )=( 1.04225, -4.67651) Y( 2,4 )=( -0.48077, 2.40385) Y( 3,1 )=( 0.00000, 3.66667) Y( 3,3 )=( 0.00000, -3.33333) Y( 4,1 )=( -0.45386, 1.89107) Y( 4,2 )=( -0.48077, 2.40385) Y( 4,4 )=( 0.93463, -4.26159) 2设定电压初值 (0)(0)0)(0)(0)(0)(jjjjfeejfejf?0?1.0;?10;11 3312123计算功率和电压偏移； 4
7、4)(00)(0)(0)(0)0)(0)0ePefBfGePfPPBG).)?27731?(0jjjssjjjjj1144(0)(0)()0)(0)(0(0(0)0BQefQGGfeeBf)?(0.05097?( jjjjsjjsjj11同理可算出 (0)(0)(0)(0)PPPQ?0?.?52596.?01960 ,ss2222222220)0()()(0(0)UUUPPP0?0?.5?.0 ,s333s3334根据求的第一次迭代时雅可比矩阵各元素的公式计算雅可比矩阵各个元素的具体值： 066667356148.0?20.992658.32104.30208.666670.?3.806298
8、03803?1.0.56148?2.30208?04.?2302080156148.826120_21.?0 00.45269130208.?256148.006629.?1?000.4?0333466667.03?20_000020.05求高斯计算后的修正量： 17 课程设计报告 课程 )0)(0(fP?0.007919-11)(0)0eQ?-0)(0(0Pf-0.108546122)(0?J ?()0)(00.023566-eQ?22?0.127652)(0)0Pf33?0.000000?(2)0e?)(0U?336计算各节点电压的一次近似值： (1)(0)(0)eee?0.?99349)
9、0)0)(eee?0.976432221)(0)(0)(eee?1.10000?333 )(0)(0)(1fff007920.111(1)(0)(0)fff0.1)(0)(0)fff?04?。经过 四返回第三步重新迭代，并校验收敛与否，令 次迭代后，收敛条件满足，停止迭代，10?求出的电压终值： Uj0.?00764?0.984591Uj0.?10910?0.95859 2Uj0.?13002?1.092293Uj0.000001.0500047计算出平衡节点4的注入功率。jS26513036787?.?0.48各节点间功率 jS00094?.00.2443712jS11804.00?.04437?14jS024320.0.2382121 Sj12694?.00?.3117924Sj1262400.04634.41Sj17564?.0321530.42Sj0.01787?0.?04487 :9网损为 2.0589610网损效率为： 18 第 5 页 共 5 页