潮流计算程序设计doc.docVIP

潮流计算程序设计 指导老师：李咸善 李欣 学生姓名：翁雪莹 学生学号：2010117141 课题内容：用Matlab语言或c语言或其他计算机语言，采用牛顿-拉夫逊法或PQ分解法，设计电力系统潮流计算程序。选择某一个题目(课本上的例题或作业题，或自己通过其他途径选择一个潮流计算的案例)，进行潮流计算相关的模型建立，包括节点导纳矩阵、功率方程，雅科比矩阵的推导；设计计算机程序，进行调试及结果分析。 潮流计算原理简介：所谓电力系统的潮流，是指系统中所有运行参数的总体。电力系统潮流计算是电力系统运行和规划中最基本和最经济的计算，其任务是要在已知（或给定）某些运行参数的情况下，计算出系统中全部运行参数，确定系统的运行状态，如各母线上的电压（幅值和相角）、网络中的功率分布及功率损耗等，是电力系统的稳态计算。 潮流计算可以用传统的手工方式进行，也可以计算机为工具通过软件完成。后者也有不同的方法来编程，我们学习的牛顿-拉夫逊法和PQ分解法就是其中的两种。 题目：网络接线如图，各支路阻抗和各节点功率均以标幺值标于图中，其中节点2连接的实际是发额定功率的发电厂。设节点1的电压保持为1.06。 基于牛顿—拉夫逊算法的计算机求解步骤： 第一步：构建节点导纳矩阵 本题有5个节点，导纳矩阵为阶，题目给出的是线路中的阻抗，首先由把线路阻抗值转换成线路导纳值。 对于对角线元素的导纳值为该行所有导纳相加，其他位置上的值是线路导纳的相反数。用计算机编程实现可以得到本题的节点导纳矩阵结果如下： 第二步：给出电压初值选取平衡节点（程序编制选取第一行为平衡节点） 本题5个节点的电压值依次为 [1.06,1,1,1,1] 第三步：利用设定的电压初值和节点导纳矩阵，求取剩下节点的功率。 由于我编程使用的是MATLAB6.5可以很方便的进行复数的计算，故功率的计算我采用的是原始的功率计算公式 第四步：构造方程中的 的行数是剩下节点数目的两倍。 本题的为 编制程序可扩展为 其中； 的行数等于的行数 本题的为 编制程序可扩展为 对于本题的则为阶，我把他分解为4个阶依次为 ， ， 鉴于MATLAB强大的矩阵运算能力使用就可得到 。 编制程序时、、、可以随着网络节点数目任意扩展形式如下： 由 其中当时 ； ； 当时 首先由，可得 ； ； 第五步 利用MATLAB inv（）命令求出 第七步 修正电压值 第六步：判断是否小于，若小于停止迭代，计算平衡节点功率，各条线路的损耗，否则回到第三步进行下一次迭代。 程序运行的结果如下 迭代次数4次 各个节点功率 （平衡节点） 0.0677 + 0.1531i 0.2000 + 0.2000i -0.4500 - 0.1500i -0.4000 - 0.0500i 0.6000 - 0.1000i 各个节点电压 （平衡节点） 1.0600 1.0572 + 0.0058i 1.0316 - 0.0269i 1.0328 - 0.0244i 1.0546 + 0.0440i 各条线路的损耗 行 1 2 3 4 5 1 0 0.0771 - 0.0758i -0.1448 - 0.0773i 0 0 2 -0.0773 + 0.0752i 0 -0.2177 - 0.0787i -0.2020 - 0.0768i 0.2970 - 0.1196i 3 0.1428 + 0.0715i 0.2148+ 0.0701i 0 0.0923 + 0.0084i 0 4 0 0.1995+ 0.0693i -0.0924 - 0.0086i 0 0.2929 - 0.0107i 5 0 -0.3006+0.1086i 0 -0.2994 - 0.0086i 0 对程序的改进 电力系统中不仅有PQ节点，还有PV节点，程序能否囊括电力系统中各种节点的潮流计算呢？答案是肯定的。 PV节点只是不知道Q，但是却给出了V，我就把对、的微分换成对、的微分；由；当时，；，当时； ；并把相应的中的换成，不变； 改动本题中的相关条件进行验证，本网络的潮流分布不变。 第一种程序检验方法 依次取5，4，3，2为PV节点，电压为最后节点电压的模值，得到的潮流分布几乎与上面的结论一样，同时取2，3，4，5同时为PV节点，电压为最后节点电压的模值，得到的潮流分布几乎与上面的结论一样。 第二种程序检验方法 利用《电力系统分析》（第三版），于永源，电力系统潮流计算的计算机算法第四章的4—13题最后潮流分布如下： 迭代次数5次 各个节点功率 （平衡节点） 3.0000 + 0.8502i -1.0000+ 0.7133i -2.0000 - 1.0000i 各个节点电压 （平衡节点） 1.000