电力系统分析实验报告
1 湖北民族学院信息工程学院实验报告 班级：0312406 姓名：王愿喜 学号：031240612成绩： 实验时间：2015年4月 日 节 实验地点：实验楼 课程名称：电力系统分析 实验类型：设计型□ 验证型□ 综合型□ 实验题目：MATLAB 在电力系统中的应用实例 实验仪器： matlab
一、实验目的 熟悉matlab在电力系统潮流计算中的应用，加深对潮流计算的了解，学会自己使用matlab进行对电力系统的计算。
二、实验原理、原理图及电路图 以书中11-5为例，用m文件编写其支路、节点、发电机和基准容量的矩阵，然后运用快速解耦法解电力系统时用runpf 函数名函数即可。
三、实验步骤及内容 首先编写11-5题中支路、节点、发电机和基准容量的参数矩阵，如下图所示； function mpc = case181181 mpc.version='2'; mpc.baseMVA = 100; mpc.bus=[ k 1 300 180 0 0 1 1 0 100 1 1.1 0.94; 2 1 550 130 0 0 1 1 0 100 1 1.1 0.94; 3 2 500 0 0 0 1 1.10 0 100 1 1.1 0.94; 4 3 0 0 0 0 1 1.05 0 100 1 1.1 0.94; ]; mpc.gen=[ 3 500 0 9990 -9999 1.10 100 1 600 0; 4 0 0 9990 -9999 1.05 100 1 600 0; ]; mpc.branch=[ 2 1 0.10 0.40 0.03056 0 0 0 0 0 1 -360 360; 3 1 0.3 0 0 0 0 0 1.10 0 1 -360 360; 4 2 0.08 0.40 0.0384 0 0 0 0 0 1 -360 360; 4 1 0.12 0.50 0.02826 0 0 0 0 0 1 -360 360; ]; return;
四、实验结果记录与分析 2
>> runpf case181181 Newton's method power did not converge in 10 iterations. Did not converge (0.22 seconds) ================================================================================ | System Summary | ================================================================================
How many? How much? P (MW) Q (MVAr) --------------------- ------------------- ------------- ----------------- Buses 4 Total Gen Capacity 1200.0 -19998.0 to 19980.0 Generators 2 On-line Capacity 1200.0 -19998.0 to 19980.0 Committed Gens 2 Generation (actual) 1675.2 738.4 Loads 3 Load 1350.0 310.0 Fixed 3 Fixed 1350.0 310.0 Dispatchable 0 Dispatchable -0.0 of -0.0 -0.0 Shunts 0 Shunt (inj) -0.0 0.0 Branches 4 Losses (I^2 * Z) 1102.94 4482.62 Transformers 1 Branch Charging (inj) - 42.3 Inter-ties 0 Total Inter-tie Flow 0.0 0.0 Areas 1
Minimum Maximum ------------------------- -------------------------------- Voltage Magnitude 1.050 p.u. @ bus 4 2.977 p.u. @ bus 2 Voltage Angle -121.58 deg @ bus 3 0.00 deg @ bus 4 P Losses (I^2*R) - 598.33 MW @ line 4-2 Q Losses (I^2*X) - 2991.63 MVAr @ line 4-2
================================================================================ | Bus Data | ================================================================================ Bus Voltage Generation Load # Mag(pu) Ang(deg) P (MW) Q (MVAr) P (MW) Q (MVAr) ----- ------- -------- -------- -------- -------- -------- 1 1.653 -120.628 - - 300.00 180.00 2 2.977 -113.355 - - 550.00 130.00 3 1.100 -121.579 500.00 9.15 500.00 0.00 3
4 1.050 0.000 1175.17 729.22 - - -------- -------- -------- -------- Total: 1675.17 738.37 1350.00 310.00
================================================================================ | Branch Data | ================================================================================ Brnch From To From Bus Injection To Bus Injection Loss (I^2 * Z) # Bus Bus P (MW) Q (MVAr) P (MW) Q (MVAr) P (MW) Q (MVAr) ----- ----- ----- -------- -------- -------- -------- -------- -------- 1 2 1 380.81 886.53 -273.06 -473.23 107.754 431.02 2 3 1 -217.73 9.15 360.19 -9.15 142.465 0.00 3 4 2 802.49 422.83 -204.16 2549.66 598.325 2991.63 4 4 1 372.68 306.39 -118.29 748.17 254.395 1059.98 -------- -------- Total: 1102.939 4482.62 4
湖北民族学院信息工程学院实验报告 班级：0312406 姓名：王愿喜 学号：031240612成绩： 实验时间：2015年4月 日 节 实验地点：实验楼 课程名称：电力系统分析 实验类型：设计型□ 验证型□ 综合型□ 实验题目：MATLAB matpower在电力系统潮流计算分析中的应用 实验仪器： matlab
一、实验目的 通过simulink建模建立电力系统潮流计算时的模型，运用powergui中的潮流计算对该系统进行分析。
二、实验原理、原理图及电路图 在模型库中的powergui中有对该系统进行潮流计算的选项，可以在里面设置后对该系统进行潮流计算分析。
三、实验步骤及内容
搭建例题的模型如下图所示：
Discrete,Ts = 5e-005 s.
Terminator4
Terminator3Terminator2Terminator1
Terminator
ABCabcT2ABCabcT1ABCabcM8ABCabcM7ABCabcM6ABCabcM5ABCabcM4ABCabcM3ABCabcM2ABCabcM1ABCABCL3ABCABCL2ABCABCL1
PmVf_mABCG2PmVf_mABCG1-C-Constant3-C-Constant2-C-Constant1-C-Constant
mABC Load3
mABC Load2mABC Load1
在powergui中对其进行相应的设置： 5 结果可得：
Machine: G2 Nominal: 100 MVA 10.5 kV rms Bus Type: Swing bus Uan phase: 0.00° Uab: 11025 Vrms [1.05 pu] 30.00° Ubc: 11025 Vrms [1.05 pu] -90.00° Uca: 11025 Vrms [1.05 pu] 150.00° Ia: 14443 Arms [2.627 pu] -28.65° Ib: 14443 Arms [2.627 pu] -148.65° Ic: 14443 Arms [2.627 pu] 91.35° P: 2.4203e+008 W [2.42 pu] Q: 1.3225e+008 Vars [1.323 pu] Pmec: 2.44e+008 W [2.44 pu] Torque: 1.5534e+007 N.m [2.44 pu] Vf: 3.9054 pu Machine: G1 Nominal: 1000 kVA 10.5 kV rms Bus Type: P & V generator Uan phase: 3.22°
电力系统分析实习报告
一、实习背景随着我国经济的快速发展和科技的不断进步，电力系统在国民经济和社会生活中的地位日益重要。
为了更好地了解电力系统的运行规律和特性，提高电力系统的运行效率和安全性，我于近期在XX电力公司进行了为期四周的电力系统分析实习。
二、实习单位简介XX电力公司是一家集电力生产、输电、配电、售电、科研、咨询、技术服务于一体的国有大型企业。
公司承担着XX地区及周边地区的电力供应任务，是保障地区电力安全稳定供应的重要力量。
三、实习目的通过本次实习，我旨在：1. 了解电力系统的基本组成和运行原理；2. 掌握电力系统分析的基本方法和技术；3. 提高电力系统运行、维护和管理的实际操作能力；4. 增强团队合作和沟通协调能力。
四、实习内容1. 电力系统基本知识学习实习期间，我首先学习了电力系统的基础知识，包括电力系统的组成、电力设备的基本原理、电力系统的运行方式等。
通过学习，我对电力系统有了初步的认识。
2. 电力系统分析软件操作为了更好地进行电力系统分析，我学习了电力系统分析软件的使用方法。
在导师的指导下，我熟练掌握了电力系统分析软件的基本操作，如潮流计算、短路计算、稳定性分析等。
3. 电力系统现场实习在导师的带领下，我参观了电力系统的现场，了解了变电站、输电线路、配电线路等设备的运行状况。
同时，我还参与了电力系统的现场操作，如变电站的设备巡检、故障处理等。
4. 电力系统分析案例研究在实习过程中，我选取了几个典型的电力系统分析案例进行深入研究。
通过对案例的分析，我掌握了电力系统分析的基本方法和技巧。
五、实习成果1. 理论知识掌握通过实习，我对电力系统的基本知识和运行原理有了更加深入的了解，为今后的学习和工作打下了坚实的基础。
2. 实际操作能力提升在实习过程中，我参与了电力系统的现场操作，提高了自己的实际操作能力。
3. 团队合作和沟通协调能力增强在实习过程中，我与同事们共同完成了多项任务，增强了团队合作和沟通协调能力。
南邮通达电力系统实验报告
单机无穷大系统稳态实验：一、整理实验数据，说明单回路送电和双回路送电对电力系统稳定运行的影响，并对实验结果进行理论分析：实验数据如下:由实验数据，我们得到如下变化规律：（1）保证励磁不变的情况下，同一回路,随着有功输出的增加，回路上电流也在增加，这是因为输出功率P=UIcos ，机端电压不变所以电流随着功率的增加而增加；(2）励磁不变情况下，同一回路，随着输出功率的增大，电压损耗在增大;这是由于电压降落△U=（PR+QX）/U，而横向分量较小，所以电压损耗也随着输出功率的增大而增大.单回路供电和双回路供电对电力系统稳定性均有一定的影响，其中双回路要稳定一些,单回路稳定性较差.二、根据不同运行状态的线路首、末端和中间开关站的实验数据、分析、比较运行状态不同时，运行参数变化的特点和变化范围。
由实验数据，我们可以得到如下结论:（1）送出相同无功相同有功的情况下:单回路所需励磁电压比双回路多,线路电流大小相等，单回路的电压损耗比双回路多；（eg。
P=1,Q=0.5时)（2)送出相同无功的条件下，双回路比单回路具有更好的静态稳定性，双回路能够输送的有功最大值要多于单回路；发生这些现象的原因是:双回路电抗比单回路小，所以所需的励磁电压小一些，电压损耗也要少一些，而线路电流由于系统电压不改变；此外，由于电抗越大，稳定性越差，所以单回路具有较好的稳定性。
三、思考题:1、影响简单系统静态稳定性的因素是哪些?答：由静稳系数S Eq=EV/X，所以影响电力系统静态稳定性的因素主要是：系统元件电抗，系统电压大小，发电机电势以及扰动的大小。
2、提高电力系统静态稳定有哪些措施？答：提高静态稳定性的措施很多，但是根本性措施是缩短"电气距离”.主要措施有：（1)、减少系统各元件的电抗:减小发电机和变压器的电抗，减少线路电抗（采用分裂导线）;(2）、提高运行电压水平;(3）、改善电力系统的结构;(4）、采用串联电容器补偿；(5）、采用自动励磁调节装置；(6）、采用直流输电。
电力各种实验报告
一、实验目的本次实验旨在使学生掌握电力系统的基本原理和实验方法，提高学生的实际操作能力，加深对电力系统理论知识的理解。
通过实验，学生应能熟练操作电力系统实验设备，掌握电力系统稳态和暂态分析的基本方法，并能够分析实验结果，得出正确的结论。
二、实验内容1. 电力系统稳态分析实验（1）实验目的：掌握电力系统稳态分析的基本方法，了解电力系统稳态运行的特点。
（2）实验内容：测量电力系统中的电压、电流、功率等参数，分析电力系统稳态运行的特点。
（3）实验步骤：①接线：按照实验电路图，正确连接实验设备。
②调试：检查实验设备是否正常，调整实验参数。
③测量：测量电力系统中的电压、电流、功率等参数。
④数据分析：对实验数据进行处理和分析，得出实验结论。
2. 电力系统暂态分析实验（1）实验目的：掌握电力系统暂态分析的基本方法，了解电力系统暂态运行的特点。
（2）实验内容：测量电力系统在故障发生时的电压、电流、功率等参数，分析电力系统暂态运行的特点。
（3）实验步骤：①接线：按照实验电路图，正确连接实验设备。
②调试：检查实验设备是否正常，调整实验参数。
③模拟故障：模拟电力系统故障，测量故障发生时的电压、电流、功率等参数。
④数据分析：对实验数据进行处理和分析，得出实验结论。
3. 电力系统保护实验（1）实验目的：掌握电力系统保护的基本原理和实验方法，提高学生的实际操作能力。
（2）实验内容：学习电力系统保护的基本原理，了解不同保护装置的工作原理。
（3）实验步骤：①学习：学习电力系统保护的基本原理，了解不同保护装置的工作原理。
②实验：按照实验电路图，正确连接实验设备，进行保护实验。
③数据分析：对实验数据进行处理和分析，得出实验结论。
4. 电力系统仿真实验（1）实验目的：掌握电力系统仿真软件的使用方法，提高学生的实际操作能力。
（2）实验内容：使用电力系统仿真软件进行电力系统稳态和暂态分析。
（3）实验步骤：①学习：学习电力系统仿真软件的使用方法。
电分实验 电力系统分析实验
电气工程学院《电力系统分析综合实验》2017年度PSASP实验报告学号：2014302540149姓名：刘玉清班级： 2014级6班一、实验目的（不超过400字介绍）潮流计算是电力系统非常重要的分析计算，用以研究系统运行和规划中提出的各种问题。
了解用PSASP进行电力系统各种计算的方法，掌握用PSASP进行电力系统潮流、稳定和短路计算的方法，对于我们学习电力系统分析有重要的意义，可以让我们更深刻地体会潮流计算的过程，方便地模拟仿真电力系统的运行状况。
对运行中的电力系统，通过潮流的计算可以分析各种负荷变化、网络结构改变等各种情况会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件（线路、变压器等）是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等；对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案（如新建变电站、线路改造、电磁环网解环等）能否满足安全稳定运行的要求。
二、潮流计算部分（1）简单介绍本次实验的潮流计算的试验内容，叙述利用PSASP软件进行潮流计算需要输入哪些数据（不超过800字）。
以上为系统常规运行方式的单线图。
由于母线 STNB-230 处负荷的增加，需对原有电网进行改造，具体方法为：在母线 GEN3-230 和 STNB-230 之间增加一回输电线，增加发电 3 的出力及其出口变压器的容量，新增或改造的元件如右图虚线所示：实验将对该系统进行计算分析。
需要输入的数据为：母线数据、交流线数据、变压器数据、发电数据、负荷数据、区域定义数据、方案定义、潮流计算作业定义。
具体数据见实验指导书。
（2）绘制仿真实验任务指导书中要求的两个潮流计算拓扑结构图，并根据两种潮流方案的潮流计算结果在拓扑结构图上标注母线节点电压和相位，利用复数功率形式标注各支路功率（参考方向选择数据表格中各支路的i侧母线至j侧），统计发电机和负荷的功率大小并统计系统网损（3）在规划潮流方式下，增加STNC-230母线负荷的有功至1.5p.u.，无功保持不变，试调试潮流，并在绘制的规划方式拓扑结构上标注母线节点电压和相位，利用复数功率形式标注各支路功率（参考方向选择数据表格中各支路的i侧母线至j侧），统计发电机和负荷的功率大小并统计系统网损。
电力系统稳态分析实验报告
电力系统稳态分析实验报告篇一：电力系统稳态分析实验指导书电力系统稳态分析实验指导书目录实验一单机－无穷大系统稳态运行方式实验 ................................................ ........................................ 2 1.1 实验目的................................................. ................................................... ........................................... 2 1.2 原理说明................................................. ................................................... ........................................... 2 1.3 实验内容与步骤 ................................................ ................................................... ................................ 3 实验二电力系统潮流计算分析实验 ................................................ ....................................................... 6 2.1 实验目的................................................. ................................................... ........................................... 6 2.2 原理说明................................................. ................................................... ........................................... 6 2.3 实验内容与步骤 ................................................ ................................................... .. (6)I实验一单机－无穷大系统稳态运行方式实验1.1 实验目的1．熟悉远距离输电的线路基本结构和参数的测试方法。
四川大学电力系统实验报告----暂态稳定实验
电力系统暂态稳定实验一、实验目的①通过实验加深对电力系统暂态稳定内容的理解，使课堂理论教学与实践结合，提高学生的感性认识。
②学生通过实际操作，从实验中观察到系统失步现象和掌握正确处理的措施③用数字式记忆示波器测出短路时短路电流的非周期分量波形图，并进行分析。
二、原理与说明电力系统暂态稳定问题是指电力系统受到较大的扰动之后，各发电机能否继续保持同步运行的问题。
在各种扰动中以短路故障的扰动最为严重。
正常运行时发电机功率特性为：P1＝（Eo×Uo）×sinδ1/X1；短路运行时发电机功率特性为：P2＝（Eo×Uo）×sinδ2/X2；故障切除发电机功率特性为：P3＝（Eo×Uo）×sinδ3/X3；对这三个公式进行比较，我们可以知道决定功率特性发生变化与阻抗和功角特性有关。
而系统保持稳定条件是切除故障角δc小于δmax，δmax可由等面积原则计算出来。
本实验就是基于此原理，由于不同短路状态下，系统阻抗X2不同，同时切除故障线路不同也使X3不同，δmax也不同，使对故障切除的时间要求也不同。
同时，在故障发生时及故障切除通过强励磁增加发电机的电势，使发电机功率特性中Eo增加，使δmax增加，相应故障切除的时间也可延长；由于电力系统发生瞬间单相接地故障较多，发生瞬间单相故障时采用自动重合闸，使系统进入正常工作状态。
这二种方法都有利于提高系统的稳定性。
三、实验电路图一次系统接线图四、实验内容及实验步骤●短路对电力系统暂态稳定的影响1．短路类型对暂态稳定的影响本实验台通过对操作台上的短路选择按钮的组合可进行单相接地短路，两相相间短路，两相接地短路和三相短路试验。
固定短路地点，短路切除时间和系统运行条件，在发电机经双回线与“无穷大”电网联网运行时，某一回线发生某种类型短路，经一定时间切除故障成单回线运行。
短路的切除时间在微机保护装置中设定，同时要设定重合闸是否投切。
武大电气《电力系统分析综合实验》2019年度PSASP实验报告
电气工程与自动化学院《电力系统分析综合实验》2019年度PSASP实验报告学号：姓名：班级：1、阐述基于PSASP的电力系统分析综合实验的目的。
实验目的：掌握用PSASP进行电力系统潮流计算，短路计算，暂态稳定计算。
（1）潮流计算可以为短路计算和暂态稳定计算提供初始状态，是电力系统计算中的基本计算，要求掌握软件的操作步骤，并对比分析牛顿拉夫逊法和PQ分解法的区别，在实验过程中体会PQ分解法相比牛顿拉夫逊法的特点。
（2）短路计算的目的要求根据数据结合对称分量法加深对于短路计算的理论知识的理解。
（3）暂态稳定计算里最关键的是故障极限切除时间的确定，加深对复杂电力系统暂态的判定的认识。
2、简要阐述本实验课程的主要实验任务（1）掌握用PSASP对电力系统进行建模。
（2）潮流计算，包括对常规方式和规划方式的电力系统进行潮流计算。
（3）短路计算，基于潮流作业1和2等5个单相接地短路、AB两相短路、复杂故障短路计算等短路计算并分析结果。
（4）暂态计算，基于潮流作业1和2的瞬时故障进行暂态稳定计算并分析结果。
3、实验方案原理图介绍。
图1（a）常规方式（b）规划方式以上为系统常规运行方式的单线图。
由于母线STNB-230 处负荷的增加，需对原有电网进行改造，具体方法为：在母线GEN3-230 和STNB-230 之间增加一回输电线，增加发电3 的出力及其出口变压器的容量，新增或改造的元件如下图虚线所示：4、计算分析用建模数据的整理表1母线数据5、按照下列作业要求，完成计算分析实验作业。
(1)基于实验二的潮流计算，对牛顿法和PQ法的原理做比较性的说明。
表6 常规方式下PQ法和NR法的潮流计算摘要信息报表表7 常规方式下PQ法和NR法的全网母线（发电、负荷）结果报表牛顿拉夫逊法每次都对电压幅值和相位进行修正，且每次计算MAX(DVR,DVI),判断是否小于允许误差0.0001，满足条件时停止迭代。
PQ分解法利用交流高压输电网中输电线路电抗远大于电阻的特点，对于牛顿拉夫逊法修正方程式的系数矩阵进行简化，节点的有功功率不平衡量只用于修正电压的相位，节点的无功功率不平衡量只用于修正电压的幅值，单次迭代计算量小，两个步骤分别轮流迭代，分别计算MAX(DP),MAX(DQ)，最终保证两者都小于允许误差0.0001。
最新电力系统分析潮流计算实验报告
电力系统分析潮流计算实验报告姓名：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. 系统经济性优化为了优化系统的经济性，我们进行了成本效益分析。
通过调整系统的运行策略和参数，我们评估了不同方案下的成本和效益。
实验结果显示，通过合理的调整发电机的输出功率和输电线路的容量，可以降低系统的运行成本，并提高系统的经济效益。
四、实验结论与建议通过本次仿真实验，我们得出了以下结论：1. 系统在大部分情况下表现出良好的稳定性和可靠性，但仍存在一些潜在的问题需要解决。
2. 加强系统的容错性和冗余性可以提高系统的可靠性。
3. 通过合理的调整运行策略和参数，可以降低系统的运行成本，并提高经济效益。
电力故障分析实验报告
一、实验目的1. 理解电力系统故障的基本概念和分类。
2. 掌握电力系统故障分析的基本方法和步骤。
3. 通过实验，加深对电力系统故障现象、故障类型、故障分析和处理方法的理解。
4. 培养实验操作能力和分析解决问题的能力。
二、实验原理电力系统故障是指电力系统中发生的各种异常现象，如短路、接地故障等。
故障分析是电力系统安全稳定运行的重要环节，主要包括故障现象的识别、故障类型的判断、故障分析计算、故障处理措施等。
三、实验仪器与设备1. 电力系统故障模拟装置2. 电流表、电压表3. 计算机及仿真软件4. 电力系统故障分析手册四、实验步骤1. 故障现象观察（1）打开电力系统故障模拟装置，观察正常运行的电力系统。
（2）模拟不同类型的故障，如三相短路、单相接地短路等，观察故障现象。
2. 故障类型判断（1）根据故障现象，判断故障类型，如短路、接地故障等。
（2）分析故障原因，如设备故障、操作失误等。
3. 故障分析计算（1）根据故障类型和故障参数，选择合适的故障分析计算方法。
（2）使用计算机和仿真软件进行故障分析计算，得到故障电流、电压等参数。
4. 故障处理措施（1）根据故障分析结果，制定故障处理措施。
（2）模拟故障处理过程，验证处理措施的有效性。
五、实验结果与分析1. 故障现象观察在实验过程中，我们模拟了三相短路、单相接地短路等故障类型。
观察到故障现象如下：（1）三相短路：故障点附近电流和电压急剧上升，保护装置动作，断路器跳闸。
（2）单相接地短路：故障点附近电流和电压发生变化，保护装置动作，断路器跳闸。
2. 故障类型判断根据故障现象，我们判断出故障类型为三相短路和单相接地短路。
3. 故障分析计算使用仿真软件对故障进行分析计算，得到以下结果：（1）三相短路：故障点电流约为额定电流的20倍，故障点电压约为额定电压的0.5倍。
（2）单相接地短路：故障点电流约为额定电流的10倍，故障点电压约为额定电压的1.7倍。
4. 故障处理措施根据故障分析结果，我们制定了以下故障处理措施：（1）三相短路：立即隔离故障点，修复设备，恢复供电。