单机无穷大算例系统说明

1、附录一算例系统I1.算例系统的数学模型算例系统采用的是文献1第12章所使用的单机无穷大系统，图fl-1为其系统单线图。系统基准频率是60Hz。下面分别介绍潮流计算和发电机初始状态计算，全部计算基于标幺值。j0.150.9950图fl-1单机无穷大系统单线图1.1. 潮流计算已知发电机机端电压幅值为Et=1.0,无穷大母线电压Eb=0.995/00发电机有功出力B?=0.9,无功出力Q?=0.3。设S为发电机机端电压相角,0X为发电机端口到无穷大母线之间的电抗之和，则根据下列公式：c?.P?=7Sin5XE0可得到8=sin-1(P-X)=sin-1(0.9x0.65=36?。0?1义0.995

2、1.2. 发电机初始状态计算：发电机参数如下表所示:表f-1发电机参数表参数（不计饱和效应）数值直轴同步电抗Xd1.81直轴暂态电抗Xd0.3直轴次暂态电抗Xd0.23交轴同步电抗Xq1.76交轴暂态电抗Xaq0.65父轴次暂态电抗Xq0.25，一，一，、一，一直轴开路暂态时间常数Td08.0直轴开路次暂态时间常数Td00.03交轴开路暂态时间常数Ta0qU1.0交轴开路次暂态时间常数Tq00.07惯性时间常数H3.5s电枢电阻Ra0.003漏抗Xl0.16极对数p2由潮流结果可知，发电机定子电流?.（P+JQ）?It一。？计算得？=0.9-J0.3=0.949Z17.571/-36?设8为发

3、电机q轴相对于无穷大母线电压的角度，机端电压、电流与发?电机内电势的关系，如图fl-2所示。图fl-2同步电机的向量图Eq是发电机等值电路中一个虚拟的计算用的电势Eq=Et+(Ra+jXq渝计算得Eq=1/36?+(0.003+j1,76)X0.949/17.57?=2.204Z81.94也就是说3=81。4机端电压Et的直轴分量和交轴分量:ed=E?sin(1-,)=1Xsin(81.94-36)=0.718eq=E?qos(Sq-1)=1Xcos(81.94-36)=0.696定子绕组出口电流R直轴分量和交轴分量id=I?sin(禽-,+?)=0.949sin(81.94-36+18.49

4、)=0.856iq=I?cos(S-S+?)=0.949cos(81.94-36+18.49)=0.411暂态电势的计算公式为?=Et+(Ra+jXd)?彳#?=1/36?+(0.003+j0.3)X0.949/17.57?=1.125Z49.84不计发电机的饱和效应，空载电势？?的计算?=Eq+I?(Xd-Xq)彳#?=2.204+0.856X(1.81-1.76)=2.24681.3. 发电机动态模型发电机转子运动方程?1?(T?-t?-k?A?)?=(?-1)?R?(,0其中T?-标幺机械转矩T?-标幺电气转矩?-机械阻尼转矩系数?-转子角速度P?-原动机功率已7-电磁功率?-转子相对于

5、同步旋转参考轴的角位移，单位为电气弧度。?同步转速，？3=2%f0=100兀/秒发电机经典模型.一.一、.、._.忽略暂态凸极效应，也就是Xd=X?7在暂态过程中，q轴阻尼绕组与励磁绕组磁链保持不变，于是？?保持不变。定子电压方程?=?-I(R?+j?发电机三绕组模型忽略定子的电磁暂态，而考虑发电机转子阻尼绕组作用的三绕组发电机模型，也就是考虑到了f绕组、D绕组、Q绕组的电磁暂态和转子运动的机电暂态的发电机模型。发电机定子电压方程2?=?-?=?？眼-?转子f,D,Q绕组电势方程分别如下:dEqT=Efq-Eq+1d（Xd-Xd）”dEq”dEqTd0/=Eq-Id（Xd-Xd）+Eq+Td0

6、端”dEd，Tq0端二-Ed+lq（Xq-Xq）其中Eq-暂态电势Efq-由励磁电压？?所决定的假想空载电势Ed-直轴次暂态电势Eq交轴次暂态电势2 .算例系统的Matlab仿真模型在Matlab环境中，从Simulink和SimPowerSystems中，选取所需元件模块，分别建立上述？?恒定的系统模型，其仿真模型如图f-1中所示。主要可以分为以下几部分模块。同步发电机模块三相输电线路模块无穷大节点的电压源模块负荷模块故障模块测量模块为了计算方便，取模型系统的额定功率？%=1MVA,额定线电压（有效值）?=v3?1000Vo这与文献1的第12章中单机无穷大系统的参数（额定功率?=2220MV

7、A,额定电压？?=24?V不同，但是两个系统的标幺值是一致的，故能够保证分析结果的一致性。2.1. 步发电机模块从SimPowerSystems的“SimPowerSystems-Machines-Synchronous路径下，分别选取经典模型的发电机模块和三绕组模型的发电机模块，用于搭建?恒定的系统模型。经典发电机模型经典发电机模型采用SimplifiedSynchronousMachine模块，如图fl-3所示。Simplifiedchronons.MachinepuLmtt图fl-3SimplifiedSynchronousMachine模块按照fl-1设置相关参数，如下图所示:图fl-

8、4经典发电机模型参数阻尼系数(dampingfactor)设为0.1的原因是：图fl-1中的单机无穷大系统中的发电机没有阻尼，而Simulink仿真系统都是按照物理元件的实际情况进行设置的，现实中的发电机都是有阻尼的。所以设置阻尼系数为0.1,既保证和图fl-1中的单机无穷大系统中的发电机近似，又能保证仿真模型能够在这个阻尼的作用下，一段时间后能够达到稳态。其中“初始条件”(Initialcondition)参数中的th(?的角度)、ia、ib、ic(发电机出口电流)，详见1.2节的计算结果。三绕组发电机模型三绕组发电机模型采用SynchronousMachine模块，如图fl-5所示Sync

9、hronous.MachinepciStandard图fl-5SynchronousMachine模块其参数设置如下图所示:图fl-6三绕组发电机模块参数三绕组发电机模块参数和经典发电机模块参数类似。2.2. 相输电线路模块图fl-7三相输电线路模块仿真系统中的线路模型采用集中参数模型。不计输电线路对地导纳和线路电阻。图fl-1中发电机端口到无穷大母线之间的电抗之和为XE=j0.5+j0.15=j0.65(pu)则线路电感的有名值为?=(?)2?-=?.Ci/(v3kV)2J0.65*1MVA=J0.65X32?x60=2?x60(H)模型参数设置如下图所示:图fl-8三相输电线路模块参数2.

10、3. 无穷大系统的仿真从SimPowerSystems的“SimPowerSystems-ElectricalSource-ThreePhaseSource”路径下，选取一个三相电压源，将其视在功率设置为100MVA远大于发电机的视在功率1MVA,故可以将其看做是一个无穷大电源。具模型如下图所示。与三派窜矣图fl-9电压源模块其端口电压的有效值=0.995XVN=0.995X西(kV)x/R=10三相电压源模块的参数设置如下图所示：图fl-10电压源模块2.4. 负荷模块1H-三季总1图fl-11负荷模块SimPowerSystems元件库中的电感元件（如三相输电线路或者是变压器）不能和电流源

11、或者是被认为是电流源的非线性元件（如发电机）直接相连，因此Matalab模型中在线路的两侧添加了负荷模块，以满足仿真环境的要求。如果没有负荷模块1,仿真系统在启动仿真时会报错，所以要在发电机和输电线路之间安放一个负荷模块。为了和fl-1中的单机无穷大系统在发电机和输电线路之间没有负荷模块，保证两者最大程度的近似，其有功功率（Activepower）设置为一个很小的数值，这里取0.001MW。图fl-12负荷模块1的参数图fl-13负荷模块2的参数2.5. 故障模块在发电机机端加三相短路故障模块，用于模拟发电机机端短路的情况。图fl-13短路故障模块故障类型设置为三相短路故障，接地电阻取0.00

12、1。设置故障起始时间200s,故障切除时间(200+E)s,这里x是故障的持续的周波数，系统承受短路故障的时间为(A)s。通过更改x的值，就可以控制故障的大小图fl-14电路模块参数2.6. 测量模块测量模块能够在模型系统仿真时，将各参量的实时数据曲线清晰直观的展现出来，并可以将数据反馈到Matlab中的Workspace中，供进一步分析或是绘制图表。功率测量模块如下图所示。VabcPQ图fl-15功率测量模块2.7. 单机无穷大系统模型组合前文中的各个模块，在Simulink仿真窗口中搭建图fl-1所示的单机无穷大模型。采用经典发电机模块所搭仿真系统如下:采用三绕组发电机模块所搭仿真系统和采

13、用经典发电机模块所搭仿真系统相近，只是发电机模块不同。2.8. 型的验证为验证所搭的仿真系统是否准确，将仿真过程中各物理量的稳态值和1.2节计算得出的稳态值相对照。0.949Z17.57表f-2仿真稳态结果对照表1.2节计算出的数据经典发电机模型系统三绕组发电机模型系统发电机机端有功功率P?发电机机端无功功率Qt发电机机端电压E?10.9670.959E的相角st49.8449.8149.7定子绕组电流？0.9490.9480.959*无穷大系统母线电压Eb0.9950.99650.9971对照可知，无论是发电机经典模型系统还是发电机三绕组模型系统，各物理量

14、稳态值与1.2节计算出的系统稳态值相近，且在容许的误差范围内，故可将。c准确性。已知发电机机端电压幅值为Et=1.0,无穷大母线电压Eb=0.995/00发电机有功出力P?=0.9,无功出力Q?=0.3。设S为发电机机端电压相角,0X为发电机端口到无穷大母线之间的电抗之和，则根据下列公式：c?3.P?=sin8?XE0可得到8=sin-1(P-X)=sin-1(0.9x0.65=36?。0?1义0.9953 .算例系统仿真在进行仿真前，要设置仿真时采用的步长算法。由于模型是带有发电机的“刚性系统”，所以选择ode-23tb算法。ode-23tb算法是在龙格库塔的第一阶段使用梯形法，第二阶段用二

15、阶的BackwardDifferentiationFormulas算法，比ode23t算法和ode15s算法精度高。3.1. 不同扰动量情况下机电振荡变化在发电机经典模型仿真系统中，在发电机机端施加三相短路故障，观察系统在不同的扰动量(通过控制故障持续时间x来实现)的条件下，系统机电振荡响应曲线：发电机功角S曲线、转速3曲线、输出电磁功率？羯线的变化。像图fl-17那样记录响应曲线前三个峰值点的数据，数据如表f-2所示图fl-17扰动量为3个周波时的功角响应曲线表f-3不同扰动量的条件下机电振荡响应曲线数据扰动情况故障时间（周波）峰值时间/幅值随min/?7in?min/?7in?J?ax加m

16、in/?nin?max/?nin?min/?%in?min?r?ax?in0.39/0.92/1.43/0.67/1.16/1.68/0.228/0.87/1.38/97.0173.9693.510.9961.0030.997-0.380.251-0.13411(23.05(0.53)(0.51)8(0.49)6(0.6328(0.51)(0.006(0.52)(0.642)0.392/0.94/1.44/0.665/1.157/1.715/0.228/0.887/1.394/22101.471.9696.750.9961.0040.997-0.4020.285-0.154(29.44(0.5

17、48(0.008(0.492(0.558(0.687(0.66)(0.507*)(0.5)0.42/0.985/1.51/0.7/1.21/1.78/0.23/0.937/1.46/108.469.72101.70.9941.0040.996-0.4770.341-0.17933(39.08(0.565(0.5259(0.51)(0.57)(0.8195(0.70(0.523)(0.009)7)0.411/1.019/1.558/0.749/1.257/1.848/0.228/0.987/1.495/113.768.32105.30.9941.0050.996-0.4870.363-0.187

18、44(45.38(0.608(0.539(0.011(0.508(0.591(0.85)(0.759(0.508)0.457/1.09/1.64/0.815/1.323/1.93/0.211/1.053/1.561/119.767.63108.90.9931.0060.995-0.4930.337-0.18455(52.07(0.633(0.55)7(0.5085(0.83)(0.842(0.508)(0.012)(0.607)0.624/1.55/2.17/1.27/1.8/2.52/0.211/0.6/1.03/125.876.63114.20.9941.0050.996-0.4160.0

19、58-0.11866(49.17(0.926(0.62)6(0.53)5(0.474(0.389(0.43)(0.010(0.72)*)76.5失稳整理得到的数据，得到扰动大小和振荡周期频率之间的关系如下表所示:表f-4不同扰动量的条件下机电振荡响应曲线数据故障持续时间（周波）初始幅值振荡周期/频率123.051.04/0.96229.441.048/0.954339.081.09/0.917445.381.147/0.872552.071.183/0.845649.171.546/0.646分析表中数据可知：系统所加扰动量不同的情况下，随着扰动量的不断加大，系统所遭受的冲击越大，相应的故障响

20、应曲线峰值越大，振荡周期越长，非线性动态电力系统的初始运行点越远离系统的稳定运行点。当扰动量增大到一定程度时，有可能超出系统的稳定运行域，导致系统失稳。通过仿真发现当扰动量为6.5个周期时，系统失去稳定。功角5失稳转速3失稳输出电磁功率？?佚稳图fl-16扰动量为6.5个周波时系统失稳3.2. 相同扰动条件下频率响应的变化研究机电振荡中发电机功角S、转速3、输出电磁功率？勺振荡周期的变化有利于我们分析系统中对于振荡期间对系统稳定性起主要所用的模式。在三绕组发电机模型仿真系统中重复上述实验，扰动量x为6个周波时,仿真结果如下：表f-5扰动量为6个周波时机电振荡响应曲线数据3序号峰值时间/幅值半周

21、期时间半周期差值峰值时间/幅值半周期时间半周期差值峰值时间/幅值半周期时间半周期差值10.59/1341.13/0.00590.551/0.048321.42/73.230.831.67/-0.00510.540.894/-0.14660.34332.040.62-0.212.3660.6960.1561.4040.510.167*/114.3/0.004/0.286842.66/80.730.6202.909/-0.00340.543-0.1531.928/-0.12780.5240.01453.25/107.40.59-0.033.552/0.00270.6430.12.638/0.170

22、60.710.18663.85/85.190.600.014.105/-0.00230.553-0.093.195/-0.08540.557-0.15374.43/103.20.58-0.024.72/0.00190.6150.0623.821/0.11530.6260.06985.01/88.080.5805.278/-0.00160.558-0.0574.394/-0.05780.573-0.05395.59/100.40.5805.88/0.00130.6020.0444.99/0.08270.5960.023106.18/89.990.590.016.444/-0.0010.564-0

23、.0385.563/-0.0370.573-0.023*17116.75/98.480.57-0.027.037/0.00090.5930.0296.148/0.06060.5850.012127.34/91.230.590.037.60/-0.00070.563-0.036.735/-0.02320.5870.002137.91/97.110.57-0.028.187/0.00060.5870.0247.332/0.04620.5970.01148.49/92.030.580.018.756/-0.00050.569-0.0187.882/-0.01290.55-0.033159.06/96

24、.10.57-0.019.335/0.00050.5790.018.465/0.03620.5830.033169.63/92.520.5709.905/-0.00040.57-0.0099.0405/-0.00530.5755-0.00751710.20/95.330.57010.483/0.000320.5780.0089.615/0.02940.5745-0.001*1810.78/92.790.580.0111.055/-0.00030.575-0.00310.202/0.00010.5870.01251911.36/94.740.58011.628/0.00020.573-0.002

25、10.748/0.02510.546-0.0412011.94/92.910.58012.191/-0.00020.563-0.0111.307/0.00370.5590.013对表中数据进行分析可知:在系统所遭受到的扰动量相同时，由于系统中存在着正阻尼，随着时间的增加，使振荡能量不断的消耗，振荡幅值减小。同时，由于系统非线性因素的影响，振荡周期也相应减小。当扰动量为6个周波的时间长度时，能够得到相同的结论。这对于我们采用轨迹辨识的分析方法(比如Prony方法)，获取系统线性化的低阶模型，得到相应传递函数，并且据此进一步设计电力系统稳定器(PSS)提供帮助。4 .无穷大节点为电压源的建模建立一

26、个三相对称的无穷大电源，可以模拟电压的幅值、频率、相角层面的扰动。这个无穷大电源子系统的ABC三相输出电压分别为:*?=V)(1+?V)sin2兀f0(1+?f)t+?0x1000624-4-?=V0(1+?V)sin2%f0(1+?f)t+?0+3兀X1000，222一、-2,?=V0(1+?V)sin2%f0(1+?f)t+?0+3兀X1000，2其中V0-单机无穷大电源节点的电压幅值。?V-电压幅值扰动量？V=Ge,sin(2%ft+0)11807，、一一花?f频率扰动量？f=Qesin(2兀ft+0)2180rJ一，一花?0-相角扰动量？0=Ge。sin(2兀ft+03)(T为扰动的衰

27、减因子，f为扰动频率，C1C2G分别为？V、？f、？0的扰动初始幅值。0、0、0是扰动量的初始相角，单位为角。公式末尾乘以1231000v2的目的是为了将标幺值换算成有名值。基于上述含扰动的无穷大电源压源模型，在Matlab平台搭建的子系统如图fl-4所不*图fl-4含扰动的无穷大电源子系统下面分别介绍各个模块的功能。输入量模块：左侧为参量输入模块，生成含10个参量的u向量u=V0,Cl,C2,C3,01,02，03,f,*t即u1=Vo,u2=G,u3=C3u10=tFcn模块Fcn模块是函数表达式模块，可对输入量进行函数运算。Fcn1模块的函数为：u1*(1+u2*exp(u9)*sin(

28、2*pi*u8*u10+u5*pi/180)*1000*sqrt(2)是对输入向量u的各个元素u1、u2、u3u10进行函数计算的结果。即Fcn1模块输出的是九工?(1+C1e;rsin(2kt+%荷)X1000，2也就是V0(1+?V)X1000V2Fcn2模块的函数为:2*pi*60*u10*(1+u3*exp(u9)*sin(2*pi*u8*u10+u6*pi/180)兀2小t知sin(2ft+”-)也就是2兀f0(1+?f)-tFcn3模块中的函数为u4*exp(u9)*sin(2*pi*u8*u10+u7*pi/180)兀Gebsin(2兀ft+0工)331807也就是？0oFcn2

29、和Fcn3输出的结果在进行求和运算（在sum模块中实现）后，在Fcn4、Fcn5、Fcn6中加入A、B、C各相的初相角并计算相应的正弦值。即Fcn4中进行的运算为sin（Fcn2中的参数+Fcn3中的参数）4Fcn5中进仃的运算为sin（Fcn2中的参数+Fcn3中的参数+；兀）3Fcn6中进行的运算为sin（Fcn2中的参数+Fcn3中的参数+兀）3也就是Fcn4输出sin2/（1+?f）t+?04Fcn5输出sin2/（1+?f）t+?0+332Fcn6输出sin2/（1+?f）t+?0+-兀3.Productl、Product2、Product3模块分别输出的就是最终反馈至Simulin

30、k主系统的VA、VB、VC的值。Productl输出的结果=Fcn4输出的结果xFcnl输出的结果Product2输出的结果=Fcn5输出的结果XFcnl输出的结果Product3输出的结果=Fcn6输出的结果XFcnl输出的结果它们的具体表达式为：VA=Vo（1+?V）sin2兀f0（1+?f）-t+?04VB=V0（1+?V）sin2/（1+?f）t+?0+w兀32VC=Vo(1+?V)sin2兀f0(1+?f)-t+?0+-%31.1.不同扰动量情况下机电振荡变化在该模型中，改变无穷大系统内扰动源的扰动频率，观察系统机电振荡响应曲线：发电机功角s曲线、转速3曲线、输出电磁功率用曲线的变化

31、。当C1=o.1时，八3、2稳态峰值如下图所示。表f-5S、3、？?峰值表扰动频率f3(X10-4)P(X10-3)6maxSminSmax-6min3max3minmmax-3minPemaxPemin星max-gmin0.1084.80680.34874.45730.8711-0.60551.4766-0.2084.94179.73255.20852.055-1.16273.2177-0.4085.92578.61757.30754.9109-3.28578.19666.9881-0.01186.99990.6086.702677.35129.35147.696-8.947616.6436

32、14.9-27.342.20.7086.114578.26997.84467.0213-8.872815.894120.9-29.850.70.8085.828877.28718.54177.8114-11.70219.513426.5-59.886.3*0.9085.507574.88310.624510.428-16.29426.72248.1-99.2147.31.084.778871.75913.019814.807-21.59836.40574.6-140.42151.183.325665.174418.151224.854-31.37256.226119.1-206.3325.41

33、.1293.02450.742842.281264.995-72.339137.334279.8-420.3700.11.1497.37345.260352.112782.209-90.625172.834363.4-515.1878.51.1696.911244.525152.386184.264-92.236176.5384.4-522.7907.11.1796.26344.58951.67483.997-91.652175.649388.4-520.3908.71.1895.54144.77750.76483.148-90.436173.584390.9-515.6906.51.293.

34、846245.51548.331280.604-87.125167.729390.1-500.98911.385.202251.100134.102161.893-65.183127.076337.5-399.7737.21.575.985258.499917.485336.48-37.31773.797231.2-269.5500.7*1.775.522865.038710.484124.544-24.83149.375173.7-209.1382.81.975.961868.79587.16618.689-18.85637.545146-177.9313870.864

35、85.44915.322-15.40730.729132.2-161293.22.376.498772.16274.33613.111-13.11826.229124.4-149.2273.62.576.79472.95393.840111.696-12.01123.707120.7-143.9264.6相角扰动量的幅值（C3=0.02）时，、9、珞随扰动频率的变化如下图所示：振幅5232129171513121811171161112111119080706040201O扰动频率图fl-5AS变化曲线图fl-63变化曲线图fl-7?变化曲线当Ci=0.05时，S、3、忆稳态峰值如下图所示。表

36、f-5S、3、？?峰值表扰动频率f3（X10-4）pe（xio-3）SmaxSminSmax-8min3max3min3max-3minPmaxRminpmax-Rmin*0.183.136480.94760.3630.2669-0.283.195380.73340.8488-0.4804-0.483.796680.12892.2612-1.87832.3361-7.2560.584.246179.52653.4073-3.29914.966-11.30.684.190979.47183.9561-4.23037.7735-14.90.883.536380.2973.4435-3.87569.

37、9383-16.41.083.209880.64533.2935-3.701313-38480.65223.431-4.013814.8-22.41.283.160380.27133.9854-4.842218.1-33.81.382.964479.5785.0664-6.079623.3-48.91.581.93177.25698.581-9.580948.6-86.7*1.780.77576.11629.9273-10.53264.4-100.81.980.389176.54929.0465-9.374565.8-11977.0917.9844-

38、8.157863.9-922.380.149377.53157.0888-7.156462.6-88.42.580.085477.82826.3532-6.369461.1-84.11.1283.195880.61423.4898-4.179615.8-83.0883.691-4.49516.6780.482312-291.183.12280.442-4.57916.-28.7363.7541391.183.0733.821-4.67016.-30.8180.4014881无扰动在不同PSS情况下，系统的稳定情况表:KL=30KI=40KH=160KL=30等幅振扬增幅振

39、扬收敛慢KI=40增幅振扬等幅振扬收敛快KH=160收敛慢收敛快等幅振扬PSS系统的Bode图如下:图fl-11系统的Bode图KL=30、KI=40、KH=160时，系统的稳定性较好。只考虑KI、KH（也就是KL=0、KI=40、KH=160时），比同时考虑KL、KI、KH（也就是KL=30、KI=40、KH=160时）的稳定性要高。当减小模型中FL、FI、FH参数时，曲线左移。反之，右移参考文献1PRABHAKUNDUR,电力系统稳定与控制.北京:中国电力出版社，20022刘取.电力系统稳定性及发电机励磁控制.北京:中国电力出版社，20073西安交通大学等,电力系统计算,北京:中国电力出版

40、社，19784动态电力系统的理论和分析,电力系统计算,北京:中国电力出版社，1978扰动量为6个周波时，S、3、巳的响应曲线5响应曲线3响应曲线已响应曲线扰动量为6个周波时，仿真结果如下:序号3pe峰值时间/幅值半周期时间半周期差值峰值时间/幅值半周期时间半周期差值峰值时间/幅值半周期时间半周期差值10.59/1341.13/0.00590.551/0.048321.42/73.230.831.67/-0.00510.540.894/-0.14660.34332.04/114.30.62-0.212.366/0.0040.6960.1561.404/0.28680.510.16742.66/8

41、0.730.6202.909/-0.00340.543-0.1531.928/-0.12780.5240.01453.250.59-0.033.5520.6430.12.6380.710.186*/107.4/0.0027/0.170663.85/85.190.600.014.105/-0.00230.553-0.093.195/-0.08540.557-0.15374.43/103.20.58-0.024.72/0.00190.6150.0623.821/0.11530.6260.06985.01/88.080.5805.278/-0.00160.558-0.0574.394/-0.0578

42、0.573-0.05395.59/100.40.5805.88/0.00130.6020.0444.99/0.08270.5960.023106.18/89.990.590.016.444/-0.00110.564-0.0385.563/-0.03770.573-0.023116.75/98.480.57-0.027.037/0.00090.5930.0296.148/0.06060.5850.012127.34/91.230.590.037.60/-0.0000.563-0.036.735/-0.0230.5870.002*72137.91/97.110.57-0.028.187/0.000

43、60.5870.0247.332/0.04620.5970.01148.49/92.030.580.018.756/-0.00050.569-0.0187.882/-0.01290.55-0.033159.06/96.10.57-0.019.335/0.00050.5790.018.465/0.03620.5830.033169.63/92.520.5709.905/-0.00040.57-0.0099.0405/-0.00530.5755-0.00751710.20/95.330.57010.483/0.000320.5780.0089.615/0.02940.5745-0.0011810.78/92.790.580.0111.055/-0.00030.575-0.00310.202/0.00010.5870.01251911.36/94.740.58011.628/0.00020.573-0.00210.748/0.02510.546-0.041*2011.94/92.910.58012.191/-0.00020.563-0.0111.307/0.00370.5590.013扰动量为4个周波时，仿真结果如