附录_潮流计算程序.doc

附录 牛顿-拉夫逊法潮流计算程序 本附录所介绍的牛顿-拉夫逊法潮流计算程序采用极坐标形式，其中所涉及的计算公式与第四章中的式(4-42)-(4-59)完全相同，计算程序框图与图4-6基本一致。程序采用C语言编制。为了便于初学者阅读，节点导纳矩阵和雅可比矩阵都用满阵存储而未采用稀疏技巧；但为了适当照顾使用的方便性，在输入数据中对节点编号次序不作任何要求。 下面先介绍输入文件的格式和要求，然后列出程序，最后说明潮流结果的输出。建议读者先从原始数据的输入中了解和熟悉它们在程序中对应的变量、结构体数组及其成员的名称，然后对照图4-6和式(4-42)-(4-59)仔细和耐心地阅读程序，最好能在计算机上亲自实现和进行调试。在调试时，可以用下面给出的对应于[例4-3]系统的输入数据，以及程序运行中得出的中间结果，逐步与[例4-3]中所给出的中间结果进行对比，从而查找错误所在并进行改正。 一、原始数据的输入 程序通过“输入数据.txt”文件输入以下5个数据段。 信息(共6个) 总节点数(变量num_node) 线路和并联电容器总数(变量num_line) 变压器支路总数(变量num_tran) 发电机节点总数(变量num_gene) 负荷节点总数(变量num_load) 节点功率不平衡量的容许误差(变量error) 线路和并联电容器数据(结构体数组line)：每一线路或并联电容器包括5个数据 线路 并联电容器 成员 I侧节点编号 所接节点编号 i J侧节点编号 同上 j Π型等值电路电阻 电容器电阻 a Π型等值电路电抗 电容器电抗(负数) b Π型等值电路一端电纳 0.0 c I侧节点编号和J侧节点编号可以对换；线路和并联电容器之间的次序可以任意，而且允许多条线路或多个电容器并联。 变压器支路数据(结构体数组tran)：每一变压器支路包括5个数据 变压器支路 成员 1侧节点编号 i 2侧节点编号 j 电阻 a 电抗 b 非标准变比 c 变压器电阻、电抗和非标准变比与两侧节点编号之间的关系服从图2-29,即电阻和电抗在1侧而非标准变比在2侧；三绕组变压器需按图2-26化成3个变压器支路，其中3侧变压器支路的非标准变比为1。变压器支路之间的次序可以任意。 发电机节点数据(结构体数组gene)：每一发电机节点包括5个数据 发电机节点 成员 所在节点编号 i 节点种类 j 发出有功功率 a 发出无功功率 b 电压 c 对于PQ节点，节点种类为1，电压可给任意值；对于PV节点，节点种类为-1，发出无功功率可给任意值；对于平衡节点，节点种类为0，发出有功功率和发出无功功率都可给任意值。发电机节点之间的次序可以任意。 负荷节点数据(结构体数组load)：每一负荷节点包括3个数据 负荷节点 成员 所在节点编号 i 吸收有功功率 a 吸收无功功率 b 负荷机节点之间的次序可以任意。 对于[例4-3]中的电力系统，输入数据如下： 4 4 1 2 2 0.00001 4 3 0.260331 0.495868 0.025864 1 4 0.173554 0.330579 0.017243 2 2 0.000000 -20.000000 0.000000 3 1 0.130165 0.247934 0.012932 1 2 0.000000 0.166667 1.128205 4 0 0.0 0.0 1.05 3 -1 0.2 0.0 1.05 2 0.5 0.3 4 0.15 0.1 二、计算程序 /* * * * * * 牛顿-拉夫逊法潮流计算程序 * * * * * */ # include stdio.h # include math.h # include stdlib.h # define pnt 1 // 1 - 输出中间结果 void read_data(); // 输入始数据 void admt_matrix(); // 形成导纳矩阵 void form_Jacobian(); // 形成雅可比矩阵, 计算功率误差 void solv_Eqn(); // 求解修正方程式 void node_flow(); // 输出节点潮流 void branch_flow(); // 输出支路潮流 double **newSpaceDouble2(int,int); void deleteSpac