电力系统暂态分析课程设计

内蒙古科技大学 电力系统稳态分析课程设计

题 目：基于MAＴLＡB的电力

系统ﻩ复杂潮流分析

学生姓名:刘建峰

学 号：11６7１3020７ 专 业:电气工程及其自动化 班 级：电气201１—２班 指导教师：刘景霞

摘要

电力系统潮流计算是电网分析的基础应用，是对复杂电力系统正常和故障条件下稳态运行状态的计算。给定电力系统的网络结构、参数和决定电力系统运行状况的边界条件,确定电力系统运行的方法之一是朝流计算。

MATLAB是一种交互式、面向对象的程序设计语言,广泛应用于工业界与学术界，主要用于矩阵运算.采用迭代法,通过建立矩阵的修正方程来依次迭代,逐步逼近真值来计算出电力网的电压，功率分布。

PQ分解法是极坐标形式牛顿-拉潮流计算的一种简化计算方法，。Ｐ—Q分解法通过对电力系统具体特点的分析，对牛顿法修正方程式的雅可比矩阵进行了有效的简化和改进。由于这些简化只涉及修正方程式的系数矩阵，并未改变节点功率平衡方程和收敛判据，因不会降低计算结果的精度。

用手算和计算机算法对其进行设计。使用MATLＡB软件进行编程,在很大程度上节省了内存,减少了计算量。通过对本题计算我们了解了一些工程计算和解决工程问题的方法。

关键词：潮流计算，PQ分解法,MATLAB

Electricａl powｅr sｙｓteｍ ｃompleｘ tidal currenｔ

aｎａlysｉs ｂasｅd ｏｎ ＭATLAＢ

Ｐoweｒ Flow Anaｌysｉｓ Gｒid ｃｏmpｕting is the basiｓ of applicatｉｏns， thｅ ｃomｐlex power system uｎdｅr normaｌ and fault cｏｎdｉtions fｏr ｔhｅ caｌｃulatｉoｎ ｏf ｓteady sｔａtｅ oｐｅratｉoｎ. Giｖｅn ｔｈe poｗer ｓyｓtｅｍ netｗoｒk sｔructuｒｅ， pａｒａmeters aｎd dｅcisions opｅrａtion of ｔhｅ power ｓystｅm bｏｕｎdary cｏnditions, to dｅtermine tｈｅ mｅthod ｏf ｏｐｅratiｏn of tｈｅ pｏwｅr system is one of Nortｈ Korea ｆlow cａｌculａｔｉoｎ.

MＡＴLＡＢ is an ｉｎｔeracｔive, object-orienteｄ progｒａmｍing ｌanguａgｅ, wｉｄely usｅd in induｓtry ａnd ａcadeｍia, maiｎｌy for ｍatrix ｃaｌculation. Ｕsing iｔeration, ｔｈe amenｄmｅｎｔ tｈrouｇh ｔhe eｓｔaｂlisｈment of matrix iteratiｖe equａｔion to turｎ, gradｕally moving towａrdｓ a true ｖａlｕｅ to calｃｕｌate the voltage ｅlecｔrｉcｉｔｙ grid, ｐoｗer distrｉbuｔiｏn.

PQ decｏmpositioｎ mｅthod is thｅ form oｆ polａr coordinates Newton － the wideｎｉｎg ｔrend of ａ siｍｐlifiｅｄ cａｌculatiｏn meｔhod. P-Q dｅcomposｉｔiｏn metｈｏd aｄｏpted by thｅ sｐecific cｈａractｅｒistiｃs of the power systｅm aｎaｌysｉｓ, Nｅwｔoｎ’s Lａｗ of the Jacｏbian mａtrｉx formula hａs effectｉvelｙ simｐlified ａnd iｍprovｅd. Aｓ a result of thｅｓe ｓimｐlifｉed formuｌａ tｈａt involｖeｓ only the coefficienｔ maｔrｉx， the baｌａｎce oｆ ｐoｗer has not cｈaｎged node eqｕatｉons anｄ ｔhe ｃoｎvｅrｇｅｎcｅ ｃriｔerｉon， becauｓe the ｒeｓｕlts wｉll noｔ rｅdｕce ｔｈe acｃuｒacｙ.

Use MATLAB softwarｅ progｒammｉｎg, saving memory to ａ ｌarｇe ｅxtent, rｅducｅ tｈe aｍount of computａtion. Ｂy thｉs calculatiｏｎ we undｅrstand thａt ａ numbeｒ ｏｆ engiｎｅerｉng calculａtiｏｎ and sｏlｖe

enginｅering ｐroblemｓ

Keｙｗordｓ ： Thｅ ｔreｎd， ｔｈｅ PＱ decomposiｔiｏｎ metｈod, MATLAB

目录

内蒙古科技大学课程设计任务书ﻩ错误!未定义书签。

第一章 引言 ....................................................................................................... 错误!未定义书签。 1.1 研究背景及意义 ............................................................................................ 错误!未定义书签。 1．２ 潮流计算的意义 ..................................................................................... 错误!未定义书签。 1.3电力系统稳态分析潮流计算总结ﻩ错误!未定义书签。

1.4 MATＬAB的概述 ......................................................................................... 错误!未定义书签。 第二章 理论计算 ................................................................................................... 错误!未定义书签。 ２.1 Ｐ-Q法潮流计算的基本步骤ﻩ错误!未定义书签。

2.２ PQ分解法潮流计算流程图 ....................................................................... 错误!未定义书签。 ２.3 两机五节点网络潮流计算 .......................................................................... 错误!未定义书签。 第三章 程序设计ﻩ错误!未定义书签。 3．1 设计流程ﻩ错误!未定义书签。

3.2 程序设计 ........................................................................................................ 错误!未定义书签。 3.3 程序运行结果 ................................................................................................ 错误!未定义书签。 第四章 设计感想ﻩ错误!未定义书签。

参考文献 ........................................................................................................ 错误!未定义书签。

内蒙古科技大学课程设计任务书

