《电力系统潮流上机》课程设计报告
院 系：电气与电子工程学院 班 级： 电气1108班 学 号： 1111180806 学生姓名： 陈子君 指导教师： 刘宝柱 设计周数： 两周
成 绩：
日期：2014年1月10日
课程课程设计报告
一、课程设计的目的与要求
培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识
二、设计正文（详细内容见附录）
1.手算： 要求应用牛顿-拉夫逊法或P-Q分解法手算求解，要求精度为0.001MW。
1.00，节点2为PQ节点，负荷功率S0.8j0.6，节节点1为平衡节点，电压U21点3是PV节点，P30.4,U31.1，两条支路分别为Z130.01j0.04，Z120.05j0.2，对地支路y30j0.33
2.计算机计算：编写潮流计算程序，要求如下：
2.1据给定的潮流计算任务书整理潮流计算的基础数据：节点的分类，线路模型，等值变压器
模型，电压等级的归算，标幺值的计算；
2.2基础数据的计算机存储：节点数据，支路数据（包括变压器）； 2.3用牛顿-拉夫逊法计算；
2.4根据所选潮流计算方法画流程图，划分出功能模块，有数据输入模块，导纳阵形成模块，
解线性方程组模块，计算不平衡功率模块，形成雅可比矩阵模块，解修正方程模块，计算线路潮流，网损，PV节点无功功率和平衡节点功率，数据输出模块； 2.5据上述模块编制程序并上机调试程序，得出潮流计算结果； 2.6源程序及其程序中的符号说明集、程序流图
简单系统如下图所示，支路数据如下：
~z120.1j0.41，z13j0.3，z140.12j0.5，z240.08j0.40
y10,2y20,1j0.01528，y10,4y40,1j0.0192，y20,4y40,2j0.01413
k1.1
节点数据如下：
0.30j0.18，S0.55j0.13， S210.5，U1.10，U41.050o S33
院系 课程课程设计报告
班级： z13S学号：1 kS学生姓名：z13 z313y10,4指导教师：k(k 1)1kz14设计周数y10,2 zy20,1y40,1 12成绩： SS4yz22440,2日期：年月日y20,4
// 2014.cpp : 定义控制台应用程序的入口点。 //
#include \"stdafx.h\"
struct Line //线路结构体 {
int Num,NumI,NumJ; //线路号 左节点名 右节点名
为变压器支路的变比）float R,X,B,K; // 电阻 电抗 电纳 变比（K等于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
不等于1
课程课程设计报告
/* Global variables */ int nL,nB,nVA,nSH; float X[NBUS]; int L;
double def[2*NBUS]; double awe[100];
int _tmain(int argc, _TCHAR* argv[]) { FILE *fp; FILE *fpout;
int i,j,l,h,n,k;
int i1=0,i2=0,i3=0,kp,kq; double a,b,k1;
double H=0,N=0,J=0,L=0,R=0,S=0;
float d1=0,d2=0,d3=0,d4=0,d5=0,d6=0,r,x,g,b1,b2,g1,g2,tt,LL,e,ps=0,qs=0,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],p[50],da,db,pp[50],ppl[50],Q[10][10],gJaccobi[50][50],sP[10][10],sQ[10][10],WP=0.0,WQ=0.0;
for(i=0;i fscanf(fp,\"%d,%d,%d\ for(i=0;i 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\ 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].Type=i2; }; 课程课程设计报告 for(i=0;i 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 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; YB[i][i]=YB[i][i]+b; YG[j][j]=YG[j][j]+g/(m*m); YB[j][j]=YB[j][j]+b/(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\ printf(\"Can not open the file named 'GGBB.txt' \\n\");exit(0);} fprintf(fp,\"---Y Matrix---\\n\"); for(i=1;i fprintf(fp,\"Y(%3d,%-3d)=(%10.5f,%10.5f)\\n\if(fp!=NULL) fclose(fp); for(i=1;i fprintf(fp,\"---电压初值---\\n\"); for(i=1;i if(sB[i-1].Type==2) { } else if(sB[i-1].Type==1) { u[i][0]=sB[i-1].Volt*cos(sB[i-1].Phase); } else if(sB[i-1].Type==0) { } u[i][0]=1.0; u[i][1]=0.0; u[i][0]=sB[i-1].Volt*cos(sB[i-1].Phase); u[i][1]=sB[i-1].Volt*sin(sB[i-1].Phase); u[i][1]=0.0; fp=fopen(\"电压初值.txt\ 课程课程设计报告 if(fp!=NULL) fclose(fp); printf(\"赋予的电压初值为(e+if):\\n\"); for(i=1;i /*循环迭代*/ { /*求偏移量*/ double P_P[5]; double P_Q[5]; double P_UU[5]; int v; for(v=1;;v++) for(i=1;i if(sB[i-1].Type==1) { double tempP=0.0; double tempQ=0.0; for(j=1;j double tempP=0.0; for(j=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]; double tempQ=0.0; P_Q[i]=(sB[i-1].GenQ-sB[i-1].LoadQ)+tempP*u[i][0]-tempQ*u[i][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][0]*u[i][0]+u[i][1]*u[i][1]); 课程课程设计报告 }} /*偏移量阵*/ double P_PQ[6];// int c=0; for(i=2;i } for(i=2;i for(i=2;i for(i=2;i /*输出偏移量*/ printf(\"输出偏移量\\n\"); if((fp=fopen(\"偏移量 .txt\ { printf(\"Can not open the file named '偏移量.txt' \\n\"); } fprintf(fp,\"--- 偏移量 ---\\n\"); for(i=0;i exit(0); 课程课程设计报告 printf(\"P_PQ[%2d]=(%10.5f)\\n\ //形成雅克比矩阵 for(i=2;i else //PV节点 { if(i!=j)//非对角元素 { H=-YB[i][j]*u[i][0]+YG[i][j]*u[i][1]; } else//对角元素 N=YG[i][j]*u[i][0]+YB[i][j]*u[i][1]; R=0; S=0; if(i!=j)//非对角元素 { } else//对角元素 { } gJaccobi[2*(i-2)][2*(j-2)]=H; gJaccobi[2*(i-2)][2*(j-2)+1]=N; gJaccobi[2*(i-2)+1][2*(j-2)]=J; gJaccobi[2*(i-2)+1][2*(j-2)+1]=L; a=0;b=0; for(k=1;k a=a+(YG[i][k]*u[k][0]-YB[i][k]*u[k][1]); b=b+(YG[i][k]*u[k][1]+YB[i][k]*u[k][0]); H=-YB[i][j]*u[i][0]+YG[i][j]*u[i][1]; N=YG[i][j]*u[i][0]+YB[i][j]*u[i][1]; J=-N; L=H; 课程课程设计报告 } } } { } gJaccobi[2*(i-2)][2*(j-2)]=H; gJaccobi[2*(i-2)][2*(j-2)+1]=N; gJaccobi[2*(i-2)+1][2*(j-2)]=R; gJaccobi[2*(i-2)+1][2*(j-2)+1]=S; a=b=0; for(k=1;k b=b+(YG[i][k]*u[k][1]+YB[i][k]*u[k][0]); a=a+(YG[i][k]*u[k][0]-YB[i][k]*u[k][1]); //output the Jaccobi fprintf(fp,\"---雅克比矩阵---\\n\"); for(i=0;i<2*(nB-1);i++)// { for(j=0;j<2*(nB-1);j++) { } fprintf(fp,\"\\n\"); fprintf(fp,\"Jaccobi(%d,%d)=(%10.5lf)\\b\雅克比矩阵从i=0开 fp=fopen(\"雅克比.txt\if(fp==NULL) { } printf(\"Can not open the file named '雅克比.txt' \\n\"); exit(0); 始 } if(fp!=NULL) fclose(fp); void Gauss(double a[50][50],double b[6], int n);// int NN=6;// Gauss(gJaccobi,P_PQ,NN);//调用高斯函数 课程课程设计报告 } for(i=0;i<6;i++)// printf(\"\\n开始计算电压新值\\n\"); for(i=2,j=0;i<=nB,j<2*(nB-1);i++,j=j+2)//i=1是平衡节点不放入 { } fp=fopen(\"电压新值.txt\路径*/ if(fp==NULL) { } fprintf(fp,\"---电压新值---\\n\"); for(i=1;i<=3;i++) // { } if(fp!=NULL) fclose(fp); double max; max=0.0; for(i=0;i<6;i++) if(max<=0.00001) printf(\"循环次数v=%d\\n\ break; if(fabs(awe[i])>max) max=fabs(awe[i]); fprintf(fp,\"ge[%d]=%10.5f,gf[%d]=%10.5f\\n\printf(\"Can not open the file named '电压新值.txt' \\n\"); exit(0); u[i][1]+=awe[j]; u[i][0]+=awe[j+1]; printf(\"u[%d][0]=%10.5lf,u[%d][1]=%10.5lf\\n\printf(\"awe[%d]=%lf\\n\ //计算平衡节点功率 double tempP=0.0; double tempQ=0.0; ps=0; qs=0; for(i=1;i 课程课程设计报告 for(j=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]; qs=-tempP*u[i][0]+tempQ*u[i][1]; } } //计算线路流动功率 for(i=1;i } } for(l=0; l b1=sL[l].B; if(fabs(sL[l].K-1.0)<0.000001) {//普通支路 sP[i][j]=u[i][0]*((u[i][0]-u[j][0])*g-(u[i][1]-u[j][1])*b)+u[i][1]*((u[i][1]-u[j][1])*g+(u[i][0]-u[j][0])*b); sQ[i][j]=-u[i][0]*((u[i][0]-u[j][0])*b+(u[i][1]-u[j][1])*g)+u[i][1]*((u[i][0]-u[j][0])*g-(u[i][1]-u[j][1])*b)-(u[i][0]*u[i][0]+u[i][1]*u[i][1])*b1; sP[j][i]=u[j][0]*((u[j][0]-u[i][0])*g-(u[j][1]-u[i][1])*b)+u[j][1]*((u[j][1]-u[i][1])*g+(u[j][0]-u[i][0])*b); sQ[j][i]=-u[j][0]*((u[j][0]-u[i][0])*b+(u[j][1]-u[i][1])*g)+u[j][1]*((u[j][0]-u[i][0])*g-(u[j][1]-u[i][1])*b)-(u[j][0]*u[j][0]+u[j][1]*u[j][1])*b1; } else {//abnormal transformer ratio k1=sL[l].K; r=sL[l].R; x=sL[l].X; d1=r*r+x*x; 课程课程设计报告 g=r/d1; g1=(k1-1)/k1*g; g2=(1-k1)/(k1*k1)*g; b1=(k1-1)/k1*b; b2=(1-k1)/(k1*k1)*b; b=-x/d1; //sP[i][j]=u[i][0]*((u[i][0]-u[j][0])*g-(u[i][1]-u[j][1])*b-u[i][1]*b1)+u[i][1]*((u[i][0]-u[j][0])*b+(u[i][1]-u[j][1])*g+u[i][0]*b1); //sP[j][i]=u[j][0]*((u[j][0]-u[i][0])*g-(u[j][1]-u[i][1])*b-u[j][1]*b2)+u[j][1]*((u[j][0]-u[i][0])*b+(u[j][1]-u[i][1])*g+u[j][0]*b2); //sQ[i][j]=u[i][1]*((u[i][0]-u[j][0])*g-(u[i][1]-u[j][1])*b-u[i][1]*b1)-u[i][0]*((u[i][0]-u[j][0])*b+(u[i][1]-u[j][1])*g+u[i][0]*b1); //sQ[i][j]=u[j][1]*((u[j][0]-u[i][0])*g-(u[j][1]-u[i][1])*b-u[j][1]*b2)-u[j][0]*((u[j][0]-u[i][0])*b+(u[j][1]-u[i][1])*g+u[i][0]*b2); sP[i][j]=u[i][0]*((u[i][0]-u[j][0])*g/k1-(u[i][1]-u[j][1])*b/k1)+u[i][1]*((u[i][1]-u[j][1])*g/k1+(u[i][0]-u[j][0])*b/k1)+(u[i][0]*u[i][0]+u[i][1]*u[i][1])*g1; sQ[i][j]=-u[i][0]*((u[i][0]-u[j][0])*b/k1+(u[i][1]-u[j][1])*g/k1)+u[i][1]*((u[i][0]-u[j][0])*g/k1-(u[i][1]-u[j][1])*b/k1)-(u[i][0]*u[i][0]+u[i][1]*u[i][1])*b1; sP[j][i]=u[j][0]*((u[j][0]-u[i][0])*g/k1-(u[j][1]-u[i][1])*b/k1)+u[j][1]*((u[j][1]-u[i][1])*g/k1+(u[j][0]-u[i][0])*b/k1)+(u[j][0]*u[j][0]+u[j][1]*u[j][1])*g2; sQ[j][i]=-u[j][0]*((u[j][0]-u[i][0])*b/k1+(u[j][1]-u[i][1])*g/k1)+u[j][1]*((u[j][0]-u[i][0])*g/k1-(u[j][1]-u[i][1])*b/k1)-(u[j][0]*u[j][0]+u[j][1]*u[j][1])*b2; } } double dsp,dsq,sumgen,dp; dsp=ps; dsq=qs; sumgen=ps; for(i=0;i dp=dsp/sumgen*100; // 输出功率情况 if((fp=fopen(\"功率情况.txt\ { printf(\"Can not open the file named '功率情况.txt' \\n\"); exit(0); } fprintf(fp,\"---各支路功率情况---\\n\"); fprintf(fp,\"平衡节点功率S=%10.5f+ i%10.5f\\n\ for(i=1;i if(fabs(sP[i][j]-0.0)>0.000001) fprintf(fp,\"S(%3d,%-3d)=(%10.5f,j%10.5f)\\n\ fprintf(fp,\"网损为%10.5f+j%10.5f,输电效率为%10.3f\\n\if(fp!=NULL) fclose(fp); //PV节点无功功率 double pP=0.0; double pQ=0.0; double tq[50]={0}; for(i=1;i for(j=1;j tq[i-1]=-pP*u[i][0]+pQ*u[i][1]; } } //输出PV无功 if((fp=fopen(\"PV无功.txt\ { printf(\"Can not open the file named 'PV无功.txt' \\n\"); exit(0); } fprintf(fp,\"---PV无功---\\n\"); for(i=1;i fprintf(fp,\"Q(%3d)=%10.5f\\n\ } if(fp!=NULL) fclose(fp); //输出最终结果 if((fp=fopen(\"最终结果.txt\ { printf(\"Can not open the file named '最终结果.txt' \\n\"); exit(0); } 课程课程设计报告 fprintf(fp,\"---最终结果---\\n\"); fprintf(fp,\"\\n\\n\\n\"); fprintf(fp,\"一共循环迭代计算%d次\\n\ fprintf(fp,\"\\n\\n\\n\"); fprintf(fp,\"---电压终值---\\n\"); for(i=1;i fprintf(fp,\"平衡节点功率S=%8.5f+j%7.5f\\n\ for(i=1;i { } fprintf(fp,\"\\n\\n\\n\"); fprintf(fp,\"---各节点间功率---\\n\"); for(i=1;i fprintf(fp,\"S(%2d,%-2d)=(%8.5f,j%8.5f)\\n\ fprintf(fp,\"\\n\\n\\n\"); fprintf(fp,\"网损为：%7.5f+j%7.5f。\\n\\n网损率为：百分之%7.5f\\n\ if(fp!=NULL) fclose(fp); } void Gauss(double a[50][50], double b[6], int n) /*定义高斯法 */ // { int JS[50]; int i,j,k; double d,t,x[50]; FILE *fp; int L=1; for(i=0;i { /*在一行中找到一个最大值赋值d,并用JS[K]记住这个最大值所在的列号*/ d=fabs(a[k][j]); JS[k]=j; } return 0; 课程课程设计报告 if(fabs(d)<0.000001) /*如果d的数值太小，做为被除数将带来很大的误差 */ L=0; else { for(i=0;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 b[k]=b[k]/a[k][k]; for(i=k+1;i a[i][j]=a[i][j]-a[i][k]*a[k][j]; /*使下三角阵的元素为0*/ for(j=k+1;j { /*用追赶法，解方程组，求未知数x*/ x[n-1]=b[n-1]; for(i=n-2;i>=0;i--) { t=0.0; for(j=i+1;j x[i]=(b[i]-t); } } if((fp=fopen(\"高斯.txt\ /*将结果写到TXT文件中*/ { } for(i=0;i if(fp!=NULL) fclose(fp); printf(\"Can not open the file named '高斯.txt'\\n\");exit(0); 课程课程设计报告 } 5）程序流图请看附录 6）得到的数据（out.txt） ---Y Matrix--- Y( 1,1 )=( 0.93463, -4.26159) Y( 1,2 )=( -0.48077, 2.40385) Y( 1,3 )=( -0.45386, 1.89107) Y( 2,1 )=( -0.48077, 2.40385) Y( 2,2 )=( 1.04225, -4.67651) Y( 2,3 )=( -0.56148, 2.30208) Y( 3,1 )=( -0.45386, 1.89107) Y( 3,2 )=( -0.56148, 2.30208) Y( 3,3 )=( 1.01534, -9.03909) Y( 3,4 )=( 0.00000, 4.43674) Y( 4,3 )=( 0.00000, 4.43674) Y( 4,4 )=( 0.00000, -4.03340) ---雅克比矩阵1--- Jaccobi(1,1)=( 4.82612)_Jaccobi(1,2)=( 1.01821)_Jaccobi(1,3)=( -2.30208)_Jaccobi(1,4)=( -0.56148)_Jaccobi(1,5)=( 0.00000)_Jaccobi(1,6)=( 0.00000)_ Jaccobi(2,1)=( -1.06629)_Jaccobi(2,2)=( 4.52691)_Jaccobi(2,3)=( 0.56148)_Jaccobi(2,4)=( -2.30208)_Jaccobi(2,5)=( 0.00000)_Jaccobi(2,6)=( 0.00000)_ Jaccobi(3,1)=( -2.30208)_Jaccobi(3,2)=( -0.56148)_Jaccobi(3,3)=( 9.16812)_Jaccobi(3,4)=( 0.99265)_Jaccobi(3,5)=( -4.43674)_Jaccobi(3,6)=( 0.00000)_ Jaccobi(4,1)=( 0.56148)_Jaccobi(4,2)=( -2.30208)_Jaccobi(4,3)=( -1.03803)_Jaccobi(4,4)=( 8.91006)_Jaccobi(4,5)=( 0.00000)_Jaccobi(4,6)=( -4.43674)_ Jaccobi(5,1)=( 0.00000)_Jaccobi(5,2)=( 0.00000)_Jaccobi(5,3)=( -4.88042)_Jaccobi(5,4)=( 0.00000)_Jaccobi(5,5)=( 4.43674)_Jaccobi(5,6)=( 0.00000)_ Jaccobi(6,1)=( 0.00000)_Jaccobi(6,2)=( 0.00000)_Jaccobi(6,3)=( 0.00000)_Jaccobi(6,4)=( 0.00000)_Jaccobi(6,5)=( 0.00000)_Jaccobi(6,6)=( 2.20000)_ ---雅克比矩阵2--- Jaccobi(1,1)=( 4.59086)_Jaccobi(1,2)=( 0.99168)_Jaccobi(1,3)=( -2.18767)_Jaccobi(1,4)=( -0.79848)_Jaccobi(1,5)=( 0.00000)_Jaccobi(1,6)=( 0.00000)_ Jaccobi(2,1)=( -2.06030)_Jaccobi(2,2)=( 4.31870)_Jaccobi(2,3)=( 0.79848)_Jaccobi(2,4)=( -2.18767)_Jaccobi(2,5)=( 0.00000)_Jaccobi(2,6)=( 0.00000)_ Jaccobi(3,1)=( -2.28417)_Jaccobi(3,2)=( -0.57680)_Jaccobi(3,3)=( 9.15929)_Jaccobi(3,4)=( 0.78329)_Jaccobi(3,5)=( -4.41095)_Jaccobi(3,6)=( -0.03580)_ Jaccobi(4,1)=( 0.57680)_Jaccobi(4,2)=( -2.28417)_Jaccobi(4,3)=( -1.38147)_Jaccobi(4,4)=( 8.79742)_Jaccobi(4,5)=( 0.03580)_Jaccobi(4,6)=( -4.41095)_ Jaccobi(5,1)=( 0.00000)_Jaccobi(5,2)=( 0.00000)_Jaccobi(5,3)=( -4.88042)_Jaccobi(5,4)=( 0.46062)_Jaccobi(5,5)=( 4.41095)_Jaccobi(5,6)=( 0.03580)_ Jaccobi(6,1)=( 0.00000)_Jaccobi(6,2)=( 0.00000)_Jaccobi(6,3)=( 0.00000)_Jaccobi(6,4)=( 0.00000)_Jaccobi(6,5)=( 0.20764)_Jaccobi(6,6)=( 2.20000)_ ---雅克比矩阵3--- Jaccobi(1,1)=( 4.57330)_Jaccobi(1,2)=( 0.96137)_Jaccobi(1,3)=( -2.14904)_Jaccobi(1,4)=( -0.79057)_Jaccobi(1,5)=( 0.00000)_Jaccobi(1,6)=( 0.00000)_ 课程课程设计报告 Jaccobi(2,1)=( -2.06173)_Jaccobi(2,2)=( 4.17946)_Jaccobi(2,3)=( 0.79057)_Jaccobi(2,4)=( -2.14904)_Jaccobi(2,5)=( 0.00000)_Jaccobi(2,6)=( 0.00000)_ Jaccobi(3,1)=( -2.26771)_Jaccobi(3,2)=( -0.57317)_Jaccobi(3,3)=( 9.09880)_Jaccobi(3,4)=( 0.77422)_Jaccobi(3,5)=( -4.37942)_Jaccobi(3,6)=( -0.03650)_ Jaccobi(4,1)=( 0.57317)_Jaccobi(4,2)=( -2.26771)_Jaccobi(4,3)=( -1.37895)_Jaccobi(4,4)=( 8.72910)_Jaccobi(4,5)=( 0.03650)_Jaccobi(4,6)=( -4.37942)_ Jaccobi(5,1)=( 0.00000)_Jaccobi(5,2)=( 0.00000)_Jaccobi(5,3)=( -4.85817)_Jaccobi(5,4)=( 0.46601)_Jaccobi(5,5)=( 4.37942)_Jaccobi(5,6)=( 0.03650)_ Jaccobi(6,1)=( 0.00000)_Jaccobi(6,2)=( 0.00000)_Jaccobi(6,3)=( 0.00000)_Jaccobi(6,4)=( 0.00000)_Jaccobi(6,5)=( 0.21007)_Jaccobi(6,6)=( 2.18997)_ ---雅克比矩阵4--- Jaccobi(1,1)=( 4.57309)_Jaccobi(1,2)=( 0.96085)_Jaccobi(1,3)=( -2.14837)_Jaccobi(1,4)=( -0.79043)_Jaccobi(1,5)=( 0.00000)_Jaccobi(1,6)=( 0.00000)_ Jaccobi(2,1)=( -2.06175)_Jaccobi(2,2)=( 4.17693)_Jaccobi(2,3)=( 0.79043)_Jaccobi(2,4)=( -2.14837)_Jaccobi(2,5)=( 0.00000)_Jaccobi(2,6)=( 0.00000)_ Jaccobi(3,1)=( -2.26751)_Jaccobi(3,2)=( -0.57315)_Jaccobi(3,3)=( 9.09805)_Jaccobi(3,4)=( 0.77422)_Jaccobi(3,5)=( -4.37903)_Jaccobi(3,6)=( -0.03657)_ Jaccobi(4,1)=( 0.57315)_Jaccobi(4,2)=( -2.26751)_Jaccobi(4,3)=( -1.37904)_Jaccobi(4,4)=( 8.72825)_Jaccobi(4,5)=( 0.03657)_Jaccobi(4,6)=( -4.37903)_ Jaccobi(5,1)=( 0.00000)_Jaccobi(5,2)=( 0.00000)_Jaccobi(5,3)=( -4.85811)_Jaccobi(5,4)=( 0.46602)_Jaccobi(5,5)=( 4.37903)_Jaccobi(5,6)=( 0.03657)_ Jaccobi(6,1)=( 0.00000)_Jaccobi(6,2)=( 0.00000)_Jaccobi(6,3)=( 0.00000)_Jaccobi(6,4)=( 0.00000)_Jaccobi(6,5)=( 0.21007)_Jaccobi(6,6)=( 2.18995)_ ---电压新值1--- ge[1]= 1.05000,gf[1]= 0.00000 ge[2]= 0.97679,gf[2]= -0.10861 ge[3]= 0.99419,gf[3]= -0.00807 --- 偏移量1 --- P_PQ( 0 )=( -0.52596) P_PQ( 1 )=( 0.01960) P_PQ( 2 )=( -0.27731) P_PQ( 3 )=( -0.05097) P_PQ( 4 )=( 0.50000) P_PQ( 5 )=( 0.00000) ---电压新值2--- ge[1]= 1.05000,gf[1]= 0.00000 ge[2]= 0.96017,gf[2]= -0.10923 ge[3]= 0.98708,gf[3]= -0.00823 --- 偏移量2 --- P_PQ( 0 )=( -0.01331) P_PQ( 1 )=( -0.05511) P_PQ( 2 )=( -0.00119) P_PQ( 3 )=( -0.00253) P_PQ( 4 )=( 0.00268) P_PQ( 5 )=( -0.01078) ---电压新值3--- ge[1]= 1.05000,gf[1]= 0.00000 ge[2]= 0.95987,gf[2]= -0.10924 ge[3]= 0.98699,gf[3]= -0.00824 --- 偏移量3 --- P_PQ( 0 )=( -0.00023) P_PQ( 1 )=( -0.00102) P_PQ( 2 )=( -0.00002) P_PQ( 3 )=( -0.00003) P_PQ( 4 )=( 0.00004) P_PQ( 5 )=( -0.00003) ---电压新值4--- ge[1]= 1.05000,gf[1]= 0.00000 ge[2]= 0.95987,gf[2]= -0.10924 ge[3]= 0.98699,gf[3]= -0.00824 --- 偏移量4 --- P_PQ( 0 )=( -0.00000) P_PQ( 1 )=( -0.00000) P_PQ( 2 )=( -0.00000) P_PQ( 3 )=( -0.00000) P_PQ( 4 )=( 0.00000) P_PQ( 5 )=( -0.00000) 课程课程设计报告 ---PV无功1--- Q( 4)= 0.07555 ---各支路功率情况1--- ---PV无功2--- Q( 4)= 0.08895 ---PV无功3--- Q( 4)= 0.08933 ---各支路功率情况2--- ---PV无功4--- Q( 4)= 0.08933 平衡节点功率S= 0.35372+ i 0.20018 S( 1,2 )=( 0.31110,j 0.11437) S( 1,3 )=( 0.04262,j 0.08581) S( 2,1 )=( -0.30285,j -0.10236) S( 2,3 )=( -0.23384,j 0.02747) S( 3,1 )=( -0.04118,j -0.11994) S( 3,2 )=( 0.23969,j -0.03336) S( 3,4 )=( -0.49732,j -0.02416) S( 4,3 )=( 0.49732,j 0.07555) 网损为 0.00372+j -0.10982,输电效率为 99.565 ---各支路功率情况3--- 平衡节点功率S= 0.36762+ i 0.25677 S( 1,2 )=( 0.32122,j 0.15676) S( 1,3 )=( 0.04639,j 0.10001) S( 2,1 )=( -0.31158,j -0.13731) S( 2,3 )=( -0.23842,j 0.00731) S( 3,1 )=( -0.04456,j -0.13225) S( 3,2 )=( 0.24456,j -0.01128) S( 3,4 )=( -0.50000,j -0.03647) S( 4,3 )=( 0.50000,j 0.08933) 网损为 0.01762+j -0.05323,输电效率为 97.970 平衡节点功率S= 0.36737+ i 0.25587 S( 1,2 )=( 0.32105,j 0.15603) S( 1,3 )=( 0.04632,j 0.09985) S( 2,1 )=( -0.31144,j -0.13672) S( 2,3 )=( -0.23834,j 0.00774) S( 3,1 )=( -0.04449,j -0.13211) S( 3,2 )=( 0.24447,j -0.01175) S( 3,4 )=( -0.49996,j -0.03612) S( 4,3 )=( 0.49996,j 0.08895) 网损为 0.01737+j -0.05413,输电效率为 97.997 ---各支路功率情况4--- 平衡节点功率S= 0.36762+ i 0.25677 S( 1,2 )=( 0.32122,j 0.15676) S( 1,3 )=( 0.04639,j 0.10001) S( 2,1 )=( -0.31158,j -0.13731) S( 2,3 )=( -0.23842,j 0.00731) S( 3,1 )=( -0.04456,j -0.13225) S( 3,2 )=( 0.24456,j -0.01128) S( 3,4 )=( -0.50000,j -0.03647) S( 4,3 )=( 0.50000,j 0.08933) 网损为 0.01762+j -0.05323,输电效率为 97.970 3.思考题（见附录） 3.1潮流计算的方法有哪些？各有何特点？ 3.2 如果交给你一个任务，请你用已有的潮流计算软件计算北京城市电网的潮流，你应该做哪些工作？（收集哪些数据，如何整理，计算结果如何分析） 3.3设计中遇到的问题和解决的办法 三、课程设计总结或结论（见附录） 四、参考文献 1. 《电力系统计算：电子数字计算机的应用》，西安交通大学等合编。北京：水利电力出版社； 2. 《现代电力系统分析》，王锡凡主编，科学出版社； 3. 《电力系统稳态分析》，陈珩，中国电力出版社，1995年，第三版； 附录（手算过程，设计流程图、程序、表格、数据等） 因篇幅问题不能全部显示，请点此查看更多更全内容