东南大学电力系统课程设计报告

东南大学电气工程学院2015年电力系统综合课程设计班级******班姓名***学号******指导老师***小组成员*****1目录一、设计任务41、原始资料41.1网络拓扑41.2发电机稳态数据41.3变电站数据41.4线路数据51.5发电机暂态数据51.6计算方法52、设计内容63、设计成品6二、设计说明书61、有功平衡62、发电厂变电所主接线设计62.1发电厂主接线选择62.2发电厂主变选择62.3变电所主接线选择73、潮流计算74、调压措施的选择95、线路5-7断路器选择115.1断路器选择的原则115.2短路电流计算过程115.3选择结果126、仿真曲线的计算与仿真136.1理论基础与分析136.1.1简单模型下的暂态稳定分析136.1.2初值计算136.1.3负荷初值的计算比较146.1.4暂态稳定分析的网络数学模型及其求解方法146.1.5用改进欧拉法求解微分方程.156.2计算过程及结果.166.2.1初值计算1626.2.2故障系统与故障后系统描述.166.2.3等效简单网络电路模型图.56.2.4实验结果6三、体会5四、附录41、变电所主接线图52、PQ解耦法计算潮流源程序53、仿真曲线计算程序53一、设计任务1、原始资料：1.1网络拓扑27893负荷2G1G2T1T256负荷1负荷34T31G31.2发电机稳态数据名称额定电压机端电压典型方式输出有功G110KV1.05平衡机G210KV1.051.60G310KV1.050.901.3变电站数据负荷额定电压为10KV变电站编号负荷1负荷2负荷3最大负荷（100MVA）1.451.11.24典型方式负荷（100MVA）最小负荷（100MVA）Tmax（小时）功率因数cosφ低压侧电压（kV）对备用要求对调压要求负荷重要性1.4线路数据名称电压电阻7-8110KV0.018-9110KV0.027-5110KV0.035-4110KV0.019-6110KV0.026-4110KV0.03负序和正序相等、零序为正序的3倍1.5发电机暂态数据1.250.91.01.050.60.75000400045000.850.80.8510101030%40%30%逆调压顺调压常调压重要重要一般电抗对地导纳（B）0.0850.0440.0920.0400.1020.0760.1620.0890.1750.0370.1830.042名称D轴电抗D轴暂态电抗阻尼系数Tj（惯性时间常数）G10.140.072.047.2G20.900.221.912.8G21.320.132.16.021.6计算方法潮流方法：直角坐标牛拉法，极坐标牛拉法，PQ解耦法任选其一；仿真方法：改进牛拉法52、设计内容1．典型方式电力系统潮流计算(小组3人共同完成）2．变电站变压器选择和主接线设计（小组其中1人完成）3．线路5-7断路器选择（小组其中1人完成）4．仿真曲线计算（E’保持恒定，三相短路，功角曲线）（小组其中1人完成）3、设计成品1．整体装订（封面、目录、正文（包括任务书、设计过程和计算过程）、感想与建议）2.附录：源代码3.同一小组内每人写各自的报告。二、设计说明书1．根据变电所需要供电的情况，做出有功功率平衡。变电站编号负荷1负荷2负荷3最大负荷（100MVA）1.451.11.2功率因数cosφ0.850.80.85由得出，所需的最大负荷为??=1.45?0.85+1.1?0.8+1.2?0.85=3.1325.??取网损7%厂用电8%热备用5%得发电厂总装机容量为PG=3.1325*PL=3.759.所以，PG1=1.259，PG2=1.6，PG3=0.9.发电厂变电所主接线设计2.1发电厂主接线选择发电厂：高压侧为110kv，功率在100MVA-200MVA之间，选择双母线分段接线，可靠性高，灵活性也较高，发电机接变压器后直接连双母线分段接线较为合适。2.2发电厂主变选择SG1=1.4，SG1*0.7=0.98，选用型号SFP7-120000/110主变两台。SG2=1.778，SG2*0.7=1.24，选用型号SFP7-120000/110主变两台。SG3=1，SG3*0.7=0.7，选用型号SF7-75000/110主变两台。2.3变电所主接线选择由于出线端负荷较高，均用10kv出线，每个变电所分配四条出线。因为进线端只有双回进6线，故采用单母线分段接线。由于各个变电所功率相近，所以均按重要变电所设计。（接线图见附图）三个变电所均采用型号SFP7-120000/110主变两台。选用的变压器参数如下：选用地点选用型号空载损耗负载损耗空载电流短路阻/KW/KW(%)抗(%)发电厂1、2主变SFP7-120000/1101064220.510.5变电所1、2、3主变发电厂3主变SF7-75000/110753000.610.5变压器型号R(Ω)X(Ω)G(10-6S)B(10-6S)SFP7-120000/1100.35459710.58758.7649.6SF7-75000/1100.64516.946.19837.19选取??=100MVA,??=110Kv,????变压器参数标幺值变压器型号RXGBSFP7-120000/1102.93*10-30.08751.06*10-36*10-3SF7-75000/1105.33*10-30.147.5*10-44.5*10-33.潮流计算我们组采用PQ解耦算法计算潮流，画出等值网络计算导纳矩阵。将平衡节点设为节点1，因此原图中节点1、2互换。在计算过程中，按照此导纳矩阵计算，迭代多次不能收敛，因此，为了简化计算，将节点4、互换，形成导纳矩阵如下：7按照程序运行，迭代十次之后收敛，得到潮流结果如下：BusVoltageGenerationLoad#Mag(pu)Ang(deg)P(MW)Q(MVAr)P(MW)Q(MVAr)----------------------------------------------------11.0500.640590.0058.58--21.0500.000*19.0177.94--31.0500.6583160.0067.74--841.06090.5839----51.06260.4370--106.2565.8561.10020.6550--85.0052.6871.02800.2096----81.04410.6265--72.0054.0091.05080.5948------------------------------------Total:269.01204.27263.25172.534．调压措施的选择根据最大最小潮流，计算出负荷端最大负荷与最小负荷的电压值。变电所1：选择变压器分接头1）最大负荷时流入功率??=123.25+j76.38MVA，??=110.033Kv??????1?????????max=P1MAXRQ1MAXX123.25?0.177+76.38?5.29=110.033=3.87KVU1MAX‘=106.163Kv??2??????‘???106.163?11????1????????=2??????2??==111.218??????2??????10.52）最小负荷时流入功率??=89.25+j55.31MVA，??=110.88Kv??????1??????RQ1MINX=89.25?0.177+55.31?5.29???=2.801Kvmin=P1MIN110.88U1MIN‘=107.287Kv??2??????‘???2??107.287?11????1????????=2??????==118.0517????10??2??????=0.5*(??)=114.617Kv??1??1????????+??1????????选择分接头（1+2.5%*2）UN=115.5KV校验大、小负荷时低压侧实际电压：??2??????=106.163?11=10.111????<10.5Kv115.5107.287?11=10.218????>10Kv??2??????=115.5所以，所选择的分接头负荷要求。变电所2选择变压器分接头1）最大负荷时流入功率??=88+j66MVA，??=107.976Kvmax1max9???=P1MAXRQ1MAXX88?0.177+66?5.29==3.378KVmaxU1MAX107.976‘=104.598Kv??2??????‘104.598?11??=??2?????????2??==112.252????1??????????2??????10.252）最小负荷时流入功率??=48+j36MVA，??=108.493Kv??????1?????????min=P1MINRQ1MINX=48?0.177+36?5.29=1.834KvU1MIN108.493‘=106.659Kv??2??????‘???2??106.659?11????1????????=2??????==109.139????10.75??2??????=0.5*(??)=110.696Kv??1??1????????+??1????????选择分接头（1+2.5%*0）??=110KV??校验大、小负荷时低压侧实际电压：106.163?11??2??????==10.616????>10.25Kv110107.287?11??==10.729????<10.75Kv2??????110所以，所选择的分接头负荷要求。变电所3选择变压器分接头2）最大负荷时流入功率??=102+j63.21MVA，??=116.039Kvmax1max???=1MAXRQ1MAXX=102?0.177+63.21?5.29=3.037KVPmaxU1MAX116.039‘=113.002Kv??2??????‘113.002?11??=??2?????????2????=10.5=118.383????1????????2??????2）最小负荷时流入功率??=59.50+j36.87MVA，??=113.916Kv??????1?????????=P1MINRQ1MINX=59.5?0.177+36.87?5.29=1.805KvminU1MIN113.916‘=112.111Kv??2??????‘???2??112.111?11????=2??????==120.904????1??????????2??????10.210??=0.5*(??+??)=119.644Kv1??1????????1????????选择分接头（1+2.5%*2）??=115.5KV??校验大、小负荷时低压侧实际电压：??=113.002?11=10.762????>10.5Kv2??????115.5??2??????=112.111?11115.5=10.677????>10.5Kv不能符合调压要求分析：如果负荷变动较小，线路上的电压损耗也较小，这种情况一般将中枢点电压保持在102%UN到105%UN的范围内，不必随负荷变化来调整中枢点的电压，仍可保证负荷点的负荷质量，这种调压方式称为恒调压或者常调压。而对于变电所3来说，最大负荷为102+j63.21MVA，最小负荷为59.50+j36.87MVA，负荷波动为42.5+j26.34MVA，波动较大，常调压并不适用。5.线路5-7断路器选择5.1断路器选择的原则高压断路器的种类和型式的选择，除满足各项技术条件和环境条件外，还应考虑便于安装调试和运行维护，并经技术经济比较后才能确定。根据我国当前生产制造情况，电压为6~220kV的电网可选用少油断路器、真空断路器和SF6断路器；330~500kV电网一般采用SF6断路器。但近年来，真空断路器在35kV及以下电力系统中得到了广泛应用，有取代有断路器的趋势。SF6断路器也已在中压10~35kV系统发展，并在城乡电网建设和改造中获得了应用。采用封闭母线的大容量机组，当需要装设断路器时，应选用发电机专用断路器。（摘自《发电厂电气部分》第六章第四节）对5-7线路，其额定电压等级为110kV，按上述规则选择额定电压为110kV的少油式断路器。其正常工作时额定电流由潮流计算可知，标幺值为1.25，有名值为1.136kA。计算其短路电流时，取最大的短路电流即三相短路情况下的短路电流。5.2短路电流计算过程将7短路和5短路之后的等值电路作出如图所示。将电路化简后求得短路电流。处发生三相短路时：11计算得线路5-7上短路电流标幺值为8.814，即有名值为8.012kA。处发生三相短路时：计算得线路5-7上短路电流标幺值为6.693，有名值为6.085kA。5.3选择结果查询书后表格（如下图）择SW4-110Ⅲ型断路器。126.仿真曲线的计算与仿真6.1理论基础与分析6.1.1简单模型下的暂态稳定分析对于地区性的电力系统来说，一般失去暂态稳定的过程发展很快，通常分析系统遭受干扰后第一摇摆周期(l-1.5s)的机电暂态过程就可以判断系统是否能够维持稳定运行。这种情况下的暂态稳定分析中，由于调速系统的惯性，使得在短时间内原动机的功率不会发生很大变化，因此可以忽略调速系统的作用而假定原动机的功率维持不变此外由于发电机励磁绕组的时间常数较大，这样在短时间内其磁链也不会发生显著变化,而对励磁调节系统的作用,可以用发电机暂态电势Eq’保持恒定来近似模拟，即认为在第一摇摆周期内，励磁绕组中自由电流分量的变化由励磁调节系统的作用所补偿，从而使励磁绕组的磁链在这段时间内保持不变。相应地，阻尼绕组的影响将略去不计。对发电机、负荷及电力网络采用不同的数学模型可以构成各种简化的暂态稳定分析程序，采用何种组合要根据所研究问题的性质而走。为说明简化暂态稳定分析程序的原理和方法，下面介绍的简化暂态稳走分析程序，来用如下的数学模型相计算方法：发电机：用发电机暂态电势Eq’保持恒定来模拟。负荷：较小的负荷用恒定阻抗模拟，较大的负荷考虑综合负荷电动机转子机械暂态过程。电力网络：用导纳矩阵描述。微分方程：用改进欧拉法求解。整个暂态稳定计算的基本过程仍如上图所示。下面给出机电暂态过程的计算的实现。6.1.2初值计算在进行暂态稳定分析前，首先应根据潮流程序算出的干扰前系统运行状态，确定微分方程求解所需的初值。在简单模型下的暂态稳定分析中，初值的计算包括干扰前瞬间发电机的暂态电势、转子角度、原动饥的帆械功率以及综合负荷电动机的滑差，等值电纳等。这些参数在系统受到干扰后的瞬间是不会发生突变的。以下各变量的下标(o)表示初值。首先，我们介绍发电机的初值计算：由潮流计算可得到干扰前各发电机的端电压V(0)VX(0)Vy(0)和各发电机注入网络的功率S(0)P(0)Q(0)，进而可以计算出发电机注入网络的电流为：I(0)IX(0)jIy(0)(S(0))U(0)这样，就可求出虚构电势EQ(0)，即：EQ(0)EQX(0)jEQy(0)U(0)(RjXq)`I(0)以此就可以确定发电机的转子角度的初值：acrtg(EQy(0))EQx(0)在系统稳定运行时，发电机转子以同步转速旋转，于是有w(0)1利用坐标变化公式可以求出发电机定子电压和电流的d，q分量：13Vd(0)sinVq(0)cosId(0)sinIq(0)cos(0)(0)(0)(0)cossincossin(0)Vx(0)(0)Vy(0)(0)Ix(0)Iy(0)然后就可以算出暂态电势的值：Eq(0)'Vq(0)RaIq(0)Xd'Id(0)另外，稳态运行时发电机的电磁功率Pe(0)等于原动机的机械功率Pm(0)，即有Pm(0)Pe(0)=P(0)(I2X(0)Iy(0)2)Ra=6.1.3负荷初值的计算比较由潮流计算结果可知干扰前各负荷节点电压V(0)和负荷所吸收的功率S(0)，据此容易得到负荷的等值导纳:Y(0)S(0)2U(0)当按恒定阻抗模拟负荷时，其等值导纳在整个暂态过程中保持不变，如前所述，可以将它包括在网络导纳矩阵里。对于考虑结合负荷电动机机械暂态过程的负荷，由于在扰动瞬间电动机的滑辈不能突变，因此负荷的等值导纳也不应突变，即扰动后瞬间负荷的等值导纳与正常运行情况的等值导纳相同。6.1.4暂态稳定分析的网络数学模型及其求解方法和潮流、短路计算中一样，电力网络的节点电压方程可用相量表示成YV=I式中I、V分别为电力网络节最注入电施和节点电压组成的列向量；Y为节点导纳矩阵。上式所描述的网络方程在形式上为线性方程组，其中的导纳矩阵Y仅由电力网络的结构和参数所决定。在以后的计算中，经常把电力网络方程写成实数形式:Vx1Ix1G11B11G1iB1iG1nB1nVy1Iy1B11.G11B1i.G1iB1n.G1nGi1Bi1GiiBiiGinBinVxiIxiVyiIyiBi1.Gi1Bii.GiiBin.GinGn1Bn1GniBniGnnBnnVxnIxnVynIynBn1.Gn1Bni.GniBnn.Gnn=式中，n为电力网络的节点数;G，B分别表示网络导纳姐阵元素Y的实部和虚部;Ii，Ij和Vi，14Vj分别表示节点注入电施和节在电压的实部和虚部。6.1.5用改进欧拉法求解微分方程简单模型下的暂态稳定分析中，全系统的微分方程式包括各发电机的转子运动方程:diws(wi1)dtdwi1PeiD(wi1))dt(PmTJ假设电力系统的机电暂态过程己经计算到t时刻，现在我们围绕微分方程的求解方法来讨论t+t时刻系统运行状态的计算过程。新的时段计算前总是先判断系统在t时刻有无故障或操作发生。若无故障或操作发生，则直接以t时刻系统状态作为初值求解微分方程，否则需要计算出故障或操作后电力网络的运行参数，用它和t时刻的状态变量作为初值求解微分方程。用改进欧拉法求解微分方理的步骤如下：(1)由t时刻各发电机的儿们相各电动机的t计算出系统所有节点的电压Vx,Vy和发电机节点的注入电流Ix,Iy。求出状态量在t时刻的导数值，即：diws(wi1)dtdwi1(PPD(wi1))dtmeiTJ其中，各发电机的电磁功率Pe按下式22Pe(t)=(VX(t)IX(t)Vy(t)Iy(t))(Ix(t)Iy(t))Ra(2)然后求出t+t此时刻状态量的初值估计：[0]dit(tt)(t)dttw(t[0]t)w(t)dwittdt(3)类似于第(1)步，由各发电机的δVx,Vy和发电机节点(tt)计算出系统所有节在的电压的注入电流Ix,Iy。类似于第（2）步，应求出在t+δt时刻状态量初值估计的导数值。最后，求出各状态量在t+δt时刻的教值，即：15[0]tdi(tt)(t)2(dtw([t0]t)w(t)t(dwi2dtdidttt)dwittt)dt6.2计算过程及结果6.2.1初值计算。计算各负荷的等值并联导纳，结果如下：S1(0)1.285883j0.796945Y1(0)2负荷（节点5）：U1(0)Y2(0)S2(0)1.040117j0.6446282负荷（节点6）：U2(0)Y3(0)S3(0)0.825351j0.6190132负荷（节点8）：U3(0)计算各发电机的暂态电势Eq'和初始功角(0)及其输入机械功率Pm(0)，结果如下表所示，各发电机角速度的初值为w1(())w2(())w3(())1。在以下计算中，不计发电机凸极效应相当于在E'q=C模型中令X=Xd',即为发电机的经典模型。发电机Eq'0Pm(0)11.128113.01390.900021.10200.65890.190131.238226.69501.60006.2.2故障系统与故障后系统描述故障期间的电力网络相当于在7号母线处并联一条阻抗为0的接地主路，这时只要将正常情况下导纳矩阵中的对角元素Y77改为无穷大(可以用大数模拟，例如在实际计算中取为10^20，即可得到故障期间系统的节点导纳矩阵YF。故障后的电力网络是切除线路5-7后的情咀。由于线路5-7对导纳矩阵的贡献为：5716...1jb...-1jx...rrjb...-1...1jbrjb...rjxYl57因此故障后革统的节点导纳短阵YPYYL57。6.2.3等效简单网络电路模型图系统稳态时的导纳矩阵Y=[-j22.857100j22.857100000;0-j22.85710000j22.857100;00-j14.285700000j14.2587;j22.8571001.2520-j34.2625-0.3796+j6.1494-0.8724+j5.3215000;000-1.2520+j34.26253.0335-j15.09020-2.6539+j9.023400;000-0.8724+j5.321501.517-j10.9225000-0.6446+j5.06406;0j22.857100-2.6539+j9.023404.0191-j43.4246-1.3652+11.60410;0000000-1.3652+11.60413.6215-j21.9412-2.2563+j10.3791;00j14.258700-0.6446+j5.06410-2.2563+j10.37912.9010-j30.2669;]程序设计框图176.2.4实验结果1.在极限切除时间前切除182.在极限切除时间后切除19三、心得体会这次课程设计是对电力系统暂态和稳态知识的综合考量，考察电力系统主接线，主变选择，潮流计算，短路电流计算等。我们小组一共三人，分工合作，我负责潮流计算和变电所主变压器选择及主接线设计，刘程才负责仿真曲线的计算，刘宇航负责断路器的计算与选择。计算潮流的时候，起初我采用的是matpower计算，之后自己按照书上流程编写程序。开始的程序因为循环条件错误陷入死循环，和同学讨论之后修改正常。不过之后将导纳矩阵带入计算的时候始终不收敛，迭代次数显示200多次。于是我按照程序流程检查，最终确定程序无误，觉得可能是导纳矩阵的原因。因为4节点的功率为0，于是将4节点与8节点调换，相当于简化矩阵计算难度，结果迭代10次以后收敛。这次综合课程设计考察了电力系统的专业知识，将以前专业课的重点知识都概括了进去，比以往上课的时候更加灵活运用所学的知识，加深了我们对专业知识的理解与运用。四、附录201、变电所主接线设计图：212、PQ解耦源程序：%%导纳矩阵%%节点1、2调换，节点4、8调换Y=[-22.85710000022.857100;0-14.28570000014.28570;00-22.85710000022.8571;000-22.59230011.7647010.8696;0000-17.918109.80398.19670;00000-11.139305.46455.7143;22.85710011.76479.80390-44.377800;014.2857008.19675.46450-27.89040;0022.857110.869605.714300-39.4145];%%B'矩阵B1=[-14.28570000014.28570;0-22.85710000022.8571;00-22.59230011.7647010.8696;000-17.918109.80398.19670;0000-11.139305.46455.7143;0011.76479.80390-44.377800;14.2857008.19675.46450-27.89040;022.857110.869605.714300-39.4145];%%B''矩阵B11=[-22.59230011.7647010.8696;220-17.918109.80398.19670;00-11.139305.46455.7143;11.76479.80390-44.377800;08.19675.46450-27.89040;10.869605.714300-39.4145];%%设置电压初值，θ初值u=[1.051.051.05111111];xita=[000000000];p(2)=0.9;p(3)=1.6;p(4)=0.72;p(5)=1.0625;p(6)=0.85;p(7)=0;p(8)=0;p(9)=0;q(4)=0.54;q(5)=0.6585;q(6)=0.5268;q(7)=0;q(8)=0;q(9)=0;%标志设置k=0;%%%%%%%dpp=[0;0;0;0;0;0;0;0;0];dp=[0;0;0;0;0;0;0;0;0];dqq=[0;0;0;0;0;0;0;0;0];dq=[0;0;0;0;0;0;0;0;0];dxita=[0;0;0;0;0;0;0;0;0];%du=[0;0;0;0;0;0;0;0;0];kk=1;%%%%%%%%whilekk>0.0000001%%%%%%%%注入有功不平衡量%%%%fori=2:9forj=1:9dpp(j)=u(i)*u(j)*Y(i,j)*sin(xita(i)-xita(j));enddp(i)=p(i)-sum(dpp);end%%%%%%%θ修正量%%%%%%dxita=-inv(B1)*[dp(2)/u(2);dp(3)/u(3);dp(4)/u(4);dp(5)/u(5);dp(6)/u(6);dp(7)/u(7);dp(8)/u(8);dp(9)/u(9);];fori=1:8dxita(i)=dxita(i)/u(i+1);xita(i+1)=xita(i+1)+dxita(i);end%%%%%%%%注入无功不平衡量%%%%23fori=4:9forj=1:9dqq(j)=-u(i)*u(j)*Y(i,j)*cos(xita(i)-xita(j));enddq(i)=q(i)-sum(dqq);end%%%%%%%U修正量%%%%%%du=-inv(B11)*[dq(4)/u(4);dq(5)/u(5);dq(6)/u(6);dq(7)/u(7);dq(8)/u(8);dq(9)/u(9);];fori=1:6u(i+3)=u(i+3)+du(i);endkk1=max(abs(dxita));kk2=max(abs(du));kk=max(kk1,kk2);k=k+1;end3、仿真曲线计算程序：%Y_Modified_EulerS1=0.9000+0.5858i;U1=1.05;delt1=7.345;Xq1=0.13;S2=0.1901+0.7794i;U2=1.05;delt2=0;Xq2=0.07;S3=1.6000+0.6774i;U3=1.05;delt3=10.986;Xq3=0.22;S5=1.0625+0.6585i;U5=0.909;delt5=-2.776;S6=0.8500+0.5268i;U6=0.904;delt6=0.676;S8=0.7200+0.5400i;U8=0.934;delt8=1.235;%计算负载导纳初值[Y5,Y6,Y8]=init_loadY_cal(S5,U5,delt5,S6,U6,delt6,S8,U8,delt8);%计算发电机暂态电动势Eq_0，初始功角delt_0，输入机械功率Pm0[Eq1_0,delt1_0,Pm1_0]=init_Eq_cal(S1,U1,delt1,Xq1);[Eq2_0,delt2_0,Pm2_0]=init_Eq_cal(S2,U1,delt2,Xq2);[Eq3_0,delt3_0,Pm3_0]=init_Eq_cal(S3,U1,delt3,Xq3);%故障系统与故障后系统描述[Y0]=init_SystemY();[Y1]=Modified_SystemY(Y5,Y6,Y8);[Y2,Y3]=Bugsystem_Y(Y5,Y6,Y8);24%用改进欧拉法计算系统，以两机之间相对功角不超过180以及时间在2s以内为运算条件t=0;dt=0.005;Tm=10;delt21=delt2_0-delt1_0;delt31=delt3_0-delt1_0;delt1_1t=delt1_0;delt2_1t=delt2_0;delt3_1t=delt3_0;Delt1=[delt1_1t;0;0];Delt2=[delt2_1t;0;0;];Delt3=[delt3_1t;0;0;];Delt21=[delt21;0;0;];Delt31=[delt31;0;0;];w1_1t=1;w2_1t=1;w3_1t=1;c=2;%做保存用的计数值whilet<2ifabs(delt21)<180&&abs(delt31)<180&&t<=Tm-0.005[delt1_1,delt2_1,delt3_1,w1_1,w2_1,w3_1]=dall_delt_w(delt1_1t,delt2_1t,delt3_1t,w1_1t,w2_1t,w3_1t,Y1);delt1_1t=delt1_1;delt2_1t=delt2_1;delt3_1t=delt3_1;w1_1t=w1_1;w2_1t=w2_1;w3_1t=w3_1;delt21=delt2_1-delt1_1;delt31=delt3_1-delt1_1;Delt1(c)=delt1_1t;Delt2(c)=delt2_1t;Delt3(c)=delt3_1t;Delt21(c)=delt21;Delt31(c)=delt31;elset1=0:dt:dt*(size(Delt21)-1);plot(t1,Delt21);holdon;plot(t1,Delt31,'-.');xlabel('时间/s');ylabel('\delta21和\delta31');legend('\delta21','\delta31');title('发电机1,2,3相对功角与时间关系');break;endt=t+dt;c=c+1;endwhilet