课程名称 设计题目 指导教师 一、教学要求 电力系统稳态分析课程设计以设计和优化电力系统的潮流分析为重点，提高学生综合能力为目标，尽可能结合实际工程进行。设计内容的安排要充分考虑学校现有的设备，设计时间及工程实际需要，并使学生初步学会运用所学知识解决工程中的实际问题。 电力系统稳态分析课程设计 两机五节点网络潮流计算— P－Q法 刘景霞 时间 1周 二、设计资料及参数 (一）设计原始资料 1、待设计电气设备系统图 2、电力系统网络各元件参数 3、电力系统电气元件的使用规范 4、电力工程电气设计手册 (二)设计参考资料 １、《电力系统稳态分析》，陈珩,中国电力出版社，2007,第三版 ２、《电力系统分析》，韩祯祥，浙江大学出版社，２005,第三版 3、《电力系统分析课程实际设计与综合实验》,祝书萍，中国电力出版社，２0０7,第一版 三、设计要求及成果 １.根据给定的参数或工程具体要求,收集和查阅资料；学习相关软件（软件自选)。 2．在给定的电力网络上画出等值电路图。 ３.运用计算机进行潮流计算。 ４.编写设计说明书。 基本要求： 1.编写潮流计算程序; ２．在计算机上调试通过(？）； 3.运行程序并计算出正确结果（？）； 4.写出课程设计报告(包括以下内容）(1份） (1)程序框图；(2)源程序;（３）符号说明表;（４）算例及计算结果 ５.编写计算说明书(1份)。 四、进度安排 根据给定的参数或工程具体要求,收集和查阅资料(半天） 学习软件（MATLAＢ或C语言等)（一天半) 编程计算复杂系统潮流计算(三天) 编写计算设计书（一天） 五、评分标准 课程设计成绩评定依据包括以下几点: 1） 工作态度（占1０%); 2） 基本技能的掌握程度（占２0%）； 3） 程序编写是否合理是否有运行结果(４0%); 4） 课程设计说明书编写水平(占30%）。 5） 分为优、良、中、合格、不合格五个等级。 考核方式:设计期间教师现场检查；评阅设计报告

G0.45+j0.150.4+j0.0510.08+j0.240.02+j0.0630.01+j0.03.180.06+j0420.04+j0.12G-(0.2+j0.2)50.6+j0.10.08+j0.24

.180j+0.06第一章 引言

1．1 研究背景及意义

电力系统在运行时，在电源电势激励作用下,电流或功率从电源通过系统各元件流入负荷，分布于电力网各处,称为潮流分布。

潮流计算是电力系统分析中的一种最基本的计算，即根据给定电力系统的网络结构和参数,在满足电力系统运行状态参数的边界条件情况下,确定电力系统稳态运行状态的最基本的方法。它的任务是在给定的接线方式和运行条件下,确定系统的运行状态,如各母线上的电压（幅值和相角）、网络中的功率分布及功率损耗等，是电力系统的稳态计算。潮流计算是对电力系统正常运行状况的分析和计算，即电力系统中的电压、电流、功率的计算，即潮流计算;潮流计算方法很多：高斯—塞德尔法、牛顿—拉夫逊法、Ｐ-Q分解法、直流潮流法,以及由高斯—塞德尔法、牛顿—拉夫逊法演变的各种潮流计算方法。

潮流计算可以用传统的手工方式进行，也可以计算机为工具通过软件完成。两种方法各有优缺点。前者物理概念清晰，可用来计算一些接线较简单的电力网,但若将其用于接线复杂的电力网则计算量过大，难于保证计算准确性。后者从数学上看可归结为用数值方法解非线性代数方程，数学逻辑简单完整,借助计算机可快速精确地完成计算,但其缺点是物理概念不明显,物理规律被埋没在循环往复的数值求解过程中。

潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出

的各种问题。对规划中的电力系统,通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求:对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各元件（线路、变压器等）是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。因此潮流计算的目的是：

① 为电力系统规划设计提供接线、电气设备选择和导线截面选择的依据。 ② 提供电力线运行方式和制定检修计划的依据。 ③ 提供继电保护、自动装置设计和整定计算的依据。

④ 为调压计算、经济运行计算、短路和稳定计算提供必要的数据。

潮流计算的研究是从20世纪50年代随着电网的产生而开始的,涌现出各种算法都是围绕着算法的可靠性或收敛性;对计算机内存量的要求和计算速度;计算的方便性和灵活性等。

1.２ 潮流计算的意义

潮流计算是研究电力系统稳态运行情况的一种基本电气计算,常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统,通过潮流计算,可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。

具体表现在以下方面：

(1)在电网规划阶段，通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案，满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。

（２)在编制年运行方式时,在预计负荷增长及新设备投运基础上，选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。

(３）正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制，指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。

(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。

总结为在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时,为了实时监控电力系统的运行状态,也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时,采用离线潮流计算;在电力系统运行状态的实时监控中，则采用在线潮流计算。

1.3电力系统稳态分析潮流计算总结

潮流计算针对电力革统各正常运行方式,而静态安全分析则要研究各种运行方式下个别系统元件退出运行后系统的状况。其目的是校验系统是否能安全运行,即是否有过负荷的元件或电压过低的母线等。原则上讲,静态安全分析也可U用潮流计算来代替。但是一般静态安全分析需要校验的状态数非常多，用严格的潮流计算来分析这些状态往往计算量过大,因此不得不寻求一些特殊的算法以满足要求。

利用电子数字计算机进行电力系统潮流计算从20 世纪５0 年代中期就己开始，此后，潮流计算曾采用了各种不同的方法,这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的，对潮流计算的要求可以归纳为下面几点:

(1）计算方法的可靠性或收敛性 (2)对计算速度和内存量的要求

（3)计算的方便性和灵活性

一般情况下，采用P-Q分解法计算时要求的迭代次数较采用牛顿-拉夫逊法时多，但每次迭代所需的时间则较采用牛顿-拉夫逊法时少，以致总的计算速度仍是P－Q分解法快。

１.4 MATＬAB的概述

目前电子计算机已广泛应用于电力系统的分析计算，潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（1)计算速度快（2)内存需要少（3)计算结果有良好的可靠性和可信性（４）适应性好,亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强(5）简单。

MATLAB是一种交互式、面向对象的程序设计语言，广泛的应用于工业界和学术界,主要用于矩阵运算,同时在数值分析、自动控制模拟、MATＬＡＢ程序设计语言结构完整,且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效地解决工业计算的问题，特别是关于矩阵和矢量的计算。MATＬＡB与C语言和FORＴRＡN语言相比更容易掌握。通过ＭATLAＢ语言，可以用类似数学公式的方式来编写算法,大大降低了程序所需的难度并节省了时间，从而可以把主要的精力集中在算法的构思而不是编程上。

另外,MATLAB提供了一种特殊的工具:工具箱（TOOＬBＯXES）.这些工具箱主要包括：信号处理（SIＧNAL PROＣESSINＧ)、控制系统(CＯNTＲOL SYSTEMS）、神经网络（ＮEURＡＬ NEＴWOＲＫS)、模糊逻辑(FUZZY LOGIC)、小波(ＷAVELETS)和模拟（ＳＩMULＡTION）等等。不同领域、不同层次的用户通过相应工具的学习和应用,可以方便地进行计算、分析及设计工作。 ＭAＴＬAB设计中,原始数据的填写格式是很关键的一个环节,它与程序使用的方便性和灵活性有着直接的关系。原始数据输入格式的设计,主要应从使用的角度出发,原则是简单明了，便于修改。

第二章 理论计算

2．1 Ｐ－Q法潮流计算的基本步骤 运用P-Q分解法计算潮流分布时的基本步骤是：

