BPA学习笔记

BPA潮流程序

【Problem】

【基础概念】

英文名词翻译：

Branch 支路

Bus 节点

Network 网络

Power flow 潮流

Base 基准

AI=area interchange 区域功率交换

名词释义：

缺省值：default value=默认值

输入、输出相关文件：

*.dat 潮流计算数据文件

*.bse潮流计算二进制结果文件（可用于潮流计算的输入或稳定计算）

*.pfo潮流计算结果文件（文本文件）

*.map供单线图格式潮流图及地理接线图格式潮流图程序使用的二进

制结果文件

*.pff，*.pfd 中间文件（正常计算结束后将自动删除。不正常时，将留在硬盘上，可随时删除）

pwrflo.dis 储存一个潮流作业计算时屏幕显示的信息。

pfcard.def 定义潮流程序卡片格式文件，用户可更改及调整该文件。该文件安装时放在与潮流程序相同的目录中。

[.dat]

潮流数据文件：1）网络数据（节点支路参数）

2）控制语句

【控制语句】

一级（……）

二级/……\

三级>……<

*必须顶格填写

*除了说明书注明外，控制语句写在网络数据前面

*空格=“，”—>名称中含空格的用“#”代替

*长度限制<100列—>超过改为多行输入：“-”表示断行，下一行继续此前内容

*“_”是无用字符，可以随意添加

*任意一行数据前面添加“.”表示该行无效（注释行）

常用控制语句

（POWERFLOW,CASEID=潮流方式名，PROJECT=工程名）

***表示潮流作业开始，必须是潮流数据文件的第一条控制语。

（END）

***表示潮流作业结束，必须是潮流数据的最后一行。该语句后面的控制语句和数据均无效。

----------------------------------------------------------------------

！计算过程的控制！用什么方法算潮流（收敛误差、带载T、边界P控制的选择）HOW

．．．？．/SOLUTION\

>AI_CONTROL=CON< ---表控制（缺省值）

MON< ---表不控制但监视

OFF< ---不控制

>LTC=ON<

ON_NV<

ON_NPS<

OFF<

ON_DCONL Y<

>SOL_ITER,DECOUPLED=(NUMBER),CURRENT=(NUMBER),NEWTOW=(NUMBER)< ***改组控制语句必须放在网络数据后面。

***上述第三级控制语句可以改写为第二级，必须放在网络数据前面。

----------------------------------------------------------------------

！用于指定【调用&输出】的文件的语句！输到哪里？WHERE?

．．．．．．

/INCLUDE_CONTROL,FILE=文件名\

/BRANCH_DATA,FILE=文件名\

/BUS_DA TA,FILE=文件名\

/NETWORK_DATA,FILE=文件名{or * },RXCHECK=ON{or OFF}\

/NEW_BASE,FILE=文件名\

/OLD_BASE,FILE=文件名\

/PF_MAP,FILE=文件名\

***可以输出的文件有①新库文件new base；②潮流图pf map；③分析报告analysis ptr ***可以输出的内容有①原始数据input；②潮流结果output

----------------------------------------------------------------------

/MV A_BASE=100\

-------------------

！输出选择（输出什么？）！WHAT?

．．．．．

/P_INPUT_LIST,NONE\

FULL,ERRORS=NO_LIST\

LIST\

ZONES=分区名 \

/PTR_SORT=BUS\

ZONE\

AREA\

/P_OUTPUT_LIST,NONE,FAILED_SOL=FULL_LIST\

FULL, PARTIAL_LIST\

ZONES=分区名，NO_LIST\

/P_ANAL YSIS_PTR,LEVEL=2，*\

1,ZONES=分区名\

3,OWNERS=所有者名\

4\

/ANL YSIS_SELECT\

>PAPER,ZONES=分区名,AREAS=区域名,OWNERS=所有者名< 输入范围指定

>USERAN< （输出用户自定义分析报告）（用户想输出的表目的缩写名称，见p19）>UNSCH< （输出系统未安排无功功率节点一览表）

……

/AI_LIST=FULL\

MATRIX\

TIELINE\

NONE\

-----------------------------------------

用户自定义控制语句

——自己定义输出的数据和格式

/USER_ANALYSIS,FILE=XXX,OUTPUT=XXX\

***该语句配有第三季控制语句——

①FILE缺省or标*，紧跟此语句后面输入

②制定了FILE文件名，则调用指定文件（含第三级控制语句）

>DEFINE_TYPE XXX<

LET 变量名=

***一般先定义物理量or分区，再定义计算，然后按自己定义的格式输出。

1、定义——物理量之间的数学运算

>DEFINE_TYPE FUCTION<

LET X=AR+BL（数学运算）

2、定义——输电线路潮流

>DEFINE_TYPE BRANCH_P<

LET 变量名=线路节点1名称，节点1基准电压，节点2名称，节点2基准电压，线路回路号

>DEFINE_TYPE BRANCH_Q<

LET 变量名=线路节点1名称，节点1基准电压，节点2名称，节点2基准电压，线路回路号

3、定义——区域联络线交换功率

4、定义——网损

5、定义——节点物理量

6、定义——分区物理量

7、输出格式

H XXX

C XXX $AAA

修改潮流数据控制语句

1、修改——数据控制语句

/CHANGES,FILE=*\

(文件名)\

数据的标志字符修改码节点名称改变后数据

