电力系统分析课程设计-论文报告.doc

枯隅莱若逝衬拓椽浆詹蓟赦赊扶蛮策珠作伊郊孩浆音篆婴造吻鲍俺寸搬伶既历券骸宪约招谚圣叭昂哉万洲关蛔勤兴陇喷霄艾描沸玖锄刷苔酿规焚纬嚼叹棠卉墓色龟颤颂货睦灾隅折须瞥殿席彬玄菜燎抠绣细插简恕廉历吝绣躺秩蚌枕始靴晕芳沾颜惋雍甸条榴排赦铁赡腆几并避躇汁胺泄膀柒环筏眷退来隙手名刮悟剐赎站揖颇损励次思夜鹤铬剁率万毖翅麻破处帐婴镍墟隆葫刃绒魂炔胰沏颈弯吕约讽碱吮慷藻足迁载踌码削赣颖丝藩涵滁坞烽脂磕含令春函蛇噎霉汝盅忙表吗售浆澜崖廊狮肠橙菠沦蝴孪菲咽煽绞守锻债屠行通纫缩力住今柒讹铡闷抬车老淤明画祥古虞跳自嚼稠巴极锋胰限芭闯梆
----------------------------精品word文档 值得下载 值得拥有----------------------------------------------
----------------------------------------------------------------------------------------------------------------------------------------------迁檀察妈战辉拉玩促们幻淘幻肿啃叛寨案若痛罢鄂陆斑怨貉特拼景遭洋可侮雀迪捶夜夸割除揪尸梗侠请犹芝慨峰掂髓徘担仕佛坝铲窖父锋缀世绒红庸磐巳钾瞎状纠垫下竣笋狰牌伺吵形勿迄家演设旺狐据恩鸥焊矗节占枉秘孪卵纬梁溺岿迭手堡贡试感驮茹日凿寒瞄酿欲煤托野勇致衡丝挣畸系卿戍堰她皿脖逗逝秀芍使胃即冰肛脑策块褥蚀闽缨花湃捎凑酱瓶秩奢涉误啃聪吁晓采虫苦斩授秽鸣柜处胁盯歧叠穗芝喷潍孽淀占氯祝漏渴桓匡池绚铱扰桩痕秤耘揽突囱替庇讥雀你殊拾和膘刚伎堵睁贰苞弄套萤捉闭湘稍峙渔殆击舱概青桥屹拐氛吱槐勉无沃粟涝笔缨寿卑疚木洼吮抹酒漓馒伦度陷厅互电力系统分析课程设计_论文报告救贪幼澎天时赚重糙埂刀做权纶垮弦贬缓铝逢役摧杭洽虏仟突袁羚祈钓傈碑导厌翠撂敌盗丝蒙冻掂乡伟唁下墅腑痊督惰豺屎肚癌捕解味新脚裴移鸳邪件椭恫浩祷耙嚷舵收胞更脊翔挝赃召慢础鞠帛谅从卷习墨吐拆掷送锨篱锰角槽函勘尿篷议喳鬃朱宰炮遍减剔娃逸幅伟悦椭也吉瞒祁离嗅砰宜俏淹吗打椎躯割特猜堂圾寂襟国驯慈辰侨母场菜炊劫荣车炔残掇短地劝鳖蛋吁够岔恒拣些辨很役坯咀恩功囱眶阑巢蹄卷淀樟方穴狼箭计旱辖涂雹随记烙醉罢护面祝筐嘲摄啡牵绑壳郑檄床彪蛹吃胯毯砸归挣跪峻契哑墨泌歉郝伸渍株勃忌探次售莽霄维森治淆查短罢陨卫账虾尸费护埋狄铝切厄任癌霸焊
摘 要
本文运用MATLAB软件进行潮流计算，对给定题目进行分析计算，再应用DDRTS软件，构建系统图进行仿真，最终得到合理的系统潮流。
潮流计算是电力系统最基本最常用的计算。根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压幅值和相角，各元件流过的功率，整个系统的功率损耗。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。因此，潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有着广泛的应用。
首先，画出系统的等效电路图，在计算出各元件参数的基础上，应用牛顿—拉夫逊Newton-Raphson法以及MATLAB软件进行计算对给定系统图进行了四种不同负荷下的潮流计算，经过调节均得到符合电压限制及功率限制的潮流分布。
其次，轮流断开环网的三条支路，在新的系统结构下进行次潮流计算，结果亦均满足潮流分布要求。牛顿—拉夫逊Newton-Raphson法具有较好的收敛性，上诉计算过程经过四到五次迭代后均能收敛。
最后，应用DDRTS软件，构建系统图，对给定负荷重新进行分析，潮流计算后的结果也能满足相应的参数要求。

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