（1） 形成系数矩阵Bˊ、Ｂ〞，并求其逆阵。

（2） 设各节点电压的初值i0(i＝1,2，…，n，i≠s)和Ui0(i=１，2，…，m，

i≠s)。

0（3） 按式Pi=pi-Ui0Uj0GijcosijBijsinij0计算有功功率的不平

0j1jn衡量Pi0,从而求出P（4） 解修正方程式 （P0U0(i=１，2，…，n≠s)。

U）=﹣BˊU0，求各节点电压相位角的变

量i0（ｉ=1，2，…,n，ｉ≠s)。

（5） 求各节点电压相位角的新值i1=i0＋i0(i=１,2，…，n，i≠ｓ)。

1（6） 按式Qi=Qi-Ui0Uj0GijsinijBijcosij1 计算无功功率的

0j1jn0Q不平衡量Qi,从而求出

0U0（i＝1,2，…,m,i≠s）。

（7） 解修正方程式(QU)=﹣B〞U（i=１,2,…，m，i≠s）。

（8） 求各节点电压大小的新值Ui1＝Ui0+Ui0（i=1,2，…，ｍ，i≠ｓ）。 （9） 运用各节点电压的新值自第三步开始进入下一次迭代。

（10）计算平衡节点功率和线路功率。

2.2 PQ分解法潮流计算流程图

启动 输入原始数据 形成节点导纳矩阵 形成BˊB＂并求其逆阵 置迭代次数k=0 k+1 计算PiQi 计算电压新值新值是否满足要求 满 停止 足 不满足

2.3 两机五节点网络潮流计算

G0.45+j0.150.4+j0.0510.08+j0.240.02+j0.0630.01+j0.03.180.06+j0420.04+j0.12G-(0.2+j0.2)系统接线图

50.6+j0.10.08+j0.24

.180j+0.06解:

1，.形成系数矩阵Bˊ、B〞并求他们的逆阵

由于节点1为平衡节点,其他的节点均为PQ节点，系数矩阵B'，B''阶数相同。又应对该等值网络，不存在除去与有功功率和电压相位或无功功率和电压大小关系较小的因素的可能性,这两个矩阵Ｂˊ、B〞完全相同。他们就由导纳矩阵的虚部部分中除第一行第一列外的各个元素所组成,即

Ｙ＝[6.２5-18.７5j -５+15j -1.2５+3．75j ０ 0

-5＋１5j 10.834－32．5j -1．６６７+5ｊ －1.6６７+5j -２.5+7.5ｊ

－1．25+3.75j -1.6６7+5j １２.９17-38.75j -10+3０j 0

0 -1．6６７+5j －10+30j 12.917-３8.７5j －1．25+３.７5ｊ

０ －2.5＋７．5j 0 -１.25+3.75j 3.75-1１．２5j]

Bˊ=B〞=

［ －３2.５00 5.000 ５．000 ７.500

5．00０ -38.7５0 30.0０0 0

５.000 3０.000 -38.750 3.7５0 ７.50０ 0 3．75０ －11.250 ]

由此可见,网络的节点导纳矩阵虽为奇异矩阵,但它的虚数部分的子阵Bˊ或Ｂ〞则是非奇异矩阵,可以求逆,其逆阵为

(B')1＝(B'')1=

[ -0．05619０ -0.０４19０５ －０.44762 -0.0５２381

-0.0419０5 -0．０99０48 －０.0８76１9 -0.057143

-0.04４762 -０．08７６19 -0.１05７14 -0.０6507９

-0.05238１ -0。0571４3 -0.065079 -０.１45５03 ]

2．计算各节点有功功率的不平衡量Pi

(0)(0) 取U1=1．06,1=0;U2=U3(0)=U4=U5(0)=１.0;2(0)=3(0)=4(0)＝5(0)=0,按下

式计算各节点有功功率不平衡量

0 Pi=pi-Ui0Uj0GijcosijBijsinij0

0j1jnP20=0．20—1.0×1.06（-5.000cos０+15.０0０ｓin0)—1.0×1.0(1０.3８４

cos0—32.500sｉn0)—1.0×１．0(-1.667ｃoｓ0＋5.000ｓin0）—1.0

×1．（－1.６６７coｓ０+5.０0０ｓin０）—１．0×１.0（-2.５０0cos０+7.500siｎ0)＝0.50000 相似地可得

P30＝-0.3７50００；P40＝-0.４０0000；P50＝-0.60000０ 3．计算各节点电压的相位角i（弧度） 由方程式 —(B')1 （P可得

U)=U0

2(0)＝—0.036９52;3(0)=—0.０８55２4；4(0)=—0．０91８１0；5(0)＝—０.108５７１; 所以

2(1)=—０.０3６9５2； 3(1)=—0.0855２4； 4(1)=—0.0918１0； 5(1)＝—0.1０85７１；

4.计算各节点无功功率不平衡量Qi 按下式计算各节点无功功率不平衡量

1 Qi=Qi-Ui0Uj0GijsinijBijcosij1