修改码含义

空白增加新数据

D 删除数据

M 修正数据

R 恢复数据

2、修改——节点类型控制语句

/CHANGE_BUS_TYPE,BQ=B,BG=BQ,BT=B,BX=B,BQ=BF,BG=BF...,AREAS=,ZONES=\ >EXCLUDE_BUS<

修改区域内不想做修改的节点名称

【数据卡片】

P68

分类1）区域控制AC A AO I

2）节点数据 B BD BM + X

3）支路数据L LD LM T R E RZ

4）节点数据修改卡P Z DZ

区域控制数据卡

a)AC——区域交换功率控制卡

功能：将区域AREA之间的交换功率控制为用户给定值

***区域输出功率由区域缓冲节点控制（每个AREA必须配有一个区域缓冲节点）全网平衡节点在本区域内，将其作为区域缓冲节点；不在选一调频电厂做（PV）

***区域交换功率=区域边界测得所有区域联络线输送功率和

***每个AREA必须各填一张AC卡-->一张AC最多20个ZONE

一个AREA包括的ZONE>20，填AC+卡，最多一张AC+,AC+最多20个ZONE

b) AO——区域输出分类卡

c) I——安排区域交换功率数据卡

节点数据卡

a)B——交流节点数据卡

*所有节点（含无源节点）都要填写

b)BD

c)BM

d)+ ——延续节点交流卡

*作用：补充制定交流节点的数据（恒阻抗负荷、恒电流负荷、恒功率负荷、发电出力etc）

e)X

Example

1、IEEE9节点【.dat】

(POWERFLOW,CASEID=IEEE9,PROJECT=IEEE_9BUS_TEST_SYSTEM)

/SOL_ITER,DECOUPLED=2,NEWTON=15,OPITM=0\-----计算过程的控制

/P_INPUT_LIST,ZONES=NONE\--------------------------------输出什么

/P_OUTPUT_LIST,ZONES=ALL\

/RPT_SORT=ZONE\

/NEW_BASE,FILE=IEEE90.BSE\-------------------------------输出到哪

/PF_MAP,FILE = IEEE90.MAP\

/P_ANAL YSIS_REPROT,LEVEL=4\---------------------------输出报告选择

/NETWORK_DATA\-----------------------------------------------网络数据来源如下

.B ----------------- Bus ----------------------------

BS 发电机1 16.501 999. 999. 1.00

BQ 发电机2 18.002 17. 8. 400.180. 120. -30. 1.00

BQ 发电机3 13.801 320. 85. 80. 10. 1.00

B 母线1 230.01 20.

B 母线2 230.02 35.0 10.0 1.0 1.0

B 母线3 230.01

.B 母线A 230.02 125. 70.0 40.

.B 母线C 230.02 100. 55.0 40.

.B 母线B 230.01 90. 40.0 20.

B 母线A 230.02 110. 70.0 20.

B 母线

C 230.02 100. 55.0 20.

B 母线B 230.01 90. 40.0 10.

.L ----------------- transmission lines ----------------------------

L 母线1 230. 母线A 230. 180. .0100 .0850 .0440

L 母线1 230. 母线B 230. 180. .0170 .0920 .0395

L 母线A 230. 母线2 230. 180. .0320 .1610 .0765

L 母线B 230. 母线3 230. 180. .0390 .1700 .0895

L 母线2 230. 母线C 230. 180. .0085 .0720 .03725

L 母线C 230. 母线3 230. 180. .0119 .1008 .05225

.T ----- transformers ---------

T 发电机1 16.5 母线1 230. .0 .0567 16.5 242.

T 发电机2 18.0 母线2 230. .0625 18.0 242.

.R 发电机2 18.02母线2 230. 母线2 230.242. 218. 12

T 发电机3 13.8 母线3 230. .001 .0586 13.8 242. ./CHANGES\

.PZ 01 1.5 1.5 1.5 1.5

.PZ 02 1.5 1.5 1.5 1.5

(END)

电力系统运行方式及潮流分析实验报告

电力系统运行方式及潮 流分析实验报告 Company Document number：WUUT-WUUY-WBBGB-BWYTT-1982GT

电力系统第一次实验报告——电力系统运行方式及潮流分析实验

实验1 电力系统运行方式及潮流分析实验 一、实验目的 1、掌握电力系统主接线电路的建立方法 2、掌握辐射形网络的潮流计算方法； 3、比较计算机潮流计算与手算潮流的差异； 4、掌握不同运行方式下潮流分布的特点。 二、实验内容 1、辐射形网络的潮流计算； 2、不同运行方式下潮流分布的比较分析 三、实验方法和步骤 1．辐射形网络主接线系统的建立 输入参数（系统图如下）： G1：300+j180MVA（平衡节点） 变压器B1：Sn=360MVA，变比=18/121，Uk％=％，Pk=230KW，P0=150KW，I0/In=1％； 变压器B2、B3：Sn=15MVA，变比=110/11 KV，Uk％=％，Pk=128KW， P0=，I0/In=％； 负荷F1：20+j15MVA；负荷F2：28+j10MVA； 线路L1、L2：长度：80km，电阻：Ω/km，电抗：Ω/km，电纳：×10-6S/km。 辐射形网络主接线图 （1）在DDRTS中绘出辐射形网络主接线图如下所示： （2）设置各项设备参数： G1：300+j180MVA（平衡节点） 变压器B1：Sn=360MVA，变比=18/121，Uk％=％，Pk=230KW，P0=150KW，I0/In=1％；