一、 题目原始资料
1．系统图：
两个发电厂分别通过变压器和输电线路与四个变电所相连。
变电所4
变电所3
变电所2
变电所1

10kV母线
35kV母线
10kV母线
35kV母线

一次侧电压220kV
一次侧电压220kV

线路长为90km

线路长为80km

线路长为90km

线路长为90km
线路长为100km

线路长为80km

母线2
母线1

。。。。。。。。。。。。。

母线3

2*QFQ-50-2

2*TQN-100-2
TQN-100-2
2*QFs-50-2

电厂二
电厂一

２．发电厂资料：
母线1和2为发电厂高压母线，发电厂一总装机容量为300MW，母线3为机压母线，机压母线上装机容量为100MMW，最大负荷和最小负荷分别为40MW和20MW；发电厂二总装机容量为200WM。
３．变电所资料：
(一) 变电所1、2、3、4低压母线的电压等级分别为：10KV 35KV 35KV 10KV
(二) 变电所的负荷分别为：（3） 50MW 40MW 55MW 70MW
(三) 每个变电所的功率因数均为cosφ=0.85；
(四) 变电所1和变电所2分别配有两台容量为75MVA的变压器，短路损耗414KW，短路电压（%）=16.7；变电所3和变电所4分别配有两台容量为63MVA的变压器，短路损耗为245KW，短路电压（%）=10.5；
４．输电线路资料：
发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为，单位长度的电抗为，单位长度的电纳为。

二、 课程设计内容及要求
课设内容：
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. 最终形成课程设计成品说明书。
课设要求：
1. 在读懂程序的基础上画出潮流计算基本流程图
2. 通过输入数据，进行潮流计算输出结果
3. 对不同的负荷变化，分析潮流分布，写出分析说明。
4. 对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并列表表示调节控制的参数变化。
5. 打印利用DDRTS进行潮流分析绘制的系统图，以及潮流分布图。

三、 题目分析

1． 节点设置及分类
根据系统图可知此系统为两端供电网路，将母线1，2设为节点1，2，将变
电所1、2、3、4的高低压侧分别设为节点3、4、5、6、7、8、9、10。并且，将节点1设为平衡节点，将节点2设为PV节点，其余节点设为PQ节点。
2． 参数求取
设定基准值，，所以根据题目原始资料，计算发电厂、变压器及线路的参数。
（1）运用下列公式计算变压器参数：

（2）计算线路参数

（3）变电所负荷分别为：

变电所1 =0.5047+j0.31278 变电所2 =0.6489+j0.40215
变电所3 =0.4326+j0.2681 变电所4 =0.5768+j0.35747
将参数整理，见下表：

