电力系统潮流计算包含牛顿拉夫逊，PQ分解法

************************************************************************** 电力系统潮流计算---- 牛顿拉夫逊法（直角坐标） ************************************************************************** No.1 ------------------------------ ** 网络中节点信息输出： 网络中节点总数：3 网络中PQ节点数：1 网络中PV节点数：1 网络中平衡节点个数：1 误差精度：0.000100 newCode oldCode type P Q U angle 0 3 0 0.0000 0.0000 1.0000 0.0000 1 1 1 -0.5000 0.4000 1.0000 0.0000 2 2 2 0.3000 0.0000 1.0200 0.0000 ** 节点信息说明： 1、newCode: 计算时用的新节点编号——第一个（0号节点）是平衡节点; 第1~m个是PQ节点; 第m+1~n个是PV节点. 2、oldCode: 用户自己定义的节点编号. 3、type: 节点类型——0 表示平衡节点，1 表示PQ节点，2 表示PV节点 4、P —— 注入母线（节点）的有功功率(MW)；Q —— 注入母线的无功功率(MVar)；U —— 母线（节点）的电压幅值(KV)；Angle —— 母线（节点）的电压辐角(弧度). No.2 --------------------------- ** 网络中平行支路信息： 网络中平行总数：2 Code， (ni,nj) (oi,oj) G+jB R+jX k 1 ( 1, 2)( 1, 2) 0.0000+j 0.0000, 0.0400+j 0.2000, 1.0000 2 ( 1, 0)( 1, 3) 0.0000+j 0.0000, 0.1000+j 0.4000, 1.0000 --------------------------- 网络中平行支路信息（各个阻抗变换成导纳后）： 网络中平行总数：2 Code， (ni,nj) (oi,oj) G+jB G+jB k 2 ( 1, 0)( 1, 3) 0.0000+j 0.0000, 0.5882+j -2.3529, 1.0000 1 ( 1, 2)( 1, 2) 0.0000+j 0.0000, 0.9615+j -4.8077, 1.0000 --------------------------- 网络中平行支路信息（计算机模型）： 网络中平行总数：2 Code， (ni,nj) (oi,oj) G+jB G+jB G+jB k 2 ( 1, 0)( 1, 3) 0.0000+j 0.0000, 0.5882+j -2.3529, 0.0000+j 0.0000, 1.0000 1 ( 1, 2)( 1, 2) 0.0000+j 0.0000, 0.9615+j -4.8077, 0.0000+j 0.0000, 1.0000 ** 普通平行支路数据说明： 1、Code —— 用户定义的普通支路编号(方便用户查询信息)，高压侧为前节点. 2、(ni,nj) —— 程序运行时新定义的支路两端节点编号，和节点新编号一致. 3、(oi,oj) —— 用户自己定义的支路两端节点编号，和节点旧编号一致，高压侧为前节点. 4、R+jX(G+jB) —— 阻抗(导纳) 5、k —— 变比，普通平行支路为1，变压器为高压侧对低压侧的电压比 No.3 --------------------------- ** 网络中母线（节点）对地支路信息： * 网络中对地支路总数：0 newCode oldCode num SignZY R+jX(G+jB) --------------------------- * 网络中母线（节点）对地支路信息（阻抗变导纳后）： newCode oldCode num G+jB * 网络中母线（节点）对地支路信息（同一个节点的导纳合并之后）: newCode oldCode num G+jB **母线（节点）对地支路信息说明： 1、newCode —— 节点新编号，同上. 2、oldCode —— 节点旧编号，同上. 3、num —— 节点上连接对地支路数目. 4、SignZY —— 节点上连接对地支路类型，阻抗为 1，导纳为 2. 5、R+jX(G+jB) —— 节点上连接对地支路阻抗或者是导纳(是前面类型而定). 6、G+jB —— 节点上连接对地支路导纳. No.4 ---------------------------------- 添加平行支路后支路导纳矩阵： G_bet: 0.0000 0.5882 0.0000 0.5882 0.0000 0.9615 0.0000 0.9615 0.0000 B_bet: 0.0000 -2.3529 0.0000 -2.3529 0.0000 -4.8077 0.0000 -4.8077 0.0000 添加平行支路后节点导纳矩阵： G: 0.5882 -0.5882 0.0000 -0.5882 1.5498 -0.9615 0.0000 -0.9615 0.9615 B: -2.3529 2.3529 0.0000 2.3529 -7.1606 4.8077 0.0000 4.8077 -4.8077 No.5 ---------------------------------- 添加节点接地支路后支路导纳矩阵： G_bet: 0.0000 0.5882 0.0000 0.5882 0.0000 0.9615 0.0000 0.9615 0.0000 B_bet: 0.0000 -2.3529 0.0000 -2.3529 0.0000 -4.8077 0.0000 -4.8077 0.0000 添加节点接地支路后节点导纳矩阵： G: 0.5882 -0.5882 0.0000 -0.5882 1.5498 -0.9615 0.0000 -0.9615 0.9615 B: -2.3529 2.3529 0.0000 2.3529 -7.1606 4.8077 0.0000 4.8077 -4.8077 ================================================== 开始迭代运算求解 ================================================== ----------- 第 1 次迭代 ----------- 输出不平衡量： PQ节点： ΔP ΔQ -0.480769 0.496154 PV节点： ΔP ΔU2 0.280385 0.000000 第二象限雅可比矩阵: 7.2568 1.5305 -1.5690 7.0645 第一象限雅可比矩阵: -4.8077 -0.9615 0.9615 -4.8077 第三象限雅可比矩阵: -4.9038 -0.9808 0.0000 0.0000 第四象限雅可比矩阵: 4.8077 1.0000 0.0000 2.0400 **雅可比增广矩阵: 7.2568 1.5305 -4.8077 -0.9615 -0.4808 -1.5690 7.0645 0.9615 -4.8077 0.4962 -4.9038 -0.9808 4.8077 1.0000 0.2804 0.0000 0.0000 0.0000 2.0400 0.0000 ===================================== 输出修正量： △e △f 0.052714 -0.097480 0.000000 -0.030356 ===================================== ** 修正： e f 1.052714 -0.097480 1.020000 -0.030356 ================================== Q:修正后电压的极坐标形式： newCode oldCode U ∠ δ 1 1 1.057217∠-0.092336 2 2 1.020452∠-0.029752 ----------- 第 2 次迭代 ----------- 输出不平衡量： PQ节点： ΔP ΔQ -0.023881 -0.072176 PV节点： ΔP ΔU2 0.009652 -0.000922 第二象限雅可比矩阵: 6.9838 1.8399 -2.8191 7.7902 第一象限雅可比矩阵: -4.9674 -1.4809 1.4809 -4.9674 第三象限雅可比矩阵: -4.8747 -1.1267 0.0000 0.0000 第四象限雅可比矩阵: 5.0965 1.4180 -0.0607 2.0400 **雅可比增广矩阵: 6.9838 1.8399 -4.9674 -1.4809 -0.0239 -2.8191 7.7902 1.4809 -4.9674 -0.0722 -4.8747 -1.1267 5.0965 1.4180 0.0097 0.0000 0.0000 -0.0607 2.0400 -0.0009 ===================================== 输出修正量： △e △f -0.010141 -0.003138 -0.000547 -0.003197 ===================================== ** 修正： e f 1.042573 -0.100618 1.019453 -0.033553 ================================== Q:修正后电压的极坐标形式： newCode oldCode U ∠ δ 1 1 1.047417∠-0.096211 2 2 1.020005∠-0.032901 ----------- 第 3 次迭代 ----------- 输出不平衡量： PQ节点： ΔP ΔQ -0.000012 -0.000762 PV节点： ΔP ΔU2 -0.000143 -0.000011 第二象限雅可比矩阵: 6.9746 1.8243 -2.8481 7.6445 第一象限雅可比矩阵: -4.9156 -1.4862 1.4862 -4.9156 第三象限雅可比矩阵: -4.8690 -1.1416 0.0000 0.0000 第四象限雅可比矩阵: 5.0446 1.4418 -0.0671 2.0389 **雅可比增广矩阵: 6.9746 1.8243 -4.9156 -1.4862 -0.0000 -2.8481 7.6445 1.4862 -4.9156 -0.0008 -4.8690 -1.1416 5.0446 1.4418 -0.0001 0.0000 0.0000 -0.0671 2.0389 -0.0000 ===================================== 输出修正量： △e △f -0.000102 -0.000036 -0.000008 -0.000083 ===================================== ** 修正： e f 1.042471 -0.100653 1.019445 -0.033637 ================================== Q:修正后电压的极坐标形式： newCode oldCode U ∠ δ 1 1 1.047319∠-0.096254 2 2 1.020000∠-0.032983 ----------- 第 4 次迭代 ----------- 输出不平衡量： PQ节点： ΔP ΔQ 0.000000 -0.000000 PV节点： ΔP ΔU2 -0.00000