電力系統(tǒng)潮流計(jì)算C語言程序及說明

1、試驗(yàn)?zāi)康男械碾娏ο到y(tǒng)計(jì)算應(yīng)用程序.通過自己設(shè)計(jì)電力系統(tǒng)計(jì)算程序使同學(xué)對電力系統(tǒng)分析有進(jìn)一步理解，同時(shí)加強(qiáng)計(jì)算機(jī)實(shí)際應(yīng)用力量的訓(xùn)練。程序計(jì)算原理1、概述合理的計(jì)算方法；第三則是選擇合適的計(jì)算機(jī)語言編制計(jì)算程序。簡單。應(yīng)要求在解算過程中占用內(nèi)存少，以利提高計(jì)算機(jī)的解題規(guī)模。算的流程圖；第三依據(jù)流程圖用選擇的語言編寫計(jì)算程序.然后上機(jī)調(diào)試，直到語法上無錯(cuò)誤。本程序接受C 語言進(jìn)行編程。所編制的程序難免存在規(guī)律錯(cuò)誤,因此先用一個(gè)已知結(jié)果的系統(tǒng)作為例題進(jìn)行計(jì)算.用程序計(jì)算的結(jié)果和已知結(jié)果相比較,假如結(jié)果相差甚遠(yuǎn)就要逐步分析程序的計(jì)算步驟，查出問題的出處；假如結(jié)果比較接近,則逐步分析誤差來源;直到結(jié)果正確

2、為止。2、電力系統(tǒng)潮流計(jì)算的程序算法潮流計(jì)算是電力系統(tǒng)分析中的一種最基本的計(jì)算統(tǒng)的運(yùn)行狀態(tài)，如母線上的電壓(幅值及相角、網(wǎng)絡(luò)中的功率分布及功率損耗等.目前計(jì)算機(jī)潮流計(jì)算的方法主要有牛頓拉夫遜算法和 PQ 分解法。牛頓-拉夫遜算法是數(shù)學(xué)上求解非線形方程組的有效方法,具有較好的收斂性,曾經(jīng)是潮流計(jì)算中應(yīng)用比較普遍 的方法。PQ 快速分解法是從牛頓拉夫遜算法演化而來的，是將純數(shù)學(xué)的牛頓拉夫遜算法與電力系統(tǒng)具體特點(diǎn)相結(jié)合并進(jìn)行簡化與改進(jìn)而得出的。PQ 快速分解法比牛頓拉夫遜算法大大提高了計(jì)算速度和節(jié)省了內(nèi)存，故而本程序以PQ快速分解法進(jìn)行潮流計(jì)算.形成節(jié)點(diǎn)導(dǎo)納矩陣（1)自導(dǎo)納的形成對節(jié)點(diǎn)i 其自導(dǎo)納Y

3、ii i 以外的全部節(jié)點(diǎn)都接地季節(jié)點(diǎn)i Yii應(yīng)等于與節(jié)點(diǎn)i 相接的各支路導(dǎo)納之和，即式中,yi0 i 與零電位節(jié)點(diǎn)之間的支路導(dǎo)納；yij 為節(jié)點(diǎn)i 與節(jié)點(diǎn)j 之間的支路導(dǎo)納?；?dǎo)納的形成對節(jié)點(diǎn)i 與節(jié)點(diǎn)k 之間的互導(dǎo)納是節(jié)點(diǎn)i、k 之間的支路導(dǎo)納的負(fù)值，即不難理解。若節(jié)點(diǎn)i和k沒有支路直接相連時(shí),便有Yik=0含變壓器支路的處理若節(jié)點(diǎn)p、q 間接有變壓器，如下圖所示，則可作出其型等值電路為：圖1變壓器型等值電路則p、q 的自導(dǎo)納和節(jié)點(diǎn)間的互導(dǎo)納分別為計(jì)算不平衡功率P、Q 并形成修正方程式對每一個(gè)PQ 節(jié)點(diǎn)或每一個(gè)PV 節(jié)點(diǎn)都可以依據(jù)下列公式計(jì)算出有功功率增量P而對于每一個(gè)PQ 節(jié)點(diǎn)還可以依據(jù)