0j1jn0Q2=０.２0—1.0×1.06[—5.0０sin(－０．０36952—0）—１５.000cos（-0.0３

６952—0)]—１.０×1.0[10.8３４sin0＋３2.5cos０］—1.0×1．0[-1.667ｓin（-０.０36952+0.０8５524）—5.０00ｃos（-0.０369５2＋０.085524)]—1．0×1.０［－１．６６7sin(—０.0３６952+0,0９1810）—5.000cｏs（－0,0369５2＋0.09１81０)]—１．0×１.0[-2.５00ｓｉｎ(—０.0369５2＋0.108571）—7.50０coｓ(-０.036952+０.108571)］=１.211930 相似地可得

0 Q30=—0.077279; Q4=—0.1９1947; Q50＝—０.３1959９；

５.计算各节点电压的大小Ui 由方程式

—(B'')1 (Q可得

U)＝U

00U2=１.039528； U30=０.008050; U4＝0.０0638６； U50=0.00００72;

所以

11U2=１.039５28； U31=0．008０５0； U4=0.00６38６； U51=0.０000

７2;

求得各节点电压的新值后,就开始第二次迭代。每次迭代所得适于表1～表３。由表1可见经6次迭代就可满足ε≤105的要求。 表1 迭代过程中各节点功率的不平衡量

Ｋ P2(k) -0.00000 ５kQ2 1.２1１930 -0．136４９7 0.014875 P3(k) -0.37５0００ ０.０4９８０2 －0．00589９ Q3k P4(k) -０.077２79 ０.018３０0 －0.0０19３6 0.000２１5 －0,00０02４ 0.000８72 -０.00０096 ０.0０0０1１ -0．００000１ 0.400０0０ 0.07６2４8 kQ4 －0.1９１947 ０.０P5(k) －0.600000 Q5k -０.319599 ０ 1 -0.４1１７２0 ０.120899 ０．040218 24849 －0.０1258-0．０04265 2 0.０43223 -0.007856 －0．002664 ８ 0．00０29５ -０.０0003３ ０.００.０0１376 －０.０00150 0.000０１6 0.0００46６ －0.00０051 0．000006 3 －０.004833 -0．０0１658 －0.0００６50 4 0.000５37 ０.000184 -0.000０72 5 -0．00006０ -０.０000２１ 0．0０0000 0.0０000８ 0.00００03 00004 ０.０0000００ -0．000002 ０．000000 6 0.0０0００７ -０．000001 0.0０0０00

表2 迭代过程中各节点电压的修正量

k 2(k) kU2 3(k) U3k 4(k) kU4 5(k) U5k 0 －0．036952 0.039５28 -０.０8５524 ０.０0８0５0 ０.０0０7５7 －0.0918１０ 0．0024６１ -０．000２85 -0．00００28 -0．０0０００3 ０．006386 －0.10８571 ０.000941 0.004597 0.０0007２ 0.00161８ 1 -０.０1０063 -0.003406 0.001８28 2 0.００1054 ０．０0０385 -０．0０0２30 －0.00００６4 0.0000０６ -０.00000１ -０.０000８2 ０.０0000８ -0.0０0001 -０.０００484 0.０0０049 -0.00005 ０-0.000149 3 －0.0０0120 -０.000043 0.０00０２３ 0．０0０0０5 -0.000002 0.00015 ０4 0.0０00１3 -０.０0０0０２ ５ -0.000０02 -０.00000１ ０．000０00 ０.00０00０ 0.000000 ０.0０0000 0.00001 ０0.000０00 ６

表3 迭代过程中各节点电压

k 2(k) 0.000000 kU2 １.000000 3(k) 0.0000０0 U3k １.00000０ 4(k) 0.000000 kU4 1.00000０ 5(k) ０.０00０００ －0，１0857１ －0,103974 U5k 1.0０00０0 １.00007２ 1．001690 0 1 －0.036952 1.039528 -０.0８５524 1．０08０５0 1.00８8０８ -0.0９１810 １．00６３８6 2 -0．０470１6 1.03612２ -0.083６96 -0,089348 1.007327 3 -0.０４59６2 １．0３6507 -０.０8392６ 1．00８７４４ １.008７50 -０.０89633 -０.08９605 1.0072４5 1．0０72５3 1.00７-0．10４458 -0.104４10 -0.10４1.00154０ 4 -0.04６0８2 1．0364６3 -0.０83903 1．001５5５ 1.00１554 5 -０.046061．０36468 -0.0８39０1.00８７５-０.0８８ ６ －0.0460７0 6 1．0364６8 -0,083906 0 9608 252 1.0072５2 415 -0.10４４14 １.00１554 １．008750 -0.089６0８

U2=１.036468-2.639610; U3=1.00８７５0-4.807433 U4=1.007252-5.134133; U5=1.０015５4-5.982486 ６．计算平衡节点功率S1和线路功率Sij

迭代收敛后,就可以计算平衡节点功率和线路功率。结果如下

S1=1.2981６2＋j0.2４４472 各线路功率如表4

表4 各线路功率Sij ﻩj i １ 0．８89５05 0.4086５7 +j0．138６+j0.1058１62 2 -0．8750７9 -j０.09５385 3 -0.39596９ 0 ０.24688４ 0.279３19 0.54８870 1 ２ 3 4 5 ~····~~~＋ｊ０.０+ｊ0.0８+ｊ0.13３28145７ -０．２43 060６ １ 0.１89０79 -j0.012１20 -j0.067７4１08 ７ -j０.07013３ 4 -0．2７4598 -0．1８8726 -j0.0664４5 +j0.013１78 ０．06332５ ＋j0.003２67 5 －0．5３６ 99０ -ｊ-０．0６30 ０8 －j0．002计算完毕。

3.１ 设计流程

0.097684 第三章３1６

程序设计

３．2 程序设计

针对两机五节点网络潮流计算基于MATLAB的程序如下：

Cｌｃ ％清除当前屏幕上显示的内容,但并不清楚工作空间的数据 Clｅar %清除工作空间内的所有变量 disp（'节点总数为:')； %显示内容 N=5

diｓp('平衡节点为:'）; 1

Ｙ=inｐut（'请输入导纳矩阵：Y=＇); pr=ｉｎｐｕt('请输入误差精度:pr=');

G＝real(Y)；Ｂ＝ｉmag（Y）; ％real实部，iｍaｇe虚部 B1=B（２:5, 2：5)； %去掉平衡节点

B2=B1;

b１=inv(B1); %取逆 ｂ２=inv（B２）;

ure=ｉnpuｔ（＇请输入平衡节点1电压初值的实部，Ｒe='); dｉm=input(＇请输入平衡节点１电压初值的虚部,Iｍ='); foｒ m=2：5

dt(m)＝0; %电压虚部初值 u（m)＝1．０； ％电压实部初值 eｎd ｆor m=2:5

p(m)=inpｕt('请输入PQ节点的有功功率：ｐ＝＇); ｑ（ｍ)＝ｉnpuｔ(＇请输入PＱ节点的无功功率：ｑ='); eｎd

k=0;wｕｃｈa＝1; %迭代次数, ｗhｉle ｗｕｃha（1）>pｒ; u（1)＝ure；

dt(1)＝dim； %平衡节点电压实部 fｏr m=2：5 %平衡节点电压虚部 fｏr ｎ=１:5

pｔ(ｎ)＝u(m)＊u(n）*(Ｇ(m,n)＊ｃos(dｔ(m)-ｄt(ｎ))+B(m,ｎ)*siｎ(ｄt(m)-ｄ

0t(ｎ))); %Ui0Uj0GijcosijBijsinij0

j1jnend

disp(＇各节点有功率不平衡量为：')

0dp(m)=p(m)-ｓum(ｐt) ％Pi=pi-Ui0Uj0GijcosijBijsinij0

0j1jnend

ｆｏr m=２：５

u1=ｄｉag（u); ％对角矩阵 ｕ２＝ｉnv（u１(2:5,2：５）); %取逆

u3=ｆｌipｕd(roｔ90（u（2：5）)); %逆时针90°上下对称 u4=ｕ3(１:5-１,：）;

ｄp1＝fliｐｕd(rot90（dp(2：5)))； %—(B')1 （Pｄp2=dp１(１：５－１,:)； ｄtｔ=(-b1*dp2／u4)*u2； ｄtt=ｄｔｔ(:,１:１）； ｄtt=fｌipｕd(ｒot90(ｄtt))

dｉｓp（'各节点电压的相位角dt为:') dｔ(ｍ）=dt(m)＋dtt(m-１) eｎd for m=２：5 for ｎ=1:5

qt(ｎ)=u(ｍ)*u(n)＊(G(m,n)*sin(dｔ（ｍ)-ｄt(ｎ））-B(m，n)*cos(dｔ(ｍ)－dt(ｎ)));

1end %Ui0Uj0GijsinijBijcosij1

j1jnU)=U0

ｄｉｓp('各节点无功率不平衡量为:') q

ｑ

0(m)=

jnj1ｑ(ｍ)-sum（qt) %