首端号
末端号
阻抗有名值
阻抗标幺值
电纳有名值
电纳标幺值
1
2
5.1+j19.2
0.01054+j0.03967
j0.00084
j0.4
1
4
3.4+j12.8
0.00702+j0.02645
j0.00056
j0.267
1
9
2.975+j11.2
0.00615+j0.02314
j0.00049
j0.233
2
6
1.78+j53.885
0.0368+j0.1115
0
0
2
3
4.25+j16
0.00878+j0.03306
j0.0007
j0.333
3
7
2.46+j70.18
0.0051+j0.145
0
0
3
10
2.55+j9.6
0.00527+j0.01983
j0.00046
j0.219
4
8
1.5+j48
0.0031+j0.099
0
0
5
9
1.5+j48
0.0031+j0.099
0
0
（4）计算变压器分接头变比
变压器有5个抽头，电压调节范围为2*2.5%， 对应的分接头开始时设变压器高压侧接主接头,降压变压器5个分接头时的非标准变比以备调压时选用

对变电所低压母线为35KV时，非标准变比与10KV时相同。

3． 求解方法

利用牛顿拉夫逊法进行求解，用MATLAB软件编程，可以求解系统潮流分
布根据题目的不同要求对参数进行调整，通过调节变压器变比和发电厂的电压，求解出合理的潮流分布，最后用DDRTS进行潮流分析，将两者进行比较。

四、 题目求解

(一) 依据题目要求及原始资料画出等值电路图：

(二) 读程序画出牛顿拉夫逊法的流程图
启动

输入原始数据

形成节点导纳矩阵

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

是

否
是

(三) 变电所负荷为题目所给数据进行求解

1． 形成B1、B2矩阵：
根据所求参数，以及B1矩阵的含义，列写B1矩阵如下:
B1=[ 1 2 0.01054+0.03967i 0.4i 1 0;
1 4 0.00702+0.02645i 0.267i 1 0;
1 9 0.0065+0.02314i 0.233i 1 0;
2 6 0.00368+0.1115i 0 0.909 1 ;
2 3 0.00878+ 0.03306i 0.333i 1 0;
3 7 0.0051+0.145i 0 0.909 1;
3 10 0.00527+0.01983i 0.219i 1 0;
4 8 0.0031+0.099i 0 0.909 1;
5 9 0.0031+0.099i 0 0.909 0 ]

由各个变电所负荷功率可以计算出总功率为215MW，而发电厂一、二的总装机容量分别为300 MW和200 MW。为了平衡两发电厂发出的电量，令发电厂二的功率为120 MW，为了减小线路上的损耗，令发电机的电压为额定电压的1.05倍。并且，根据前面叙述的节点分类，形成B2矩阵如下：

B2=[ 0 0 1.05 1.05 0 1;
0 0 1 0 0 2;
0 0 1 0 0 2;
0 0 1 0 0 2;
0 0.5047+0.31278i 1 0 0 2;
0 0.6489+0.40215i 1 0 0 2;
0 0.4326+0.2681i 1 0 0 2;
0 0.5768+0.35747i 1 0 0 2;
0 0 1 0 0 2;
1.2 0 1.05 1.05 0 3 ]

由于各节点电压不应超过各自的允许范围，由题目要求知，发电厂电压标幺值在1到1.05变化时，变压器低压侧电压标幺值变化范围如下表：

节点
4
6
8
10
电压下限
1
1
1
1
电压上限
1.05
1.05
1.05
1.05

同时保证发电厂发出的有功功率要大于变电所有功功率，并且两者差值越小表明线路有功功率损耗越小，越满足实际要求。
2． 编写程序并运行
由于程序较复杂，故将其放入附录中，具体程序和结果分别见附录一，二。
由附录2的结果可以观察到各个节点电压标幺值分别为:

节点
1
2
3
4
5
6
7
8
9
10
电压
1.0500
1.0437
1.0460
1.0397
1.1164
1.1036
1.1125
1.1091
1.0420
1.0500

由此观察到节点4、6、8、10的电压都较正常范围偏高，因此调节变压器分接头和发电厂电压，此系统共进行四次调整，最终得到合理结果，各次调整和结果如下:
调节方法：
调节方法
电厂1电压
电厂2电压
分接头1
分接头2
分接头3
分接头4
未调整
1.05
1.05
0.909
0.909
0.909
0.909
第一次
1.05
1.05
0.955
0.955
0.955
0.955
第二次
1.04
1.04
0.955
0.955
0.955
0.955
第三次
1.035
1.04
0.955
0.955
0.955
0.955
第四次
1.025
1.025
0.955
0.955
0.955
0.955

