电力系统潮流计算C语言程序及说明

1、实验目的根据所给的电力系统，编制潮流计算程序，通过计算机进行调试，最后完成一个切实可 行的电力系统计算应用程序。通过自己设计电力系统计算程序使同学对电力系统分析有进一 步理解，同时加强计算机实际应用能力的训练。程序计算原理1、概述应用计算机进行电力系统计算，首先要掌握电力系统相应计算的数学模型；其次是运用 合理的计算方法；第三则是选择合适的计算机语言编制计算程序。建立电力系统计算的相关数学模型，就是建立用于描述电力系统相应计算的有关参数间 的相互关系的数学方程式。该数学模型的建立往往要突出问题的主要方面，即考虑影响问题 的主要因素，而忽略一些次要因素，使数学模型既能正确地反映实际问题，又使计算

2、不过于 复杂。运用合理的计算方法，就是要求所选用的计算方法能快速准确地得出正确结果，同时还 应要求在解算过程中占用内存少，以利提高计算机的解题规模。选择合适的语言编写程序，就是首先确定用什么计算机语言来编制程序；其次是作出计 算的流程图；第三根据流程图用选择的语言编写计算程序。然后上机调试，直到语法上无错 误。本程序采用C语言进行编程。所编制的程序难免存在逻辑错误，因此先用一个已知结果的系统作为例题进行计算。用 程序计算的结果和已知结果相比较，如果结果相差甚远就要逐步分析程序的计算步骤，查出 问题的出处；如果结果比较接近，则逐步分析误差来源；直到结果正确为止。2、电力系统潮流计算的程序算法潮流

3、计算是电力系统分析中的一种最基本的计算，它的任务是对给定的运行条件确定系 统的运行状态，如母线上的电压（幅值及相角）、网络中的功率分布及功率损耗等。目前计算机潮流计算的方法主要有牛顿-拉夫逊算法和PQ分解法。牛顿-拉夫逊算法是 数学上求解非线形方程组的有效方法，具有较好的收敛性，曾经是潮流计算中应用比较普遍 的方法。PQ快速分解法是从牛顿-拉夫逊算法演变而来的，是将纯数学的牛顿-拉夫逊算法 与电力系统具体特点相结合并进行简化与改进而得出的PQ快速分解法比牛顿-拉夫逊算法 大大提高了计算速度和节省了内存，故而本程序以PQ快速分解法进行潮流计算。1）形成节点导纳矩阵（1）自导纳的形成对节点1其自导

4、纳Yii是节点i以外的所有节点都接地时节点i对地的总导纳。显然，Yii应等于与节点i相接的各支路导纳之和，即七=七.0+z七j式中，yi0为节点i与零电位节点之间的支路导纳；y为节点i与节点j之间的支路导纳。（2）互导纳的形成对节点i与节点k之间的互导纳是节点i、k之间的支路导纳的负值，即丫 =-） ik ik不难理解 =Yk。若节点i和k没有支路直接相连时，便有Yik=0（3）含变压器支路的处理若节点p、q间接有变压器，如下图所示，则可作出其n型等值电路为:z l:kP _QQ_ qkzP t! 1q质 I图1变压器n型等值电路则p、q的自导纳和节点间的互导纳分别为Y = + 三=-pp k

5、z kz z1 1 - k 1Yqq一 += kz k 2 z k 2 zpq=Yqp1kz2）计算不平衡功率 P.AQ并形成修正方程式对每一个PQ节点或每一个PV节点都可以根据下列公式计算出有功功率增量勘 P = P -P = P -VILv (G cos6 + B sin6 )(i = 1,2,，n-1)i is i is i j ij ij ij ijj=1而对于每一个PQ节点还可以根据下面的公式计算出无功功率增量8 Q = Q -Q = Q -V8v (G sin6 -B cos6 )(i = 1,2,.,m)i is i is i j ij ij ij ijj=1在有功功率增量和无功

