電力系統(tǒng)潮流計算c程序

1、電力系統(tǒng)潮流計算注：這是一個基于N－R法的潮流計算通用程序，僅提供了子程序，需要做些處理才能成為一個可運(yùn)行的計算程序！此程序非我原創(chuàng)，僅與大家共享?。?！這里提供的是電力系統(tǒng)潮流計算機(jī)解法的五個子程序采用的方法是Newton_Raphson法.程序中所用的變量說明如下:N:網(wǎng)絡(luò)節(jié)點總數(shù).M:網(wǎng)絡(luò)的PQ節(jié)點數(shù).L:網(wǎng)絡(luò)的支路總數(shù).N0:雅可比矩陣的行數(shù).N1:N01K:打印開關(guān).K=1則打印否則不打印.K1:子程序PLSC中判斷輸入電壓的形

2、式.K1=1則為極座標(biāo)形式.否則為直角坐標(biāo)形式.D:有功及無功功率誤差的最大值.G(IJ):Ybus的電導(dǎo)元素(實部).B(IJ):Ybus的電納元素(虛部).G1(I):第I支路的串聯(lián)電導(dǎo).B1(I):第I支路的串聯(lián)電納.C1(I):第I支路的pie型對稱接地電納.C(IJ):第I節(jié)點J支路不對稱接地電納.CO(I):第I節(jié)點的接地電納.S1(I):第I節(jié)點的起始節(jié)點號.E1(I):第I節(jié)點的終止節(jié)點號.P(I):第I節(jié)點的注入有功功

3、率.Q(I):第I節(jié)點的注入無功功率.P0(I):第I節(jié)點有功功率誤差.Q0(I):第I節(jié)點無功功率誤差.V0(I):第I節(jié)點(PV節(jié)點)的電壓誤差(平方誤差).intijioi0intpos1pos2intstenif(file4==NULL)fp=stdoutelsefp=file4輸出到文件初始化矩陣GBf(i=1i=ni)f(j=1j=nj)pos2=f2(ijn)g[pos2]=0b[pos2]=0計算支路導(dǎo)納f(i=1i=l

4、i)計算對角元pos1=f1(i)st=s1[pos1]en=e1[pos1]pos2=f2(ststn)g[pos2]=g1[pos1]b[pos2]=b1[pos1]c1[pos1]pos2=f2(enenn)g[pos2]=g1[pos1]b[pos2]=b1[pos1]c1[pos1]計算非對角元pos2=f2(stenn)g[pos2]=g1[pos1]b[pos2]=b1[pos1]g[f2(enstn)]=g[f2(ste