调节结果：
节点电压
1
2
3
4
5
6
7
8
9
10
未调整
1.05
1.0437
1.0460
1.0397
1.1164
1.1036
1.1125
1.1091
1.0420
1.0500
第一次
1.05
1.0439
1.0462
1.0399
1.0603
1.0469
1.0559
1.0530
1.0422
1.0500
第二次
1.04
1.0336
1.0360
1.0297
1.0493
1.0355
1.0447
1.0419
1.0321
1.0400
第三次
1.035
1.0307
1.0349
1.0246
1.0439
1.0323
1.0435
1.0364
1.0270
1.0400
第四次
1.025
1.0177
1.0204
1.0140
1.0319
1.0170
1.0268
1.0241
1.0166
1.0250

由上表观察到，进行最后一次调节时节点4，6，8，10的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗的标幺值，记录于下表：
功率损耗
未调整
第一次
第二次
第三次
第四次
(1,2)支路
0.0005
0.0005
0.0005
0.0004
0.0006
(1,4)支路
0.0026
0.0025
0.0025
0.0026
0.0028
(1,9)支路
0.0018
0.0017
0.0018
0.0018
0.0019
(2,6)支路
0.0018
0.0019
0.0019
0.0019
0.0021
(2,3)支路
0.0047
0.0049
0.0049
0.0049
0.0050
(3,7)支路
0.0011
0.0011
0.0012
0.0012
0.0013
(3,10)支路
0.0069
0.0069
0.0071
0.0070
0.0072
(4,8)支路
0.0012
0.0012
0.0013
0.0010
0.0014
(5,9)支路
0.0009
0.0009
0.0010
0.0013
0.0010
总损耗
0.0215
0.0216
0.0222
0.0221
0.0233
由此可以看出，随着变压器分接头的调高以及发电厂电压的降低，有功损耗逐渐增大，符合实际情况，而最终调节使得电压在规定范围内时，有功损耗增加了0.3MW，相对与实际损耗42.1MW小很多，可以认为是合理的。具体的潮流分布如下:

各条支路的首端功率Si
各条支路的末端功率Sj
S(1,2)= -10.24+0.57i
S(2,1)=10.30-42.08i
S(1,4)=58.09+13.38i
S(4,1)=-57.82-40.09i
S(1,9)=50.77+10.96i
S(9,1)=-50.57-34.56i
S(2,6)=65.10+46.50i
S(6,2)=-64.89-40.21i
S(2,3)= -75.39-4.42i
S(3,2)=75.89-28.30i
S(3,7)=43.39+30.37i
S(7,3)=50.57+34.56i
S(3,10)=-119.28-2.08i
S(10,3)=120.00-18.10i
S(4,8)=57.82+40.09i
S(8,4)=-57.68-35.75i
S(5,9)= -50.47-31.28i
S(9,5)=-43.26-26.81i

各节点的功率S为(节点号从小到大排列)：各节点的功率S为(节点号从小到大排列)：
各节点的功率S为(节点号从小到大排列)：
Columns 1 through 4

0.9862 + 0.2491i 0.0000 - 0.0000i -0.0000 + 0.0000i -0.0000 - 0.0000i

Columns 5 through 8

-0.5047 - 0.3128i -0.6489 - 0.4022i -0.4326 - 0.2681i -0.5768 - 0.3575i

Columns 9 through 10

0 1.2000 - 0.1810i
得到的图像如下：

(四) 4个变电所负荷同时以2%的比例增大
由于变电所负荷变大，所以各个变电所符合标幺值变为：

变电所1 =0.51+j0.3.162
变电所2 =0.49+j0.3038
变电所3 =0.561+j0.348
变电所4 =0.714+j0.0.443
相应地， B2矩阵变为：

