基于matlab的电力系统潮流计算课程设计(毕业论文).docx

课程设计论文

基于MATLAB的电力系统潮流计算

目录：

一、电力系统潮流计算概述 3
二、潮流计算方法概述 3
2．1节点分类 3
2.2牛顿—拉夫逊法概要 3
2.2.1牛顿—拉夫逊法迭代原理 3
2．2.1 牛顿法的框图及求解过程 3
三、课程设计任务 3
3.1题目原始数据 3
3.2课程设计要求 3
四、初步分析 3
4.1节点设置及分类 3
4.2 参数计算 3
4.3等值电路绘制 3
五、潮流计算 3
5.1 给定负荷下的潮流计算 3
5.1.1 B1/B2矩阵的形成 3
5.1.2 潮流调整及分析 3
5.1.3 仿真比较 3
5.2 变电所负荷变化时的潮流计算 3
5.2.1 4个变电所负荷同时以2%的比例增大 3
5.2.2 4个变电所负荷成比例2%下降 3
5.2.3 1、4号负荷下降，2、3号负荷上升 3
5.2.4 仿真比较 3
5.3断线潮流计算 3
5.5.1 断开1、3节点间的一条支路 3
5.3.2 断开1、5支路的一条线 3
5.3.3 断开1、7支路的一条线 3
5.3.4 断开2、9支路的一条线 3
5.3.5 断开7、9支路的一条线 3
5.3.6 仿真比较 3
六、潮流分析总结 3
七、心得体会 3

参考文献： 3
附录： 3

潮流计算课程设计
摘要
电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压。各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量的分析比较供电方案或运行方式的合理性、可靠性和经济性。此外，电力系统的潮流计算也是计算机系统动态稳定和静态稳定的基础，所以潮流计算是研究电力系统的一种重要的计算。
潮流计算在数学上是多元非线性方程组的求解问题，求解的方法有很多种，牛顿—拉夫逊法是数学上解非线性方程组的有效方法，有较好的收敛性。本文应用了电力系统潮流计算仿真软件DDTR与利用程序计算的结果进行比较，使计算的结果更加准确。利用成形的程序对系统中出现的各种情况，例如负荷的变化以及线路上所发生的变化进行计算，并对母线上不满足范围的数据进行调控，使得系统处于一个较稳定的状态。

关键词：牛顿—拉夫逊法 MATLAB DDRTS 潮流计算

一、电力系统潮流计算概述
在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。为了保证电力系统的稳定运行，要进行潮流调节。随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。它的发展主要围绕这样几个方面：计算方法的收敛性、可靠性；计算速度的快速性；对计算机存储容量的要求以及计算的方便、灵活等。 
牛顿-拉夫逊法是数学中解决非线性方程式的典型方法，有较好的收敛性。。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高潮流计算的效率。
二、潮流计算方法概述
2．1节点分类
常规的电力系统潮流计算中一般具有三种类型的节点：PQ、PV及平衡节点。一个节点有四个变量，即注入有功功率、注入无功功率，电压大小及相角。常规的潮流计算一般给定其中的二个变量：PQ节点（注入有功功率及无功功率），PV节点（注入有功功率及电压的大小），平衡节点（电压的大小及相角）。
第一类称PQ节点：等值负荷功率PLi、QLi和等值电源功率PGi、QGi是给定的，从而注入功率P、Q是给定的，待求的则是节点电压的大小Ui和相位角。属于这类节点的有按给定有功、无功率发电的发电厂母线和没有其他电源的变电所母线。 
第二类称PV节点：等值负荷和等值电源的有功功率PLi、PGi是给定的，从而注入有功功率Pi是给定的。等值负荷的无功功率QLi和节点电压的大小Ui也是给定的。待求的则是等值电源的无功功率QGi，从而注入无功功率Qi和节点电压的相位角。有一定无功功率储备的发电厂和有一定无功功率电源的变电所母线都可以作为PV节点。
第三类平衡节点：潮流计算时一般只设一个平衡节点。等值负荷功率PLs、QLs是给定的，节点电压的大小和相位也是给定的。担负调整系统频率任务的发电厂母线往往被选作为平衡节点。
2.2牛顿—拉夫逊法概要
2.2.1牛顿—拉夫逊法迭代原理
已知一个变量X函数为： ，由适当的近似值出发，根据：

反复进行计算，当满足适当的收敛条件就是上面方程的根。这样的方法就是所谓的牛顿—拉夫逊法。
这一方法还可以做下面的解释，设第次迭代得到的解语真值之差，即的误差为时，则：

