吴青峰语录 1.曾经得到和失去的，只要我们愿意，都还来得及的，来得及对我们生命做一个交代。我们 可以也应该把自己的灵魂装在皮箱里面到处流浪， 我们应该多了解自己一点。 我们不应该一 味想着别人怎么看待我们，别人希望我们怎么样，我们应该自己对自己负责，不能让自己的 生命一点都不值得。我们都应该飞翔，发光，即便是断翅和黑暗的，勇敢爱自己，做自己的 样子都是飞翔和发光。2.看穿或厌倦了太多虚伪的轻薄风景， 请让自己回到自己安静的内心， 诚实接受自己的样子。3.我只想唱来一些温暖，让我们心灵不会腐坏。4.别害怕去喜欢你自己的样子。5.有时候我们在生活中，面对不开心的事常常不去解决，而是搁着等待时间风化，或是自己 骗自己遗忘甚至闪躲，但事实上一旦发现了才知道原来从不曾真正放下。6.我们都想为对方发光，为对方飞翔，即使自己不能像风筝般飞翔，我们也有彼此，不用苦 苦追赶。7.如果你选错了时机和舞台，就别怪别人忽略你。8.在这样的城市，我还是愿意微笑并真诚，我宁愿傻傻地相信自己可以融化那些孤独岛。9.我们最需要的有时候不是快乐，不是梦想，而是相信事情，相信自己的勇气。敢勇敢地相 信自己的勇气，敢勇敢地相信某些信念，是最坚强的。10.不管我们身在多么纷扰的世界当中，我们还得为自己找到相信的东西和理由，不然我们 的生命将会一点意义也没有。11.请让我在你身边，一起记念这一天12.为什么喜欢苏打绿的多数是年轻人，因为他们的心单纯、开放，对我们没有成见，可以 接受不一样的东西。13.没有想到你下手如此狠毒...XD14. 如果有不切实际的梦想，才是实际的根本。15.有时候我会觉得，你们好象都有跟我们一样的梦想，而我们只是代表你们来实现这个梦 想的。所以，每次只要想到你们，就觉得你们一定站在我们身边和我们连成一片，谢谢你们 给我们的这些力量！如果你们有一些梦想，苏打绿今天站在这边，我们都可以做到，你们一 定也可以。 这些话我讲了很多遍，但是今天还是想要再说一遍。 请你一定要相信自己， 一定要接受喜欢自己的样子， 一定要让自己变成你真心会喜欢的 样子。如果你想要做的不是长辈所控制你的样子，不是社会所规定你的样子，请你一定要为 自己勇敢的站出来，温柔的推翻这个世界，然后把世界变成我们的。16.我会让自己站的好好的，就算摇摇晃晃也绝不倒下。17.如果我们不能理解自己虚构的世界,就无法理解真实的世界。18.是你浪费在我身上的时间、 使我变得如此珍贵。19. 不被他人知道， 也不被

自己了解的落寞， 连落寞本身想必也因为连主人都无法理解而落 寞吧。

 

第二篇：031140635向青峰

信息工程学院课程设计（论文）

信息工程学院

课程设计报告书

题目: 电力系统潮流计算

专 业：电气工程及其自动化

班 级： 0311406

学 号： 031140635

学生姓名： 向青峰

指导教师： 耿东山

2014 年 5 月 30 日

1

信息工程学院课程设计（论文）

信息工程学院课程设计任务书

年 月 日

2

信息工程学院课程设计（论文）

信息工程学院课程设计成绩评定表

3

信息工程学院课程设计（论文）

摘 要

在如今的社会，电力已经成为人们必不可少的需求，而建立结构合理的大型电力系统不仅便于电能生产与消费的集中管理、统一调度和分配，减少总装机容量，节省动力设施投资，且有利于地区能源资源的合理开发利用，更大限度地满足地区国民经济日益增长的用电需要。电力系统建设往往是国家及地区国民经济发展规划的重要组成部分。

电力系统的出现，使高效、无污染、使用方便、易于调控的电能得到广泛应用，推动了社会生产各个领域的变化，开创了电力时代，发生了第二次技术革命。电力系统的规模和技术水准已成为一个国家经济发展水平的标志之一。

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

潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。是电力系统研究人员长期研究的一个课题。它既是对电力系统规划设计和运行方式的合理性、 可靠性及经济性进行定量分析的依据 ， 又是电力系统静态和暂态稳定计算的基础。