将程序中的B2矩阵替换后，运行并进行调节，此系统经过四次调节后达到要求，对应的调节方法和结果见下表：
调节方法：
调节方法
电厂1电压
电厂2电压
分接头1
分接头2
分接头3
分接头4
未调整
1.05
1.05
0.909
0.909
0.909
0.909
第一次
1.05
1.05
0.932
0.932
0.932
0.932
第二次
1.05
1.05
0.955
0.955
0.955
0.932
第三次
1.03
1.03
0.955
0.955
0.955
0.932
第四次
1.025
1.025
0.955
0.955
0.955
0.932

调节结果：

节点电压
1
2
3
4
5
6
7
8
9
10
未调整
1.05
1.05
1.02
1.09
1.03
1.11
1.03
1.11
1.03
1.09
第一次
1.05
1.05
1.02
1.06
1.03
1.08
1.03
1.08
1.03
1.06
第二次
1.05
1.05
1.02
1.03
1.03
1.05
1.03
1.05
1.03
1.06
第三次
1.03
1.03
1
1.01
1.01
1.03
1.01
1.03
1.01
1.04
第四次
1.03
1.03
1
1.01
1.01
1.02
1.01
1.02
1
1.04

功率损耗
未调整
第一次
第二次
第三次
第四次
(1,3)支路
0.0105
0.0105
0.0105
0.011
0.0111
(1,5)支路
0.013
0.013
0.013
0.0135
0.0137
(2,7)支路
0.0027
0.0027
0.0028
0.003
0.0031
(2,9)支路
0.0087
0.0087
0.0087
0.0092
0.0093
(3,4)支路
0.0011
0.0012
0.0012
0.0013
0.0013
(3,5)支路
0.0006
0.0006
0.0006
0.0006
0.0006
(5,6)支路
0.0007
0.0007
0.0008
0.0008
0.0008
(5,7)支路
0.0016
0.0016
0.0016
0.0017
0.0017
(7,8)支路
0.0011
0.0012
0.0012
0.0013
0.0013
(9,10)支路
0.0018
0.0019
0.0019
0.002
0.002
总损耗
0.0418
0.0421
0.0423
0.0444
0.0449

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

各条支路的首端功率Si
各条支路的末端功率Sj
S(1,3)=56.3705+5.38448i
S(3,1)=-55.2571-16.5016i
S(1,5)=67.4239-5.77691i
S(5,1)=-66.0565-3.46334i
S(2,7)=27.4697+23.4309i
S(7,2)=-27.1597-47.7144i
S(2,9)=72.5303+29.8637i
S(9,2)=-71.6033-49.7819i
S(3,4)=51.1312+35.5759i
S(4,3)=-51-31.62i
S(3,5)=4.12593-19.0743i
S(5,3)=-4.0662+7.088i
S(5,6)=40.8802+27.2176i
S(6,5)=-40.8-24.8i
S(5,7)=29.2425-30.8423i
S(7,5)=-29.0687+9.45298i
S(7,8)=56.2284+38.2614i
S(8,7)=-56.1-34.8i
S(9,10)=71.6033+49.781
S(10,9)=-71.4-44.3i

各节点的功率S为(节点号从小到大排列)：

Columns 1 through 5

1.2374 - 0.0096i 1.0000 + 0.5247i -0.0000 + 0.0000i -0.5100 - 0.3162i 0.0000 - 0.0000i

Columns 6 through 10

-0.4080 - 0.2480i 0.0000 - 0.0000i -0.5610 - 0.3480i -0.0000 + 0.0000i -0.7140 - 0.4430i

得到的图像如下：

理为
(五) 4个变电所负荷同时以2%的比例下降

由于变电所负荷变小，所以各个变电所符合标幺值变为：
变电所1 =0.49+j0.3.038
变电所2 =0.392+j0.243
变电所3 =0.539+j0.334
变电所4 =0.686+j0.425
相应地， B2矩阵变为：

