电力系统暂态上机计算课程设计报告

1、课程设计报告( 20102011年度第二学期)名 称： 电力系统暂态上机 院 系： 电气与电子工程学院 班 级： 电气0803班 学 号： 1081140222 学生姓名： 张野 指导教师： 姜彤 设计周数： 两周 成 绩： 日期： 2011年 7月6一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。二、 主要内容1电磁暂态过程计算1.1 三相短路计算1.1.1 计算机计算三相短路流程图（程序见附录）输入数据形成网络节点导纳矩阵应用式（3-17）计

2、算if应用式（3-20）计算任一支路电流结束选择短路点f应用式（3-19）计算任一节点电压附：式3-17 if = uf|0| / ( zff + zf )1/( zff + zf ) 1/zff式3-19 u1= u1|0| +u1 = u1|0| z1fif 1 z1fifuf = uf|0| +uf = zfifun= un|0|+un = un|0| znfif 1 znfif式3-20 iij = (ui uj) / zij (ui uj ) / zij = (ui uj) yij1.1.2 计算教材电力系统暂态分析(第三版)p68例（3-2）例3-2： 下图所示一环形网络，已知各元

3、件参数为：发电机：g1g3 100mw,10.5kv,cosn=0.86，xd”=0.183。变压器：t1t3 120mva,115/10.5kv,us(%)=10.5。线路：三条线路完全相同，长50km,电抗0.44欧姆/试计算母线三相短路后时刻的各节点电压、各支路电流以及发电机端电压。电路图和等值电路图如下：故障分量网络：1.1.3节点发生三相短路接地计算结果得：（下面各数值均为标幺值）各节点的三相电压节点标号123abcabcabc各节点三相电压向量0.4507-0.2254- 0.3903i-0.2254+ 0.3903i0.5352-0.2676- 0.4635i-0.2676+ 0

4、.4635i000各节点三相电压幅值0.45070.45070.45070.53520.53520.5352000 各支路电流支路编号1、2支路2、3支路1、3支路abcabcabc各支路三相电流向量0.8451i0.7319- 0.4225i-0.731- 0.4225i-5.3521i-4.6351+ 2.6761i4.6351+ 2.6761i-4.5070i-3.9032+ 2.2535i3.9032+ 2.2535i各支路三相电流幅值0.84510.84510.84515.35215.35215.35214.5074.5074.507发电机的端电压发电机节点节点4节点5abcabc发

5、电机端点三相电压向量0.2625 + 0.5000i-0.0000- 1.0000i-0.2625+ 0.5000i0.3567 + 0.5000i-0.0000- 1.0000i-0.3567+ 0.5000i发电机端点三相电压幅值0.56471.00000.56470.61421.00000.61421.2.不对称短路计算过程流程图和计算结果1.2.1不对称短路计算过程流程图：输入数据潮流计算形成三个序网的导纳矩阵选择故障点k根据故障类型选择不同公式计算故障处各序电流、电压，最后合成得相电流、电压应用式（5-39）、式（5-58）计算任一节点各序电压，应用式（5-40）计算任一支路各序电流