关键词：潮流计算，等值电路，MATLAB仿真，电力系统分析

4

信息工程学院课程设计（论文）

目 录

1 任务提出与方案论证 ....................................................... 6

2 总体设计 ................................................................. 7

2.1潮流计算等值电路 .................................................... 7

2.2建立电力系统模型 .................................................... 7

2.3模型的调试与运行 .................................................... 7

3 详细设计 ................................................................. 8

3.1 计算前提 ........................................................... 8

3.2手工计算 ........................................................... 11

4设计图及源程序 .......................................................... 15

4.1MATLAB仿真 ......................................................... 15

4.2潮流计算源程序 ..................................................... 15

5 总结 .................................................................... 23

参考文献 .................................................................. 23

5

信息工程学院课程设计（论文）

1 任务提出与方案论证

潮流计算是在给定电力系统网络结构、参数和决定系统运行状态的边界条件的情况下确定系统稳态运行状态的一种基本方法,是电力系统规划和运营中不可缺少的一个重要组成部分。可以说,它是电力系统分析中最基本、 最重要的计算,是系统安全、 经济分析和实时控制与调度的基础。常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。是电力系统研究人员长期研究的一个课题。它既是对电力系统规划设计和运行方式的合理性、 可靠性及经济性进行定量分析的依据 ， 又是电力系统静态和暂态稳定计算的基础。

潮流计算经历了一个由手工到应用数字电子计算机的发展过程，现在的潮流算法都以计算机的应用为前提用计算机进行潮流计算主要步骤在于编制计算机程序，这是一项非常复杂的工作。对系统进行潮流分析，本文利用 MATLAB中的SimpowerSystems工具箱设计电力系统，在Simulink 环境下，不仅可以仿真系统的动态过程，还可以对系统进行稳态潮流分析。

6

信息工程学院课程设计（论文）

2 总体设计

SimpowerSystems使用Simulink环境，可以将该系统中的发电机、变压器，线路等模型联结起来，形成电力系统仿真模拟图。在加人测量模块，并对各元件的参数进行设置后，用measurement和sink中的仪器可以观察各元件的电压、电流、功率的大小。

2.1潮流计算等值电路

YN,d11kV

,d11?10MWVAkV

YN,2?4?15MW"?xd?0.136?x2?0.16?x?0.073?0

?cos?N?0.8

3?12MW

"?xd?0.128?x2?0.154?x?0.054?0

?cos?N?0.85

"?xd?0.128?x2?0.157?x?0.0591?0

?cos?N?0.8

?x?0.136?x2?0.161?x?0.075?0

?cos?N?0.8

"d

2.2建立电力系统模型

在Simulink中按照电力系统原型选择元件进行建模。所建立的模型和建立的方法在详细设计中详述。

在电力系统模型的建立工程中主要涉及到的是：元器件的选择及其参数的设置；发电机选型；变压器选择；线路的选择；负荷模型的选择；母线选择。

2.3模型的调试与运行

建立系统模型，并设置好参数以后，就可以在Simulink环境下进行仿真运行。运行的具体结果和分析也在详细设计中详述。

7

kV

信息工程学院课程设计（论文）

3 详细设计

3.1 计算前提

首先是发电机的参数计算，先对5个发电厂简化为5台发电机来计算。 发电机G1：

P1?4?15?60MW

Q1?60?tan(arccos0.8)?45MVar

发电机G2：

P2?4?63?252MW

Q2?252?tan(arccos0.85)?156MVar

发电机G3：

P3?3?12?36MW

Q3?36?tan(arccos0.8)?27MVar

发电机G4：

P4?1?50?50MW

Q4?50?tan(arccos0.85)?31MVar

发电机G5：

P5?2?25?50MW

Q5?50?tan(arccos0.8)?37.5MVar

变压器T1： 其次是变电站的参数计算，我们还是对7个变电站简化为7台变压器来计算。

2?ps?VN73?11023RT1??10??103?3.450?232SN(16?10)

2Vs%?VN10.5?1102

XT1??10??10?79.406?SN16?103

?S01??p0?j

变压器T2：（双并联） I0%?SN?(0.0157?j0.0800)MVA 100

RT2