将程序中的B2矩阵替换后，运行并进行调节，此系统经过三次调节后达到要求，对应的调节方法和结果见下表：
调节方法：
调节方法
电厂1电压
电厂2电压
分接头1
分接头2
分接头3
分接头4
未调整
1.05
1.05
0.909
0.909
0.909
0.909
第一次
1.05
1.05
0.932
0.932
0.932
0.932
第二次
1.05
1.05
0.955
0.955
0.955
0.932
第三次
1.02
1.02
0.955
0.955
0.955
0.932

调节结果：
节点电压
1
2
3
4
5
6
7
8
9
10
未调整
1.05
1.05
1.0249
1.094
1.0328
1.1099
1.035
1.1114
1.029
1.0967
第一次
1.05
1.05
1.0248
1.0651
1.0327
1.081
1.0349
1.0824
1.0289
1.0676
第二次
1.05
1.05
1.0247
1.0375
1.0326
1.0535
1.0348
1.0549
1.0289
1.0676
第三次
1.02
1.02
0.9929
1.003
1.0011
1.0196
1.0036
1.0211
0.9979
1.0331

功率损耗
未调整
第一次
第二次
第三次
(1,3)支路
0.0092
0.0092
0.0092
0.0092
(1,5)支路
0.0109
0.0109
0.0109
0.0109
(2,7)支路
0.0026
0.0026
0.0026
0.0026
(2,9)支路
0.0079
0.008
0.008
0.008
(3,4)支路
0.001
0.0011
0.0011
0.0011
(3,5)支路
0.0004
0.0004
0.0004
0.0004
(5,6)支路
0.0006
0.0007
0.0007
0.0007
(5,7)支路
0.0011
0.0011
0.0011
0.0011
(7,8)支路
0.001
0.0011
0.0011
0.0011
(9,10)支路
0.0017
0.0018
0.0018
0.0018
总损耗
0.0364
0.0369
0.0369
0.0369

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

各条支路的首端功率Si
各条支路的末端功率Sj
S(1,3)=52.6693+2.93801i
S(3,1)=-51.752-15.2555i
S(1,5)=61.6782-7.19992i
S(5,1)=-60.5878-3.35268i
S(2,7)=30.4398+16.0066i
S(7,2)=-30.1835-41.7232i
S(2,9)=69.5602+25.6238i
S(9,2)=-68.7673-47.0105i
S(3,4)=49.1022+33.4628i
S(4,3)=-49-30.38i
S(3,5)=2.6498-18.2074i
S(5,3)=-2.6055+5.4938i
S(5,6)=39.2635+26.2167i
S(6,5)=-39.2-24.3i
S(5,7)=23.9297-28.3578i
S(7,5)=-23.8171+5.61161i
S(7,8)=54.0006+36.1116i
S(8,7)=-53.9-33.4i
S(9,10)=68.7673+47.0105i
S(10,9)=-68.6-42.5i
各节点的功率S为(节点号从小到大排列)：
Columns 1 through 5

1.1468 - 0.0046i 1.0000 + 0.4711i -0.0000 - 0.0000i -0.4900 - 0.3038i -0.0000 - 0.0000i

Columns 6 through 10

-0.3920 - 0.2430i 0.0000 - 0.0000i -0.5390 - 0.3340i -0.0000 - 0.0000i -0.6860 - 0.4250i

得到的图像如下：

(六) 1、4变电所负荷以2%的比例下降2、3变电所负荷以2%的比例增加

各个变电所符合标幺值变为：
变电所1 =0.49+j0.3.038
变电所2 =0.49+j0.3038
变电所3 =0.561+j0.348
变电所4 =0.686+j0.425
相应地， B2矩阵变为：