4、下面的公式計(jì)算出無功功率增量Q在有功功率增量和無功功率增量不滿足如下約束條件時(shí)利用PQ 分解法則可以形成如下修正方程3）利用因子表法求解修正方程 數(shù)項(xiàng) F，系數(shù)矩陣保持不變。依據(jù)一般的高斯消去法，對每一轉(zhuǎn)變的常數(shù)項(xiàng)，形成包括常數(shù) PAGE PAGE 13項(xiàng)及系數(shù)矩陣在內(nèi)的增廣矩陣,A 矩陣元素的消元都是重復(fù)的,為了避開這種重復(fù)，我們把對相同的系數(shù)矩陣重復(fù)進(jìn)行的消去與對 中將對常數(shù)項(xiàng)進(jìn)行消去運(yùn)算的運(yùn)算因子保存下來,形成所謂因子表,這就是因子表法。由于因子表記錄了高斯消去法對常數(shù)項(xiàng)進(jìn)行消去的全部信息，利用它便可對不同常數(shù)項(xiàng)進(jìn)行消去, 形成上三角矩陣，最終求出全部未知數(shù).在使用PQ 分解法時(shí)，其系數(shù)

5、矩陣是在迭代過程中保持不變的，所以為了節(jié)省內(nèi)存和縮 是我們可以實(shí)行只保存系數(shù)矩陣的上三角陣來使運(yùn)算更為簡化。若線性方程組一般形式如下：其中稱為系數(shù)矩陣,稱為未知數(shù)向量， 稱為常數(shù)項(xiàng)向量。將矩陣A 的元素進(jìn)行如下處理:得到因子表其中 ；再利用因子表進(jìn)行前代過程，求出每次迭代后的常數(shù)項(xiàng)。其前代公式是: 求得向量；再由因子表與前代得到的向量F，得到方程組求解出此方程即可得到線性方程組的解向量。4）V 和以及全線路功率利用上面所介紹的方法求解修正方程組可以求得和。V VPV Q 和平衡節(jié)點(diǎn)的功率以及網(wǎng)絡(luò)中的功率分布.PV 節(jié)點(diǎn)及平衡節(jié)點(diǎn)無功功率計(jì)算公式為:平衡節(jié)點(diǎn)有功功率計(jì)算公式為：以下圖所標(biāo)示的正方

6、向，輸電線路功率的計(jì)算公式如下:圖2支路功率計(jì)算對其進(jìn)行實(shí)部虛部進(jìn)行分解可得P、Q 計(jì)算公式為：程序及說明1、主要變量說明結(jié)構(gòu)體類型說明節(jié)點(diǎn)功率結(jié)構(gòu)體struct Nodetypefloat P，Q；;其中，P 為節(jié)點(diǎn)的有功功率,Q 為正,。線路參數(shù)結(jié)構(gòu)體struct Linetypefloat G,B，B0，k;其中，GB 為線路的導(dǎo)納和容納；B0 為線路的考慮變壓器型等值電路后的對地充電容納的一半Bc/2；k 為折算到標(biāo)準(zhǔn)變壓器支路后的變壓器變比。具體說明參見下一章“算例及結(jié)果的其次節(jié)“源數(shù)據(jù)格式說明變量說明表2程序主要變量說明主要變量Node主要變量NodeNP NQNumNo VDlt

7、aYBP、BQcount eP、eQkp、kqdP、dQ類型含義int int int int*struct Nodetype *float *float struct Linetype *float *unsigned intconstintfloat *系統(tǒng)總節(jié)點(diǎn)數(shù)PV+PQ 節(jié)點(diǎn)數(shù),即非平衡節(jié)點(diǎn)數(shù)PQ 節(jié)點(diǎn)數(shù)原始節(jié)點(diǎn)編號與程序表示編號映射數(shù)組節(jié)點(diǎn)功率數(shù)組節(jié)點(diǎn)電壓有效值數(shù)組節(jié)點(diǎn)電壓相角數(shù)組 線路參數(shù)矩陣有功、無功簡化雅克比矩陣B、B”PQ 迭代次數(shù)有功、無功迭代精度把握有功、無功迭代結(jié)束標(biāo)志有功、無功不平衡量數(shù)組開頭開頭讀取節(jié)點(diǎn)數(shù)據(jù)讀取線路數(shù)據(jù)并形成Y 矩陣形成矩陣B及B并求取因子表設(shè)定P