XT221?ps?VN189?110233???10???10?1.346?2322SN2(20?10)21Vs%?VN110.5?1102???10???10?31.7625? 32SN220?10

?S02?2?(?p0?j

变压器T3：（四并联） I0%?SN)?(0.0372?j0.2000)MVA100

8

信息工程学院课程设计（论文）

21?ps?VN1121?1102

33RT3???10???10?0.092?2324SN4(63?10)

XT321Vs%?VN110.5?1102???10???10?5.042? 34SN463?10?S03?4?(?p0?j

变压器T4：（双并联） I0%?SN)?(0.1760?j0.8820)MVA100

1RT1?1.7250?2

1XT4?XT1?39.7030? 2

?S04?2?S01?(0.0314?j0.1600)MVART4?

变压器T5：

RT5?4RT3?0.3680?

XT5?4XT3?20.168?

?S05?1?S03?(0.0440?j0.2205)MVA4

变压器T6：（两个三绕组变压器并联）

RT6?1?RT6?2?RT6?3

Vs1%?163?352???103?0.386? 322(10?10)1?[Vs(1?2)%?Vs(1?3)%?Vs(2?3)%]?10.752

1Vs2%??[Vs(1?2)%?Vs(2?3)%?Vs(1?3)%]??0.25 2

1Vs3%??[Vs(1?3)%?Vs(2?3)%?Vs(1?2)%]?6.752

21Vs1%?VNXT6?1???10?6.584?2SN

XT6?221Vs2%?VN???10??0.153?2SN

21Vs3%?VNXT6?3???10?4.134? 2SN

I%?S06?2?(?P06?j0?10)?(0.0264?j0.1100)MVA 100

变压器T7：（双并联）

21?ps?VN150?3523RT7???10???103?0.306?2322SN2(10?10)

XT721Vs%?VN110.5?352???10???10?6.431? 32SN210?10?S07?2?(?p0?jI0%?SN)?(0.0220?j0.1100)MVA100

9 再次是传输线参数计算，5条传输线的具体计算如下。

信息工程学院课程设计（论文）

根据教材查得r0?0.21?/km x0?0.4?/km b0?2.8?10?6S/km 线路L1：

RL1?r0?l1?0.21?40?8.4?XL1?x0?l1?0.4?40?16?BL1?b0?l1?2.8?10?6?40?1.12?10?4S

12?QL1??BL1VN??0.6776MVar2

线路L2：

RL2?r0?l2?0.21?130?27.3?XL2?x0?l2?0.4?130?52?BL2?b0?l2?2.8?10?6?130?3.64?10?4S

12?QL2??BL2VN??2.2022MVar2

线路L3：（双回路）

11RL3??r0?l3??0.21?70?7.35?22

11XL3??x0?l3??0.4?70?14?22 BL3?2?b0?l3?2?2.8?10?6?70?3.92?10?4S

1?QL3??BL3VN2??2.3716MVar2

线路L4：

RL4?r0?l4?0.21?60?12.6?XL4?x0?l4?0.4?60?24?BL4?b0?l4?2.8?10?6?60?1.68?10?4S

12?QL4??BL4VN??1.0164MVar2

线路L5：（双回路）

11RL5??r0?l5??0.21?20?2.1?22

11XL5??x0?l5??0.4?20?4?22 BL5?2?b0?l5?2?2.8?10?6?20?1.12?10?4S

1?QL5??BL3VN2??0.0686MVar2

10

信息工程学院课程设计（论文）

3.2手工计算

FLR1：

P2102

?ST1?2(RT1?jXT1)?(3.450?j74.406)?(0.0285?j0.6562)MVAVN1102

Sa?10MW??ST1??S01?j?QL1?(10.0442?j0.1142)MVA

P2?Q210.04422?0.11422?SL1?(RL1?jXL1)?(8.4?j16)?(0.070?j0.1334)MVA22VN110?ST2P2?Q2402?452?(RT2?jXT2)?(1.346?j31.7625)?(0.4032?j9.5156)MVA22VN110

Sb?SG1?20??ST2?60?j45?20?0.4032?j9.5156?(39.5968?j35.4844)MVASc?Sb?Sa?25?jQL1??SL1?(4.4826?j35.9144)MVAFLR2：

?ST3P2?Q22522?1562?(RT3?jXT3)?(0.092?j5.042)?(0.6679?j36.6024)MVAVN21102