变压器B2、B3：Sn=15MVA，变比=110/11 KV，Uk％=％，Pk=128KW， P0=，I0/In=％； 负荷F1：20+j15MVA；负荷F2：28+j10MVA； 线路L1、L2：长度：80km，电阻：Ω/km，电抗：Ω/km，电纳：×10-6S/km。2．辐射形网络的潮流计算 （1）调节发电机输出电压，使母线A的电压为115KV，运行DDRTS进行系统潮流计算，在监控图页上观察计算结果 项目DDRTS潮流计算结果 变压器B2输入功率+ 变压器B2输出功率+ 变压器B3输入功率+ 变压器B3输出功率+ 线路L1输入功率+ 线路L1输出功率+ 线路L2输入功率+ 线路L2输出功率+ （2）手算潮流： （3）计算比较误差分析 通过比较可以看出，手算结果与计算机仿真结果相差不大。产生误差原因：手算时是已知首端电压、末端功率的潮流计算，计算过程中要将输电线路对地电容吸收的功率以及变压器励磁回路吸收的功率归算到运算负荷中，并且在每一轮的潮流计算中都用上一轮的电压或功率的值（第一轮电压用额定电压）。 3．不同运行方式下潮流比较分析 （1）实验网络结构图如上。由线路上的断路器切换以下实验运行方式： ①双回线运行（L1、L2均投入运行） ②单回线运行（L1投入运行，L2退出）将断路器断开 对上述两种运行方式分别运行潮流计算功能，将潮流计算结果填入下表：

关于电力系统经济调度的潮流计算分析

关于电力系统经济调度的潮流计算分析 发表时间：2016-05-24T15:57:29.347Z 来源：《电力设备》2016年第2期作者：秦先威 [导读] （国网山东省电力公司烟台市牟平区供电公司山东烟台 264100）随着经济的快速发展和科技的不断进步，社会各行业对电力资源的需求量越来越大，我国的电力系统建设规模也越来越大。 （国网山东省电力公司烟台市牟平区供电公司山东烟台 264100） 摘要:潮流计算是电力调度中最重要也是最基本的计算之一,它应用于电力系统中实时电价计算、输电权分配、网络阻塞管理等多方面。 关键词:电力系统；经济调度；潮流计算 前言 随着经济的快速发展和科技的不断进步，社会各行业对电力资源的需求量越来越大，我国的电力系统建设规模也越来越大。电力调度对电力系统的正常运行有很大的影响，而潮流计算则是电力调度中最重要的基本计算方法，潮流计算对电价计算、输电分配、电网线路管理有十分重要的影响。随着经济的快速发展，我国的电力企业得到了飞速的发展，与此同时，人们对供电质量的要求也越来越高，为满足人们的用电需求，电力系统在运行过程中，必须保证电力调度的合理性、科学性，潮流计算是电力系统经济调度最重要的计算方法之一，潮流计算的结果准确性很高，科学性很强，潮流计算对电力系统经济调度有十分重要的作用。 一、潮流计算的概述 1.1 潮流计算的概述 潮流计算是指利用已知的电网接线方式、参数、运行条件，将电力系统的各个母线电压、支路电流、功率、网损计算出来。通过潮流计算能判断出正在运行的电力系统的母线电压、支路电流、功率是否在允许范围内运行，如果超出允许范围，就需要采用合理的措施，对电力系统的进行方式进行调整。在电力系统规划过程中，采用潮流计算，能为电网供电方案、电气设备的选择提供科学的依据，同时潮流计算还能为自动装置定整计算、继电保护、电力系统稳定计算、故障计算提供原始数据。 1.2 潮流计算的电气量 潮流计算是根据电力系统接线方式、运行条件、参数等已知条件，将稳定状态下电力系统的电气量计算出来。一般情况下，给出的条件有电源、负荷节点的功率、平衡节点的电压、相位角、枢纽点的电压，需要计算的电气量有各节点的电压、相位角、各支路通过的电流、功率、网络的功率损耗等。 1.3 传统的潮流计算方法 传统的潮流计算方法，包括很多不同的内容，具有一定的优点和缺点。例如，传统的潮流计算方法，包括非线性规划法、二次规划法和线性规划法等。在电力系统经济调度的过程中，应用传统的潮流计算方法，优点是：可以根据目标函数的导数信息，确定需要进行搜索的方向，因此在计算的时候，具有较快的速度和清晰的计算过程。而且，可信度比较高。 1.5 智能的潮流计算方法 潮流计算中人工智能方法的优点是：随机性：属于全局优化算法，跳出局部极值点比较容易；与导数无关性：在工程中，一些优化问题的目标函数处于不可导状态。如果进行近似和假设，会对求解的真实性造成影响；内在并行性：操作对象为一组可行解，在一定程度上可以克服内在并发性开放中性能的不足。而其缺点，主要是：需要按照概率进行操作，不能保证可以完全获取最优解；算法中的一些控制参数需要根据经验人文地给出，对专家经验和一定量的试验要求比较高；表现不稳定，在同一问题的不同实例中应用算法会出现不同的效果。 二、潮流计算的分类 根据电力系统的运行状态，潮流计算可以分为离线计算和在线计算两种方法，离线计算主要用于电力系统规划设计和电力系统运行方式安排中；在线计算主要用于电力系统运行监控和控制中；根据潮流计算的发展，潮流计算可以分为传统方法和人工智能方法两种情况，下面分别对这两种方法进行分析。 2.1 潮流计算的传统方法 潮流计算的传统方法有非线性规划法、线性规划法、二次规划法等几种情况，潮流计算的传统方法具有计算速度快、解析过程清晰、结果真实可靠等优点，但传统方法对目标函数有一定的限制，需要简化处理，这样求出来的值有可能不是最优值。 2.2 潮流计算的人工智能方法 潮流计算的人工智能方法是一种新兴的方法，人工智能方法不会过于依赖精确的数学模型，它有粒子群优化算法、遗传法、模拟退火法等几种情况，人工智能方法的计算结果和导数没有关系，其操作对象是一组可行解，能克服内在并行性存在的问题，但人工智能方法表现不太稳定，在计算过程中，有的控制参数需要根据经验得出，因此，采用人工智能方法进行计算时，需要计算人员有丰富的经验。 三、潮流计算在电力系统经济调度中的应用 3.1 在输电线路线损计算的应用 在进行输电线路线损计算过程中，通过潮流计算能得出经济潮流数据。潮流程度能根据线路的功率因数、有功负荷、无功负荷等参数，计算出潮流线损，例如一条长为38.1km，型号为LGJ—150的导线，当潮流为20MW、功率因数为0.9时，该线路线损为0.24MW，线损率为1.18%；当潮流为30MW、功率因数为0.9时，该线路线损为0.57MW，线损率为1.91%；潮流为50MW、功率因数为0.9时，该线路线损为1.95MW，线损率为3.90%；由此可以看出，潮流小于30MW时，线损率小于2%，潮流超过50MW时，线损率将超过4%，因此，该输电线路的经济输送潮流为30MW以下。调度人员可以根据计算结果，编制线路经济运行方案，从而实现节能调度。 3.2 在变压器变损中的应用 调度人员可以利用潮流计算程序，将变压器在不同负荷下的损耗、变损率计算出来，从而为变压器控制提供依据。例如一台40MVA双

