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

1、電力系統(tǒng)潮流計(jì)算注：這是一個(gè)基于N－R法的潮流計(jì)算通用程序，僅提供了子程序，需要做些處理才能成為一個(gè)可運(yùn)行的計(jì)算程序！此程序非我原創(chuàng)，僅與大家共享?。?！這里提供的是電力系統(tǒng)潮流計(jì)算機(jī)解法的五個(gè)子程序采用的方法是Newton_Raphson法.程序中所用的變量說明如下:N:網(wǎng)絡(luò)節(jié)點(diǎn)總數(shù).M:網(wǎng)絡(luò)的PQ節(jié)點(diǎn)數(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)元素(實(shí)部).B(IJ):Ybus的電納元素(虛部).G1(I):第I支路的串聯(lián)電導(dǎo).B1(I):第I支路的串聯(lián)電納.C1(I):第I支路的pie型對稱接地電納.C(IJ):第I節(jié)點(diǎn)J支路不對稱接地電納.CO(I):第I節(jié)點(diǎn)的接地電納.S1(I):第I節(jié)點(diǎn)的起始節(jié)點(diǎn)號.E1(I):第I節(jié)點(diǎn)的終止節(jié)點(diǎn)號.P(I):第I節(jié)點(diǎn)的注入有功功

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

4、i)計(jì)算對角元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]計(jì)算非對角元pos2=f2(stenn)g[pos2]=g1[pos1]b[pos2]=b1[pos1]g[f2(enstn)]=g[f2(ste