1Bijcosij1 Qi=Qi-Ui0Uj0Gijsinijend

fｏr m=２:５

ｑq１=ｆlipｕｄ(rot９0(qｑ（２:5）)); ｑq2=qq1(１:5-1,:)； ut=-b2*qq2/u4;

disｐ( ＇各节点电压的修正量：') ut

dｉｓp( '各节点电压的大小为:'） u（m）=u(m)+ut(m－1)

eｎd ｆｏr n=1:5

U(n)=u(n)*(cos（dｔ（ｎ))+ｊ*siｎ(dt(n))); %电压 end

foｒ ｎ=１：5

Ｉ(n)＝Y(1,n)*U（n)； end

ｄisp（'各节点功率s为:')

Ｓ1＝U(1)＊sum（conj(I)) %电流 foｒ ｍ=1：5 for ｎ=１:5

S(m，ｎ）＝U(m)*（conj(U（m))-coｎj(Ｕ（ｎ）)）*cｏnj（-Y(ｍ，n）); 功率 end end

ｄｉsp('各节点功率sij为：'） S

ｗucha=ｍａx(abs(ut)）

ｋ=k+1 %迭代+１ ｅnd

% ３.3 程序运行结果

节点总数为: N = 5 平衡节点为： aｎs = １ 请输入导纳矩阵:

Ｙ=［6.25-18.75i -5+1５ｉ -１.25+3.７5i 0 0

-５+15i １０.８34-3２.５ｉ -１.667+5i -1．6６7＋5i －2．５＋７.５i

-1.２5＋３.７5i -1.667+5i １2.917-3８.75i -１0+３0i 0

0 -1．667+5ｉ -10+30ｉ 12．91７－38．75i -1.25+3.7５i

0 -２．５+7.5i 0 -１.25+3.75i ３.75-11．25ｉ]

请输入误差精度：pr=0.00001

请输入平衡节点1电压初值的实部，Re＝１.06 请输入平衡节点１电压初值的虚部，Im=0 请输入PＱ节点的有功功率：p=0.2 请输入PQ节点的无功功率：ｑ=0.２ 请输入ＰＱ节点的有功功率：p=－０．45 请输入PQ节点的无功功率：q＝-0.15 请输入PQ节点的有功功率:p=-0．4 请输入PQ节点的无功功率:ｑ=－0.05 请输入ＰQ节点的有功功率：ｐ=-０.6 请输入PQ节点的无功功率:q=-０.1

各节点有功率不平衡量为: dｐ =

０ 0.５000 各节点有功率不平衡量为: ｄｐ ＝

０ 0.５０00 -0.３750 各节点有功率不平衡量为: dp =

０ 0.50０0 -0.3750 -0.４000 各节点有功率不平衡量为: dp ＝

0 0.5000 －0．3750 -0.40０0 -0.６000 dｔt =

-0．03７0 －0.08５5 -0.0918 -0.1086 各节点电压的相位角dｔ为： dt =

０ -０.0370 0 0 ｄtｔ =

－0.0370 －０.０855 -0.０918 -0.1086 各节点电压的相位角dt为： dｔ =

0 －０．0３70 －0.0８５5 ０ dtt =

-0.0370 －0.0855 －0.0918 -0.108６ 各节点电压的相位角ｄｔ为： ｄt =

0 -0.０37０ －0.0８55 －0.0９18 dｔt =

0 0 ０ -0.0３70 －0.0８５5 -0.０918 －0.1０86 各节点电压的相位角dt为： dｔ ＝

0 -０.０３70 -0.0８5５ -0.09１8 -0.1086

各节点无功率不平衡量为: qq =

0 １.2119 各节点无功率不平衡量为： qq =

０ 1.2119 -0.0773 各节点无功率不平衡量为: qq =

0 1.21１9 -0.0７73 -0.19１９ 各节点无功率不平衡量为: qq =

0 1．21１9 -0.０７7３ -0.19１9 -0.31９6 各节点电压的修正量: ut =

0.0395 0 0 0 0.008１ 0 ０ 0 0.０064 0 0 0 0.0001 0 0 0 各节点电压的大小为： u ＝

1.0600 1.０3９5 １.０00０ １．0000 1．000０ 各节点电压的修正量: ut ＝

0.0395 0 0 0

0.０081 0 0 0 ０.0064 ０ 0 0 0.０001 0 ０ 0 各节点电压的大小为: u =

１.0６0０ 1．０395 1.０081 １.0000 1.0００0 各节点电压的修正量: ut =

0．03９5 0 0 0 0.0081 0 0 0 ０.０0６4 ０ ０ ０ 0.0001 0 0 0 各节点电压的大小为： u ＝

1.060０ 1.0３95 １.0０81 1.0064 1．0000 各节点电压的修正量： ut ＝

０．0３９５ ０ 0 0 ０．００81 ０ ０ 0 ０.0064 0 ０ 0 0.０00１ 0 0 0 各节点电压的大小为: u =

１.0600 １.0395 １.0081 1．0064 1.00０１ 各节点功率s为: S1 ＝

1.13８9 + 0.24０３i 各节点功率sij为: S =

0 0.722９ + ０.１332ｉ 0．４16０ + ０.1071i 0 ０ -０.７1３3 - 0.10４４i 0 0.3１10 ＋ ０.085０i 0.3469 + 0.０８4５i 0.６671 + ０.1416i

-０.４０29 - 0.0６７6i －０.3052 － 0.0677i 0 0.2083 - 0.0128i 0

0 -０.3398 － 0．０633ｉ -0．2０7８ ＋ ０.０141i ０ 0．0714 + 0.0０３3ｉ

0 -0.6４99 - 0．0900ｉ -0.0710 － ０．002１i 0 wucha =

0．03９5 ０ ０ ０ k = １

各节点有功率不平衡量为: ｄｐ =

0 -0.4１17 -0.375０ -0.４000 -0.6000 各节点有功率不平衡量为： dp =

0 -0.41１7 0．0４９8 -０.４000 -０．６000 各节点有功率不平衡量为: dp =

0 -0.411７ ０.0４98 0.０762 -0．600０ 各节点有功率不平衡量为: dp =

0 －0．4１17 0.049８ ０.０76２ 0．1２09 ｄtt =

-0.０105 0.００1２ 0.０017 0．０0３5 各节点电压的相位角ｄｔ为：

0 dt =

0 -0.０４7４ -0.0８5５ -0．091８ -0.1０86 dtｔ =

－0.0105 0.0012 ０．00１７ 0.0０35 各节点电压的相位角ｄt为： ｄｔ =

0 -0.0474 －0.0843 -０．0９1８ -0.１086 ｄｔt ＝

－0.０10５ 0.0012 0.0017 0.003５ 各节点电压的相位角dt为： dｔ =

０ -０.０４7４ -0.0843 －0.０９01 －0.10８6 ｄtt =

-0.0105 0．００1２ ０．00１7 0.0035 各节点电压的相位角dt为: ｄt =

０ -0．0４7４ －0.０８4３ -0.０9０1 －０.10５0 各节点无功率不平衡量为: qq =

0 -０.1３67 -0.0773 -0.1919 -0.3196 各节点无功率不平衡量为： qq =

０ －0．1３67 0.017５ -0．１919 -0.319６ 各节点无功率不平衡量为： qｑ =

0 -０．1３67 0．01７5 ０.0237 -0.31９６ 各节点无功率不平衡量为: qｑ =

0 －０.１3６7 0.0１75 0.０２37 ０.0３７8

各节点电压的修正量： ut ＝

－0.0038 0 0 ０ ０.０0０2 ０ 0 0 0.000４ 0 0 0 0．０0０9 0 0 0 各节点电压的大小为： u =

1.060０ １.０3５8 各节点电压的修正量： ut =

-0.00３8 0 ０.0０0２ 0 0.00０4 0 0.0０09 0 各节点电压的大小为： u ＝

1.0600 １.０35８ 各节点电压的修正量: ut =

-0.0038 0 ０．0０0２ 0 0.0００4 ０ 0.0009 0 各节点电压的大小为: u =

1.06０0 1.０3５8 各节点电压的修正量: ut =

1．0081 0 0 0 0 １.0０8３ ０ ０ 0 0 1．0０8３ 1.０064 0 ０ 0 0 1.006４ 0 ０ 0 1.００68 1.0０01 1.０00１0 1.0００1 －０.0０３８ 0 0 0 0.000２ 0 ０ 0 0.0００４ 0 0 0 0.0009 0 0 0 各节点电压的大小为: ｕ ＝

1．０６0０ 1.0３58 1.００83 1.0068 1.０009 各节点功率s为： S1 =

1.3２60 + ０．250８i 各节点功率ｓiｊ为： S ＝

0 0.９15１ + 0.143５ｉ 0．4109 + 0.1０7２ｉ ０ 0 -０.89９8 - 0.0977i 0.２4１5 + 0.0８16ｉ ０.2741 + 0.０8０9ｉ ０.542３ + ０.1３4４ｉ

-0.3９８１ － 0.０６8７i －0.2３78 － 0．0７０7i 0 ０.1903 - 0.0１14ｉ 0

0 －0.2695 － 0．０6７2i -０．18９9 + 0.０12５i 0 0.0６３9 ＋ 0.0０3６i

0 －0.5３06 - 0.０995i ０ -０．０6３6 - 0．０02６i 0 wucha =

0.０038 0 0 ０ ｋ = 2

各节点有功率不平衡量为： ｄp =

0 0．0419 0.0４98 0.07６2 0.１209 各节点有功率不平衡量为：

ｄｐ ＝

０ 0.04１９ -0.0０４4 0．076２ 0.1209 各节点有功率不平衡量为: dｐ ＝

0 ０.0419 －0．００４４ －０．004４ 0.1209 各节点有功率不平衡量为: dp =

0 0.041９ -0．004４ -0.0044 -0.0058 dtt ＝

0.0016 0.0006 0.00０6 0.0008 各节点电压的相位角ｄt为: ｄt ＝

０ -0．0459 -0.0843 -0.0901 -0.10５0 dtt =

0.00１６ ０.0006 0．0０06 ０.0008 各节点电压的相位角ｄt为： dｔ ＝

0 －0.0459 －0.０838 -０.0901 －0．１05０ dtt =

0.001６ 0．0０0６ ０.0０0６ 0.0０0８ 各节点电压的相位角dt为: dｔ ＝

0 -0.０459 -０.08３8 －0．0895 －0.1050 dtt =

0.0016 0.00０6 0.０0０6 0.000８ 各节点电压的相位角dt为: ｄｔ ＝

0 -0.0４5９ -０.0838 -0.0895 －0.10４3 各节点无功率不平衡量为：

qq =

０ 0.0153 ０．0175 ０.0237 0.037８ 各节点无功率不平衡量为: qq ＝

0 0．01５3 -0．00０5 ０.０237 0.0378 各节点无功率不平衡量为： qq =

0 0.01５３ -0.00０５ －0.000６ 0．0３７８ 各节点无功率不平衡量为： ｑq =

0 0.0153 －０.000５ -0.0006 －0.０001 各节点电压的修正量： uｔ =

1.０e-0０3 *

0.7７76 ０ 0 0 0.51１8 ０ 0 0 0．５４9１ 0 ０ 0 ０.６921 0 0 ０ 各节点电压的大小为: ｕ =

1．０600 1.０366 1.00８３ 1.０06８ 1.００0９ 各节点电压的修正量： ut =

1.0e-０03 *

0．７776 0 0 0 0.5１18 0 ０ 0 0.5491 ０ 0 ０ 0.6921 0 0 ０ 各节点电压的大小为：

u =

1.０600 1.０366 1.0０８8 １.0068 1．０0０9 各节点电压的修正量： uｔ =

1.０e-００3 *

0.７７76 0 0 0 0.5118 0 ０ ０ 0.5４91 0 ０ ０ 0.6921 0 0 ０ 各节点电压的大小为: u =

1．０600 1.03６6 1．０088 1.0073 1.0009 各节点电压的修正量: ut =

1．0ｅ-003 *

0.7７７6 0 0 0 0.５１18 0 0 0 0.５491 0 0 0 0．6９21 ０ 0 0 各节点电压的大小为: ｕ ＝

１.06０0 １.０366 1.00８8 1．0073 1.0016 各节点功率ｓ为： Ｓ1 =

１.２９36 + 0.2441i 各节点功率sij为： S =

0 0.８8５5 + 0.1３83ｉ 0.４０81 + 0．１057ｉ 0 -0.８71２ - ０.0955i 0 0.2475 + 0.0８15ｉ 0 0.27９

9 + 0.0８06i ０.５４95 + 0.1３33ｉ

-0.3954 － ０.０678i -0．2437 － 0.０701i 0 0.1889 - 0.０１22i 0

0 －０.２75１ － ０.0664ｉ －０.1８85 + ０．０１32i 0 ０.063２ + 0.００３3i

0 -0.5376 - 0.0９７６i 0 －０.0629 - 0.0０23i 0 ｗucha = 1.0e-00３ *

０.7776 ０ 0 0 k ＝ 3

各节点有功率不平衡量为： dｐ ＝

0 －０.００57 -0.０04４ -０.０044 -０.0058 各节点有功率不平衡量为: dｐ =

0 -０.005７ ０.000３ -0.0０4４ －0．００58 各节点有功率不平衡量为: dp =

0 -0.0０57 0.００03 0．0００４ －0.０05８ 各节点有功率不平衡量为: dp =

0 -０.00５7 0.０003 0.0０04 0.00０５ dtt =

1.0ｅ－0０３ *

－０.２４42 －0.1346 －0.141３ -0.1６74 各节点电压的相位角ｄt为: dｔ =

0 -0.0461 －0.083８ －0.0８95 -０.１043 dｔｔ ＝ 1.0e-003 *

-0.244２ -０.1346 -0.1413 -0.1674 各节点电压的相位角dt为: dt =

0 －０．0４６1 -0．083９ －0.089５ -0.1043 dtt =

１.0ｅ－003 ＊

-0.２442 －0.1３46 -0.14１3 -0.167４ 各节点电压的相位角dt为: dt =

0 -0.0461 -0.0839 -0.0８9６ -0.1043 dｔt ＝ 1.0e-00３ *

-0.2４42 -0.１346 －0.１413 -0.１６7４ 各节点电压的相位角dt为: dt =

０ -0.0461 -０.０839 -0.0８96 -0.1044 各节点无功率不平衡量为: qｑ =

0 -0.0０20 -0.０0０5 -０．０006 -0．00０1 各节点无功率不平衡量为： qq =

０ -0．0０20 0.００0１ -0.0006 -0.000１ 各节点无功率不平衡量为： qq ＝

0 －0.0020 ０.０001 ０.0０01 -0．00０1

各节点无功率不平衡量为: qq ＝

0 -0.002０ 0.000１ 0.0００1 ０.00０1 各节点电压的修正量: uｔ =

１.0e－０04 ＊

-0.9362 ０ -0．5617 0 －0.６0０2 0 -0．7１1９ ０ 各节点电压的大小为： u =

1．060０ １．03６５ 各节点电压的修正量： ut =

1.0e-004 *

-0.9３6２ 0 －0.５６17 ０ -０.600２ 0 -0.７１19 0 各节点电压的大小为： u =

１.0600 1.03６5 1.0087 各节点电压的修正量： uｔ = 1.0ｅ-０04 *

-0．９３６2 0 －0.5617 0 －0.6002 0 ０ 0 0 0 0 0

0 0 .0088 1.0073 1．０0１6 0 0 0 0 0 0 ０ 0 1.0０7３ １.0０１６ 0 0 0 0 0 0

１ -0.７119 ０ 0 0 各节点电压的大小为: u =

1.0600 1.0365 １.０087 1.0072 １.0０16 各节点电压的修正量： ut = 1.0e-004 *

-0.9362 0 0 0 －0.5617 0 0 0 －0.6002 ０ 0 0 -0.７119 0 ０ 0 各节点电压的大小为: u ＝

1.0600 １.0365 １.0087 1.０0７2 １.0０1５ 各节点功率s为: S１ =

１.2987 + 0.２445ｉ 各节点功率sｉj为： S =

０ 0．8900 + ０.1３８7i 0.4087 + 0.1０５8i 0 0 －0．８75５ - 0．０954i 0 0．246８ + 0．08１5i 0.２7９２ ＋ ０．08０6ｉ 0.５488 + ０.13３3ｉ

-0.396０ - 0．0677i -0．243０ - ０.0701ｉ 0 ０．1891 - 0.0121i 0

0 －０.2745 - 0.0６６5i －0．1887 + 0.0132i ０ 0．0６３3 + 0.003３i

0 -０．5369 - 0.097７ｉ 0 -0．06３０ － 0.002３i ０ ｗuchａ =

1．0ｅ-00４ *

0.９3６2 0 ０ 0 k = 4

各节点有功率不平衡量为: dｐ ＝

１.0ｅ-00３ *

0 0.7371 0.２61８ 0．4425 0.5138 各节点有功率不平衡量为: dp =

1.0e-00３ *

０ 0.7３７1 －0.0４82 0.４４2５ 0．5138 各节点有功率不平衡量为: dｐ ＝

1．0e－００３ *

0 0．７371 －0.04８2 －0．0６39 0．5138 各节点有功率不平衡量为： dp =

1．0e-0０3 ＊

0 0．７371 －0.048２ -0．０639 -0．1022 dtt =

1.0ｅ-0０4 *

0．2903 0.13６6 ０.１431 0.1５67 各节点电压的相位角ｄt为： dt =

0 -0.０461 -０.0839 -0.08９6 -0.1044 dtt =

1.０e-００4 ＊

0.２９0３ 0.1３66 0.14３1 0.1567 各节点电压的相位角dt为: dｔ =

０ -０.０461 -0．０839 -0.089６ -0.1044 dtｔ = 1.0e-0０4 *

0.2903 ０．1３66 ０.1431 0.１567 各节点电压的相位角ｄt为： dt =

0 -0.0461 －０．０8３９ -0.０８９6 -0.1０４4 ｄtt ＝ 1.０e-0０4 *

0．２９０3 0．1366 0.1431 ０.1567 各节点电压的相位角dt为： ｄt =

0 -0.0461 -0．０839 -0．0８96 －0.10４4 各节点无功率不平衡量为： qｑ = 1.0e-003 *

0 0．2564 0．０953 0.０９８０ 0.1310 各节点无功率不平衡量为: qq =

1.0ｅ-00３ *

0 ０．2564 －0．0１6４ 0．０980 ０.1310 各节点无功率不平衡量为: ｑq =

1.０e-00３ *

0 0.２564 -0.0１64 -0.0149 0.131０ 各节点无功率不平衡量为:

qq ＝

１.０ｅ-003 *

0 0．２５64 -0．0１６４ －０.01４９ -0.0256 各节点电压的修正量： ut =

１.0ｅ-004 *

０.1１30 ０ 0 ０ 0.0613 0 0 0 0.0656 0 0 ０ 0.0753 0 0 0 各节点电压的大小为: u =

1.０６00 １.03６5 1．0０87 1．0072 1.00１５ 各节点电压的修正量： ut =

１.0ｅ-００4 *

０．1130 0 0 0 0．０613 0 0 ０ 0．０６56 ０ 0 ０ 0.０7５３ 0 0 0 各节点电压的大小为: u =

１.0600 １.036５ 1．００88 1．0072 1.００1５ 各节点电压的修正量: ut =

1．０ｅ-０04 *

０．1１３0 0 ０ 0 ０.０６13 0 0 0 0.０656 0 ０ 0

0.0753 0 0 ０ 各节点电压的大小为: u =

１.０600 1.０３65 1.0０88 1.0073 1.0０15 各节点电压的修正量: ｕt =

1．０e－０04 *

0.１１30 0 0 0 ０.0６１3 ０ 0 ０ 0．0６56 ０ 0 0 0．0753 0 0 ０ 各节点电压的大小为： u =

1.0600 1.0365 １．00８８ 1．0０73 1．0０１6 各节点功率s为: S1 =

1.２981 + 0．24４5i 各节点功率sij为： S ＝

0 0.88９４ + 0.138７ｉ 0．4０８7 ＋ ０.１０５8i 0 0 -０.8７５０ － ０．0954ｉ 0 0.2469 ＋ ０.0815i 0.2７９3 + 0．0806i 0.548９ + 0.1333ｉ

－0.3９６0 - ０．0677i -0.2431 - ０.0701i ０ ０.1891 - 0.0１21i 0

0 -０.27４6 － 0.０664ｉ －０.18８７ ＋ 0.0132ｉ 0 ０.0633 + 0．0033ｉ

0 -0.5370 - ０.０9７7i ０ -0.06３0 - 0．0０2３i 0

wucｈa = １.0e－0０4 *

0．1130 0 0 0 ｋ = 5

各节点有功率不平衡量为: dp =

１.0e－003 *

0 -0.０9２4 -0.０482 -0.0６39 -０．1022 各节点有功率不平衡量为： dｐ ＝ 1.０ｅ-003 *

0 -０．０９2４ 0.0０6７ －０.06３9 -0.102２ 各节点有功率不平衡量为: dp =

1.0e-0０3 *

０ -0.0924 0.006７ 0.00８0 -0.1022 各节点有功率不平衡量为： dp =

1.0ｅ-004 *

0 -0.923７ 0.0667 0.0７97 ０.125５ dtｔ ＝ 1.0e-0０５ *

-０.3６27 -0.1670 －0.1760 －０.1966 各节点电压的相位角dt为： dt =

０ －0.０4６1 －０．0839 -0.089６ －0.1044 dtｔ = 1.0ｅ-005 *

-0.３６2７ －0.1６7０ -0．1７６0 -0.1９6６ 各节点电压的相位角dt为： dt =

０ -０.0４6１ -0.08３９ -0．０896 -0.1０44 ｄｔｔ = 1.0e-00５ *

-０.３６27 -0.16７0 -0.1７60 -０.1９66 各节点电压的相位角dｔ为： ｄt ＝

0 -０.046１ -０.0839 -0.089６ －0．1044 dtt ＝

1．０e-005 ＊

-0.3627 －0.1670 -０.1760 -0.1966 各节点电压的相位角dt为: dt =

0 －0.04６１ －0．０８39 -0.0896 -0．1０４4 各节点无功率不平衡量为： qq =

1.０e-004 *

0 -0.323５ -0.16３６ -0．149２ -０.2558、 各节点无功率不平衡量为： qｑ ＝ 1.0e-004 *

0 -0．3２35 0.0199 -0．1492 -0.2５5８ 各节点无功率不平衡量为： qq ＝ 1.0e-004 *

0 -0．32３5 0.0199 0.0180 -0.2558 各节点无功率不平衡量为:

qq =

1．０e-0０4 *

０ －０.3２35 0.019９ 0.0１80 0.0276 各节点电压的修正量： uｔ = 1.0e-00５ ＊

-0.１４56 ０ ０ 0 -0．0813 0 0 ０ -0.0872 0 0 0 -0.１0２5 0 ０ 0 各节点电压的大小为： u =

1．06００ 1.0365 １.00８8 １.00７3 1．００1６ 各节点电压的修正量： ut =

１.０e-005 *

-０．１4５6 0 0 0 -0．0８13 0 0 0 －０．0８７2 0 0 0 －0.10２5 ０ 0 0 各节点电压的大小为: u =

1.0600 １.0365 １.00８7 １.00７3 １.０016 各节点电压的修正量： ｕt =

1．０e-005 *

-0.1456 ０ ０ 0 -０.0８１３ 0 0 ０ -0.08７2 0 0 0

-０.10２５ 0 0 0 各节点电压的大小为： u =

1.06０0 １.０3６５ 1．0087 1．0０73 1．0０16 各节点电压的修正量： ut =

1.0e－005 *

-0.14５6 0 0 ０ -0．0813 0 0 0 -0．0872 0 0 ０ -０.10２5 ０ ０ 0 各节点电压的大小为: u =

１.0６0０ 1.036５ １.００８7 １.０073 １.0016 各节点功率s为: S1 =

1．2982 + 0．2445i 各节点功率siｊ为： S =

0 ０.８895 + 0.1387i 0.40８7 + 0.1０58i ０ 0

-0.8751 - 0.095４ｉ 0 0．2４６9 + 0.0815i 0.2７93 + 0.0８０6i

0．5489 + 0．１3３3ｉ －0.39６0 - 0.０６７7i -0.24３1 - ０.０7０1ｉ 0 0．1891 - 0．0121i 0

0 －0.274６ - ０.0６64i -0.18８７ + 0.013２ｉ 0 0.0633 + 0.00３３i

0 -０.5３70 - 0．09７7i ０ －0.０６3０ - 0.0023i 0

wuｃha = １.０ｅ-0０5 *

0.1４5６ 0 ０ 0 k = 6 ＞>

第四章 设计感想

疯狂忙碌了近一周后，课程设计也进入了尾声，着实有很多不舍。

首先,此次课程设计，是我们大学以来第一次严格进行的课程设计,这是我们进入社会前的一次小小的磨练,为我们今后参加工作打下一点基础。所以收获还是蛮多的。

一周内,我们初步的掌握了对MAＴLAB的应用。这是一个很优秀的软件，在生活,学习,航天,军工等方面都有很大的用处。

其次，我们对稳态潮流分析,掌握的更加深刻，不单单可以进行简单的计算，还可以凭借计算机,对各种复杂网络进行潮流计算。

无论是收获到的知识，还是技能,我觉得都比不上这次课程设计给我们带来的态度上的改变。俗话说得好，不逼自己，你永远不知道自己有多强。当拿到课程设计任务书的时候,觉得很不可思议：要用自己从来没接触过的一门计算机语言,去独立完成一个高难度的计算。军令如山。一周内,我们从寝室到图书馆，从食堂到课设教室,过着几点一线的生活。早上８点多到教室，12点下课了还浑然不知,很多次,晚上10点半时都是门卫打着手电过来催我们走。晚上回去睡觉还在想哪做的不对，为什么不对,第二天早起就去图书馆疯补不会的知识。 什么是幸福?忙碌着，收获着，简单着,快乐着,这就是幸福。 感谢这次课设计带给我的幸福。 最后感谢老师对我的悉心教导。

生无所息,奋斗不止。

参考文献

1、陈珩.电力系统稳态分析.北京；中国电力出版社,2００7

2、原思聪．MAＴLAＢ语言与应用技术，北京;国防工业出版社2０11 3、陈泽，占海明．详解ＭＡTLAB在科学计算中的应用，北京;电子工业出版社,2０11

４、祝淑萍.电力系统分析课程设计与综合试验.北京：中国电力出版社，2007