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