多平衡节点+普通PQ分解

1、/多平衡节点+普通PQ分解#include "fstream.h"#include "math.h"#include "stdlib.h"#include "iomanip.h"#include "time.h"#define LINEMAX 5000/最大线路数#define GENERATORMAX 500/最大发电机数#define LOADMAX 2000/最大负荷数#define NODEMAX 2000/最大节点数#define SWINGMAX 20/最大平衡节点数#define

2、PVMAX 500/最大PV节点数#define NODEFACTOR10/导纳矩阵中非零非对角元素的个数相对于/最大节点数(NODEMAX)的倍数.#define Deg_to_Rad 0.017453292/度到弧度的转换系数#define Rad_to_Deg 57.29577951/弧度到度的转换系数#define SinglePai 3.14159265/圆周率#define DoublePai 6.2831853/两倍的圆周率struct Line/线路参数结构定义int Node_No2;/线路两端节点名(号):0-左节点;1-右节点int Flag;/线路类型标志:0-普通支路

3、;1、2-变压器支路。其中,/1-非标准变比在左侧节点,2-非标准变比在右侧节点。/对接地支路，左右节点名（号）相同。double RXBK3;/0-R;1-X;2-Bc/2 或 KLLineLINEMAX;struct Generator/发电机参数结构定义int Node_No;/发电机节点名(号)int Flag;/发电机节点类型标志:0-平衡节点;1-PQ节点;2-PV节点double PQV2;/对平衡节点，0-V,1-Angle;对PQ节点,0-P,1-Q;对PV节/点,0-P,1-V。GGenGENERATORMAX;struct Load/负荷参数结构定义int Node_No

4、;/负荷节点名(号)int Flag;/负荷节点静特性标志:0-不计静特性;1-计静特性。double ABC6;/PL=a1*V*V+b1*V+c1,QL=a2*V*V+b2*V+c2。/0-a1;1-a2;2-b1;3-b2;4-c1;5-c2。LLoadLOADMAX;int Node_Name_NewtoOldNODEMAX;/新节点名（号）->旧节点名（号）int Node_FlagNODEMAX;/节点类型标志：0-平衡,1-PQ,2-PVint Line_NodeNameLINEMAX2;/线路的左、右节点新名（号）int Line_No_NewtoOldLINEMAX;/

5、新线路号->旧线路号int Line_FlagLINEMAX;/新线路的类型标志:0,1,2说明同Line结构int Gen_NodeNameGENERATORMAX;/发电机节点的新节点名（号）int Gen_No_NewtoOldGENERATORMAX;/新发电机顺序号->旧发电机顺序号int Gen_SWNodeSWINGMAX2;/平衡节点数据：0-新节点名（号）;/1-对应的旧发电机顺序号int Gen_PVNodePVMAX2;/发电机PV节点数据：0-新节点名（号）;/1-对应的旧发电机顺序号int Gen_PQNodeGENERATORMAX2;/发电机PQ节点数

6、据：0-新节点名（号）;/1-对应的旧发电机顺序号int Load_NodeNameLOADMAX;/负荷节点的新节点名（号）int Load_No_NewtoOldLOADMAX;/新负荷顺序号->旧负荷顺序号/读数据子程void Data_Input(int &Num_Line,int &Num_Gen,int &Num_Load,double &Eps,int &Iter_Max,int &VolIni_Flag,int &VolRes_Flag)int i,j;/ 从键盘读数据文件名char argv20;cout<&

