电力系统潮流计算仿真实验报告
电力仿真算法实验报告
一、实验目的本次实验旨在通过电力系统仿真软件对电力系统进行仿真分析，验证电力系统仿真算法的有效性，并进一步了解电力系统在不同运行条件下的稳定性和性能。
实验内容包括电力系统潮流计算、暂态稳定分析、短路电流计算等。
二、实验内容1. 电力系统潮流计算（1）实验背景：以某地区110kV电网为例，分析该电网在不同运行方式下的潮流分布。
（2）实验步骤：① 利用电力系统仿真软件建立110kV电网模型；② 设置电网参数，包括各节点电压、线路参数等；③ 运行潮流计算程序，得到潮流分布结果；④ 分析潮流分布结果，判断电网的稳定性。
2. 电力系统暂态稳定分析（1）实验背景：以某地区110kV电网为例，分析该电网在发生单相接地故障时的暂态稳定性。
（2）实验步骤：① 利用电力系统仿真软件建立110kV电网模型；② 设置电网参数，包括各节点电压、线路参数等；③ 设置故障参数，包括故障类型、故障位置等；④ 运行暂态稳定分析程序，得到暂态稳定结果；⑤ 分析暂态稳定结果，判断电网的稳定性。
3. 电力系统短路电流计算（1）实验背景：以某地区110kV电网为例，计算电网在发生短路故障时的短路电流。
（2）实验步骤：① 利用电力系统仿真软件建立110kV电网模型；② 设置电网参数，包括各节点电压、线路参数等；③ 设置故障参数，包括故障类型、故障位置等；④ 运行短路电流计算程序，得到短路电流结果；⑤ 分析短路电流结果，判断电网的短路容量。
三、实验结果与分析1. 电力系统潮流计算结果通过潮流计算，得到110kV电网在不同运行方式下的潮流分布。
结果表明，在正常运行方式下，电网的潮流分布合理，节点电压满足要求。
在故障运行方式下，电网的潮流分布发生较大变化，部分节点电压超出了允许范围。
2. 电力系统暂态稳定分析结果通过暂态稳定分析，得到110kV电网在发生单相接地故障时的暂态稳定结果。
结果表明，在故障发生初期，电网暂态稳定，但故障持续一段时间后，电网发生暂态失稳。
仿真潮流实验报告
一、实验目的1. 理解电力系统潮流计算的基本原理和方法。
2. 掌握MATLAB/Simulink在电力系统仿真中的应用。
3. 通过仿真实验，验证潮流计算的正确性和实用性。
二、实验原理与内容1. 潮流计算的基本原理潮流计算是电力系统分析的重要手段，用于计算电力系统各节点的电压、相角、功率等参数。
其基本原理如下：（1）根据电力系统的网络结构和参数，建立节点方程和支路方程。
（2）利用节点方程和支路方程，求解节点电压和相角。
（3）根据节点电压和相角，计算各节点的有功功率和无功功率。
2. 仿真实验内容本次仿真实验采用MATLAB/Simulink搭建一个简单的2机5节点电力系统模型，并利用PowerGUI进行潮流计算。
（1）建立电力系统模型首先，在MATLAB/Simulink中搭建电力系统模型，包括发电机、负荷、线路等元件。
根据实验要求，设置发电机参数、负荷参数和线路参数。
（2）潮流计算利用PowerGUI进行潮流计算，设置求解器参数，如迭代次数、收敛精度等。
运行潮流计算，得到各节点的电压、相角、有功功率和无功功率等参数。
（3）结果分析对潮流计算结果进行分析，验证潮流计算的正确性和实用性。
比较不同运行方式下的潮流计算结果，分析系统稳定性。
三、实验方法1. 利用MATLAB/Simulink搭建电力系统模型。
2. 利用PowerGUI进行潮流计算。
3. 对潮流计算结果进行分析。
四、实验步骤1. 启动MATLAB/Simulink，新建一个仿真模型。
2. 在仿真模型中，添加发电机、负荷、线路等元件，设置相应参数。
3. 将搭建好的电力系统模型连接起来，形成一个完整的系统。
4. 打开PowerGUI，选择潮流计算模块。
5. 在潮流计算模块中，设置求解器参数，如迭代次数、收敛精度等。
6. 运行潮流计算，得到各节点的电压、相角、有功功率和无功功率等参数。
7. 对潮流计算结果进行分析，验证潮流计算的正确性和实用性。
五、实验结果与分析1. 潮流计算结果本次仿真实验中，潮流计算结果如下：（1）节点电压：U1=1.02p.u., U2=1.05p.u., U3=1.03p.u., U4=1.00p.u., U5=1.01p.u.（2）节点相角：δ1=0.5°, δ2=1.0°, δ3=0.7°, δ4=0.0°, δ5=0.6°（3）有功功率：P1=100MW, P2=100MW, P3=100MW, P4=100MW, P5=100MW（4）无功功率：Q1=20Mvar, Q2=20Mvar, Q3=20Mvar, Q4=20Mvar, Q5=20Mvar2. 结果分析（1）节点电压和相角在合理范围内，说明潮流计算正确。
武汉大学电力系统分析仿真实验报告
PSASP电力系统分析仿真实验报告姓名：学号：实验二基于PSASP的电力系统潮流计算实验一、实验目的掌握用PSASP进行电力系统潮流计算方法。
二、实验内容以上为系统常规运行方式的单线图。
由于母线STNB-230处负荷的增加，需对原有电网进行改造，具体方法为：在母线GEN3-230和STNB-230之间增加一回输电线，增加发电3的出力及其出口变压器的容量，新增或改造的元件如下图虚线所示：母线数据:交流线数据其基础数据如下：数 I 测 据 母线 J 测母 线 编 所 单 正序 号 属 位 阻抗正序 电抗正序 充电 零序 电阻 零序 电抗零序充电 电纳的 组区电纳1/2变压器数据续上表发电数据续上表负荷数据区域定义数据方案定义潮流计算作业定义三、实验步骤（1）点击|“电力系统分析综合程序（PSASP）”；（2）点击“创建”，创建文件；（3）点击“图形支持环境”；（4）点击“编辑模式”，可进行绘图和参数录入：a、绘制出所有母线，输入母线数据；b、添加发电机、负荷、交流线、变压器、支路，输入该元件数据；（5）关闭“编辑模式”窗口；（6）点击“运行模式”：（7）点击“作业”菜单项，执行“方案定义”命令（例如方案为1，数据组选择BASIC），点击“确定”。
（8）点击“作业”菜单项，执行“潮流”命令，定义作业；（9）点击“视图”菜单项，执行“潮流数据”命令，作业选择。
（10）点击“计算”菜单项，执行“潮流”命令；（11）点击“格式”菜单项，进行元件参数格式选择；（12）点击“报表”菜单项，执行“潮流”命令,计算结果输出有图示、报表输出两种方。
四、实验注意事项（1）本系统文件请存入D：\PSASP\学号\潮流\下；电cosθg总有功负荷cosθl总有功总无功损（2）严禁删除或更改计算机中除上述目录以外的一切内容。
五、实验结果（1）将实验结果采用图示、报表列出。
（2）需要完成的计算任务:1、查看作业号1的计算结果潮流计算摘要信息报表PSASP(Load Flow)EPRI,China计算日期：2014/04/02时间：15:47:50作业号：1作业描述：计算方法：Newton(Power Equation)基准容量：100.0000(MW)允许误差：0.000100本系统上限母线：910000发电机：33000负荷：36000交流线：610000直流线：010两绕组变压器：37000三绕组变压器：02000移相变压器：0200UD模型调用次数：0200UP调用个数：010结果综述报表作业号：1计算日期：2014/04/02时间：15:47:50单位：p.u.区域名区域号总有功发电总无功发荷总无功负损耗耗区域-11 2.48-0.042060.9998610.350.943860.01917-0.44282区域20.716410.270460.93555 2.150.80.937220.02724-0.47878号I侧有功I侧无功I侧充电功率J侧有功J侧无-2全网 3.196410.22840.99746 3.15 1.150.939360.04641-0.9216全网母线(发电、负荷)结果报表作业号：1计算日期：2014/04/02时间：15:47:50单位：p.u.母线名电压幅值电压相角GEN1-230 1.02579-2.2168GEN2-230 1.02577 3.7197GEN3-230 1.03235 1.9667STNA-2300.99563-3.9888STNB-230 1.01265-3.6874STNC-230 1.015880.7275发电1 1.040发电2 1.0259.28发电3 1.025 4.6648全网交流线结果报表作业号：1计算日期：2014/04/02时间：15:47:50单位：p.u.I侧母线名J侧母线名编功J侧充电功率GEN1-230STNA-23010.409370.228930.09260.40680.386870.08723 GEN2-230STNC-23030.7638-0.007970.078390.759050.107040.07689 GEN3-230STNB-23050.60817-0.180750.190770.594630.134570.18356 STNA-230GEN2-2302-0.8432-0.113130.15167-0.86620.083810.16099 STNB-230GEN1-2306-0.30537-0.165430.08101-0.30704-0.01030.08313 STNC-230GEN3-2304-0.24095-0.242960.10785-0.24183-0.03120.11137全网两绕组变压器结果报表作业号：1计算日期：2014/04/02时间：15:47:50单位：p.u.I侧母线名J侧母线名编号I侧有功I侧无功J侧有功J侧无功发电1GEN1-23070.716410.270460.716410.23923发电2GEN2-2308 1.630.06654 1.63-0.09178发电3GEN3-23090.85-0.10860.85-0.14955电cosθg总有功负荷cosθl总有功总无功损改用PQ分解法重复计算作业1，查看计算结果，与牛顿法结果做比较。
电力系统分析仿真实验报告
电力系统分析仿真实验报告****名目实验一电力系统分析综合程序PSASP概述一、实验目的了解用PSASP进行电力系统各种计算的方法。
二、PSASP简介1．PSASP是一套功能强大，使用方便的电力系统分析综合程序，是具有我国自主知识产权的大型软件包。
2．PSASP的体系结构：第一层是：公用数据和模型资源库，第二层是应用程序包，第三层是计算结果和分析工具。
3．PSASP的使用方法：〔以短路计算为例〕1).输进电网数据，形成电网根底数据库及元件公用参数数据库，〔后者含励磁调节器，调速器，PSS等的固定模型〕，也可使用用户自定义模型UD。
在此，可将数据合理组织成假设干数据组，以便下一步形成不同的计算方案。
✧文本支持环境：点击“数据〞菜单项，执行“根底数据〞和“公用参数〞命令，可依次输进各电网元件的参数。
✧图形支持环境：在“编辑模式下〞，利用工具箱，输进电网接线图。
作图时，假设元件参数尚未输进，会自动弹出相关数据录进窗口，如今输进数据即可。
注重：两种环境下，均应先输进母线数据，再处理其他元件！！！2).方案定义：从根底数据库中抽取数据组，组合成不同方案，以确定电网的规模，结构和运行方式。
✧文本支持环境：点击“计算〞菜单项，执行“方案定义〞命令。
✧图形支持环境：“运行模式〞下，点击“作业〞菜单项，执行“方案定义〞命令。
3〕数据检查：对确定的电网结构进行检查，检查网架结构的合理性，计算规模是否超出范围。
✧文本支持环境：点击“计算〞菜单项，执行“数据检查〞命令。
✧图形支持环境：“运行模式〞下，点击“作业〞菜单项，执行“数据检查〞命令。
4〕作业定义：给出计算操纵信息，明确具体的计算任务。
✧文本支持环境：点击“计算〞菜单项，执行“短路〞命令。
✧图形支持环境：“运行模式〞下，点击“作业〞菜单项，执行“短路〞命令。
5〕执行计算：✧文本支持环境：在上述“短路计算信息〞窗口，完成作业定义之后，点击“计算〞按钮即可。
✧图形支持环境：“运行模式〞下，a.点击“视图〞菜单项，执行“短路〞命令，选择作业；b.点击“计算〞菜单项，执行“短路〞命令，执行计算；c.点击“格式〞菜单项，执行“短路结果〞命令，确定计算结果在图上的显示方式。
电力系统仿真实训报告电力系统仿真实训
电力系统仿真实训报告1 前言电力系统潮流计算是研究电力系统稳态运行情况的一种计算,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态。
在电力系统规划设计和现有电力系统运行方式的研究中,都需要利用潮流计算来定量分析、比较供电方案或运行方式的合理性、可靠性和经济性。
本次课程设计任务是闭环网络的潮流计算，用到的方法为PQ分解法潮流计算。
2 实训目的与要求2.1实训目的电力系统分析的潮流计算是电力系统分析的一个重要的部分。
通过对电力系统潮流分布的分析和计算，可进一步对系统运行的安全性，经济性进行分析、评估，提出改进措施。
电力系统潮流的计算和分析是电力系统运行和规划工作的基础。
潮流计算是指对电力系统正常运行状况的分析和计算。
通常需要已知系统参数和条件，给定一些初始条件，从而计算出系统运行的电压和功率等；潮流计算方法很多：高斯-塞德尔法、牛顿-拉夫逊法、P-Q分解法、直流潮流法，以及由高斯-塞德尔法、牛顿-拉夫逊法演变的各种潮流计算方法。
本实验采用P-Q分解法进行电力系统分析的潮流计算程序的编制与调试，获得电力系统中各节点电压，为进一步进行电力系统分析作准备。
通过实验教学加深学生对电力系统潮流计算原理的理解和计算，初步学会运用计算机知识解决电力系统的问题，掌握潮流计算的过程及其特点。
熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。
2.2实训要求编制调试电力系统潮流计算的计算机程序。
程序要求根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，完成该电力系统的潮流计算，要求计算出节点电压、功率等参数。
3 实训内容1 基于PSASP的电力系统潮流计算仿真1.1 实验要求要求在掌握电力系统稳态分析知识的基础上，根据PSASP中电力系统潮流计算的步骤，利用该软件实现电力系统的潮流计算，并能根据潮流计算结果，对电力系统进行运行情况分析。
电力系统潮流计算数字仿真实验实验报告（PASAP软件的使用）
电力系统潮流计算数字仿真实验实验报告（PASAP软件的使用）电力系统潮流计算数字仿真实验实验报告班级姓名学号五节点潮流计算（原图） (3)不同负荷水平时刻的潮流计算——负荷加倍 (6)各节点负荷减半作为谷时 (8)两种不同运行方式的潮流计算 (10)两种不同调度结果的潮流计算 (12)2一、五节点潮流计算（原图）结果导出：①物理母线单位：kA\kV\MW\Mvar区域分区厂站全网全网全网母线名称电压幅值电压相角-------- -------- --------* b 103.48250 -8.0217g 217.28740 -5.1551x 106.31940 -10.2418Ⅰ 242.00000 0Ⅱ 116.39210 -2.35793②发电机单位：kA\kV\MW\Mvar区域分区厂站全网全网全网发电机名称母线名类型有功发电无功发电功率因数---------- ------ ---- -------- -------- -------- Gen_1 Ⅰ Vθ 202.4800 155.9800 0.79220 Gen_2 Ⅱ PQ 40.0000 30.0000 0.80000 Gen_3 g PQ0.0000 10.0000 0.00000③负荷结果报表单位：kA\kV\MW\Mvar区域分区厂站全网全网全网负荷名称母线名类型有功负荷无功负荷功率因数-------- ------ ---- -------- -------- -------- Load_1 b PQ 50.0000 30.0000 0.85749 Load_3 x PQ 180.0000 100.00000.87416④交流线结果报表单位：kA\kV\MW\Mvar区域分区全网全网交流线名称 I侧母线 J侧母线 I侧电压 I侧有功 I侧无功---------- ------- -------- ------- ------- ------- AC_1 Ⅰ g 242.00000 180.4800 162.8100 AC_3 Ⅱ b 116.39210 61.9700 22.1800 AC_9 bx 103.48250 6.7700 -15.8100 J侧电压 J侧有功 J侧无功------- ------- -------217.28740 174.5300 131.0300103.48250 56.7700 14.1900106.31940 6.3200 -16.50004⑤两绕组变压器结果报表单位：kA\kV\MW\Mvar区域分区厂站全网全网全网两绕组变压器名称 I侧母线 J侧母线 I侧有功---------------- ------- ------- -------T2w_1 ⅠⅡ 21.9900 T2w_2 g x 174.5300I侧无功 J侧有功 J侧无功------- ------- --------6.8200 21.9600 -7.8200141.0300 173.6800 116.50005二、不同负荷水平时刻的潮流计算——负荷加倍问题：如果不增加无功功率，负荷在正常运行基础上加倍时潮流计算不成功，因为根据电压损耗公式，负荷加倍后电压损耗变大，部分母线电压过低，所以潮流计算不成功。
电力系统分析仿真实验报告
电力系统分析仿真实验报告电力系统分析仿真实验报告一、实验目的本次实验的目的是通过使用电力系统仿真软件进行电力系统模拟分析，掌握电力系统运行特点及原因、掌握电力系统基本传输线路的参数，以及了解电力系统的潮流分布计算和短路分析流程。
二、实验原理电力系统仿真软件是针对电力系统运行及其各种故障情况下的仿真软件。
仿真软件将电力系统进行模拟分析，可以让使用者对电力系统进行检测修正，达到保证电网质量的目的。
仿真软件主要采用数学模型进行计算，本次实验中使用的仿真软件为PSASP。
第一，电力系统线路模型：电力系统线路模型是电力系统分析的核心内容，此模型通过研究电力线路本身的运行特点，为电力系统计算和仿真打下基础。
电力系统线路模型假设电力系统线路为非常规线路，包括零序电感和阻抗、平衡、非平衡模型的相间电路等，具体包括电感、电容、电阻三部分。
第二，电力系统模拟分析：电力系统的仿真分析，就是对电力系统进行计算、仿真，从而得出电力系统的各种参数或特性。
模拟分析主要包括电力系统的潮流计算、电力系统的短路分析等两个方面。
（1）电力系统潮流计算：电力系统潮流计算是指通过对电力系统进行数学建模，来分析电力系统中电流、电压等各种状态量的分布规律。
具体的计算过程采用功率系统仿真软件进行计算。
（2）电力系统短路分析：电力系统短路分析是针对电力系统在遭受外部灾害时计算其在各种短路状态下的可能损伤程度，在电力系统建设过程中非常重要。
同时也是保障电网电力质量安全的必要手段。
三、实验内容实验的主要内容分为两个部分，第一部分是电力系统潮流计算实验，第二部分是电力系统短路分析实验。
（1）潮流计算实验这部分实验的主要内容是计算电力系统的电流分布以及电压分布等参数，实验过程如下：1. 打开PSASP软件，新建项目档案。
根据实际需求设置主进程，建立相应关系文件，并完成电力系统初始化操作。
2. 添加仿真数据。
根据实验要求，添加相应的电力系统数据。
其中包括节点数据、主变和传输线路数据、变压器等数据。
电力系统计算机辅助潮流计算 实验报告(DOC)
（二 〇 一 一 年 十 二 月（2011-2012学年第一学期）电力系统计算机辅助潮流计算实验报告学 号： ************学生姓名：学 院：电力学院 系 别：电力系专 业：电气工程及其自动化 班 级： 授课教师：1、实验目的：了解计算机潮流分析的基本原理、主要步骤；掌握节点导纳矩阵形成和修改的方法，掌握数据处理的基本方法；熟悉Matlab运行环境，了解Matlab基本编程语句和语法；运用潮流分析程序对给定网络的运行方式做潮流分析，并初步分析计算结果2、实验要求：通过预习，对计算机潮流分析基本理论有深入了解；为程序准备必要的、准确的原始数据；熟悉Matlab运行环境，输入潮流程序，上机独立完成程序的调试，给出潮流分析的结果并按要求绘制潮流分布图3、实验内容：输入网络参数，包括节点号、节点导纳矩阵、节点功率等；输入潮流程序、调试并输出结果，绘制潮流分布图4、实验步骤：1、熟悉原始资料：根据计算要求，整理数据，包括：计算网络中线路、变压器的参数、形成节点导纳矩阵；表示各节点的注入功率。
（以上数据均采用有名值计算）2、读通潮流程序：完成程序的解释和说明，必要时附加对应的公式和程序语言的说明3、上机调试：熟悉Matlab的运行环境，准确输入原始数据、节点编号、节点注入功率等信息4、整理计算结果：根据计算结果作电网潮流分布图原始网络：5、实验数据及处理：一、实验程序：clearG(1,1)=3.75;B(1,1)=-11.25;G(1,2)=-2.5;B(1,2)=7.5;G(1,3)=0;B(1,3)=0; G(1,4)=-1.25;B(1,4)=3.75;G(1,5)=0;B(1,5)=0;G(2,1)=-2.5;B(2,1)=7.5;G(2,2)=10.834;B(2,2)=-32.5;G(2,3)=-1.667;B(2,3)=5;G(2,4)=-1.667;B(2,4)=5;G(2,5)=-5;B(2,5)=15;G(3,1)=0;B(3,1)=0;G(3,2)=-1.667;B(3,2)=5;G(3,3)=12.917;B(3,3)=-38.75;G(3,4)=-10;B(3,4)=30;G(3,5)=-1.25;B(3,5)=3.75;G(4,1)=-1.25;B(4,1)=3.75;G(4,2)=-1.667;B(4,2)=5;G(4,3)=-10;B(4,3)=30;G(4,4)=12.917;B(4,4)=-38.750;G(4,5)=0;B(4,5)=0;G(5,1)=0;B(5,1)=0;G(5,2)=-5;B(5,2)=15;G(5,3)=-1.25;B(5,3)=3.75;G(5,4)=0;B(5,4)=0;G(5,5)=6.25;B(5,5)=-18.75;Y=G+j*B %形成节点导纳矩阵delt(1)=0;delt(2)=0;delt(3)=0; delt(4)=0;u(1)=1.0;u(2)=1.0;u(3)=1.0;u(4)=1.0;ps(1)=-0.6;qs(1)=-0.10;ps(2)=0.2;qs(2)=0.2;ps(3)=-0.45;qs(3)=-0.15; ps(4)=-0.4;qs(4)=-0.05; %设迭代初值k=1;precision=1 %设迭代次数和精度N1=4; %PQ节点数while precision>0.0000001 %判断是否满足精度要求delt(5)=0;u(5)=1.06; %给定平衡节点编号for m=1:N1for n=1:N1+1pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt (n)));qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt (n)));endpi(m)=sum(pt);qi(m)=sum(qt); %计算PQ节点的注入功率 dp(m)=ps(m)-pi(m);dq(m)=qs(m)-qi(m); %计算PQ节点的功率不平衡量endfor m=1:N1for n=1:N1if m==nH(m,m)=-qi(m)-u(m)^2*B(m,m); N(m,m)=pi(m)+u(m)^2*G(m,m);J(m,m)=pi(m)-u(m)^2*G(m,m); L(m,m)=qi(m)-u(m)^2*B(m,m);JJ(2*m-1,2*m-1)=H(m,m); JJ(2*m-1,2*m)=N(m,m);JJ(2*m,2*m-1)=J(m,m); JJ(2*m,2*m)=L(m,m);elseH(m,n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-del t(n)));J(m,n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-de lt(n)));N(m,n)=-J(m,n);L(m,n)=H(m,n);JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=N(m,n);JJ(2*m,2*n-1)=J(m,n); JJ(2*m,2*n)=L(m,n);Endendend %计算jocbi各项,并放入统一矩阵JJ中，对JJ下标统一编号JJfor m=1:N1PP(2*m-1)=dp(m);PP(2*m)=dq(m);End %按统一矩阵形成功率不平衡 uu=inv(JJ)*PP';precision=max(abs(uu)); %判断是否收敛 for n=1:N1delt(n)=delt(n)+uu(2*n-1);u(n)=u(n)+uu(2*n)*u(n); %将结果分解为电压幅值和角度 end %求解修正方程，得电压幅值变化量（标幺值）和角度变化量 k=k+1; end for n=1:N1+1U(n)=u(n)*(cos(delt(n))+j*sin(delt(n))); endfor m=1:N1+1I(m)=Y(5,m)*U(m); %求平衡节点的注入电流5551j j i I Y U ==∑endS5=U(5)*sum(conj(I)) %求平衡节点的注入功率*555S V I = for m=1:N1+1 for n=1:N1+1S(m,n)=U(m)*(conj(U(m))-conj(U(n)))*conj(-Y(m,n));% end endend %求节点i,j 节点之间的功率，方向为由i 指向j, *ij ij i S V I = S %显示支路功率二、实验结果：1、节点导纳矩阵Y =3.7500 -11.2500i -2.5000 + 7.5000i 0 -1.2500 + 3.7500i 0 -2.5000 + 7.5000i 10.8340 -32.5000i -1.6670 + 5.0000i -1.6670 + 5.0000i -5.0000 +15.0000i 0 -1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i -1.2500 + 3.7500i -1.2500 + 3.7500i -1.6670 + 5.0000i -10.0000 +30.0000i 12.9170 -38.7500i 0 0 -5.0000 +15.0000i -1.2500 + 3.7500i 0 6.2500 -18.7500i节点导纳矩阵特点：1.节点导纳矩阵的对角元就等于各该节点所连接导纳的总和2.节点导纳矩阵是稀疏矩阵3.节点导纳矩阵一般是对称矩阵4.节点导纳矩阵的非对角元Yij等于连接节点i、j支路导纳的负值5. 节点导纳矩阵是方阵2、迭代过程数据：电压变化量du = 0.0094 - 0.0010i 0.0448 - 0.0021i 0.0156 - 0.0013i 0.0142 - 0.0013idu =-0.0077 - 0.0000i -0.0064 - 0.0000i -0.0065 - 0.0000i -0.0068 -0.0000idu = -0.6060 - 0.0000i -0.4686 - 0.0000i -0.5040 - 0.0000i -0.5213 -0.0000idu = -0.3934 - 0.0000i -0.2741 + 0.0000i -0.3107 - 0.0000i -0.3238 -0.0000i功率不平衡量dS = -0.6000 - 0.1000i 0.5000 + 1.1000i -0.3750 + 0.0750i -0.4000 - 0.0500idS = 0.0208 - 0.0206i -0.0468 - 0.0876i 0.0047 - 0.0226i 0.0155 - 0.0096i dS= 1.0e-003 *0.2233 - 0.2193i -0.4990 - 0.5087i -0.0052 - 0.1658i 0.1556 - 0.0997i dS =1.0e-007 *0.2011 - 0.1958i -0.4104 - 0.2076i -0.0402 - 0.0900i 0.1273 - 0.0836i雅可比矩阵JJ =11.2500 3.7500 -7.5000 -2.5000 0 0 -3.7500 -1.2500-3.7500 11.2500 2.5000 -7.5000 0 0 1.2500 -3.7500-7.5000 -2.5000 33.4000 10.5340 -5.0000 -1.6670 -5.0000 -1.66702.5000 -7.5000 -11.1340 31.6000 1.6670 -5.0000 1.6670 -5.00000 0 -5.0000 -1.6670 38.9750 12.8420 -30.0000 -10.00000 0 1.6670 -5.0000 -12.9920 38.5250 10.0000 -30.0000-3.7500 -1.2500 -5.0000 -1.6670 -30.0000 -10.0000 38.7500 12.91701.2500 -3.7500 1.6670 -5.0000 10.0000 -30.0000 -12.9170 38.7500JJ =11.5406 3.1996 -7.7223 -3.1029 0 0 -3.8183 -1.3384-4.4412 11.3818 3.1029 -7.7223 0 0 1.3384 -3.8183-8.0396 -2.1511 35.0648 12.0317 -5.3599 -1.5576 -5.3622 -1.52382.1511 -8.0396 -11.5380 35.6400 1.5576 -5.3599 1.5238 -5.36220 0 -5.2222 -1.9705 40.0793 12.8630 -30.9519 -10.11360 0 1.9705 -5.2222 -13.7724 39.8246 10.1136 -30.9519-3.8576 -1.2203 -5.2038 -1.9989 -30.8297 -10.4802 39.8912 12.86851.2203 -3.8576 1.9989 -5.2038 10.4802 -30.8297 -13.6994 39.8104JJ =11.3861 3.1619 -7.6217 -3.0453 0 0 -3.7644 -1.3171-4.3623 11.1866 3.0453 -7.6217 0 0 1.3171 -3.7644-7.9246 -2.1368 34.7163 11.8401 -5.2904 -1.5440 -5.2913 -1.51162.1368 -7.9246 -11.4391 35.1173 1.5440 -5.2904 1.5116 -5.29130 0 -5.1585 -1.9397 39.5849 12.6953 -30.5425 -9.98760 0 1.9397 -5.1585 -13.5953 39.2852 9.9876 -30.5425-3.8018 -1.2050 -5.1397 -1.9661 -30.4266 -10.3354 39.3681 12.70621.2050 -3.8018 1.9661 -5.1397 10.3354 -30.4266 -13.5065 39.2683JJ =11.3850 3.1617 -7.6210 -3.0448 0 0 -3.7640 -1.3169-4.3617 11.1850 3.0448 -7.6210 0 0 1.3169 -3.7640-7.9237 -2.1368 34.7136 11.8386 -5.2899 -1.5439 -5.2907 -1.51152.1368 -7.9237 -11.4386 35.1136 1.5439 -5.2899 1.5115 -5.29070 0 -5.1580 -1.9394 39.5811 12.6940 -30.5394 -9.98670 0 1.9394 -5.1580 -13.5940 39.2811 9.9867 -30.5394-3.8013 -1.2049 -5.1392 -1.9659 -30.4235 -10.3343 39.3641 12.70501.2049 -3.8013 1.9659 -5.1392 10.3343 -30.4235 -13.5050 39.26413、收敛后数据：支路功率S =0 -0.5370 - 0.0977i 0 -0.0630 - 0.0023i 0 0.5489 + 0.1333i 0 0.2469 + 0.0815i 0.2793 + 0.0806i -0.8751 - 0.0954i 0 -0.2431 - 0.0701i 0 0.1891 - 0.0121i -0.3960 - 0.0677i 0.0633 + 0.0033i -0.2746 - 0.0664i -0.1887 + 0.0132i 0 0 0 0.8895 + 0.1387i 0.4087 + 0.1058i 0 0平衡节点功率 S5 =1.2982 + 0.2445i潮流分布图GG1.250-j3.750 10.000-j30.0002.500-j7.5005.000-j15.001.250-j3.7501.667-j5.0001.667-j5.000534210.45+j0.150.4+j0.050.2+j0.20.6+j0.10.4087+j0.1058 0.3960+j0.06770.1891-j0.0121 0.1887-j0.01320.5489+j0.150.5370+j0.09770.2431+j0.07010.2469+j0.08150.2746+j0.06640.2793+j0.0806 0.8895+j0.1387 0.8751+j0.09540.0633+0.00330.0630+j0.00231.2982+j0.2445。
电力系统仿真实验实验报告
电力系统仿真实验实验报告
一、实验目的
1. 通过电力系统仿真实验，了解电力系统组成及其分析的知识；
2. 掌握电力系统仿真的基本方法，并能熟练使用电力系统仿真软件；
3. 了解电力系统受外界因素引起的各种失稳的原因，熟悉失稳的分析和处理方法；
4. 了解电力系统异步机的调节原理，掌握调节器参数选择和设置的知识。
二、实验内容
1. 利用电力系统仿真软件，建立一个单机电力系统模型，熟悉软件操作。
2. 观察失稳现象，绘制各类型失稳的时域曲线和频率曲线，分析其原因，探讨失稳的处理方法。
3. 分析不同调速器的调速性能，并应用比较几种调节器的性能，比较几种调节器的参数优化方法。
4. 在单机电力系统中调整相位调节器，使系统负荷增加时能保持稳定。
三、实验结果
1. 通过仿真实验，我们建立了一个单机电力系统的模型，并熟悉了电力系统仿真软件的操作。
2. 观察到系统可能会出现的失稳现象，绘制了各类型失稳的时
域曲线和频率曲线，分析系统失稳的原因；采取合理的处理措施，改善系统的稳定性。
3. 对比了不同调速器的调速性能，探讨了几种调节器的参数优化方法；在单机电力系统中，调整了相位调节器，使系统负荷增加时能保持稳定。
四、总结
本次实验中，我们熟悉了电力系统的组成及其分析，学习了电力系统仿真的基本方法和电力系统软件的操作，分析了电力系统受外界因素所引起的失稳现象，研究了失稳的分析和处理方法。
对于调节器的调节原理和参数选择设置等也有了深入的了解。
通过本次实验，我们掌握了电力系统仿真技术的基本知识，为今后的研究和发展打下了良好的基础。
电力系统分析仿真实验报告
电力系统分析仿真实验报告一、实验目的本实验的目的是通过电力系统分析仿真来研究电力系统的稳态和暂态运行特性，并通过实验结果分析电力系统中存在的问题和改进方案。
二、实验原理1.电力系统稳态分析电力系统稳态分析是指在电力系统稳定运行条件下，对电力系统进行负荷流量和节点电压的计算和分析。
稳态分析的目的是确定电力系统的潮流分布、负荷特性和节点电压，从而评估系统的稳定性和能量传输效率。
2.电力系统暂态分析电力系统暂态分析是指在电力系统出现故障或突发负荷变化时，对系统暂时的电压、电流和功率进行计算和分析。
暂态分析的目的是研究系统在故障或负荷突变时的动态响应和稳定性，以便采取相应措施保障系统的安全稳定运行。
三、实验过程1.电力系统稳态分析实验（1）建立电力系统模型：根据实际情况，建立包含发电机、变电站、输电线路和负荷的电力系统模型。
（2）潮流计算：通过潮流计算方法，对电力系统的负荷流量、节点电压和功率分布进行计算。
（3）结果分析：分析潮流计算结果，评估系统的稳定性和能量传输效率，检查是否存在过负荷或电压偏差等问题。
2.电力系统暂态分析实验（1）建立电力系统模型：在稳态模型的基础上，引入系统故障或负荷突变事件，如短路故障、突发负荷增加等。
（2）暂态计算：通过暂态计算方法，对系统的电压、电流和功率在故障或负荷突变时的动态变化进行计算。
（3）结果分析：分析暂态计算结果，评估系统在故障或负荷突变时的动态响应和稳定性，检查是否存在电压暂降或过载等问题。
四、实验结果与分析1.电力系统稳态分析结果分析：根据潮流计算结果，评估系统的稳定性和能量传输效率，检查系统是否存在过负荷或电压偏差等问题。
如果存在问题，可以通过调整发电机发电功率、变压器变比或线路容量来改善系统运行状况。
2.电力系统暂态分析结果分析：根据暂态计算结果，评估系统在故障或负荷突变时的动态响应和稳定性，检查是否存在电压暂降或过载等问题。
如果存在问题，可以通过引入自动重启装置、电力调度系统等措施来提高系统的恢复能力和稳定性。
电力系统仿真计算报告
电力系统仿真计算报告
目录
一、潮流计算 (2)
1.1计算条件及基础数据 (2)
1.2 常规方式潮流运算 (4)
1.3 规划方式潮流运算 (5)
二、短路计算 (6)
2.1三相短路 (6)
2.2 单相接地 (9)
2.3 两相短路 (12)
2.4 两相接地短路 (14)
三、暂态稳定计算 (15)
3.1 基于常规方式 (16)
3.2 基于规划暂稳计算 (20)
一、潮流计算
1.1计算条件及基础数据
1.2 常规方式潮流运算
图1.1 常规潮流单线图
图1.2 常规潮流计算结果1.3 规划方式潮流运算
图1.3 规划潮流单线图
图1.4 规划潮流计算结果
二、短路计算
2.1三相短路
图2.1 三相短路计算条件
图2.2 三相短路单线图
图2.3 三相短路部分计算结果2.2 单相接地
图2.4 单相接地计算条件
图2.5 单相接地单线图
图 2.6 单相接地部分母线计算电压
2.3 两相短路
图2.6 两相短路计算条件
图2.7 两相短路单线图
图2.8 两相短路部分计算结果
2.4 两相接地短路
图2.9 两相接地计算条件
图2.10 两相接地单线图
图2. 11 两相接地部分计算结果三、暂态稳定计算
3.1 基于常规方式
图3.1常规暂稳计算条件
图3.2 常规暂稳单线图
图3.3 常规暂稳发电机功角
图3.4 常规暂稳部分母线电压
3.2 基于规划暂稳计算
图3.5 规划暂稳计算条件
图3.6 规划暂稳单线图
3.7 规划暂稳发电机功角
图3.8 规划暂稳部分母线电压。
潮流计算实验
PSASP 潮流计算一、实验目的理解电力系统分析中潮流计算的相关概念，掌握用PSASP 软件对系统潮流进行计算的过程。
学会在文本方式下和图形方式下的对潮流计算结果进行分析。
二、预习要求复习《电力系统分析》中有关潮流计算的内容，了解有关潮流计算的功能，理解常用潮流计算方法，了解PQ、PV和Vθ（平衡节点，在PSASP中称为Slack节点）的设置。
三、实验内容（一） PSASP 潮流计算概述潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。
通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。
待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。
PSASP 潮流计算的流程和结构如下图所示：潮流计算各种计算公共部分文本方式图形方式结果的编辑和输出计算结果库执行计算图形方式文本方式计算作业的定义(运行方式和计算控制)方案定义(电网结构)文本方式图形方式用户自定义模型库电网基础数据库图形方式文本方式数据录入和编辑以一个图所示9节点系统为例，计算其在常规、规划两种运行方式下的潮流。
规划运行方式即在常规运行方式下，其中接于一条母线（STNB-230）处的负荷增加，对原有电网进行改造后的运行方式，具体方法为：在母线GEN3-230和STNB-230之间增加一回输电线，增加发电机3的出力及其出口变压器的容量，新增或改造的元件如下图虚线所示。
（二） 数据准备1. 指定数据目录及基准容量双击PSASP图标，弹出PSASP封面后，按任意键，即进入PSASP主画面：在该画面中，要完成的工作如下：(1) 指定数据目录第一次可通过“创建数据目录” 按钮，建立新目录，如：C:\CLJS。
以后可通过“选择数据目录”按钮，选择该目录。
(2) 给定系统基准容量系统基准容量项中，键入该系统基准容量，如100MVA。
电力系统分析综合实验一：潮流计算实验
电力系统分析综合实验一：潮流计算实验课程名称：电力系统分综合实验指导老师：成绩：实验名称：潮流分析实验实验类型：一、实验目的和要求（必填）二、实验内容和原理（必填）三、主要仪器设备（必填）四、操作方法和实验步骤五、实验数据记录和处理六、实验结果与分析（必填）七、讨论、心得一、实验目的和要求通过本实验熟悉PSCAD仿真软件，了解电力系统有功、无功的概念，掌握简单的潮流计算二、原理与说明（1）PSCAD/EMTDC是加拿大马尼托巴高压直流研究中心（ManitobaHVDCResearchCenter）推出的一款电力系统电磁暂态仿真软件。
EMTDC （ElectroMagneticTransientinDCSystem）是一种电力系统仿真分析软件，既可以研究交直流电力系统问题，又能完成电力电子仿真及其非线性控制的多功能工具。
PSCAD（PowerSystemComputerAidedDesign）是EMTDC的用户图形界面，也是EMTDC的图形用户接口，它的存在是为了更方便用户使用EMTDC。
PSCAD的开发成功，使得用户能更方便地使用EMTDC进行电力系统分析，使电力系统复杂部分可视化成为可能，而且软件可以作为实时数字仿真器的前置端。
用PSCAD进行潮流计算的时候要注意：①选取元件的时候，须辨别元件的属性，例如电源是三相还是单相；②设置元件参数时，要注意参数的物理意义③注意测量信号的方向④将多个曲线绘制在一张图表上时，须选取合适的范围，将单位标注清楚。
（2）三节点单相实验系统介绍在PSCAD界面搭建一个系统为50Hz的单相系统，如下图所示参数如下：G1：理想电压电源，100kVRMS，相角为0，频率为50HzG2：理想电压电源，100kVRMS，相角为60，频率为50HzG3：理想电压电源，100kVRMS，相角为30，频率为50HzX1：频率为50Hz的条件下，感抗为2X2：频率为50Hz的条件下，感抗为2X3：频率为50Hz的条件下，感抗为3X4：频率为50Hz的条件下，感抗为3X5：频率为50Hz的条件下，感抗为4X6：频率为50Hz的条件下，感抗为4Xline1：频率为50Hz的条件下，感抗为40Xline2：频率为50Hz的条件下，感抗为20Xline3：频率为50Hz的条件下，感抗为20仿真设置：电源：电压源输出坡至1.0p.u.的时间：0.05s仿真过程总时长：1s （确保仿真过程达到稳态）仿真步长：=50s（3）并联电容器并联电容器又称为移相电容器，是电力系统中一种重要的无功功率补偿设备，广泛地应用于改善负荷的功率因数。
电力系统-潮流计算仿真报告
Beijing Jiaotong University电力系统潮流计算仿真报告…姓名：TYP班级：电气0906学号：指导老师：吴俊勇完成日期：一、实验内容电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。
它的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。
电力系统潮流计算的结果是电力系统稳定计算和故障分析的基础。
对于简单系统，可以将其分为开式网络和闭式网络手工计算。
对于复杂电力系统，根据定解条件，应用牛顿—拉夫逊法进行计算，在手工计算中，由于涉及大量变量、微分方程、矩阵计算，求解很烦琐，而且容易出错，计算不同系统时需要重新计算。
故而我们可以借助计算机来进行潮流计算，方便快捷且准确率高。
二、计算机潮流计算方法@我们常用牛顿—拉夫逊法来进行潮流计算。
牛顿—拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方法，其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程，即通常所称的逐次线性化过程。
1、基本原理从几何意义上，牛顿—拉夫逊法实质上就是切线法，是一种逐步线性化的方法。
2、牛顿—拉夫逊法潮流求解过程以下讨论的是用直角坐标形式的牛顿—拉夫逊法潮流的求解过程。
当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量，由于平衡节点的电压向量是给定的，因此待求量共2(n-1)需要2(n-1)个方程式。
事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。
求解过程大致可以分为以下步骤：！(1)形成节点导纳矩阵；(2)将各节点电压设初值；(3)将节点初值代入相关求式，求出修正方程式的常数项向量；(4)将节点电压初值代入求式，求出雅可比矩阵元素；(5)求解修正方程，求修正向量；(6)求取节点电压的新值；(7)检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步；(8)计算支路功率分布，PV节点无功功率和平衡节点功率。
电力网潮流计算实训报告
一、实训目的本次电力网潮流计算实训旨在使学生掌握电力系统潮流计算的基本原理和方法，提高学生对电力系统运行状态的分析能力，培养学生运用所学知识解决实际问题的能力。
通过实训，使学生能够熟练运用MATLAB等工具进行电力系统潮流计算，并能够根据计算结果对电力系统运行状态进行分析和评价。
二、实训内容1. 电力系统潮流计算的基本原理电力系统潮流计算是一种分析电力系统稳态运行状态的方法，其主要任务是确定电力系统中各节点的电压、电流和功率分布。
电力系统潮流计算的基本原理如下：（1）节点电压方程：根据节点电压和注入功率之间的关系，建立节点电压方程。
（2）支路功率方程：根据支路两端电压和电流之间的关系，建立支路功率方程。
（3）潮流计算方法：采用迭代方法求解节点电压方程和支路功率方程，得到电力系统中各节点的电压、电流和功率分布。
2. 电力系统潮流计算实训本次实训以MATLAB为工具，进行电力系统潮流计算。
实训步骤如下：（1）建立电力系统模型：根据实训要求，建立电力系统模型，包括节点、支路、变压器等。
（2）设置电力系统参数：根据电力系统模型，设置各节点的注入功率、支路参数、变压器参数等。
（3）进行潮流计算：运用MATLAB中的电力系统分析工具箱（Power System Toolbox）进行潮流计算。
（4）分析计算结果：对潮流计算结果进行分析，包括节点电压、电流、功率分布、网络损耗等。
三、实训结果与分析1. 潮流计算结果通过潮流计算，得到以下结果：（1）节点电压：各节点电压的幅值和相角。
（2）电流：各支路电流的幅值和相角。
（3）功率分布：各支路的有功功率和无功功率。
（4）网络损耗：电力系统中的网络损耗。
2. 结果分析（1）节点电压分析：分析各节点电压是否满足要求，是否存在过电压或欠电压现象。
（2）电流分析：分析各支路电流是否满足要求，是否存在过载现象。
（3）功率分布分析：分析电力系统中功率分布是否合理，是否存在不平衡现象。
武汉大学电力系统分析仿真实验报告
PSASP电力系统分析仿真实验报告姓名：学号：实验二基于PSASP的电力系统潮流计算实验一、实验目的掌握用PSASP进行电力系统潮流计算方法。
二、实验内容以上为系统常规运行方式的单线图。
由于母线 STNB-230 处负荷的增加，需对原有电网进行改造，具体方法为：在母线 GEN3-230 和STNB-230 之间增加一回输电线，增加发电 3 的出力及其出口变压器的容量，新增或改造的元件如下图虚线所示：其基础数据如下：母线数据:母线名基准电压区域号电压上限电压下限单相短路容量三相短路容量发电1 16.5002 18.15014.8500.00000.0000发电2 18.0001 19.80016.2000.00000.0000发电3 13.8001 15.18012.4200.00000.0000GEN1-23 0 230.00002 0.0000 0.0000 0.00000.0000GEN2-23 0 230.00001 0.0000 0.0000 0.00000.0000GEN3-23 0 230.00001 0.0000 0.0000 0.00000.0000STNA-23 0 230.00002 0.0000 0.0000 0.00000.0000STNB-23 0 230.00002 0.0000 0.0000 0.00000.0000STNC-23 0 230.00001 0.0000 0.0000 0.00000.0000交流线数据数据组I测母线J测母线编号所属区单位正序阻抗正序电抗正序充电电纳零序电阻零序电抗零序充电电纳的1/2变压器数据续上表发电数据续上表负荷数据区域定义数据方案定义潮流计算作业定义三、实验步骤（1）点击|“电力系统分析综合程序（PSASP）”；（2）点击“创建”，创建文件；（3）点击“图形支持环境”；（4）点击“编辑模式”，可进行绘图和参数录入：a、绘制出所有母线，输入母线数据；b、添加发电机、负荷、交流线、变压器、支路，输入该元件数据；（5）关闭“编辑模式”窗口；（6）点击“运行模式”：（7）点击“作业”菜单项，执行“方案定义”命令（例如方案为1，数据组选择BASIC），点击“确定”。
电力系统潮流计算实验报告
１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.2000231）计算节点导纳矩阵由2000.00500.012j Z 71.418.112j y ;2000.00500.013j Z71.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;形成导纳矩阵BY ：71.418.10071.418.10071.418.171.418.171.418.171.418.142.936.2j j j j j j j j j Y B2）计算各PQ、PV 节点功率的不平衡量，及PV 节点电压的不平衡量：取:000.0000.1)0(1)0(1)0(1j jf e U000.0000.1)0(2)0(2)0(2j jf e U节点3是平衡节点，保持000.0000.1333j jf e U为定值。
nj j jij jij ijij jij i ieB fG f fB eG e P1)0()0()0()0()0()0()0(；２nj j jij jij ijij jij i ie B fG e f B eG f Q 1)0()0()0()0()0()0()0(；);(2)0(2)0(2)0(iiif e U)0.142.90.036.2(0.0)0.042.90.136.2(0.1)0(1P)0.171.40.018.1(0.0)0.071.40.118.1(0.1 )0.171.40.018.1(0.0)0.071.40.118.1(0.1 0.0 ；)0.142.90.036.2(0.1)0.042.90.136.2(0.0)0(1Q)0.171.40.018.1(0.1)0.071.40.118.1(0.0 )0.171.40.018.1(0.1)0.071.40.118.1(0.0 0.0 ；)0.171.40.018.1(0.0)0.071.40.118.1(0.1)0(2P)0.171.40.018.1(0.0)0.071.40.118.1(0.1 )0.00.00.00.0(0.0)0.10.00.10.0(0.1 0.0 ；101)(222)0(22)0(22)0(2f e U；于是：;)0()0(iiiP P P ;)0()0(iiiQQ Q);(2)0(2)0(22)0(iiiif e UU5.00.05.0)0(11)0(1P P P ;35.00.035.0)0(11)0(1QQ Q;4.00.04.0)0(22)0(2P P P ;1025.0)01(05.1)(2222)0(22)0(2222)0(2f e UU3）计算雅可比矩阵中各元素雅可比矩阵的各个元素分别为：３ji ij ji ij j i ij j i ij ji ij j i ij e U S f U R e Q L f Q J e P N f P H 22;;; 又： nj j jij jij i jij jij i ieB fG f fB eG e P1)0()0()0()0()0()0()0(； nj j jij jij ijij jij iieB fG e fB eG f Q 1)0()0()0()0()0()0()0(；);(2)0(2)0(2)0(iiif e U)0(1P )0(111)0(111)0(1)0(111)0(111)0(1e Bf G f f B e G e)0(212)0(212)0(1)0(212)0(212)0(1e B fG f f B e G e313313)0(1313313)0(1e Bf G f f B e G e ；)()()0(111)0(111)0(1)0(111)0(111)0(1)0(1e Bf Ge f B e G f Q)()()0(212)0(212)0(1)0(212)0(212)0(1e Bf G e f B e G f)()(313313)0(1313313)0(1e Bf G e f B e G f；)0(2P )0(121)0(121)0(2)0(121)0(121)0(2e Bf G f f B e G e)0(222)0(222)0(2)0(222)0(222)0(2eB fG f fBeG e323323)0(2323323)0(2e Bf G f f B e G e ；)(2)0(22)0(22)0(2f e U42.90.171.40.171.4313)0(212)0(1)0(1)0(11e B e Bf P H ； 36.20.118.10.118.10.136.222313)0(212)0(111)0(1)0(1)0(11 e G e G e G e P N 36.20.118.10.118.1313)0(212)0(1)0(1)0(11 e G e G f Q J４42.90.171.40.171.40.142.922313)0(212)0(111)0(1)0(1)0(11 e B e B e B e Q L 71.40.171.4)0(112)0(2)0(1)0(12 e B f P H ； 18.10.118.1)0(112)0(2)0(1)0(12 e G e P N ； 18.10.118.1)0(112)0(2)0(1)0(12 e G f Q J ；71.40.171.4)0(112)0(2)0(1)0(12 e B e Q L ； 71.40.171.4)0(221)0(1)0(2)0(21 e B f P H ； 11.40.111.4)0(221)0(1)0(2)0(21 e G e P N ； 0)0(12)0(2)0(21 f U R ； 0)0(12)0(2)0(21 e U S ； 71.40.10.00.171.4323)0(121)0(2)0(2)0(22 e B e B f P H ； 18.10.10.00.118.10.118.122323)0(121)0(222)0(2)0(2)0(22 e G e G e G e P N ；02)0(2)0(22)0(2)0(22 f f U R ； 0.20.122)0(2)0(22)0(2)0(22 e e U S ； 得到K=0时的雅可比矩阵：0.200018.171.418.171.471.418.142.936.218.171.436.242.9)0(J4）建立修正方程组：５)0(2)0(2)0(1)0(10.200011.4959.1011.4959.10959.1011.4918.2122.811.4959.1022.8918.210975.04.035.08.0e f e f 解得：04875.001828.00504.00176.0)0(2)0(2)0(1)0(1e f e f 因为 )0()0()1(iiie e e ； )0()0()1(iiif f f ；所以 9782.00218.00.1)0(1)0(1)1(1e e e ； 0158.00158.00)0(1)0(1)1(1f f f ；05125.105125.00.1)0(2)0(2)1(2e e e ；05085.005085.00)0(2)0(2)1(2f f f ；5）运用各节点电压的新值进行下一次迭代：即取: 0158.09782.0)1(1)1(1)1(1j jf e U05085.005125.1)1(2)1(2)1(2j jf e U节点3时平衡节点，保持000.0000.1333j jf e U为定值。
最新电力系统分析潮流计算实验报告
电力系统分析潮流计算实验报告姓名：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、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。