P2?Q24.49312?34.10482?(RL2?jXL2)?(27.3?j52)?(2.67?j5.0854)MVA22VN110Sc'?(4.4931?j34.1048)MVA?SL2

Sd?SG2?Sc'?120??ST3??S03?jQL2??SL2?(132.9792?j149.229)MVAFLR3：

?ST4P2?Q262?272?(RT4?jXT4)?(1.725?j39.703)?(0.1091?j2.5101)MVA22VN110

'Sd?(133.5955?j149.9956)MVA

P2?Q2133.59552?149.99562?SL3?(RL3?jXL3)?(7.35?j14)?(24.51?j46.682)MVAVN21102

'Se?SG3?Sd?30?25??ST4??S04?jQL3??SL3?(89.945?j130.0151)MVA

FLR4：

?ST5P2?Q2502?312?(RT5?jXT5)?(0.368?j20.168)?(0.1052?j5.7687)MVA22VN110

P2?Q292.74872?133.99372?(RL4?jXL4)?(12.6?j24)?(27.654?j52.674)MVAVN21102Se'?(92.7481?j133.9937)MVA?SL4

Sf?SG4?Se'?80??ST5??S05?jQL4??SL4?(34.9449?j107.3469)MVAFLR5:

11

信息工程学院课程设计（论文）

152

?ST7?2?(0.306?j6.431)?(0.0562?j1.1812)MVA35

Sh?15??ST7??S07?j?QL5?(15.0782?j0.3422)MVA

15.07822?0.34222

?SL5??(2.1?j4)?(0.3899?j0.743)MVA352

Si?Sh??SL5??S06?j?QL5?5?(20.4945?j1.1266)MVA

152?37.52

?ST6?3??(0.386?j4.34)?(0.514?j5.7793)MVA352

20.65052?0.54512

?ST6?2??(0.386?j0.153)?(0.1345?j0.0533)MVA235

26.3362?98.73692

?ST6?1??(0.386?j6.584)?(3.2905?j56.1256)MVA352

Sg?Sf??ST6?1?SG5??ST6?2??ST6?3?Si?35?(25.5114?j194.12)MVA

计算每一个FLR的功率分布和电压分布计算如下： FLR1：

?VT2?PR?QX40?1.346?45?31.7625??12.8970kVVN115

Vb?115??VT2?102.1030kV

PR?QX10.0442?8.4?0.1442?16?VL1???0.8489kVVb102.1030

Va?Vb??VL1?101.2541kV

FLR2：

功率分布：

SL2?Z?Z*T3*

L2*Sd?T3(Vb?VN)Z?ZL2

****?VN?(0.092?j5.042)?(132.9792?j149.229)?1418.6727.392?j57.042T3?(4.8812?j13.8097)MVAST3?Z?Z*L2*

L2*Sd?T3(Vb?VN)Z?ZL2?VN?(27.3?j52)?(132.9792?j149.229)?1418.6727.392?j57.042T3

?(108.687?j122.62)MVA

电压分布：

Sc1?SL2??SL2?(4.8812?j13.8097)?(2.67?j5.0854)?(7.5512?j8.7243)MVA

7.5512?27.3?8.7243?52??2.424kV102.1030

Vd?Vb??VL2?102.103?(?2.424)?104.527kV?VL2?FLR3：

功率分布：

12

信息工程学院课程设计（论文）

*

SL3?Z?Z*T4*

L3Se?T4(VG3?Vd)Z?ZL3

****?VN?(1.725?j39.703)?(89.945?j130.0151)?1037.9279.075?j53.73T4?(59.444?j16.846)MVAST4?Z?Z*L3*

L3*Se?T4(Vb?VN)Z?ZL3?VN?(7.35?j14)?(89.945?j130.0151)?1037.9279.075?j53.73T4

?(31.811?j60.1256)MVA

电压分布：

Se1?SL3??SL3?(59.444?j19.846)?(24.51?j46.682)?(83.954?j26.836)MVA

83.954?7.35?26.836?14?9.404kV105.5643

Ve?Vd??VL3?96.16kV?VL3?

FLR4：

功率分布：

SL4?Z?Z*T5*

L4*Sf?T5(VG3?Vd)Z?ZL4

****?VN=(0.368?j20.168)?(34.9449?j107.3469)?1037.92712.968?j44.168T5?(20.843?j19.689)MVAST4?Z?Z*L4*

