5.-电力系统分析.pptVIP

⑶ 请输入各节点参数形成的矩阵B2 矩阵B2的每行是由下列参数构成的： 节点所接发电机的功率SG； 节点负荷的功率SL； 节点电压的初始值。 PU节点电压U的给定值。 节点所接的无功补偿设备的容量。 节点分类标号igl。 1 －－平衡节点 igl= 2－－PQ节点 3－－PU节点 （4）请输入由节点号及其容抗形成的矩阵X 3.4实例 2.先形成节点导纳矩阵。 3.根据式（19.15）和（19.16）求出修正方程式的常数项向量 4.据式（19.18）和（19.19）求出雅可比矩阵各元素值，即可得到第一次迭代时的修正方程式。 5.解方程式（19.16），求 6.修正各节点电压，即得出第一次迭代后各节点的电压值。 7.按以上计算步骤迭代下去，当收敛精度取 （即PR=0.0001）时，需要进行四次迭代。求出了各节点电压后，即可求各支路的潮流分布。 3.4实例 针对例19.1输入数据如下： 请输入节点数:n=5 请输入支路数:nl=5 请输入平衡母线节点号:isb=1 请输入误差精度:pr=0.00001 请输入由支路参数形成的矩阵:B1=[1 2 0.03i 0 1.05 0; 2 3 0.08+0.3i 0.5i 1 0;2 4 0.1+0.35i 0 1 0; 3 4 0.04+0.25i 0.5i 1 0;3 5 0.015i 0 1.05 1] 请输入各节点参数形成的矩阵:B2=[0 0 1.05 1.05 0 1; 0 3.7+1.3i 1 0 0 2;0 2+1i 1 0 0 2;0 1.6+0.8i 1 0 0 2; 5 0 1.05 1.05 0 3] 请输入由节点号及其对地阻抗形成的矩阵:X=[1 0;2 0;3 0; 4 0;5 0] 3.4实例 结果如下： 迭代次数 5 没有达到精度要求的个数 7 8 8 6 0 各节点的实际电压标么值E为(节点号从小到大排列): 1.0500 1.0335 - 0.0774i 1.0260 + 0.3305i 0.8592 - 0.0718i 0.9746 + 0.3907i 各节点的电压大小U为(节点号从小到大排列): 1.0500 1.0364 1.0779 0.8622 1.0500 各节点的电压相角为θ(节点号从小到大排列): 0 -4.2819 17.8535 -4.7785 21.8433 各节点的功率S为(节点号从小到大排列): 2.5794 + 2.2994i -3.7000 - 1.3000i -2.0000 - 1.0000i -1.6000 - 0.8000i 5.0000 + 1.8131i 3.4实例 各条支路的首端功率Si为(顺序同您输入B1时一样): 2.5794 + 2.2994i -1.2774 + 0.2032i 0.1568 + 0.4713i 1.5845 + 0.6726i 5.0000 + 1.8131i 各条支路的末端功率Sj为(顺序同您输入B1时一样): -2.5794 - 1.9745i 1.4155 - 0.2443i -0.1338 - 0.3909i -1.4662 - 0.4091i -5.0000 - 1.4282i 各条支路的功率损耗DS为(顺序同您输入B1时一样): -0.0000 + 0.3249i 0.1381 - 0.0412i 0.0230 + 0.0804i 0.1184 + 0.2635i 0.0000 + 0.3849i 3.4实例 每次迭代后各节点的电压值如图 3.4实例 4 PQ分解法潮流计算 4.1 PQ分解法的基本方程式 4.2 计算步骤和程序框图 4.3 实例 4.1 PQ分解法的基本方程式 （19.24） （19.25） 快速分解法（又称P—Q分解法）是从简化牛顿法极坐标形式计算潮流程序的基础上提出来的。 节点功率方程在用极坐标形式表示时，它的修正方程式可简化为： 在一般情况下，线路两端电压的相角是不大的（不超过10°～20 ° ）。因此， 与系统各节点无功功率相应的导纳B 远远小于该节点自导纳的虚部，即 因而 （19.26） 式（19.24）的系数矩阵中的各元素可表示为： （i,j=1,2,……，n-1） （19.27） （i,j=1,2,……，m） (19.28) 4.1 PQ分解法的基本方程式 而系数矩阵H和