電力潮流計(jì)算C源程序

1、電力潮流計(jì)算 C 源程序* FLOW.C */* 這里提供的是電力系統(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:N0+1 K: 打印開關(guān) .K=1, 則打印 ;否則 , 不打印 .* K1: 子程序 PLSC 中判斷輸入電壓的形式 .K1=1, 則為極座標(biāo)形式 .否則 * 為直角坐標(biāo)形式 . * D: 有功及無功功率誤差的最大值 . * G(I,J):Ybus 的電導(dǎo)元素 ( 實(shí)部 ).

2、* B(I,J):Ybus 的電納元素 ( 虛部 ). * G1(I) : 第 I 支路的串聯(lián)電導(dǎo) . B1(I): 第 I 支路的串聯(lián)電納 . * C1(I) : 第 I 支路的 pie 型對稱接地電納 . * C(I,J): 第 I 節(jié)點(diǎn) J 支路不對稱接地電納 . * CO(I) : 第 I 節(jié)點(diǎn)的接地電納 . * S1(I) : 第 I 節(jié)點(diǎn)的起始節(jié)點(diǎn)號(hào) . E1(I): 第 I 節(jié)點(diǎn)的終止節(jié)點(diǎn)號(hào) . * P(I) : 第 I 節(jié)點(diǎn)的注入有功功率 . Q(I): 第 I 節(jié)點(diǎn)的注入無功功率 .* P0(I) : 第 I 節(jié)點(diǎn)有功功率誤差 . Q0(I): 第 I 節(jié)點(diǎn)無功功率誤差 .

3、* V0(I) : 第 I 節(jié)點(diǎn) (PV 節(jié)點(diǎn) ) 的電壓誤差 ( 平方誤差 ). * V(I) : 第 I 節(jié)點(diǎn)的電壓誤差幅值 . * * E(I) : 第 I 節(jié)點(diǎn)的電壓的實(shí)部 . F(I): 第 I 節(jié)點(diǎn)的電壓的虛部 . * JM(I,J):Jacoby 矩陣的第 I 行 J 列元素 . * * A(I,J): 修正方程的增廣矩陣 ,三角化矩陣的第 I 行 J 列元素 , 運(yùn)算結(jié) * 束后 A 矩陣的最后一列存放修正的解 . * P1(I) : 第 I 支路由 S1(I) 節(jié)點(diǎn)注入的有功功率 . * Q1(I) : 第 I 支路由 S1(I)節(jié)點(diǎn)注入的無功功率 . * P2(I) : 第

4、 I 支路由 E1(I)節(jié)點(diǎn)注入的有功功率 . * Q2(I) : 第 I 支路由 E1(I) 節(jié)點(diǎn)注入的無功功率 . * P3(I) : 第 I 支路的有功功率損耗 . * Q3(I) : 第 I 支路的無功功率損耗 . * ANGLE(I): 第 I 節(jié)點(diǎn)電壓的角度 . * */#include #include #define f1(i) (i-1)/* 把習(xí)慣的一階矩陣的下標(biāo)轉(zhuǎn)化為 C 語言數(shù)組下標(biāo) */#define f2(i,j,n) (i-1)*(n)+j-1)/* 把習(xí)慣的二階矩陣的下標(biāo)轉(zhuǎn)化為 C 語言數(shù)組下標(biāo) */ /* 本子程序根據(jù)所給的支路導(dǎo)納及有關(guān)信息,形成 結(jié)點(diǎn) *

5、導(dǎo)納矩陣 ,如打印參數(shù) K=1, 則輸出電導(dǎo)矩陣 G 和電納矩 B * */void ybus(int n,int l,int m,float *g,float *b,float *g1,float *b1,float *c1,float *c,float *co,int k,int *s1,int *e1)extern FILE *file4;FILE *fp;int i,j,io,i0;int pos1,pos2;int st,en;if(file4=NULL)fp=stdout;elsefp=file4; /* 輸出到文件 */* 初始化矩陣 G,B */for(i=1;i=n;i+)for(j=1;j=n;j+)pos2=f2(i,j,n);gpos2=0;bpos2=0;/* 計(jì)算支路導(dǎo)納 */ for(i=1;i=l;i+)/* 計(jì)算對角元 */pos1=f1(i);st=s1pos1;en=e1pos1;pos2=f2(st,st,n);gpos2+=g1pos1;bpos2+=b1pos1+c1pos1;pos2=f2(en,en,n);gpos2+=g1pos1; bpos2+=b1pos1+c1pos1;/* 計(jì)算非對角元 */pos2=f2(st,en,n);gpos2-=g1pos1;bpos