潮流计算设计.doc

毕业设计（论文）题 目 配电网时尚计算与程序设计完毕日期年月日目录1绪论4 1.1时尚计算目和意义4 1.2时尚计算发展4 1.3本文重要工作62 配电网络模型6 2.1线路模型6 2.2变压器模型9 2.3 负荷模型133基于前推回代法时尚计算144基于matlab程序设计15 4.1 MATLAB简介15 4.2 程序设计16 4.3 算例205 总结24配电网时尚计算与程序设计摘要：在电力系统正常运行中，伴随用电负荷变化和系统运行方式变化，网络中损耗也将发生变化。要严格保证所有顾客在任何时刻均有额定电压是不也许，因此系统运行中个节点出现电压偏移是不可防止。为了保证电力系统稳定运行，要进行时尚调整。Abstract: In the normal operation of the power system, with the change of power load and the change of system operation mode, the network losses will also change. To strictly ensure that all users at any time have a rated voltage is impossible, so the node appears in the operation of the system voltage deviation is inevitable. In order to guarantee the stable operation of power system, must carry on the trend of the regulation.关键词: 电网时尚计算 前推后裔法 MATLABKeywords: Power flow calculation Push the offspring before MATLAB序言电力是以电能作为动力能源。发明于19世纪70 年代，电力发明和应用掀起了第二次工业化高潮。成为人类历史18世纪以来，世界发生三次科技革命之一，从此科技变化了人们生活。既是是当今互联网时代我们仍然对电力有着持续增长需求，由于我们发明了电脑、家电等更多使用电力产品。不可否认新技术不停出现使得电力成为人们必需品。20世纪出现大规模电力系统是人类工程科学史上最重要成就之一，是由发电、输电、变电、配电和用电等环节构成电力生产与消费系统。它功能是将自然界一次能源通过发电动力装置转化成电能，再经输电、变电和配电将电能供应到各顾客。为实现这一功能，电力系统在各个环节和不一样层次还具有对应信息与控制系统，对电能生产过程进行测量、调整、控制、保护、通信和调度，以保证顾客获得安全、经济、优质电能。电力工业发展初期，电能是直接在顾客附近发电站（或称发电厂）中生产，各发电站孤立运行。伴随工农业生产和都市发展，电能需要量迅速增长，而热能资源和水能资源丰富地区又往往远离用电比较集中都市和工矿区，为了处理这个矛盾，就需要在动力资源丰富地区建立大型发电站，然后将电能远距离输送给电力顾客。同步，为了提高供电可靠性以及资源运用综合经济性，又把许多分散多种形式发电站，通过送电线路和变电所联络起来。这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接起来整体，即称为电力系统。伴随电力系统及在线应用发展，计算机网络已经形成，为电力系统时尚计算提供了物质基础。电力系统时尚计算是电力系统分析计算中最基本内容，也是电力系统运行及设计中必不可少工具。根据系统给定运行条件、网络接线及元件参数，通过时尚计算可以确定各母线电压幅值及相角、各元件中流过功率、整个系统功率损耗等。时尚计算是实现电力系统安全经济发供电必要手段和重要工作环节，因此时尚计算在电力系统规划设计、生产运行、调度管理及科学研究中均有着广泛应用。它发展重要围绕这样几种方面：计算措施收敛性、可靠性；计算速度迅速性；对计算机存储容量规定以及计算以便、灵活等。现代电力系统提出了“灵活交流输电和新型直流输电”概念。灵活交流输电技术是指运用固态电子器件与现代自动控制技术对交流电网电压、相位角、阻抗、功率以及电路通断进行实时闭环控制，从而提高高压输电线路诉讼能力和电力系统稳态水平。新型直流输电技术是指应用现电力电子技术最新成果，改善和简化变流站造价等。运行方式管理中，时尚是确定电网运行方式基本出发点：在规划领域，需要进行时尚分析验证规划方案合理性；在实时运行环境，调度员时尚提供了电网在预想操作预想下电网时尚分布以及校验运行可靠性。在电力系统调度运行多种领域都波及到电网时尚计算。时尚是确定电力网咯运行状态基本原因，时尚问题是研究电力系统稳态问题基础和前提。1 绪论1.1时尚计算目和意义电力系统时尚计算是研究电力系统稳态运行状况一种基本电气计算。它目是根据给定运行条件和网路构造确定整个系统运行状态，如各母线上电压（幅值及相角）、网络中功率分布以及功率损耗等。电力系统时尚计算成果是电力系统稳定计算和故障分析基础。时尚计算在电网中具有重要意义： (1)在电网规划阶段,通过时尚计算,合理规划电源容量及接入点,合理规划网架,选择无功赔偿方案,满足规划水平大、小方式下时尚互换控制、调峰、调相、调压规定。(2)在编制年运行方式时,在估计负荷增长及新设备投运基础上,选择经典方式进行时尚计算,发现电网中微弱环节,供调度员平常调度控制参照,并对规划、基建部门提出改善网架构造,加紧基建进度提议。(3)正常检修及特殊运行方式下时尚计算,用于日运行方式编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定规定及电压质量规定。(4)预想事故、设备退出运行对静态安全影响分析及作出预想运行方式调整方案。总结为在电力系统运行方式和规划方案研究中，都需要进行时尚计算以比较运行方式或规划供电方案可行性、可靠性和经济性。同步，为了实时监控电力系统运行状态，也需要进行大量而迅速时尚计算。因此，时尚计算是电力系统中应用最广泛、最基本和最重要一种电气运算。1.2时尚计算发展史运用电子计算机进行时尚计算从20世纪50年代中期就已经开始。此后，时尚计算曾采用了多种不一样措施，这些措施发展重要是围绕着对时尚计算某些基本规定进行。对时尚计算规定可以归纳为下面几点：（1）算法可靠性或收敛性（2）计算速度和内存占用量（3）计算以便性和灵活性电力系统时尚计算属于稳态分析范围，不波及系统元件动态特性和过渡过程。因此其数学模型不包括微分方程，是一组高阶非线性方程。非线性代数方程组解法离不开迭代，因此，时尚计算措施首先规定它是能可靠收敛，并给出对答案。伴随电力系统规模不停扩大，时尚问题方程式阶数越来越高，目前已到达几千阶甚至上万阶，对这样规模方程式并不是采用任何数学措施都能保证给出对答案。这种状况促使电力系统研究人员不停寻求新更可靠计算措施。在用数字计算机求解电力系统时尚问题开始阶段，人们普遍采用以节点导纳矩阵为基础高斯-赛德尔迭代法（一下简称导纳法）。这个措施原理比较简朴，规定数字计算机内存量也比较小，适应当时电子数字计算机制作水平和电力系统理论水平，于是电力系记录算人员转向以阻抗矩阵为主逐次代入法（如下简称阻抗法）。20世纪60年代初，数字计算机已经发展到第二代，计算机内存和计算速度发生了很大飞跃，从而为阻抗法采用发明了条件。阻抗矩阵是满矩阵，阻抗法规定计算机储存表征系统接线和参数阻抗矩阵。这就需要较大内存量。并且阻抗法每迭代一次都规定顺次取阻抗矩阵中每一种元素进行计算，因此，每次迭代计算量很大。阻抗法改善了电力系统时尚计算问题收敛性，处理了导纳法无法处理某些系统时尚计算，在当时获得了广泛应用，曾为我国电力系统设计、运行和研究作出了很大奉献。不过，阻抗法重要缺陷就是占用计算机内存很大，每次迭代计算量很大。当系统不停扩大时，这些缺陷就愈加突出。为了克服阻抗法在内存和速度方面缺陷，后来发展了以阻抗矩阵为基础分块阻抗法。这个措施把一种大系统分割为几种小地区系统，在计算机内只需存储各个地区系统阻抗矩阵及它们之间联络线阻抗，这样不仅大幅度节省了内存容量，同步也提高了计算速度。克服阻抗法缺陷另一途径是采用牛顿-拉夫逊法（如下简称牛顿法）。牛顿法是数学中求解非线性方程式经典措施，有很好收敛性。处理电力系统时尚计算问题是以导纳矩阵为基础，因此，只要在迭代过程中尽量保持方程式系数矩阵稀疏性，就可以大大提高牛顿时尚程序计算效率。自从20世纪60年代中期采用了最佳次序消去法后来，牛顿法在收敛性、内存规定、计算速度方面都超过了阻抗法，成为直到目前仍被广泛采用措施。在牛顿法基础上，根据电力系统特点，抓住重要矛盾，对纯数学牛顿法进行了改造，得到了P-Q分解法。P-Q分解法在计算速度方面有明显提高，迅速得到了推广。牛顿法特点是将非线性方程线性化。20世纪70年代后期，有人提出采用更精确模型，即将泰勒级数高阶项也包括进来，但愿以此提高算法性能，这便产生了保留非线性时尚算法。此外，为了处理病态时尚计算，出现了将时尚计算表达为一种无约束非线性规划问题模型，即非线性规划时尚算法。近20数年来，时尚算法研究仍然非常活跃，不过大多数研究都是围绕改善牛顿法和P-Q分解法进行。此外，伴随人工智能理论发展，遗传算法、模糊算法也逐渐被引入时尚计算。不过，到目前为止这些新模型和算法还不能取代牛顿法和P-Q分解法地位。由于电力系统规模不停扩大，对计算速度规定不停提高，计算机并行计算技术也将在时尚计算中得到广泛应用，成为重要研究领域。1.3本文重要工作(1)熟悉电力系统时尚计算（前推后裔法）有关理论。(2)运用所学专业知识，提出一种合理时尚计算算法。(3)熟悉MATLAB软件使用，并根据算法编程(4)对照经典系统，对程序有效性进行验证2 配电网络模型图2-1 电力线路单相等值电路图配电网中基本元件诸多，如变压器、线路、电容器、调相机、电容器等等，配电网中基本元件诸多，如变压器、线路、电容器、调相机、电容器等等，本章重要简介线路模型、变压器模型以及负荷模型。122.1线路模型电力系统中线路模型就是以电阻、电抗、电纳、电导来表达它们等值电路。按式 求得单位长度导线电阻、电抗、电纳、电导后，就可作最原始电力线路等值电路。以单相等值电路代表三相虽已经简化了不少计算，但由于电力线路长度往往有数十乃至数百公里，如将每公里电阻、电抗、电纳、电导都一一绘于图上，所得等值电路仍十分复杂。何况，严格说来，电力线路参数是十分均匀分布，虽然是极短一段线路，均有对应大小电阻、电抗、电纳、电导。换言之虽然如此复杂等值电路，也不能认为精确。但好在电力线路一般不长，需分析又往往只是它们端点状况两端电压、电流、功率，一般可不考虑线路这种分布参数，只是在个别状况下才要用双曲函数研究具有均匀分布参数特性线路。如下，先讨论一般线路等值电路。所谓一般线路，指中等及中等如下长度线路。对架空线路，这长度大概为300km；对电缆线路，大概为100km。线路长度不超过这些数值时，可不考虑它们分布参数特性，而只用将线路参数简朴地集中起来电路表达。在如下讨论中，R（），X（），G（S），B（S）分别表达全线路每相总电阻、电抗、电纳、电导。显然线路长度为l(km)时 (2-1)一般线路中，又有短线路和中等长度线路之分。 图2-2 短线路等值电路图Z所谓短线路，是指长度不超过100km架空线路。线路电压不高时，这种线路导纳B影响一般不大，可略去。从而，这种线路等值电路最简朴。只有一串联总电抗，如图2-2所示。 显然，假如电缆线路不长，电纳影响不大时，也可采用这种等值电路。由图2-2可得 (2-2)将式(2-2)与电路理论课程中简介过两端口或四端网络方程式 (2-3)相比较,可得这种等值电路通用常数A、B、C、D （2-4）所谓中等长度电路，是指长度在100km 300km之间架空线路和不超过100km电缆线路。这种线路电纳B一般不能略去。这种线路等值电路有二种形等值电路和T形等值电路，如图2-3（a）（b）所示。其中，一般是形等值电路。在形等值电路中，除串联线路总阻抗外，还将线路总导纳分为两半，分别并联在线路始末端。在T形等值电路中，线路总导纳集中在中间，而线路总阻抗则分为两半，分别串联在它两侧。因此，这两种电路都是近似等值电路，并且，互相间并不等值，即它们不能用-Y变换公式互相变换。由图2-3（a）可得，流过串联阻抗Z电流为，从而 (2-5)流入始端导纳电流为，从而 (2-6)由此又可得 (2-7)将式(2-7)与式(2-3)相比较,可得这种等值电路通用常数 (2-8)Z/2 Z/2ZY (a) (b)(a) 形等值电路(b) T形等值电路图2-3 中等长线路等值电路图相似可以到图2-3(b)所示等值电路通用常数 (2-9)2.2变压器模型（1）双绕组变压器参数和数学模型由于变压器短路损耗近似等于额定电流流过变压器时高下压绕组中总铜耗，即而铜耗与电阻之间有如下关系 可得式中，、以V、VA为单位，以W为单位。如改以kV，改以MVA为单位，则可得 （2-10）式中 变压器高下压侧绕组总阻抗（）； 变压器短路损耗（kW）； 变压器额定容量（MVA）； 变压器额定电压（kV）；在电力系记录算中，求取变压器电抗措施和电机学课程中简介略有不一样。由于大容量变压器阻抗中以电抗为主，亦即变压器电抗和阻抗数值上靠近相等，可大体认为变压器短路电压百分值%与变压器电抗有如下关系%从而 （2-11）式中 变压器高下压侧绕组总电抗（）； %变压器短路电压百分值； 、代表意义与上面相似。变压器励磁支路有两种表达方式，即以阻抗表达和以导纳表达。前者在电机学课程中常用，后者则在电力系记录算中常用。它们分别示于图2-4（a）（b）。而与之对应空载运行时电压、电流相量图则示于图2-5（a）（b）。 变压器励磁支路一导纳表达时，其对应是变压器铁耗。因变压器铁耗近似与变压器空载损耗相等，电导也可于空载损耗相对应。而由图2-4（b）可见，两者之间有如下关系 （2-12）式中 变压器电导（S）； 变压器空载损耗（kW）； 变压器额定电压（kV）（a）励磁支路以阻抗表达时（b）励磁支路以导纳表达图2-4 双绕组变压器等值电路图由图2-5（b）可见，变压器空载电流中流经电纳部分占很大比重，从而，它和空载在电流在数值上靠近相等，可以替代求取变压器电纳。亦即，由于 （2-13）而 （2-14）而得 （2-15）将代入，最终得 （2-16）式中 变压器电纳（S）； 变压器空载电流百分数； 、代表意义与（2-10）同。求得变压器阻抗、导纳后，即可作变压器等值电路。变压器等值电路有两种，即形等值电路和T形等值电路。在电力系记录算中，一般用形等值电路，且将励磁支路接在电源侧。这种等值电路就如图2-4（b）所示。（2） 三绕组变压器参数和数学模型计算三绕组变压器各绕组阻抗措施虽与计算双绕组变压器措施没有本质区别，不过由于三绕组变压器各绕组容量比有不一样组合，而各绕组在铁芯上排列又有不一样方式，计算时需注意。 三绕组变压器按三个绕组容量比不一样有三种不一样类型.第种为100/100/100，即三个绕组容量都等于变压器额定容量；第种100/100/50，即第三绕组容量仅为变压器额定容量 50%；第种为100/50/100，即第二绕组容量仅为变压器额定容量50%。目前已在系统中使用三绕组变压器，从制造厂搜集到往往是它三个绕组两两作短路试验时测得短路损耗。如该变压器属第一类型，可由提供短路损耗、直接按下式求取各绕组短路损耗 （2-17）然后按与双绕组变压器相似公式计算各绕组电阻 （2-18）如该变压器属第、第种类型，则制造厂提供短路损耗数据是一对绕组中容量较小一方到达它自身额定容量，即时值。这时，应首先将各绕组间短路损耗数据归算为额定电流值，在运用上列公式求取各绕组短路损耗和电阻。例如，对100/50/100类型变压器，制造厂提供短路损耗、都是第二绕组中流过它自身额定电流，即二分之一变压器额定电流时测得数据。因此，应首先将它们归算到对应于变压器额定电流 （2-19）然后再按式（2-17）（2-18）计算。求取三绕组变压器导纳措施和求取双绕组变压器导纳措施相似。2.3 负荷模型在电力系统稳态分析中，负荷数学模型最简朴，就是以给定有功功率和无功功率表达。只有在对计算精度规定较高时，才需计及负荷静态特性。负荷静态特性可以用函数或多项式表达，如静态电压特性可为也可为:式中 在额定电压下有功功率、无功功率负荷； P、Q电压偏离额定值时有功功率、无功功率负荷；待定系数，它们数值可通过拟合对应特性曲线而得。一般可将与节点有关负荷模型描述为 （2-23）式中，U为节点实际电压；为节点参照电压。假如式（3.23）中，S为恒功率负荷；假如，S为恒电流负荷；假如，S为恒阻抗负荷。为了讨论以便，假定S为恒阻抗负荷，则有 （2-24）因此，可以将节点恒阻抗表达为 （2-25）式中为节点电压。一般认为节点负荷为恒功率，对于运行在正常工作条件下配电系统，其节点电压变化幅度在5%以内，可以认为节点电压是恒定，此时恒功率负荷可以作为恒阻抗负荷来处理。3基于前推回代法时尚计算已知配电网始端电压和末端负荷，以馈线为基本计算单位。最初假设全网电压都为额定电压，根据负荷功率由末端向始端逐段推算，仅计算各元件中功率损耗而不计算节点电压，求得各支路上电流和功率损耗，并据此获得始端功率，这是回代过程；再根据给定始端电压和求得始端功率，由始端向末端逐段推算电压降落，求得各节点电压，这是前推过程。如此反复上述过程，直至各个节点功率偏差满足容许条件为止。前推回代法基于支路电流进行，首先假定各节点电压幅值为1，幅角为0，详细计算环节为：（1）从第一层节点开始，根据基尔霍夫电流定律，求支路上电流： （3-1）式中，是节点 j 功率，是节点 j 电压。（2）从第二层开始逐层计算非末梢节点注入电流，根据基尔霍夫电流定律应等于（1-1）式与该节点流出电流之和： （3-2）（3）由环节1）和2）可求出所有支路支电流，再运用已知根节点电压，从根节点向后顺次求得各个负荷节点电压 （3-3）其中i为父节点，j为子节点，Z为i、j间支路阻抗。（4）计算各节点电压幅值修正量： （3-4）（5）计算节点电压修正量最大值；（6）鉴别收敛条件： （3-5）其中k为迭代次数，若最大电压修正量不不小于阈值，则跳出循环，输出电压计算成果；否则反复环节（1）（6），直到满足（1-5）式条件为止。（7）在得到各个节点电压电流后，就可以计算线路时尚S： （3-6）总之，（3-1）（3-3）式为一组递归方程，对树进行前向遍历，从树末梢节点出发，运用已知负荷功率，逐一计算（3-1）（3-2）式，即可求得根节点处电流，再从根节点出发，对树进行后向遍历，用（3-3）式可求各节点电压。这样就完毕一次前推后裔计算，迭代反复进行，直至满足收敛原则为止。4 基于matlab程序设计4.1 MATLAB简介MATLAB，是美国MathWorks企业出品商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算高级技术计算语言和交互式环境，重要包括MATLAB和Simulink两大部分。MATLAB是矩阵试验室（Matrix，Laboratory）简称，和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创立顾客界面、连接其他编程语言程序等，重要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。MATLAB基本数据单位是矩阵，它指令体现式与数学、工程中常用形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完毕相似事情简捷得多，并且mathwork也吸取了像Maple等软件长处，使MATLAB成为一种强大数学软件。在新版本中也加入了对C，FORTRAN，C+，JAVA支持。可以直接调用，顾客也可以将自己编写实用程序导入到MATLAB函数库中以便自己后来调用，此外许多MATLAB爱好者都编写了某些经典程序，顾客可以直接进行下载就可以用。在科学研究和工程应用中，往往要进行大量数学计算，其中包括矩阵运算。这些运算一般来说难以用手工精确和快捷地进行，而要借助计算机编制对应程序做近似计算。目前流行用Basic、Fortran和c语言编制计算程序,，既需要对有关算法有深刻理解，还需要纯熟地掌握所用语言语法及编程技巧。对多数科学工作者而言，同步具有这两方面技能有一定困难。一般，编制程序也是繁杂，不仅消耗人力与物力,并且影响工作进程和效率。为克服上述困难，美国Mathwork企业于1967年推出了“Matrix，Laboratory”（缩写为Matlab）软件包，并不停更新和扩充。MATLAB产品族可以用来进行如下多种工作：数值分析数值和符号计算工程与科学绘图控制系统设计与仿真数字图像处理数字信号处理通讯系统设计与仿真财务与金融工程MATLAB应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加工具箱（单独提供专用MATLAB函数集）扩展了MATLAB环境，以处理这些应用领域内特定类型问题。MATLAB可进行如下运算：（1）与常数运算常数与矩阵运算即是同该矩阵每一元素进行计算。但需注意进行数除时，常数一般只能做除量。基本函数运算中，矩阵函数运算是矩阵预算中最实用部分，常用重要有如下几种： det(a) 求矩阵a行列式 eig(a) 求矩阵a特性值 inv(a)或a(-1) 求矩阵a逆矩阵 rank(a) 求矩阵a秩 trace(a) 求矩阵a迹（对角线元素之和） 我们进行工程计算时常常碰到矩阵对应元素之间运算。这种运算不一样于前面讲数学运算，为有所区别，我们称之为数组运算。（2）基本数学运算 数组加、减与矩阵加、减运算完全相似。而乘除法运算有相称大区别，数组乘除法是指两同维数组对应元素之间乘除法，它们运算符为“.*”和“./”或“.。”前面讲过常数与矩阵除法运算中常数只能做除数。在数组运算中有了“对应关系”规定，数组与常数之间除法运算没有任何限制。 此外，矩阵数组运算中尚有幂运算（运算符.）、指数运算（exp）、对数运算(log)、和开方运算（sqrt）、等，有了“对应元素”规定，数组运算实质上就是针对数组内部每个元素进行。矩阵幂运算与数组幂运算有很大区别。（3）逻辑关系运算逻辑运算是MATLAB中数组运算所特有一种运算形式，也是几乎所有高级语言普遍合用一种运算。4.2 程序设计(1)程序运算流程图如图4-1所示图4-1(2)程序代码clc;clear all;% 第一列为首节点，第二列为末节点，第三列为支路电阻，第四列为支路电抗BranchData = 0 1 0.933 0.655; 1 2 1.203 0.453; 2 3 1.279 0.521; 3 4 2.145 0.964; 4 5 3.022 1.477; 5 6 4.5 1.223; 6 7 2.89 0.818; 7 8 1.249 0.384; %第一列为节点号，第二列为负荷有功，第三列为负荷无功NodeData = 1 98 98; 2 70.00 60.00; 3 85.00 80.00; 4 55.00 55.00; 5 40.00 40.00; 6 50.00 50.00; 7 40.00 30.00; 8 20.00 15.00; UB = 10; % 基准电压 kVNN = 9; % 节点数SB = 100; % 基准功率 MVAZB = UB2/SB; % 基准阻抗 ohm nb=size(BranchData,1); %支路数BranchData(:,3,4) = BranchData(:,3,4) / ZB; % 阻抗标幺化NodeData(:,2,3) = NodeData(:,2,3) / SB / 1000; % 功率标幺化BranchData(:,1,2) = BranchData(:,1,2) + 1; % 节点从1开始A0 = zeros(NN); for n = 1:(NN-1) A0(BranchData(n,1),BranchData(n,2) = 1;end % 形成 A0i=sqrt(-1);A0T = A0; % 形成 A0 转置S = 0;-NodeData(:,2) - i*NodeData(:,3); % 形成 SZL = 0;BranchData(:,3) + i*BranchData(:,4); % 形成 ZLV = ones(NN,1);V(1) = 1; % 各个节点电压赋初值IL(NN,1) = -conj(S(NN) / V(NN); % 最末支路电流赋初值Delta = 1; % 收敛判据赋初值TempV = V; % 赋初值，用于记忆上次迭代成果k=0;tic %计时while Delta 1e-6 IN = conj(S ./ V); % 节点注入电流 for n = 1:NN-1 IL(NN-n) = A0(NN-n,NN-n+1:end) * IL(NN-n+1:end) - IN(NN-n); end % 电流回代过程 for n = 2:NN V(n) = A0T(n,1:n-1) * V(1:n-1) - ZL(n) * IL(n); end % 电压前推过程 Delta = max(abs(V-TempV); % 更新收敛判据 TempV = V; % 记忆迭代成果 k=k+1; if k=30 break; endendtocv=abs(V); %取电压幅值theta=angle(V)*180/pi; %取电压相角Sij=zeros(nb,1); %用来存储损耗for n=1:nb %计算支路损耗 Sij(n)=(V(BranchData(n,1)-V(BranchData(n,2)*. conj(V(BranchData(n,1)-V(BranchData(n,2)/ZL(n+1);endPloss=real(Sij); %有功损耗Qloss=imag(Sij); %无功损耗if k=30disp( 节点 电压幅值 电角度); for m=1:NN fprintf(%6d%12.4f%13.4fn,m,v(m),theta(m);end fprintf(n); disp( 首节点 末节点 有功损耗 无功损耗); for m=1:nb fprintf(%6d%8d%12.4f%13.4fn,BranchData(m,1),. BranchData(m,2),Ploss(m),Qloss(m);End fprintf(n); fprintf(迭代次数 ：%dn,k);else disp(不收敛！);end4.3 算例验证（1）算例1图4-2给出了简朴配电网示意图，网络构造参数由表1和表2分别给出，运行成果见表4，运用此措施进行迭代运算。图4-2将运行成果表3与文献数据表4做对比，误差在合理范围内，阐明此算法有效。（2）算例2以IEEE33节点系统为例，其网络如图4-3所示，其支路参数见表5，运行成果见表6。图4-3表5 IEEE33节点支路参数将运行成果表6与文献表7做比较，误差在合理范围内，阐明此算法有效。5 总结致 谢参 考 文 献1 马瑾.辐射状配电网时尚计算. 2月.沈阳工业大学学报2 张荣.改善带二阶项配电网迅速时尚算法.7月.电工技术学报3 陈旺虎.配电网时尚前推回代法一种简朴实现措施.7月.江苏电机工程