6、功率增量不满足如下约束条件时max U P( k)J 8max，Q (k) 8iI利用PQ分解法则可以形成如下修正方程11121,-1P2P Q121222, n-1V Sn - 1,1n-1,2n - 1,n -1V S11121,m.BV21222, m.Bm,1m,2m,mVV3）利用因子表法求解修正方程在电网计算中经常遇到这样的问题，对方程组需要反复多次求解，而每次求解仅改变常 数项F,系数矩阵保持不变。按照一般的高斯消去法，对每一改变的常数项，形成包括常数 项及系数矩阵在内的增广矩阵，然后消去回代求出其解。可以看出，每次对增广矩阵中A 矩阵元素的消元都是重复的，为了避免这种重复，我们

7、把对相同的系数矩阵重复进行的消去 与对不同的常数项进行的消去分开进行，因此对系数矩阵的消去只需进行一次，并在消去的 过程中将对常数项进行消去运算的运算因子保存下来，形成所谓因子表，这就是因子表法。 因为因子表记录了高斯消去法对常数项进行消去的全部信息，利用它便可对不同常数项进行 消去，形成上三角矩阵，最后求出全部未知数。在使用PQ分解法时，其系数矩阵是在迭代过程中保持不变的，所以为了节省内存和缩 短运算时间我们采取了因子表法。同时由于电网的节点导纳矩阵矩阵是稀疏阵和对称阵，于 是我们可以采取只保存系数矩阵的上三角阵来使运算更为简化。若线性方程组一般形式如下：其中A =a11a12a22a13a

8、23a33a aaa1112131naaa22232naa333nannnaa2na3nannx 1f 1xf22x一f33:x x_ f _nn称为系数矩阵，X=*X3七P未知数向量，Tf2 f称为常数项向量。将矩阵A的元素进行如下处理:aj）二 aj-1）- a（k-i）a（幻（k -1,2, . . ., i -1; j = k +1, k + 2, . . ., n）a （i）ija（ i-1）ija（ i-1）ii（j =i +1,i + 2,n）得到因子表D11U12D22U13U23D33U1nU2nU3nDnn其中 D = 1/aa-1）Uj = a i（）i j ；）再利用因

9、子表进行前代过程，求出每次迭代后的常数项。其前代公式是:f（i-1）- Ufi-1）求得向量F = ffff （n） V； L 123 n -1再由因子表与前代得到的向量F，得到方程组1 U U U12131n1UU232n1U3n. !1工1X2X.3:=，j f（2）2琴Xnf（ n） n求解出此方程即可得到线性方程组的解向量X =气气七卜4）多次迭代最终求得V和5以及全线路功率利用上面所介绍的方法求解修正方程组0 PB B.B111121,n - 10PBBB黄:2122:2, n-1:2:BB.Bn-1,1n-1,2n - 1,n -10Pn 1iVn-1V 05 一V10512: 2

10、.V 05n - 1n - 1Q1V1QX2V2:QmVmB11B12B1,mB B B21222, m:Bm,1Bm,2- Bm,mV1V2:Vm=角 5 . 5 1。123 n-1可以求得V = 匕V2 V.W 1和i和max pQ（k）e己时，迭代结束。迭代结束后即可得再利用5（5 =5成）+口5（k），V（5 = V（Q +V（Q求得每次迭代后的结果。多次迭代当 iii ii其满足约束条件max 4 P（k） e和max 到各节点的V和5，再根据V、5来计算PV节点的无功功率Q和平衡节点的功率以及网络 中的功率分布。PV节点及平衡节点无功功率计算公式为：P = V Zv (G cos

11、5 + B sin 5 ) i i j ij ij ij ij j=1平衡节点有功功率计算公式为：Q = V ZV (G sin 5 - B cos 5 )i i j ij ij ij ij j=1以下图所标示的正方向，输电线路功率的计算公式如下：S = P + jQ = V I = V 2 y + V (V - V)ij ij ij i ij iij 0i i j ij5皿)图2支路功率计算对其进行实部虚部进行分解可得P、Q计算公式为：P = V2G - VV (G cos 5 + B sin 5 )ij i ij i j ij ij ij ijQ =-V2(B + B ) - VV (G

12、sin 5 + B cos 5 )ij i ijij 0i j ij ij ij ij程序及说明1、主要变量说明1）结构体类型说明（1）节点功率结构体struct Nodetypefloat P,Q;；其中，P为节点的有功功率，Q为无功功率。节点功率不区分负荷功率和发电机功率， 其值为本节点连接的各支路输入功率及节点所接负荷、发电机功率之和，且规定功率流入节 点为正，流出为负。详细说明参见下一章“算例及结果”的第二节“源数据格式说明”。（2）线路参数结构体struct Linetypefloat G,B,B0,k;其中，G、B为线路的导纳和容纳；B0为线路的考虑变压器n型等值电路后的对地充电