把在附近对用泰勒级数展开

上式省略去以后部分

的误差可以近似由上式计算出来。

比较两式，可以看出牛顿—拉夫逊法的休整量和的误差的一次项相等。
用同样的方法考虑，给出个变量的个方程：

对其近似解得修正量可以通过解下边的方程来确定：

式中等号右边的矩阵都是对于的值。这一矩阵称为雅可比（JACOBI）矩阵。按上述得到的修正向量后，得到如下关系

这比更接近真实值。这一步在收敛到希望的值以前重复进行，一般要反复计算满足

为预先规定的小正数，是第n次迭代的近似值。
2．2.1 牛顿法的框图及求解过程
（1）用牛顿法计算潮流时，有以下的步骤：
①给这各节点电压初始值；
②将以上电压初始值代入公式，求修正方程的常数项向量 ；
③将电压初始值在带入上述公式，求出修正方程中系数矩阵的各元素。
④解修正方程式；
⑤修正各节点电压，；
⑥将，在带入方程式，求出；
⑦检验是否收敛，即
如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功率，并打印
输出结果。如果不收敛，转回②进行下次迭代计算，直到收敛为止。
（2） 牛顿—拉夫逊法计算程序框图

启动

输入原始数据

形成节点导纳矩阵

分解各节点初始电压的实部和虚部
迭代次数K=0
求PQ节点的,，求PV节点的，
置节点号i=0
雅克比矩阵是否形成，i>n
求得雅克比矩阵各元素
增大节点号i=i+1
把雅克比矩阵单位化
求解修正方程，得，
求解最大修正量，
是否收敛
回带各电压新值，K=K+1
计算输出电压大小及相角，节电功率及支路损耗
停止
否
是

否
是

三、课程设计任务
3.1题目原始数据
1.系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。

变电所1
变电所2
变电所3
变电所4
35kV母线
10kV母线
35kV母线
10kV母线
一次侧电压220kV
一次侧电压220kV
母线1
母线3
母线2
线路长为100km
线路长为50km
线路长为90km
线路长为100km
2*QFQ-50-2
2*QFS-50-2
TQN-100-2
2*TQN-100-2
。。。。。。。。。。。。。
线路长为80km
。。。。。。。。。。。。。
电厂一
电厂二

2、发电厂资料：
母线1和2为发电厂高压母线，发电厂一总装机容量为（ 300MW ），母线3为机压母线，机压母线上装机容量（100MW），最大负荷和最小负荷为50MW和20MW；发电厂二总装机容量（200MW）。
3、变电所资料：
（一）变电所1、2、3、4低压母线的电压等级分别为：35KV 10KV 35KV 10KV（二）变电所的负荷分别为：60MW 50MW 40MW 60MW
（三） 每个变电所的功率因数均为cosφ=0.85；
（四） 变电所1和变电所3分别配有两台容量为75MVA的变压器，短路损耗414KW，短路电压（%）=16.7；变电所2和变电所4分别配有两台容量为63MVA的变压器，短路损耗为245KW，短路电压（%）=10.5；
4、输电线路资料：
发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为，单位长度的电抗为，单位长度的电纳为。
3.2课程设计要求
1. 对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。
2. 输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷情况下的潮流计算，并对计算结果进行分析。
3. 跟随变电所负荷按一定比例发生变化，进行潮流计算分析。
1) 4个变电所的负荷同时以2%的比例增大；
2) 4个变电所的负荷同时以2%的比例下降
3) 1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的负荷同时以2%的比例上升；
4. 在不同的负荷情况下，分析潮流计算的结果，如果各母线电压不满足要求，进行电压的调整。（变电所低压母线电压10KV要求调整范围在9.5-10.5之间；电压35KV要求调整范围在35-36之间）
5. 轮流断开支路双回线中的一条，分析潮流的分布。（几条支路断几次）
6. 利用DDRTS软件，进行绘制系统图进行上述各种情况潮流的分析，并进行结果的比较。
7. 最终形成课程设计成品说明书。
四、初步分析
4.1节点设置及分类
根据系统图可知此系统为两端供电网路，为了潮流计算的完整性，将母线1，2设为节点1，2，将变电所1、2、3、4的高低压侧分别设为节点3、4、5、6、7、8、9、10。并且，将节点1设为平衡节点，将节点2设为PV节点，其余节点设为PQ节点。则网络共具有10个节点和9条支路。
4.2 参数计算
根据题目原始资料，计算发电厂、变压器及线路的有名值参数。
（1）运用下列公式计算变压器低压侧参数：