7、lt;"Please input diskette data-file name:"<<endl; cin.getline(argv,sizeof(argv);/ 打开数据文件 ifstream infile(argv);if(infile.fail()cerr<<"Error opening the diskette data-file:"<<argv<<endl;exit(0);else/读记数器变量值infile>>Num_Line>>Num_Gen>>Num_Lo

8、ad>>Eps>>Iter_Max>>VolIni_Flag>>VolRes_Flag;/*cout<<setw(5)<<Num_Line<<setw(5)<<Num_Gen<<setw(5)<<Num_Load<<setw(10)<<Eps<<Iter_Max<<VolIni_Flag<<VolRes_Flag<<endl;*/if(Num_Line>LINEMAX)cout<<&qu

9、ot;Lines Number > LINEMAX!"<<endl;exit(0);if(Num_Gen>GENERATORMAX)cout<<"Generators Number > GENERATORMAX!"<<endl;exit(1);if(Num_Load>LOADMAX)cout<<"Loads Number > LOADMAX!"<<endl;exit(2);/读线路参数for(i=0;i<Num_Line;i+)infile>&

10、gt;LLinei.Node_No0>>LLinei.Node_No1>>LLinei.Flag;for(j=0;j<3;j+)infile>>LLinei.RXBKj;/*cout<<setw(5)<<LLinei.Node_No0<<setw(5)<<LLinei.Node_No1<<setw(5)<<LLinei.Flag;cout<<setw(10)<<LLinei.RXBK0<<setw(10)<<LLinei.RXBK1&

11、lt;<setw(10)<<LLinei.RXBK2<<endl;*/读发电机参数for(i=0;i<Num_Gen;i+)infile>>GGeni.Node_No>>GGeni.Flag>>GGeni.PQV0>>GGeni.PQV1;/*cout<<setw(5)<<GGeni.Node_No<<setw(5)<<GGeni.Flag<<setw(10)<<GGeni.PQV0<<setw(10)<<GGeni

12、.PQV1<<endl;*/读负荷参数for(i=0;i<Num_Load;i+)infile>>LLoadi.Node_No>>LLoadi.Flag;for(j=0;j<6;j+)infile>>LLoadi.ABCj;/*cout<<setw(5)<<LLoadi.Node_No<<setw(5)<<LLoadi.Flag;for(j=0;j<6;j+)cout<<setw(10)<<LLoadi.ABCj;cout<<endl;*/inf

13、ile.close();/序号处理子程void Node_Sequen(int &Num_Node,int Num_Line,int Num_Gen,int Num_Load, int &Num_Swing,int &Num_GPV,int &Num_GPQ)int i,j,Flag,temp,np;int Node_NameNODEMAX2;/0-节点名(号);1-节点出线数 /统计各节点的出线数for (i=0;i<NODEMAX;i+)Node_Namei1=0; /节点出线数初始化为0for(i=0;i<Num_Line;i+)if(LLin

14、ei.Node_No0=LLinei.Node_No1)continue;/接地支路(左右节点相同)不在出线统计之内Flag=0;/左节点出线数分析开始for(j=0;j<Num_Node;j+)if(LLinei.Node_No0=Node_Namej0)/该节点已经在节点/数组中出现,只需Node_Namej1+;/出线数加1。Flag=1;if(Flag=1)break;if(Flag=0)/该节点还没在/节点数组中出Node_NameNum_Node0=LLinei.Node_No0;/现,需将该节Node_NameNum_Node1+;/点名(号)添加Num_Node+;/到节

15、点数组中,if(Num_Node>NODEMAX)/然后该节点的/出线数加1,并cout<<"Nodes Number > NODEMAX!"<<endl;/将节点数也加exit(4);/1。Flag=0;/右节点出线数分析开始for(j=0;j<Num_Node;j+)if(LLinei.Node_No1=Node_Namej0)/该节点已经在节点/数组中出现,只需Node_Namej1+;/出线数加1Flag=1;if(Flag=1)break;if(Flag=0)/该节点还没在/节点数组中出Node_NameNum_Node0

16、=LLinei.Node_No1;/现,需将该节Node_NameNum_Node1+;/点名(号)添加Num_Node+;/到节点数组中,if(Num_Node>NODEMAX)/然后该节点的/出线数加1,并cout<<"Nodes Number > NODEMAX!"<<endl;/将节点数也加exit(5);/1。/节点出线数统计完毕，屏幕输出。其中，Num_Node为总节点数。/*cout<<endl<<"Node: Number, Name and out-line number"&l

17、t;<endl;for(i=0;i<Num_Node;i+)cout<<setw(5)<<i<<setw(5)<<Node_Namei0<<setw(5)<<Node_Namei1<<endl;*/根据出线数由小到大的顺序对节点进行排序(冒泡算法)for(i=0;i<Num_Node-1;i+)np=i;for(int j=i+1;j<Num_Node;j+)if(Node_Namenp1>Node_Namej1)np=j;temp=Node_Namei0;Node_Namei0=

18、Node_Namenp0;Node_Namenp0=temp;temp=Node_Namei1;Node_Namei1=Node_Namenp1;Node_Namenp1=temp;/*cout<<endl<<"Node Sequence: Number, Name and out-line number"<<endl;for(i=0;i<Num_Node;i+)cout<<setw(5)<<i<<setw(5)<<Node_Namei0<<setw(5)<<N

19、ode_Namei1<<endl;*/平衡节点统计:总数及各节点的名(号)int Node_Name_SwingSWINGMAX;for(i=0;i<Num_Gen;i+)if(GGeni.Flag=0)Node_Name_SwingNum_Swing=GGeni.Node_No;Num_Swing+;if(Num_Swing>SWINGMAX)cout<<"Swinging Generators Number > SWINGMAX!"<<endl;exit(6);/*cout<<endl<<&q

20、uot;Swing: Number and Name "<<Num_Swing<<endl;for(i=0;i<Num_Swing;i+)cout<<setw(5)<<i<<setw(5)<<Node_Name_Swingi<<endl;*/根据出线数由小到大的顺序对节点排序,并将平衡节点排在最后(序号最大)int Nswing=0,Nnode=0;for(i=0;i<Num_Node;i+)Flag=0;for(j=0;j<Num_Swing;j+)if(Node_Namei0=N

21、ode_Name_Swingj)Flag=1;if(Flag=1)break;/Flag=1时,表示该节点为平衡节点,/需排在靠后的位置上。if(Flag=0)Node_Name_NewtoOldNnode=Node_Namei0;Nnode+;else/最后的各平衡节点间也仍然按出线数由小到大的顺序排列Node_Name_NewtoOldNum_Node-Num_Swing+Nswing=Node_Namei0;Nswing+;/*cout<<endl<<"Node Sequence with swing nodes being in the end: &q

22、uot;<<"Number, Name"<<endl;for(i=0;i<Num_Node;i+)cout<<setw(5)<<i<<setw(5)<<Node_Name_NewtoOldi<<endl;*/新线路类型标志赋初值for(i=0;i<Num_Line;i+)Line_Flagi=LLinei.Flag;/线路名(号)处理:变成新的节点名(号)且左节点的绝对值小于右节点的绝对值for(i=0;i<Num_Line;i+)Flag=0;for(j=0;j<N

23、um_Node;j+)if(LLinei.Node_No0=Node_Name_NewtoOldj)/左节点处理Line_NodeNamei0=j;/赋新名(号)Flag=1;if(Flag=1)break;Flag=0;for(j=0;j<Num_Node;j+)if(LLinei.Node_No1=Node_Name_NewtoOldj)/右节点处理Line_NodeNamei1=j;/赋新名(号)Flag=1;if(Flag=1)break;if(Line_NodeNamei0>Line_NodeNamei1)/左节点的绝对值小于 /右节点的绝对值处理if(LLinei.Fl

24、ag=1)Line_Flagi=2;/变压器的非标准if(LLinei.Flag=2)Line_Flagi=1;/变比侧发生变化temp=Line_NodeNamei0;Line_NodeNamei0=Line_NodeNamei1;Line_NodeNamei1=temp;/*cout<<endl<<"Line:line number,type,left node,right node"<<endl;for(i=0;i<Num_Line;i+)cout<<setw(5)<<i<<setw(5)&

25、lt;<Line_Flagi<<setw(5)<<Line_NodeNamei0<<setw(5)<<Line_NodeNamei1<<endl;*/线路排序:按照左节点的绝对值由小到大、若左节点的绝对值相等则按照右节/点的绝对值由小到大顺序排序（双排序冒泡算法）for(i=0;i<Num_Line;i+)Line_No_NewtoOldi=i;for(i=0;i<Num_Line-1;i+)np=i;for(j=i+1;j<Num_Line;j+)if(Line_NodeNamej0<Line_Node

26、Namenp0 |(Line_NodeNamej0=Line_NodeNamenp0 &&Line_NodeNamej1<Line_NodeNamenp1)np=j;temp=Line_NodeNamenp0;Line_NodeNamenp0=Line_NodeNamei0;Line_NodeNamei0=temp;temp=Line_NodeNamenp1;Line_NodeNamenp1=Line_NodeNamei1;Line_NodeNamei1=temp;temp=Line_No_NewtoOldnp;Line_No_NewtoOldnp=Line_No_New

27、toOldi;Line_No_NewtoOldi=temp;temp=Line_Flagnp;Line_Flagnp=Line_Flagi;Line_Flagi=temp;/*cout<<endl<<"Line sequencing: new line number,new type, new left"<<" node, new right node, old line number"<<endl;for(i=0;i<Num_Line;i+)cout<<setw(5)<<i&

28、lt;<setw(5)<<Line_Flagi<<setw(5)<<Line_NodeNamei0<<setw(5)<<Line_NodeNamei1<<setw(5)<<Line_No_NewtoOldi<<endl;*/发电机节点名(号)处理:变成新的节点名(号)for(i=0;i<Num_Gen;i+)Flag=0;for(j=0;j<Num_Node;j+)if(GGeni.Node_No=Node_Name_NewtoOldj)Gen_NodeNamei=j;/赋新名(号

29、)Flag=1;if(Flag=1)break;/*cout<<endl<<"Generator node new name:number, new node name"<<endl;for(i=0;i<Num_Gen;i+)cout<<setw(5)<<i<<setw(5)<<Gen_NodeNamei<<endl;*/发电机排序：按照新节点名（号）由小到大的顺序排序，并找出新发电机序号/对应的旧发电机序号for(i=0;i<Num_Gen;i+)Gen_No_Ne

30、wtoOldi=i;for(i=0;i<Num_Gen-1;i+)np=i;for(j=i+1;j<Num_Gen;j+)if(Gen_NodeNamej<Gen_NodeNamenp)np=j;temp=Gen_NodeNamenp;Gen_NodeNamenp=Gen_NodeNamei;Gen_NodeNamei=temp;temp=Gen_No_NewtoOldnp;Gen_No_NewtoOldnp=Gen_No_NewtoOldi;Gen_No_NewtoOldi=temp;/*cout<<endl<<"Generator seq

31、uencing: new gen number, new node, old"<<" gen number"<<endl;for(i=0;i<Num_Gen;i+)cout<<setw(5)<<i<<setw(5)<<Gen_NodeNamei<<setw(5)<<Gen_No_NewtoOldi<<endl;*/负荷节点名(号)处理:变成新的节点名(号)for(i=0;i<Num_Load;i+)Flag=0;for(j=0;j<Num_

32、Node;j+)if(LLoadi.Node_No=Node_Name_NewtoOldj)Load_NodeNamei=j;/赋新名(号)Flag=1;if(Flag=1)break;/*cout<<endl<<"Load node new name: number, new node name"<<endl;for(i=0;i<Num_Load;i+)cout<<setw(5)<<i<<setw(5)<<Load_NodeNamei<<endl;*/负荷排序：按照新节点

33、名（号）由小到大的顺序排序，并找出新负荷序号/对应的旧负荷序号for(i=0;i<Num_Load;i+)Load_No_NewtoOldi=i;for(i=0;i<Num_Load-1;i+)np=i;for(j=i+1;j<Num_Load;j+)if(Load_NodeNamej<Load_NodeNamenp)np=j;temp=Load_NodeNamenp;Load_NodeNamenp=Load_NodeNamei;Load_NodeNamei=temp;temp=Load_No_NewtoOldnp;Load_No_NewtoOldnp=Load_No_

34、NewtoOldi;Load_No_NewtoOldi=temp;/*cout<<endl<<"Load sequencing: new Load number, new node name, "<<"old Load number"<<endl;for(i=0;i<Num_Load;i+)cout<<setw(5)<<i<<setw(5)<<Load_NodeNamei<<setw(5)<<Load_No_NewtoOldi&l

35、t;<endl;*/从发电机节点数据中归纳出平衡节点、PV节点、PQ节点的新节点名（号）和对/应的旧发电机序号，并对平衡节点和PV节点修正其节点类型标志for(i=0;i<Num_Node;i+)Node_Flagi=1;/节点类型赋初值1(PQ节点)Nswing=0;for(i=0;i<Num_Gen;i+)j=Gen_No_NewtoOldi;/发电机节点旧顺序号if(GGenj.Flag=0)Gen_SWNodeNswing0=Gen_NodeNamei;/发电机节点名称Gen_SWNodeNswing1=j;Node_FlagGen_NodeNamei=0;Nswin

36、g+;else if(GGenj.Flag=1)Gen_PQNodeNum_GPQ0=Gen_NodeNamei;/发电机节点名称Gen_PQNodeNum_GPQ1=j;(Num_GPQ)+;else if(GGenj.Flag=2)Gen_PVNodeNum_GPV0=Gen_NodeNamei;/发电机节点名称Gen_PVNodeNum_GPV1=j;Node_FlagGen_NodeNamei=2;(Num_GPV)+;if(Num_GPV>PVMAX)cout<<"PV Generators Number > PVMAX!"<<

37、endl;exit(7);double Y_DiagNODEMAX2;/节点导纳阵的对角元:0-实部;/1-虚部。double Y_UpTriNODEMAX*NODEFACTOR2;/节点导纳阵上三角的非零元:/0-实部;1-虚部。int Foot_Y_UpTriNODEMAX*NODEFACTOR;/上三角按行压缩存储的非零元的/列足码。int Num_Y_UpTriNODEMAX;/上三角各行非零元素的个数int No_First_Y_UpTriNODEMAX;/上三角各行第一个非零元素在/Y_UpTri中的顺序号。int Foot_Y_DownTriNODEMAX*NODEFACTOR;

38、/下三角按行压缩存储的非零元的/列足码。int Num_Y_DownTriNODEMAX;/下三角各行非零元素的个数int No_First_Y_DownTriNODEMAX;/下三角各行第一个非零元素在按/行压缩存储序列中的顺序号int No_Y_DownTri_RowtoColNODEMAX*NODEFACTOR;/下三角某行非零元所对/应的按列压缩存储序列/中的序号/形成节点导纳矩阵1(不包括线路充电容纳及非标准变比的影响)子程void Y_Bus1(int Num_Node,int Num_Line,/*int Num_Load,*/int Num_Swing)int i,j,k,k_

39、old,Flag,l;double X,B;/线路参数工作单元l=0;for(i=0;i<Num_Node-Num_Swing;i+)/初始化Y_Diagi1=0.0;Num_Y_UpTrii=0;for(k=0;k<Num_Line;k+)i=Line_NodeNamek0;/线路左节点j=Line_NodeNamek1;/线路右节点if(i>=Num_Node-Num_Swing)/左右节点均为平衡节点，对导纳阵无break;/影响。k_old=Line_No_NewtoOldk;/对应的旧线路顺序号X=LLinek_old.RXBK1;/取线路电抗值B=-1.0/X;/

40、不计线路电阻后的线路支路电纳if(j>=Num_Node-Num_Swing)/左为普通节点，右为平衡节点Y_Diagi1=Y_Diagi1+B;else/左、右节点均为普通节点Flag=0;if(k>0&&(i=Line_NodeNamek-10)&&(j=Line_NodeNamek-11)Flag=1;/多回线Y_Diagi1=Y_Diagi1+B;if(i!=j)/非接地支路Y_Diagj1=Y_Diagj1+B;if(Flag=0)/第一回线Y_UpTril1=-B;Foot_Y_UpTril=j;Num_Y_UpTrii+;l+;if(l

41、>NODEMAX*NODEFACTOR)cout<<"Number of none-zero elements of ""up_triangle > NODEMAX*NODEFACTOR!"<<endl;exit(8);else/多回线Y_UpTril-11=Y_UpTril-11-B;No_First_Y_UpTri0=0;for(i=0;i<Num_Node-Num_Swing;i+)No_First_Y_UpTrii+1=No_First_Y_UpTrii+Num_Y_UpTrii;/稀疏导纳矩阵上三角及对

42、角元的结果输出/*cout<<endl<<"对角元素的新节点名(号)、旧节点名(号)及元素值"<<endl;for(i=0;i<Num_Node-Num_Swing;i+)cout<<setw(5)<<i<<setw(5)<<Node_Name_NewtoOldi<<setw(10)<<Y_Diagi0<<setw(10)<<Y_Diagi1<<endl;cout<<endl<<"上三角 &q

43、uot;<<l<<" 个非对角元素的顺序号、列足码及元素值"<<endl;for(k=0;k<l;k+)cout<<setw(5)<<k<<setw(5)<<Foot_Y_UpTrik<<setw(10)<<Y_UpTrik0<<setw(10)<<Y_UpTrik1<<endl;cout<<endl<<"导纳阵上三角每行非对角元素的个数：行号及个数"<<endl;for

44、(i=0;i<Num_Node-Num_Swing;i+)cout<<setw(5)<<i<<setw(5)<<Num_Y_UpTrii<<endl;cout<<endl<<"上三角每行第一个非对角元素的顺序号：行号及顺序号"<<endl;for(i=0;i<Num_Node-Num_Swing+1;i+)cout<<setw(5)<<i<<setw(5)<<No_First_Y_UpTrii<<endl;

45、*/导纳矩阵下三角按行压缩存储时各行非零元的个数、每行第一个非零元的序/号、按行压缩存储时非零元的列足码、某一非零元所对应的下三角阵按列压/缩存储时相同非零元的序号，上述这些值的求取目的在于快速地处理下面迭/代过程中修正方程Jacobi矩阵的按行压缩形式下的取值、消去和回代运算。int Row_DownNODEMAX;/下三角某列非零元序号的下限工作单元int Row_UpNODEMAX;/下三角某列非零元序号的上限工作单元int li;for(i=0;i<Num_Node-Num_Swing;i+)/下三角各行非零元个数数组清零Num_Y_DownTrii=0;for(j=0;j<

46、;Num_Node-Num_Swing;j+)/该循环统计下三角各行非零元个数for(k=No_First_Y_UpTrij;k<No_First_Y_UpTrij+1;k+)/针对下三角第j列非零元作处理i=Foot_Y_UpTrik;/行足码Num_Y_DownTrii+;/下三角第i行非零元个数增1Row_Downj=No_First_Y_UpTrij;Row_Upj=No_First_Y_UpTrij+1;No_First_Y_DownTri0=0;for(i=0;i<Num_Node-Num_Swing;i+)/下三角各行第一个非零元的序号No_First_Y_DownT

47、rii+1=No_First_Y_DownTrii+Num_Y_DownTrii;for(i=1;i<Num_Node-Num_Swing;i+)/该循环确定下三角各行非零元的/列足码。li=No_First_Y_DownTrii;/下三角第i行第一个非零元序号for(j=0;j<i;j+)/该循环搜寻下三角第0i-1列中/行号为i的非零元。for(k=Row_Downj;k<Row_Upj;k+)if(Foot_Y_UpTrik=i)Foot_Y_DownTrili=j;/记录i行第li个非零元的列足码No_Y_DownTri_RowtoColli=k;/记录该元素在下三角

48、按/列压缩存储序列中序号li+;/序号计数器增1，备下次使用Row_Downj+;break;/稀疏导纳矩阵下三角的结果输出/*cout<<endl<<"下三角 "<<l<<" 个非对角元素的顺序号、列足码及对应"<<"的按列压缩存储时的顺序号"<<endl;for(k=0;k<l;k+)cout<<setw(5)<<k<<setw(5)<<Foot_Y_DownTrik<<setw(5)<&

49、lt;No_Y_DownTri_RowtoColk<<endl; cout<<endl<<"导纳阵下三角每行非对角元素的个数:行号及个数"<<endl;for(i=0;i<Num_Node-Num_Swing;i+)cout<<setw(5)<<i<<setw(5)<<Num_Y_DownTrii<<endl;cout<<endl<<"下三角每行第一个非对角元的顺序号:行号及顺序号"<<endl;for(i

50、=0;i<Num_Node-Num_Swing+1;i+)cout<<setw(5)<<i<<setw(5)<<No_First_Y_DownTrii<<endl; */形成节点导纳矩阵2(包括线路充电容纳及非标准变比的影响)子程void Y_Bus2(int Num_Node,int Num_Line,int Num_Load,int Num_Swing)int i,j,k,k_old,Flag,l;double R,X,Z,G,B,BK;l=0;for(i=0;i<Num_Node-Num_Swing;i+)/初始化Y

51、_Diagi0=0.0;Y_Diagi1=0.0;for(k=0;k<Num_Line;k+)i=Line_NodeNamek0;/线路左节点j=Line_NodeNamek1;/线路右节点if(i>=Num_Node-Num_Swing)/左右节点均为平衡节点，对导纳阵无break;/影响。k_old=Line_No_NewtoOldk;/对应的旧线路顺序号R=LLinek_old.RXBK0;/取线路电阻值X=LLinek_old.RXBK1;/取线路电抗值BK=LLinek_old.RXBK2;/取线路容纳半值或变压器变比值Z=R*R+X*X;G=R/Z;/电导B=-X/Z;/电纳if(j>=Num_Node-Num_Swing)/左为普通节点,右为平衡节点if(Line_Flagk=0)/普通支路Y_Diagi0=Y_Diag