华电电力系统暂态上机报告书(程序完整无误,带注释)
华电潮流上机课程设计报告+程序
院系：电气与电子工程学院班级：学号：学生姓名：指导教师：刘宝柱设计周数：成绩：日期：2012年1月5日一、课程设计的目的与要求培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识二、设计正文（详细内容见附录）1．手算2．计算机计算3．思考题3.1潮流计算的方法有哪些？各有何特点？答：潮流计算的方法主要有高斯-赛德尔迭代法、牛顿-拉夫逊迭代法和P-Q分解法。
它们各自的特点如下：（1）高斯-赛德尔迭代法分为以节点导纳矩阵为基础的高斯-赛德尔迭代法和以以节点阻抗矩阵为基础的高斯-赛德尔迭代法的原理比较简单，要求的数字计算机的内存量也比较小，但收敛性差，当系统规模变大时，迭代次数急剧上升，往往出现迭代不收敛的情况；而阻抗法改善了电力系统潮流计算导纳法德收敛性问题，在当时获得了广泛的应用，但是，阻抗法的主要缺点是占用计算机的内存很大，每次迭代的计算量很大。
当系统不断扩大时，这些缺点就更加突出。
（2）牛顿-拉夫逊法是数学中求解非线性方程式的典型方法，有较好的收敛性。
只要在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿潮流计算程序的计算效率，牛顿法在收敛性、内存要求、计算速度方面都超过了阻抗法，成为知道目前仍被广泛采用的方法。
（3）P-Q分解法潮流计算派生于以极坐标表示时的牛顿-拉夫逊法，它根据电力系统的特点，抓住主要矛盾，对纯数学的牛顿法进行了改造。
与牛顿法相比，P-Q分解法的修正方程的系数矩阵B’和B”分别是（n-1）和（m-1）的方阵，替代了原有的（n+m-2）阶系数矩阵J；B’、B”在迭代过程中可以保持不变且为对称的系数矩阵，提高了计算速度，降低了对存储容量的要求。
P-Q分解法在计算速度方面有显著地提高，迅速得到了推广。
3.2如果交给你一个任务，请你用已有的潮流计算软件计算北京城市电网的潮流，你应该做哪些工作？（收集哪些数据，如何整理，计算结果如何分析）答：(1)在进行北京城市电网的潮流计算之前需要了解北京城市电网中所有的节点支路的相关数据，并对节点和支路分类。
电力系统暂态上机计算教学指导书(适用于matlab)
电力系统暂态上机计算教学指导书(适用于Ｍaｔlab）一、基础知识1．为完成本次设计，需要掌握的Maｔｌab的基础知识有：数组的创建和计算,循环语句,条件选择语句,条件判断语句,绘制二维图的语句;2．提供了Ｍaｔｌab的学习参考《精通Matlab6.5》光盘版；3．Ｍａtlab是区分大小写的;4．本次设计需要学习和掌握的专业知识:《电力系统暂态分析》相关知识;5．ｍatlaｂ在三角函数计算时默认是弧度,而不是度;6．Maｔlａb 数组（矩阵)的下标从1开始;7．Ｍatlab 默认小数点后保留4位;8．如果不想显示变量的数值,则在该行末尾加“;”符号。
三、程序设计指导提供了两个程序,（1)不对称短路,名字是buｄuichen.m;(２）稳定计算的程序,名字是wendinｇ.m。
所有程序里“%”表示该行为注释行。
1．不对称短路程序fｕｎction buｄｕｉchen%数据来源于教材《电力系统暂态分析》P12９例（5－1)%给学生提供的不对称短路故障的部分程序。
%要求的工作（１）自动形成导纳阵;（2)完善短路类型，两相相间短路和两相短路接地；（３）求短路后发电机的电压cleaｒ %清除所有变量的值％输入支路各序阻抗,ZZ１表示支路正序阻抗,ZZ２表示支路负序阻抗，ZZ０表示支路零序阻抗ZZ1(1,2)=j*0.1；ZZ1(1,３)=j＊0.1；ＺZ１（2,3)=j＊0.1;ZZ2(1,2）=j＊0.１;ZZ2（2,３)=j*0.１;ZZ2（1，3)=j*０．1；ZZ0（1,2)=ｊ*0．2;ZZ0（2,3)=j*0．2;ZZ０(1,3）=j*0．2;％输入节点导纳矩阵,Y１表示不计发电机节点的正序网络节点导纳阵,Y2表示不计发电机节点的负序网络节点导纳阵,Y0表示不计发电机节点的零序网络节点导纳阵，Y1＝[-j*２6.62６6 j*1０ j*10j＊10 -j*33.293３ j*１０j＊1０ j*1０-j＊19.96];Y２=[-ｊ*26．62６6 j*10 ｊ＊10ｊ*10 -j＊33.2933 j*10j*10 j＊10 -j*19.96];Ｙ０=[－ｊ*30 ｊ*５ j*5ｊ*５－j*50 j*５ｊ＊5 j*５ -j＊１0];ＹY１=[-j＊39.９6 ｊ*10 j*10 j*2０ 0 %YＹ1表示计及发电机节点的正序网络节点导纳阵，YＹ2表示计及发电机节点的负序网络节点导纳阵j＊１０ -ｊ＊59.９６ j*10 0 j*4０j＊10 ｊ＊10 -j*１9．96 0 0 j＊２００0 -ｊ*３0 0 ０ j*4０ 0 ０-j*60］；YY２＝YY1；N１＝3; %输入网络的节点数N2＝5；%输入包括所有发电机节点的网络的节点数k＝3；％输入短路点的节点号fault=1； %输入短路类型 ;f（1)＝１单相短路接地；f（2)=２两相相间短路;f（１，１)＝3两相短路接地%第一部分:计算所有节点的a,b,c三相电压for p=1：N1ｉf p==ｋI（p)=1;ｅlseI(p)=0;eｎdenｄ%正序Z1=eye(N１)／Y1;% eye(N１)是3阶的单位矩阵,求节点导纳矩阵Y1的逆阵,即节点阻抗矩阵ZZ1(：,k)=Ｚ１*I';%短路点k的自阻抗和互阻抗Zk1＝Z１(:,k）;%负序Z2=eｙe（N１)／Y２;Z２(:，k)=Ｚ2*I';Zk２=Z2(：，k);%零序Z0=eye(N１)/Ｙ0;Z0(:,k)=Ｚ0*I＇;Ｚk0=Z0(:,k);％计根据故障类型选择不同的计算公式，算故障点各序电流if ｆaulｔ==1Iｋ１=1／（Z１（k,ｋ）+Z2（k,k)+Z0(k，k));Ｉk2＝Ik1；Ｉｋ0＝Ik1；elｓeiｆfａult==2elｓeiｆ fault==3ｅndeｎdendend％计算各序电流的故障分量ｆoｒp=1:Ｎ１pif p==kＩ１(p)=-1*Ik1；I2(ｐ）=－１＊Ik2；I0（p)=－１*Ik０;elseI1(ｐ)=0;I2(p)=０;I０(p)=0;eｎdend%计算故障点各序电压的故障分量P９6公式（４-１3)ｕu1（:,k)＝Ｙ1\I１.＇; ％“\”表示左除 A\Ｂ=inv(Ａ)*B ｉnv是求逆阵。
电力系统静态稳定暂态稳定实验报告
电力系统静态、暂态稳定实验报告一、实验目的1．了解和掌握对称稳定情况下，输电系统的各种运行状态与运行参数的数值变化范围；2．通过实验加深对电力系统暂态稳定内容的理解3．通过实际操作，从实验中观察到系统失步现象和掌握正确处理的措施二、原理与说明实验用一次系统接线图如图1所示：图1. 一次系统接线图实验中采用直流电动机来模拟原动机，原动机输出功率的大小，可通过给定直流电动机的电枢电压来调节。
实验系统用标准小型三相同步发电机来模拟电力系统的同步发电机，虽然其参数不能与大型发电机相似，但也可以看成是一种具有特殊参数的电力系统的发电机。
发电机的励磁系统可以用外加直流电源通过手动来调节，也可以切换到台上的微机励磁调节器来实现自动调节。
实验台的输电线路是用多个接成链型的电抗线圈来模拟，其电抗值满足相似条件。
“无穷大”母线就直接用实验室的交流电源，因为它是由实际电力系统供电的，因此，它基本上符合“无穷大”母线的条件。
为了进行测量，实验台设置了测量系统，以测量各种电量（电流、电压、功率、频率）。
为了测量发电机转子与系统的相对位置角（功率角），在发电机轴上装设了闪光测角装置。
此外，台上还设置了模拟短路故障等控制设备。
电力系统静态稳定问题是指电力系统受到小干扰后，各发电机能否不失同步恢复到原来稳定状态的能力。
在实验中测量单回路和双回路运行时，发电机不同出力情况下各节点的电压值，并测出静态稳定极限数值记录在表格中。
电力系统暂态稳定问题是指电力系统受到较大的扰动之后，各发电机能否过渡到新的稳定状态，继续保持同步运行的问题。
在各种扰动中以短路故障的扰动最为严重。
正常运行时发电机功率特性为：P1＝（Eo×Uo）×sinδ1/X1；短路运行时发电机功率特性为：P2＝（Eo×Uo）×sinδ2/X2；故障切除发电机功率特性为：P3＝（Eo×Uo）×sinδ3/X3；对这三个公式进行比较，我们可以知道决定功率特性发生变化与阻抗和功角特性有关。
新版电力系统暂态上机计算课程设计报告完整无误版附程序
课程设计报告( 2014—2015年度第二学期)名称：电力系统暂态上机计算院系：电气与电子工程学院班级：电气1211学号：1121210205学生姓名：郝阳指导教师：陈艳波设计周数：两周成绩：日期：2015年7月4日一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。
二、设计正文（详细内容见附录,用A4纸，页数不限）1．对称短路计算过程流程图和计算结果2．不对称短路计算过程流程图和计算结果3．静态稳定性计算过程流程图和计算结果4．暂态稳定性计算过程流程图和计算结果5．思考题三、课程设计总结或结论本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。
通过本次课程设计，本人对电力系统故障分析有了更深刻的理解，电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。
在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。
在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。
同时我对用来分析电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。
这与手算系统短路电流时使用的网络化简方法大大不同。
在学习中，参照潮流程序，我加深了对节点导纳矩阵建立方法的理解与学习，巩固了不同类型短路的短路电流计算方法，和序电压、序电流，相电压、相电流基于matlab软件的计算。
利用小干扰分析法判别静态稳定可以使用劳斯判据，也可以使用特征根判别的方法。
二者都能够判断，在书中的例题中使用的是劳斯判据，但是，特征根判别放法更适合在matlab软件的环境下使用，因为该软件提供了计算矩阵特征根的函数，使用起来简便易懂。
华北电力大学电力系统第二章暂态作业
暂态第七章
7-1 什么是静态稳定问题？静稳储备系数的定义是什么？ 7-2 电力系统静态稳定的判据是什么？（Eq 恒定，单机对无穷大系统） 7-3 发电机的自动励磁调节系统改善静态稳定性的原理是什么？ 7-4 如下图所示，写出发电机的二阶动态方程（不考虑阻尼）。其中发电机为隐极机，电磁
暂态第 2 章
2-1 什么是 Park 变换？ 2-2 写出同步发电机定子、转子绕组正电流产生正磁链的有阻尼绕组发电机 d、q、0 系统
的磁链方程。
2-3 请给出同步发电机的各种电抗参数的等值电路和表达式，并排出大小顺序。 xd ， xq ，
xd′ ， xd′′ ， xq′′ ， x f ， xD ， xQ
2-4 请给出同步发电机的各种电势的等值电路、表达式及相量图。 Eq ， EQ ， E′ ， Eq′ ，
E′′ ， Eq′′ ， Ed′′
2-5 试写出同步发电机（隐极机和凸极机）稳态运行时的基本方程，并画出等值电路和相 量图。
2-6 某同步发电机的参数为： f N= 50Hz ， xσ = 0.15 ， xd = 1.1 ， xq = 0.7 ， x f = 1.0 ，
什么？
6-3 什么是电力系统的稳定问题？P154
6-4 简单电力系统接线如图所示，
I
U =1
其中单机无穷大母线电压为U = 1 ，功率因数为 0.98，发电机输出功率为 1.0，单根线路电
抗为 0.6。变压器电抗为 0.25。假设发电机为隐极机，同步电抗为 1.0，暂态电抗为 0.2。试 计算：
电力系统分析电力系统暂态实验报告
电力系统分析电力系统暂态实验报告备注：序号（一）、（二）、（三）为实验预习填写项v+-voutv +-vin Continuous pow erguii +-iline Vin1PI Section LineCurrentBreakerclose at 0.02 sec.2.572e5 v60 Hz10 ohms图6-14 Powergui模块参数设置（2）Ac Voltage Source模块：电源电压2.57kV，60Hz，参考实验一，参数设置如图6-15所示。
图6-15 AC Voltage Source模块参数设置（3）电阻：参考实验一，参数设置10欧姆。
（4）Breaker模块：断路器，设置内阻0.001欧姆，0.02秒闭合。
（a）（b）图6-16 Breaker模块参数设置（5）输电线П模型模块（Pi Section Line）：100公里，60Hz,串联电阻0.2568Ω串联感抗2mH ，并联容抗8.6nF 。
（a）(b)图6-17 Pi Section Line（输电线П模型）模块参数设置（6）其余量测元件参考实验一【参考波形】U1与U2电压：图6-18 U1与U2电压波形图输电线电流I line电流响应与局部放大响应：（a）（b）图6-19 输电线电流I line电流响应与局部放大响应五、程序调试及实验总结程序调试：在实验过程中，电阻模块没有注意看实验二电阻值的设定值，把设置成了实验一的电阻值，导致输出波形图跟实验指导给出参考的波形图相差甚大。
后面问了同学更改过来后能输出正确的波形了。
实验总结：通过这一次的实验，让我加深对于电力系统暂态稳定内容的理解，使得我能把课堂理论教学的知识与实践相结合起来，进而提高我对电力系统暂态表现的认识。
也让我通过进行实际操作的过程中，从实验中观察到系统暂态响应发生时的现象和掌握正确处理的措施，并使用用MATLAB/Simulink来观测输出的波形图，并进行分析。
新版电力系统暂态上机计算课程设计报告完整无误版附程序
课程设计报告( 2014—2015年度第二学期)名称：电力系统暂态上机计算院系：电气与电子工程学院班级：电气1211学号：1121210205学生姓名：郝阳指导教师：陈艳波设计周数：两周成绩：日期：2015年7月4日一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。
二、设计正文（详细内容见附录,用A4纸，页数不限）1．对称短路计算过程流程图和计算结果2．不对称短路计算过程流程图和计算结果3．静态稳定性计算过程流程图和计算结果4．暂态稳定性计算过程流程图和计算结果5．思考题三、课程设计总结或结论本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。
通过本次课程设计，本人对电力系统故障分析有了更深刻的理解，电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。
在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。
在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。
同时我对用来分析电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。
这与手算系统短路电流时使用的网络化简方法大大不同。
在学习中，参照潮流程序，我加深了对节点导纳矩阵建立方法的理解与学习，巩固了不同类型短路的短路电流计算方法，和序电压、序电流，相电压、相电流基于matlab软件的计算。
利用小干扰分析法判别静态稳定可以使用劳斯判据，也可以使用特征根判别的方法。
二者都能够判断，在书中的例题中使用的是劳斯判据，但是，特征根判别放法更适合在matlab软件的环境下使用，因为该软件提供了计算矩阵特征根的函数，使用起来简便易懂。
华北电力大学电力系统暂态分析教学大纲
华北电力大学电力系统暂态分析教学大纲 2007-07-30 考研论坛卡希·闹【字体：大中小】课程的目的和任务本课程是电气工程及其自动化专业的专业必修课也是本专业学习后续其它专业课及专业选修课的基础1 使学生全面了解由短路故障引起的电力系统电磁暂态过程及各种扰动引发的电力系统机电暂态过程2 使学生了解同步发电机短路的物理过程数学模型以及发电机短路无限大系统短路不对称短路的计算方法3 使学生了解电力系统基本元件模型了解静态稳定暂态稳定电压稳定的基本概念分析方法和提高稳定性的措施课程的基本要求1 建立起由短路引起的电力系统电磁暂态过程及各种扰动引起的机电暂态过程的完整概念2 关于电力系统电磁暂态过程应掌握同步发电机短路物理过程同步发电机参数数学模型发电机三相短路无限大系统短路各种不对称短路的计算方法以及励磁调节对短路的影响3 关于机电暂态过程部分应了解电力系统基本元件模型建立的方法掌握基本元件模型静态稳定暂态稳定电压稳定的基本概念及基本分析方法和提高稳定性的措施课程的基本内容本课程包含两部分电磁暂态部分1 概述短路成因类型危害标幺值2 无限大容量系统短路计算3 同步发电机三相短路基本假设及物理过程暂态次暂态电势暂态次暂态电抗的概念4 同步发电机数学模型稳态下同步发电机方程和相量图5 同步发电机空载及负载下三相短路等值电路时间常数及短路电流计算6 励磁调节对短路的影响7 三相短路实用计算转移阻抗异步机对短路电流影响8 对称分量法9 发电机变压器及输电线的各序参数计算及等值电路10 各种不对称短路计算方法正序等效定则非故障点电流电压计算11 断线计算电力系统稳定部分1 概述系统稳定性及分类2 同步发电机转子运动方程同步发电机功角特性方程电势方程励磁调节器的结构模型异步发电机转矩及运动方程3电力系统静态稳定物理过程电压稳定性的概念4 小扰动法分析电力系统静态稳定阻尼对静态稳定的影响5 励磁调节对静态稳定的影响6 多机系统静态稳定分析7 各种提高静态稳定性的措施8 电力系统暂态稳定物理过程分析等面积定则9 分段法和改进欧拉法求解系统运动方程10 励磁调节对暂态稳定影响11 多机系统暂态稳定分析12各种提高暂态稳定性的措施考研论坛-考研加油站/。
电力系统暂态上机计算课程设计报告
-0.5489+ 0.3169i
5
0.6091
0
0
0.5977+ 0.3451i
0.0000 - 0.6901i
-0.5977+ 0.3451i
c.各支路三序电流、三相电流：
支路
I1
I2
I0
Ia
Ib
Ic
1-2
0 + 0.8451i
0
0
0 + 0.8451i
0.7319- 0.4225i
3.《电力系统故障的计算机辅助分析》 大学 米麟书等
4.《电力系统潮流计算》 大学 宋文南等
5.《电力系统故障分析》 清华大学 周荣光
6.《短路电流实用计算方法》 电力工业 交通大学等
7.《精通Matlab6.5》航空航天大学，志涌
附录（设计流程图、计算结果、思考题答案）
1、对称短路计算流程图及计算结果
-0.7319- 0.4225i
1-3
0 - 4.5070i
0
0
0 - 4.5070i
-3.9032+ 2.2535i
3.9032+ 2.2535i
2-3
0 - 5.3521i
0
0
0 - 5.3521i
-4.6351+ 2.6761i
4.6351+ 2.6761i
1-4
0 + 3.6620i
0
0
0+3.6620i
另外，通过此次课程设计，我学习并掌握了一定的Matlab知识，对简单的选择语句、循环语句、数组的使用等有了一些了解，为今后学习、工作中再次学习使用此软件打下了基础，获益匪浅。
电力系统暂态上机暂态稳定
%归算后的发电机惯性时间常数%正常运行时发电机向无穷大系统传输的有功功率 %故障存在时发电机的最大功率%故障切除后发电机的最大功率%下面是利用改进欧拉法进行逐段计算 %需要注意故障切除前后电磁功率有跃变h=0.01;%设置步长 0.01s Duration=20 %设置计算时段长度 20sCutTime=input(' 输入故障切除时间 \n');%故障发生时的功角变化过程Delta(1)=33.92;%初始功角 Omega(1)=1; % 初始转速t(1)=0;for i=1:round(CutTime/h)d(i)=(Omega(i)-1)*360*f;o(i)=(PT-P2M*sin(Delta(i)*pi/180))/Tj;% 时段初变化率 Delta(i+1)=Delta(i)+d(i)*h;Omega(i+1)=Omega(i)+o(i)*h;% 时段末估计植dd(i)=(Omega(i+1)-1)*360*f;oo(i)=(PT-P2M*sin(Delta(i+1)*pi/180))/Tj;% 时段末变化率 ddd(i)=(d(i)+dd(i))/2;ooo(i)=(o(i)+oo(i))/2;% 平均变化率Delta(i+1)=Delta(i)+ddd(i)*h;Omega(i+1)=Omega(i)+ooo(i)*h;% 时段末计算值 t(i+1)=i*h;end;%显示此时功角 s=sprintf('Delta=%f\n',Delta(i+1)); disp(s);%故障切除后的功角变化过程for i=round(CutTime/h)+1:round(Duration/h)d(i)=(Omega(i)-1)*360*f;clear;clc;f=50; Tj=8.47; PT=1; P2M=0.48; P3M=1.38; %系统额定频率o(i)=(PT-P3M*sin(Delta(i)*pi/180))/Tj;% 时段初变化率Delta(i+1)=Delta(i)+d(i)*h;Omega(i+1)=Omega(i)+o(i)*h;% 时段末估计植dd(i)=(Omega(i+1)-1)*360*f;oo(i)=(PT-P3M*sin(Delta(i+1)*pi/180))/Tj;% 时段末变化率ddd(i)=(d(i)+dd(i))/2;ooo(i)=(o(i)+oo(i))/2;% 平均变化率Delta(i+1)=Delta(i)+ddd(i)*h;Omega(i+1)=Omega(i)+ooo(i)*h;% 时段末计算值t(i+1)=i*h;end;%画出摇摆曲线plot(t,Delta);。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
课程设计报告( 2013—2014年度第二学期)名称：电力系统暂态上机计算院系：电气与电子工程学院班级：电气1108学号：1111550112学生姓名：龙日尚指导教师：刘宝柱设计周数：两周成绩：日期：2014年7月11日一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。
二、设计正文（详细内容见附录,用A4纸，页数不限）1．对称短路计算过程流程图和计算结果2．不对称短路计算过程流程图和计算结果3．静态稳定性计算过程流程图和计算结果4．暂态稳定性计算过程流程图和计算结果5．思考题三、课程设计总结或结论本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。
通过本次课程设计，本人对电力系统故障分析有了更深刻的理解，电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。
在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。
在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。
同时我对用来分析电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。
这与手算系统短路电流时使用的网络化简方法大大不同。
另外，通过此次课程设计，我学习并掌握了一定的Matlab知识，对简单的选择语句、循环语句、数组的使用等有了一些了解，为今后学习、工作中再次学习使用此软件打下了基础，获益匪浅。
四、参考文献1.《电力系统暂态分析》，李光琦，中国电力出版社，2006年，第三版；2.《电力系统分析》（上、下），何仰赞，华中科技大学出版社，1998年，第二版。
3.《电力系统故障的计算机辅助分析》重庆大学出版社米麟书等4.《电力系统潮流计算》天津大学出版社宋文南等5.《电力系统故障分析》清华大学出版社周荣光6.《短路电流实用计算方法》电力工业出版社西安交通大学等7.《精通Matlab6.5》北京航空航天大学出版社，张志涌附录（设计流程图、计算结果、思考题答案）1.对称短路计算过程流程图和计算结果流程图计算结果：（1）导纳矩阵：Y1 =0 -26.6667i 0 +10.0000i 0 +10.0000i 0 +10.0000i 0 -33.3333i 0 +10.0000i 0 +10.0000i 0 +10.0000i 0 -20.0000i ⎛⎫ ⎪ ⎪ ⎪⎝⎭Y2 =0 -26.6667i 0 +10.0000i 0 +10.0000i 0 +10.0000i 0 -33.3333i 0 +10.0000i 0 +10.0000i 0 +10.0000i 0 -20.0000i ⎛⎫ ⎪⎪ ⎪⎝⎭Y0 = 0 -29.0476i 0 + 5.0000i 0 + 5.0000i 0 + 5.0000i 0 -48.0228i 0 + 5.0000i 0 + 5.0000i 0 + 5.0000i 0 -10.0000i ⎛⎫ ⎪ ⎪ ⎪⎝⎭YY1=0 -40.0000i 0 +10.0000i 0 +10.0000i 0 +20.0000i 0 0 +10.0000i 0 -60.0000i 0 +10.0000i 0 0 +40.0000i 0 +10.0000i 0 +10.0000i 0 -20.0000i 0 0 0 +20.0000i 0 0 0 -30.0000i 0 0 0 +40.0000i 0 0 0 -60.0000i ⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭YY2 = 0 -40.0000i 0 +10.0000i 0 +10.0000i 0 +20.0000i 0 0 +10.0000i 0 -60.0000i 0 +10.0000i 0 0 +40.0000i 0 +10.0000i 0 +10.0000i 0 -20.0000i 0 0 0 +20.0000i 0 0 0 -30.0000i 0 0 0 +40.0000i 0 0 0 -60.0000i ⎛⎫ ⎪ ⎪⎪ ⎪ ⎪ ⎪⎝⎭(2) 节点3发生三相短路故障a.故障点三序电流、三相电流：I1= 0 - 9.8592i I2=0 I3=0 ; Ia= 0 - 9.8592i Ib= -8.5383 + 4.9296i Ic= 8.5383 + 4.9296ib.各节点三序电压、三相电压：2．不对称短路计算过程流程图和计算结果计算结果：(1) 节点3发生A相短路接地故障a. 故障点三序电流、三相电流：I1= 0 - 3.1080i I2= 0 - 3.1080i I3= 0 0 - 3.1080i ; Ia= 0 - 9.3241i Ib= 0 - 0.0000i Ic= 0 - 0.0000i(2) 节点3发生A相经10Ω电阻接地故障a. 故障点三序电流、三相电流：I1= 0.0333 - 0.0004i I2= 0.0333 - 0.0004i I3= 0.0333 - 0.0004i ; Ia= 0.1000 - 0.0011i Ib= 0 Ic= 0(3) 节点3发生b、c两相短路故障a. 故障点三序电流、三相电流：I1= 0 - 4.9296i I2= 0 + 4.9296i I3= 0 ; Ia= 0 Ib= -8.5383 - 0.0000i Ic= 8.5383 + 0.0000i(4) 节点3发生b、c两相短路接地故障a. 故障点三序电流、三相电流：I1= 0 - 6.4036i I2=0 + 3.4556i I3= 0 + 2.9480i ; Ia= 0 Ib= -8.5383 + 4.4220i Ic= 8.5383 + 4.4220i3.静态稳定性计算过程流程图和计算结果流程图计算结果：Ke=0.100000, delta=91, P=1.345117 Ke=0.200000, delta=92, P=1.363000 Ke=0.300000, delta=93, P=1.380693 Ke=0.400000, delta=95, P=1.398409 Ke=0.500000, delta=96, P=1.415747 Ke=0.600000, delta=97, P=1.432842 Ke=0.700000, delta=98, P=1.449679 Ke=0.800000, delta=99, P=1.466244 Ke=0.900000, delta=100, P=1.482520 Ke=1.000000, delta=101, P=1.498487 Ke=1.100000, delta=101, P=1.513985 Ke=1.200000, delta=102, P=1.529380 Ke=1.300000, delta=103, P=1.544429 Ke=1.400000, delta=103, P=1.558878 Ke=1.500000, delta=104, P=1.573336 Ke=1.600000, delta=105, P=1.587402 Ke=1.700000, delta=105, P=1.600902 Ke=1.800000, delta=105, P=1.613924 Ke=1.900000, delta=102, P=1.621480 Ke=2.000000, delta=99, P=1.622734 Ke=2.100000, delta=97, P=1.623379 Ke=2.200000, delta=94, P=1.614349 Ke=2.300000, delta=92, P=1.607801 Ke=2.400000, delta=89, P=1.589402 Ke=2.500000, delta=87, P=1.576380 Ke=2.600000, delta=85, P=1.560855 Ke=2.700000, delta=83, P=1.542968 Ke=2.800000, delta=81, P=1.522859 Ke=2.900000, delta=80, P=1.514020 Ke=3.000000, delta=78, P=1.490598 Ke=3.100000, delta=76, P=1.465314 Ke=3.200000, delta=75, P=1.453440 Ke=3.300000, delta=73, P=1.425420 Ke=3.400000, delta=72, P=1.411962 Ke=3.500000, delta=70, P=1.381541 Ke=3.600000, delta=69, P=1.366704 Ke=3.700000, delta=68, P=1.351397 Ke=3.800000, delta=66, P=1.318171 Ke=3.900000, delta=65, P=1.301745 Ke=4.000000, delta=64, P=1.284942 Ke=4.100000, delta=63, P=1.267785 Ke=4.200000, delta=62, P=1.250299 Ke=4.300000, delta=61, P=1.232503 Ke=4.400000, delta=60, P=1.214420 Ke=4.500000, delta=59, P=1.196069 Ke=4.600000, delta=58, P=1.177469 Ke=4.700000, delta=57, P=1.158638 Ke=4.800000, delta=56, P=1.139594 Ke=4.900000, delta=55, P=1.120353 Ke=5.000000, delta=54, P=1.100930 Ke=5.100000, delta=53, P=1.081340 Ke=5.200000, delta=52, P=1.061599 Ke=5.300000, delta=51, P=1.041718 Ke=5.400000, delta=51, P=1.041792 Ke=5.500000, delta=50, P=1.021746 Ke=5.600000, delta=49, P=1.001589 Ke=5.600000, delta=49, P=1.001589最终选择放大倍数Ke=2.14．暂态稳定性计算过程流程图和计算结果流程图结果：0.15s时切除故障的摇摆曲线0.25s时切除故障的摇摆曲线根据摇摆曲线判断：0.15s时切除故障系统暂态稳定，0.25s时切除故障系统失稳。