（2）计算线路参数

将参数整理成如下表格：
首端号
末端号
阻抗有名值
电纳有名值
1
3
8.5+j20.1
j2.780e-4
1
5
4.25+j10.05
j1.390e-4
1
7
7.65+j18.09
j2.502e-4
2
9
8.5+j20.1
j2.780e-4
3
4
0.05455+j1.65024
0
5
6
3.375e-3+j0.10083
0
7
8
0.05455+j1.65024
0
7
9
6.8+j16.08
j2.224e-4
9
10
3.375e-3+j0.10083
0

（3）变电所负荷分别为：
变电所1 =60+j37.185
变电所2 =50+j30.987
变电所3 =40+j24.790
变电所4 =60+j37.185
(4）计算变压器分接头变比
变压器有5个抽头，电压调节范围为2*2.5%， 对应的分接头开始时设变压器高压侧接主接头,降压变压器5个分接头时变比K以备调压时选用。

4.3等值电路绘制
根据题目给定的参数和要求，手绘形成了网络的等值电路图,并且在图中标注出了各参数和节点编号。
60+37.18466j

发电厂一
1:k2
1:k3
1
3
5
7
4
6
8

δ=0
U=235
P=200MW
8.5+20.1j
4.25+10.05j
7.65+18.09j
0.05455+1.65024j
0.05455+1.65024j
0.00373+0.10084j
1:k1
50+30.98722j
40+24.78977j
0.000139j
0.00025 j
0.00025j
0.000139j
0.01405+0.03322j
0.00028j
0.00028j
0.00022j
0.00022 j
8.5+20.1j
发电厂二
U=235
60+37.18466j

图4.3 等值电路图
4.4计算方法
利用牛顿拉夫逊法进行求解，用MATLAB软件编程，可以求解系统潮流分
布。根据题目的不同要求对参数进行调整，通过调节变压器变比和发电厂的电压，求解出合理的潮流分布，最后用DDTRS进行潮流分析，将两者进行比较。
五、潮流计算
5.1 给定负荷下的潮流计算
5.1.1 B1/B2矩阵的形成
根据所求参数，以及B1、B2矩阵的含义，列写B1、B2矩阵如下：

1
3
8.5+j20.1
j5.560e-4
1
0
0
1
5
4.25+j10.05
j2.780e-4
1
0
0
B1=
1
7
7.65+j18.09
j5.004e-4
1
0
0
2
9
8.5+j20.1
j5.560e-4
1
0
0
3
4
0.05455+j1.65024
0
231/38.5
1
1
5
6
3.375e-3+j0.10083
0
231/11
1
1
7
8
0.05455+j1.65024
0
231/38.5
1
1
7
9
6.8+j16.08
j4.448e-4
1
0
0
9
10
3.375e-3+j0.10083
0
231/11
1
1

0 0 235 0 1
200 0 235 0 3
0 0 220 0 2
B2=
0 60+j37.18466 35 0 2
0 0 220 0 2
0 50+j30.98722 10 0 2
0 0 220 0 2
0 40+j24.78977 35 0 2
0 0 220 0 2
0 60+j37.18466 10 0 2

5.1.2 潮流调整及分析
将编好的的程序(程序见附录)输入MATLAB中运行，可以得到如下运行结果:
节点
1
2
3
4
5
电压
235.0000
235.0000
230.2417
36.5113
232.9385
节点
6
7
8
9
10
电压
10.7756
232.1691
37.5076
231.3951
10.6310

由表观察到节点4,6,8,10的电压值都较正常范围偏高，因此调节发电厂电压和变压器分接头，经过多次调整，使节点电压均满足要求，调节结果如下：
调整过程：
调节方法
发电厂一
发电厂二
分接头一
分接头二
分接头三
分接头四
未调整
235
235
231
231
231
231
第一次
230
230
231
231
231
231
第二次
225
225
231
231
231
231
第三次
225
225
225.5
231
231
231
调整结果：
电压
1
2
3
4
5
6
7
8
9
10
未调整
235.00
235.00
230.24
36.51
232.93
10.77
232.16
37.50
231.39
10.63
第一次
230.00
230.00
225.04
35.59
227.87
10.52
226.96
36.60
226.15
10.37
第二次
225.00
225.00
219.83
34.66
222.80
10.27
221.75
35.70
220.91
10.11
第三次
225.00
225.00
219.86
35.62
222.80
10.27
221.75
35.70
222.91
10.11

由上表观察到，先后进行发电机机端电压调节和变压器分接头的调节使得节点4,6,8,10的节点电压均在允许范围内。
调节前后各支路损耗对比如下：
功率损耗
未调整
调整后
(1,3)支路
0.712464-28.4049i
0.798149-25.6253i
(1,5)支路
0.251484-14.6238i
0.279321-13.2764i
(1,7)支路
1.90622-22.796i
2.13376-19.9238i
(2,9)支路
6.33378-15.2602i
6.8637-11.41i
(3,4)支路
0.20389+6.1682i
0.21416+6.4789i
(5,6)支路
0.1113+3.0048i
0.12238+3.3037i
(7,8)支路
0.085869+2.5977i
0.094746+2.8663i
(7,9)支路
2.77942-17.3235i
3.02607-14.6339i
(9,10)支路
0.16467+4.4453i
0.18207+4.915i
总损耗
12.5491-82.1924i
13.7144-67.3056i
由此可以看出，随着发电机端电压的降低，网络有功损耗逐渐增大，符合实际情况，而最终调节使得电压在规定范围内。调整后的各支路的功率如下表所示：
各条支路的首端功率Si
各条支路的末端功率Sj
S(1,3)=61.0123+18.0382i
S(3,1)=-60.2142-43.6635i
S(1,5)=50.4017+21.0145i
S(5,1)=-50.1224-34.2909i
S(1,7)=-87.69965+67.51637i
S(7,1)=89.83341-87.4402i
S(2,9)=200-43.72836i
S(9,2)=-193.1363+32.31837i
S(3,4)=60.2142+43.6635i
S(4,3)=-60-37.1847i
S(5,6)=50.1224+34.2909i
S(6,5)=-50-30.9872i
S(7,8)=40.0947+27.656i
S(8,7)=-40-24.7898i
S(7,9)=-129.9282+59.78417i
S(9,7)=132.9542-74.41806i
S(9,10)=60.1821+42.0997i
S(10,9)=-60-37.1847i
支路的末端功率加上支路损耗等于该支路的首段功率，符合实际情况。
且得到各节点的功率 (节点号从小到大排列)：
0.2371+1.0657i 2.0000-0.4373i -0.0000+0.0000i -0.6000-0.3718i -0.0000-0.0000i -0.5000-0.3099i -0.0000+0.0000i -0.4000-0.2479i 0.0000+0.0000i -0.6000-0.3718i
5.1.3 仿真比较
由MATLAB得到电压、功角和具体的潮流分布图，调整好以后的数据如图：

在DDTRS中绘制出电路图，仿真并进行调节以后可以得到潮流分布图：

在用DDTRS绘制的系统潮流分布图中，各个节点电压在调试之后均能满足要求，均达到系统电压等级标准。
将PSAT仿真得到的潮流图中个点电压的有名值与MATLAB计算得到的各节点电压进行比较：
电压
1
2
3
4
5
6
7
8
9
10
DDTRS
224.99
224.99
219.85
35.61
222.79
10.27
221.74
35.70
221.90
10.11
MATLAB
225.00
225.00
219.86
35.62
222.80
10.27
221.75
35.70
222.91
10.11
由上述比较可观察到: 应用DDTRS仿真和MATLAB计算的到的结果大致相同，可以认为此潮流计算是正确的。
5.2 变电所负荷变化时的潮流计算
5.2.1 4个变电所负荷同时以2%的比例增大
由于各个变电所的负荷同时以2%的比例增大，所以此时的各负荷变为： 变电所1 =61.2+j37.92835
变电所2 = 51+j31.60696
变电所3 =40.8+j25.28557
变电所4 =61.2+j37.92835
相应的，B2矩阵变为：

0 0 235 0 1
200 0 235 0 3
0 0 220 0 2
B2=
0 61.2+j37.92835 35 0 2
0 0 220 0 2
0 51+j31.60696 10 0 2
0 0 220 0 2
0 40.8+j25.28557 35 0 2
0 0 220 0 2
0 61.2+j37.92835 10 0 2

将原程序中的B2矩阵替换后，运行并进行调节，经过调节后达到要求，对应的调节方法和结果见下表：
调节方法：
调节方法
发电厂一
发电厂二
分接头一
分接头二
分接头三
分接头四
未调整
235
235
231
231
231
231
调整后
225
225
225.5
231
231
231
调整结果：
电压
1
2
3
4
5
6
7
8
9
10
未调整
235.00
235.00
230.10
36.44
232.88
10.76
232.08
37.46
231.29
10.61
调整后
225.00
225.00
219.71
35.55
222.74
10.26
221.66
35.66
220.80
10.09
由上表观察到，调节后的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：
功率损耗
未调整
调整后
(1,3)支路
0.746602-28.3064i
0.836494-25.5166i
(1,5)支路
0.262717-14.5939i
0.291806-13.2435i
(1,7)支路
1.86419-22.8856i
2.08903-20.0197i
(2,9)支路
6.32203-15.275i
6.8524-11.4236i
(3,4)支路
0.2129+6.4407i
0.22369+6.767i
(5,6)支路
0.116+3.1315i
0.12756+3.4437i
(7,8)支路
0.089529+2.7084i
0.098809+2.9891i
(7,9)支路
2.738-17.4024i
2.98121-14.7207i
(9,10)支路
0.17175+4.6366i
0.18995+5.128i
总损耗
12.5237-81.5461i
13.6909-66.5963i
由电压和功率损耗可以观察到，最后一次的电压和有功损耗均符合题目要求，可认为是合理的。具体的潮流分布如下：
各条支路的首端功率Si
各条支路的末端功率Sj
S(1,3)=62.2602+19.1788i
S(3,1)=-61.4237-44.6954i
S(1,5)=51.4194+21.8072i
S(5,1)=-51.1276-35.0507i
S(1,7)=-85.78861+67.73745i
S(7,1)=87.87764-87.75718i
S(2,9)=200-42.57042i
S(9,2)=-193.1476+31.14682i
S(3,4)=61.4237+44.6954
S(4,3)=-61.2-37.9284i
S(5,6)=51.1276+35.0507i
S(6,5)=-51-31.607i
S(7,8)=40.8988+28.2747i
S(8,7)=-40.8-25.2856i
S(7,9)=-128.7764+59.48247i
S(9,7)=131.7577-74.20318i
S(9,10)=61.39+43.0564i
S(10,9)=-61.2-37.9284i
各节点的功率S为(节点号从小到大排列)：
0.278 +1.0872i 2.0000-0.4257i 0.0000-0.0000i -0.6120-0.3793i 0.0000+0.0000i -0.5100 0.3161i -0.000-0.0000i -0.4080-0.2529i 0.0000+0.0000i -0.6120-0.3793i
5.2.2 4个变电所负荷成比例2%下降
由于各个变电所的负荷同时以2%的比例下降，所以此时的各负荷变为： 变电所1 =58.8+j36.44097
变电所2 = 49+j30.36747
变电所3 =39.2+j24.29398
变电所4 =58.8+j36.44097
相应的，B2矩阵变为：

0 0 235 0 1
200 0 235 0 3
0 0 220 0 2
B2=
0 58.8+j36.44097 35 0 2
0 0 220 0 2
0 49+j30.36747 10 0 2
0 0 220 0 2
0 39.2+j24.29398 35 0 2
0 0 220 0 2
0 58.8+j36.44097 10 0 2

将程序中的B2矩阵替换后，运行并进行调节，此系统经过调节后达到要求，对应的调节方法和结果见下表：
调节方法：
调节方法
发电厂一
发电厂二
分接头一
分接头二
分接头三
分接头四
未调整
235
235
231
231
231
231
调整后
225
225
225.5
231
231
231
调整结果：
电压
1
2
3
4
5
6
7
8
9
10
未调整
235.00
235.00
230.38
36.57
232.98
10.78
232.25
37.54
231.49
10.64
调整后
225.00
225.00
220.01
35.69
222.85
10.28
221.83
35.74
221.01
10.12
由上表观察到，调节后的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：
功率损耗
未调整
调整后
(1,3)支路
0.679291-28.5011i
0.760896-25.7314i
(1,5)支路
0.24052-14.653i
0.267137-13.3086i
(1,7)支路
1.94936-22.7038i
2.17972-19.825i
(2,9)支路
6.34587-15.2445i
6.87539-11.3954i
(3,4)支路
0.19512+5.9027i
0.20489+6.1982i
(5,6)支路
0.10671+2.8808i
0.11731+3.1668i
(7,8)支路
0.082294+2.4896i
0.09078+2.7463i
(7,9)支路
2.82114-17.2439i
3.07126-14.5462i
(9,10)支路
0.15775+4.2586i
0.17437+4.7072i
总损耗
12.5781-82.8145i
13.7417-67.988i
由电压和功率损耗可以观察到，最后一次的电压和有功损耗均符合题目要求，可认为是合理的。具体的潮流分布如下：
各条支路的首端功率Si
各条支路的末端功率Sj
S(1,3)=59.7658+16.9078i
S(3,1)=-59.0049-42.6392i
S(1,5)=49.3844+20.2257i
S(5,1)=-49.1173-33.5343i
S(1,7)=-89.60849+67.30186i
S(7,1)=91.7882-87.12685i
S(2,9)=200-44.88005i
S(9,2)=-193.1246+33.48465i
S(3,4)=59.0049+42.6392i
S(4,3)=-58.8-36.441i
S(5,6)=49.1173+33.5343i
S(6,5)=-49-30.3675i
S(7,8)=39.2908+27.0402i
S(8,7)=-39.2-24.294i
S(7,9)=-131.079+60.08661i
S(9,7)=134.1502-74.63283i
S(9,10)=58.9744+41.1482i
S(10,9)=-58.8-36.441i
各节点的功率S为(节点号从小到大排列)：
0.1954+1.0444i 2.0000-0.4488i -0.0000+0.0000i -0.5880-0.3644i 0.0000+0.0000i -0.4900-0.3037i -0.0000+0.0000i -0.3920-0.2429i -0.0000-0.0000i - 0.5880-0.3644i
5.2.3 1、4号变电所负荷成2%下降，2、3号变电负荷成2%上升
此时的各负荷变为：
变电所1 =58.8+j36.44097
变电所2 = 51+j31.60696
变电所3 =40.8+j25.28557
变电所4 =58.8+j36.44097
相应的，B2矩阵变为：

0 0 235 0 1
200 0 235 0 3
0 0 220 0 2
B2=
0 58.8+j36.44097 35 0 2
0 0 220 0 2
0 51+j31.60696 10 0 2
0 0 220 0 2
0 40.8+j25.28557 35 0 2
0 0 220 0 2
0 58.8+j36.44097 10 0 2

将程序中的B2矩阵替换后，运行并进行调节，此系统经过调节后达到要求，对应的调节方法和结果见下表：
调节方法：
调节方法
发电厂一
发电厂二
分接头一
分接头二
分接头三
分接头四
未调整
235
235
231
231
231
231
调整后
225
225
225.5
231
231
231
调整结果：
电压
1
2
3
4
5
6
7
8
9
10
未调整
235.00
235.00
230.38
36.57
232.88
10.76
232.15
37.48
231.44
10.64
调整后
225.00
225.00
220.01
35.69
222.74
10.26
221.73
35.67
220.96
10.12
由上表观察到，调节后的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：

功率损耗
未调整
调整后
(1,3)支路
0.679291-28.5011i
0.760896-25.7314i
(1,5)支路
0.262717-14.5939i
0.291806-13.2435i
(1,7)支路
1.92183-22.7577i
2.15108-19.8814i
(2,9)支路
6.33938-15.253i
6.8691-11.4033i
(3,4)支路
0.19512+5.9027i
0.20489+6.1982i
(5,6)支路
0.116+3.1315i
0.12756+3.4437i
(7,8)支路
0.08947+2.7066i
0.098736+2.987i
(7,9)支路
2.81284-17.248i
3.06221-14.552i
(9,10)支路
0.15783+4.2608i
0.17446+4.7098i
总损耗
12.5745-82.352i
13.7407-67.4728i

由电压和功率损耗可以观察到，最后一次的电压和有功损耗均符合题目要求，可认为是合理的。具体的潮流分布如下：

各条支路的首端功率Si
各条支路的末端功率Sj
S(1,3)=59.7658+16.9078i
S(3,1)=-59.0049-42.6392i
S(1,5)=51.4194+21.8072i
S(5,1)=-51.1276-35.0507i
S(1,7)=-88.04441+67.85231i
S(7,1)=90.19549-87.7337i
S(2,9)=200-44.26567i
S(9,2)=-193.1309+32.86236i
S(3,4)=59.0049+42.6392i
S(4,3)=-58.8-36.441i
S(5,6)=51.1276+35.0507i
S(6,5)=-51-31.607i
S(7,8)=40.8987+28.2725i
S(8,7)=-40.8-25.2856i
S(7,9)=-131.0942+59.46118i
S(9,7)=134.1564-74.01315i
S(9,10