电力系统剖析潮流计算例题[整理版]
电力系统课程设计潮流计算
电力系统课程设计潮流计算潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。
对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求；对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(线路、变压器等)是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。
潮流计算是电力系统分析最基本的计算。
除它自身的重要作用之外，潮流计算还是网损计算、静态安全分析、暂态稳定计算、小干扰静态稳定计算、短路计算、静态和动态等值计算的基础。
实际电力系统的潮流计算主要采用牛顿-拉夫逊法。
按电压的不同表示方法，牛顿-拉夫逊潮流计算分为直角坐标形式和极坐标形式两种。
本次计算采用直角坐标形式下的牛顿-拉夫逊法，牛顿-拉夫逊法有很好的收敛性，但要求有合适的初值。
传统的潮流计算程序缺乏图形用户界面，结果显示不直接难与其他分析功能集成。
网络原始数据输入工作大量且易于出错。
本文采用MATLAB 语言运行WINDOWS操作系统的潮流计算软件。
目前MATLAB已成为国际控制界最流行、使用最广泛的语言了。
它的强大的矩阵处理功能给电力系统的分析、计算带来很多方便，而且采用MATLAB界面直观，运行稳定，计算准确。
所以本次课程设计程序设计采用MATLAB计算。
1.1.2设计要求1.程序源代码；2.给定题目的输入，输出文件；3.程序说明；4.给定系统的程序计算过程；5.给定系统的手算过程（至少迭代2次）。
1.2设计题目电力系统潮流计算（牛顿-拉夫逊法、P-Q分解法）1.3设计内容1.根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；2.赋予各节点电压变量（直角坐标系形式）初值后，求解不平衡量；3.形成雅可比矩阵；4.求解修正量后，重新修改初值，从2开始重新循环计算；5.求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和平衡节点功率；6.上机编程调试；7.计算分析给定系统潮流分析并与手工计算结果做比较分析；8.书写课程设计说明书。
输配电系统潮流计算与分析考核试卷
3.针对复杂网络，如何选择合适的潮流分布计算方法？请结合实际案例分析。
4.在潮流分布计算中，如何分析电网的稳定性？请结合实际案例进行说明。
六、案例题（本题共2小题，每题5分，共10分）
1.案例题：某地区电网由10kV和35kV两个电压等级组成，包含10个节点和15条线路。已知各节点负荷和发电机出力，以及线路的阻抗参数。请利用潮流分布计算方法，计算电网在正常运行状态下的潮流分布，并绘制潮流分布图。
A.确定电网运行状态
B.分析电网运行稳定性
C.优化电网运行方式
D.评估电网设备性能
6.以下哪些是潮流分布计算中的节点类型？（）
A. PQ节点
B. PV节点
C.负荷节点
D.发电机节点
7.潮流分布计算中，以下哪些是影响潮流分布的因素？（）
A.线路阻抗
B.节点负荷
C.发电机出力
D.网络拓扑结构
8.以下哪些是潮流分布计算中的收敛判据？（）
A.是固定的
B.是变化的
C.是对称的
D.是非对称的
10.下列哪种方法在潮流计算中用于求解节点电压（）。
A.牛顿-拉夫逊法
B.迭代法
C.直接法
D.数值积分法
11.潮流分布图中，负荷节点的电压幅值（）。
A.是固定的
B.是变化的
C.是对称的
D.是非对称的
12.在潮流计算中，若网络中存在开路，则该开路为（）。
A.迭代误差
B.计算时间
C.迭代次数
D.网络潮流分布图
9.潮流分布计算中，以下哪些是牛顿-拉夫逊法的优点？（）
A.计算速度快
B.精度高
电力系统潮流计算与分析考核试卷
2.潮流计算的基本方程包括功率方程、电压方程和______。
（）
3.牛顿-拉夫逊法在潮流计算中属于______求解法。
（）
4.电力系统的基准功率通常选择为系统最大______的______倍。
（）（）
5.在电力系统中，PV节点是指有功功率______，电压幅值已知的节点。
（）
6.电力系统的潮流计算主要目的是确定系统各节点的______和______。
D.有功功率和无功功率可以通过变压器相互转换
18.以下哪些设备在电力系统中通常用于控制无功功率？（）
A.电容器
B.电感器
C.静止无功发生器（SVG）
D.晶闸管控制串联电容器（TSC）
19.在电力系统潮流计算中，以下哪些情况可能导致计算不精确？（）
A.系统包含大量的分布式发电
B.系统包含多个不同电压等级的网络
C.系统运行在重载状态
D.系统中的线路参数变化较大
20.以下哪些是电力系统潮流计算的重要应用？（）
A.系统优化
B.故障分析
C.安全评估
D.经济调度
（以下继续其他题型内容）
三、填空题（本题共10小题，每小题2分，共20分，请将正确答案填到题目空白处）
1.在电力系统潮流计算中，平衡节点的电压幅值通常设定为______，相角为______。
A. PQ节点、PV节点、平衡节点
B. PQ节点、PV节点、阻抗节点
C. PQ节点、平衡节点、电流节点
D. PV节点、平衡节点、电压节点
7.以下哪一种方法属于潮流计算的直接法？（）
A.牛顿-拉夫逊法
B.高斯-赛德尔法
C.快速分解法
D.矩阵分解法
8.在潮流计算中，如果某一节点的有功功率和无功功率均已知，则该节点属于以下哪种类型？（）
电力系统基础第三章习题
第三章 电力系统潮流计算一、填空题1. 输电线路始末两端电压的向量差称为 。
2. 输电线路始末两端电压的数值差称为 。
3. 输电线路始端电压或末端电压与线路额定电压的差值称为 。
4. 电力系统的潮流计算是求解电力网 、 以及 。
5. 所谓 是指输电线首端与末端电压的(相量)之差。
是指输电线某点的实际电压与额定电压的(数值)的差。
6. 两端供电网络的两端电压相等时其 。
二、判断题 1. 电压降落的表达式为 （ ）2. 电压损耗的表达式为 （ ）3. 所谓线损率系指线路上损耗的电能与线路始端输入电能的比值。
（ ）4. 在环型网络的潮流中，其有功分点与无功分点总是重合的（ ）。
5. 线路首端的电压一定高于末端的电压。
( )6. 高压电网中无功功率分点的电压最低。
（ ）7. 任何多电压等级环网中都存在循环功率。
（ ）三、选择题1． 输电线路单位长度的电阻主要决定于（ D ）。
A 材料与对地高度B 电晕损耗与电压C 几何均距与半径D 材料与截面大小2． 线路电压在220KV 及以上时，采用分裂导线的目的是（ ）。
A 减少导线重量B 增大导线应力C 增大线路电容D 减少电晕与线路电抗3． 架空输电线路全换位的目的是（ ）。
A 使三相线路的电阻参数相等；B 使三相线路的电抗与电纳参数相等；C 减小线路电抗；D 减小线路电阻。
4． 架空输电线路的电抗与导线之间几何平均距离的关系为（ A ）。
A 几何平均距离越大，电抗越大；B 几何平均距离越大，电抗越小；C 输电线路的电抗与几何平均距离无关；D 改变导线之间的几何平均距离可以明显改变线路的电抗。
5． 架空输电线路的电纳与导线之间几何平均距离的关系为（ ）。
A 几何平均距离越大，电纳越大；B 几何平均距离越大，电纳越小；U j U U U δ+∆=-.2.1U U U∆≈-21C 输电线路的电纳与几何平均距离无关；D 改变导线之间的几何平均距离可以明显改变线路的电纳。
潮流计算 例题
例题：如图1所示的简单电力网中，已知变压器的参数为S N =31.5MV A ，0S S 031kW,190kW,%=10.5,%=0.7P P U I ∆=∆=；线路单位长度的参数为61110.21/km,0.416/km, 2.7410S/km r x b -=Ω=Ω=⨯。
如图所示的简单电力网中，当线路首端电压U A =120kV 时，试求：（1）线路和变压器的电压损耗；（2）变压器运行在额定变比时的低压侧电压及电压偏移。
说明：以上计算忽略电压降落的横分量。
图1解：如题画等值电路图如下：线路参数为：0.21408.40.4164016.64l l l l R rl X x l ==⨯=Ω==⨯=Ω变压器参数为Ω=⨯⨯⨯=⨯∆=317.210)105.31(110190103232322N N S T S U P R Ω=⨯⨯⨯⨯⨯=33.4010105.311001105.1010100%33232N N S T S U U X （1） 变压器的功率损耗和励磁功率为222T S 02N 2515()19031193.760.1937631.5S P P P kW MW S +∆=∆+∆=⨯+==222S 0N T N %%10.5(2515)0.731.5 3.0538var 10010010031.5100U S I S Q M S ⨯+⨯∆=+=+=⨯1点处线路的充电功率var 66308.01104074.22121222M lU b Q N l B =⨯⨯⨯==计算L S 2 为：MVAj j Q Q Q j P P S B T LD T LD L 39.1719.25)66308.00538.315(19376.025)(22+=-+++=-∆++∆+=线路阻抗中的功率损耗为：MW R U Q P P l L L l 65044.0104021.011017390251901032223222222=⨯⨯⨯+=⨯+=∆-- v a r 2885.11040416.011017390251901032223222222M X U Q P Q l L L l =⨯⨯⨯+=⨯+=∆-- 计算功率1S 为 M V Aj j Q Q j P P S l L l I 68.1884.25)2885.139.17(65044.019.25)(221+=+++=∆++∆+=线路电压损耗（忽略电压降落的横分量） 1125.848.418.6816.64 4.40120=l l l A PR Q X U kV U +⨯+⨯∆== 1点电压为：1120 4.40115.60-A l U U U kV =∆=-=计算功率TS 2 为 M V Aj j Q Q j P P S T LD T LD T 833.17163.25)833.215(163.025)(2+=+++='∆++'∆+= 变压器电压损耗 22125.163 2.3217.83340.33 6.73115.60=T T T T T P R Q X U kV U +⨯+⨯∆== （2） 变压器低压侧折算到高压侧的电压为21115.60 6.73108.87=-T U U U kV '∆=-= 变压器低压侧的实际电压 22108.8710.8910=U U kV k '== 电压偏移为2210.8910%100%8.9%10N N U U m U --=⨯==。
潮流计算实例计算
潮流计算实例计算潮流例题：根据给定的参数或⼯程具体要求（如图），收集和查阅资料；学习相关软件（软件⾃选：本设计选择Matlab进⾏设计）。
2.在给定的电⼒⽹络上画出等值电路图。
3.运⽤计算机进⾏潮流计算。
4.编写设计说明书。
⼀、设计原理1．⽜顿-拉夫逊原理⽜顿迭代法是取x0 之后，在这个基础上，找到⽐x0 更接近的⽅程的跟，⼀步⼀步迭代，从⽽找到更接近⽅程根的近似跟。
⽜顿迭代法是求⽅程根的重要⽅法之⼀，其最⼤优点是在⽅程f(x) = 0 的单根附近具有平⽅收敛，⽽且该法还可以⽤来求⽅程的重根、复根。
电⼒系统潮流计算，⼀般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据⽹络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率⽅程，由于功率⽅程⾥功率是已知的，电压的幅值和相⾓是未知的，这样潮流计算的问题就转化为求解⾮线性⽅程组的问题了。
为了便于⽤迭代法解⽅程组，需要将上述功率⽅程改写成功率平衡⽅程，并对功率平衡⽅程求偏导，得出对应的雅可⽐矩阵，给未知节点赋电压初值，⼀般为额定电压，将初值带⼊功率平衡⽅程，得到功率不平衡量，这样由功率不平衡量、雅可⽐矩阵、节点电压不平衡量（未知的）构成了误差⽅程，解误差⽅程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带⼊原来的功率平衡⽅程，并重新形成雅可⽐矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，⼀般迭代三到五次就能收敛。
⽜顿—拉夫逊迭代法的⼀般步骤：（1）形成各节点导纳矩阵Y。
（2）设个节点电压的初始值U和相⾓初始值e 还有迭代次数初值为0。
（3）计算各个节点的功率不平衡量。
（4）根据收敛条件判断是否满⾜，若不满⾜则向下进⾏。
（5）计算雅可⽐矩阵中的各元素。
（6）修正⽅程式个节点电压（7）利⽤新值⾃第（3）步开始进⼊下⼀次迭代，直⾄达到精度退出循环。
（8）计算平衡节点输出功率和各线路功率2．⽹络节点的优化１）静态地按最少出线⽀路数编号这种⽅法由称为静态优化法。
电力系统课程设计潮流计算
摘要潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。
对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求；对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(线路、变压器等)是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。
潮流计算是电力系统分析最基本的计算。
除它自身的重要作用之外，潮流计算还是网损计算、静态安全分析、暂态稳定计算、小干扰静态稳定计算、短路计算、静态和动态等值计算的基础。
实际电力系统的潮流计算主要采用牛顿-拉夫逊法。
按电压的不同表示方法，牛顿-拉夫逊潮流计算分为直角坐标形式和极坐标形式两种。
本次计算采用直角坐标形式下的牛顿-拉夫逊法，牛顿-拉夫逊法有很好的收敛性，但要求有合适的初值。
传统的潮流计算程序缺乏图形用户界面，结果显示不直接难与其他分析功能集成。
网络原始数据输入工作大量且易于出错。
本文采用MATLAB语言运行WINDOWS操作系统的潮流计算软件。
目前MATLAB已成为国际控制界最流行、使用最广泛的语言了。
它的强大的矩阵处理功能给电力系统的分析、计算带来很多方便，而且采用MATLAB界面直观，运行稳定，计算准确。
所以本次课程设计程序设计采用MATLAB计算。
关键词：电力系统潮流计算牛顿—拉夫逊法潮流计算 MATLAB目录一、概述1.1设计目的与要求 (3)1.1.1 设计目的 (3)1.1.2 设计要求 (3)1.2 设计题目 (3)1.3 设计内容 (3)二电力系统潮流计算概述 (4)2.1 电力系统简介 (4)2.2 潮流计算简介 (4)2.3 潮流计算的意义及其发展..................... . (5)三潮流计算设计题目 (6)3.1 潮流计算题目 (6)3.2 对课题的分析及求解思路 (7)四潮流计算算法及手工计算 (7)4.1 极坐标下P-Q法的算法 (7)4.2 节点电压方程 (8)4.3节点导纳矩阵 (9)4.4 导纳矩阵在潮流计算 (10)4.5 潮流计算的手工计算 (12)五 Matlab概述 (13)5.1 Matlab简介 (14)5.2 Matlab的应用 (14)5.3 矩阵的运算 (14)5.3.1 与常数的运算 (14)5.3.2 基本数学运算 (14)5.3.3 逻辑关系运算 (14)5.4 Matlab中的一些命令 (15)六潮流计算流程图及源程序 (18)6.1 潮流计算流程图 (18)6.2 潮流计算源程序图 (19)6.3 运行计算结果 (27)七总结 (29)八参考文献 (29)第一章系统概述1.1 设计目的与要求1.1.1设计目的1.掌握电力系统潮流计算的基本原理；2.掌握并能熟练运用一门计算机语言（MATLAB语言或C语言或C++语言）；3.采用计算机语言对潮流计算进行计算机编程。
5、开式电力系统的潮流计算_修改
将理想变压器后面的参数折算到一次侧，再前推回代； 注意，变压器后面的电压是折算后的电压。 麻烦！
5.3 多电压等级开式电网潮流计算
Q：系统中的变压器应该如何处理？
5.1 潮流计算的概念 5.2 开式电力系统的潮流计算 5.3多电压等级开式电网潮流计算 5.4 开式潮流计算例题
按照图c进行前推回代； 经过理想变压器时，功率不变，电压比等于变比； 手算方便。
5.1 潮流计算的概念 5.2 开式电力系统的潮流计算 5.3多电压等级开式电网潮流计算 5.4 开式潮流计算例题
5.1 潮流计算的概念
潮流计算与电路计算的区别？
已知条件的变化：由电流到复功率。 模型物理基础的变化：由欧姆定律到功率守恒定律。 计算模型的变化：由线性方程组到非线性方程组。
潮流计算的两个基础——电压降落和功率损耗
U D
∆S0
U
S3:从电源点开始，顺功率传输方向，计算各段线路的
电压降落和各节点电压（不算功率损耗），完成回代；
S4:用新电压反复迭代，直到满足精度要求，收敛。
两步计算（近似）
Q：系统中的变压器应该如何处理？
5.3 多电压等级开式电网潮流计算
Q：系统中的变压器应该如何处理？
5.1 潮流计算的概念 5.2 开式电力系统的潮流计算 5.3多电压等级开式电网潮流计算 5.4 开式潮流计算例题
5.3多电压等级开式电网潮流计算
D
5.4 开式潮流计算例题
SD
U D
已知： SD、UD 已知： SA、UA 但实际中，如果已知不同点的电压和功率，怎么办？
5.2 开式电力系统的潮流计算
一般都是已知不同点的电压和功率：
5.1 潮流计算的概念 5.2 开式电力系统的潮流计算
(完整版)电力系统分析试题及参考答案
作业题一参考答案一、填空题1．降压变压器高压侧的主分接头电压为220kv ，若选择＋2×2.5%的分接头，则该分接头电压为 231KV 。
2．电力系统中性点有效接地方式指的是 中性点直接接地 。
3．输电线路的电气参数包括电抗、电导、电纳和 电阻 。
4．输电线路的电压偏移是指线路始端或末端母线的实际运行电压与线路 额定电压 的数值差。
5．电力系统的潮流分布一般是用各节点的电压和 功率 表示。
6．调整发电机组输出的有功功率用来调整电力系统运行的 频率 。
7．复合故障一般是指某一时刻在电力系统 二个及以上地方 发生故障。
8．用对称分量法计算不对称故障，当三相阻抗完全对称时，则其序阻抗矩阵Zsc 的非对角元素为 零 。
9．系统中发生单相接地短路时故障点短路电流的大小是零序电流的 3 倍。
10．减小输出电元件的电抗将 提高（改善） 系统的静态稳定性。
二、单项选择题在每小题列出的四个选项中只有一个选项是符合题目要求的，请将正确选项前的字母填在题后的括号内。
11．同步发电机的转速和系统频率之间是否有严格的关系（ ② ）①否 ②是 ③不一定 ④根据发电机的形式定12．三绕组变压器的结构、通常将高压绕组放在（ ③ ）①内层 ②中间层 ③外层 ④独立设置13．中性点以消弧线圈接地的电力系统，通常采用的补偿方式是（ ③ ）①全补偿 ②欠补偿 ③过补偿 ④有时全补偿，有时欠补偿14．三相导线的几何均距越大，则导线的电抗（ ② ）①越大 ②越小 ③不变 ④无法确定15．变压器的电导参数G T ，主要决定于哪一个实验数据（ ① ）①△P O ②△P K ③U K % ④I O %16.当功率的有名值为s ＝P ＋jQ 时（功率因数角为ϕ）取基准功率为S n ，则有功功率的标么值为（ ③ ） ①ϕcos S P n ⋅ ②ϕsin S P n ⋅ ③n S P ④nS cos P ϕ⋅17.环网中功率的自然分布是（ ④ ）①与电阻成正比分布 ②与电抗成正比分布③与阻抗成正比分布 ④与阻抗成反比分布18．电力系统中PQ 节点的数量（ ② ）①全都是 ②大量的 ③少量的 ④必有且一般只设一个19．潮流计算中，要求某些节点之间电压的相位差应满足的约束条件是（④ ）①|-j i δδ｜>|-j i δδ｜min ②|-j i δδ｜<|-j i δδ｜min③|-j i δδ｜>|-j i δδ｜max ④|-j i δδ｜<|-j i δδ｜max20．在同一时间内，电力网的电能损耗与供电量之比的百分值称为（ ② ）①负载率 ②网损率 ③供电率 ④厂用电率21．电力系统的频率主要决定于（ ① ）①有功功率的平衡 ②无功功率的平衡③电压质量 ④电流的大小22．关于顺调压电压调整方式的描述，错误的是（ ② ）①高峰负荷时允许中枢点电压略低②低谷负荷时允许中枢点电压略低③适用于用户对电压要求不高的场合④适用于供电线路不长的场合23．通过改变变压器变比，实质上（ ③ ）①改变了电压损耗的数值 ②改变了负荷变化时次级电压的变化幅度③改变了电力网的无功功率分布 ④增加了整个电力系统的无功功率容量24．三相短路时，非周期分量极小值，只能是在某种情况（① ）①一相中出现 ②同时在两相中出现③三相均不出现 ④只有故障相出现其它相不出现25．同步机的各种电抗间关系为(③ )①'x "x x x d d q d >>> ②"x 'x x x d d d q >>>③ "x 'x x x d d q d >>> ④"x x 'x x d q d d >>>26.若ac 两相发生短路，则基准相选择为( ② )①a 相 ②b 相 ③c 相 ④a 相和c27．下网K 点发生两相短路接地，其复合序网为图所示( ③ )(其中，1,2，0分别为正序、负序、零序阻抗)28．越靠近电源点负序电压越( ① )①低 ②高 ③不变 ④无法确定29．作为判据0d dP E >δ主要应用于分析简单系统的（ ③ ） ①暂态稳定 ②故障计算 ③静态稳定 ④调压计算30．分析简单系统的暂态稳定性可以应用（ ② ）①等耗量微增率准则 ②等面积定则 ③小干扰法 ④对称分量法三、简答题31．电力变压器的主要作用是什么？答：电力变压器的主要作用是升高或降低电压，另外还起到将不同电压等级电网相联系的作用。
《电力系统分析》习题第3-6章(1)
3 简单电力系统潮流计算3．1 思考题、习题1）电力线路阻抗中的功率损耗表达式是什么？电力线路始、末端的电容功率表达式是什么？2）电力线路阻抗中电压降落的纵分量和横分量的表达式是什么？3）什么叫电压降落、电压损耗、电压偏移、电压调整及输电效率？5）对简单开式网络、变电所较多的开式网络和环形网络潮流计算的内容及步骤是什么？6）变压器在额定状况下，其功率损耗的简单表达式是什么？9）为什么要对电力网络的潮流进行调整控制？调整控制潮流的手段主要有哪些？10）欲改变电力网络的有功功率和无功功率分布，分别需要调整网络的什么参数？16）110kV 双回架空线路，长度为150kM ，导线型号为LGJ-120，导线计算外径为15.2mm ，三相导线几何平均距离为5m 。
已知电力线路末端负荷为30+j15MVA ，末端电压为106kV ，求始端电压、功率，并作出电压向量图。
17）220kV 单回架空线路，长度为200kM ，导线型号为LGJ-300，导线计算外径为24.2mm ，三相导线几何平均距离为7.5m 。
已知电力线路始端输入功率为120+j50MVA ，始端电压为240kV ，求末端电压、功率，并作出电压向量图。
18）110kV 单回架空线路，长度为80kM ，导线型号为LGJ-95，导线计算外径为13.7mm ，三相导线几何平均距离为5m 。
已知电力线路末端负荷为15+j10MVA ，始端电压为116kV ，求末端电压和始端功率。
19）220kV 单回架空线路，长度为220kM ，电力线路每公里的参数分别为：kM S b kM x kM r /1066.2,/42.0,/108.06111-⨯=Ω=Ω=、线路空载运行，当线路末端电压为205kV ，求线路始端的电压。
20）有一台三绕组变压器，其归算至高压侧的等值电路如图3-1所示，其中,68~,45~,8.3747.2,5.147.2,6547.232321MVA j S MVA j S j Z j Z j Z T T T +=+=Ω+=Ω-=Ω+=当变压器变比为110/38.5（1+5%）/6.6kV ，U 3=6kV 时，试计算高压、中压侧的实际电压。
电力系统分析 第三章例题
第三章 电力系统潮流分布计算3-2 已知图3-2所示输电线路始末端电压分别为248kV 、220kV ，末端有功功率负荷为220MW ，无功功率负荷为165 MV AR 。
试求始端功率因数。
3-2 解：62.26105.5220422=⨯⨯=∆-y Q (MV AR) 83.33105.5248421=⨯⨯=∆-y Q (MV AR)38.13822062.261652202j j j S +=-+='∙(MV A)求Z 12中的功率损耗：21.194165.23183.55165.1138.13822083.55165.11)408(22038.138220122212j j j S j j S +=+++='+=++=∆∙∙38.160165.23183.3321.194165.2311j j j S +=-+=∙(MV A)8216.038.160165.231165.231cos 22=+=ϕ3-8 额定电压110 kV 的辐射形电网各段阻抗及负荷如图3-8所示。
已知电源A 的电压为121 kV ，求功率分布和各母线电压。
（注：考虑功率损耗，可以不计电压降落的横分量U δ）。
3-8 解：设︒∠=︒∠=∙01100N C U U220kVA习题解图3-8P 2=220MW2=165MV AR8+j40Ω 习题图 3-2 A习题图 3-8083.27545.32676.4338.2407.22207.30676.4338.2)4020(110407.22207.30407.22207.30953.7793.93040593.7793.9407.0271.0)810(407.0271.0)3020(110810222222j j j S S S j j S j j j S S S j j j S S S j j S AB B A ABB B B BC C BBC+=+++=∆+'-=+=++=∆+=--+=''+='--=++--=∆+-=''+=++=∆∙∙∙∙∙∙∙∙∙∙∙已知U A =121kV332.1412140083.2720545.32=⨯+⨯=∆AB U kV668.106332.14121=-=∆-=AB A B U U U kV972.3668.10630)593.7(20)793.9(-=⨯-+⨯-=∆BC U kV64.110972.3668.106=-=∆-=BC B C U U U kV3-13 由A 、B 两端供电的电力网，其线路阻抗和负荷功率等如图3-13示。
电力系统分析潮流计算最终完整版
电力系统分析潮流计算实验报告姓名：XXXXXX 学号：XXXXXXXXXX 班级：XXXXXXXX一、实验目的掌握潮流计算计算机算法的方法，熟悉MATLAB的程序调试方法。
二、实验准备根据课程内容，熟悉MATLAB软件的使用方法，自行学习MATLAB程序的基础语法，并根据所学知识编写潮流计算牛顿拉夫逊法(或PQ分解法) 的计算程序，用相应的算例在MATLAB上进行计算、调试和验证。
三、实验要求每人一组，在实验课时内，调试和修改运行程序，用算例计算输出潮流结果。
四、程序流程五、实验程序%本程序的功能是用牛拉法进行潮流计算%原理介绍详见鞠平著《电气工程》%默认数据为鞠平著《电气工程》例8.4所示数据%B1是支路参数矩阵%第一列和第二列是节点编号。
节点编号由小到大编写%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点编号%第三列为支路的串列阻抗参数，含变压器支路此值为变压器短路电抗%第四列为支路的对地导纳参数，含变压器支路此值不代入计算%第五烈为含变压器支路的变压器的变比，变压器非标准电压比%第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，“0”为不含有变压器%B2为节点参数矩阵%第一列为节点注入发电功率参数%第二列为节点负荷功率参数%第三列为节点电压参数%第四列%第五列%第六列为节点类型参数，“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数%X为节点号和对地参数矩阵%第一列为节点编号%第二列为节点对地参数%默认算例% n=4;% n1=4;% isb=4;% pr=0.00001;% B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0];% B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1];% X=[1 0;2 0.05i;3 0;4 0];clear;clc;num=input('是否采用默认数据？(1-默认数据;2-手动输入)');if num==1n=4;n1=4;isb=4;pr=0.00001;B1=[1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0];B2=[0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1];X=[1 0;2 0.05i;3 0;4 0];elsen=input('请输入节点数:n=');n1=input('请输入支路数:n1=');isb=input('请输入平衡节点号:isb=');pr=input('请输入误差精度:pr=');B1=input('请输入支路参数:B1=');B2=input('请输入节点参数:B2=');X=input('节点号和对地参数:X=');endTimes=1; %迭代次数%创建节点导纳矩阵Y=zeros(n);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)-B1(i,5)/B1(i,3);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+B1(i,5)/B1(i,3)+(1-B1(i,5))/B1(i,3);Y(q,q)=Y(q,q)+B1(i,5)/B1(i,3)+(B1(i,5)*(B1(i,5)-1))/B1(i,3);endendfor i=1:n1Y(i,i)=Y(i,i)+X(i,2); %计及补偿电容电纳enddisp('导纳矩阵为：');disp(Y); %显示导纳矩阵%初始化OrgS、DetaSOrgS=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)==2 %不是平衡点&是PQ点h=h+1;for j=1:n%公式8-74%Pi=ei*(Gij*ej-Bij*fj)+fi*(Gij*fj+Bij*ej)%Qi=fi*(Gij*ej-Bij*fj)-ei*(Gij*fj+Bij*ej)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)));endendendfor i=1:n %对PV节点的处理，注意这时不可再将h初始化为0if i~=isb&B2(i,6)==3 %不是平衡点&是PV点h=h+1;for j=1:n%公式8-75-a%Pi=ei*(Gij*ej-Bij*fj)+fi*(Gij*fj+Bij*ej)%Qi=fi*(Gij*ej-Bij*fj)-ei*(Gij*fj+Bij*ej)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)));endendend%创建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);endend%创建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); %delPiDetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1); %delQiendendt=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,1))-OrgS(2*h-1,1); %delPiDetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2; %delUi endend% DetaS%创建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));%conj求共轭endend%创建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;endendendendenddisp('初始雅可比矩阵为：');disp(Jacbi);%求解修正方程，获取节点电压的不平衡量DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS; %inv矩阵求逆% 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);endend% B2%开始循环**********************************************************************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))*imag(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))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j ))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));endendend% OrgS%创建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,1)=real(B2(i,1))-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;endend% DetaS%创建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));endend% I%创建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;endendendend% JacbiDetaU=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);endend% B2Times=Times+1; %迭代次数加1enddisp('迭代次数为：');disp(Times);disp('收敛时电压修正量为：：');disp(DetaU);for k=1:nE(k)=B2(k,3);e(k)=real(E(k));f(k)=imag(E(k));V(k)=sqrt(e(k)^2+f(k)^2);sida(k)=atan(f(k)./e(k))*180./pi;end%=============== 计算各输出量=========================== disp('各节点的实际电压标幺值E为(节点号从小到大排列)：'); disp(E); %显示各节点的实际电压标幺值E用复数表示disp('-----------------------------------------------------')disp('各节点的电压大小V为(节点号从小到大排列)：');disp(V); %显示各节点的电压大小V的模值disp('-----------------------------------------------------');disp('各节点的电压相角sida为(节点号从小到大排列)：');disp(sida); %显示各节点的电压相for p=1:nfor q=1:nC(p)=C(p)+conj(Y(p,q))*conj(E(q)); %计算各节点的注入电流的共轭值endS(p)=E(p)*C(p); %计算各节点的功率S = 电压X 注入电流的共轭值enddisp('各节点的功率S为(节点号从小到大排列)：');disp(S); %显示各节点的注入功率Sline=zeros(n1,5);disp('-----------------------------------------------------');disp('各条支路的首端功率Si为(顺序同您输入B1时一致)：');for i=1:n1p=B1(i,1);q=B1(i,2);Sline(i,1)=B1(i,1);Sline(i,2)=B1(i,2);if B1(i,6)==0Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))));Siz(i)=Si(p,q);elseSi(p,q)=E(p)*(conj(E(p))*((1-B1(i,5))/B1(i,3))+(conj(E(p))-conj(E(q)))*(B1(i,5)/B1(i,3)));Siz(i)=Si(p,q);endSSi(p,q)=Si(p,q);Sline(i,3)=Siz(i);ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))];disp(ZF);enddisp('-----------------------------------------------------');disp('各条支路的末端功率Sj为(顺序同您输入B1时一致)：');for i=1:n1p=B1(i,1);q=B1(i,2);if B1(i,6)==0Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))));Sjy(i)=Sj(q,p);elseSj(q,p)=E(q)*(conj(E(q))*((B1(i,5)*(B1(i,5)-1))/B1(i,3))+(conj(E(q))-conj(E(p)))*(B1(i,5)/B1(i,3)));Sjy(i)=Sj(q,p);endSSj(q,p)=Sj(q,p);Sline(i,4)=Sjy(i);ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))];disp(ZF);enddisp('-----------------------------------------------------');disp('各条支路的功率损耗DS为(顺序同您输入B1时一致)：');for i=1:n1p=B1(i,1);q=B1(i,2);DS(i)=Si(p,q)+Sj(q,p);DDS(i)=DS(i);Sline(i,5)=DS(i);ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))];disp(ZF);enddisp('-----------------------------------------------------');disp('各支路首端编号末端编号首端功率末端功率线路损耗');disp(Sline);六、运行结果及其分析是否采用默认数据？(1-默认数据;2-手动输入)1导纳矩阵为：2.9056 -11.5015i 0.0000 + 5.3173i -1.6606 +3.1617i -1.2450 + 2.3710i0.0000 + 5.3173i 0.0000 - 4.6633i 0.0000 + 0.0000i 0.0000 + 0.0000i-1.6606 + 3.1617i 0.0000 + 0.0000i 2.4904 - 4.7039i -0.8298 + 1.5809i-1.2450 + 2.3710i 0.0000 + 0.0000i -0.8298 + 1.5809i 2.0749 - 3.9089i初始雅可比矩阵为：11.1267 2.7603 -5.3173 0 -3.1617 -1.6606-3.0509 11.8762 0 -5.3173 1.6606 -3.1617-5.3173 0 5.3173 0 0 00 -5.3173 0 4.0092 0 0-3.3198 -1.7436 0 0 4.8217 2.69800 0 0 0 0 2.1000迭代次数为：4收敛时电压修正量为：：1.0e-05 *0.0349-0.2445-0.0101-0.5713-0.0931-0.0073各节点的实际电压标幺值E为(节点号从小到大排列)：0.9673 - 0.0655i 1.0252 - 0.1666i 1.0495 - 0.0337i 1.0500 + 0.0000i -----------------------------------------------------各节点的电压大小V为(节点号从小到大排列)：0.9695 1.0387 1.0500 1.0500-----------------------------------------------------各节点的电压相角sida为(节点号从小到大排列)：-3.8734 -9.2315 -1.8419 0各节点的功率S为(节点号从小到大排列)：-0.0000 + 0.0000i -0.5000 - 0.3000i 0.2000 + 0.1969i 0.3277 + 0.0443i -----------------------------------------------------各条支路的首端功率Si为(顺序同您输入B1时一致)：S(1,2)=-0.5-0.30713iS(1,3)=-0.24266-0.197iS(1,4)=-0.25734-0.11013iS(3,4)=-0.055551+0.0017528i-----------------------------------------------------各条支路的末端功率Sj为(顺序同您输入B1时一致)：S(2,1)=0.5+0.24606iS(3,1)=0.25555+0.1952iS(4,1)=0.2712+0.1014iS(4,3)=0.056496-0.057061i-----------------------------------------------------各条支路的功率损耗DS为(顺序同您输入B1时一致)：DS(1,2)=0-0.06107iDS(1,3)=0.012892-0.0018014iDS(1,4)=0.013863-0.0087295iDS(3,4)=0.00094545-0.055308i-----------------------------------------------------各支路首端编号末端编号首端功率末端功率线路损耗1.0000 + 0.0000i2.0000 + 0.0000i -0.5000 - 0.3071i 0.5000 + 0.2461i 0.0000 - 0.0611i 1.0000 + 0.0000i3.0000 + 0.0000i -0.2427 - 0.1970i 0.2556 + 0.1952i 0.0129 - 0.0018i 1.0000 + 0.0000i4.0000 + 0.0000i -0.2573 - 0.1101i 0.2712 + 0.1014i 0.0139 - 0.0087i3.0000 + 0.0000i4.0000 + 0.0000i -0.0556 + 0.0018i 0.0565 - 0.0571i 0.0009 - 0.0553i七、实验体会及感悟通过这次实验，首先让我对matlab软件有了初步的了解，对它强大的矩阵运算能力有了更深的体会，同时掌握了设置断点和断点调试的一般方法，结合课本上的程序流程图和参考资料上的例子单步跟踪调试，再一次的熟悉了牛顿拉夫逊法潮流计算的一般方法和步骤，对计算机计算潮流计算有了更进一步的认识，在学习潮流计算时，虽然依次学习了节点导纳矩阵，功率方程、雅可比矩阵，但不能将它们联系起来，更不知道其中的原委，通过程序的编写，知道了其中的联系，也知道了每个方程、矩阵在计算中的作用。
电力系统正常运行潮流计算
~
1、某发电厂高压 110kV 母线上有 S L = 20 + J15MVA 负荷，试以（1）串联电阻
和电抗（2）并联电导和电纳来表示该负荷。
*
~
解： S
.*
=V I
.
=V
V
*
=
V2
*
ZZ
→Z
=
V2
*
S
ZL
=
VL
*
2
SL
= 1102 20 − j15
=
1102 25e− j36.9
110kV
110kV
RL
GL
jBL
jXL
阻抗表示的负荷形式
导纳表示的负荷形式
２、已知 10kV 系统某条线路的等值电路图如下，其首端母线电压为 11.05kV，试 做潮流分布，包括支路功率及节点电压。
V1 2+ j 6(Ω) V2 1+ j 3(Ω)
V3
2+ j 3(MVA) 解= 386 +
j290(Ω )
RL = 386(Ω ) ， X L = 290(Ω )
YL
=
1 ZL
=
1 484e j36.9
= 0.0027e− j36.9
= 0.00166 −
j0.00124(S )
GL = 0.00166(S ) ， BL = 0.00124(S ) 等值电路分别为：
= 9.86 × 2 +15.59 ×6 11.05
= 10.249kV
Create PDF with GO2PDF for free, if you wish to remove this line, click here to buy Virtual PDF Printer
电力系统分析(潮流计算)
电力系统分析（一）：电力系统的基本概念No.1电力系统的组成和接线方式1、电力系统的四大主要元件：发电机、变压器、电力线路、负荷。
2、动力系统包括动力部分（火电厂的锅炉和汽轮机、水电厂的水库和水轮机、核电厂的核反应堆和汽轮机）和电力系统。
3、电力网包括变压器和电力线路。
4、用户只能从一回线路获得电能的接线方式称为无备用接线方式。
No.2电力系统的运行特点1、电能的生产、传输、分配和消费具有：①重要性、②快速性、③同时性。
2、电力系统运行的基本要求：①安全可靠持续供电（首要要求）、②优质、③经济3、根据负荷的重要程度（供电可靠性）将负荷分为三级。
4、电压质量分为：①电压允许偏差、②三相电压允许不平衡度、③公网谐波、④电压允许波动与闪变5、衡量电能质量的指标：①电压、②频率、③波形（电压畸变率）6、10kV公用电网电压畸变率不超过4%。
7、抑制谐波的主要措施：①变压器星三角接线、②加装调谐波器、③并联电容/串联电抗、④增加整流器的脉冲次数8、衡量电力系统运行经济性的指标：①燃料损耗率、②厂用电率、③网损率9、线损包括：①管理线损、②理论线损、③不明线损10、线损计算方法：①最大负荷损耗时间法②最大负荷损失因数法③均方根电流法No.3电力系统的额定频率和额定电压1、电力线路的额定电压（也称电力网的额定电压）与用电设备的额定电压相同。
2、正常运行时电力线路首端的运行电压常为用电设备额定电压的105%，末端电压为额定电压。
3、发电机的额定电压比电力网的额定电压高5%。
4、变压器的一次绕组相当于用电设备，其额定电压与电力线路的额定电压相同；但变压器直接与发电机相连时，其额定电压与发电机额定电压相同，即为该电压级额定电压的105%。
5、变压器的二次绕组相当于电源，其输出电压应较额定电压高5%，但因变压器本身漏抗的电压损耗在额定负荷时约为5%，所以变压器二次侧的额定电压规定比额定电压高10%。
6、降压变压器二次侧连接10kV线路，当短路电压百分比小于7.5%（变压器本身漏抗的电压损耗较小）时，比线路额定电压高5%。
现代电力系统分析潮流计算2
现代电力系统分析潮流计算2现代电力系统分析潮流计算2潮流计算是电力系统分析中的一项重要工作，其目的是确定电力系统中各节点的电压和功率信息。
在现代电力系统中，潮流计算是实现电力系统的稳态分析和规划的基础工作。
潮流计算可以帮助系统操作员确定电网传输能力、系统稳定性等参数，对电力系统的运行和设计进行优化和改进。
潮流计算是通过求解潮流方程来得到各个节点的电压和功率。
潮流方程描述了电力系统中各个节点的电压和功率之间的关系。
潮流计算的基本原理是功率守恒原理，即系统输入功率等于输出功率。
通过潮流计算，可以得到电力系统中各节点的电压和功率信息，进而分析电力系统的稳定性、传输能力等指标。
直流潮流计算是最早出现的潮流计算方法，它采用直流模型对电力系统进行建模，忽略了电压相位的影响。
直流潮流计算的基本假设是电压相位小，即各节点之间的相位差非常小，因此可以忽略相位差，只考虑电压的模值变化。
直流潮流计算方法简单、迭代速度快，适用于简化的电力系统模型。
但由于忽略了相位差的影响，直流潮流计算在分析稳定性和传输能力等方面存在一定的局限性。
交流潮流计算是一种精确的潮流计算方法，它采用交流模型对电力系统进行建模，考虑了电压相位的影响。
交流潮流计算要求求解非线性方程组，通常使用牛顿-拉夫逊法、高斯-赛德尔法等迭代方法进行求解。
交流潮流计算方法适用于复杂的电力系统模型，可以对系统的稳定性和传输能力等进行精确分析。
为了提高潮流计算的速度和精度，现代电力系统还采用了一些改进的潮流计算方法，如快速潮流计算、修正潮流计算等。
快速潮流计算方法通过选择合适的系数矩阵进行近似计算，可以大大提高计算速度。
修正潮流计算方法通过修正潮流方程，可以减小误差，并提高计算结果的精度。
总之，潮流计算是现代电力系统分析中的一项重要工作。
通过潮流计算，可以得到电力系统中各节点的电压和功率信息，进而分析电力系统的稳定性、传输能力等指标。
潮流计算方法包括直流潮流计算和交流潮流计算等，还有一些改进的方法。
例题-第四章 电力系统潮流的计算机计算
第4章复杂电力系统的潮流计算一、填空题1.用计算机进行潮流计算时，按照给定量的不同，可将电力系统节点分为节点、节点、节点三大类，其中, 节点数目最多，节点数目很少、可有可无，节点至少要有一个.二、选择题1。
若在两个节点i、j之间增加一条支路，则下列关于节点导纳矩阵的说法正确的是（）A。
阶数增加1B.节点i的自导纳不变C。
节点i、j间的互导纳发生变化D。
节点j的自导纳不变2.若从节点i引出一条对地支路，则下列关于节点导纳矩阵的说法正确的是（)A.阶数增加1B.节点i的自导纳发生变化C。
节点i和其余节点间的互导纳均发生变化D。
节点导纳矩阵的所有元素均不变3。
若从两个节点i、j之间切除掉一条支路，则下列关于节点导纳矩阵的说法正确的是（)A.阶数减少1B.节点i、j间的互导纳一定变为0C。
节点i、j间的互导纳发生变化，但不一定变为0D.节点i、j的自导纳均不变4.若网络中增加一个节点k，且增加一条节点i与之相连的支路，则下列关于节点导纳矩阵的说法正确的是（）（1）阶数增加1(2)节点k的自导纳等于题干中所述支路的导纳(3)节点i的自导纳等于题干中所述支路的导纳（4）节点i、k间的互导纳等于题干中所述支路的导纳A.（1）（2）B。
(2)(3） C.（1）（4）D。
（2）（4）三、简答题1.什么是潮流计算?潮流计算的主要作用有哪些？潮流计算,电力学名词，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下,计算有功功率、无功功率及电压在电力网中的分布。
潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题.对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求;对运行中的电力系统,通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(线路、变压器等)是否会出现过负荷,以及可能出现过负荷时应事先采取哪些预防措施等。
电力系统潮流计算问答题
电力系统潮流计算问答题潮流计算数学模型与数值方法1. 什么是潮流计算？潮流计算的主要作用有哪些？潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。
对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。
对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。
潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。
2. 潮流计算有哪些待求量、已知量？（已知量：1、电力系统网络结构、参数2、决定系统运行状态的边界条件待求量：系统稳态运行状态例如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等）通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。
待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。
3. 潮流计算节点分成哪几类？分类根据是什么？（分成三类：PQ 节点、PV 节点和平衡节点，分类依据是给定变量的不同）PU 节点（电压控制母线）有功功率Pi 和电压幅值Ui 为给定。
这种类型节点相当于发电机母线节点，或者相当于一个装有调相机或静止补偿器的变电所母线。
PQ 节点注入有功功率Pi 和无功功率Qi 是给定的。
相当于实际电力系统中的一个负荷节点，或有功和无功功率给定的发电机母线。
平衡节点用来平衡全电网的功率。
平衡节点的电压幅值Ui 和相角δi 是给定的，通常以它的相角为参考点，即取其电压相角为零。
一个独立的电力网中只设一个平衡节点。
4. 教材牛顿-拉夫逊法及有功-无功分解法是基于何种电路方程？可否采用其它类型方程？基于节点电压方程，还可以采用回路电流方程和割集电压方程等。
但是后两者不常用。
5. 教材牛顿-拉夫逊法是基于节点阻抗方程、还是基于节点导纳方程进行迭代计算的？试阐述这两种方程的优点与缺点。
电力系统潮流计算问答题
电力系统潮流计算问答题潮流计算数学模型与数值方法1. 什么是潮流计算？潮流计算的主要作用有哪些？潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。
对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。
对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。
潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。
2. 潮流计算有哪些待求量、已知量？（已知量：1、电力系统网络结构、参数2、决定系统运行状态的边界条件待求量：系统稳态运行状态例如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等）通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。
待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。
3. 潮流计算节点分成哪几类？分类根据是什么？（分成三类：PQ节点、PV 节点和平衡节点，分类依据是给定变量的不同）PU节点（电压控制母线）有功功率Pi和电压幅值Ui为给定。
这种类型节点相当于发电机母线节点，或者相当于一个装有调相机或静止补偿器的变电所母线。
PQ节点注入有功功率Pi和无功功率Qi是给定的。
相当于实际电力系统中的一个负荷节点，或有功和无功功率给定的发电机母线。
平衡节点用来平衡全电网的功率。
平衡节点的电压幅值Ui和相角δi 是给定的，通常以它的相角为参考点，即取其电压相角为零。
一个独立的电力网中只设一个平衡节点。
4. 教材牛顿-拉夫逊法及有功-无功分解法是基于何种电路方程？可否采用其它类型方程？基于节点电压方程，还可以采用回路电流方程和割集电压方程等。
但是后两者不常用。
5. 教材牛顿-拉夫逊法是基于节点阻抗方程、还是基于节点导纳方程进行迭代计算的？试阐述这两种方程的优点与缺点。