基于内点法的最优潮流计算

基于内点法的最优潮流计 算 Prepared on 24 November 2020

摘要 内点法是一种能在可行域内部寻优的方法，即从初始内点出发，沿着中心路径方向在可行域内部直接走向最优解的方法。其中路径跟踪法是目前最具有发展潜力的一类内点算法，该方法鲁棒性强，对初值的选择不敏感，在目前电力系统优化问题中得到了广泛的应用。本文采用路径跟踪法进行最优求解，首先介绍了路径跟踪法的基本模型，并且结合具体算例，用编写的Matlab程序进行仿真分析，验证了该方法在最优潮流计算中的优越性能。 关键词：最优潮流、内点法、路径跟踪法、仿真

目次

0、引言 电力系统最优潮流，简称OPF（Optimal Power Flow）。OPF问题是一个复杂的非线性规划问题，要求满足待定的电力系统运行和安全约束条件下，通过调整系统中可利用控制手段实现预定目标最优的系统稳定运行状态。针对不同的应用，OPF模型课以选择不同的控制变量、状态变量集合，不同的目标函数，以及不同的约束条件，其数学模型可描述为确定一组最优控制变量u，以使目标函数取极小值，并且满足如下等式和不等式。 {min u f(x,u) S.t.?(x,u)=0 g(x,u)≤0 （0-1）其中min u f(x,u)为优化的目标函数，可以表示系统运行成本最小、或者系统运行网损最小。S.t.?(x,u)=0为等式约束，表示满足系统稳定运行的功率平衡。g(x,u)≤0为不等式约束，表示电源有功出力的上下界约束、节点电压上下线约束、线路传输功率上下线约束等等。 电力系统最优潮流算法大致可以分为两类：经典算法和智能算法。其中经典算法主要是指以简化梯度法、牛顿法、内点法和解耦法为代表的基于线性规划和非线性规划以及解耦原则的算法，是研究最多的最优潮流算法, 这类算法的特点是以一阶或二阶梯度作为寻找最优解的主要信息。智能算法主要是指遗传算法和模拟退火发等，这类算法的特点是不以梯度作为寻优信息，属于非导数的优化方法。 因此经典算法的优点是能按目标函数的导数信息确定搜索方向，计算速度快，算法比较成熟，结果可信度高。缺点是对目标函数及约束条件有一定的限

武汉大学电力系统分析实验报告

电气工程学院 《电力系统分析综合实验》2017年度PSASP实验报告 学号： 姓名： 班级：

实验目的： 通过电力系统分析的课程学习，我们都对简单电力系统的正常和故障运行状态有了大致的了解。但电力系统结构较为复杂，对电力系统极性分析计算量大，如果手工计算，将花费 大量的时间和精力，且容易发生错误。而通过使用电力系统分析程序PSASP，我们能对电 力系统潮流以及故障状态进行快速、准确的分析和计算。在实验过程中，我们能够加深对电力系统分析的了解，并学会了如何使用计算机软件等工具进行电力系统分析计算，这对我们以后的学习和工作都是有帮助的。 潮流计算部分： 本次实验潮流计算部分包括使用牛顿法对常规运行方式下的潮流进行计算，以及应用PQ分解法规划运行方式下的潮流计算。在规划潮流运行方式下，增加STNC-230母线负荷的有功至1.5.p.u，无功保持不变，计算潮流。潮流计算中，需要添加母线并输入所有母线 的数据，然后再添加发电机、负荷、交流线、变压器、支路，输入这些元件的数据。对运行方案和潮流计算作业进行定义，就可以定义的潮流计算作业进行潮流计算。 因为软件存在安装存在问题，无法使用图形支持模式，故只能使用文本支持模式，所以 无法使用PSASP绘制网络拓扑结构图，实验报告中的网络拓扑结构图均使用Visio绘制， 请见谅。 常规潮流计算： 下图是常规模式下的网络拓扑结构图，并在各节点标注电压大小以及相位。 下图为利用复数功率形式表示的各支路功率（参考方向选择数据表格中各支路的i侧母

线至j侧），因为无法使用图形支持模式，故只能通过文本支持环境计算出个交流线功率，下图为计算结果。

电力系统分析潮流实验报告(优选.)

最新文件---------------- 仅供参考--------------------已改成-----------word文本 --------------------- 方便更改 赠人玫瑰，手留余香。 南昌大学实验报告 学生姓名：学号：专业班级： 实验类型：□验证□综合■设计□创新实验日期：实验成绩： 电力系统潮流计算实验 一、实验目的： 本实验通过对电力系统潮流计算的计算机程序的编制与调试，获得对复杂电力系统进行潮流计算的计算机程序，使系统潮流计算能够由计算机自行完成，即根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，由计算程序运行完成该电力系统的潮流计算。通过实验教学加深学生对复杂电力系统潮流计算计算方法的理解，学会运用电力系统的数学模型，掌握潮流计算的过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。 二、实验内容： 编制调试电力系统潮流计算的计算机程序。程序要求根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，完成该电力系统的潮流计算，要求计算出节点电压、功率等参数。 1、在各种潮流计算的算法中选择一种，按照计算方法编制程序。 2、将事先编制好的电力系统潮流计算的计算程序原代码由自备移动存储设备导入计算机。

3、在相应的编程环境下对程序进行组织调试。 4、应用计算例题验证程序的计算效果。 三、实验程序： function [e,f,p,q]=flow_out(g,b,kind,e,f) %计算潮流后efpq的终值 s=flow(g,b,kind,e,f); k=0; while max(abs(s))>10^-5 J=J_out(g,b,kind,e,f); J_ni=inv(J); dv=J_ni*s; l=length(dv)/2; for i=1:l e(i)=e(i)-dv(2*i-1); f(i)=f(i)-dv(2*i); end s=flow(g,b,kind,e,f); end l=length(e); for i=1:l s1=0; s2=0; for j=1:l s1=s1+g(i,j)*e(j)-b(i,j)*f(j); s2=s2+g(i,j)*f(j)+b(i,j)*e(j); end p(i)=e(i)*s1+f(i)*s2; q(i)=f(i)*s1-e(i)*s2;

潮流计算的基本算法及其使用方法

潮流计算的基本算法及使用方法 一、 潮流计算的基本算法 1. 牛顿－拉夫逊法 1．1 概述 牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。这种方法的特点就是把对非线 性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。 牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏 导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。 1．2 一般概念 对于非线性代数方程组 ()0=x f 即 ()0,,,21=n i x x x f Λ ()n i Λ,2,1= (1－1) 在待求量x 的某一个初始计算值() 0x 附件，将上式展开泰勒级数并略去二阶及以上的高 阶项，得到如下的线性化的方程组 ()()()() ()0000=?'+x x f x f (1－2) 上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 ()() ()[]()()0 1 00x f x f x -'-=? (1－3) 将() 0x ?和() 0x 相加，得到变量的第一次改进值()1x 。接着再从() 1x 出发，重复上述计算 过程。因此从一定的初值() 0x 出发，应用牛顿法求解的迭代格式为 ()()()()() k k k x f x x f -=?' (1－4) ()()()k k k x x x ?+=+1 (1－5) 上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代

潮流计算问答题

1.什么是潮流计算？潮流计算的主要作用有哪些？ 潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。 对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置整定计算、电力系统故障计算和稳定计算等提供原始数据。 2.潮流计算有哪些待求量、已知量？ （已知量： 电力系统网络结构、参数； 决定系统运行状态的边界条件 待求量：系统稳态运行状态 例如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等）通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。 待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。 3.潮流计算节点分成哪几类？分类根据是什么？ （分成三类：PQ节点、PV节点和平衡节点，分类依据是给定变量的不同） PV节点（电压控制母线）：有功功率Pi和电压幅值Ui为给定。这种类型节点相当于发电机母线节点，或者相当于一个装有调相机或静止补偿器的变电所母线。 PQ节点：注入有功功率Pi和无功功率Qi是给定的。相当于实际电力系统中的一个负荷节点，或有功和无功功率给定的发电机母线。 平衡节点：用来平衡全电网的功率。平衡节点的电压幅值Ui和相角δi是给定的，通常以它的相角为参考点，即取其电压相角为零。 一个独立的电力网中只设一个平衡节点。 4.教材牛顿-拉夫逊法及有功-无功分解法是基于何种电路方程？可否采用其它类型方程？ 基于节点电压方程，还可以采用回路电流方程和割集电压方程等。但是后两者不常用。

潮流计算(matlab)实例计算

潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新

的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图

潮流计算实验

电力系统分析实验报告 实验一：潮流计算的计算机算法 >> clear; n=10; nl=10; isb=1; pr=0.00001; B1=[120.03512+0.08306i0.13455i10; 230.0068+0.18375i0 1.023811; 140.05620+0.13289i0.05382i10; 450.00811+0.24549i0 1.023811; 160.05620+0.13289i0.05382i10; 460.04215+0.09967i0.04037i10; 670.0068+0.18375i0 1.023811; 680.02810+0.06645i0.10764i10; 8100.00811+0.24549i011; 890.03512+0.08306i0.13455i10] B2=[00 1.1 1.101; 001002; 00.343+0.21256i1002; 001002; 00.204+0.12638i1002; 001002; 00.306+0.18962i1002; 001002; 0.50 1.1 1.103; 00.343+0.21256i1002] Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); for i=1:nl if B1(i,6)==0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;