8、Q 節(jié)點(diǎn)電壓初值及各節(jié)點(diǎn)電壓相角的初值迭代計(jì)算計(jì)算平衡節(jié)點(diǎn)功率及全部線路功率結(jié)束圖3程序主流程圖開頭開頭Kp=1，Kq=1 Kq=0Count+Kp=0計(jì)算P，P/NKp=0 ?求+= Kq=0Kq=0計(jì)算Q、Q/NKq=0？求V結(jié)束V+=V圖4迭代部分流程圖4、程序代碼/*FUNCTION : POWER FLOW*/*/WRITTEN BY : HUANGYANG&TONGLAST EDITED : 20081124*/*/include #include math。h/* 宏定義 */#define eP 0.00001#define eQ 0.00001#define Y_(i,j)

9、(*（*(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.Gcos(tmp）+Yij。Bsin(tmp) #define Qji Yji。G*sin（tmp）-Yji。B*cos(tmp）define Qij Yij。Gsin（tmp）Yij。Bcos（tmp)/*結(jié)構(gòu)體變量定義 */ struct Nodetype / 節(jié)點(diǎn)功率 /float P，Q;struct Linetype /* 線路類型 */float G，B,

10、B0,k;/*子函數(shù)聲明 */void in_node()； / 讀節(jié)點(diǎn)功率 /void in_line（）;/ 讀線路參數(shù) /void B_Form（）;/* BP、BQ 矩陣 */void factor（); / 求因子表 /void solve(float *B，float *X,int N); / 解方程組 */void PrtNod（； */void ErrorMsg(int Flag）; / 錯(cuò)誤提示信息 /* 全局變量聲明 */int Node；/節(jié)點(diǎn)數(shù)*/int *Num；/*保存原始節(jié)點(diǎn)序號*/ int NP,NQ=0;/PV+PQ、PQ 節(jié)點(diǎn)數(shù)*/struct Nodety

11、pe No； /*節(jié)點(diǎn)數(shù)據(jù)/struct Linetype *Y； /線路參數(shù)/float *BP,*BQ; /有功、無功簡化雅克比矩陣B*/float V;/節(jié)點(diǎn)電壓有效值*/float Dlta；/節(jié)點(diǎn)電壓相角值*/unsigned int count=0； /*迭代計(jì)數(shù)/ int i，j,k,m；/*通用下標(biāo)值*/float tmp；/臨時(shí)數(shù)據(jù)暫存*/char Type； /*節(jié)點(diǎn)類型*/FILE in，*out；/*輸入、輸出文件指針/*主函數(shù) */*主函數(shù) */int main（void）int kp=1，kq=1; /*P、Q 精度推斷:1不滿足，0滿足精度要求*/float *d

12、P，*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(2); in_node(）;/*讀取節(jié)點(diǎn)參數(shù)并統(tǒng)計(jì)PQ、PV 節(jié)點(diǎn)數(shù)*/ in_line； 讀取線路參數(shù)并形成Y 矩*/B_Form(）;/形成B（BP&BQ)矩陣*/ factor；/求B 因子仍保存于for（i=0;iNQ；i+) (V+i）=1；

13、 /*對PQ 節(jié)點(diǎn)電壓V 賦初值*/dP=（float )malloc（sizeof（float）NP);/P/dQ=dP； /Q/*P、Q 不同時(shí)存在，故而可共用空間*/loop：/*開頭迭代*/if(kp=0&kq=0） goto loopEnd；count+; /*迭代次數(shù)加一*/if（count=65535) ErrorMsg（99)； /不收斂，退出*/kp=0；for（i=0;iNP;i+）dPi=dP+i; Yi=*（Y+i)i； Dltai=（Dlta+；*dPi=0; for(j=0；jj) *dPi+=(V+j)*(Pji）;else *dPi+=*（V+j)*(Pij)；