13、容纳的一半Bc/2； k为折算到标准变压器支路后的变压器变比。详细说明参见下一章“算例 及结果”的第二节“源数据格式说明”。2）变量说明表2程序主要变量说明主要变量类型含义Nodeint系统总节点数NPintPV+PQ节点数，即非平衡节点数NQintPQ节点数Numint *原始节点编号与程序表示编号映射数组Nostruct Nodetype *节点功率数组Vfloat *节点电压有效值数组Dltafloat *节点电压相角数组Ystruct Linetype *线路参数矩阵BP、BQfloat *有功、无功简化雅克比矩阵B、B”countunsigned intPQ迭代次数eP、eQcons

14、t有功、无功迭代精度控制kp、kqint有功、无功迭代结束标志dP、dQfloat *有功、无功不平衡量数组3、程序流程图开始图3程序主流程图开始Kp=0 且 Kq=0.Kp=0 ?Kp=0 且 Kq=0Kq=0?结束6 +=A6计算 Q、Q/NKq=0求+=Kp=0Count+计算 P,AP/NKp=1, Kq=1求逐图4迭代部分流程图4、程序代码/* FUNCTION : POWER FLOW*/* WRITTEN BY : HUANG&YANG&TONG */* LAST EDITED : 2008-11-24*/#include #include /* 宏定义 */#define eP

15、 0.00001#define eQ 0.00001#define Y_(i,j) (*(*(Y+i)-i+j)#define Yij (*(Yi+j)#define Yji (*(*(Y+j)-j+i)#define Pji Yji.G*cos(tmp)+Yji.B*sin(tmp)#define Pij Yij.G*cos(tmp)+Yij.B*sin(tmp)#define Qji Yji.G*sin(tmp)-Yji.B*cos(tmp)#define Qij Yij.G*sin(tmp)-Yij.B*cos(tmp)/*结构体变量定义*/struct Nodetype /* 节点功率

16、 */float P,Q;struct Linetype /* 线路类型 */float G,B,B0,k;/*子函数声明 */void in_node(); /* 读节点功率 */void in_line(); /* 读线路参数 */void B_Form(); /* 生成 BP、BQ 矩阵 */void factor(); /* 求因子表 */void solve(float *B,float *X,int N); /* 解方程组 */void PrtNode(); /*打印节点参数*/void ErrorMsg(int Flag); /* 错误提示信息 */* 全局变量声明 */int

17、Node; /* 节点数 */int *Num;/*保存原始节点序号*/ int NP,NQ=0; /*PV+PQ、PQ 节点数*/ struct Nodetype *No; /* 节点数据*/ struct Linetype *Y; /*线路参数*/float *BP,*BQ; /*有功、无功简化雅克比矩阵B*/ float *V;/*节点电压有效值*/float *Dlta;/*节点电压相角值*/unsigned int count=0; /* 迭代计数*/int i,j,k,m;/*通用下标值*/float tmp;/*临时数据暂存*/char *Type; /*节点类型*/FILE *