L4*Sf?T5(VG3?Vd)Z?ZL4?VN=(12.6?j24)?(34.9449?j107.3469)?1037.92712.968?j44.168T5

?(1.398?j44.389)MVA

电压分布：

Se1?SL3??SL3?(59.444?j16.846)?(24.51?j46.682)?(83.954?j63.528)MVA

83.954?12.6?63.528?24?24.464kV105.5643

Ve?Vd??VL3?81.10kV?VL4?

FLR5:

这里我们先将f点和发电机G5当做电源，经过ZT61和ZT63构成两端供电网络以g点作为运算负荷进行计算。

ST6?

ST4(0.386?j4.134)?(20.2656?j70.9293)?(22.0938?37)?35?(3.900?j25.1175)MVA0.772?j10.718(0.386?j6.584)?(20.2656?j70.9293)?(22.0938?37)?35??(16.5061?j91.7905)MVA0.772?j10.718电压分布：

13

信息工程学院课程设计（论文）

ST631?ST63??ST63?(16.6421?j97.5698)MVA

16.6421?0.386?97.5698?4.134?10.9186kV37

Vg?37??VT63?26.0814V?VT63?

20.2656?0.386?70.9293?(?0.153)??0.1162kV 26.0814

Vi?Vg??VT62?26.1976?VT62?

20.4945?2.1?1.1266?4?1.815kV26.1976

Vh?Vi??VL5?24.3826?VL5?

14

信息工程学院课程设计（论文）

4设计图及源程序

4.1MATLAB仿真

相关的原始数据输入格式如下：

1、B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写。

2、对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点编号，将变压器的串联阻抗置于低压侧处理，第三列为支路的串列阻抗参数，第四列为支路的对地导纳参数，第五烈为含变压器支路的变压器的变比，第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，“0”为不含有变压器。

3、B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。

4、X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地参数。

4.2潮流计算源程序

close all

clear all

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:n1

if 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));

end

15

信息工程学院课程设计（论文）

end

Y

OrgS=zeros(2*n-2,1);

DetaS=zeros(2*n-2,1); %将OrgS、DetaS初始化

%创建OrgS，用于存储初始功率参数

h=0;

j=0;

for i=1:n %对PQ节点的处理

if i~=isb&B2(i,6)==2

h=h+1;

for j=1:n

OrgS(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))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));

end

end

end

for i=1:n %对PV节点的处理，注意这时不可再将h初始化为0

if i~=isb&B2(i,6)==3

h=h+1;

for j=1:n

OrgS(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))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));

end

end

end

OrgS

%创建PVU 用于存储PV节点的初始电压

PVU=zeros(n-h-1,1);

t=0;

for i=1:n

if B2(i,6)==3

t=t+1;

PVU(t,1)=B2(i,3);

end

end

PVU

16

信息工程学院课程设计（论文）

%创建DetaS，用于存储有功功率、无功功率和电压幅值的不平衡量

h=0;

for i=1:n %对PQ节点的处理

if i~=isb&B2(i,6)==2

h=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);

end

end

t=0;

for i=1:n %对PV节点的处理，注意这时不可再将h初始化为0

if i~=isb&B2(i,6)==3

h=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; end

end

DetaS

%创建I，用于存储节点电流参数

i=zeros(n-1,1);

h=0;

for i=1:n

if i~=isb

h=h+1;

I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));

end

end

I

%创建Jacbi(雅可比矩阵)

Jacbi=zeros(2*n-2);

h=0;

k=0;

for i=1:n %对PQ节点的处理

if B2(i,6)==2

h=h+1;

for j=1:n

if j~=isb

k=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));

17

信息工程学院课程设计（论文）

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);

end

if k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行 k=0;

end

end

end

end

end

k=0;

for i=1:n %对PV节点的处理

if B2(i,6)==3

h=h+1;

for j=1:n

if j~=isb

k=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;

end

if k==(n-1) %将用于内循环的指针置于初始值，以确保雅可比矩阵换行

k=0;

end

end

end

end

end

Jacbi

%求解修正方程，获取节点电压的不平衡量

18

信息工程学院课程设计（论文）

DetaU=zeros(2*n-2,1);

DetaU=inv(Jacbi)*DetaS;

DetaU

%修正节点电压

j=0;

