电力系统潮流计算
注这是一个基于NR法的潮流计算通用程序仅提供了子程序需要做些处理才能成为一个可运行的计算程序此程序非我原创仅与大家共享这里提供的是电力系统潮流计算机解法的五个子程序采用的方法是Newto
_Raphso
法程序中所用的变量说明如下N网络节点总数M网络的PQ节点数L网络的支路总数N0雅可比矩阵的行数N1N01K打印开关K1则打印否则不打印K1子程序PLSC中判断输入电压的形式K11则为极座标形式否则为直角坐标形式D有功及无功功率误差的最大值GIJYbus的电导元素实部BIJYbus的电纳元素虚部G1I第I支路的串联电导B1I第I支路的串联电纳C1I第I支路的pie型对称接地电纳CIJ第I节点J支路不对称接地电纳COI第I节点的接地电纳S1I第I节点的起始节点号E1I第I节点的终止节点号PI第I节点的注入有功功率QI第I节点的注入无功功率P0I第I节点有功功率误差Q0I第I节点无功功率误差V0I第I节点PV节点的电压误差平方误
f差VI第I节点的电压误差幅值EI第I节点的电压的实部FI第I节点的电压的虚部JMIJJacoby矩阵的第I行J列元素AIJ修正方程的增广矩阵三角化矩阵的第I行J列元素运算结束后A矩阵的最后一列存放修正的解P1I第I支路由S1I节点注入的有功功率Q1I第I支路由S1I节点注入的无功功率P2I第I支路由E1I节点注入的有功功率Q2I第I支路由E1I节点注入的无功功率P3I第I支路的有功功率损耗Q3I第I支路的无功功率损耗ANGLEI第I节点电压的角度i
cludemathhi
cludestdiohdefi
ef1ii1把习惯的一阶矩阵的下标转化为C语言数组下标defi
ef2ij
i1
j1把习惯的二阶矩阵的下标转化为C语言数组下标本子程序根据所给的支路导纳及有关信息形成结点导纳矩阵如打印参数K1则输出电导矩阵G和电纳矩Bvoidybusi
t
i
tli
tmfloatgfloatbfloatg1floatb1floatc1floatcfloatcoi
tki
ts1i
te1exter
FILEfile4FILEfp
fi
tijioi0i
tpos1pos2i
tste
iffile4NULLfpstdoutelsefpfile4输出到文件初始化矩阵GBfori1i
iforj1j
jpor