14、/留意到Y(jié) 矩陣為上三角矩陣,ij 時(shí)要交換下標(biāo)*/*dPi=*（V+i）;*dPi=（*（No+i） .P-*dPi； /求得Pi*/ if（fabs(dPi）0 x8fffffff） ErrorMsg（99）;/不收斂，退出*/if(fab（dPeP) kp=1; /有不滿足精度的P 即令kp=1*/dPi/=*（V+i)； /求得常數(shù)項(xiàng)Pi/Vi/if(kp=0） goto loopQ;solve(BP，dP，NP）;for(i=0；iNP;i+) （Dlta+i）+=(（dP+i）/(*(V+i); /修正相角+=*/loopQ:if（kp=0&kq=0） goto loopEnd;k

15、q=0;for（i=0；iNQ；i+）dQi=dQ+i； Yi=*（Y+i)i; Dltai=（Dlta+i)；dQi=0；for(j=0;jj 時(shí)要交換下標(biāo)*/*dQi*=（V+i)；*dQi=（*（No+i） .QdQi; /*求得Qi*/ if(fabs（dQi)0 x8fffffff) ErrorMsg(99)； /*不收斂，退出/ if（fabs（dQi)eQ） kq=1； /*有不滿足精度的Q 即令kq=1/*dQi/=(V+i）; /*求得常數(shù)項(xiàng)Qi/Vi*/if(kq=0) goto loop；solve(BQ,dQ,NQ);for（i=0；ij） (*Noi).Q+=*(V+

16、j)*（Qji)； else (*Noi）.Q+=*（V+j）(Qij）;/*Y 矩陣為上三角矩陣，ij 時(shí)要交換下標(biāo)*/NoiQ*=（V+i；/*計(jì)算平衡節(jié)點(diǎn)的有功功率P*/ i=NP；Noi=No+i; Dltai=*（Dlta+i); for（j=0；jNode；j+）tmp=Dltai*（Dlta+j);/*tmp 即ij/（*NoiP+V+jPji；/留意到Y(jié) 矩陣為上三角矩陣，ij 時(shí)要交換下標(biāo)/（NoiP*V+i);/*輸出最終結(jié)果/fprintf（out,”nn【 潮流計(jì)算結(jié)果（節(jié)點(diǎn)) 】 （ 迭代次數(shù)k=3d )n，count1)；PrtNode()；/*計(jì)算全部線路功率*/

17、fprintf(ou”n【 潮流計(jì)算結(jié)果（線路）】 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)if(Yji。B=0） continue; /*Bij=0，則節(jié)點(diǎn)i、j 無直接聯(lián)系/tP=*(V+j)*(Pji）;tP=*(V+i)Yji.GtP;tP*（V+；tQ=-(V+j)(Qji)； tQ=*（V+i）

18、*（Yji.BYji.B0/Yji.k)；tQ*=（V+i）;elseif（Yij.B=0） continue; /若Bij=0,則節(jié)點(diǎn)i、j 無直接聯(lián)系/tP=*(V+j)*(Pij)；tP=(V+i)Yij。GtP；tP=（V+i）; tQV+j)*(Qi；tQ-=(V+i)（Yij.B-Yij.B0);tQ=(V+i)；fprintf（out，”Sd，%d= （10.6f，%10。6f)n”,k+1,m+1，tP，-tQ);fclos（ou；d c t out.txt；return（0);/*主函數(shù) */* 主函數(shù) */* 子函數(shù)：讀節(jié)點(diǎn)數(shù)據(jù) */void in_node()struct

19、 Nodetype *Noi; /*臨時(shí)變量/fscanf(in,%d %d”,&Node，&k);/讀取節(jié)點(diǎn)數(shù)Node*/NP=Node-1; /*PV+PQ 節(jié)點(diǎn)數(shù)，即非平衡節(jié)點(diǎn)數(shù)目/Num=(int ）malloc(sizeof（int)*Node）; /*Node 個(gè)空間,每節(jié)點(diǎn)一個(gè)/V=(float )malloc(sizeof(float)Node）; /電壓*/Dlta=（float ）malloc(sizeof（float)Node); /電壓相角/No=（struct Nodetype *）malloc（sizeof（struct Nodetype）*Node)；/節(jié)點(diǎn)功率*

