統的運行狀態，如母線上的電壓(幅值及相角)、網絡中的功率分布及功率損耗等。1)形成節點導納矩陣(1)自導納的形成(2)互導納的形成3)利用因子表法求解修正方程在電網計算中經常遇到這樣的問題，對方程組需要反復多次求解，而每次求解僅改變常數項F,系數矩陣保持不變。按照一般的高斯消去法，對每一改變的常數項，形成包括常數項及系數矩陣在內的增廣矩陣，然后消去回代求出其解。可以看出，每次對增廣矩陣中A矩陣元素的消元都是重復的，為了避免這種重復，我們把對相同的系數矩陣重復進行的消去與對不同的常數項進行的消去分開進行，因此對系數矩陣的消去只需進行一次，并在消去的過程中將對常數項進行消去運算的運算因子保存下來，形成所謂因子表，這就是因子表法。因為因子表記錄了高斯消去法對常數項進行消去的全部信息，利用它便可對不同常數項進行消去，形成上三角矩陣，最后求出全部未知數。在使用PQ分解法時，其系數矩陣是在迭代過程中保持不變的，所以為了節省內存和縮短運算時間我們采取了因子表法。同時由于電網的節點導納矩陣矩陣是稀疏陣和對稱陣，于是我們可以采取只保存系數矩陣的上三角陣來使運算更為簡化。其中X=[x?x?x?…x,]稱為未知數向量，T4)多次迭代最終求得V和δ以及全線路功率o再利用δ(+1)=8)+△8),V{+1)=V()+△V)求得每次迭代后的結果。多次迭代當*i程序及說明(1)節點功率結構體(2)線路參數結構體類型含義系統總節點數PQ節點數原始節點編號與程序表示編號映射數組節點功率數組V節點電壓有效值數組節點電壓相角數組Y線路參數矩陣有功、無功簡化雅克比矩陣B'、B”PQ迭代次數有功、無功迭代精度控制有功、無功迭代結束標志有功、無功不平衡量數組開始開始結束YNYN求△δYNYN求△V結束/*WRITTENBYHUANGYANGTONGdefinePjiYjicostmpYji{{if((infopenDatatxtNULLErrorMsgYi=*(Y+i)-i;3Noi=No+i;Noi=No+i;【潮流計算結果(節點)】(迭代次數k=%3d)\n",count-1);Noi=No+i;j=*(Num+m);}}j=1;Noi=No+k;fscanf(in,"%d%f%f%fNoiNoiDlta*(Num+i)=k;/*第i個Num元素中存放i節點在No中的下標*/PV節點：%d\n",Node,NQ,NP-NQ);for(i=0;i<Node;i++)/*再在第i行上面開辟Node-ifscanf(in,"%d%d%f%f%f%f',&(*(*(Y+i)).B+=tmp;/*將變壓器的對地充電容納累加到自導納*1if(i>j){tmp=i;i=jj=tmp;kYij.k=k;/*且有B0ji=B0ij/k*/Yij.G=-R;/*生成互導納*1(*(*(Y+i)).G+=R;/*將線路互導納累加到自導納*/3i=*(Num+i);/*查取第i節點在程序中存儲序號*/j=*(Num+m);/*查取第j節點在程序中存儲****************子函數：生成BP、BQ矩陣****************Yi=*(Y+i)-i;)BQ=BP;/*BP包含BQ,BP左上角的NQ*NQ子陣即BQ*{子函數：求因子表for(j=ij<NPj++BPitmpBPkPQVj=*(Num+i);/*查取第i節點在程序中存儲序號*/Noi=No+j;++i,Type,(*Noi).P,(*Noi).Q,*(V+j),*(Dlta+j)//****************子函數：錯誤信息****************/caseprintftErrorItsEman程序編寫過程中，采用了《電力系統分析》教材的4節點、IEEE9節點參數進行了驗證。本說明書中僅列出4節點驗算過程。4節點測試數據來源于參考文獻中的《電力系統分析》(下冊)第61頁例11-5,用于檢驗程序正確性的參考計算結果見67頁。對其基本情況總結如下。32)節點參數電壓相角δ110012002300340003)線路參數節點1121141241310040(1)節點數Node第一列數字2表明該節點的類型為2-PV節點；第二列數字3表明該行數據為節點3的(3)節點類型的判斷3)平衡節點：沒有給定P、Q初始值，僅給定V的初始值，用數字代碼3表示；(1)線路阻抗Z前兩列代表線路兩端的節點編號，亦即線路編號14(即41);第三四列代表線路阻抗(2)變壓器支路2)變壓器線路的編號31有特定含義：3對應節點p,1對應節點q。即變壓器支路的編【節點參數表】節點類型PQUδ總節點：4PQ節點：2pu節點：1【節點導納矩陣Y】(1-042093,-8.242867)(-9.588235,2.352941)(-0.900900,3-666666)(-0.453858,1.891074)【潮流計算結果(節點)】(迭代次數k=4)【、潮流計算結果(線路)】線路