电力系统稳态分析牛顿拉夫逊法
- 格式:doc
- 大小:999.00 KB
- 文档页数:31
摘要本文,首先简单介绍了基于在MALAB中行潮流计算的原理、意义,然后用具体的实例,简单介绍了如何利用MALAB去进行电力系统中的潮流计算。
众所周知,电力系统潮流计算是研究电力系统稳态运行情况的一种计算,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态:各线的电压、各元件中流过的功率、系统的功率损耗等等。
在电力系统规划的设计和现有电力系统运行方式的研究中,都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性、可靠性和经济性。
此外,在进行电力系统静态及暂态稳定计算时,要利用潮流计算的结果作为其计算的基础;一些故障分析以及优化计算也需要有相应的潮流计算作配合;潮流计算往往成为上述计算程序的一个重要组成部分。
以上这些,主要是在系统规划设计及运行方式安排中的应用,属于离线计算范畴。
牛顿-拉夫逊法在电力系统潮流计算的常用算法之一,它收敛性好,迭代次数少。
本文介绍了电力系统潮流计算机辅助分析的基本知识及潮流计算牛顿-拉夫逊法,最后介绍了利用MTALAB程序运行的结果。
关键词:电力系统潮流计算,牛顿-拉夫逊法,MATLABABSTRACTThis article first introduces the flow calculation based on the principle of MALAB Bank of China, meaning, and then use specific examples, a brief introduction, how to use MALAB to the flow calculation in power systems.As we all know, is the study of power flow calculation of power system steady-state operation of a calculation, which according to the given operating conditions and system wiring the entire power system to determine the operational status of each part: the bus voltage flowing through the components power, system power loss and so on. In power system planning power system design and operation mode of the current study, are required to quantitatively calculated using the trend analysis and comparison of the program or run mode power supply reasonable, reliability and economy.In addition, during the power system static and transient stability calculation, the results of calculation to take advantage of the trend as its basis of calculation; number of fault analysis and optimization also requires a corresponding flow calculation for cooperation; power flow calculation program often become the an important part. These, mainly in the way of system design and operation arrangements in the application areas are off-line calculation.Newton - Raphson power flow calculation in power system is one commonly used method, it is good convergence of the iteration number of small, introduce the trend of computer-aided power system analysis of the basic knowledge and power flow Newton - Raphson method, introduced by the last matlab run results.Keywords:power system flow calculation, Newton – Raphson method, matlab目录1 绪论 (1)1.1 课题背景 (1)1.2 电力系统潮流计算的意义 (2)1.3 电力系统潮流计算的发展 (2)1.4 潮流计算的发展趋势 (4)2 潮流计算的数学模型 (5)2.1 电力线路的数学模型及其应用 (5)2.2 等值双绕组变压器模型及其应用 (6)2.3 电力网络的数学模型 (8)2.4 节点导纳矩阵 (9)2.4.1 节点导纳矩阵的形成 (9)2.4.2 节点导纳矩阵的修改 (10)2.5 潮流计算节点的类型 (11)2.6 节点功率方程 (12)2·7 潮流计算的约束条件 (13)3 牛顿-拉夫逊法潮流计算基本原理 (14)3.1 牛顿-拉夫逊法的基本原理 (14)3.2 牛顿-拉夫逊法潮流计算的修正方程 (17)3.3 潮流计算的基本特点 (20)3.4 节点功率方程 (21)4牛顿-拉夫逊法分解潮流程序 (22)4·1 牛顿-拉夫逊法分解潮流程序原理总框图 (22)4.2 形成节点导纳矩阵程序框图及代码 (23)4.2。
牛顿拉夫逊法计算潮流步骤牛顿拉夫逊法(Newton-Raphson Method)是一种常用于计算潮流的数值求解方法。
它是基于潮流计算的功率流方程的非线性特性而设计的,通过迭代求解来逼近潮流计算的稳态解。
下面将介绍牛顿拉夫逊法计算潮流的基本步骤。
首先,我们需要明确潮流计算的目标,即确定电力系统中各节点的电压相角和幅值。
这些节点是电力系统中的发电机、负荷和交流输电线路的连接点。
通过潮流计算,我们可以得到各节点的电压相角和幅值,从而分析系统的功率分布、电压稳定性等运行特性。
接下来,我们需要建立电力系统的潮流计算模型。
这个模型中,我们需要考虑发电机的注入功率、负荷的吸收功率、线路的传输损耗等因素。
通过利用功率流方程,我们可以将这些因素表示为电压、功率和导纳之间的方程。
然后,我们需要进行初始化操作。
在进行牛顿拉夫逊法迭代计算之前,我们需要对电力系统的各节点进行初始电压值的设定。
这些初始值可以根据经验或者历史数据来得到,但需要满足物理约束条件,如一致性、电压幅值在合理范围内等。
接下来,我们进入迭代计算的过程。
首先,我们需要对系统的节点进行编号,然后选择某一节点作为基准节点,其他节点相对于基准节点的电压相角进行计算。
然后,我们根据节点注入功率和导纳矩阵的关系,得到节点注入电流。
接着,我们根据节点注入电流和电压相角的关系,计算各节点的电压相角和幅值的改变量。
在计算改变量后,我们需要对节点电压进行更新。
更新后,我们判断系统是否达到收敛条件。
如果满足收敛条件,则停止迭代,得到最终的潮流计算结果;如果不满足收敛条件,则继续进行下一轮迭代计算。
最后,我们对潮流计算结果进行分析和验证。
通过比较计算得到的结果和实际运行数据进行对比,我们可以评估潮流计算的准确性。
同时,我们还可以通过故障分析、电压稳定性评估等手段对电力系统进行优化和改进。
总而言之,牛顿拉夫逊法是一种常用的求解潮流计算问题的方法。
它通过迭代求解潮流计算的功率流方程,逼近潮流计算的稳态解。
电力系统牛顿拉夫逊法频率电力系统中的牛顿拉夫逊法是一种常用的频率分析方法,用于研究电力系统在不同频率下的稳定性和响应特性。
该方法基于牛顿第二定律和拉夫逊法则,结合电力系统的动态方程,可以求解系统的频率响应和稳定性。
在电力系统中,频率是指电压和电流信号中的周期性变化。
电力系统中的频率通常为50Hz或60Hz,代表每秒钟电压和电流信号的周期数。
频率的稳定性对于电力系统的正常运行至关重要,因为频率的偏离会导致电力设备的失效和电网的不稳定。
牛顿拉夫逊法是一种通过求解系统的微分方程来研究系统动态行为的方法。
在电力系统中,牛顿拉夫逊法可以用于求解系统的振荡频率和振荡模态。
通过分析系统的振荡频率和振荡模态,可以评估系统的稳定性和抗干扰能力。
牛顿拉夫逊法的基本思想是将电力系统的动态方程转化为一组一阶微分方程组,然后通过数值求解方法求解这组微分方程。
这组微分方程描述了电力系统中各个节点的电压和相角随时间的变化规律。
通过求解这组微分方程,可以得到系统在不同频率下的响应特性。
在应用牛顿拉夫逊法进行频率分析时,首先需要建立电力系统的动态模型。
这个模型包括系统的节点、传输线、发电机、负荷以及其他与系统运行相关的参数。
根据这个模型,可以得到系统的微分方程组。
然后,需要选择适当的求解方法来求解微分方程组。
常用的求解方法包括欧拉法、龙格-库塔法和改进的欧拉法等。
这些方法可以通过迭代计算的方式得到系统在不同时刻的电压和相角值。
在求解微分方程组之后,可以通过对系统的电压和相角进行频谱分析来得到系统的频率响应。
频谱分析是一种将信号分解成不同频率分量的方法,通过计算每个频率分量的幅值和相位,可以得到系统在不同频率下的响应特性。
通过牛顿拉夫逊法进行频率分析可以帮助电力系统工程师评估系统的稳定性和可靠性。
通过分析系统的振荡频率和振荡模态,可以发现系统中存在的潜在问题,并采取相应的措施来提高系统的稳定性。
牛顿拉夫逊法是一种常用的电力系统频率分析方法,通过求解系统的微分方程组,可以得到系统在不同频率下的响应特性。
电力系统稳态分析摘要电力系统潮流计算是研究电力系统稳态运行情况的一种重要的分析计算,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态:各母线的电压,各元件中流过的功率,系统的功率损耗。
所以,电力系统潮流计算是进行电力系统故障计算,继电保护整定,安全分析的必要工具。
本文介绍了基于MATLAB软件的牛顿—拉夫逊法和P—Q分解法潮流计算的程序,该程序用于计算中小型电力网络的潮流。
在本文中,采用的是一个5节点的算例进行分析,并对仿真结果进行比较,算例的结果验证了程序的正确性和迭代法的有效性。
关键词:电力系统潮流计算;MATLAB;牛顿—拉夫逊法;P-Q分解法;目次1 绪论 01.1背景及意义 01.2相关理论 01。
3本文的主要工作 (1)2 潮流计算的基本理论 (2)2。
1节点的分类 (2)2。
2基本功率方程式(极坐标下) (2)2.3本章小结 (3)3 潮流计算的两种算法 (4)3。
1牛顿—拉夫逊算法 (4)3.2PQ分解算法 (10)3。
3本章小结 (14)4 算例 (15)4.1系统模型 (15)4.2结果分析 (15)4。
3本章小结 (18)结论 (19)参考文献 (20)附录 (21)1 绪论1。
1背景及意义电力系统稳态分析是研究电力系统运行和规划方案最重要和最基本的手段。
电力系统稳态分析根据给定的发电运行方式和系统接线方式来确定系统的稳态运行状态,其中潮流计算针对电力系统的各种正常的运行方式进行稳态分析.潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件,确定电力系统各部分稳态运行状态参数的计算.通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。
待求的运行状态参量包括电网各母线节点的电压幅值和相角,以及各支路的功率分布、网络的功率损耗等.电力系统潮流计算问题在数学上是一组多元非线性方程式求解问题,其解法都离不开迭代.潮流计算方法的改进过程中,经历了高斯-赛德尔迭代法、阻抗法、分块阻抗法、牛顿-拉夫逊法、改进牛顿法、P—Q分解法等。
电力系统稳态分析课程设计电力系统稳态分析课程设计题目:两机五节点网络潮流运算—牛拉法姓名:朱润民学号:1167130230学院:信息工程学院专业:电气工程及其自动化班级:11级电气2班指导教师:刘景霞名目摘要潮流运算的目的在于:确定电力系统的运行方式;检查系统中各元件是否过电压或过载;为电力系统继电爱护的整定提供依据;为电力系统的稳固运算提供初值,为电力系统规划和经济运行提供分析基础。
牛顿迭代法(Newton's method)又称为牛顿-拉夫逊方法(Newton-Raphson method),牛顿--拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方法。
其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程。
MATLAB 是一种交互式、面向对象的程序设计语言,广泛应用于工业界与学术界,要紧用于矩阵运算.采纳迭代法,通过建立矩阵的修正方程来依次迭代,逐步靠近真值来运算出电力网的电压,功率分布。
采纳迭代法,通过建立矩阵的修正方程来依次迭代,逐步靠近真值来运算出电力网的电压,功率分布。
本文采纳牛顿-拉夫逊法解算电力稳态潮流,用手算和运算机算法对其进行设计。
关键词:电力系统潮流运算;牛顿—拉夫逊法潮流运算;程序;ABSTRACTThe Power Flow computation's goal lies in: Definite electrical power system's movement way; In checkout system various parts whether overvoltage or overload; Provides the basis for the electrical power system relay protection's installation; Provides the starting value for electrical power system's stable computation, is the electrical power system plan and the economical movement provides the analysis foundation.The Newton iteration method (Newton's method) is called Newton - Rough to abdicate the method (Newton - Rough method), Newton--Rough abdicates the law (i.e. Newton law) is solves the misalignment algebraic equation in mathematics the efficacious device. Its main point is turns the misalignment equation solution process carries on repeatedly to the corresponding linear equation the solution the process. MATLAB is one kind interactive, the object-oriented programming language, widely applies in the industrial world and the academic circle, mainly uses in the matrix operation. Uses the repetitive process, iterates in turn through the establishment matrix's modified equation ,approaches the true value to calculate electric power network's voltage gradually, the power distribution.key word:Electrical power flow computation;Newton - Rough abdicates the law tidal current computation;Procedurekey word:Electrical power flow computation;Newton - Rough abdicates the law tidal current computation;Procedure内蒙古科技大学课程设计任务书12系统接线图其中节点1为平稳节点,节点2、3、4、5为PQ节点。
电力系统的稳态与暂态分析方法稳态和暂态是电力系统分析中两个重要的概念。
稳态分析主要用于评估电力系统在正常运行情况下的性能和稳定性,而暂态分析则关注电力系统在发生故障或其他异常情况下的响应和恢复过程。
本文将介绍电力系统中的稳态与暂态分析方法,并探讨其在电力系统规划、运行和故障处理中的应用。
一、稳态分析方法稳态是指电力系统在正常运行情况下,各电压、电流和功率等参数保持在稳定状态的能力。
稳态分析主要涉及电压、功率、功率因数等参数的计算和评估。
常用的稳态分析方法包括潮流计算、负荷流计算、电压稳定性评估等。
1. 潮流计算潮流计算是稳态分析中最基础的方法之一,用于计算电力系统中各节点的电压、电流和功率等参数。
通过潮流计算,可以确定电力系统中各节点的电压稳定程度,评估传输能力和合理分配负载等。
常用的潮流计算方法包括高斯-赛德尔法、牛顿-拉夫逊法等。
2. 负荷流计算负荷流计算是潮流计算的一种特殊形式,用于分析电力系统中负载的分布和负载对系统潮流的影响。
负荷流计算可以帮助确定合理的负载分配方案,提高系统的稳定性和经济性。
3. 电压稳定性评估电压稳定性是一个评估电力系统稳定性的重要指标,特别是在大规模电力系统中。
电压稳定性评估主要通过计算稳态电压变化范围和电压裕度等参数来判断系统的电压稳定性,并采取相应的调整措施。
二、暂态分析方法暂态是指电力系统在出现故障或其他异常情况下,系统中各参数发生瞬时变化并逐渐恢复到正常状态的过程。
暂态分析主要关注电力系统在故障发生后的动态响应和恢复。
常用的暂态分析方法包括短路分析、稳定性分析和电磁暂态分析等。
1. 短路分析短路分析主要用于分析电力系统中发生短路故障时的电流和电压等参数的变化。
通过短路分析,可以确定故障点、故障类型和故障电流等信息,为故障处理和保护设备的选择提供依据。
2. 稳定性分析稳定性分析是评估电力系统在故障发生后是否能够保持稳定运行的一项重要工作。
稳定性分析主要关注系统的动态行为和振荡特性,通过模拟故障后系统的响应来判断系统的稳定性和选择合适的控制策略。
牛顿拉夫逊法计算潮流步骤牛顿拉夫逊法(Newton-Raphson method)是一种用于求解非线性方程组的迭代方法,它可以用来计算电力系统潮流的解。
潮流计算是电力系统规划和运行中的重要任务,它的目标是求解电力系统中各节点的电压幅值和相角,以及线路的功率流向等参数,用于分析电力系统的稳定性和安全性,以及进行电力系统规划和调度。
下面是使用牛顿拉夫逊法计算潮流的一般步骤:步骤1:初始化首先,需要对电力系统的各个节点(包括发电机节点和负荷节点)的电压幅值和相角进行初始化,一般可以使用其中一种估计值或者历史数据作为初始值。
步骤2:建立潮流方程根据电力系统的潮流计算模型,可以建立节点电压幅值和相角的平衡方程,一般采用节点注入功率和节点电压的关系来表示。
潮流方程一般是一个非线性方程组,包含了各个节点之间的复杂关系。
步骤3:线性化方程组将潮流方程组进行线性化处理,一般采用泰勒展开的方法,将非线性方程组变为线性方程组。
线性化的过程需要计算雅可比矩阵,即方程组中的系数矩阵。
步骤4:求解线性方程组利用线性方程组的求解方法,比如高斯消元法或LU分解法等,求解线性方程组,得到电压幅值和相角的修正量。
步骤5:更新节点电压根据线性方程组的解,更新各个节点的电压幅值和相角,得到新的节点电压。
步骤6:检查收敛性判断节点电压的修正量是否小于设定的收敛阈值,如果满足收敛条件,则停止迭代;否则,返回步骤3,循环进行线性化方程组和线性方程组的求解。
步骤7:输出结果当潮流计算收敛时,输出最终的节点电压幅值和相角,以及线路的功率流向等参数。
牛顿拉夫逊法是一种高效、快速且收敛性良好的方法,广泛应用于电力系统潮流计算。
在实际应用中,可能会遇到多次迭代或者收敛性不好的情况,此时可以采用退火技术或其他优化算法进行改进。
此外,牛顿拉夫逊法的计算也可以并行化,利用多核处理器或者分布式计算集群来加速计算过程。
总之,牛顿拉夫逊法是一种重要的潮流计算方法,通过迭代计算逼近非线性方程组的解,可以得到电力系统中各节点的电压幅值和相角,用于分析电力系统的稳定性和安全性。
0 引言潮流是配电网络分析的基础,用于电网调度、运行分析、操作模拟和设计规划,同时也是电压优化和网络接线变化所要参考的内容.潮流计算通过数值仿真的方法把电力系统的详细运行情况呈现给工作人员,从而便于研究系统在给定条件下的稳态运行特点。
随着市场经济的发展,经济利益是企业十分看重的,而线损却是现阶段阻碍企业提高效益的一大因素.及时、准确的潮流计算结果,可以给出配电网的潮流分布、理论线损及其在网络中的分布,从而为配电网的安全经济运行提供参考.从数学的角度来看,牛顿—拉夫逊法能有效进行非线性代数方程组的计算且具有二次收敛的特点,具有收敛快、精度高的特点,在输电网中得到广泛应用.随着现代计算机技术的发展,利用编程和相关软件,可以更好、更快地实现配电网功能,本文就是结合牛顿—拉夫逊法的基本原理,利用C++程序进行潮流计算,计算结果表明该方法具有良好的收敛性、可靠性及正确性。
1 牛顿-拉夫逊法基本介绍1。
1 潮流方程对于N个节点的电力网络(地作为参考节点不包括在内),如果网络结构和元件参数已知,则网络方程可表示为:YV I (1—1)=式中,Y为N*N阶节点导纳矩阵;V为N*1维节点电压列向量;I为N*1维节点注入电流列向量。
如果不计网络元件的非线性,也不考虑移相变压器,则Y为对称矩阵。
电力系统计算中,给定的运行变量是节点注入功率,而不是节点注入电流,这两者之间有如下关系:ˆˆ=EI S(1—2)式中,S为节点的注入复功率,是N*1维列矢量;ˆS为S的共轭;ˆˆi diag ⎡⎤=⎢⎥⎣⎦E V 是由节点电压的共轭组成的N*N 阶对角线矩阵。
由(1-1)和(1-2),可得:ˆˆ=S EYV上式就是潮流方程的复数形式,是N 维的非线性复数代数方程组.将其展开,有:ˆi i iij j j iP jQ V Y V ∈-=∑ j=1,2,….,N (1—3)式中, j i ∈表示所有和i 相连的节点j ,包括j i =。
目录1任务书 (2)2.模型简介及等值电路 (3)3.设计原理 (4)4.修正方程的建立 (6)5.程序流程图及MATLAB程序编写 (8)6.结果分析 (14)7.设计总结 (17)8.参考文献 (17)《电力系统分析》课程设计任务书题目极坐标表示的牛顿拉夫逊法潮流计算程序设计学生姓名学号专业班级电气工程及其自动化设计内容与要求1. 设计要求掌握MATLAB语言编程方法;理解和掌握运用计算机进行潮流计算的基本算法原理;针对某一具体电网,进行潮流计算程序设计。
其目的在于加深学生对电力系统稳态分析中课程中基本概念和计算方法的理解,培养学生运用所学知识分析和解决问题的能力。
2. 内容1)学习并掌握MATLAB语言。
2)掌握变压器非标准变比概念及非标准变比变压器的等值电路。
掌握节点导纳矩阵的概念及导纳矩阵的形成和修改方法。
3)掌握电力系统功率方程、变量和节点分类。
4)掌握利用极坐标表示的牛-拉法进行潮流计算的方法和步骤。
5)选择一个某一具体电网,编制程序流程框图。
6)利用MATLAB语言编写该模型的潮流计算程序,并上机调试程序,对计算结果进行分析。
7)整理课程设计论文。
起止时间2013 年7 月 4 日至2013 年7月10日指导教师签名年月日系(教研室)主任签名年月日学生签名年月日2 模型简介及等值电路 2.1 课程设计模型:模型3电力网络接线如下图所示,各支路阻抗标幺值参数如下:Z12=0.02+j0.06,Z13=0.08+j0.24, Z23=0.06+j0.18, Z24=0.06+j0.12, Z25=0.04+j0.12, Z34=0.01+j0.03, Z45=0.08+j0.24, k=1.1。
该系统中,节点1为平衡节点,保持11.060V j =+&为定值;节点2、3、4都是PQ 节点,节点5为PV 节点,给定的注入功率分别为:20.200.20S j =+,3-0.45-0.15S j =,40.400.05S j =--,50.500.00S j =-+,5 1.10V =&。
电力系统的稳定性分析电力系统的稳定性分析是电力工程中的重要课题之一,它涉及到电力系统运行的可靠性和安全性。
稳定性分析主要考虑电力系统在各种外界扰动下的稳定性能,例如电力负荷突然增加或减小、供电故障等。
一、电力系统稳定性的概念电力系统稳定性指的是电力系统在负荷变化或外界扰动下,能够保持正常运行而不发生系统级别的不稳定或系统崩溃。
主要包括功率稳定性和动态稳定性两个方面。
1. 功率稳定性功率稳定性是指系统在负荷变化或供电故障的情况下,能够保持电压和频率稳定的能力。
这是电力系统必须具备的基本稳定性。
2. 动态稳定性动态稳定性是指电力系统在负荷突然变化或供电故障等大干扰下,能够在一定时间内恢复到稳定工作状态的能力。
这是保证系统能够迅速恢复到正常供电状态的重要指标。
二、电力系统稳定性分析方法在电力系统稳定性分析中,常用的方法主要有牛顿-拉夫逊法、潮流灵敏度法、方程迭代法和直接解法等。
1. 牛顿-拉夫逊法牛顿-拉夫逊法是常用的潮流计算方法,它通过迭代求解电流、电压和功率等参数来判断系统的稳定性。
通过计算节点电压和功率的变化情况,可以得出系统是否稳定以及稳定的程度。
2. 潮流灵敏度法潮流灵敏度法是一种通过计算电力系统中各个参数的灵敏度来评估系统稳定性的方法。
它可以分析发电机输出功率、传输线路电流和变压器负载等参数对系统稳定性的影响,有助于识别出系统中薄弱环节。
3. 方程迭代法方程迭代法是通过建立电力系统的状态方程,并利用迭代计算的方法来得出系统的稳定性。
通过不断迭代求解状态方程,得到系统的稳定情况。
4. 直接解法直接解法是指通过求解系统非线性方程组的方法来得到系统的稳定性。
这种方法常用于小规模系统或者用于求解系统的特定问题。
三、电力系统稳定性分析的应用电力系统稳定性分析在电力工程中有着广泛的应用,主要用于以下几个方面:1. 发电机组调度稳定性分析可以帮助电力系统运营人员制定合理的发电机组调度策略,以保证系统在负荷变化下的稳定运行。
牛顿—拉夫逊法在电力系统潮流计算的应用与分析潮流计算是电力系统进行稳定计算和故障分析的基础,可以得到各电网各节点的电压,并求得网络的潮流及网络中各元件的电力损耗,进而求得电能损耗。
随着现代电力系统的不断扩大,需要对传统的牛顿-拉夫逊法进行改进,降低初值选取的敏感性和提高收敛速度。
经典的牛顿法给定潮流计算时各节点的类型,确定导纳矩阵、修正方程和迭代收敛条件,将非线性方程组线性化为修正方程组反复迭代求解,因此收敛范围依赖电压的初值;同时求解雅克比矩阵计算量较大,影响计算速度。
1 算法原理1.1 原理介绍牛顿迭代法是取之后,找更接近的方程根,一步一步迭代,找到更接近方程根的近似根。
牛顿迭代法是求方程根的重要方法之一,最大优点是在方程的单根附近具有平方收敛,且还可用来求方程的重根、复根。
电力系统潮流计算,各个母线所供负荷的功率是已知的,各个平衡节点外的节点电压是未知的,可以根据网络结构形成节点导纳矩阵,由节点导纳矩阵列写功率方程,由于功率方程里功率是已知的,电压的幅值和相角是未知的,这样潮流计算的问题就转为求解非线性方程组的问题。
为便于用迭代法解方程组,需将上述功率方程改写成功率平衡方程,并对功率平衡方程求偏导,得出对应的雅可比矩阵,给未知节点赋电压初值,将初值带入功率平衡方程,得到功率不平衡量,这样由功率不平衡量、雅可比矩阵、未知节点电压不平衡量构成误差方程,解方程得到节点电压不平衡量,节点电压加上其不平衡量构成新的节点电压初值,将其带入原功率平衡方程,重新形成雅可比矩阵,计算新的电压不平衡量,这样不断迭代,一般迭代三到五次就能收敛。
1.2 牛顿—拉夫逊迭代法的一般步骤:(1)形成各节点导纳矩阵 Y。
(2)设节点电压的初始值 U、相角初始值 e、迭代次数初值 0。
(3)计算各节点的功率不平衡量。
(4)根据收敛条件判断是否满足,若不满足则向下进行。
(5)计算雅可比矩阵中的各元素。
(6)修正方程式节点电压。
牛顿拉夫逊法潮流计算牛顿-拉夫逊法(Newton-Raphson method)是一种用于求解非线性方程的数值方法。
它通过迭代逼近根的方式,将非线性方程转化为一系列的线性方程来求解。
在电力系统中,潮流计算用于确定电力网中节点的电压幅值和相角。
潮流计算是电力系统分析的重要基础,可以用于计算电力系统的潮流分布、功率损耗、节点电压稳定度等参数,为电力系统的规划、运行和控制提供参考依据。
牛顿-拉夫逊法是一种常用的潮流计算方法,它的基本思想是通过不断迭代来逼近电网的潮流分布,直到满足一定的收敛条件。
下面将对牛顿-拉夫逊法的具体步骤进行详细介绍。
首先,我们需要建立电力网络的节点潮流方程,即功率方程。
对于每一个节点i,其节点功率方程可以表示为:Pi - Vi * (sum(Gij * cos(θi - θj)) - sum(Bij * sin(θi -θj))) = 0Qi - Vi * (sum(Gij * sin(θi - θj)) + sum(Bij * cos(θi -θj))) = 0其中,Pi和Qi分别为节点i的有功功率和无功功率,Vi和θi分别为节点i的电压幅值和相角,Gij和Bij分别为节点i和节点j之间的导纳和电纳。
接下来,我们需要对每个节点的电压幅值和相角进行初始化。
一般情况下,可以将电压幅值设置为1,相角设置为0。
然后,我们可以开始进行迭代计算。
在每一轮迭代中,我们需要计算每个节点的雅可比矩阵和功率残差,然后更新电压幅值和相角。
雅可比矩阵可以通过对节点功率方程进行求导得到,具体如下:dPi/dVi = -sum(Vj * (Gij * sin(θi - θj) + Bij * cos(θi - θj)))dPi/dθi = sum(Vj * (Gij * Vi * cos(θi - θj) - Bij * Vi * sin(θi - θj)))dQi/dVi = sum(Vj * (Gij * cos(θi - θj) - Bij * sin(θi - θj)))dQi/dθi = sum(Vj * (Gij * Vi * sin(θi - θj) + Bij * Vi * cos(θi - θj)))功率残差可以通过将节点功率方程代入得到,如下:RPi = Pi - Vi * (sum(Gij * cos(θi - θj)) - sum(Bij *sin(θi - θj)))RQi = Qi - Vi * (sum(Gij * sin(θi - θj)) + sum(Bij *cos(θi - θj)))最后,我们可以使用牛顿-拉夫逊法的迭代公式来更新电压幅值和相角,具体如下:Vi(new) = Vi(old) + ΔViθi(new) = θi(old) + Δθi其中,ΔVi和Δθi分别为通过求解线性方程组得到的电压幅值和相角的增量。
安徽工程大学本科生课程设计说明书目录安徽工程大学课程设计任务书 (3)摘要 (5)Abstract (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)安徽工程大学课程设计任务书12系统接线图其中节点1为平衡节点,节点2、3、4、5为PQ节点。
摘要潮流计算,指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下,计算有功功率、无功功率及电压在电力网中的分布。
潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件,确定电力系统各部分稳态运行状态参数的计算。
通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。
待求的运行状态参量包括电网各母线节点的电压幅值和相角,以及各支路的功率分布、网络的功率损耗等。
它是基于配电网络特有的层次结构特性,论文提出了一种新颖的分层前推回代算法。
该算法将网络支路按层次进行分类,并分层并行计算各层次的支路功率损耗和电压损耗,因而可大幅度提高配电网潮流的计算速度。
论文在MATLAB环境下,利用其快速的复数矩阵运算功能,实现了文中所提的分层前推回代算法,并取得了非常明显的速度效益。
另外,论文还讨论发现,当变压器支路阻抗过小时,利用Π型模型会产生数值巨大的对地导纳,由此会导致潮流不收敛。
电力系统稳态潮流计算上机实验报告一、问题如下图所示的电力系统网络,分别用牛顿拉夫逊法、PQ解耦法、高斯赛德尔法、保留非线性法计算该电力系统的潮流。
发电机的参数如下,*表示任意值负荷参数如下,如上图所示的电力系统,可以看出,节点1、2、3是PQ节点,节点4是PV节点,而将节点5作为平衡节点。
根据问题所需,采用牛顿拉夫逊法、PQ解耦法、高斯赛德尔法、保留非线性法,通过对每次修正量的收敛判据的判断,得出整个电力系统的潮流,并分析这四种方法的收敛速度等等。
算法分析1.牛顿拉夫逊法节点5为平衡节点,不参加整个的迭代过程,节点1、2、3为PQ节点,节点4为PV 节点,计算修正方程中各量,进而得到修正量,判断修正量是否收敛,如果不收敛,迭代继续,如果收敛,算出PQ节点的电压幅值以及电压相角,得出PV节点的无功量以及电压相角,得出平衡节点的输出功率。
潮流方程的直角坐标形式,()()∑∑∈∈++-=ij j ij j ij i ij j ij j ij i i e B f G f f B e G e P()()∑∑∈∈+--=ij j ij j ij i ij j ij j ij i i e B f G e f B e G f Q直角坐标形式的修正方程式,11112n n n m n m -----∆⎡⎤⎡⎤∆⎡⎤⎢⎥⎢⎥∆=-⎢⎥⎢⎥⎢⎥∆⎣⎦⎢⎥⎢⎥∆⎣⎦⎣⎦PHN e Q M L f UR S修正方程式中的各量值的计算,()()][∑∑∈∈++--=∆ij j ij j ij i ij j ij j ij i is i e B f G f f B e G e p P()()][∑∑∈∈+---=∆ij j ij j ij i ij j ij j ij i is i e B f G e f B e G f Q Q)(2222i i is i f e U U +-=∆Jacobi 矩阵的元素计算,()()()ij i ij i i ijij j ij j ii i ii i jj iB e G f j i Q M G f B e B e G f j i e ∈-⎧≠∂∆⎪==⎨++-=∂⎪⎩∑()()()ij i ij i i ijij j ij j ii i ii i jj iG e B f j i Q L G e B f G e B f j i f ∈+⎧≠∂∆⎪==⎨--++=∂⎪⎩∑)()(202i j i j e e U R ijij i =≠⎩⎨⎧-=∂∆∂=)()(202i j i j f f U S ijij i =≠⎩⎨⎧-=∂∆∂=牛顿拉夫逊法潮流计算的流程图如下,2.PQ 解耦法如同牛顿拉夫逊法,快速解耦法的前提是,输电线路的阻抗要比电阻大得多,并且输电线路两端的电压相角相差不大,此时可利用PQ 快速解耦法,来计算整个电力系统网络的潮流。
第三节牛顿拉夫逊法潮流计算牛顿-拉夫逊法(Newton-Raphson method)是一种数值计算方法,用于求解非线性方程和潮流计算问题。
它是基于牛顿迭代法和拉夫逊迭代法的结合,可高效地求解电力系统潮流计算问题。
潮流计算是电力系统运行分析中的重要环节,其目标是确定系统中每个节点的电压和相角,并计算各个支路的电流,以评估系统的功率传输和稳定性。
在传统的高压电力系统中,由于负荷、发电机和传输线等元件的非线性特性,潮流计算问题呈现为非线性的数学方程组,通常采用迭代方法求解。
牛顿-拉夫逊法的基本思想是通过对方程组的线性化近似,迭代求解线性方程组的解,以接近方程组的精确解。
它通过将非线性方程组转化为以下形式进行迭代:F(x)=0其中,F(x)是非线性方程组的向量函数,x是未知向量。
牛顿-拉夫逊法的迭代过程可通过以下步骤进行:1.初始化变量:根据系统的初始状态进行节点电压和相角的初始化。
2.计算雅可比矩阵:通过对非线性方程组进行偏导,得到雅可比矩阵。
雅可比矩阵描述了各个节点潮流量与节点电压和相角之间的关系。
3.迭代计算:通过牛顿迭代法进行迭代计算,直到达到指定的收敛条件。
具体步骤为:a.解线性方程组:根据雅可比矩阵和当前节点电压和相角,求解线性方程组,得到修正量。
b.更新变量:根据修正量和当前节点电压和相角,更新节点电压和相角的值。
c.判断收敛:判断修正量是否满足收敛条件,如果满足则结束迭代计算,否则返回步骤a。
牛顿-拉夫逊法的优点是收敛速度快,精度高。
然而,它的缺点是对于方程组的收敛性和初始值的选择要求较高,存在收敛到局部最小值的问题。
为了克服这些问题,可以采用改进的牛顿-拉夫逊法,如增加松弛因子或采用多起点迭代法等。
总之,牛顿-拉夫逊法是一种高效的求解非线性方程组和潮流计算问题的数值方法。
它在电力系统潮流计算中广泛应用,帮助分析和评估电力系统的稳定性和功率传输能力。
随着电力系统的规模和复杂性的增加,牛顿-拉夫逊法的进一步改进和优化仍然是一个研究的热点问题。
电力系统牛顿拉夫逊潮流计算电力系统是现代社会不可或缺的基础设施之一、为了确保电力系统的安全和稳定运行,需要对电力系统的潮流进行计算和分析。
牛顿拉夫逊潮流计算是一种常用的潮流计算方法,下面将详细介绍该方法及其应用。
牛顿拉夫逊潮流计算是一种基于潮流方程的数学模型,用于计算电力系统中节点电压和线路功率等参数。
该方法是按照功率平衡和节点电压平衡原理建立的,可以通过迭代的方式求解电力系统的潮流分布。
首先,我们需要了解电力系统的基本元件和参数。
电力系统包括发电机、变压器、输电线路和负荷等。
发电机产生的电功率通过变压器输送到负荷上,同时输电线路的电阻和电抗对电功率的传输也起到了一定的影响。
负荷是电力系统的终端用户,需要提供稳定的电能供应。
在牛顿拉夫逊潮流计算中,我们首先需要确定电力系统中各个节点的功率注入值和节点电压大小。
节点功率注入包括发电机的有功功率和无功功率注入以及负荷的有功功率和无功功率消耗。
节点电压大小是指各个节点之间的电压差距。
其次,我们需要建立潮流方程。
潮流方程是通过节点电压和线路阻抗得到的,用于计算各个节点的电压大小和线路的功率输送。
潮流方程通常是一个非线性方程组,需要通过迭代的方式求解。
接下来,我们需要选择一个合适的起始值进行计算。
起始值可以通过经验值或者实测值确定,然后根据潮流方程进行迭代计算,直到满足一定的收敛条件为止。
迭代计算的具体过程是,首先将起始值代入到潮流方程中,计算得到新的节点电压和线路功率。
然后,根据新的节点电压和线路功率重新计算潮流方程,并对比上一次计算结果,判断是否满足收敛条件。
如果满足收敛条件,则计算结束,否则继续迭代计算,直到满足收敛条件为止。
牛顿拉夫逊潮流计算方法具有以下优点:计算结果准确,收敛速度快,适用于大范围的电力系统。
然而,该方法也存在一些缺点,比如计算复杂度较高,计算过程需要重复迭代。
牛顿拉夫逊潮流计算方法在电力系统中有着广泛的应用。
首先,该方法可以用于电力系统的规划和设计。
电力系统网络潮流计算—牛顿拉夫逊法牛顿拉弗逊法(Newton-Raphson Method)是一种常用的电力系统网络潮流计算方法,用于求解复杂电力系统中的节点电压和支路潮流分布。
本文将对牛顿拉弗逊法进行详细介绍,并讨论其优缺点及应用范围。
牛顿拉弗逊法的基本原理是通过迭代计算,将电力系统网络潮流计算问题转化为一个非线性方程组的求解问题。
假设电力系统有n个节点,则该方程组的节点电压和支路潮流分布可以通过以下公式表示:f(x)=0其中,f为非线性函数,x为待求解的节点电压和支路潮流分布。
通过泰勒展开,可以将f在其中一点x_k处展开为:f(x)≈f(x_k)+J_k(x-x_k)其中,J_k为f在x_k处的雅可比矩阵,x_k为当前迭代步骤的解。
通过令f(x)≈f(x_k)+J_k(x-x_k)=0,可以求解方程J_k(x-x_k)=-f(x_k),得到下一步的迭代解x_{k+1}。
通过不断迭代,可以逐步接近真实的解,直到满足收敛条件为止。
牛顿拉弗逊法的迭代公式如下:x_{k+1}=x_k-(J_k)^{-1}f(x_k)其中,(J_k)^{-1}为雅可比矩阵J_k的逆矩阵。
牛顿拉弗逊法的优点之一是收敛速度快。
相比其他方法,如高斯赛德尔法,牛顿拉弗逊法通常需要更少的迭代次数才能达到收敛条件。
这是因为牛顿拉弗逊法利用了函数的一阶导数信息,能够更快地找到接近解的方向。
然而,牛顿拉弗逊法也存在一些缺点。
首先,该方法要求求解雅可比矩阵的逆矩阵,计算量较大。
尤其是在大型电力系统网络中,雅可比矩阵往往非常大,计算逆矩阵的复杂度高。
其次,如果初始猜测值不合理,可能会导致算法无法收敛,需要选择合适的初始值,否则可能陷入局部极小值。
牛顿拉弗逊法在电力系统网络潮流计算中有广泛的应用。
该方法可以用于计算节点电压和支路潮流分布,提供电力系统分析和设计的重要数据。
它可以用于稳态分析、短路分析、负荷流分析等多种电力系统问题的求解。
这些问题在电力系统规划、运行和控制等方面都具有重要意义。
0引言潮流是配电网络分析的基础,用于电网调度、运行分析、操作模拟和设计规划,同时也是电压优化和网络接线变化所要参考的内容。
潮流计算通过数值仿真的方法把电力系统的详细运行情况呈现给工作人员,从而便于研究系统在给定条件下的稳态运行特点。
随着市场经济的发展,经济利益是企业十分看重的,而线损却是现阶段阻碍企业提高效益的一大因素。
及时、准确的潮流计算结果,可以给出配电网的潮流分布、理论线损及其在网络中的分布,从而为配电网的安全经济运行提供参考。
从数学的角度来看,牛顿-拉夫逊法能有效进行非线性代数方程组的计算且具有二次收敛的特点,具有收敛快、精度高的特点,在输电网中得到广泛应用。
随着现代计算机技术的发展,利用编程和相关软件,可以更好、更快地实现配电网功能,本文就是结合牛顿-拉夫逊法的基本原理,利用C++程序进行潮流计算,计算结果表明该方法具有良好的收敛性、可靠性及正确性。
1牛顿-拉夫逊法基本介绍1.1潮流方程对于N 个节点的电力网络(地作为参考节点不包括在内),如果网络结构和元件参数已知,则网络方程可表示为:=&&YV I (1-1) 式中,Y 为N*N 阶节点导纳矩阵;&V 为N*1维节点电压列向量;&I 为N*1维节点注入电流列向量。
如果不计网络元件的非线性,也不考虑移相变压器,则Y 为对称矩阵。
电力系统计算中,给定的运行变量是节点注入功率,而不是节点注入电流,这两者之间有如下关系:ˆˆ=&&&EI S (1-2)式中,&S 为节点的注入复功率,是N*1维列矢量;ˆ&S 为&S 的共轭;ˆˆi diag ⎡⎤=⎢⎥⎣⎦&&E V 是由节点电压的共轭组成的N*N 阶对角线矩阵。
由(1-1)和(1-2),可得:ˆˆ=&&&S EYV上式就是潮流方程的复数形式,是N 维的非线性复数代数方程组。
将其展开,有:ˆi i i ij jj iP jQ V Y V ∈-=∑&& j=1,2,….,N (1-3)式中,j i ∈表示所有和i 相连的节点j ,包括j i =。
将节点电压用极坐标表示,即令ii i V V θ=∠&,代入式(1-3)中则有: ()i i i i ij ij j j j iP jQ V G jB V θθ∈-=∠-+∠∑()()cos sin i j ij ij ij ij j iV V G jB j θθ∈=+-∑故有:()()cos sin sin cos i i j ij ij ij ij j ii i j ij ij ij ij j i P V V G B Q V V G B θθθθ∈∈⎧=+⎪⎨=-⎪⎩∑∑ i=1,2,…,N (1-4) 式(1-4)是用极坐标表示的潮流方程。
而节点功率误差:(cos sin )θθ∈∆=-+∑SPi ii j ij ij ij ij j iP P V V G B (1-5) (cos sin )θθ∈∆=--∑SP i i i j ij ij ij ij j iQ Q V V G B (1-6)式中:SP i P ,SP i Q 为节点i 给定的有功功率及无功功率。
1.2牛顿-拉夫逊法基本原理1.2.1牛拉法的一般描述牛拉法是把非线性方程式的求解过程变成反复对相应的线性方程式的求解过程,即非线性问题通过线性化逐步近似,这就是牛拉法的核心。
下面以非线性方程式的求解过程来进行说明。
设电力网络的节点功率方程一般形式如下:()SP =y y x (1-7)式中,SP y 为节点注入功率给定值;y 为SP y 对应的物理量和节点电压之间的函数表达式;x 为节点电压。
写成功率偏差的形式:()()SP =-=0f x y y x (1-8)应用牛拉法求解如下。
在给定的初值()0x 处将式(1-8)作一阶泰勒展开:()()()00T x ∂+∆=∂f f x x x定义T∂=∂f J x为潮流方程的雅克比矩阵,0J 为J 在()0x 处的值,则有: ()()10-∆=-x J f x用∆x 修正()0x 就得到x 的新值。
如果用k 表示迭代次数,写成一般的表达式,有:()()()()()()()()11k k k k k k -+⎧∆=-⎪⎨⎪=+∆⎩x J x f x x x x (1-9) 对于潮流收敛的情况,()1k +x 应比()kx 更接近于解点。
收敛条件为:()()max ki f x ε<由简单迭代法收敛性分析的结论知,越接近解点,牛顿-拉夫逊法收敛越快,它具有二阶收敛速度。
由图1.1可以直观地了解牛拉法的步骤:图1.1 牛顿-拉夫逊法的几何解释1.2.2极坐标的牛顿-拉夫逊法在极坐标中,()f x 有如下的形式:()()()()(),,,,nSP SP n r-⎡⎤⎡⎤∆-⎢⎥⎢⎥==⎢⎥⎢⎥∆-⎣⎦⎣⎦P V P P V f x Q V Q Q V θθθθ(1-10) 共2n-r 个方程,状态变量为:1212T T Tn n r V V V θθθ-⎡⎤⎡⎤==⎣⎦⎣⎦L L x V θ 共2n-r 个待求量。
r 个PV 节点的电压幅值给定,不需求解。
潮流雅克比矩阵的维数是(2n-r )*(2n-r),结构如下:nT T TTTn r-⎡⎤∂∆∂∆⎢⎥∂∂∂⎢⎥==∂∂∆∂∆⎢⎥⎢⎥∂∂⎣⎦PPfV J x Q Q V θθ 上式右侧的对电压幅值的偏导数项中的电压幅值的阶数减少了1,为使雅克比矩阵的各部分子矩阵具有一致的形式,在实际计算中,常将该项乘以电压幅值,并选取1122////Tn r n r V V V V V V --⎡⎤⎡⎤∆=∆∆∆⎣⎦⎣⎦L V V 作为待求的修正量,则雅克比矩阵可写成:nT T TT Tn r-⎡⎤∂∆∂∆⎢⎥∂∂∂⎢⎥==∂∂∆∂∆⎢⎥⎢⎥∂∂⎣⎦P P V f V J x Q Q V V θθ(1-11)将式(1-10)和(1-11)代入式(1-9)的修正方程即可求得x 的修正量x ∆,用它修正x 直到()()max ki f x ε<为止。
将式(1-11)用下式表示:⎡⎤=⎢⎥⎣⎦H N J ML 其中每个字块的计算公式如下:对角元素: ''2''2''2''2,,,,i i ii i ii i ii ii i i i i ii i i ii i ii ii i i i i ii i ii i ii ii i i i iii i i ii i ii ii i i P Q H V H V H B V P P N V V N V N G V V Q P M V M V M G V Q Q L V V L V L B V V θθ⎧∂∆===+⎪∂⎪⎪∂∆⎪===--∂⎪⎨∂∆⎪===-⎪∂⎪⎪∂∆===-⎪∂⎩(1-12)非对角元素: '''''''''',cos sin ,cos sin ,,i ij i ij j ij ij ij ij ijji ij j i ij j ij ij ij ij ijji ij i ij j ij ijji ij j i ij j ij ijjP H V H V H B G P N V V N V N G B V Q M V M V M N Q L V V L V L H V θθθθθθ⎧∂∆===-⎪∂⎪⎪∂∆⎪===--∂⎪⎨∂∆⎪===-⎪∂⎪⎪∂∆===⎪∂⎩(1-13)2牛顿法潮流计算步骤2.1程序流程图在了解了牛拉法的原理之后,明确程序编写思路,如图 2.1、2.2所示。
其中图2.1中的“计算电压幅值和角度”步骤较多,单独用图2.2表示出来。
图2.1 牛顿法计算潮流的程序框图 图2.2 电压幅值和角度求解步骤框图当不符合收敛的条件“amontk>1”时,即认为计算不收敛。
具体程序见附录。
2.2计算步骤下面讨论的是极坐标形式的牛顿法求解过程,大致分为以下几个步骤: ① 形成节点导纳矩阵;② 给各节点电压设初值((0)(0),i i V θ);③根据式(1-12)、(1-13)生成雅克比矩阵(H 、N 、M 、L );④ 将节点电压初值代入式(1-5)、式(1-6),求出修正方程式的常数项向量,i i P Q ∆∆;⑤ 求解修正方程,得到电压幅值和角度;⑥判断是否收敛,若收敛,计算平衡节点和线路功率; ⑦输出结果,并结束。
3算例3.1系统模型本文以图3.1所示电力网络为例,调用基于牛顿-拉夫逊法的C++程序。
1.1:10.3j0.12+0.5j0.1+0.4j0.08+0.4j①②③④图3.1 系统模型其中节点4设为平衡节点,电压标幺值为1.05,计算误差为0.000001。
3.2输入与输出将图3.1所示模型的相关数据放在data.dat文件中图3.2 输入节点和支路数据对各个数字含义的解释如下:网络模型有四个节点,四条支路,编号见图3.1。
第一个零下面三行数为支路参数,分别表示三条支路的起始和终止节点编号,后面的为电阻、电抗和电纳,电导均为0,例如:1 2 0.1 0.4 0.01528。
第二个零下面的为变压器支路,各数字意义同支路参数。
接下去三行均为节点参数,分别表示注入有功功率和无功功率。
调用text.cpp文件,得到运行结果,见图3.3和图3.4。
图3.3 运行结果1图3.4运行结果23.3结果分析将上述仿真结果整理为表格3.1、3.2,其中“+”表示节点i输出功率给节点j,“-”表示节点j输出功率给i(纵向为i,横向为j)。
表3.1节点有功功率输入与输出节点号 1 2 3 41 0 +0.245981 -0.5 -0.0465632 -0.24316 0 0 -0.3129493 0.5 0 0 04 0.0482143 0.319671 0 0表3.2节点无功功率输入与输出节点号 1 2 3 41 0 -0.014708 -0.029001 -0.1361872 0.0110505 0 0 -0.140363 0.097016 0 0 04 0.10464 0.160255 0 0根据表格计算:节点1有功功率:0+0.245981-0.5-0.046563=-0.300582无功功率:0-0.014708-0.029001-0.136187=-0.179896节点2有功功率:-0.24316+0+0-0.312949=-0.556109无功功率:0.0110505+0+0-0.14036=-0.1293095 节点3有功功率:0.5+0+0+0=0.5无功功率:0.097016+0+0+0=0.097016节点4有功功率:0.0482143+0.319671+0+0=0.3678853无功功率:0.10464+0.160255+0+0=0.264895根据已知条件,两个PQ 节点的注入有功、无功分别为:P1=0.3,Q1=0.18; P2=0.55,Q2=0.13 潮流计算误差:10.3005820.3100%0.194%0.3iP P ∆-=⨯=110.180.179896100%0.0578%0.18Q Q ∆-=⨯=220.5561090.55100% 1.11%0.55P P ∆-=⨯=220.130.1293095100%0.531%0.13Q Q ∆-=⨯=可见,误差均在允许范围内。