有关电力系统三种潮流计算方法的比较.docx

电 力 系 统 三 种 潮 流 计 算 方 法 的 比 较 一、高斯 -赛德尔迭代法： 以导纳矩阵为基础， 并应用高斯 -- 塞德尔迭代的算法是在电力系统中最早得到应用的潮流计算方法，目前高斯一塞德尔法已很少使用。 将所求方程 f ( x ) 0 改写为 x ( x ) 不能直接得出方程的根，给一个猜测值 x 0 得 x 1( x 0 ) 又可取 x1 为猜测值，进一步得： x 2 ( x 1 ) 反复猜测 x k 1 迭代 则方程的根 ( x k ) 优点： 1. 原理简单，程序设计十分容易。 2. 导纳矩阵是一个对称且高度稀疏的矩阵，因此占用内存非常节省。 3. 就每次迭代所需的计算量而言，是各种潮流算法中最小的，并且和网络所包 含的节点数成正比关系。 缺点： 1. 收敛速度很慢。 2. 对病态条件系统，计算往往会发生收敛困难：如节点间相位角差很大的重负 荷系统、包含有负电抗支路 (如某些三绕组变压器或线路串联电容等 )的系统、具有较长的辐射形线路的系统、长线路与短线路接在同一节点上，而且长短 线路的长度比值又很大的系统。 3. 平衡节点所在位置的不同选择，也会影响到收敛性能。 二、牛顿 -拉夫逊法： 求解 f ( x ) 0 设 x x 0 x ，则 按牛顿二项式展开： 当 △x 不大，则取线性化（仅取一次项） 则可得修正量 对 得： 作变量修正： x k 1x k x k ，求解修正方程 20 世纪 牛顿法是数学中求解非线性方程式的典型方法，有较好的收敛性。自从 60 年代中期采用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、计算速度方面都超过了其他方法，成为直到目前仍被广泛采用的方法。 优点： 1. 收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭 代 4—5 次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。 2. 具有良好的收敛可靠性， 对于前面提到的对以节点导纳矩阵为基础的高斯一 塞德尔法呈病态的系统，牛顿法均能可靠地收敛。 3. 牛顿法所需的内存量及每次迭代所需时间均较前述的高斯一塞德尔法为多， 并与程序设计技巧有密切关系。 缺点：