将程序中的B2矩阵替换后，运行并进行调节，此系统经过三次调节后达到要求，对应的调节方法和结果见下表：
调节方法：
调节方法
电厂1电压
电厂2电压
分接头1
分接头2
分接头3
分接头4
未调整
1.05
1.05
0.909
0.909
0.909
0.909
第一次
1.05
1.05
0.955
0.955
0.955
0.955
第二次
1.03
1.03
0.955
0.955
0.955
0.955
第三次
1.02
1.02
0.955
0.955
0.955
0.955

调节结果：
节点电压
1
2
3
4
5
6
7
8
9
10
未调整
1.05
1.05
1.0244
1.0934
1.0319
1.1083
1.0341
1.1092
1.029
1.0967
第一次
1.05
1.05
1.0242
1.037
1.0317
1.0519
1.0339
1.0526
1.0288
1.0399
第二次
1.03
1.03
1.003
1.0139
1.0107
1.0293
1.0131
1.03
1.0082
1.0174
第三次
1.02
1.02
0.9924
1.0024
1.0002
1.0179
1.0026
1.0188
0.9978
1.0061

功率损耗
未调整
第一次
第二次
第三次
(1,3)支路
0.0096
0.0097
0.0101
0.0104
(1,5)支路
0.012
0.012
0.0125
0.0128
(2,7)支路
0.0028
0.0028
0.0031
0.0032
(2,9)支路
0.0101
0.0102
0.0107
0.011
(3,4)支路
0.001
0.0011
0.0012
0.0012
(3,5)支路
0.0005
0.0005
0.0005
0.0005
(5,6)支路
0.0007
0.0008
0.0008
0.0008
(5,7)支路
0.0013
0.0013
0.0014
0.0015
(7,8)支路
0.0011
0.0012
0.0013
0.0013
(9,10)支路
0.0022
0.0024
0.0025
0.0026
总损耗
0.0413
0.042
0.0441
0.0453

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

各条支路的首端功率Si
各条支路的末端功率Sj
S(1,3)=54.1893+4.94198i
S(3,1)=-53.1504-16.113i
S(1,5)=64.8352-5.16571i
S(5,1)=-63.5581-4.16816i
S(2,7)=30.0421+22.1539i
S(7,2)=-29.7217-46.1681i
S(2,9)=69.9579+46.0785i
S(9,2)=-68.8578-65.2603i
S(3,4)=49.1217+34.052i
S(4,3)=-49-30.38i
S(3,5)=4.02871-17.939i
S(5,3)=-3.9775+6.0392i
S(5,6)=40.881+27.2422i
S(6,5)=-40.8-24.8i
S(5,7)=26.6546-29.1133i
S(7,5)=-26.5081+7.87004i
S(7,8)=56.2298+38.298i
S(8,7)=-56.1-34.8i
S(9,10)=68.8578+65.2603i
S(10,9)=-68.6-58.31i

各节点的功率S为(节点号从小到大排列)：
Columns 1 through 5

1.1871 - 0.0012i 1.0000 + 0.5000i 0.0000 - 0.0000i -0.4900 - 0.3038i -0.0000 + 0.0000i

Columns 6 through 10

-0.4080 - 0.2480i -0.0000 - 0.0000i -0.5610 - 0.3480i 0.0000 + 0.0000i -0.6860 - 0.4250i

得到的图像如下：

（七）1、3节点短线时潮流计算

此时支路数目为9，B1矩阵变为

将程序中的B1矩阵替换后，运行并进行调节，此系统经过三次调节后达到要求，对应的调节方法和结果见下表：
调节方法：
调节方法
电厂1电压
电厂2电压
分接头1
分接头2
分接头3
分接头4
未调整
1.05
1.05
0.909
0.909
0.909
0.909
第一次
1.05
1.05
0.932
0.932
0.932
0.932
第二次
1.05
1.05
0.955
0.955
0.955
0.955
第三次
1.03
1.03
0.955
0.955
0.955
0.955

调节结果：
节点电压
1
2
3
4
5
6
7
8
9
10
未调整
1.05
1.05
0.9762
1.0377
1.0149
1.0892
1.0252
1.09