电力系统潮流的计算机算法.pptVIP

极座标法矩阵表示二、极坐标形式的N-R潮流算法第29页，共59页，星期日，2025年，2月5日极座标法系数推导展开式计及第30页，共59页，星期日，2025年，2月5日极座标法系数推导当i≠j，对特定的j，只有特定节点的δj，从而δij=δi-δj是变量对特定的j，只有该特定节点的Uj是变量第31页，共59页，星期日，2025年，2月5日极座标法系数推导当i=j，由于δi是变量，从而所有δij=δi-δj都是变量，可得相似地，由于Ui是变量，可得第32页，共59页，星期日，2025年，2月5日极座标法潮流计算的基本步骤：1.输入原始数据和信息：y、C、Pis、Qis、Uis、约束条件2.形成节点导纳矩阵YB=CTyC设置各节点电压初值Ui(0),δi(0)将初始值代入(4-48)求不平衡量?Pi(0),?Qi(0)5.计算雅可比矩阵各元素（Hij、Lij、Nij、Jij）6.解修正方程(4-53)，求?Ui(k),?δi(k)7.求节点电压新值Ui(k+1)=Ui(k)-?Ui(k),δi(k+1)=δi(k)-?δi(k)第33页，共59页，星期日，2025年，2月5日极座标法潮流计算的基本步骤：8.判断是否收敛：Max|?Ui(k)|≤ε,Max|?δi(k)|≤ε9.重复迭代第4、5、6、7步，直到满足第8步的条件求平衡节点的功率和PV节点的Qi及各支路的功率第34页，共59页，星期日，2025年，2月5日三、直角坐标形式的N-R潮流算法1.直角座标法：PQ节点PU节点第35页，共59页，星期日，2025年，2月5日直角坐标的缩写形式：第36页，共59页，星期日，2025年，2月5日直角座标法矩阵表示第37页，共59页，星期日，2025年，2月5日雅可比矩阵元素值非对角元素(i≠j)对角元素(i=j)第38页，共59页，星期日，2025年，2月5日雅可比矩阵的特点：各元素是各节点电压的函数不是对称矩阵∵Yij=0，∴Hij=Nij=Jij=Lij=0，另Rij=Sij=0，故稀疏直角坐标潮流计算的基本步骤：1.输入原始数据和信息：y、C、Pis、Qis、Uis、约束条件2.形成节点导纳矩阵YB=CTyC设置各节点电压初值ei(0),fi(0)4.将初始值代入(4-38)求不平衡量?Pi(0),?Qi(0),?Ui2(0)第39页，共59页，星期日，2025年，2月5日5.计算雅可比矩阵各元素（Hij、Lij、Nij、Jij、Rij、Sij）6.解修正方程(4-37)，求?ei(k),?fi(k)7.求节点电压新值ei(k+1)=ei(k)-?ei(k),fi(k+1)=fi(k)-?fi(k)8.判断是否收敛：Max|?fi(k)|≤ε,Max|?ei(k)|≤ε9.重复迭代第4、5、6、7步，直到满足第8步的条件求平衡节点的功率和PV节点的Qi及各支路的功率返回第40页，共59页，星期日，2025年，2月5日第四节N-R潮流计算中的收敛性和稀疏技术一、稀疏矩阵的存储151081154207149123222171316181、按坐标存储的方案顺序号L123456DIAG151120122218对角元素第41页，共59页，星期日，2025年，2月5日15108115420714912322217131618非对角元素顺序号m123456789101112OFFD10485714932171316IROW112333445566ICOL231146354635特点：按坐标位置存储，简单、直观，便于检索，但不便于运算。需n+3N个存储单元（n为对角元数，N为非零非对角元数）。第42页，共59页，星期日，2025年，2月5日第1页，共59页，星期日，2025年，2月5日第四章复杂电力系统潮流的计算机算法第二节潮流计算的节点功率方程和节点分类第三节潮流计算的牛顿-拉夫逊法第一节网络方程式第四节牛顿-拉夫逊潮流计算中的收敛性和稀疏技术第五节其他潮流计算方法简介第2页，共59页，星期日，2025年，2月5日第一节网络方程式

网络方程2用节点导纳矩阵表示的网络方程