题目19-根据MATLAB的电力系统复杂潮流计算.doc

/. 编号 课 程 设 计 （ 2013级本科） 题 目： 基于MATLAB的复杂潮流计算 系 院： 物理与机电工程学院 专 业： 电气工程及其自动化 作者姓名： 指导教师： 田 娜 职称： 助 教 完成日期： 2016 年 6 月 30 日 二○一六 年 六 月 河西学院本科生课程设计任务书 设 计 题 目 基于MATLAB的复杂潮流计算 作 者 姓 名 学院、专业、年级 物电学院电气工程及其自动化专业13级 指导教师姓名、职称 田娜 助教 任务下达日期 2016年6月1日 1.设计内容 原始参考图 其中节点1为平衡节点，节点2、3、4、5为PQ节点。根据原始参考图编写MATLAB程序进行潮流分析，计算出正确的潮流结果。 2.设计的基本要求 2.1设计及计算说明书 （1）说明书要求书写整齐，条理分明，表达正确、语言正确。 （2）计算书内容：为各设计内容最终成果、确定提供依据进行的技术分析、论证和定量计算，如。 （3）计算书要求：计算无误，分析论证过程简单明了，各设计内容列表汇总。 2.2图纸 （1）绘制分析所需的必要图纸 （2）图纸要求：用标准符号绘制，布置均匀，设备符号大小合适，清晰美观。 (3)说明书后应附录MATLAB程序。 3.论文(设计)进度安排 阶段 论文（设计）各阶段名称 起止日期 1 熟悉设计任务书、设计题目及设计背景资料 2 查阅有关资料 3 阅读设计要求必读的参考资料 4 书写设计说明书 5 上交设计成果 4.需收集和阅读的资料及参考文献（指导教师指定） [1]: 陈珩.电力系统稳态分析（第三版）[M]，北京，中国电力出版社，2007 [2]：何仰赞，温增银.《电力系统分析》第三版[M]，武汉,华中科技大学出版社，2002 [3]：陈悦.《电气工程毕业设计指南电力系统分册》[M]，北京，中国水利水电出版社，2008 教 研 室 意 见 负责人签名： 年 月 日 学 院 意 见 负责人签名： 年 月 日 目录 摘要5 第一章 电力系统潮流计算概述6 1.1电力系统概述6 1.2 电力系统潮流概述7 1.3 潮流计算的目的8 1.4电力系统的发展和分析计算9 1.5、MATLAB软件的应用10 第二章 牛顿—拉夫逊法潮流计算基本原理11 2.1牛顿—拉夫逊法潮流计算简介11 2.2牛顿—拉夫逊法潮流计算计算公式11 2.3牛顿—拉夫逊法解题的一般步骤14 第三章 网络潮流计算15 3.1电力系统设计图15 3.2网络潮流计算的手工算法15 3.3牛拉法潮流计算的流程图17 3.4 MATLAB算法的计算程序18 3.5 MATLAB的计算结果23 总结及感想37 参考文献及资料37 摘 要 潮流计算，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。 它是基于配电网络特有的层次结构特性，论文提出了一种新颖的分层前推回代算法。该算法将网络支路按层次进行分类，并分层并行计算各层次的支路功率损耗和电压损耗，因而可大幅度提高配电网潮流的计算速度。论文在MATLAB环境下，利用其快速的复数矩阵运算功能，实现了文中所提的分层前推回代算法，并取得了非常明显的速度效益。另外，论文还讨论发现，当变压器支路阻抗过小时，利用Π型模型会产生数值巨大的对地导纳，由此会导致潮流不收敛。为此，论文根据理想变压器对功率和电压的变换原理，提出了一种有效的电压变换模型来处理变压器支路，从而改善了潮流算法的收敛特性。 关键词：电力系统；潮流分析；MATLAB 第一章 电力系统潮流计算概述 1.1电力系统概述 1831年法拉第发现了电磁感应定律。在此基础上，很快出现了原始的交流发电机、直流发电机和直流电动机。由于当时发电机发出的电能仅用于电化学工业和电弧灯，而电动机所需的电能又来自蓄电池，电机制造和电力输送技术的发展最初集中于直流电。原始的 电力线路使用的就是100～400V低压直流电。由于输电电压低，输送的距离不可能远，输送的功率也不可能大。与100余年前电力系统的雏形相比，近代电力系统不仅在输电电压、输送距离、输送功率等方面有了千百倍的增长，而且电源构成、负荷成分等方面也有了很大变化。 电力系统加上发电机的原动机(如汽轮机、水轮机)，原动机的力能部分(如热力锅炉、水库、原子能电站的反应堆)、供热和用热设备，则称为动力系统。 电力工业发展初期，电能是直接在用户附近的发电站(或称发电厂)中生产的，各发电站孤立运行。随着工农业生产和城市的发展，电能的需要量迅速增加，而热能资源(如煤田)和水能资源丰富的地区又往往远离用电比较集中的城市和工矿区，为了解决这个矛盾，就需要在动力资源丰富的地区建立大型发电站，然后将电能远距离输送给电力用户。同时，为了提高供电可靠性以及资源利用的综合经济性，又把许多分散的各种形式的发电站，通过送电线路和变电所联系起来。这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接起来的整体，即称为电力系统。 现代电力系统提出了“灵活交流输电与新型直流输电”的概念。灵活交流输电技术是指运用固态电子器件与现代自动控制技术对交流电网的电压、相位角、阻抗、功率以及电路的通断进行实时闭环控制，从而提高高压输电线路的输送能力和电力系统的稳定水平。新型直流输电技术是指应用现电力电子技术的最新成果，改善和简化变流站的造价等。 运行方式管理中，潮流是确定电网运行方式的基本出发点；在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了电网在预想操作情况下电网的潮流分布以校验运行可靠性。在电力系统调度运行的多个领域都涉及到电网潮流计算。潮流是确定电力网络运行状态的基本因素，潮流问题是研究电力系统稳态问题的基础和前经验主义的研究方法 1.2 电力系统潮流概述 电力系统在运行时，在电源电势激励作用下，电流或功率从电源通过系统各元件流入负荷，分布于电力网各处，称为潮流分布。 什么是电力系统的潮流？ 什么是潮流方向？ 指电力系统上网.或是输电的发展方向.用电力系统上网不会成为注流.用远程在功率输电,肯定的发展方向是直流输电. 潮流计算是电力系统分析中的一种最基本的计算，它的任务是在给定的接线方式和运行条件下，确定系统的运行状态，如各母线上的电压（幅值和相角）、网络中的功率分布及功率损耗等，是电力系统的稳态计算。潮流计算是对电力系统正常运行状况的分析和计算，即电力系统中的电压、电流、功率的计算，即潮流计算；潮流计算方法很多：高斯—塞德尔法、牛顿—拉夫逊法、P-Q分解法、直流潮流法，以及由高斯—塞德尔法、牛顿—拉夫逊法演变的各种潮流计算方法。 潮流计算可以用传统的手工方式进行，也可以计算机为工具通过软件完成。两种方法各有优缺点。前者物理概念清晰，可用来计算一些接线较简单的电力网，但若将其用于接线复杂的电力网则计算量过大，难于保证计算准确性。后者从数学上看可归结为用数值方法解非线性代数方程，数学逻辑简单完整，借助计算机可快速精确地完成计算，但其缺点是物理概念不明显，物理规律被埋没在循环往复的数值求解过程中。 潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求：对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各元件（线路、变压器等）是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。因此潮流计算的目的是： ① 为电力系统规划设计提供接线、电气设备选择和导线截面选择的依据。 ② 提供电力线运行方式和制定检修计划的依据。 ③ 提供继电保护、自动装置设计和整定计算的依据。 ④ 为调压计算、经济运行计算、短路和稳定计算提供必要的数据。 1.3 潮流计算的目的 电力系统的潮流计算最主要的目的是为了让电力系统能够安全稳定运行的同时做到经济运行。所以考留到经及调度、电网规划、电力系统可靠性分析。 具体表现在以下方面： ①在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。 ②在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。 ③正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。 ④预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。 总结为在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。 1.4电力系统的发展和分析计算 潮流计算针对电力系统各种正常运行方式，而静态安全分析则要研究各种运行方式下个别系统元件退出运行后系统的状况。其目的是校验系统是否能安全运行，即是否有过负荷的元件或电压过低的母线等。原则上讲，静态安全分析也可用潮流计算来代替。但是一般静态安全分析需要校验的状态数非常多，用严格的潮流计算来分析这些状态往往计算量过大，因此不得不寻求一些特殊的算法以满足要求。 利用电子数字计算机进行电力系统潮流计算从20 世纪50 年代中期就己开始，此后，潮流计算曾采用了各种不同的方法，这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的，对潮流计算的要求可以归纳为下面几点： ①计算方法的可靠性或收敛性 ②对计算速度和内存量的要求 ③计算的方便性和灵活性 现在主要的研究都是围绕着牛顿拉夫逊法来进行的，牛顿法是数学中解决非线性方程式的典型方法，有较好的收敛性。解决电力系统潮流计算问题是以导纳矩阵为基础的，因此，只要在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的放率。自从20 世纪60 年代中期利用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、速度方面都超过了阻抗法，成为直到目前仍在广泛采用的优秀方法。同样优秀的方法还有很多，fast decoupled, PQ分解法等等。 1.5、MATLAB软件的应用 MATLAB Compiler是一种编译工具，它能够将M编写的函数文件生成函数库或者可执行文件COM组件等，以提供给其他高级语言如C++、C#等进行调用由此扩展MATLAB的应用范围，将MATLAB的开发效率与其他高级语言的运行结合起来，取长补短，丰富程序开发的手段。 目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（1）计算速度快（2）内存需要少（3）计算结果有良好的可靠性和可信性（4）适应性好，即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强（5）简单。 MATLAB是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。 MATLAB程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。MATLAB与C语言和FORTRAN语言相比更容易被掌握。通过M语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上。 另外，MATLAB提供了一种特殊的工具：工具箱（TOOLBOXES）.这些工具箱主要包括：信号处理（SIGNAL PROCESSING）、控制系统（CONTROL SYSTEMS）、神经网络（NEURAL NETWORKS）、模糊逻辑(FUZZY LOGIC)、小波(WAVELETS)和模拟（SIMULATION）等等。不同领域、不同层次的用户通过相应工具的学习和应用，可以方便地进行计算、分析及设计工作。 MATLAB设计中，原始数据的填写格式是很关键的一个环节，它与程序使用的方便性和灵活性有着直接的关系。原始数据输入格式的设计，主要应从使用的角度出发，原则是简单明了，便于修改。 第二章 牛顿—拉夫逊法潮流计算基本原理 2.1牛顿—拉夫逊法潮流计算简介 牛顿迭代法（Newtons method）又称为牛顿-拉夫逊方法，它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。多数方程不存在求根公式，因此求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要。方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。 设r是f(x) = 0的根，选取x0作为r初始近似值，过点（x0,f(x0)）做曲线y = f(x)的切线L，L的方程为y = f(x0) f(x0)(x-x0)，求出L与x轴交点的横坐标 x1 = x0-f(x0)/f(x0)，称x1为r的一次近似值。过点（x1,f(x1)）作曲线y = f(x)的切线，并求该切线与x轴的横坐标 x2 = x1-f(x1)/f(x1)，称x2为r的二次近似值。重复以上过程，得r的近似值序列，其中x(n+1)=x(n)－f(x(n))/f(x(n))，称为r的n+1次近似值，上式称为牛顿迭代公式。 解非线性方程f(x)=0的牛顿法是把非线性方程线性化的一种近似方法。把f(x)在x0点附近展开成泰勒级数 f(x) = f(x0)+(x－x0)f(x0)+(x－x0)^2*f(x0)/2! +… 取其线性部分，作为非线性方程f(x) = 0的近似方程，即泰勒展开的前两项，则有f(x0)+f(x0)(x－x0)=f(x)=0 设f(x0)≠0则其解为x1=x0－f(x0)/f(x0) 这样，得到牛顿法的一个迭代序列：x(n+1)=x(n)－f(x(n))/f(x(n))。 2.2牛顿——拉夫逊法潮流计算计算公式 把牛顿法用于潮流计算，采用直角坐标形式表示的如式（2-2）所示的形式。其中电压和支路导纳可表示为： （2-1） 将上述表示式（2-1）代入功率方程，展开并分出 实部和虚部，便得： （2-2） 按照以上的分类，PQ节点的输出有功功率和无功功率是给定的，则第i节点的给定功率设为和（称为注入功率）。 假定系统中的第1、2、…、m节点为PQ节点，对其中每一个节点的N-R法表达式 F(x)=0[如、、]形式有些下列方程： （2-3） =（1、2、…、m） PV节点的有功功率和节点电压幅值是给定的。假定系统中的第m+1、m+2、…、n-1节点为PV节点，则对其中每一PV节点可以列写方程： （2-4） =（m+1、m+2、…、n-1） 形成雅可比矩阵。对多维变量求偏导（、、、、、、、…），并以矩阵的形式表达称为雅可比矩阵。 当j=i时，对角元素为： （2-5） 当时,矩阵非对角元素为： （2-6） 由上式不难看出，雅可比矩阵有以下特点。 ① 雅可比矩阵中的诸元素都是节点电压的函数，因此在迭代过程中，它们将随着节点电压的变化而不断的变化。 ② 雅可比矩阵具有结构对称性，数据不对称。如非对角，，。 ③ 由式（2-6）可以看出，当导纳矩阵中非对角元素为零时，。雅可比矩阵中相应的元素也为零，即矩阵是非常稀疏的。因此，修正方程的求解同样可以应用稀疏矩阵的求解技巧。正是由于这一点才使N-R法获得广泛的应用。 2.3牛顿—拉夫逊法解题的一般步骤 以上讨论的是用直角坐标形式的牛顿—拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量由于平衡节点的电压向量是给定的，因此待求共需要2(n-1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。 （2-3-0） 对PQ节点来说，是给定的，因而可以写出 （2-3-1） 对PV节点来说，给定量是，因此可以列出式（2-3-2） （2-3-2） 求解过程大致可以分为以下步骤： （1）形成节点导纳矩阵 （2）将各节点电压设初值U， （3）将节点初值代入相关求式，求出修正方程式的常数项向量 （4）将节点电压初值代入求式，求出雅可比矩阵元素 （5）求解修正方程，求修正向量 （6）求取节点电压的新值 （7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步 （8）计算支路功率分布，PV节点无功功率和平衡节点注入功率。 第三章 复杂网络潮流计算 3.1 电力系统设计图 系统接线图 （其中节点1为平衡节点，节点2、3、4、5为PQ节点。） 3.2复杂网络潮流计算的手工算法 解：依题意，可知其等值阻抗电路图为 节点1为平衡节点,U1=1.06+J0为一值,其它四个节点都是PQ节点给定的注入功率为:S2 =0.20+J0.20,S3=-0.45-J0.15,S4=-0.40-J0.05,S5=-0.60-J0.10. 由上图可得相应的节点导纳矩阵 Y= 6.250-j18.750-5.000+j15.000-1.250+j3.75000-5.000+j15.00010.834-j32.500-1.667+j5.000-1.667+j5.000-2.500+j7.500-1.250+j3.750-1.667+j5.00012.917-j38.750-10.000+j30.00000-1.667+j5.000-10.000+j30.00012.917-j38.750-1.250+j3.7500-2.500+j7.5000-1.250+j3.7503.750-j11.250 计算各节点功率的不平衡量： 取=1.06，=0；====1.0;= ===0, 根据式（2-2）计算各节点初始功率、得： =-0.300；=-0.900 =-0.0750；=-0.2250 =0.0；=0.0 =0.0；=0.0 根据式（2-3）得各节点功率的不平衡量为： =0.5000；=1.1000 =-0.3750；=0.0750 =-0.4000；=-0.0500 =-0.6000；=-0.1000 根据公式： 3.3牛拉法潮流计算的流程图 3.4 MATLAB算法的计算程序 %开始 clc clear disp(节点总数为：); N=5 disp(平衡节点为：); 1 disp(PQ节点为：); JD=[2,3,4,5] e=[1.06 1 1 1 1]; f=[0 0 0 0 0]; P1=0; Q1=0; P2=-0.2; Q2=-0.2; P3=0.45; Q3=0.15; P4=0.4; Q4=0.05; P5=0.6; Q5=0.1; G=[6.2500,-5.0000,-1.2500,0,0;-5.0000,10.8340,-1.6670,-1.6670,-2.5000;-1.2500,-1.6670,12.9170,-10.0000,0;0,-1.6670,-10.0000,12.9170,-1.2500;0,-2.5000,0,-1.2500,3.7500];%形成电导矩阵。 B=[-18.75,15.0000,3.7500,0,0;15.0000,-32.5000,5.0000,5.0000,7.5000;3.7500,5.0000,-38.7500,30.0000,0;0,5.0000,30.0000,-38.7500,3.7500;0,7.5000,0,3.7500,-11.2500];%形成电纳矩阵。 disp(节点电导矩阵G为:); disp(G) disp(节点电纳矩阵B为:); disp(B) k=0; for v=1:7 I=[0,0;0,0;0,0;0,0;0,0]; for n=1:5 I(1,1)=I(1,1)+G(1,n)*e(n)-B(1,n)*f(n); I(1,2)=I(1,2)+G(1,n)*f(n)+B(1,n)*e(n); end for n=1:5 I(2,1)=I(2,1)+G(2,n)*e(n)-B(2,n)*f(n); I(2,2)=I(2,2)+G(2,n)*f(n)+B(2,n)*e(n); end for n=1:5 I(3,1)=I(3,1)+G(3,n)*e(n)-B(3,n)*f(n); I(3,2)=I(3,2)+G(3,n)*f(n)+B(3,n)*e(n); end for n=1:5 I(4,1)=I(4,1)+G(4,n)*e(n)-B(4,n)*f(n); I(4,2)=I(4,2)+G(4,n)*f(n)+B(4,n)*e(n); end for n=1:5 I(5,1)=I(5,1)+G(5,n)*e(n)-B(5,n)*f(n); I(5,2)=I(5,2)+G(5,n)*f(n)+B(5,n)*e(n); end H=[]; N=[]; M=[]; L=[]; J=[]; P2=P2-e(2)*I(2,1)-f(2)*I(2,2); %有功功率的不平衡量 Q2=Q2-f(2)*I(2,1)+e(2)*I(2,2); %无功功率的不平衡量 P3=P3-e(3)*I(3,1)-f(3)*I(3,2); Q3=Q3-f(3)*I(3,1)+e(3)*I(3,2); P4=P4-e(4)*I(4,1)-f(4)*I(4,2); Q4=Q4-f(4)*I(4,1)+e(4)*I(4,2); P5=P5-e(5)*I(5,1)-f(5)*I(5,2); Q5=Q5-f(5)*I(5,1)+e(5)*I(5,2); for m=2:5 for n=2:5 if(m==n) H(m,m)=-B(m,m)*e(m)+G(m,m)*f(m)+I(m,2); N(m,m)=G(m,m)*e(m)+B(m,m)*f(m)+I(m,1); M(m,m)=-G(m,m)*e(m)-B(m,m)*f(m)+I(m,1); L(m,m)=-B(m,m)*e(m)+G(m,m)*f(m)-I(m,2); else H(m,n)=-B(m,n)*e(m)+G(m,n)*f(m); N(m,n)=G(m,n)*e(m)+B(m,n)*f(m); M(m,n)=-N(m,n); L(m,n)=H(m,n); end end end J=[H(2,2),N(2,2),H(2,3),N(2,3),H(2,4),N(2,4),H(2,5),N(2,5);M(2,2),L(2,2),M(2,3),L(2,3),M(2,4),L(2,4),M(2,5),L(2,5);H(3,2),N(3,2),H(3,3),N(3,3),H(3,4),N(3,4),H(3,5),N(3,5);M(3,2),L(3,2),M(3,3),L(3,3),M(3,4),L(3,4),M(3,5),L(3,5);H(4,2),N(4,2),H(4,3),N(4,3),H(4,4),N(4,4),H(4,5),N(4,5);M(4,2),L(4,2),M(4,3),L(4,3),M(4,4),L(4,4),M(4,5),L(4,5);H(5,2),N(5,2),H(5,3),N(5,3),H(5,4),N(5,4),H(5,5),N(5,5);M(5,2),L(5,2),M(5,3),L(5,3),M(5,4),L(5,4),M(5,5),L(5,5)]; disp(雅克比矩阵J:); disp(J); A=[]; C=[P2;Q2;P3;Q3;P4;Q4;P5;Q5] A=J\C;%解修正方程式 disp(第M次修正方程的解A:); disp(A); f(2)=f(2) +A(1,1); e(2)=e(2) +A(2,1); %计算新值 f(3)=f(3) +A(3,1); e(3)=e(3) +A(4,1); f(4)=f(4) +A(5,1); e(4)=e(4) +A(6,1); f(5)=f(5) +A(7,1); e(5)=e(5) +A(8,1); disp(各点的电压实部e(单位:V)为(节点号从小到大排列):); disp(e) disp(各点的电压虚部f单位:V)为(节点号从小到大排列):); disp(f); u=e+f*i; disp(节点电压的第C(k)次近似值:); disp(u); k=k+1; disp(迭代次数:); disp(k); end for m=1:5 I(m)=(G(1,m)+B(1,m)*i) *u(m); end disp(平衡节点的功率); S1=u(1)*sum(conj(I))%计算平衡节点的功率 for m=1:5 for n=1:5 S(m,n)=u(m)*(conj(u(m))-conj(u(n)))*conj(-(G(m,n)+B(m,n)*i));%计算各支路功率 end end disp(各支路功率);disp(S) %结束 3.5 MATLAB的计算结果 节点总数为： N = 5 平衡节点为： ans = 1 PQ节点为： JD = 2 3 4 5 节点电导矩阵G为: 6.2500 -5.0000 -1.2500 0 0 -5.0000 10.8340 -1.6670 -1.6670 -2.5000 -1.2500 -1.6670 12.9170 -10.0000 0 0 -1.6670 -10.0000 12.9170 -1.2500 0 -2.5000 0 -1.2500 3.7500 节点电纳矩阵B为: -18.7500 15.0000 3.7500 0 0 15.0000 -32.5000 5.0000 5.0000 7.5000 3.7500 5.0000 -38.7500 30.0000 0 0 5.0000 30.0000 -38.7500 3.7500 0 7.5000 0 3.7500 -11.2500 雅克比矩阵J: Columns 1 through 6 33.4000 10.5340 -5.0000 -1.6670 -5.0000 -1.6670 -11.1340 31.6000 1.6670 -5.0000 1.6670 -5.0000 -5.0000 -1.6670 38.9750 12.8420 -30.0000 -10.0000 1.6670 -5.0000 -12.9920 38.5250 10.0000 -30.0000 -5.0000 -1.6670 -30.0000 -10.0000 38.7500 12.9170 1.6670 -5.0000 10.0000 -30.0000 -12.9170 38.7500 -7.5000 -2.5000 0 0 -3.7500 -1.2500 2.5000 -7.5000 0 0 1.2500 -3.7500 Columns 7 through 8 -7.5000 -2.5000 2.5000 -7.5000 0 0 0 0 -3.7500 -1.2500 1.2500 -3.7500 11.2500 3.7500 -3.7500 11.2500 C = 0.1000 0.7000 0.5250 0.3750 0.4000 0.0500 0.6000 0.1000 第M次修正方程的解A: 0.0473 0.0847 0.0863 0.1123 0.0922 0.1136 0.1076 0.1183 各点的电压实部e(单位:V)为(节点号从小到大排列): 1.0600 1.0847 1.1123 1.1136 1.1183 各点的电压虚部f单位:V)为(节点号从小到大排列): 0 0.0473 0.0863 0.0922 0.1076 节点电压的第C(k)次近似值: Columns 1 through 4 1.0600 1.0847 + 0.0473i 1.1123 + 0.0863i 1.1136 + 0.0922i Column 5 1.1183 + 0.1076i 迭代次数: 1 雅克比矩阵J: Columns 1 through 6 36.7747 9.6974 -5.5023 -1.5717 -5.5023 -1.5717 -10.7317 34.7555 1.5717 -5.5023 1.5717 -5.5023 -5.7052 -1.4227 44.2583 11.3874 -34.2310 -8.5339 1.4227 -5.7052 -10.6594 44.1718 8.5339 -34.2310 -5.7215 -1.3952 -34.3289 -8.3688 44.2915 11.2101 1.3952 -5.7215 8.3688 -34.3289 -10.4101 44.3915 -8.6564 -1.9888 0 0 -4.3282 -0.9944 1.9888 -8.6564 0 0 0.9944 -4.3282 Columns 7