實驗目的根據所給的電力系統(tǒng)，編制潮流計算程序，通過計算機進行調試，最后完成一個切實可行的電力系統(tǒng)計算應用程序。通過自己設計電力系統(tǒng)計算程序使同學對電力系統(tǒng)分析有進一步理解，同時加強計算機實際應用能力的訓練。程序計算原理1、概述應用計算機進行電力系統(tǒng)計算，首先要掌握電力系統(tǒng)相應計算的數學模型；其次是運用合理的計算方法；第三則是選擇合適的計算機語言編制計算程序。建立電力系統(tǒng)計算的相關數學模型，就是建立用于描述電力系統(tǒng)相應計算的有關參數間的相互關系的數學方程式。該數學模型的建立往往要突出問題的主要方面，即考慮影響問題的主要因素，而忽略一些次要因素，使數學模型既能正確地反映實際問題，又使計算不過于復雜。運用合理的計算方法，就是要求所選用的計算方法能快速準確地得出正確結果，同時還應要求在解算過程中占用內存少，以利提高計算機的解題規(guī)模。選擇合適的語言編寫程序，就是首先確定用什么計算機語言來編制程序；其次是作出計算的流程圖；第三根據流程圖用選擇的語言編寫計算程序。然后上機調試，直到語法上無錯誤。本程序采用C語言進行編程。所編制的程序難免存在邏輯錯誤，因此先用一個已知結果的系統(tǒng)作為例題進行計算。用程序計算的結果和已知結果相比較，如果結果相差甚遠就要逐步分析程序的計算步驟，查出問題的出處；如果結果比較接近，則逐步分析誤差來源；直到結果正確為止。2、電力系統(tǒng)潮流計算的程序算法潮流計算是電力系統(tǒng)分析中的一種最基本的計算，它的任務是對給定的運行條件確定系統(tǒng)的運行狀態(tài)，如母線上的電壓（幅值及相角）、網絡中的功率分布及功率損耗等。目前計算機潮流計算的方法主要有牛頓-拉夫遜算法和PQ分解法。牛頓-拉夫遜算法是數學上求解非線形方程組的有效方法，具有較好的收斂性，曾經是潮流計算中應用比較普遍的方法。PQ快速分解法是從牛頓-拉夫遜算法演變而來的，是將純數學的牛頓-拉夫遜算法與電力系統(tǒng)具體特點相結合并進行簡化與改進而得出的。PQ快速分解法比牛頓-拉夫遜算法大大提高了計算速度和節(jié)省了內存，故而本程序以PQ快速分解法進行潮流計算。1）形成節(jié)點導納矩陣（1）自導納的形成對節(jié)點i其自導納Yii是節(jié)點i以外的所有節(jié)點都接地時節(jié)點i對地的總導納。顯然，Yii應等于與節(jié)點i相接的各支路導納之和，即Y二y+Xyiii0 ijj式中，yi0為節(jié)點i與零電位節(jié)點之間的支路導納；yij為節(jié)點i與節(jié)點j之間的支路導納。2）互導納的形成對節(jié)點i與節(jié)點k之間的互導納是節(jié)點i、k之間的支路導納的負值，即y二-yikik不難理解Y二Y。若節(jié)點i和k沒有支路直接相連時，便有Yik=°kiik（3）含變壓器支路的處理若節(jié)點p、q間接有變壓器，如下圖所示，則可作出其n型等值電路為：圖1變壓器口型等值電路則p、q的自導納和節(jié)點間的互導納分別為2） 計算不平衡功率AP'AQ并形成修正方程式對每一個PQ節(jié)點或每一個PV節(jié)點都可以根據下列公式計算出有功功率增量AP而對于每一個PQ節(jié)點還可以根據下面的公式計算出無功功率增量AQ在有功功率增量和無功功率增量不滿足如下約束條件時利用PQ分解法則可以形成如下修正方程3） 利用因子表法求解修正方程在電網計算中經常遇到這樣的問題，對方程組需要反復多次求解，而每次求解僅改變常數項F,系數矩陣保持不變。按照一般的高斯消去法，對每一改變的常數項，形成包括常數項及系數矩陣在內的增廣矩陣，然后消去回代求出其解?？梢钥闯觯看螌υ鰪V矩陣中A矩陣元素的消元都是重復的，為了避免這種重復，我們把對相同的系數矩陣重復進行的消去與對不同的常數項進行的消去分開進行，因此對系數矩陣的消去只需進頁腳內容2行一次，并在消去的過程中將對常數項進行消去運算的運算因子保存下來，形成所謂因子表，這就是因子表法。因為因子表記錄了高斯消去法對常數項進行消去的全部信息，利用它便可對不同常數項進行消去，形成上三角矩陣，最后求出全部未知數。在使用PQ分解法時，其系數矩陣是在迭代過程中保持不變的，所以為了節(jié)省內存和縮短運算時間我們采取了因子表法。同時由于電網的節(jié)點導納矩陣矩陣是稀疏陣和對稱陣，于是我們可以采取只保存系數矩陣的上三角陣來使運算更為簡化。aaaLa111213L1naaa22232naLa333nOM其中A=nn稱為系數矩陣，其中A=nn稱為系數矩陣，aX=[xxxLxA稱為未知數向量，123nF=[fffLf｝稱為常數項向量。將矩陣A的元素進行如下處理:123n得到因子表其中D=1a(i-i),U=a(i)(i<j)；ii ii ijij再利用因子表進行前代過程，求出每次迭代后的常數項。其前代公式是求得向量F=rf⑴f⑵f⑶Lf(n)T；1 2 3 n再由因子表與前代得到的向量F，得到方程組求解出此方程即可得到線性方程組的解向量X=[xxxLxb。1 2 3n4）多次迭代最終求得V和6以及全線路功率利用上面所介紹的方法求解修正方程組可以求得VV二[vvVVVVLVV]T和V6=[v6V6V6LV6]T。12 3 m 1 23 n-1再利用6（k+1）=6（k）+V6仏），V（k+1）=V（k）+VV（k）求得每次迭代后的結果。多次迭代當其滿足iiiiii約束條件max｛P（k）｝＜e和max[vQ（k）｝＜e時，迭代結束。迭代結束后即可得到各節(jié)點的Vi P i Q和6，再根據V、6來計算PV節(jié)點的無功功率Q和平衡節(jié)點的功率以及網絡中的功率分布。PV節(jié)點及平衡節(jié)點無功功率計算公式為：平衡節(jié)點有功功率計算公式為：以下圖所標示的正方向，輸電線路功率的計算公式如下圖2支路功率計算對其進行實部虛部進行分解可得P、Q計算公式為：程序及說明1、主要變量說明1）結構體類型說明1）節(jié)點功率結構體structNodetype{floatP,Q}其中，P為節(jié)點的有功功率，Q為無功功率。節(jié)點功率不區(qū)分負荷功率和發(fā)電機功率，其值為本節(jié)點連接的各支路輸入功率及節(jié)點所接負荷、發(fā)電機功率之和，且規(guī)定功率流入節(jié)點為正，流出為負。詳細說明參見下一章“算例及結果”的第二節(jié)“源數據格式說明”。（2）線路參數結構體structLinetype{floatG,B,B0,k}其中，G、B為線路的導納和容納；B0為線路的考慮變壓器n型等值電路后的對地充電容納的一半Bc/2；k為折算到標準變壓器支路后的變壓器變比。詳細說明參見下一章“算例及結果”的第二節(jié)“源數據格式說明”。2）變量說明

