电力系统分析课程设计——牛顿拉夫逊潮流计算

电力系统编程潮流计算1 设计任务及初步分析1。1 设计任务条件：节点数:3 支路数:3 计算精度:0。00010支路 1: 0.0300+j0.0900 1┠—————□——-——┨2支路 2: 0。0200+j0。0900 2┠———--□-————┨3支路 3: 0。0300+j0。0900 3┠——--—□—-——-┨1节点 1:PQ 节点,S（1）=—0.5000—j0.2000节点 2：PQ 节点,S（2）=-0.6000—j0。2500节点 3：平衡节点，U（3）=1。0000∠0。0000要求:编写程序计算潮流1。2 初步分析潮流计算在数学上可归结为求解非线性方程组，其数学模型简写如下： 2 牛顿-拉夫逊法简介2。1 概述牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程,通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解.因为越靠近解,偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性.而所谓“某一邻域"是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压（相角为0,幅值为1)启动即在此邻域内.2．2 一般概念对于非线性代数方程组即 （2－1)在待求量的某一个初始计算值附件,将上式展开泰勒级数并略去二阶及以上的高阶项，得到如下的线性化的方程组 （2－2)上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 （2－3)将和相加，得到变量的第一次改进值。接着再从出发，重复上述计算过程。因此从一定的初值出发，应用牛顿法求解的迭代格式为 (2－4) (2－5）上两式中:是函数对于变量的一阶偏导数矩阵，即雅可比矩阵；为迭代次数.由式（2－4)和式子（2－5）可见，牛顿法的核心便是反复形成求解修正方程式。牛顿法当初始估量值和方程的精确解足够接近时，收敛速度非常快 ,具有平方收敛特性.2。3 潮流计算的修正方程运用牛顿－拉夫逊法计算潮流分布时,首先要找出描述电力系统的非线性方程.这里仍从节点电压方程入手，设电力系统导纳矩阵已知，则系统中某节点（节点)电压方程为从而得 进而有 （2－6)式（2－6）中，左边第一项为给定的节点注入功率，第二项为由节点电压求得的节点注入功率。他们二者之差...