电力系统分析基础(第四章).ppt

第三节高斯—塞德尔迭代法潮流计算三、高斯-塞德尔迭代法潮流计算4.潮流计算：平衡节点的功率：支路功率：支路功率损耗：第29页，共75页，星期日，2025年，2月5日第三节高斯—塞德尔迭代法潮流计算四、高斯-塞德尔迭代法潮流计算流程图见书上P155第30页，共75页，星期日，2025年，2月5日第三节高斯—塞德尔迭代法潮流计算四、例题：用G-S计算潮流分布解：网络的节点导纳距阵为：~~1231.17-j4.71y135.88-j23.5j0.33y12y30平衡节点U1=1.00°PQ节点S2=-0.8-j0.6PU节点P3=0.4,U3=1.1第31页，共75页，星期日，2025年，2月5日第三节高斯—塞德尔迭代法潮流计算设，代入式（6-24）求第32页，共75页，星期日，2025年，2月5日第三节高斯—塞德尔迭代法潮流计算修正U3为，再用式（6-25）计算：然后开始第二次迭代：第33页，共75页，星期日，2025年，2月5日第三节高斯—塞德尔迭代法潮流计算再修正U3为：因此，第二次迭代结束时节点2的电压为节点3的电压相位角为δ3=2.940o,与之对应的节点3的无功功率为Q3=0.0596.再计算第34页，共75页，星期日，2025年，2月5日第四节牛顿—拉夫逊法潮流计算一、N-R原理1.非线性方程的求解：f(x)=0 设：x(0)为的初始近似解，?x(0)为与真实解的偏差则：x=x(0)—?x(0) f(x(0)—?x(0))=0 按Taylor’s展开f(x(0)—?x(0))= f(x(0))-f(x(0))?x(0)+…..+(-1)nfn(x(0))(?x(0))n/n!+….=0第35页，共75页，星期日，2025年，2月5日第四节牛顿—拉夫逊法潮流计算由于?x(0）较小，故忽略高次项后：f(x(0)—?x(0))=f(x(0))-f(x(0))?x(0)=0?x(0)=f(x(0))/f(x(0)) x(1)=x(0)-?x(0)=x(0)-f(x(0))/f(x(0))k次迭代时修正方程为：f(x(k))-f(x(k))?x(k)=0 ?x(k)=f(x(k))/f(x(k)) x(k+1)=x(k)-f(x(k))/f(x(k)) 第36页，共75页，星期日，2025年，2月5日第四节牛顿—拉夫逊法潮流计算结束迭代的条件（收敛）：|f(x(k))|ε1或|?x(k)|ε2物理意义初值不当不收敛第37页，共75页，星期日，2025年，2月5日第四节牛顿—拉夫逊法潮流计算2.非线性方程组的求解：推广于（6-16）表示的多变量非线性方程组(4-31)第38页，共75页，星期日，2025年，2月5日第四节牛顿—拉夫逊法潮流计算(4-32)式中为函数fi(x1,x2,….xn)对自变量xj的偏导在初始值处的值第39页，共75页，星期日，2025年，2月5日第四节牛顿—拉夫逊法潮流计算(4-33)用矩阵表示得到新的近似解：第40页，共75页，星期日，2025年，2月5日第四节牛顿—拉夫逊法潮流计算(4-33a)更一般的表示第k+1次迭代后的解为：(4-33b)第41页，共75页，星期日，2025年，2月5日第四节牛顿—拉夫逊法潮流计算式(4-33a)可简写为：F(x(k))=J(k)?x(k)(4-34)J(k)为n?n阶雅可比矩阵,其元素为函数fi(x1,x2,….xn)对自变量xj的偏导在点（x(k)）的值式(4-33b)可简写为：x(k+1)=x(k)—?x(k)(4-34b)第k次迭代后用下面的公式检查是否收敛 (4-35a)(4-35b)第42页，共75页，星期日，2025年，2月5日第四节牛顿—拉夫逊法潮流计算二、潮流计算时的修正方程式（直角座标法）(4-38a)(4-38b)1.直角座标法：PQ节点PU节点第43页，共75页，星期日，2025年，2月5日