表2程序主要變量說明主要變量類型含義Nodeint系統(tǒng)總節(jié)點數NPintPV+PQ節(jié)點數，即非平衡節(jié)點數NQintPQ節(jié)點數Numint*原始節(jié)點編號與程序表示編號映射數組NostructNodetype*節(jié)點功率數組Vfloat*節(jié)點電壓有效值數組Ditafloat*節(jié)點電壓相角數組YstructLinetype**線路參數矩陣BP、BQfloat**有功、無功簡化雅克比矩陣B、B”countunsignedintPQ迭代次數eP、eQconst有功、無功迭代精度控制kp、kqint有功、無功迭代結束標志dP、dQfloat*有功、無功不平衡量數組3、程序流程圖圖3程序主流程圖圖4迭代部分流程圖4、程序代碼/*FUNCTION:POWERFLOW*//*WRITTENBY:HUANG&YANG&TONG*//*LASTEDITED:2008-11-24*/#include<stdio.h>#include<math.h>/***宏定義***/#defineeP0.00001#defineeQ0.00001#defineY_(i,j)(*(*(Y+i)-i+j))#defineYij(*(Yi+j))#defineYji(*(*(Y+j)-j+i))#definePjiYji.G*cos(tmp)+Yji.B*sin(tmp)#definePijYij.G*cos(tmp)+Yij.B*sin(tmp)#defineQjiYji.G*sin(tmp)-Yji.B*cos(tmp)#defineQijYij.G*sin(tmp)-Yij.B*cos(tmp)/***結構體變量定義***/structNodetype/*節(jié)點功率*/{floatP,Q}structLinetype/*線路類型*/{floatG,B,B0,k}/***子函數聲明***/voidin_node()/*讀節(jié)點功率*/voidin_line()/*讀線路參數*/voidB_Form()/*生成BP、BQ矩陣*/voidfactor()/*求因子表*/voidsolve(float**B,float*X,intN)/*解方程組*/voidPrtNode()/*打印節(jié)點參數*/voidErrorMsg(intFlag)/*錯誤提示信息*//***全局變量聲明***/intNode/*節(jié)點數*/int*Num/*保存原始節(jié)點序號*/intNP,NQ=0/*PV+PQ、PQ節(jié)點數*/structNodetype*No/*節(jié)點數據*/structLinetype**Y/*線路參數*/float**BP,**BQ/*有功、無功簡化雅克比矩陣B*/float*V/*節(jié)點電壓有效值*/float*Dlta/*節(jié)點電壓相角值*/unsignedintcount=0/*迭代計數*/inti,j,k,m/*通用下標值*/floattmp/*臨時數據暫存*/char*Type/*節(jié)點類型*/FILE*in,*out/*輸入、輸出文件指針*//**I****I****I**主函數**I****I****I**/intmain(void){intkp=1,kq=1/*P、Q精度判斷：1-不滿足，0-滿足精度要求*/float*dP,*dPi,*dQ,*dQi/*AP、AQ*/floatDltaistructLinetype*YistructNodetype*NoifloattP,tQif((in=fopen("Data.txt","r"))==NULL)ErrorMsg(1)if((out=fopen("out.txt","w"))==NULL)ErrorMsg(2)in_node()/*讀取節(jié)點參數并統(tǒng)計PQ、PV節(jié)點數*/in_line()/*讀取線路參數并形成Y矩陣*/B_Form()/*形成B(BP&BQ)矩陣*/factor。/*求B因子式(仍保存于BP&BQ)*/for(i=Oi<NQi++)*(V+i)=l/*對PQ節(jié)點電壓V賦初值*/dP=(float*)malloc(sizeof(float)*NP)/*AP*/dQ=dP/*AQ*//*AP、AQ不同時存在，故而可共用空間*/loop：/****開始迭代****/if(kp==O&&kq==O)gotoloopEndcount++/*迭代次數加一*/if(count==65535)ErrorMsg(99)/*不收斂，退出*/kp=Ofor(i=Oi<NPi++){dPi=dP+iYi=*(Y+i)-iDltai=*(Dlta+i)*dPi=0for(j=0j<Nodej++){tmp=Dltai-*(Dlta+j)/*tmp即§ij*/if(i>j)*dPi+=*(V+j)*(Pji)else*dPi+=*(V+j)*(Pij)}/*注意到Y矩陣為上三角矩陣，i>j時要交換下標*/*dPi*=*(V+i)*dPi=(*(No+i)).P-*dPi/*求得APi*/if(fabs(*dPi)>0x8fffffff)ErrorMsg(99)/*不收斂，退出*/if(fabs(*dPi)>eP)kp=l/*有不滿足精度的AP即令kp=1*/*dPi/=*(V+i)/*求得常數項APi/Vi*/}if(kp==0)gotoloopQsolve(BP,dP,NP)for(i=Oi<NPi++)*(Dlta+i)+=(*(dP+i)/(*(V+i)))/*修正相角6+=AS*/loopQ:if(kp==O&&kq==O)gotoloopEndkq=Ofor(i=Oi<NQi++){dQi=dQ+iYi=*(Y+i)-iDltai=*(Dlta+i)*dQi=Ofor(j=Oj<Nodej++){tmp=Dltai-*(Dlta+j)/*tmp即6ij*/if(i>j)*dQi+=*(V+j)*(Qji)else*dQi+=*(V+j)*(Qij)}/*注意到y(tǒng)矩陣為上三角矩陣，i>j時要交換下標*/*dQi*=*(V+i)*dQi=(*(No+i)).Q-*dQi/*求得AQi*/if(fabs(*dQi)>0x8fffffff)ErrorMsg(99)/*不收斂，退出*/f(fabs(*dQi)>eQ)kq=l/*有不滿足精度的AQ即令kq=l*/*dQi/=*(V+i)/*求得常數項AQi/Vi*/}if(kq==0)gotoloopsolve(BQ,dQ,NQ)for(i=0i<NQi++)*(V+i)+=*(dQ+i)/*修正PQ節(jié)點電壓V+=AV*/gotoloop/*無功迭代，則必定需要下一輪回迭代判斷*/loopEnd：/****迭代結束****/free(dP)/*釋放內存空間*//****計算PV節(jié)點和平衡節(jié)點的無功功率Q****/for(i=NQi<Nodei++){Noi=No+iYi=*(Y+i)-iDltai=*(Dlta+i)for(j=0j<Nodej++){tmp=Dltai-*(Dlta+j)/*tmp即§ij*/if(i>j)(*Noi).Q+=*(V+j)*(Qji)else(*Noi).Q+=*(V+j)*(Qij)}/*注意到Y矩陣為上三角矩陣，i>j時要交換下標*/(*Noi).Q*=*(V+i)}/****計算平衡節(jié)點的有功功率P****/i=NPNoi=No+iDltai=*(Dlta+i)for(j=0j<Nodej++){tmp=Dltai-*(Dlta+j)/*tmp即§ij*/(*Noi).P+=*(V+j)*(Pji)}/*注意到Y矩陣為上三角矩陣，i>j時要交換下標*/(*Noi).P*=*(V+i)/****輸出最終結果****/fprintf(out,"\n\n【潮流計算結果(節(jié)點)】(迭代次數k=%3d)\n",count-l)PrtNode()/****計算全部線路功率****/fprintf(out,"\n\n【潮流計算結果(線路)】\n")fprintf(out,"線路PQ\n")for(k=0k<Nodek++){i=*(Num+k)Yi=*(Y+i)-iDltai=*(Dlta+i)Noi=No+ifor(m=0m<Nodem++){j=*(Num+m)if(j==i)continuetmp=Dltai-*(Dlta+j)/*tmp即§ij*/if(j<i){f(Yji.B==0)continue/*若Bij=O,則節(jié)點i、j無直接聯(lián)系*/tP=*(V+j)*(Pji)tP=*(V+i)*Yji.G-tPtP*=*(V+i)tQ=-*(V+j)*(Qji)tQ-=*(V+i)*(Yji.B-Yji.B0/Yji.k)tQ*=*(V+i)}else{f(Yij.B==0)continue/*若Bij=O,則節(jié)點i、j無直接聯(lián)系*/tP=*(V+j)*(Pij)tP=*(V+i)*Yij.G-tPtP*=*(V+i)tQ=-*(V+j)*(Qij)tQ-=*(V+i)*(Yij.B-Yij.B0)tQ*=*(V+i)fprintf(out,"S[%d,%d]=(%10.6f,%10.6f)\n",k+1,m+1,-tP,-tQ)}}fclose(out)system("cmd/cstartout.txt")return(0)/**f****f****f**主函數**f****f****f**//******************主函數******************//****************子函數：讀節(jié)點數據****************/voidin_node(){structNodetype*Noi/*臨時變量*/fscanf(in,"%d%d",&Node,&k)/*讀取節(jié)點數Node*/NP=Node-1/*PV+PQ節(jié)點數，即非平衡節(jié)點數目*/Num=(int*)malloc(sizeof(int)*Node)/*開Node個空間，每節(jié)點一個*/V=(float*)malloc(sizeof(float)*Node)/*電壓*/Dlta=(float*)malloc(sizeof(float)*Node)/*電壓相角*/No=(structNodetype*)malloc(sizeof(structNodetype)*Node)/*節(jié)點功率*/j=1while(k!=O)/*若k=0,表明節(jié)點數據讀取完畢*/{switch(k){case1:k=NQNQ++break/*NQ統(tǒng)計PQ節(jié)點個數*/case2:k=NP-jj++break/*從NP-1個空間倒著保存PV節(jié)點*/case3:k=NPbreak/*平衡節(jié)點*/default:ErrorMsg(3)}Noi=No+kfscanf(in,"%d%f%f%f%f",&i,&(*Noi).P,&(*Noi).Q,V+k,Dlta+k)i--/*節(jié)點編號減一，以和程序表達方式兼容*/*(Num+i)=k/*第i個Num元素中存放i節(jié)點在No中的下標*/fscanf(in,"%d",&k)/*讀取節(jié)點類型*/}if(NQ+j!=Node)ErrorMsg(4)/*檢驗節(jié)點數據個數是否夠Node個*/fprintf(out,"【節(jié)點參數表】\n")PrtNode()fprintf(out,"總節(jié)點：％d\nPQ節(jié)點：％d\nPV節(jié)點：％d\n",Node,NQ,NP-NQ)}/************子函數：讀線路數據，并形成節(jié)點導納矩陣Y************/voidin_line(){structLinetype*YifloatR,X,k,Bm=sizeof(structLinetype)Y=(structLinetype**)malloc(m*Node)/*先開Node行,每一個節(jié)點一行*/for(i=0i<Nodei++)/*再在第i行上面開辟Node-i列*/{/*即以上三角存儲Y矩陣*/*(Y+i)=(structLinetype*)malloc(m*(Node-i))Yi=*(Y+i)-ifor(j=ij<Nodej++){Yij.G=Yij.B=Yij.B0=Yij.k=0}/*初始化*/}while(feof(in)==0)/*文件指針到文件末*/{fscanf(in,"%d%d%f%f%f%f",&i,&j,&R,&X,&k,&B)i--j--i=*(Num+i)/*轉換節(jié)點號為該節(jié)點在程序中的儲存編號*/j=*(Num+j)(*(*(Y+i))).B+=B/*將對地充電導納累加到自導納*/(*(*(Y+j))).B+=Bif(k!=1.0){X*=kR=0tmp=(1-k)/X(*(*(Y+i))).B+=tmp/*將變壓器的對地充電容納累加到自導納*/(*(*(Y+j))).B+=-(tmp/k)B=tmpk=-k}if(i>j){tmp=ii=jj=tmpk=1/kB*=k}Yi=*(y+i)-i/*以Yi代替*(Y+i)-i,簡化表達式并避免重復計算*/Yij.BO二B/*保存ijO、jiO對地充電電容到BijO*/Yij.k=k/*且有B0ji=B0ij/k*/tmp=R*R+X*XR/=tmpX/=tmpYij.G=-R/*生成互導納*/Yij.B=X(*(*(Y+i))).G+=R/*將線路互導納累加到自導納*/(*(*(Y+i))).B+=-X(*(*(Y+j))).G+=R(*(*(Y+j))).B+=-X}fclose(in)fprintf(out,"\n【節(jié)點導納矩陣Y】\n")for(k=0k<Nodek++){i=ki=*(Num+i)/*查取第i節(jié)點在程序中存儲序號*/for(j=0j<kj++)fprintf(out,"\t\t\t")for(m=km<Nodem++){j=*(Num+m)/*查取第j節(jié)點在程序中存儲序號*/if(i<j)fprintf(out,"(%10.6f,%10.6f)",Y_(i,j).G,Y_(i,j).B)elsefprintf(out,"(%10.6f,%10.6f)",Y_(j,i).G,Y_(j,i).B)}fprintf(out,"\n")}/****************子函數：生成BP、BQ矩陣****************/voidB_Form()float*BPi,*BQistructLinetype*Yiintsize=sizeof(float)BP=(float**)malloc(size*NP)/*以上三角存儲*/for(i=0i<NPi++)*(BP+i)=(float*)malloc(size*(NP-i))for(i=0i<NPi++){BPi=*(BP+i)-i/*以BPi代替*(BP+i)-i,避免重復計算*/Yi=*(Y+i)-ifor(j=ij<NPj++)*(BPi+j)=Yij.B/*(BPi+j)即相當于BP[i][j]*/}BQ=BP/*BP包含BQ,BP左上角的NQ*NQ子陣即BQ*//****************子函數：求因子表****************/voidfactor(){float*BPi,*BPk,*BQifor(i=0i<NPi++){BPi=*(BP+i)-ifor(k=0k<ik++){BPk=*(BP+k)-ktmp=(*(BPk+i))/(*(BPk+k))for(j=ij<NPj++)(*(BPi+j))-=tmp*(*(BPk+j))}*(BPi+i)=1/(*(BPi+i))for(j=i+1j<NPj++)*(BPi+j)*=*(BPi+i)}/****************子函數：解方程組****************/voidsolve(float**B,float*X,intN){float*Bi,*Xifor(i=0i<Ni++)*(X+i)=-*(X+i)/*對常數項取負*//****對常數列進行前代****/for(i=0i<Ni++){Bi=*(B+i)-iXi=X+ifor(j=i+1j<Nj++)*(X+j)-=*(Bi+j)**Xi*Xi*=*(Bi+i)}/****回代以求解方程組****/for(i=N-1i>=0i--){Bi=*(B+i)-iXi=X+ifor(j=N-1j>ij--)*Xi-=*(Bi+j)**(X+j)}}/****打印節(jié)點參數****/voidPrtNode(){structNodetype*Noifprintf(out,"節(jié)點類型PQV6\n")for(i=0i<Node){j=*(Num+i)/*查取第i節(jié)點在程序中存儲序號*/Noi=No+jif(j<NQ)Type="PQ"elseType="PV"if(j==NP)Type="BS"fprintf(out,"%3d%s%10.6f%10.6f%10.6f%10.6f\n",++i,Type,(*Noi).P,(*Noi).Q,*(V+j),*(Dlta+j)/0.017453)}/****************子函數：錯誤信息****************/voidErrorMsg(intFlag){switch(Flag){case1:printf("\n\tError(1):FailedtoOpenFile\"Data.txt\"!")breakcase2:printf("\n\tError(2):FailedtoCreatFile\"out.txt\"!")breakcase3:printf("\n\tError(3):NodeDataError,PleaseCheck!")breakcase4:printf("\n\tError(4):LackNodeData,PleaseCheck!")breakcase99:printf("k=%d\n\tError(99):It'sEmanative!",count)break}getch()fclose(out)exit(Flag)}算例及結果1、原始數據程序編寫過程中，采用了《電力系統(tǒng)分析》教材的4節(jié)點、IEEE9節(jié)點參數進行了驗證。本說明書中僅列出4節(jié)點驗算過程。4節(jié)點測試數據來源于參考文獻中的《電力系統(tǒng)分析》(下冊)第61頁例11-5，用于檢驗程序正確性的參考計算結果見67頁。對其基本情況總結如下。網絡連接圖圖6算例電力系統(tǒng)網絡圖節(jié)點參數表4算例節(jié)點參數表