20、/j=1；while(k!=0） /*k=0，表明節(jié)點(diǎn)數(shù)據(jù)讀取完畢*/switch（k）case 1:k=NQ;NQ+；break;/NQ 統(tǒng)計(jì)PQ 節(jié)點(diǎn)個(gè)數(shù)/case 2:k=NP-j；j+;break; /從NP-1 個(gè)空間倒著保存PV 節(jié)點(diǎn)/case 3:k=NP;break； /平衡節(jié)點(diǎn)/default:ErrorMsg(3);Noi=No+k;fscan（d %f %f %f %f”,&，(*Noi).，（Noi）.V+Dlta+；i-;/*節(jié)點(diǎn)編號減一，以和程序表達(dá)方式兼容/*（Num+i)=k； /第i Num 元素中存放i 節(jié)點(diǎn)在No 中的下標(biāo)/fscanf（in，d,&k)；

21、 /讀取節(jié)點(diǎn)類型/if(NQ+j！=Node) ErrorMsg(4); /*Node 個(gè)/fprint（out,【 節(jié)點(diǎn)參數(shù)表】）;PrtNod（；fprintf(out，總節(jié)點(diǎn):%dnPQ 節(jié)點(diǎn)：%dnPV 節(jié)點(diǎn)：dn”,Node,NQ，NPNQ）;/* 子函數(shù)：讀線路數(shù)據(jù)，并形成節(jié)點(diǎn)導(dǎo)納矩陣Y*/void in_line（)struct Linetype Yi；float R,X，k,B; m=sizeof（struct Linetype)；Y=(struct Linetype ）malloc（m*Node); /*先開Node 行，每一個(gè)節(jié)點(diǎn)一行/ for（i=0；iNode；i+）/

22、*再在第i 行上面開拓Nodei 列*/*即以上三角存儲 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，簡化表達(dá)式并避開重復(fù)計(jì)算 */Yij.B0=B； /*ij0、ji0 對地充電電容到Bij0/Yij。k=k；/且有B0ji=B0ij/k/tmp=RR+X*X;R/=tmp;X/=tmp；Yij.G=-R； / 生成互導(dǎo)納 */Yij。B=X；（(*(Y+i） ).G+=R; /* 將線路互導(dǎo)

23、納累加到自導(dǎo)納 /（*(Y+i)）.B+=-X;(*(*(Y+j)） .G+=R；(*（Y+j).B+=X；fclose(in）;fprintf(ou，【 節(jié)點(diǎn)導(dǎo)納矩陣Yn；for(k=0；kNode;k+）i=k;i=（Num+i）; /查取第i 節(jié)點(diǎn)在程序中存儲序號*/for(j=0;jk;j+) fprintf（out，ttt”)； for（m=k;mNode；m+）j=（Num+m）; /查取第j 節(jié)點(diǎn)在程序中存儲序號/if(） fprintf(ou，（10.6，%16f) ，Y）.Y（i,jB； else fprintf(out，”(10。6f,10.6f） ，Y_（j，i).G,Y

24、_（j,i).B)；fprintf(out,n)；/*子函數(shù)：生成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，避開重復(fù)計(jì)算 /Yi=*(Y+i）i；for(j=i;jNP;j+） *(BPi+j）=Yij

25、。B； /*(BPi+j)即相當(dāng)于BPij/BQ=BP; / BP 包含BQ,BP 左上角的NQ*NQ 子陣即BQ */*子函數(shù)：求因子表 */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+j）*=(BPi+i）;/*子函數(shù)：解方程組 */vo

26、id solve（float *B，float X,int N）float *Bi,Xi；for（i=0;iN;i+） *(X+i）=-*（X+i）;/*對常數(shù)項(xiàng)取負(fù)/*對常數(shù)列進(jìn)行前代*/ 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=N1；ji;j-) *Xi-=(Bi+j）*(X+j）;/*打印節(jié)點(diǎn)參數(shù)/void PrtNode()struct Nodetype Noi；fprint（ou”節(jié)點(diǎn) 類型PQV for（i=0;iNode；)j=*(Num+； /查取第i Noi=No+j;if

27、(jNQ） Type=PQ; else Type=”PV；if(j=NP) Type=”BS”;fprintf（out,%3ds%10。6f %10.6f %10。6f %10.6fn”,+iType,(*NoP,Noi，*(V+j)(Dlta+）/0.01745）;/*子函數(shù)：錯(cuò)誤信息 */void ErrorMsg（int Flag）switch（Flag）case 1:printf(ntError（1)： Failed to Open File Data.txt”！”)；break； e 2:printf(ntErro2): Failed to Creat File outxt！brea

28、k; 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(ou； exit(Fla；算例及結(jié)果1、原始數(shù)據(jù)4 IEEE9 節(jié)點(diǎn)參數(shù)進(jìn)行了驗(yàn)證。4 節(jié)點(diǎn)驗(yàn)算過程.4 節(jié)點(diǎn)測試數(shù)據(jù)來源于參考文獻(xiàn)中的電力系統(tǒng)分析(下冊）61 115，用于67

29、 頁.對其基本狀況總結(jié)如下。網(wǎng)絡(luò)連接圖2）節(jié)點(diǎn)參數(shù)表4算例節(jié)點(diǎn)參數(shù)表圖6算例電力系統(tǒng)網(wǎng)絡(luò)圖節(jié)點(diǎn)類型節(jié)點(diǎn)編號PQV電壓相角110。30-0.1800120.550。1300230.501。10034001。0503）線路參數(shù)5算例線路參數(shù)表12RXkBc/2120.10.410。01528140.120.510.01920240.080.410.014133100.30.90909104)data。txt434001.050230.501.10011-0.30-0.180012-0.55-0.130002、源數(shù)據(jù)格式說明圖7數(shù)據(jù)文件格式131221440.100.120.080.41131221

30、440.100.120.080.410.015280.30.90909100.510.019200.410.01413包含兩類參數(shù)：節(jié)點(diǎn)參數(shù)和線路參數(shù)。節(jié)點(diǎn)數(shù)據(jù)塊與線路數(shù)據(jù)塊之間用數(shù)字 0 作為間隔,即在節(jié)點(diǎn)數(shù)據(jù)塊結(jié)束后，另起一行輸入0，然后再在后面按格式要求錄入線路參數(shù)1）節(jié)點(diǎn)參數(shù)節(jié)點(diǎn)參數(shù)包括：系統(tǒng)節(jié)點(diǎn)數(shù)Node、節(jié)點(diǎn)功率（有功P、無功Q）,節(jié)點(diǎn)電壓（有效值V、相角Delta）參數(shù)組織格式:Node節(jié)點(diǎn)數(shù)Node 寫在參數(shù)文件的開頭,如：4 表明為四節(jié)點(diǎn)系統(tǒng)。P/Q/V/Delta首先給出節(jié)點(diǎn)參數(shù)示例：230。501.1002 表明該節(jié)點(diǎn)的類型為2-PV 節(jié)點(diǎn);其次列數(shù)字3 表明該行數(shù)據(jù)為節(jié)點(diǎn)3 的參數(shù)；后面三列依次為P、Q、V 的給定值，給定值為0,表明該項(xiàng)參數(shù)未知；第六列為相角平衡節(jié)點(diǎn)的即為PQ 迭代的初始相角值，平衡節(jié)點(diǎn)的即為給定的相角值。節(jié)點(diǎn)功率為各支路輸入功率之和,且規(guī)定功率流入節(jié)點(diǎn)為正，流出為負(fù)。故而負(fù)荷功率為負(fù)值，發(fā)電機(jī)功率為正值。節(jié)點(diǎn)類型的推斷依據(jù)節(jié)點(diǎn)的給定參數(shù)可以將節(jié)點(diǎn)分為三種類型：1）PQ 節(jié)點(diǎn)：給定P、Q 1 表示；PV 節(jié)點(diǎn)：給定P、V