课程设计指导方案计算潮流计算程序设计1.doc

电力系统分析课程设计讲义 一、 课程设计目的： 其目的在于加深学生对电力系统稳态分析中课程中基本概念和计算方法的理解，培 养学生运用所学知识分析和解决问题的能力。具体目标如下： 1）理解和掌握运用计算机进行潮流计算的基本算法原理； 2）针对某一具体电网，能够建立网络的计算模型，并可根据该等值电路模型进行 潮流计算程序设计。 3）学习并掌握 MATLAB 语言编程方法； 4）利用 MATLAB 语言实现潮流计算程序的编译，并对运行结果进行验证分析。 二、设计主要内容要求： 1）学习并掌握 MATLAB 语言。 2）掌握变压器非标准变比概念及非标准变比变压器的等值电路。掌握节点导纳矩 阵的概念及导纳矩阵的形成和修改方法。 3）掌握电力系统功率方程、变量和节点分类。 4）选择一种计算方法，设计潮流计算的方法和步骤。 5）选择一个某一具体电网，编制程序流程框图。 6）利用 MATLAB 语言???制潮流计算的程序，并上机调试程序，对计算结果进行 验证分析。 7）整理课程设计论文。 三、课程设计相关算法原理 1 节点导纳矩阵的形成和修改 1.1 节点导纳矩阵的形成 节点导纳矩阵的形成既可根据自导纳和互导纳的定义直接求取,也可运用支路节点 关联矩阵计算。 根据定义直接求取节点导纳矩阵时,需注意以下几点: 1) 节点导纳矩阵是方阵,其阶数就等于网络中初参考节点外的节点数 n。 2) 节点导纳矩阵是稀疏矩阵，其各行非零非对角元素等于该行相对应节点所连接 的不接地支路数。 3) 节点导纳矩阵的对角元就等于连接节点支路导纳的和。 4) 节点导纳矩阵的非对角元Yij 等于连接节点 i , j 支路导纳的负值。 5) 节点导纳矩阵一般是对称矩阵，这是网络的互易特性所决定的。从而，一般只 要求取这个矩阵的上三角或下三角部分就可以了。 6) 网络中的变压器，仍可按上述原则进行计算。 1.2 节点导纳矩阵的修改 在电力系统潮流计算中,往往要计算不同结线方式下的运行状况,例如,某电力线路 或变压器投入前后的状况以及某些元件参数变更前后的运行状况.由于改变一个支路的 参数或它的投入、退出状态只影响该支路两端节点的自导纳和它们之间的互导纳,可不 必重新形成与新运行状况相对应的节点导纳矩阵,仅需对原有的矩阵作某些修改以下介 绍几种典型的修改方法： 1) 从原有网络引出一支路,同时增加一节点,如图 1(a) 所示。 图 1 设 i 为原有节点，j 为新增加节点，新增加支路导纳 yij 。则因新增加一节点，节点 导纳矩阵将增加一阶。新增加的对角元 Y jj ，由于在节点 j 只有一个支路 yij ，将为 Y jj = yij 。 则因新增加一节点，节点导纳矩阵将增加一阶。 新增的对角元 Y jj ，由于在节点 j 只有一个支路 yij ，将为 Y jj = yij ；新增的非对角元 Yij = Y ji 则为 Yij = Y ji =- yij ；原有矩阵中的对角元Yii 将增加 ΔYii = yij 。 2) 在原有网络的节点 i、j 之间增加一支路,如图 1(b)所示。 这时由于仅增加支路不增加节点，节点导纳矩阵阶数不变，但与节点 i、j 有关元素 应作如下修改 ΔYii = yij ； ΔY jj = yij ； ΔYij = ΔY ji = ? yij 3) 在原有网络的节点 i、j 之间切除一支路,如图 1(c)所示。 切除一导纳为 yij 的支路相当于增加一导纳为- yij 的支路，从而与节点 i、j 有关元素 应作如下修改：  ΔYii = ? yij ；  ΔY jj  = ? yij ；  ΔYij  = ΔY ji = yij  ij ij ij ij 4) 原有网络节点 i、j 之间的导纳由 y 改变为 y ,如图 1(d)所示。这种情况相当于切 除一导纳为 y 的支路并增加一导纳为 y 的支路, 从而与节点 i、j 有关元素应作如下修 改： ΔYii = yij ? yij ； ΔY jj = yij ? yij ； ΔYij = ΔY ji = yij ? yij ；  2 迭代法求解非线性方程 2.1 高斯--塞德尔迭代法 高斯--塞德尔迭代法既可用以解线性方程组,也可用以解非线性方程组,其标准形 式如下。 设有方程组 a11 x1 + a12 x2 + a13 x3 = y1 ?  a21 x1 + a22 x2 + a23 x3 ? = y2 ? ?  (2.1)  它可改写为 a31 x1 + a32 x2 + a33 x3 = y3 ?  x1 = 1 a11  ( y1 a12 x2 ? a13 x3