for i=1:n %对PQ节点处理

if B2(i,6)==2

j=j+1;

B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);

end

end

for i=1:n %对PV节点的处理

if B2(i,6)==3

j=j+1;

B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);

end

end

B2

%开始循环********************************************************************** while abs(max(DetaU))>pr

OrgS=zeros(2*n-2,1); %!!!初始功率参数在迭代过程中是不累加的，所以在这里必须将其初始化为零矩阵

h=0;

j=0;

for i=1:n

if i~=isb&B2(i,6)==2

h=h+1;

for j=1:n

OrgS(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))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));

end

end

end

for i=1:n

if i~=isb&B2(i,6)==3

h=h+1;

for j=1:n

OrgS(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))* 19

信息工程学院课程设计（论文）

imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));

end

end

end

OrgS

%创建DetaS

h=0;

for i=1:n

if i~=isb&B2(i,6)==2

h=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);

end

end

t=0;

for i=1:n

if i~=isb&B2(i,6)==3

h=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; end

end

DetaS

%创建I

i=zeros(n-1,1);

h=0;

for i=1:n

if i~=isb

h=h+1;

I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));

end

end

I

%创建Jacbi

Jacbi=zeros(2*n-2);

h=0;

k=0;

for i=1:n

if B2(i,6)==2

h=h+1;

for j=1:n

if j~=isb

k=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));

20

信息工程学院课程设计（论文）

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);

end

if k==(n-1)

k=0;

end

end

end

end

end

k=0;

for i=1:n

if B2(i,6)==3

h=h+1;

for j=1:n

if j~=isb

k=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;

end

if k==(n-1)

k=0;

end

end

end

end

end

Jacbi

21

信息工程学院课程设计（论文）

DetaU=zeros(2*n-2,1);

DetaU=inv(Jacbi)*DetaS;

DetaU

%修正节点电压

j=0;

for i=1:n

if B2(i,6)==2

j=j+1;

B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1); end

end

for i=1:n

if B2(i,6)==3

j=j+1;

B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1); end

end

B2

Times=Times+1; %迭代次数加1

end

Times

22

信息工程学院课程设计（论文）

5 总结

通过本次课程设计让我有复习了一次潮流计算的相关知识，跟家清晰了什么事潮流计算以及潮流计算的在电力系统的重要性。电力系统的稳定运行状况即是正常运行状况，是指电力系统在稳定运行条件下电压、功率的分布，也称为潮流分布。电力系统分析的潮流计算是电力系统分析的一个重要的部分。通过对电力系统潮流分布的分析和计算，可进一步对系统运行的安全性，经济性进行分析、评估，提出改进措施。同时潮流分布也是电力系统规划设计的一项基础工作。

整个计算过程的模型建立并不是十分复杂，但计算过程十分繁琐、计算量相当的大，而且由于枝节太多很容易算错。不过在计算潮流计算的过程中却对以往学过的电力系统分析的相关知识进行了一次较为深入的复习。而且整个计算对计算量的要求很大，锻炼了我们的计算能力。而且对细节的把握也得到了锻炼，做题的精细程度得到了提高。

最后，我想衷心的感谢我们的指导老师耿东山老师。对于手工计算的每个问题，老师都为我们进行最详细的指导，对于程序的设计，老师也全力帮助我们进行程序语言的学习，谢谢老师。

23

信息工程学院课程设计（论文）

参考文献

[1]何仰赞, 温增银《电力系统分析》（第三版）［M]. 华中科技大学，2002

[2]/view/627420.

[3]王守相，刘玉田 电力系统潮流计算研究现状--《山东电力技术》19xx年05期

[4]何仰赞，温增银. 电力系统分析（上册）第三版［M]. 湖北：华中科技大学出版社，2002

[5] 刘同娟. MATLAB在电路分析中的应用.电气电子教学学报.2002

[6] 西安交通大学等.电力系统计算[M].北京：水利电力出版社，1993.12

[7] 李光琦.电力系统暂态分析[M].北京： 水利电力出版社，2002.5

[8]何仰赞，温增银. 电力系统分析（下册）第三版［M]. 湖北：华中科技大学出版社，2002

[9]韦化，李滨，杭乃善，等．大规模水一火电力系统最优潮流的现代内点算法实现[J]．中

国电机工程学报，2003．23(6)：13一l8．

24