电力系统暂态上机计算课程设计报告
-0.5489+ 0.3169i
5
0.6091
0
0
0.5977+ 0.3451i
0.0000 - 0.6901i
-0.5977+ 0.3451i
c.各支路三序电流、三相电流：
支路
I1
I2
I0
Ia
Ib
Ic
1-2
0 + 0.8451i
0
0
0 + 0.8451i
0.7319- 0.4225i
3.《电力系统故障的计算机辅助分析》 大学 米麟书等
4.《电力系统潮流计算》 大学 宋文南等
5.《电力系统故障分析》 清华大学 周荣光
6.《短路电流实用计算方法》 电力工业 交通大学等
7.《精通Matlab6.5》航空航天大学，志涌
附录（设计流程图、计算结果、思考题答案）
1、对称短路计算流程图及计算结果
-0.7319- 0.4225i
1-3
0 - 4.5070i
0
0
0 - 4.5070i
-3.9032+ 2.2535i
3.9032+ 2.2535i
2-3
0 - 5.3521i
0
0
0 - 5.3521i
-4.6351+ 2.6761i
4.6351+ 2.6761i
1-4
0 + 3.6620i
0
0
0+3.6620i
另外，通过此次课程设计，我学习并掌握了一定的Matlab知识，对简单的选择语句、循环语句、数组的使用等有了一些了解，为今后学习、工作中再次学习使用此软件打下了基础，获益匪浅。
四、参考文献
1.《电力系统暂态分析》，光琦，中国电力，2006年，第三版；
2.《电力系统分析》（上、下），何仰赞，华中科技大学，1998年，第二版。
0.0000 - 1.0000i
-0.6968 + 0.5000i
c.各支路三序电流、三相电流：
支路
I1
3.1714-
1.8310i
-3.1714-
1.8310i
2-5
0 + 6.1972i
0
0
0+6.1972i
5.3669-
3.0968i
-5.3669-
3.0968i
发电机1
-1.8310 + 3.1714i
0
0
-1.8310 + 3.1714i
3.6620
-1.8310 -3.1714i
发电机23
b.各节点三序电压、三相电压：
节点
U1
U2
U0
Ua
Ub
Uc
1
0.8268
-0.1732
-0.0715
0.5822
-0.3983 - 0.8660i
-0.3983 + 0.8660i
2
0.8535
-0.1465
-0.0459
0.6610
-0.3994 - 0.8660i
-0.3994 + 0.8660i
二、设计正文（详细容见附录,用A4纸，页数不限）
1．对称短路计算过程流程图和计算结果
2．不对称短路计算过程流程图和计算结果
3．静态稳定性计算过程流程图和计算结果
4．暂态稳定性计算过程流程图和计算结果
5、思考题
6、实验程序
三、课程设计总结或结论
本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。在短路电流的计算过程中，先通过计算系统各个节点的导纳矩阵，再通过判断系统的故障节点号和故障类型代入相应的公式进行计算。其中，将故障后的电气量分成故障分量和正常分量，通过故障类型判断故障分量和正常分量的值，相加后通过对称分量法的矩阵转换即可得出所求的电气量的数值。在静态稳定性的计算中，计算出发电机中的各个参数值和放大倍数，根据相应的对应关系通过特征值判据判断稳定性，进而计算Pe的值，从而判断其最大功率。在暂态稳定性的计算中，设置了相应的切除时间和分段时间，每一个分段时间进行一次计算，在切除故障前和后分别进行计算，进而得出相应曲线，并通过多次尝试的除最大切除时间大约为2s。通过此次学习，了解到电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。在分析过程常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。
b.各节点三序电压、三相电压：
节点
U1
U2
U0
Ua
Ub
Uc
1
0.4507
0
0
0.4507
-0.2254- 0.3903i
-0.2254+ 0.3903i
2
0.5352
0
0
0.5352
-0.2676- 0.4635i
-0.2676+ 0.4635i
3
0
0
0
0
0
0
4
0.6338
0
0
0.5489+ 0.3169i
-3.0986 + 5.3669i
0
0
-3.0986 + 5.3669i
6.1972 - 0.0000i
-3.0986 -5.3669i
2、不对称短路计算的流程图及计算结果
流程图
计算结果
(1) 节点3发生A相短路接地故障
a.故障点三序电流、三相电流：
I1=0 - 3.1080iI2=0 - 3.1080iI0=0 - 3.1080i; Ia=0 - 9.3241iIb=0 - 0.0000iI3年度第二学期)
名 称：电力系统暂态上机
院 系：电气与电子工程学院
班 级：
学 号：
学生：
指导教师：
设计周数：两周
成 绩：
日期： 2013年 7月1日
一、课程设计的目的与要求
巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。
过程流程图
计算结果节点导纳矩阵
图1 Y1、Y2的计算结果
图2 Y0的计算结果
图3 YY1、YY2的结算结果
(2) 节点3发生三相短路故障
a、故障点三序电流、三相电流：
I1=0 - 9.8592iI2=0 I0=0 ; Ia=0 - 9.8592iIb=-8.5383+ 4.9296iIc=8.5383+ 4.9296i
3
0.6848
-0.3152
-0.3695
0.0000
-0.5543 - 0.8660i
-0.5543 + 0.8660i
4
0.8846
-0.1154
0
0.6661 + 0.5000i
0.0000 - 1.0000i
-0.6661 + 0.5000i
5
0.9023
-0.0977
0
0.6968 + 0.5000i
电力系统暂态的课程设计
电力系统暂态的课程设计一、课程目标知识目标：1. 学生能够理解电力系统暂态过程的基本概念，掌握暂态过程中的关键参数及其相互影响。
2. 学生能够描述常见电力系统暂态现象，如短路、故障、暂态稳定性问题等，并分析其原因及危害。
3. 学生能够掌握电力系统暂态稳定性分析的基本方法，并运用相关公式进行简单计算。
技能目标：1. 学生能够运用所学知识，分析实际电力系统暂态问题，提出解决策略。
2. 学生能够通过查阅资料，了解电力系统暂态过程的最新研究动态和发展趋势。
3. 学生能够利用仿真软件对电力系统暂态过程进行模拟，验证所学理论和方法。
情感态度价值观目标：1. 培养学生对电力系统暂态问题的关注和兴趣，激发其探索精神。
2. 增强学生的团队合作意识，使其在分析和解决电力系统暂态问题时，能够积极与他人交流与合作。
3. 培养学生严谨、务实的科学态度，使其在分析问题时，能够充分考虑到实际工况，注重实际应用。
课程性质：本课程为电力系统专业课程，旨在帮助学生深入了解电力系统暂态过程，掌握相关分析方法和解决策略。
学生特点：学生具备一定的电力系统基础知识，具有较强的学习能力和动手能力。
教学要求：结合学生特点，注重理论与实践相结合，提高学生的实际操作能力和解决问题的能力。
在教学过程中，注重启发式教学，引导学生主动思考，培养其创新意识。
二、教学内容1. 电力系统暂态过程基本概念：包括暂态过程的定义、分类及其在电力系统中的重要性。
- 教材章节：第二章 电力系统暂态过程概述- 内容列举：暂态过程的产生原因、暂态过程的基本类型及其特点。
2. 常见电力系统暂态现象及分析：如短路、故障、暂态稳定性问题等。
- 教材章节：第三章 常见电力系统暂态现象- 内容列举：短路故障的分类及影响、暂态稳定性问题的分析方法及其防治措施。
3. 电力系统暂态稳定性分析：介绍暂态稳定性分析的基本原理、方法和应用。
- 教材章节：第四章 电力系统暂态稳定性分析- 内容列举：暂态稳定性分析的基本理论、时域仿真法、直接法等分析方法。
电力系统暂态课程设计
电力系统暂态课程设计一、课程目标知识目标：1. 掌握电力系统暂态过程的基本概念、分类及特点；2. 理解电力系统暂态稳定性分析的基本原理和方法；3. 了解电力系统暂态过程中的主要影响因素及其作用机理。
技能目标：1. 能够运用所学知识对电力系统暂态过程进行正确分析；2. 能够运用暂态稳定性分析的方法，评估电力系统的稳定性；3. 能够根据实际案例，提出改善电力系统暂态稳定性的措施。
情感态度价值观目标：1. 培养学生对电力系统暂态过程研究的兴趣，激发学生探索精神；2. 增强学生的团队协作意识，培养合作解决问题的能力；3. 提高学生关注电力系统安全、稳定的意识，树立社会责任感。
本课程针对高年级学生，结合电力系统暂态过程的学科特点，注重理论与实践相结合。
通过本课程的学习，使学生能够深入理解电力系统暂态过程的本质，掌握相关分析方法和技能，为今后从事电力系统运行、维护和管理等工作奠定基础。
同时，培养学生关注社会热点问题，提高解决实际问题的能力。
课程目标的设定旨在使学生在知识、技能和情感态度价值观方面取得全面发展。
二、教学内容1. 电力系统暂态过程基本概念：暂态过程的定义、分类及其在电力系统中的重要性；教材章节：第一章第一节。
2. 电力系统暂态稳定性分析原理：暂态稳定性分析的基本理论、方法及其应用；教材章节：第二章。
3. 电力系统暂态过程中的影响因素：分析电力系统暂态过程中的主要因素，如发电机、变压器、线路等；教材章节：第三章。
4. 暂态稳定性分析的方法与案例：介绍暂态稳定性分析的方法，结合实际案例进行分析；教材章节：第四章。
5. 改善电力系统暂态稳定性的措施：探讨提高暂态稳定性的措施，如励磁系统、PSS、直流输电等；教材章节：第五章。
6. 电力系统暂态过程实验与仿真：组织学生进行暂态过程实验，利用仿真软件进行模拟分析；教材章节：第六章。
教学内容按照课程目标进行科学、系统地组织和安排，注重理论与实践相结合。
教学进度根据教材章节进行制定，确保学生能够在掌握基本理论和方法的基础上，通过实际案例和实验，深入理解电力系统暂态过程及其稳定性问题。
华电潮流上机计算实验报告全解
院系：电气与电子工程学院班级：电气1205学号：1121181223学生姓名：王城指导教师：孙英云设计周数：两周成绩：日期：2015年7月7日一、课程设计的目的与要求培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识二、设计正文（详细内容见附录）1．手算2．计算机计算3．思考题三、课程设计总结或结论（详细内容见附录）四、参考文献1.《电力系统计算：电子数字计算机的应用》，西安交通大学等合编。
北京：水利电力出版社；2.《现代电力系统分析》，王锡凡主编，科学出版社；3.《电力系统稳态分析》，陈珩，中国电力出版社，1995年，第三版；附录（设计流程图、程序、表格、数据等）4.机算潮流程序及结果// dierti.cpp : 定义控制台应用程序的入口点。
//#include "stdafx.h"struct Line //线路结构体{int Num,NumI,NumJ; //线路号左节点名右节点名float R,X,B,K; //电阻电抗电纳变比（K等于1为普通支路，不等于1为变压器支路的变比） };struct Bus //节点结构体{int Num ;float Volt,Phase,GenP,GenQ,LoadP,LoadQ;int Type;};#include"stdio.h"#include"string.h"#include"math.h"#include"stdlib.h"#define NBUS 4#define NLINE 4/* Global variables */int nL,nB,nVA,nSH;float X[NBUS];int L;double def[2*NBUS];double mn[50];void Gauss(double a[50][50],double b[50], int n) /*定义高斯法 */{int JS[50];int i,j,k;float d,t,x[50];FILE *fp;int L=1;for(i=0;i<50;i++) JS[i]=0;for(k=0;k<n;k++){d=0.0;for(j=k;j<n;j++)if(fabs(a[k][j])>d){ /*在一行中找到一个最大值赋值d,并用JS[K]记住这个最大值所在的列号*/ d=fabs(a[k][j]);JS[k]=j;}if(fabs(d)<0.000001) /*如果d的数值太小，做为被除数将带来很大的误差 */L=0;else {if(JS[k]!=k)for(i=0;i<n;i++){t=a[i][k];a[i][k]=a[i][JS[k]]; /*进行列交换，让最大值始终在对角元上*/a[i][JS[k]]=t;}}if(L==0)break;for(j=k+1;j<n;j++)a[k][j]=a[k][j]/a[k][k]; /*对角元上的元素消为1*/b[k]=b[k]/a[k][k];for(i=k+1;i<n;i++){for(j=k+1;j<n;j++)a[i][j]=a[i][j]-a[i][k]*a[k][j]; /*使下三角阵的元素为0*/b[i]=b[i]-a[i][k]*b[k];}}if(fabs(a[n-1][n-1])>0.00001){ /*用追赶法，解方程组，求未知数x*/ x[n-1]=b[n-1];for(i=n-2;i>=0;i--){t=0.0;for(j=i+1;j<n;j++)t=t+a[i][j]*x[j];x[i]=(b[i]-t);}}if((fp=fopen("gauss.txt","w"))==NULL) /*将结果写到TXT文件中*/{printf("err");exit(0);}for(i=0;i<n;i++){fprintf(fp,"%lf",x[i]);mn[i]=x[i];fprintf(fp,"\n");}fclose(fp);if(fp!=NULL) fclose(fp);}int _tmain(int argc, _TCHAR* argv[]){FILE *fp;FILE *fpout;int i,j,k,l,h,n,v;int i1,i2,i3,kp,kq;float d1,d2,d3,d4,d5,d6,r,x,g,b,tt,LL,e,ps,qs,shsh,m;struct Line sL[NLINE];struct Bus sB[NBUS];float YG[NBUS+1][NBUS+1],YB[NBUS+1][NBUS+1];double u[50][2];i1=i2=i3=0;d1=d2=d3=d4=d5=d6=ps=qs=0.0;for(i=0;i<NBUS;i++)if((fp=fopen("in.txt","r"))==NULL){ printf("Can not open the file named 'in.txt' \n");exit(0);}fscanf(fp,"%d,%d,%d",&nB,&nL,&nSH);for(i=0;i<nB;i++){sB[i].Num=sB[i].Type=0;sB[i].Volt=1.0;sB[i].Phase=sB[i].GenP=sB[i].GenQ=sB[i].LoadP=sB[i].LoadQ=0.0;fscanf(fp,"%d,%f,%f,%f,%f,%f,%f,%d",&i1,&d1,&d2,&d3,&d4,&d5,&d6,&i2);sB[i].Num=i1;sB[i].Volt=d1;sB[i].Phase=d2;sB[i].GenP=d3;sB[i].GenQ=d4;sB[i].LoadP=d5,sB[i].LoadQ=d6;sB[i].T ype=i2;};for(i=0;i<nL;i++){sL[i].Num=sL[i].NumI=sL[i].NumJ=0;sL[i].R=sL[i].X=sL[i].B=0.0;sL[i].K=1.0;fscanf(fp,"%2d %3d %3d %f %f %f %f",&i1,&i2,&i3,&d1,&d2,&d3,&d4);sL[i].Num=i1;sL[i].NumI=i2;sL[i].NumJ=i3;sL[i].R=d1;sL[i].X=d2;sL[i].B=d3;sL[i].K=d4;}if(fp!=NULL) fclose(fp);/*Make Y Matrix*/for(i=1;i<nB+1;i++)for(j=1;j<nB+1;j++){YG[i][j]=0.0;YB[i][j]=0.0;};for(l=0; l<nL; l++){i=sL[l].NumI;j=sL[l].NumJ;r=sL[l].R;x=sL[l].X;d1=r*r+x*x;g=r/d1;b=-x/d1;m=sL[l].K;if(fabs(sL[l].K-1.0)<0.000001) //普通支路 {YG[i][i]=YG[i][i]+g;YG[j][j]=YG[j][j]+g;YB[i][i]=YB[i][i]+b+sL[l].B;YB[j][j]=YB[j][j]+b+sL[l].B;YG[i][j]=YG[i][j]-g;YG[j][i]=YG[j][i]-g;YB[i][j]=YB[i][j]-b;YB[j][i]=YB[j][i]-b;}else //变压器支路{YG[i][i]=YG[i][i]+g/m+g*(m-1)/m;YG[j][j]=YG[j][j]+g/m+g*(1-m)/m/m;YB[i][i]=YB[i][i]+b/m+b*(m-1)/m;YB[j][j]=YB[j][j]+b/m+b*(1-m)/m/m;YG[i][j]=YG[i][j]-g/m;YG[j][i]=YG[j][i]-g/m;YB[i][j]=YB[i][j]-b/m;YB[j][i]=YB[j][i]-b/m; }}/* Check the Y matrix */if((fp=fopen("GGBB.txt","w"))==NULL){printf("Can not open the file named 'GGBB.txt' \n");exit(0);}fprintf(fp,"---Y Matrix---\n");for(i=1;i<nB+1;i++)for(j=1;j<nB+1;j++)if(fabs(YB[i][j]-0.0)>0.000001) fprintf(fp,"Y(%3d,%-3d)=(%10.5f,%10.5f)\n",i,j,YG[i][j],YB[i][j]);if(fp!=NULL) fclose(fp);/* 节点电压附初值 */for(i=1;i<nB+1;i++){if(sB[i-1].Type==0){u[i][0]=0.0;u[i][1]=1.0;}else if(sB[i-1].Type==1){u[i][1]=sB[i-1].Volt;u[i][0]=0.0;}else if(sB[i-1].Type==2){u[i][1]=sB[i-1].Volt;u[i][0]= sB[i-1].Phase;}}for(v=1;;v++)/* 迭代次数可以无限大 */{/* 节点电压附初值 */printf("迭代第%d次赋予的电压初值为e+jf:\n",v); for(i=1;i<nB+1;i++)printf("%lf,%lf\n",u[i][1],u[i][0]);printf("\n");printf("\n");/* 求偏移量 */double P_P[10];double P_Q[10];double P_UU[10];for(i=1;i<nB+1;i++){if(sB[i-1].Type==2){P_P[i]=0.0;P_Q[i]=0.0;P_UU[i]=1.05;}if(sB[i-1].Type==0){double tempP=0.0;double tempQ=0.0;for(j=1;j<nB+1;j++){tempP+=YG[i][j]*u[j][1]-YB[i][j]*u[j][0];tempQ+=YG[i][j]*u[j][0]+YB[i][j]*u[j][1];}P_P[i]=(sB[i-1].GenP-sB[i-1].LoadP)-tempP*u[i][1]-tempQ*u[i][0]; P_Q[i]=(sB[i-1].GenQ-sB[i-1].LoadQ)-tempP*u[i][0]+tempQ*u[i][1]; P_UU[i]=0.0;}if(sB[i-1].Type==1){double tempP=0.0;double tempQ=0.0;for(j=1;j<nB+1;j++){tempP+=YG[i][j]*u[j][1]-YB[i][j]*u[j][0];tempQ+=YG[i][j]*u[j][0]+YB[i][j]*u[j][1];P_P[i]=(sB[i-1].GenP-sB[i-1].LoadP)-tempP*u[i][1]-tempQ*u[i][0]; }P_UU[i]=sB[i-1].Volt*sB[i-1].Volt-u[i][1]*u[i][1]-u[i][0]*u[i][0]; P_Q[i]=0.0;}}/* 偏移量阵 */double P_PQ[6];int a=0;for(i=1;i<3;i++){P_PQ[a]=P_P[i];a=a+2;}a=1;for(i=1;i<3;i++){P_PQ[a]=P_Q[i];a=a+2;P_PQ[4]=P_P[3];P_PQ[5]=P_UU[3];printf("迭代第%d次的偏移量为：\n",v);for(i=0;i<6;i++){printf("%f",P_PQ[i]);printf("\n");}printf("\n");printf("\n");/* 雅可比矩阵 */double H[6][6],N[6][6],J[6][6],L[6][6],R[6][6],S[6][6],aa[6],bb[6]; for(i=1;i<5;i++){ if(fabs(sB[i-1].Type-2.0)<0.000001)continue;else{for(j=1;j<5;j++)if(i!=j){H[i][j]=-YB[i][j]*u[i][1]+YG[i][j]*u[i][0];N[i][j]=YG[i][j]*u[i][1]+YB[i][j]*u[i][0];J[i][j]=-N[i][j];L[i][j]=H[i][j];R[i][j]=0;S[i][j]=0;}else{aa[i]=bb[i]=0.0;aa[i]+=YG[i][n]*u[n][1]-YB[i][n]*u[n][0];bb[i]+=YG[i][n]*u[n][0]+YB[i][n]*u[n][1];}H[i][i]=-YB[i][i]*u[i][1]+YG[i][i]*u[i][0]+bb[i]; N[i][i]=YG[i][i]*u[i][1]+YB[i][i]*u[i][0]+aa[i]; J[i][i]=-YG[i][i]*u[i][1]-YB[i][i]*u[i][0]+aa[i]; L[i][i]=YG[i][i]*u[i][0]-YB[i][i]*u[i][1]-bb[i]; R[i][i]=2*u[i][0];S[i][i]=2*u[i][1];}}}double ss[50][50];for(i=0;i<6;i++)for(j=0;j<6;j++)ss[i][j]=0.0;for(i=1;i<3;i++)for(j=1;j<4;j++){ss[2*i-2][2*j-2]=H[i][j];ss[2*i-2][2*j-1]=N[i][j];ss[2*i-1][2*j-2]=J[i][j];ss[2*i-1][2*j-1]=L[i][j];}i=3;for(j=1;j<4;j++){ss[2*i-2][2*j-2]=H[i][j];ss[2*i-2][2*j-1]=N[i][j];ss[2*i-1][2*j-2]=R[i][j];ss[2*i-1][2*j-1]=S[i][j];}printf("迭代第%d次的雅可比矩阵为：\n",v);for(i=0;i<6;i++){for(j=0;j<6;j++)printf("%10f",ss[i][j]);printf("\n");}printf("\n");printf("\n");Gauss(ss,P_PQ,6);for(i=1;i<nB;i++){u[i][0]=u[i][0]+mn[2*(i-1)];u[i][1]=u[i][1]+mn[2*i-1];}double max;max=fabs(P_PQ[0]);for(i=0;i<=5;i++)if (max<fabs(P_PQ[i]))max=fabs(P_PQ[i]);if(fabs(max)<0.0001){printf("满足精度要求，迭代终止，迭代次数为%d\n",v); printf("\n");}/* 叠代循环的括号 */printf("最终求得的节点电压值为e+jf:\n");for(i=1;i<nB+1;i++)printf("%lf,%lf\n",u[i][1],u[i][0]);printf("\n");printf("\n");double uu[5],Phase[5];for(i=1;i<nB+1;i++){uu[i]=sqrt(u[i][1]*u[i][1]+u[i][0]*u[i][0]); Phase[i]=atan(u[i][0]/u[i][1]);}for(i=1;i<nB+1;i++)printf("%lf,%lf\n",uu[i],Phase[i]);*计算线路功率和平衡节点 PV节点功率*/double P[5],Q[5];double tempP=0.0;double tempQ=0.0;for(i=1;i<nB+1;i++){for(j=1;j<nB+1;j++){tempP+=YG[i][j]*u[j][1]-YB[i][j]*u[j][0]; tempQ+=YG[i][j]*u[j][0]+YB[i][j]*u[j][1];}P[i]=tempP*u[i][1]+tempQ*u[i][0];Q[i]=tempP*u[i][0]-tempQ*u[i][1];}for(i=1;i<nB+1;i++)printf("节点%d注入功率为%lf,%lf\n",i,P[i],Q[i]);/* 支路功率 */double V[4][2];for(i=1;i<5;i++)for(j=0;j<3;j++)V[i][j]=u[i][j];double sP[5][5],sQ[5][5];double dsq,dsp,dp,sumgen;for(i=1;i<NBUS+1;i++){for(j=1;j<NBUS+1;j++){sP[i][j]=0.0;sQ[i][j]=0.0;}}for(l=0; l<nL; l++){i=sL[l].NumI;j=sL[l].NumJ;r=sL[l].R;x=sL[l].X;d1=r*r+x*x;if(fabs(sL[l].K-1.0)<0.000001){/*Normal lines or transformers*/sP[i][j]=V[i][1]*V[i][1]*g-V[i][1]*V[j][1]*(g*cos(V[i][0]-V[j][0])+b*sin(V[i][0]-V[j][0]));sQ[i][j]=-(V[i][1]*V[i][1]*sL[l].B+V[i][1]*V[i][1]*b+V[i][1]*V[j][1]*(g*sin(V[i][0]-V[j][0])-b*cos(V[i ][0]-V[j][0])));sP[j][i]=V[j][1]*V[j][1]*g-V[i][1]*V[j][1]*(g*cos(V[j][0]-V[i][0])+b*sin(V[j][0]-V[i][0]));sQ[j][i]=-(V[j][1]*V[j][1]*sL[l].B+V[j][1]*V[j][1]*b+V[i][1]*V[j][1]*(g*sin(V[j][0]-V[i][0])-b*cos(V[j ][0]-V[i][0])));}else{/*abnormal transformer ratio*/sP[i][j]=V[i][1]*V[i][1]*g/sL[l].B/sL[l].B-V[i][1]*V[j][1]*(g*cos(V[i][0]-V[j][0])/sL[l].B+b*sin(V[i][ 0]-V[j][0])/sL[l].B);sQ[i][j]=-(V[i][1]*V[i][1]*b/sL[l].B/sL[l].B+V[i][1]*V[j][1]*(g*sin(V[i][0]-V[j][0])/sL[l].B-b*cos(V[i ][0]-V[j][0])/sL[l].B));sP[j][i]=V[j][1]*V[j][1]*g-V[i][1]*V[j][1]*(g*cos(V[j][0]-V[i][0])/sL[l].B+b*sin(V[j][0]-V[i][0])/sL[l ].B);sQ[j][i]=-(V[i][1]*V[i][1]*b+V[i][1]*V[j][1]*(g*sin(V[j][0]-V[i][0])/sL[l].B-b*cos(V[j][0]-V[i][0])/sL [l].B));}}/* 输电效率 */dsp=P[4];sumgen=P[4];for(i=0;i<NBUS;i++){dsp+=sB[i].GenP-sB[i].LoadP;dsq+=sB[i].GenQ-sB[i].LoadQ;sumgen+=sB[i].GenP;}dp=dsp/sumgen*100;/* 输出功率情况 */if((fp=fopen("功率情况.txt","w"))==NULL){printf("Can not open the file named '功率情况.txt' \n");exit(0);}fprintf(fp,"---功率情况---\n");fprintf(fp,"平衡节点功率S=%10.5f+ j%10.5f\n",P[4],Q[4]);for(i=1;i<NBUS+1;i++)for(j=1;j<NBUS+1;j++)if(fabs(sP[i][j]-0.0)>0.000001)fprintf(fp,"S(%3d,%-3d)=(%10.5f,j%10.5f)\n",i,j,sP[i][j],sQ[i][j]); fprintf(fp,"网损为%10.5f+j%10.3f,输电效率为%10.3f\n",dsp,dsq,100-dp);if(fp!=NULL) fclose(fp);return 0;}结果：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 )=( -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．设定电压初值01.1;01;01)0(3)0(3)0(2)0(2)0(1)0(1j jf e j jf e j jf e +=++=++=+ 3．计算功率和电压偏移；27731.0])()([41)0(11)0(1)0(141)0(11)0(1)0(11)0(11)0(1-=++--=-=∆∑∑==j j j j j j jj s s e B f G ff B eG e P P P P0.05097])()([41)0(11)0(1)0(141)0(11)0(1)0(11)0(11)0(1-=----=-=∆∑∑==jj j j jj jj s s e B f G ef B e G f Q Q Q Q同理可算出52596.0)0(22)0(2-=-=∆P P P s ,0196.0)0(22)0(2=-=∆Q Q Q s 5.0)0(33)0(3=-=∆P P P s ,0.02)0(3232)0(3=-=∆U U U s 4．根据求的第一次迭代时雅可比矩阵各元素的公式计算雅可比矩阵各个元素的具体值：⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤-----------⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡20000.200066667.300003334.40052691.406629.130208.256148.00001821.182612.456148.030208.266667.3030208.256148.006298.803803.1066667.356148.030208.299265.032104.85．求高斯计算后的修正量：⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∆∆∆∆∆∆-0.0000000.1276520.023566-0.108546-0.006511-0.007919-2)0(3)0(3)0(2)0(2)0(1)0(11)0()0(3)0(3)0(2)0(2)0(1)0(1U P Q P Q P e f e f e f J 6．计算各节点电压的一次近似值：12765.010855.000792.010000.197643.099349.0)0(3)0(3)1(3)0(2)0(2)1(2)0(1)0(1)1(1)0(3)0(3)1(3)0(2)0(2)1(2)0(1)0(1)1(1=∆+=-=∆+=-=∆+==∆+==∆+==∆+=f f f f f f fffe e e e e e e e e返回第三步重新迭代，并校验收敛与否，令410-=ε。
(完整word版)电力系统暂态分析课程设计
电力系统的短路计算能够帮助我们避免很多不必要的损失，随着科学技术的发展，电力已经和人们的生活密切相关，而建立结构合理的大型电力系统不仅便于电能生产与消费的集中管理、统一调度和分配，减少总装机容量，节省动力设施投资，且有利于地区能源资源的合理开发利用，更大限度地满足地区国民经济日益增长的用电需要。
电力系统的规模和技术水准已成为一个国家经济发展水平的标志之一。
电力系统的短路故障是严重的，而又是发生几率最多的故障，当发生短路时，其短路电流可达数万安以至十几万安，它们所产生的热效应和电动力效应将使电气设备遭受严重破环。
为使非故障部分从不正常运行情况下解脱出来，这要求电气设备必须有足够的机械强度和热稳定度，开关电气设备必须具备足够的开断能力。
因此，电力系统短路电流计算是电力系统运行分析，设计计算的重要环节，许多电业设计单位和个人倾注极大精力从事这一工作的研究。
由于电力系统结构复杂，随着生产发展，技术进步系统日趋扩大和复杂化，短路电流计算工作量也随之增大，采用计算机辅助计算势在并行。
我们这次课程设计使用PSCAD软件对电力系统进行线路建模，从而计算出不同短路类型时的短路电流及其波形。
关键字：短路计算；pscad软件使用；电力系统建模。
目录摘要1 课题概述 (1)1.1课题要求 (1)1.2课题内容 (5)1.3 课题目的 (6)2短路故障分析 (6)2.1不对称故障的分析 (6)2.2三相短路故障分析 (7)2.3短路电流计算步骤 (9)3方案设计 (10)3.1方案概述 (10)3.2 课题设计图 (10)3.3 课题步骤图 (11)4两相短路的仿真分析 (12)4.1 PSCAD简介 (12)4.2 两相短路故障的仿真 (13)5总结 (15)参考文献 (16)1课题概述1.1 课题要求（1）通过课程设计是学生掌握电力系统三相短路计算的基本原理和方法；（2）掌握并能熟练运用PSCAD/MATLAB仿真软件；（3）建立系统三相接线图的仿真过程；（5）编写短路计算流程图；（4）得出仿真结果。
新版电力系统暂态上机计算课程设计报告完整无误版附程序
课程设计报告( 2014—2015年度第二学期)名称：电力系统暂态上机计算院系：电气与电子工程学院班级：电气1211学号：1121210205学生姓名：郝阳指导教师：陈艳波设计周数：两周成绩：日期：2015年7月4日一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。
二、设计正文（详细内容见附录,用A4纸，页数不限）1．对称短路计算过程流程图和计算结果2．不对称短路计算过程流程图和计算结果3．静态稳定性计算过程流程图和计算结果4．暂态稳定性计算过程流程图和计算结果5．思考题三、课程设计总结或结论本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。
通过本次课程设计，本人对电力系统故障分析有了更深刻的理解，电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。
在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。
在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。
同时我对用来分析电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。
这与手算系统短路电流时使用的网络化简方法大大不同。
在学习中，参照潮流程序，我加深了对节点导纳矩阵建立方法的理解与学习，巩固了不同类型短路的短路电流计算方法，和序电压、序电流，相电压、相电流基于matlab软件的计算。
利用小干扰分析法判别静态稳定可以使用劳斯判据，也可以使用特征根判别的方法。
二者都能够判断，在书中的例题中使用的是劳斯判据，但是，特征根判别放法更适合在matlab软件的环境下使用，因为该软件提供了计算矩阵特征根的函数，使用起来简便易懂。
故障上机计算 电力系统故障分析课程设计
课程设计(综合实验)报告( 2010 -- 2011年度第二学期)名称：电力系统故障分析题目：故障上机计算院系：电力工程系班级：电自学号：2008学生姓名：李博指导教师：余洋设计周数： 2 周成绩：日期：年月一、概论：电力系统故障分析课程设计是“电力系统故障分析”课程的实践性教学环节，具有很强的实践性和现实意义。
通过这次的课程设计可以使我们对用计算机计算电力系统故障的原理和过程有更为清晰和深刻的理解，同时进一步掌握电力系统故障计算的基本知识和方法，提高动手编制计算程序的能力，为今后进行电力系统故障计算打下初步的基础。
而随着计算机技术在电力系统中应用的日益普及，这次的实践对自己在日后的工作和学习中运用计算机解决实际问题有很大的帮助。
经过约一周时间的编程和调试，在指导老师和同学的帮助下，编制出了电力系统短路故障计算程序，完成了课程设计指定的各项基本要求任务。
二、任务：1、 对于给定的系统接线图，按照要求手工计算出各种简单故障的电气量值；2、 用C 语言或相关的语言编制一个电力系统短路故障的计算程序并上机调试，要求该程序能处理各种简单的故障，并将计算结果以文件形式输出。
三、任务设计阶段：（一） 计算的假设1、 不计各元件的电阻，即R=0；2、 各元件的正序和负序电抗相等，即X1=X2；3、 不计线路电容电流的影响；4、 计算简单短路故障时，各发电机电势的标幺值均为E=j1。
5、 各种短路均为金属性短路，纵向不对称故障指一相断相和两相断相；6、 故障形式取A 相为特殊相，故障计算以A 相为基准相。
（二） 准备工作1、 始数据计算根据给定的系统原始数据计算系统各元件的标幺值参数，取av B B U U MVA S ==,1000。
（1）发电机参数①计算公式：X "d = X "dNB S S ⎪⎪⎭⎫ ⎝⎛B N U U 2(2)变压器①计算公式：X T=⋅100k U N B S S , )(21g kj g ki j ki ki U U U U ----+= (i,j,k,=1,2,3)②计算结果：(3) 线路① 计算公式：114.0L X = ，211BBU S X X ⋅=*，KV U B 115= ，103X X =② 计算结果：（4）系统系统Ⅰ：X1=0.36 X0=0.87 系统Ⅱ：X1=0.38 X0=0.512、序网络的制定根据上面的基本数据可以制定出系统的正序、负序、零序网络图。
新版电力系统暂态上机计算课程设计报告完整无误版附程序
课程设计报告( 2014—2015年度第二学期)名称：电力系统暂态上机计算院系：电气与电子工程学院班级：电气1211学号：1121210205学生姓名：郝阳指导教师：陈艳波设计周数：两周成绩：日期：2015年7月4日一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。
二、设计正文（详细内容见附录,用A4纸，页数不限）1．对称短路计算过程流程图和计算结果2．不对称短路计算过程流程图和计算结果3．静态稳定性计算过程流程图和计算结果4．暂态稳定性计算过程流程图和计算结果5．思考题三、课程设计总结或结论本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。
通过本次课程设计，本人对电力系统故障分析有了更深刻的理解，电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。
在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。
在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。
同时我对用来分析电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。
这与手算系统短路电流时使用的网络化简方法大大不同。
在学习中，参照潮流程序，我加深了对节点导纳矩阵建立方法的理解与学习，巩固了不同类型短路的短路电流计算方法，和序电压、序电流，相电压、相电流基于matlab软件的计算。
利用小干扰分析法判别静态稳定可以使用劳斯判据，也可以使用特征根判别的方法。
二者都能够判断，在书中的例题中使用的是劳斯判据，但是，特征根判别放法更适合在matlab软件的环境下使用，因为该软件提供了计算矩阵特征根的函数，使用起来简便易懂。
电力系统暂态上机计算报告
电力系统暂态上机计算报告系别：班级：姓名：FREEMAN学号：2011年1月三、给定题目的程序计算结果：1．电磁暂态输入文件：input.txt2 100,10.5,0.18 120,0.11,1 200,10.5,0.14 200,0.10,231,2,110,0.41,3,60,0.42,3,60,0.42,3,40100115 输出文件：output.txt------------------------------------原始数据-----------------------------------发电机G1: SGN1= 100 MV A, UGN1=10.50 KV , xd1=0.180变压器T1: STN1= 120 MV A, 变比为10.50 / 115.00 KV , Uk1%=11.00000 ,所连节点号为1电源支路数发电机的SN 、发电机的UN 、发电机的次暂态电抗 发电机相连的变压器的SN 、变压器的电抗百分数、发电机支路所接节点 节点总数 支路所接的两个节点号,线路长度,支路单位电抗 故障所在支路的两节点号，距第一个点的距离百分数 功率基准值 线路电压的基准值发电机G2: SGN2= 200 MV A, UGN2=10.50 KV, xd2=0.140变压器T2: STN2= 200 MV A, 变比为10.50 / 115.00 KV, Uk2%=10.00000 ,所连节点号为2输电线路1-2: 长110 KM, 单位电抗0.400 Ω/km输电线路1-3: 长60 KM, 单位电抗0.400 Ω/km输电线路2-3: 长60 KM, 单位电抗0.400 Ω/km故障情况：在线路2-3上距2节点40%处发生三相短路---------------------------------------------------------------------------------------------------------------本程序计算结果为:--------------------------------计算后线路电抗为：线路1-2 :0.332703, 线路1-3 :0.181474, 线路2-4 :0.072590, 线路3-4 :0.108885,节点导纳矩阵为：j(-12.197080) j( 3.005682) j( 5.510417) j( 0.000000)j( 3.005682) j(-25.115057) j( 0.000000) j(13.776042)j( 5.510417) j( 0.000000) j(-14.694444) j( 9.184028) j( 0.000000) j(13.776042) j( 9.184028) j(-22.960069)求得的R和D因子为:j(0.081987) j(-0.246426) j(-0.451782) j(0.000000) 0.00000 j(0.041027) j(-0.055711) j(-0.565185) 0.00000 0.00000 j(0.082445) j(-0.820452) 0.00000 0.00000 0.00000 j(0.142667)所求得电压、电流结果为：短路电流:相量值为if =j( -7.009324)有名值为If = 3.518987 KA各节点电压为:相量值为u1 = 0.478795, 有名值为U1 = 55.061387 KV相量值为u2 = 0.389107, 有名值为U2 = 44.747275 KV相量值为u3 = 0.179548, 有名值为U3 = 20.648020 KV各线路故障电流为:相量值i1-2 = j(-0.269573) , 有名值I1-2 = 0.135338 KA相量值i1-3 = j(-1.648974) , 有名值I1-3 = 0.827857 KA相量值i2-f = j(-5.360351) , 有名值I2-f = 2.691130 KA相量值i3-f = j(-1.648974) , 有名值I3-f = 0.827857 KA相量值△iG1 = j(-1.918547) , 有名值△IG1 = 10.549274 KA, 相量值△iG2 = j(-5.090777) , 有名值△IG2 = 27.992015 KA, -----------------------------------END----------------------------------------2．机电暂态a.输入文件：input.txt200,115240,10.5,0.80,1.1,0.9,0.3,0.44,1.54300,242,10.5,12.0280,220,121,14.0230,0.4,0.4,1.2250,0.98,1152,0.65,0容量基准值，电压基准值发电机容量，电压，功率因数，直轴电抗，交轴电抗，暂态电抗，负序电抗，惯性时间变压器T1容量，高压侧电压，低压侧电压，Uk%变压器T2容量，高压侧电压，低压侧电压，Uk%输电线路长度，单位长度正序、负序、零序电抗正常情况下向无穷大系统输出的功率、功率因数、无穷大系统电压故障类型，故障点距线路始端百分数，接地电抗b.输出文件：output.txt******************* 原始数据****************容量基准值:200.000000 电压基准值:115.000000----------------------------------------------------------------------------------- 发电机参数：额定有功功率:240.000000 额定电压:10.500000 功率因数:0.800000 直轴电抗:1.100000 交轴电抗:0.900000 直轴暂态电抗:0.300000 负序电抗:0.440000 Tj:1.540000 ----------------------------------------------------------------------------------- 变压器参数：变压器1容量:300.000000 高压侧额定电压:242.000000 低压侧额定电压电压:10.500000 电压百分数:12.000000变压器2容量:280.000000 高压侧额定电压:220.000000 低压侧额定电压电压:121.000000 电压百分数:14.000000----------------------------------------------------------------------------------- 输电线路参数：线路长度:230.000000 单位正序电抗:0.400000 单位负序电抗:0.400000 单位零序电抗=1.200000----------------------------------------------------------------------------------- 无穷大系统参数：故障前吸收有功:250.000000 功率因数:0.980000 恒定电压:115.000000----------------------------------------------------------------------------------- 故障地点：距线路始端0.650000, 故障类型：两相短路，故障点接地电抗为：j0.000000*************** 原始数据标么值归算************** 容量基准值:200.000000 电压基准值:115.000000 220KV电压等级电压基准值:209.090912----------------------------------------------------------------------------------- 发电机参数：额定有功功率:1.200000 额定电压:1.157391 功率因数:0.800000 直轴电抗:0.982340 交轴电抗:0.803733直轴暂态电抗:0.267911 负序电抗:0.392936 Tj:2.310000----------------------------------------------------------------------------------- 变压器参数：变压器T1容量:1.500000 高压侧额定电压:242.000000 低压侧额定电压电压:10.500000 电抗:0.107164 变压器T2容量:1.400000 高压侧额定电压:220.000000 低压侧额定电压电压:121.000000 电抗:0.110707----------------------------------------------------------------------------------- 输电线路参数：线路长度:230.000000 单位正序电抗:0.420870 单位负序电抗:0.420870 单位零序电抗=1.262609-----------------------------------------------------------------------------------无穷大系统参数：故障前吸收有功:1.250000 功率因数:0.980000恒定电压:1.000000-----------------------------------------------------------------------------------故障点接地电抗：0.000000****************** 正常运行时参数***************直轴暂态电抗和=0.696217交轴电抗和=1.232039虚构电动势E'=1.463569正常运行时候的功角=36.485753度。
新版电力系统暂态上机计算课程设计报告完整无误版
课程设计报告（2014 —2015年度第二学期）名称：电力系统暂态上机计算院系：电气与电子工程学院班级：电气1211学号：1121210205学生姓名：_________ 郝阳_________指导教师：陈艳波设计周数：_________ 两周成绩：__________________________日期：2015年7月4日、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。
二、设计正文（详细内容见附录，用A4纸，页数不限）1．对称短路计算过程流程图和计算结果2．不对称短路计算过程流程图和计算结果3．静态稳定性计算过程流程图和计算结果4．暂态稳定性计算过程流程图和计算结果5．思考题三、课程设计总结或结论本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。
通过本次课程设计，本人对电力系统故障分析有了更深刻的理解，电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。
在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。
在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。
同时我对用来分析电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。
这与手算系统短路电流时使用的网络化简方法大大不同。
在学习中，参照潮流程序，我加深了对节点导纳矩阵建立方法的理解与学习，巩固了不同类型短路的短路电流计算方法，和序电压、序电流，相电压、相电流基于matlab 软件的计算。
利用小干扰分析法判别静态稳定可以使用劳斯判据，也可以使用特征根判别的方法。
电力系统潮流上机课程设计报告
课程设计报告一、课程设计的目的与要求培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识二、设计正文（详细内容见附录）总体步骤如下：读入相关数据 -> 构造节点导纳阵Y -> 计算ΔP,ΔQ，构造Δf -> 构造雅克比矩阵 -> 计算潮流分布（包括迭代修正数据） -> 输出结果3 2 1母线数据1 3 1.0 0 0 0 0 02 1 1.0 0 0.8 0.6 0 03 2 1.1 0 0 0 0.4 0支路数据1 1 3 1 0.0497 0.2000 0 02 1 2 1 0.0100 0.0400 0 0接地支路1 3 0 0.33迭代数量：iteration = 3计算结果：电压相角（弧度）相角1 0 00.96653 -0.0269036 -1.541461.1 0.05205522.98255母线上的P和Q注入：编号 P Q1 0.425894 0.2458942 -0.8 -0.63 0.4 0.0587502线路潮流分布：编号 I J Sij Sji1 1 3 -0.38481+j-0.396925 -0.139061+j0.5920072 1 2 0.810705+j0.642819 0.423529+j-0.905883功率损耗：Ploss = 0.0258945 Qloss = -0.2953569 9 0母线数据1 3 1.04 0. 0. 0. 71.641 27.045922 2 1.025 0. 0. 0. 163. 6.653663 2 1.025 0. 0. 0. 85. -10.85974 1 1. 0. 0. 0. 0. 0.5 1 1. 0. 125. 50. 0. 0.6 1 1. 0. 90. 30. 0. 0.7 1 1. 0. 0. 0. 0. 0.8 1 1. 0. 100. 35. 0. 0.9 1 1. 0. 0. 0. 0. 0. 支路数据1 4 12 0. 0.0576 0. 1.2 7 2 2 0. 0.0625 0. 1.3 9 3 2 0. 0.0586 0. 1.4 7 8 1 0.0085 0.072 0.149 05 9 8 1 0.0119 0.1008 0.209 06 7 5 1 0.032 0.161 0.306 07 9 6 1 0.039 0.17 0.358 08 5 4 1 0.01 0.085 0.176 09 6 4 1 0.017 0.092 0.158 0迭代数量：iteration = 4计算结果：电压相角（弧度）相角1.04 0 01.025 0.165578 9.486941.025 0.0833249 4.774171.0105 -0.039445 -2.260050.9727 -0.070866 -4.060350.9890 -0.064734 -3.709031.0113 0.067145 3.847140.9972 0.013673 0.783421.0180 0.0355712.03811母线上的P和Q注入：编号 P Q1 0.71641 0.5468462 1.63 0.3040183 0.85 0.1424344 0 05 -1.25 -0.56 -0.9 -0.37 0 08 -1 -0.359 0 0线路潮流分布：编号 I J Sij Sji1 4 1 -0.719499+j-0.503352 -0.719499+j0.5468462 7 2 -1.63+j-0.140465 -1.63+j0.3040183 9 3 -0.85+j-0.101004 -0.85+j0.1424344 7 8 0.763816+j0.0521904 0.717905+j-0.3345885 9 8 0.242231+j0.0755265 0.193826+j-0.3301066 7 5 0.866184+j-0.0280722 0.775227+j-0.4739837 9 6 0.607769+j-0.121426 0.533117+j-0.436418 5 4 -0.407805+j-0.461515 -0.305617+j0.4004949 6 4 -0.306286+j-0.248276 -0.222811+j0.200074功率损耗：Ploss = 0.0494989 Qloss = -0.156702三、课程设计总结或结论本次课程实验有如下的收获：1.学会使用基本的c++语法进行程序编写，掌握了指针与数组的关系与其应用。