电力系统分析潮流计算例题

电力系统的潮流计算 西安交通大学自动化学院 2012.10 3．1 电网结构如图3—11所示，其额定电压为10KV 。已知各节点的负荷功率及参数： MVA j S )2.03.0(2 +=， MVA j S )3.05.0(3+=， MVA j S )15.02.0(4+= Ω+=)4.22.1(12j Z ，Ω+=)0.20.1(23j Z ，Ω+=)0.35.1(24j Z 试求电压和功率分布。 解：（1）先假设各节点电压均为额定电压，求线路始端功率。 0068.00034.0)21(103.05.0)(2 2223232232323j j jX R V Q P S N +=++=++=?0019.00009.0)35.1(10 15.02.0)(2 2 224242242424j j jX R V Q P S N +=++=++=?

则： 3068.05034.023323j S S S +=?+= 1519.02009.024424j S S S +=?+= 6587.00043.122423' 12 j S S S S +=++= 又 0346 .00173.0)4.22.1(106587.00043.1)(2 2 212122'12'1212j j jX R V Q P S N +=++=++=? 故： 6933.00216.112'1212 j S S S +=?+= （2） 再用已知的线路始端电压kV V 5.101 =及上述求得的线路始端功率 12 S ，求出线 路 各 点 电 压 。

kV V X Q R P V 2752.05 .104.26933.02.10216.1)(11212121212=?+?=+=? kV V V V 2248.101212=?-≈ kV V V V kV V X Q R P V 1508.100740.0) (24242 2424242424=?-≈?=+=? kV V V V kV V X Q R P V 1156.101092.0) (23232 2323232323=?-≈?=+=? （3）根据上述求得的线路各点电压，重新计算各线路的功率损耗和线路始端功率。 0066.00033.0)21(12.103.05.02 2 223j j S +=++=? 0018.00009.0)35.1(15 .1015.02.02 2 224j j S +=++=? 故 3066.05033.023323j S S S +=?+= 1518.02009.024424j S S S +=?+= 则 6584.00042.122423' 12 j S S S S +=++= 又 0331.00166.0)4.22.1(22 .106584.00042.12 2 212j j S +=++=? 从而可得线路始端功率 6915.00208.112 j S +=

电力系统潮流计算实验报告

1． 手算过程 已知： 节点1：PQ 节点， s(1)= -0.5000-j0.3500 节点2：PV 节点， p(2)=0.4000 v(2)=1.0500 节点3：平衡节点，U(3)=1.0000∠0.0000 网络的连接图： 0.0500+j0.2000 1 0.0500+j0.2000 2 3 1）计算节点导纳矩阵 由2000.00500.012j Z += ? 71.418.112j y -=; 2000.00500.013j Z += ? 71.418.113j y -=; ∴导纳矩阵中的各元素： 42.936.271.418.171.418.1131211j j j y y Y -=-+-=+=; 71.418.11212j y Y +-=-=; 71.418.11313j y Y +-=-=; =21Y 71.418.11212j y Y +-=-=; 71.418.12122j y Y -==; 002323j y Y +=-=; =31Y 71.418.11313j y Y +-=-=; =32Y 002323j y Y +=-=; 71.418.13133j y Y -==; ∴形成导纳矩阵B Y ： ?? ?? ? ?????-++-+-+-+-+--=71.418.10071.418.10071.418.171.418.171.418.171.418.142.936.2j j j j j j j j j Y B 2）计算各PQ 、PV 节点功率的不平衡量，及PV 节点电压的不平衡量： 取:000.0000.1)0(1)0(1)0(1j jf e U +=+= 000.0000.1)0(2) 0(2)0(2j jf e U +=+= 节点3是平衡节点，保持000.0000.1333j jf e U +=+=为定值。 ()()[] ∑==++-=n j j j ij j ij i j ij j ij i i e B f G f f B e G e P 1 )0()0()0()0()0()0() 0(；

完整word版潮流计算方法

由于本人参加我们电气学院的电气小课堂，主讲的是计算机算法计算潮流这章，所以潜心玩了一个星期，下面整理给大家分享下。 本人一个星期以来的汗水，弄清楚了计算机算法计算潮流的基础，如果有什么不懂的可以发信息到邮箱：zenghao616@https://www.wendangku.net/doc/ba16787551.html, 接下来开始弄潮流的优化问题，吼吼！ 电力系统的潮流计算的计算机算法：以MATLAB为环境 这里理论不做过多介绍，推荐一本专门讲解电力系统分析的计算机算法的书籍---------《电力系统分析的计算机算法》—邱晓燕、刘天琪编著。 这里以这本书上的例题【2-1】说明计算机算法计算的过程，分别是牛顿拉弗逊算法的直角坐标和极坐标算法、P-Q分解算法。主要是简单的网络的潮流计算，其实简单网络计算和大型网络计算并无本质区别，代码里面只需要修改循环迭代的N即可，这里旨在弄清计算机算法计算潮流的本质。代码均有详细的注释. 其中简单的高斯赛德尔迭代法是以我们的电稳教材为例子讲，其实都差不多，只要把导纳矩阵Y 给你，节点的编号和分类给你，就可以进行计算了，不必要找到原始的电气接线图。 理论不多说，直接上代码： 简单的高斯赛德尔迭代法： 这里我们只是迭代算出各个节点的电压值，支路功率并没有计算。 S_ij=P_ij+Q_ij=V_i(V_i* -V_j*) * y_ij* 可以计算出各个线路的功率 在显示最终电压幅角的时候注意在MATLAB里面默认的是弧度的形式，需要转化成角度显示。clear;clc; %电稳书Page 102 例题3-5 %计算网络的潮流分布 --- 高斯-赛德尔算法 %其中节点1是平衡节点 %节点2、3是PV节点，其余是PQ节点 % 如果节点有对地导纳支路 %需将对地导纳支路算到自导纳里面 %------------------------------------------------% %输入原始数据，每条支路的导纳数值，包括自导和互导纳； y=zeros(5,5); y(1,2)=1/(0.0194+0.0592*1i); y(1,5)=1/(0.054+0.223*1i); y(2,3)=1/(0.04699+0.198*1i); y(2,4)=1/(0.0581+0.1763*1i); 由于电路网络的互易性，导纳矩阵为对称的矩阵%． for i=1:1:5 for j=1:1:5

2021年潮流计算的基本算法及使用方法

潮流计算的基本算法及使用方法 一、 欧阳光明（2021.03.07） 二、 潮流计算的基本算法 1. 牛顿－拉夫逊法 1．1 概述 牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。这种方法的特点就是把对非线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程，就是牛顿－拉夫逊法的核心。 牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。因为越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。而所谓“某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点，一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。 1．2一般概念 对于非线性代数方程组 即 ()0,,,21=n i x x x f ()n i ,2,1= (1－1) 在待求量x 的某一个初始计算值()0x 附件，将上式展开泰勒级数并

略去二阶及以上的高阶项，得到如下的线性化的方程组 ()()()() ()0000=?'+x x f x f (1－2) 上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 ()() ()[]()()0 1 00x f x f x -'-=? (1－3) 将()0x ?和()0x 相加，得到变量的第一次改进值()1x 。接着再从()1x 出 发，重复上述计算过程。因此从一定的初值()0x 出发，应用牛顿法求解的迭代格式为 ()()()() ()k k k x f x x f -=?' (1－4) ()()()k k k x x x ?+=+1 (1－5) 上两式中：()x f '是函数()x f 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ；k 为迭代次数。 由式（1－4）和式子（1－5）可见，牛顿法的核心便是反复形成 求解修正方程式。牛顿法当初始估计值()0x 和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。 1．3潮流计算的修正方程 运用牛顿－拉夫逊法计算潮流分布时，首先要找出描述电力系统的非线性方程。这里仍从节点电压方程入手，设电力系统导纳矩阵已知，则系统中某节点（i 节点）电压方程为 从而得 ∑=* *? ?=n j j ij i i U Y U S 1 进而有 ()01 =-+* =* ? ∑j n j ij i i i U Y U jQ P （1 －6）

PQ分解法潮流计算实验

xxxx实验报告 学生姓名：学号：专业班级： 实验类型：□验证□综合■设计□创新实验日期：实验成绩： 一、实验项目名称 P-Q分解法潮流计算实验 二、实验目的与要求： 目的：电力系统分析的潮流计算是电力系统分析的一个重要的部分。通过对电力系统潮流分布的分析和计算，可进一步对系统运行的安全性，经济性进行分析、评估，提出改进措施。电力系统潮流的计算和分析是电力系统运行和规划工作的基础。 潮流计算是指对电力系统正常运行状况的分析和计算。通常需要已知系统参数和条件，给定一些初始条件，从而计算出系统运行的电压和功率等；潮流计算方法很多：高斯-塞德尔法、牛顿-拉夫逊法、P-Q分解法、直流潮流法，以及由高斯-塞德尔法、牛顿-拉夫逊法演变的各种潮流计算方法。 本实验采用P-Q分解法进行电力系统分析的潮流计算程序的编制与调试，获得电力系统中各节点电压，为进一步进行电力系统分析作准备。通过实验教学加深学生对电力系统潮流计算原理的理解和计算，初步学会运用计算机知识解决电力系统的问题，掌握潮流计算的过程及其特点。熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。 要求：编制调试电力系统潮流计算的计算机程序。程序要求根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，完成该电力系统的潮流计算，要求计算出节点电压、功率等参数。 三、主要仪器设备及耗材 每组计算机1台、相关计算软件1套 四、实验内容： 1.理论分析： P-Q分解法是从改进和简化牛顿法潮流程序的基础上提出来的，它的基本思想是：把节点功率表示为电压向量的极坐标方程式，抓住主要矛盾，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功功率和无功功率迭代分开来进行。 牛顿法潮流程序的核心是求解修正方程式，当节点功率方程式采取极坐标系统时，

潮流计算简答题

潮流计算数学模型与数值方法 1. 什么是潮流计算？潮流计算的主要作用有哪些？ 潮流计算，电力学名词，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。 潮流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求；对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(线路、变压器等)是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。 2. 潮流计算有哪些待求量、已知量？ （已知量：1、电力系统网络结构、参数 2、决定系统运行状态的边界条件 待求量：系统稳态运行状态 例如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等） 3. 潮流计算节点分成哪几类？分类根据是什么？ （分成三类：PQ 节点、PV 节点和平衡节点，分类依据是给定变量的不同） 4. 教材牛顿-拉夫逊法及有功-无功分解法是基于何种电路方程？可否采用其它类型方程？ 答:基于节点电压方程，还可以采用回路电流方程和割集电压方程等。但是后两者不常用。 5. 教材牛顿-拉夫逊法是基于节点阻抗方程、还是基于节点导纳方程进行迭代计算的？试阐述这两种方程的优点与缺点。 1.不能由等值电路直接求出 2.满秩矩阵内存量大 3.对角占优矩阵。。 节点导纳矩阵的特点：1.直观容易形成2.对称阵3.稀疏矩阵（零元素多）：每一行的零元素个数=该节点直接连出的支路数。 6. 说出至少两种建立节点导纳矩阵的方法，阐述其中一种方法的原理与过程。 方法：1.根据自导纳和互导纳的定义直接求取2.运用一节点关联矩阵计算3.阻抗矩阵的逆矩阵 节点导纳矩阵的形成：1.对角线元素ii Y 的求解)1,,0(=≠==i j I i ii U i j U U I Y 【除i 外的其他节点接地，0=j U ，只在i 节点加单位电压值】解析ii Y 等于与i 节点直接相连的的所有支路导纳和2.互导纳),0,1(j k U U U I Y k j j i ij ≠===，ji ij Y Y =（无源网络导纳之间是对称的）解析：ij Y 等于j i ,节点之间直接相连的支路导纳的负值。 7. 潮流计算需要考虑哪些约束条件？ 答: 为了保证系统的正常运行必须满足以下的约束条件：

matlab潮流计算

附录1 使用牛顿拉夫逊法进行潮流计算的Matlab程序代码 % 牛拉法计算潮流程序 %----------------------------------------------------------------------- % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量% 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；3为PV节点； %------------------------------------------------------------------------ clear all; format long; n=input('请输入节点数:nodes='); nl=input('请输入支路数:lines='); isb=input('请输入平衡母线节点号:balance=');

pr=input('请输入误差精度:precision='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); %------------------------------------------------------------------ for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); %非对角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4); %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4); %对角元1侧 end %求导纳矩阵