18、in,*out;/*输入、输出文件指针*/ /* 主函数 */ /* | * | * | * 主函数 * | * | * | */ int main(void)int kp=1,kq=1;/*P、Q精度判断：1-不满足，0-满足精度要求*/float *dP,*dPi,*dQ,*dQi;阻 P、 Q*/float Dltai;struct Linetype *Yi;struct Nodetype *Noi;float tP,tQ;if(in=fopen(Data.txt”,r)=NULL) ErrorMsg(1);if(out=fopen(out.txt”,w)=NULL) ErrorMsg(

19、2);in_node(); /*读取节点参数并统计PQ、PV节点数*/in_line(); /*读取线路参数并形成Y矩阵*/B_Form();/*形成 B(BP&BQ)矩阵*/factor();/*求B因子式(仍保存于BP&BQ)*/for(i=0;iNQ;i+) *(V+i)=1; /*对 PQ 节点电压 V 赋初值*/ dP=(float *)malloc(sizeof(float)*NP); /*AP*/dQ=dP; /* Q*/*A P、 Q不同时存在，故而可共用空间*/loop:/* 开始迭代 */if(kp=0&kq=0) goto loopEnd;退出*/count+; /*迭代

20、次数加一*/if(count=65535) ErrorMsg(99); /*不收敛， kp=0;for(i=0;iNP;i+)dPi=dP+i;Yi=*(Y+i)-i;Dltai=*(Dlta+i);*dPi=0;for(j=0;jj) *dPi+=*(V+j)*(Pji);else *dPi+=*(V+j)*(Pij);/*注意到Y矩阵为上三角矩阵，ij时要交换下标*/*dPi*=*(V+i);*dPi=(*(No+i).P-*dPi; /*求得 Pi*/if(fabs(*dPi)0 x8fffffff) ErrorMsg(99); /*不收敛，退出*/if(fabs(*dPi)eP) kp

21、=1; /*有不满足精度的 P 即令 kp=1*/ *dPi/=*(V+i); /* 求得常数项 Pi/Vi*/if(kp=0) goto loopQ;solve(BP,dP,NP);for(i=0;iNP;i+) *(Dlta+i)+=(*(dP+i)/(*(V+i); /*修正相角0 += 0 */loopQ:if(kp=0&kq=0) goto loopEnd;kq=0;for(i=0;iNQ;i+)dQi=dQ+i;Yi=*(Y+i)-i;Dltai=*(Dlta+i);*dQi=0;for(j=0;jj) *dQi+=*(V+j)*(Qji);else *dQi+=*(V+j)*(Q

22、ij);/*注意到Y矩阵为上三角矩阵，ij时要交换下标*/*dQi*=*(V+i);*dQi=(*(No+i).Q-*dQi; /*求得 Qi*/if(fabs(*dQi)0 x8fffffff) ErrorMsg(99); /*不收敛，退出*/if(fabs(*dQi)eQ) kq=1; /*有不满足精度的 Q 即令 kq=1*/ *dQi/=*(V+i); /* 求得常数项 Qi/Vi*/if(kq=0) goto loop;solve(BQ,dQ,NQ);for(i=0;iNQ;i+) *(V+i)+=*(dQ+i); /*修正 PQ 节点电压 V+=A V*/goto loop; /*

23、无功迭代，则必定需要下一轮回迭代判断*/loopEnd:/* 迭代结束 */free(dP); /*释放内存空间*/*计算PV节点和平衡节点的无功功率Q*/for(i=NQ;iNode;i+)Noi=No+i;Yi=*(Y+i)-i;Dltai=*(Dlta+i);for(j=0;jj) (*Noi).Q+=*(V+j)*(Qji);else (*Noi).Q+=*(V+j)*(Qij);/*注意到Y矩阵为上三角矩阵，ij时要交换下标*/(*Noi).Q*=*(V+i);/*计算平衡节点的有功功率P*/i=NP;Noi=No+i;Dltai=*(Dlta+i);for(j=0;jj时要交换下标

24、*/(*Noi).P*=*(V+i);/*输出最终结果*/fprintf(out,nn【潮流计算结果(节点)】(迭代次数k=%3d )n,count-1); PrtNode();/*计算全部线路功率*/fprintf(out,nn【潮流计算结果(线路)】n);fprintf(out,” 线路PQn);for(k=0;kNode;k+)i=*(Num+k);Yi=*(Y+i)-i;Dltai=*(Dlta+i);Noi=No+i;for(m=0;mNode;m+)j=*(Num+m);if(j=i) continue;tmp=Dltai-*(Dlta+j); /*tmp 即。ij*/if(ji)

25、if(Yji.B=0) continue; /*若 Bij=0,则节点 i、j 无直接联系*/ tP=*(V+j)*(Pji);tP=*(V+i)*Yji.G-tP;tP*=*(V+i);tQ=-*(V+j)*(Qji);tQ-=*(V+i)*(Yji.B-Yji.B0/Yji.k);tQ*=*(V+i);elseif(Yij.B=0) continue; /喏 Bij=0,则节点 i、j 无直接联系*/ tP=*(V+j)*(Pij);tP=*(V+i)*Yij.G-tP;tP*=*(V+i);tQ=-*(V+j)*(Qij);tQ-=*(V+i)*(Yij.B-Yij.B0);tQ*=*(

26、V+i);fprintf(out,”S%d,%d= (%10.6f,%10.6f)n”,k+1,m+1,-tP,-tQ);fclose(out);system(cmd /c start out.txt);return(0);/* f * f * f * 主函数 * f * f * f */* 主函数 */ /* 子函数：读节点数据 */ void in_node()struct Nodetype *Noi; /* 临时变量*/fscanf(in,%d %d”,&Node,&k);/*读取节点数 Node*/NP=Node-1; /*PV+PQ节点数，即非平衡节点数目*/Num=(int *)m

27、alloc(sizeof(int)*Node); /*开 Node 个空间，每节点一个*/V=(float *)malloc(sizeof(float)*Node); /* 电压 */Dlta=(float *)malloc(sizeof(float)*Node); /* 电压相角*/No=(struct Nodetype *)malloc(sizeof(struct Nodetype)*Node);/*节点功率 */ j=1;while(k!=0) /*若k=0,表明节点数据读取完毕*/switch(k)case 1:k=NQ;NQ+;break; /*NQ 统计 PQ 节点个数*/case

28、 2:k=NP-j;j+;break; /*从 NP-1 个空间倒着保存 PV 节点*/case 3:k=NP;break; /* 平衡节点*/default:ErrorMsg(3);Noi=No+k;fscanf(in,%d %f %f %f %f”,&i,&(*Noi).P,&(*Noi).Q,V+k,Dlta+k);i-; /*节点编号减一，以和程序表达方式兼容*/*(Num+i)=k; /*第i个Num元素中存放i节点在No中的下标*/ fscanf(in,”%d”,&k); /* 读取节点类型 */if(NQ+j!=Node) ErrorMsg(4); /*检验节点数据个数是否够No

29、de个*/fprintf(out,”【节点参数表】n);PrtNode();fprintf(out,” 总节点：%dn PQ 节点：%dn PV 节点：dn”,Node,NQ,NP-NQ);/* 子函数：读线路数据，并形成节点导纳矩阵Y */void in_line()struct Linetype *Yi;float R,X,k,B;m=sizeof(struct Linetype);Y=(struct Linetype *)malloc(m*Node); /*先开 Node 行，每一个节点一行*/for(i=0;iNode;i+)/*再在第 i 行上面开辟 Node-i 列*/*即以上三角

30、存储Y矩阵*/*(Y+i)=(struct Linetype *)malloc(m*(Node-i);Yi=*(Y+i)-i;for(j=i;jj) tmp=i;i=j;j=tmp;k=1/k;B*=k;Yi=*(Y+i)-i; /*以Yi代替*(Y+i)-i,简化表达式并避免重复计算*/Yij.B0=B; /*保存ij0、ji0对地充电电容到Bij0*/Yij.k=k; /*且有 B0ji=B0ij/k*/tmp=R*R+X*X;R/=tmp;X/=tmp;Yij.G=-R; /*生成互导纳*/Yij.B=X;(*(*(Y+i).G+=R; /*将线路互导纳累加到自导纳*/(*(*(Y+i)

31、.B+=-X;(*(*(Y+j).G+=R;(*(*(Y+j).B+=-X;fclose(in);fprintf(out,n【节点导纳矩阵Y】n);for(k=0;kNode;k+)i=k;i=*(Num+i); /*查取第i节点在程序中存储序号*/for(j=0;jk;j+) fprintf(out,ttt);for(m=k;mNode;m+)j=*(Num+m); /*查取第j节点在程序中存储序号*/if(ij) fprintf(out,”(%10.6f,%10.6f) ”，丫_(i,j).G,Y_(i,j).B);else fprintf(out,”(%10.6f,%10.6f) ”,Y

32、_(j,i).GY_(j,i).B);fprintf(out,n);3 “/*子函数：生成 BP、BQ 矩阵 */void B_Form()float *BPi,*BQi;struct Linetype *Yi;int size=sizeof(float);BP=(float *)malloc(size*NP); /* 以上三角存储 */for(i=0;iNP;i+) *(BP+i)=(float *)malloc(size*(NP-i);for(i=0;iNP;i+)BPi=*(BP+i)-i; /* 以 BPi 代替*(BP+i)-i,避免重复计算 */ Yi=*(Y+i)-i;for(j

33、=i;jNP;j+) *(BPi+j)=Yij.B; /*(BPi+j)即相当于 BPij*/ BQ=BP; /* BP包含BQ, BP左上角的NQ*NQ子阵即BQ */* 子函数：求因子表 */ void factor()float *BPi,*BPk,*BQi;for(i=0;iNP;i+)BPi=*(BP+i)-i;for(k=0;ki;k+)BPk=*(BP+k)-k;tmp=(*(BPk+i)/(*(BPk+k);for(j=i;jNP;j+) (*(BPi+j)-=tmp*(*(BPk+j);*(BPi+i)=1/(*(BPi+i);for(j=i+1;jNP;j+) *(BPi+

34、j)*=*(BPi+i);/* 子函数：解方程组 */ void solve(float *B,float *X,int N)float *Bi,*Xi;for(i=0;iN;i+) *(X+i)=-*(X+i);/*对常数项取负 */*对常数列进行前代*/for(i=0;iN;i+)Bi=*(B+i)-i;Xi=X+i;for(j=i+1;j=0;i-)Bi=*(B+i)-i;Xi=X+i;for(j=N-1;ji;j-) *Xi-=*(Bi+j)*(X+j);/*打印节点参数*/void PrtNode()struct Nodetype *Noi;fprintf(out,”节点 类型 PQ

35、V(5 n);for(i=0;iNode;)j=*(Num+i); /*查取第i节点在程序中存储序号*/Noi=No+j;if(jNQ) Type=PQ”; else Type=PV”;if(j=NP) Type=BS”;fprintf(out,%3d %s %10.6f %10.6f %10.6f %10.6fn,+i,Type,(*Noi).P,(*Noi).Q,*(V+j),*(Dlta+j)/0.017453);/* 子函数：错误信息 */void ErrorMsg(int Flag)switch(Flag)case 1:printf(ntError(1): Failed to Ope

36、n File Data.txt!”);break;case 2:printf(ntError(2): Failed to Creat File out.txt!”);break;case 3:printf(ntError(3): Node Data Error, Please Check!);break;case 4:printf(ntError(4): Lack Node Data, Please Check!);break;case 99:printf(k=%dntError(99): Its Emanative!”,count);break;getch();fclose(out);exi

37、t(Flag);算例及结果1、原始数据程序编写过程中，采用了电力系统分析教材的4节点、IEEE9节点参数进行了验证。本说明书中仅列出4节点验算过程。4节点测试数据来源于参考文献中的电力系统分析（下册）第61页例11-5，用于检 验程序正确性的参考计算结果见67页。对其基本情况总结如下。1）网络连接图图6算例电力系统网络图2）节点参数表4算例节点参数表节点类型节点编号有功功率P无功功率Q电压V电压相角511-0.30-0.180012-0.55-0.1300230.501.10034001.0503）线路参数表5算例线路参数表节点1节点2阻抗R感抗X变压器变比k充电容纳Bc/2120.10.41

38、0.01528140.120.510.01920240.080.410.014133100.30.90909104）数据文件data.txt34001.050230.501.10011-0.30-0.180012-0.55-0.13000120.10.410.015283100.30.9090910140.120.510.01920240.080.410.01413图7数据文件格式2、源数据格式说明源数据有功功率、无功功率、电压、阻抗、感抗、对地充电容纳均以标幺值表示。数据文件必须命名为data.txt且与潮流计算程序放置于同一个文件中。数据文件data.txt 包含两类参数：节点参数和线路参

39、数。节点数据块与线路数据块之间用数字0作为间隔，即 在节点数据块结束后，另起一行输入0，然后再在后面按格式要求录入线路参数1）节点参数节点参数包括：系统节点数Node、节点功率（有功P、无功Q），节点电压（有效值V、 相角Delta）参数组织格式：（1）节点数Node节点数Node写在参数文件的开头，如：4表明为四节点系统。（2）功率和电压P/Q/V/Delta首先给出节点参数示例：230.501.100第一列数字2表明该节点的类型为2-PV节点；第二列数字3表明该行数据为节点3的 参数；后面三列依次为P、Q、V的给定值，给定值为0，表明该项参数未知；第六列为相 角0，非平衡节点的0即为PQ迭

40、代的初始相角值，平衡节点的即为给定的相角值。节点功率为各支路输入功率之和，且规定功率流入节点为正，流出为负。故而负荷功率 为负值，发电机功率为正值。（3）节点类型的判断根据节点的给定参数可以将节点分为三种类型：1）PQ节点：给定P、Q初始值的节点，用数字代码1表示；2）PV节点：给定P、V初始值的节点，用数字代码2表示；3）平衡节点：没有给定P、Q初始值，仅给定V的初始值，用数字代码3表示；节点本身不分先后，不必刻意按编号大小来列写各节点，只需将每一个节点的参数按格 式列写即可。2）线路参数线路参数包括：线路阻抗（R,X）.变压器变比k和对地充电（接地支路）容纳b。每一条线 路包括节点号在内共

41、有6个参数，6个参数缺一不可。（1）线路阻抗Z首先给出线路参数示例：1 4 0.12 0.510.01920前两列代表线路两端的节点编号，亦即线路编号14（即41）；第三四列代表线路阻抗 Z=0.12+j0.5 ；第五列为变压器变比k；最后一列为线路的对地充电电容的一半，即B/2。k不为0。当k=1则表明该支路为普通支路；否则该支路为变压器支路。普通支路没有 任何特殊要求，但对于变压器支路，有以下注意事项。（2）变压器支路首先给出变压器线路示例：3 1 0 0.3 0.909091 01）变比k为转换为变压器支路的标准等值电路（如下图）后的变比。图8变压器n型等值电路2）变压器线路的编号31有特定含义：3对应节点p，1对应节点q。即变压器支路的编 号对节点顺序有要求，p节点编号在前，q节点在后。3）q节点为全系统参考电压侧。程序智能化读取线路数据，用户无需刻意对线路参数进行归类、排序，可随意输入，