6、结束附：式（5-39）ui(1)=ui|0|-zif(1)if(1) ui(2)=-zif(2)if(2) ui(0)=-zif(0)if(0) 式（5-40）iij(1) = ( ui(1) - uj(1) ) / zij(1) iij(2) = ( ui(2) - uj(2) ) / zij(2) iij(0) = ( ui(0) - uj(0) ) / zij(0)式（5-58）ui(1)=ui|0|-(ziq(1)-zik(1)i(1) ui(2)=-(ziq(2)-zik(2)i(2)ui(0)=-(ziq(0)-zik(0)i(0)1.2.2计算教材电力系统暂态分析(第三版)p12

7、7例（5-1）例5-1：在例3-2的系统中又已知三台发电机中性点均不接地；三台变压器均为ynd11接线（发电机侧为三角形）；经试验得三条输电线路的零序电抗均为0.20（以60mva为基准值）。要求计算节点分别发生a相接地短路、a相经10欧姆的电阻接地短路、bc两相相间短路、bc两相接地短路时各节点的电压、各支路电流以及发电机端电压。三个序网图：正序：负序： 零序：3.2.1 计算结果（标幺值） （1）a相接地短路各节点电压：节点标号节点1节点2节点3abcabcabc各节点三相电压向量0.5828-0.3838- 0.8660i-0.3838+ 0.8660i0.6938-0.3549- 0.

8、8660i-0.3549+ 0.8660i0.0000-0.5289- 0.8660i-0.5289+ 0.8660i各节点三相电压幅值0.58280.94730.94730.69380.93590.93590.00001.01481.0148各支路电流：支路编号1、2支路2、3支路1、3支路abcabcabc各支路三相电流向量0.0000 + 0.8286i0.0000+ 0.0079i-0.0000+ 0.0079i0.0000 - 5.2014i-0.0000- 0.0040i-0.0000- 0.0040i0.0000 - 4.3728i0.0000+ 0.0040i-0.0000+

9、0.0040i各支路三相电流幅值0.82860.00790.00795.20140.00400.00404.37280.00400.0040发电机端电压发电机节点节点4节点5abcabc发电机端点三相电压向量0.6707 + 0.5000i-0.0000- 1.0000i-0.6707+ 0.5000i0.7011 + 0.5000i-0.0000- 1.0000i-0.7011+ 0.5000i发电机端点三相电压幅值0.83651.00000.83650.86121.00000.8612（2）a相经10欧姆的电阻接地短路各节点电压：节点标号节点1节点2节点3abcabcabc各节点三相电压向

10、量0.6490- 0.1525i-0.4022- 0.8236i-0.4022+ 0.9085i0.7424- 0.1119i-0.3779-0.8130i-0.3779-0.9191i0.1587 - 0.3654i-0.5243- 0.8766i-0.5243+ 0.8555i各节点三相电压幅值0.66670.91650.99360.75080.89650.99370.39841.02141.0034各支路电流支路编号1、2支路2、3支路3、1支路abcabcabc各支路三相电流向量-0.3028+ 0.6971i-0.0029+ 0.0067i-0.0029+ 0.0067i1.9007

11、- 4.3758i0.0015- 0.0033i0.0015- 0.0033i1.5979- 3.6787i-0.0015+ 0.0033i-0.0015+ 0.0033i各支路三相电流幅值0.76000.00730.00734.77080.00360.00364.01080.00360.0036发电机端电压发电机节点节点4节点5abcabc发电机端点三相电压向量0.7017 + 0.4286i0.0000- 1.0000i-0.7017+ 0.5714i0.7273 + 0.4397i0.0000- 1.0000i-0.7273+ 0.5603i发电机端点三相电压幅值0.82221.0000

12、0.90490.84991.00000.9181（3）bc两相相间短路各节点电压节点标号节点1节点2节点3abcabcabc各节点三相电压向量1.0000-0.5000- 0.3903i-0.5000+ 0.3903i1.0000-0.5000- 0.4635i-0.5000+ 0.4635i1.0000-0.5000-0.5000各节点三相电压幅值1.00000.63430.63431.00000.68180.68181.0000-0.5000-0.5000各支路电流支路编号1、2支路2、3支路3、1支路abcabcabc各支路三相电流向量0.00000.7319+ 0.0000i-0.73

13、19-0.0000i0.0000- 0.0000i-4.6351+ 0.0000i4.6351+ 0.0000i0.0000 - 0.0000i-3.9032+ 0.0000i3.9032+ 0.0000i各支路三相电流幅值0.00000.73190.73194.77084.63514.63510.00003.90323.9032发电机端电压发电机节点节点4节点5abcabc发电机端点三相电压向量0.5643 + 0.5000i-0.0000- 1.0000i-0.5643+ 0.5000i0.6113 + 0.5000i-0.0000- 1.0000i-0.6113+ 0.5000i发电机端

14、点三相电压幅值0.75391.00000.75390.78981.00000.7898（4）bc两相接地短路各节点电压节点标号节点1节点2节点3abcabcabc各节点三相电压向量-0.8871-0.3537- 0.3903i-0.3537+ 0.3903i0.8589-0.4217- 0.4635i-0.4217+ 0.4635i1.02810.00000.0000各节点三相电压幅值0.88710.52680.52680.85890.62670.62671.02810.00000.0000各支路电流支路编号1、2支路2、3支路3、1支路abcabcabc各支路三相电流向量0.0000 - 0

15、.0077i0.7319 - 0.4065i-0.7319- 0.4065i0.0000+ 0.0039i-4.6351+ 2.5296i4.6351+ 2.5296i0.0000 - 0.0039i-3.9032+ 2.1231i3.9032+ 2.1231i各支路三相电流幅值0.00770.83720.83720.00395.28045.28040.00394.44334.4433发电机端电压发电机节点节点4节点5abcabc发电机端点三相电压向量0.4693 + 0.5000i-0.0000- 1.0000i-0.4693+ 0.5000i0.5312 + 0.5000i0.0000-

16、1.0000i0.5312 + 0.5000i发电机端点三相电压幅值0.68581.00000.68580.72951.00000.72952静态稳定计算2.1静态稳定性计算过程流程图：由0和eq0算各个运行变量eq0,ugd0,ugq0,ug0以及k1k6根据正常运行的参数算出kemax和kemin得到劳斯稳定判据预选ke（ke=0.1）逐步增加，从49到110变化算出eq以及pe计算eq,ugd,ugq,ug,k1k6算出kemax和kemin打印输出此时的功率p和功角ke是否小于5.6结束是否是否满足劳斯判据ke=ke+0.1 =49否是2.2 计算教材电力系统暂态分析(第三版)p191

17、例（7-2）中，采用试探法选择放大倍数例7-2：一简单电力系统参数为：xd=xq=0.982,xd=0.344,xe=0.504,xd=1.486, xd=0.48。已知正常运行时pe0=1，eq0=1.972, 0=49,u=1。若发电机装有按端电压偏差调节的比例式励磁调节器，采用试探法求不同励磁放大倍数情况下，系统的静态稳定性。2.3 计算结果：ke delta pke=0.100000, delta=91, p=1.345117ke=0.200000, delta=92, p=1.363000ke=0.300000, delta=93, p=1.380693ke=0.400000, de

18、lta=95, p=1.398409ke=0.500000, delta=96, p=1.415747ke=0.600000, delta=97, p=1.432842ke=0.700000, delta=98, p=1.449679ke=0.800000, delta=99, p=1.466244ke=0.900000, delta=100, p=1.482520ke=1.000000, delta=101, p=1.498487ke=1.100000, delta=101, p=1.513985ke=1.200000, delta=102, p=1.529380ke=1.300000, d

19、elta=103, p=1.544429ke=1.400000, delta=103, p=1.558878ke=1.500000, delta=104, p=1.573336ke=1.600000, delta=105, p=1.587402ke=1.700000, delta=105, p=1.600902ke=1.800000, delta=105, p=1.613924ke=1.900000, delta=102, p=1.621480ke=2.000000, delta=99, p=1.622734ke=2.100000, delta=97, p=1.623379ke=2.20000

20、0, delta=94, p=1.614349ke=2.300000, delta=92, p=1.607801ke=2.400000, delta=89, p=1.589402ke=2.500000, delta=87, p=1.576380ke=2.600000, delta=85, p=1.560855ke=2.700000, delta=83, p=1.542968ke=2.800000, delta=81, p=1.522859ke=2.900000, delta=80, p=1.514020将三者关系两两绘成图线：用试探法可知，随着励磁放大倍数ke的增大，功角delta先增大再减小

21、，功率极限p也是先增大再减小。分析数据可知：当ke=1.8时，系统可达静态稳定极限，即有功角deltasl105、p=1.613924。3.暂态稳定计算3.1暂态稳定性计算过程流程图：根据正常运行方式计算所需初值：计算e及n=0判断此时刻系统有无扰动，根据扰动性质选择功率特性曲线(短路时为，切除故障后为)打印n要求计算时段n=n+1结束是否3.2 已知：系统的额定频率f=50，归算后的发电机惯性时间常数tj=8.47，正常运行时发电机向无穷大系统传输的有功功率pt=1，发生故障时发电机的最大功率p2m=0.48，故障切除后发电机的最大功率 p3m=1.38。用改进欧拉法编程，分析切除时间为0.

22、15s和0.25s两种情况下摇摆曲线以及分析系统的暂态稳定性。3.3 计算结果对于不同的切除时间其结果是不一样的。（1） 切除时间为0.15秒时摇摆曲线：可以看出由于及时切除故障，系统会保持暂态稳定。（2） 切除时间为0.25秒时摇摆曲线：可以看出由于切除故障时间过晚，系统无法保持暂态稳定。4思考题4.1计算短路电流，书中给出的手算方法与计算机编程方法有何区别？答：手算方法不用形成导纳矩阵，通过星角变换化简，求得最后的等值阻抗，但是如果要求非故障点的电压或电流量还得倒回去求。计算机编程方法要形成导纳矩阵，主要是利用公式来求故障点和非故障点的量。4.2用计算机方法，进行电力系统潮流计算和短路电流

23、计算，有哪些区别？答：电力系统常规的潮流计算的任务是根据给定的电网结构，发电规划及负荷分布情况，求出整个电网的运行状态，如各母线上的电压（幅值及相角），网络中的功率分布以及功率损耗等，潮流计算的结果是电力系统稳定计算和故障分析的基础，是正常时候，系统里功率和电流的分布。而短路电流计算则是在系统发生故障时计算短路后系统里各点电压，各支路流过的电流。它要先算出正常时电压和电流，即潮流分布。再算故障分量，两者相加才是最后系统里各点电压，流过的电流。可以说潮流计算是短路电流计算的一部分。4.3如果交给你一个任务，请你用已有的短路电流计算软件计算某地区电网某点发生短路后的短路电流，你应该做哪些工作？（收

24、集哪些数据，如何整理，计算结果如何分析）答：首先要收集该地区电网的电气连接图，收集各元件的参数，有发电机、变压器和线路的正序、负序和零序的阻抗参数。将电气连接图画出电网的等值电路的图，将有名值参数装换成在同一基准下的标幺值。利用保护装置可以检测出是哪里出现何种故障，根据故障点位置以及故障类型可以计算出短路电流。看短路电流值的大小，如果较小，电网可以持续运行12个小时，但应发出报警信号，让工作人员去把故障切除；如果短路电流较小，应迅速切除故障，可以采用自动重合闸来减低故障带来的危害。4.4静态稳定的概念，暂态稳定的概念？答：静态稳定：是指电力系统受到小扰动后，不发生非周期性失步，自动恢复到初始运

25、行状态的能力。暂态稳定：是指电力系统受到大扰动后，各同步电机保持同步运行并过渡到新的或是恢复到原来稳定运行方式的能力。4.5提高静态稳定、暂态稳定的措施？答：提高静态稳定措施：采用自动调节励磁装置；减小元件的电抗：采用分裂导线、提高线路的额定电压等级、采用串联电容补偿；改善系统的结构和采用中间补偿设备。提高暂态稳定措施：故障的快速切除和自动重合闸装置的应用；提高发电机输出的电磁功率：对发电机施行强行励磁、电气制动、变压器中性点经小电阻接地；减少原动机输出的机械功率。三.课程设计总结或结论为期两周的课程设计使我获益良多：首先初步掌握了matlab的基本使用方法；之后对自己学过的稳态暂态知识有了新

26、的理解和感悟，并在整个过程中发现了很多自己对知识掌握的不足之处；另外，基本掌握了采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固了计算机编程的能力。对系统导纳矩阵和阻抗矩阵：首先一定要理解好自导纳、互导纳以及自阻抗、互阻抗的物理意义，以便其更好地服务于对短路电流计算；其次要掌握形成节点导纳矩阵的正确方法，对节点导纳矩阵求逆就可以形成阻抗矩阵；此外注意嵌套循环的逻辑顺序，尽量使循环清楚明朗便于执行。对于短路电流计算部分：首先要理解好用对称分量法计算不对称故障短路电流的基本原理；其次，要会画三序网图，以及不同故障类型系统复合序网图；最后，除了会计算故障点的各参数外，还

27、要会算非故障处电压、电流值。对称故障是不对称故障的一种特殊情况，可以一起讨论。在编程中药注意的是，计算发电机端电压时，对称分量经变压器后的相位变化。静态稳定部分：首先要熟练掌握并使用公式，了解公式的物理意义；其次要学会利用所学的知识，解决实际的静态稳定问题，了解影响静态问的相关因素，并对改善静态稳定性的措施有所了解；另外，还需理解试探法的精神，学会使用这一方法解决其他类似问题。暂态稳定部分：学会使用改进的欧拉法，通过迭代算出下时段的delta和omega，绘制摇摆曲线，并通过摇摆曲线能判断出在相应的切除时间下，系统能否保持暂态稳定，可求取极限切除时间；另外，需掌握影响暂态稳定的基本因素，并了解

28、改善暂态稳定的基本措施，学会解决实际问题。感谢老师的指导和同学的热心帮助，才得以解决一个又一个的编程问题、补足知识的缺陷，顺利完成这次课程设计。四.参考文献1. 电力系统暂态分析，李光琦，中国电力出版社，2006年，第三版2. 电力系统分析（上、下），何仰赞，华中科技大学出版社，1998年，第二版3. 电力系统故障的计算机辅助分析 重庆大学出版社 米麟书等 4. 电力系统潮流计算 天津大学出版社 宋文南等 5. 电力系统故障分析 清华大学出版社 周荣光6. 短路电流实用计算方法 电力工业出版社 西安交通大学等7. 精通matlab6.5北京航空航天大学出版社，张志涌附录：1.matrix.m，

29、网络原始数据以及节点导纳矩阵、节点阻抗矩阵的形成程序clc;%清空屏幕显示clear;%清空所有变量n1=3; %输入网络的节点数n2=5; %输入包括所有发电机节点的网络的节点数k=3; %输入短路点的节点号%输入支路各序阻抗，zz1表示支路正序阻抗，zz2表示支路负序阻抗，zz0表示支路零序阻抗zz1(1,1)=j*0.15;zz1(1,2)=j*0.1;zz1(1,3)=j*0.1;zz1(2,1)=j*0.1;zz1(2,2)=j*0.075;zz1(2,3)=j*0.1;zz1(3,1)=j*0.1;zz1(3,2)=j*0.1;zz1(3,3)=j*0.0;zz2=zz1;yyy1

30、=zeros(n1);yyy2=zeros(n1);yyy0=zeros(n1);for p=1:n1 for q=1:n1 if zz1(p,q)=0 yyy1(p,q)=0; else yyy1(p,q)=1/zz1(p,q); end endendyyy2=yyy1;zz0(1,1)=j*0.0525;zz0(1,2)=j*0.2;zz0(1,3)=j*0.2;zz0(2,1)=j*0.2;zz0(2,2)=j*0.00263;zz0(2,3)=j*0.2;zz0(3,1)=j*0.2;zz0(3,2)=j*0.2;zz0(3,3)=j*0.0;for p=1:n1 for q=1:n1

31、 if zz0(p,q)=0 yyy0(p,q)=0; else yyy0(p,q)=1/zz0(p,q); end endend%求节点导纳矩阵，y1表示不计发电机节点的正序网络节点导纳阵，y2表示不计发电机节点的负序网络节点导纳阵，y0表示不计发电机节点的零序网络节点导纳阵，y1=zeros(n1);for p=1:n1 for q=1:n1 if p=q for a=1:n1 y1(p,q)=y1(p,q)+yyy1(p,a);%? end else y1(p,q)=(-1)*yyy1(p,q); end endendy2=zeros(n1);y2=y1;y0=zeros(n1);for

32、 p=1:n1 for q=1:n1 if p=q for a=1:n1 y0(p,q)=y0(p,q)+yyy0(p,a);%? end else y0(p,q)=(-1)*yyy0(p,q); end endendy1y2y0%节点导纳矩阵参考结果，yy1表示计及发电机节点的正序网络节点导纳阵，yy2表示计及发电机节点的负序网络节点导纳阵%输入支路各序阻抗，zz1表示支路正序阻抗，zz2表示支路负序阻抗，zz0表示支路零序阻抗zz1(1,1)=0;zz1(1,2)=j*0.1;zz1(1,3)=j*0.1;zz1(1,4)=j*0.0525;zz1(1,5)=0;zz1(2,1)=j*0.

33、1;zz1(2,2)=0;zz1(2,3)=j*0.1;zz1(2,4)=0;zz1(2,5)=j*0.0263;zz1(3,1)=j*0.1;zz1(3,2)=j*0.1;zz1(3,3)=j*0.0;zz1(3,4)=0;zz1(3,5)=0;zz1(4,1)=zz1(1,4);zz1(4,2)=zz1(2,4);zz1(4,3)=zz1(3,4);zz1(4,4)=j*0.095;zz1(4,5)=0;zz1(5,1)=zz1(1,5);zz1(5,2)=zz1(2,5);zz1(5,3)=zz1(3,5);zz1(5,4)=0;zz1(5,5)=j*0.0475;zz2=zz1;yyy

34、1=zeros(n2);yyy2=zeros(n2);for p=1:n2 for q=1:n2 if zz1(p,q)=0 yyy1(p,q)=0; else yyy1(p,q)=1/zz1(p,q); end endendyyy2=yyy1;yy1=zeros(n2);for p=1:n2 for q=1:n2 if p=q for a=1:n2 yy1(p,q)=yy1(p,q)+yyy1(p,a); end else yy1(p,q)=(-1)*yyy1(p,q); end endendyy2=zeros(n2);yy2=yy1;yy1yy2%正序节点阻抗z1=eye(n1)/y1;%

35、 eye(n1)是3阶的单位矩阵,求节点导纳矩阵y1的逆阵,即节点阻抗矩阵z%负序节点阻抗z2=eye(n1)/y2;%零序节点阻抗z0=eye(n1)/y0;2.shortcircuit.m，短路电流计算的程序clc;%清空屏幕显示clear;%清空所有变量matrix %计算节点导纳矩阵、节点阻抗矩阵，形成全局变量fault=input(输入短路类型 ;(0)fault=0为三相短路;n(1)fault=1为a相接地短路;n(2)fault=2为a相经10欧姆电阻接地短路;n(3)fault=3为bc两相相间短路;n(4)faul=4为bc两相短路接地;nfault=)%第一部分：计算所有

36、节点的a,b,c三相电压%计根据故障类型选择不同的计算公式，计算故障点各序电流if fault=0 ik1=1/z1(k,k);%fault=0为三相短路 ik2=0; ik0=0;elseif fault=1 ik1=1/(z1(k,k)+z2(k,k)+z0(k,k);%fault=1为a相接地短路 ik2=ik1; ik0=ik1;elseif fault=2 zf=10*60/1152;%换成标幺值 ik1=1/(z1(k,k)+z2(k,k)+z0(k,k)+3*zf);%fault=2为a相经10欧姆电阻接地短路 ik2=ik1; ik0=ik1;elseif fault=3 ik

37、1=1/(z1(k,k)+z2(k,k);%fault=3为bc两相相间短路 ik2=(-1)*ik1; ik0=0; elseif fault=4 ik1=1/(z1(k,k)+z2(k,k)*z0(k,k)/(z2(k,k)+z0(k,k); %fault=4为bc两相短路接地 ik2=(-1)*ik1*z0(k,k)/(z0(k,k)+z2(k,k); ik0=(-1)*ik1*z2(k,k)/(z0(k,k)+z2(k,k);end;%计算各序电流的故障分量for p=1:n1 if p=k i1(p)=-1*ik1; i2(p)=-1*ik2; i0(p)=-1*ik0; else

38、i1(p)=0; i2(p)=0; i0(p)=0; end;end;%计算各点各序电压的故障分量uu1(:,k)=y1i1.; %“”表示左除 ab=inv(a)*b inv是求逆阵。 即ax=b。 i.表示转置uu2(:,k)=y2i2.; %uu(:,k)取uu2的第k列 p89uu0(:,k)=y0i0.;%假设故障前空载，各节点电压正序分量为1 for p=1:n1 u1(p)=1;end;%计算发生故障后的电压值:正常运行的分量+故障分量u1=u1+uu1(:,k); %计算所有节点正序电压u2=uu2(:,k); %计算所有节点负序电压u0=uu0(:,k); %计算所有节点零序

39、电压a=-0.5+j*sqrt(3)/2;t=1 1 1 % t为对称分量法的变换矩阵 a2 a 1 a a2 1;for p=1:n1 u=u1(p) u2(p) u0(p); p uabc=t*u. %uabc表示p节点的a,b,c三相电压向量 uuabc=abs(uabc) %uuabc表示p节点的a,b,c三相电压有效值end;%第二部分：计算支路电流for m=1:n1 for n=1:n1 if mn mn(1)=m; %mn是什么？ mn(2)=n; mn i1(m,n)=(u1(m)-u1(n)/zz1(m,n);%支路正序电流 i2(m,n)=(u2(m)-u2(n)/zz2

40、(m,n);%支路负序电流 i0(m,n)=(u0(m)-u0(n)/zz0(m,n);%支路零序电流 iabc=t*i1(m,n) i2(m,n) i0(m,n).;%由序分量合成相分量 iabc abs(iabc) end; end;end;%第三部分：计算发电机的端电压，只写了部分程序，余下的工作由同学们完成。要考虑变压器绕组是11点接线的for p=1:n2 if p=k ii(p)=-ik1; else ii(p)=0; end;end;vv1(:,k)=yy1ii.;vv2(:,k)=yy2ii.; %假设故障前空载，各节点电压正序分量为1 for p=1:n2 v1(p)=1;e

41、nd;vv0=zeros(n2);%计算故障点发生故障后的电压值:正常运行的分量+故障分量v1=v1+vv1(:,k); %计算所有节点正序电压v2=vv2(:,k); %计算所有节点负序电压v0=vv0(:,k); %计算所有节点零序电压 for p=1:n2 if p=4|p=5 v=v1(p)*(exp(j*pi/6) v2(p)*(exp(-j*pi/6) v0(p) ;%对称分量经变压器后的相位变化 p vabc=t*v. %vabc表示p节点的a,b,c三相电压向量 vvabc=abs(vabc) %vvabc表示p节点的a,b,c三相电压有效值 end;end;3.stabili

42、ty_smallsignal.m，静态稳定计算的程序clear;clc;u=1; %系统电压eq0=1.972; %空载电动势ug0=1.21; %机端电压xe=0.504; %线路电抗值xd=0.982; %同步电抗xdd=0.344; %暂态电抗xd_all=1.486; %系统电抗(包含同步电抗)xdd_all=0.848; %系统电抗(包含暂态电抗)ke=0.1;p_result=1:56;k=0.1:0.1:5.6;delta_result=1:56;i=1;q=1;while(ke0)&(ke=kemin)&(ke=kemax) ke_result=ke; delta_result1=delta; p_result1=eq*u*sin(delta*pi/180)/xd_all; else s=sprintf(ke=%f, delta=%d, p=%fn,ke_result,delta_result1,p_result1); p_result(i)=p_result1; delta_result(i)=delta_result1; i=i+1; disp(s); break; end; end; k(q)=ke; q=q+1; ke=ke+0.1; end;subplot(2