电力系统分析潮流计算的计算机算法
- 格式：doc
- 大小：298.00 KB
- 文档页数：20
下载文档原格式
合集下载
电力系统稳态分析-第四章 复杂电力系统潮流的计算机算法
第二节 功率方程及其迭代解法
在实际电力系统中，已知的运行条件往往不是节点的注 入电流而是负荷和发电机的功率，而且这些功率一般不随节 点电压的变化而变化。同时，在节点功率不变的情况下，节 点的注入电流随节点电压的变化而变化。而在节点电压未知 的情况下，节点注入电流是无法求得的。因此，不能用上节 介绍的网络方程进行潮流计算。必须在已知节点导纳矩阵的 情况下，用已知的节点功率来代替未知的节点注入电流，建 立起潮流计算用的节点功率方程，才能求出节点电压，进而 求出整个系统的潮流分布。
非线性方程组没有直接的解析方法，只能用迭代求解发方法。
第二节 功率方程及其迭代解法
为了更好的理解功率方程的意义，先以两母线系统为例， 然后推广到n母线系统 1、两母线系统的功率方程 G 1
~ ~ SG1 PG1 jQG1 SG 2 PG 2 jQG 2
U1
等值电源功率
SG1和SG2 ？
3) 输电线模型：是一个分布参数的电路，可用一个集中 参数的∏型等值电路表示；
4) 变压器模型：通常用集中参数的г型等值电路表示。
第一节 电力网络方程
要进行复杂系统的潮流计算，借助计算机程序进行计 算时，需要建立电力网络的网络方程。它是反映系统中 电流与电压之间相互关系的数学方程。需要对电力网进 行数学的抽象。
i
Yii Y jj yij Yij Y ji yij
yij
电力网
j
Yii Y jj Yii Yij Y ji Yij
(0) (0)
Yii Yij
第一节 电力网络方程
（4）在原有网络的节点i、j之间的导纳由yij改变为y'ij
i
Yii y yij ij Yij Y ji yij y ij
第四章电力系统潮流的计算机算法
1 z ij
(4) 原有节点ij之间阻抗由Zij变为Zij’
i j
－Zij
Yii
Yj
j
y i' jyi
j
1 z'ij
1 zij
Z’ij
Yij＝Yji
yi
j
y
i'
j＝z1ij
1 z'ij
(4) 原有节点ij之间变压器的变比由K*变为K*’时。
i j
返回
-ZT K*:1
ZT K’*:1
Z1 Y T(k-1 )/k
(2)节点导纳矩阵是稀疏矩阵，非对角非零 元素的个数等于对应节点所连的不接地 支路数。
(3)对角元素（自导纳）等于相应节点所连 支路的导纳之和。
(4)非对角元素（互导纳）等于两节点间支 路导纳的负值。
(5)节点导纳矩阵是对称方阵，只需求上三 角或是下三角元素。
标准变比：在采用有名值时，是指归算参数时所 取的变比。采用标么值时，是指折算参数时所 取各基准电压之比。
•
I1
Z 1 U 1 k :1
I1
•
I2
ZT
U2
Z2
U 1/k
I2
～～
S1 ＝ S 2
U1I 1 U1I2 k
I1 I2 / k U 1/kU 2I 2ZT
I1
U1 ZT k 2
U2 ZT k
I2
U1 ZT k
U2 ZT
I 1(y10y12)U 1y12 U 2 I 2 y2U 1 1(y20y21)U 2
2n个扰动变量是已知的，给定2（n－1）个控制变量， 给定2个状态变量，要求确定2（n－1）个状态变量。 已知：4n个变量，待求：2n个变量
电力系统中的潮流计算与分析
电力系统中的潮流计算与分析摘要本文介绍了电力系统中的潮流计算与分析，潮流计算是电力系统计算的基础，通过对电力系统中的电流、电压和功率进行计算和分析，可以有效地评估电力系统的稳定性和安全性。
在本文中，我们讨论了潮流计算的原理和方法，并介绍了一种基于改进的高斯-赛德尔迭代算法的潮流计算方法。
同时，我们还介绍了一种基于Python语言的潮流计算程序的设计和实现，该程序可以对电力系统进行潮流计算和分析，并生成相关的报告和图表。
最后，我们利用该程序对IEEE 14节点测试系统进行了潮流计算和分析，并分析了系统的稳定性和安全性。
关键词：电力系统；潮流计算；高斯-赛德尔迭代算法；Python语言AbstractThis paper introduces the load flow calculation and analysis in power system. Load flow calculation is the basis of power system calculation. By calculating and analyzing the current, voltage and power in the power system, the stability and safety of the power system can be effectively evaluated. In this paper, we discuss the principles and methods of load flow calculation, and introduce an improved Gauss-Seidel iterative algorithm based load flow calculation method. At the same time, we also introduce the design and implementation of a load flow calculation program based on the Python language. The program can perform load flow calculation and analysis on the power system, and generate relevant reports and charts. Finally, we use the program to perform load flow calculation and analysis on the IEEE 14-bus test system, and analyze the stability and safety of the system.Keywords: power system; load flow calculation; Gauss-Seidel iterative algorithm; Python language一、引言电力系统是现代工业和生活的基础设施之一，它承担着输送和分配电能的重要任务。
电力系统中的潮流计算与优化方法
电力系统中的潮流计算与优化方法潮流计算是电力系统运行和规划中的重要环节，它用于计算电力系统中各节点的电压、相角、有功、无功功率以及线路、变压器等的潮流分布情况。
对电力系统进行潮流计算可以帮助电力系统运行人员了解系统的稳定性、可靠性以及容载能力，也可以为电力系统规划提供数据支持。
本文将介绍电力系统潮流计算的基本方法与优化技术。
一、潮流计算的基本方法1.1 普通潮流计算方法潮流计算的基本方法是牛顿-拉夫逊迭代法（Newton-Raphson Iteration Method）和高尔顿法（Gauss-Seidel Method）。
牛顿-拉夫逊迭代法主要是通过不断迭代求解雅可比矩阵的逆，直到迭代误差小于给定阀值时停止迭代；高尔顿法则是逐一更新所有节点的电压与相角，直至所有节点的迭代误差都小于给定阀值。
1.2 快速潮流计算方法在大型电力系统中，普通的潮流计算方法计算速度较慢。
因此，研究人员提出了一些针对快速潮流计算的方法，如快速牛顿-拉夫逊法（Fast Newton-Raphson Method）和DC潮流计算方法。
快速牛顿-拉夫逊法通过简化牛顿-拉夫逊法的迭代公式，减少计算量，提高计算速度；DC潮流计算方法则是将潮流计算问题转化为一个线性方程组的求解问题，进一步提升计算效率。
二、潮流计算的优化技术2.1 改进的潮流计算算法为了提高潮流计算的准确性和收敛速度，研究人员提出了一些改进的潮流计算算法。
其中，改进的牛顿-拉夫逊法（Improved Newton-Raphson Method）是一种结合牛顿-拉夫逊法和割线法的算法，通过混合使用这两种方法，实现在减小迭代误差的同时加快计算速度。
此外，基于粒子群优化算法（Particle Swarm Optimization）和遗传算法（Genetic Algorithm）的潮流计算算法也得到了广泛研究和应用。
2.2 潮流优化潮流计算不仅可以用于分析电力系统的工作状态，还可以作为优化问题的约束条件。
潮流计算的计算机方法
一、潮流计算的计算机方法对于复杂网络的潮流计算，一般必须借助电子计算机进行。
其计算步骤是：建立电力网络的数学模型，确定计算方法、制定框图和编制程序。
本章重点介绍前两部分，并着重阐述在电力系统潮流实际计算中常用的、基本的方法。
1，电力网络的数学模型电力网络的数学模型指的是将网络有关参数相变量及其相互关系归纳起来所组成的．可以反映网络性能的数学方程式组。
也就是对电力系统的运行状态、变量和网络参数之间相互关系的—种数学描述。
电力网络的数学模型有节点电压方程和回路电流方程等，前者在电力系统潮流计算中广泛采用。
节点电压方程又分为以节点导纳矩阵表示的节点电压方程和以节点阻抗矩阵表示的节点电压方程。
（1）节点导纳矩阵在电路理论课中。
已讲过了用节点导纳矩阵表示的节点电压方程：对于n个节点的网络其展开为：上式中，I是节点注入电流的列向量。
在电力系统计算中，节点注入电流可理解为节点电源电流与负荷电流之和，并规定电源向网络节点的注人电流为正。
那么，只有负荷节点的注入电流为负，而仅起联络作用的联络节点的注入电流为零。
U是节点电压的列向量。
网络中有接地支路时，通常以大地作参考点，节点电压就是各节点的对地电压。
并规定地节点的编号为0。
y是一个n×n阶节点导纳矩阵，其阶数n就等于网络中除参考节点外的节点数。
物理意义：节点i单位电压，其余节点接地，此时各节点向网络注入的电流就是节点i 的自导纳和其余节点的与节点i之间的互导纳。
特点：对称矩阵，稀疏矩阵，对角占优(2) 节点阻抗矩阵对导纳阵求逆，得：其中称为节点阻抗矩阵，是节点导纳矩阵的逆阵。
物理意义：节点i注入单位电流，其余节点不注入电流，此时各节点的电压就是节点i 的自阻抗和其余节点的与节点i之间的互阻抗。
特点：满阵，对称，对角占优2，功率方程、变量和节点分类（1）功率方程已知的是节点的注入功率，因此，需要重新列写方程： **==B B B B B U S I U Y其展开式为： i i i nj j ij U jQ P U Y ~1-=∑= 所以：∑=**=+nj jij i i i U Y U jQ P 1 展开写成极坐标方程的形式：)cos sin ()sin cos (11ij ij ij ij n j j i i ij ij ij ij n j j i i B G U U Q B G U U P δδδδ-=+=∑∑==所以节点的功率方程为：)cos sin ()sin cos (11ij ij ij ij n j j i di Gi i ij ij ij ij nj j i di Gi i B G U U Q Q Q B G U U P P P δδδδ---=∆+--=∆∑∑==（2） 变量分类负荷消耗的有功、无功功率取决于用户，因而是无法控制的，故称为不可控变量或扰动变量。
电力系统分析第三章-新
是已知的，每个节点
•3.2 功率方程
•变量的分类： ① 不可控变量（扰动变量）：PLi，QLi――由用户决定，无
法由电力系统控制； • ② 控制变量：PGi，QGi――由电力系统控制； ③ 状态变量：Ui，δi――受控制变量控制；其中Ui 主要受 ④ QGi 控制，δi 主要受PGi 控制。 • ☆ 若电力系统有n个节点，则对应共有6n个变量，其中不可 • 控变量、控制变量、状态变量各2n个； • ☆ 每个节点必须已知或给定其中的4个变量，才能求解功率 • 方程。
•
待求的是等值电源无功功率 QGi和节点电压相位角 δi 。
•3.2 功率方程
•选择：通常可以将有一定无功储备的发电厂母线和有一定无
•
功电源的变电所母线看作PV节点。
•3、平衡节点：
• 特点：进行潮流计算时通常只设一个平衡节点。给定平衡节
•
点的是等值负荷功率PLs 、QLs和节点电压的幅值Us 和
。
•⑦ 计算平衡节点功率和线路功率。
•3.3 潮流分布计算的计算机算法
•潮流计算流程 图（极坐标）
•3.3 潮流分布计算的计算机算法
•三、PQ分解法潮流计算：
•
也称牛顿－拉夫逊法快速解耦法潮流计算
•1、问题的提出:牛顿-拉夫逊法分析
•(1) 雅可比矩阵 J 不对称；
•(2) J 是变化的，每一步都要重新计算，重新分析；
；
• ⑤ 利用x (1) 重新计算∆f (1)和雅可比矩阵J (1)，进而得到∆x (1)
；
• 如此反复迭代：
；直至解出精确解或
• 得到满足精度要求的解。
•3.3 潮流分布计算的计算机算法
•二、牛顿－拉夫逊法潮流计算：迭代求解非线性功率方程
电力系统潮流计算方法分析
电力系统潮流计算方法分析电力系统潮流计算是电力系统运行中的基础性分析方法之一，它用于求解电力系统中各个节点的电压、相角以及线路的功率、电流等变量。
潮流计算是电力系统规划、运行和控制等方面的重要工具。
本文将对电力系统潮流计算方法进行分析。
电力系统潮流计算方法主要有两种，即直接法和迭代法。
直接法又分为解析法和数值法，迭代法包括高斯赛德尔迭代法、牛顿-拉夫逊迭代法等。
解析法是通过电力系统各个节点之间的网络拓扑关系和节点电压平衡条件的方程式，直接求解节点电压和线路功率等参数。
解析法的优点是计算速度快，但其适用范围较窄，主要适用于小型简单电力系统，对于大型复杂电力系统的潮流计算会出现计算量庞大的问题。
数值法是通过将连续变量离散化，将微分方程转化为差分方程，并利用数值解法求解离散的方程组来得到电力系统潮流计算结果。
数值法的优点是适用范围广，能够处理大型复杂电力系统的潮流计算，但其缺点是计算速度相对较慢。
在迭代法中，高斯赛德尔迭代法是一种经典的迭代法，它通过先假设节点电压的初值，然后利用节点注入功率与节点电压之间的关系不断迭代计算，最终达到收敛条件为止。
高斯赛德尔迭代法的优点是收敛速度快，计算精度高，但其缺点是收敛性有时不易保证，并且计算速度会随着系统规模的增大而变慢。
牛顿-拉夫逊迭代法是一种基于牛顿迭代法的改进方法，它引入雅可比矩阵，通过牛顿迭代法的迭代过程来求解节点电压和线路功率等参数。
牛顿-拉夫逊迭代法的优点是收敛性好，计算速度快，但其缺点是在实际应用中需要预先计算雅可比矩阵，会增加计算的复杂度。
综上所述，电力系统潮流计算方法有直接法和迭代法两种，其中直接法包括解析法和数值法，迭代法包括高斯赛德尔迭代法和牛顿-拉夫逊迭代法。
在实际应用中，根据电力系统的规模和复杂程度选择合适的方法进行潮流计算，以得到准确可靠的计算结果。
此外，随着计算机技术的不断发展，还可以利用并行计算和分布式计算等方法来提高潮流计算的效率。
潮流计算的主要方法
潮流计算的主要方法
最近几年，随着计算机仿真技术和复杂系统全面发展，潮流计算也受到越来越多的重视。
潮流计算是研究不同电力网络的物理特性和操作规律的一项重要工作。
针对潮流计算的主要方法，总结如下：
一、基于动力学的方法
1. 碰撞模型：根据动力学方法，计算电力系统的运行稳定性。
基于动力学的碰撞模型能够快速而精确地预测两个潮流的变化情况。
2. 时变快速收敛：在碰撞模型的基础上，为快速求解电力系统潮流，提出了时变快速收敛算法。
可以更快地获得潮流解。
二、基于牛顿迭代法的方法
1.牛顿迭代潮流计算方法：根据牛顿迭代法，采用迭代算法，求解电力系统潮流运行状态。
2. 功率流计算方法：计算机基于牛顿迭代法，快速求解节点电能的功率流公式。
可以有效的缩短潮流计算的时间，提高计算效率。
三、基于模糊聚类算法的方法
1. 基于模糊聚类的潮流计算方法：采用模糊聚类算法，对潮流计算进行多维度分析，可以得出最优的潮流结果。
2. 基于模糊划分的多目标模糊控制：根据模糊聚类理论，对潮流算法进行最佳控制，以满足电力网不同优化目标。
四、基于期望最大化的方法
1、基于粒子群优化的潮流计算方法：采用粒子群优化算法，将电力网潮流计算定义为多目标最优化问题，以期望最大化来求解潮流值，提高计算效率。
2、基于遗传算法的潮流计算方法：遗传算法利用进化过程来搜索全局最优解，使用遗传变异原则来改变候选解，以期望最大化来求解潮流计算问题。
简单电力系统分析潮流计算
简单电力系统分析潮流计算电力系统潮流计算是电力系统分析中的一项重要任务。
其目的是通过计算各个节点的电压、电流、有功功率、无功功率等参数，来确定系统中各个元件的运行状态和互相之间的相互影响。
本文将介绍电力系统潮流计算的基本原理、计算方法以及应用。
潮流计算的基本原理是基于电力系统的节点电压和支路功率之间的网络方程。
通过对节点电压进行迭代计算，直到满足所有支路功率平衡方程为止，得到系统的运行状态。
潮流计算的基本问题可以表示为以下方程组：P_i = V_i * (G_i * cos(θ_i - θ_j ) + B_i * sin(θ_i -θ_j )) - V_j * (G_i * cos(θ_i - θ_j ) - B_i * sin(θ_i -θ_j )) (1)Q_i = V_i * (G_i * sin(θ_i - θ_j ) - B_i * cos(θ_i -θ_j )) - V_j * (G_i * sin(θ_i - θ_j ) + B_i * cos(θ_i -θ_j )) (2)其中，P_i为节点i的有功功率注入；Q_i为节点i的无功功率注入；V_i和θ_i分别为节点i的电压幅值和相角；V_j和θ_j分别为节点j的电压幅值和相角；G_i和B_i分别为支路i的导纳的实部和虚部。
对于一个电力系统，如果知道了节点注入功率和线路的导纳，就可以通过潮流计算求解出各节点的电压和功率。
这是一种不断迭代的过程，直到系统达到平衡状态。
潮流计算的方法有多种，常见的有高斯-赛德尔迭代法、牛顿-拉夫逊迭代法等。
其中，高斯-赛德尔迭代法是最常用的一种方法。
高斯-赛德尔迭代法的思想是从已知节点开始，逐步更新其他节点的电压值，直到所有节点的电压值收敛为止。
具体步骤如下：1.初始化所有节点电压的初始值；2.根据已知节点的注入功率和节点电压，计算其他节点的电压值；3.判断节点电压是否收敛，如果收敛则结束计算，否则继续迭代；4.更新未收敛节点的电压值，返回步骤2高斯-赛德尔迭代法的优点是简单有效，但其收敛速度较慢。
电力系统分析课件：第4章 电力系统潮流的计算机算法
Yjj Yjj Yjj yij
Yji Yji Yji yij
Y11 Y12 Y1i Y1 j Y1n
Y21 Y22 Y2i Y2 j Y2n
Y i1 Y i2 Y ii Y ij Y in
Y j1 Y j2 Y ji Y jj Y jn
nn
Yn1 Yn2 Yni Ynj Ynn
•
I2
I1 I2
Y11 Y21
0 Y31
Y12 Y22 Y32
Y13 Y23 Y33
UU12 U 3
一般化表示：
节点电压 列向量
n*n阶，阶 数不包括参
考节点
节点注入电 流列向量
I1 I2
Y11 Y21
Ii In
Yi1 Yn1
矩阵向量表示：
Y12 Y22
1) k
4 功率方程
前文已述，若已知： ① 网络的导纳矩阵YB（ YB体现了网络的各个参
数） ② 各母线注入电流IB 则：通过节点电压方程求逆：UB YB1I B 即可：求出各节点电压：U i
进而：求得各支路电流： Iij yij (Ui U j )
然而：
在电力系统中，通常已知的不是各节点注入电
Yii Yii Yii yij yij Yjj Yjj Yjj yij yij
Yij Yij Yij yij yij
Yij Yij Yij yij yij
（5）若在原有的i,j节点的变压器的变比由 k 变为 k
结果：相当于在i、j节点之间切除一台变压器变比为 k 的变
ZB
Z i1
Z n1
Z12 Z1i Z1n Y11
Z22
Z2i
Z2n
Y21
Zi 2
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
潮流计算的计算机算法实验报告姓名：学号：班级：一、实验目的掌握潮流计算的计算机算法。
熟悉MATLAB，并掌握MATLAB程序的基本调试方法。
二、实验准备根据课程内容，熟悉MATLAB软件的使用方法，自行学习MATLAB程序的基础语法，并根据所学知识编写潮流计算牛顿拉夫逊法(或PQ分解法) 的计算程序，用相应的算例在MATLAB上进行计算、调试和验证。
三、实验要求每人一组，在实验课时内，用MATLAB调试和修改运行程序，用算例计算输出潮流结果。
四、实验程序clear; %清空内存n=input('请输入节点数:n=');n1=input('请输入支路数:n1=');isb=input('请输入平衡节点号:isb=');pr=input('请输入误差精度:pr=');B1=input('请输入支路参数:B1=');B2=input('请输入节点参数:B2=');X=input('节点号和对地参数:X=');Y=zeros(n);Times=1;%一：创建节点导纳矩阵for i=1:n1if B1(i,6)==0 %不含变压器的支路p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/B1(i,3);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);else %含有变压器的支路p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5));Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/B1(i,3);Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3));endendY;%将OrgS、DetaS初始化OrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1);%二：创建OrgS，用于存储初始功率参数h=0;j=0;for i=1:n %对PQ节点的处理if i~=isb&B2(i,6)==2h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j ))*Imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*i mag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendend%三：对PV节点的处理，注意这时不可再将h初始化为0for i=1:nif i~=isb&B2(i,6)==3h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j ))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*i mag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendendOrgS;%四：创建PVU 用于存储PV节点的初始电压PVU=zeros(n-h-1,1);t=0;for i=1:nif B2(i,6)==3t=t+1;PVU(t,1)=B2(i,3);endendPVU;%五：创建DetaS,用于存储有功功率、无功功率和电压幅值的不平衡量h=0;for i=1:n %对PQ节点的处理if i~=isb&B2(i,6)==2h=h+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);endendt=0;for i=1:n%六：对PV节点的处理，注意这时不可再将h初始化为0if i~=isb&B2(i,6)==3h=h+1;t=t+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag( B2(i,3))^2;endendDetaS;%七：创建I，用于存储节点电流参数i=zeros(n-1,1);h=0;for i=1:nif i~=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));endendI;%八：创建Jacbi(雅可比矩阵)Jacbi=zeros(2*n-2);h=0;k=0;for i=1:n %对PQ节点的处理if B2(i,6)==2h=h+1;for j=1:nif j~=isbk=k+1;if i==j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I (h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;endendendendendk=0;for i=1:n %对PV节点的处理if B2(i,6)==3h=h+1;for j=1:nif j~=isbk=k+1;if i==j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I (h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1)) ;Jacbi(2*h,2*k-1)=2*imag(B2(i,3));Jacbi(2*h,2*k)=2*real(B2(i,3));else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;endendendendendJacbi;%九：求解修正方程，获取节点电压的不平衡量DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU;%修正节点电压j=0;for i=1:n %对PQ节点处理if B2(i,6)==2j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfor i=1:n %对PV节点的处理if B2(i,6)==3j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendB2;%十：开始循环**********************************************************************while abs(max(DetaU))>prOrgS=zeros(2*n-2,1); %!!!初始功率参数在迭代过程中是不累加的，所以在这里必须将其初始化为零矩阵h=0;j=0;for i=1:nif i~=isb&B2(i,6)==2h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j ))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*i mag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendendfor i=1:nif i~=isb&B2(i,6)==3h=h+1;for j=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j ))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j, 3)));OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*i mag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendendOrgS;%创建DetaSh=0;for i=1:nif i~=isb&B2(i,6)==2h=h+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);endendt=0;for i=1:nif i~=isb&B2(i,6)==3h=h+1;t=t+1;DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i, 3))^2;endendDetaS;%创建Ii=zeros(n-1,1);h=0;for i=1:nif i~=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));endendI;%创建JacbiJacbi=zeros(2*n-2);h=0;k=0;for i=1:nif B2(i,6)==2h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I (h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1)) ;Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endif k==(n-1)k=0;endendendendendk=0;for i=1:nif B2(i,6)==3h=h+1;for j=1:nif j~=isbk=k+1;if i==jJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I (h,1));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1)) ;Jacbi(2*h,2*k-1)=2*imag(B2(i,3));Jacbi(2*h,2*k)=2*real(B2(i,3));elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endif k==(n-1)k=0;endendendendendJacbi;DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU;%修正节点电压j=0;for i=1:nif B2(i,6)==2j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfor i=1:nif B2(i,6)==3j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendB2;Times=Times+1; %迭代次数加1endTimes;五、实验流程六、实验结果参数输入：运行结果：七、实验体会通过这次实验，让我第一次接触到了MATLAB,并深切体会到了它的强大之处；潮流计算的计算机算法的实现不仅巩固了我的学过的知识，还让我学到一些MATLAB的编程，虽然在实验的过程中出现了很多的错误，但在老师的细心指导下，问题都解决啦；计算机为我们省去了大量的人工计算，希望在以后的学习中能接触到更多的软件，学习到更多的知识。