節(jié)點類型節(jié)點編號有功功率P無功功率Q電壓V電壓相角811-0.30-0.180012-0.55-0.1300230.501.10034001.0503）線路參數表5算例線路參數表節(jié)點1節(jié)點2阻抗R感抗X變壓器變比k充電容納Bc/2120.10.410.01528140.120.510.01920240.080.410.014133100.30.90909104）數據文件data.txt434001.050230.501.10011-0.30-0.1800圖7數據文件格式2、源數據格式說明源數據有功功率、無功功率、電壓、阻抗、感抗、對地充電容納均以標幺值表示。數據文件必須命名為data.txt且與潮流計算程序放置于同一個文件中。數據文件data.txt包含兩類參數：節(jié)點參數和線路參數。節(jié)點數據塊與線路數據塊之間用數字0作為間隔，即在節(jié)點數據塊結束后，另起一行輸入0，然后再在后面按格式要求錄入線路參數1）節(jié)點參數節(jié)點參數包括：系統(tǒng)節(jié)點數Node、節(jié)點功率（有功P、無功Q）,節(jié)點電壓（有效值頁腳內容35V、相角Delta）參數組織格式：（1）節(jié)點數Node節(jié)點數Node寫在參數文件的開頭，如：4表明為四節(jié)點系統(tǒng)。（2） 功率和電壓P/Q/V/Delta首先給出節(jié)點參數示例：230.501.100第一列數字2表明該節(jié)點的類型為2-PV節(jié)點；第二列數字3表明該行數據為節(jié)點3的參數；后面三列依次為P、Q、V的給定值，給定值為0，表明該項參數未知；第六列為相角非平衡節(jié)點的§即為PQ迭代的初始相角值，平衡節(jié)點的即為給定的相角值。節(jié)點功率為各支路輸入功率之和，且規(guī)定功率流入節(jié)點為正，流出為負。故而負荷功率為負值，發(fā)電機功率為正值。（3） 節(jié)