单机无穷大算例系统说明1214

1、附录一算例系统I1 .算例系统的数学模型算例系统采用的是文献第12章所使用的单机无穷大系统，图fl-1为其系 统单线图。系统基准频率是60Hz。下面分别介绍潮流计算和发电机初始状态 计算，全部计算基于标幺值。图fll单机无穷大系统单线图1.1. 潮流计算已知发电机机端电压幅值为叁=1.0,无穷人母线电压备=0.995Z00, 发电机有功出力R = 0.9,无功出力Qt=O,3。设说为发电机机端电压相角,X2为发电机端口到无穷人母线之间的电抗之和，则根据卜列公式:_ E*b .Pt 可得到bo = siL （篙=siL （粽）=36 -1.2. 发电机初始状态计算:发电机参数如下表所示:表fl发

2、电机参数表参数(不计饱和效应)数值直轴同步电抗Xd1.81直轴暂态电抗X；0.3直轴次暂态电抗X；0.23交轴同步电抗Xq1.76交轴暂态电抗X；0.65交轴次暂态电抗X；0.25直轴开路暂态时间常数T；o8.0直轴开路次暂态时间常数Td00.03交轴开路暂态时间常数T；。1.0交轴开路次暂态时间常数Tq00.07惯性时间常数H3.5s电枢电阻Ra0.003漏抗X0.16极对数p2由潮流结果可知，发电机定子电流Y (P + jQ)* it = 耳计算得 = 0.949/17.57。IN36°设即为发电机q轴相对于无穷人母线电压的角度，机端电压、电流与发电机内电势的关系，如图fl-2所

3、示。fl-2同步电机的向量Eq是发电机等值电路中一个虚拟的计算用的电势EQ = Et +(Ra+jXq)/t计算得=2.204Z81.940 。Eq = 1N36。4- (0.003 4-jl.76) X 0.949Z17.570 也就是说,=81.94°机端电压国的直轴分量和交轴分量:ed = Etsin (6q-6e) = 1 X sin(81.94° - 36°) = 0.718eq = Et cos (8q 6e ) = 1 X cos (81.94° 36°) = 0.696定子绕组出口电流it直轴分量和交轴分量id = Itsin(

4、8q - 6e + 0)= 0.949 sin(81.94° - 36° + 18.49°) = 0.856iq = Itcos(8q - 5e + 0)= 0.949 cos(81.94° - 36° + 18.49°) = 0.411暂态电势的计算公式为彦=展 + (Ra + jXd)It得户' = 1N36。+ (0.003 + j0.3) X 0.949Z17.57* = 1.125Z49.84°不计发电机的饱和效应，空载电势0的计算Eq =EQ + l(Xd-Xq)得Eq = 2.204 + 0.856 X

5、 ( 1.81- 1.76 ) = 2.246813.发电机动态模型发电机转子运动方程dA 0)1k二瓦(以一一KM”)其中以-标幺机械转矩T标幺电气转矩kd一一机械阻尼转矩系数如一一转子角速度Pm-原动机功率Pe 电磁功率5转子相对于同步旋转参考轴的角位移，单位为电气弧度。oo0同步转速，= 27Tf0 = IOOtt/秒发电机经典模型忽略暂态凸极效应，也就是X； = X在暂态过程中，q轴阻尼绕组与励磁绕 组磁链保持不变，于是百'保持不变。定子电压方程Et = Ed I(R0+jXd)发电机三绕组模型忽略定子的电磁暂态，而考虑发电机转子阻尼绕组作用的三绕组发电机模 型，也就是考虑到了

6、f绕组、D绕组、Q绕组的电磁暂态和转子运动的机电暂态 的发电机模型。发电机定子电压方程= Xq1 -R/dq = Eq - Xdid Raiq转子f, D, Q绕组电势方程分别如下:,dE；,丁亚不 = Efq - Eq + 1d(Xd XQTd0 学=-Eq - Id(xd - xd) + E； + Td0 端TqO j = -Ed +- X；)其中E；暂态电势Efq一一-由励磁电压4所决定的假想空我电势E；一直轴次暂态电势E；交轴次暂态电势2.算例系统的Matlab仿真模型在Mat lab环境中，从Simulink和SimPowerSystems中，选取所需元件模块， 分别建立上述于恒定的

7、系统模型，其仿真模型如图f-1中所示。主要可以分为 以下几部分模块。 同步发电机模块 三相输电线路模块 无穷大节点的电压源模块 负荷模块 故障模块测量模块为了计算方便，取模型系统的额定功率Sn = 1MVA,额定线电压（有效值）= V3 * lOOOVo这与文献的第12章中单机无穷大系统的参数（额定功率Sn = 2220MVA额定电压匕=24AV）不同，但是两个系统的标幺值是一致的，故能够保证分析结果的一致性。1.1. 步发电机模块从 SimPowerSystems 的 “SimPowerSystems - Machines-Synchronous ” 路径 下，分别选取经典模型的发电机模块和

8、二绕组模型的发电机模块，用于搭建百'恒 定的系统模型。经典发电机模型经典发电机模型采用Simplified Synchronous Machine模块,如图fl-3 所示。Simplified Synchronous Machine pj Units图 fl-3Simplified Synchronous Machine 模块按照fl-1设置相关参数，如下图所示：图件4经典发电机模型参数阻尼系数（damping factor）设为0.1的原因是：图fl-1中的单机无穷大系统中 的发电机没有阻尼，而Simulink仿真系统都是按照物理元件的实际情况进行设 置的，现实中的发电机都是有阻尼的

9、。所以设置阻尼系数为0.1,既保证和图fl-1 中的单机无穷大系统中的发电机近似，又能保证仿真模型能够在这个阻尼的作 用下，一段时间后能够达到稳态。其中“初始条件"（Initial condition惨数中的th停 的角度）、ia、%、%（发电 机出口电流），详见1.2节的计算结果。三绕组发电机模型三绕组发电机模型采用Synchronous Machine模块，如图fl-5所示。Synchronous Mach ne pu Standardfl-5 Synchronous Machine 模块其参数设置如下图所示:图d6三绕组发电机模块参数三绕组发电机模块参数和经典发电机模块参数类似

10、。1.2. 相输电线路模块« C C R三相输电皴路图fl7三相输电线路模块仿真系统中的线路模型采用集中参数模型。不计输电线路对地导纳和线路电阻。图fl-1中发电机端口到无穷人母线之间的电抗之和为xz = J0.5+J0.15 = j0.65 (pu)则线路电感的有名值为0)j0.65 X(同 V)21MVA27r x 60j0.65 x 327r x 60(H)模型参数设置如卜图所示:件8三相输电线路模块参数1.3. 无穷大系统的仿真从 SimPowerSystems 的 u SimPowerSystems-Electrical Source-Three Phase Source”

11、路径下，选取一个三相电压源，将其视在功率设置为100MVA 远大于发电机的视在功率1MVA,故可以将其看做是一个无穷大电源。其模型 如下图所示。图fl-9电压源模块其端口电压的有效值=0.995 XVN = 0.995 xV3(kV)x/R=10三相电压源模块的参数设置如下图所示:图fl-10电压源模块1.4. 负荷模块力商模央1图fl-11负荷模块SimPowerSystems元件库中的电感兀件（如三相输电线路或者是变压器）不能和电流源或者是被认为是电流源的非线性元件(如发电机)直接相连，因此Matalab模型中在线路的两侧添加了负荷模块，以满足仿真环境的要求。如果没有负荷模块1,仿真系统在

12、启动仿真时会报错，所以要在发电机和 输电线路之间安放一个负荷模块。为了和fl-l中的单机无穷大系统在发电机和 输电线路之间没有负荷模块，保证两者最大程度的近似，其有功功率(Active power)设置为一个很小的数值，这里取O.OOlMWo图fl-12负荷模块1的参数fl-13负荷模块2的参数1.5. 故障模块在发电机机端加三相短路故障模块，用于模拟发电机机端短路的情况。fl-13短路故障模块故障类型设置为三相短路故障，接地电阻取0.001。设置故障起始时间200s,故障切除时间（200 +六）s,这里x是故障的持续的周波数，系统承受短路故障的时间为（1）s。通过更改x的值，就可以控制故60

13、障的大小。fl-14电路模块参数1.6. 测量模块测量模块能够在模型系统仿真时，将各参量的实时数据曲线清晰直观的展 现出来，并可以将数据反馈到Matlab中的Workspace中，供进一步分析或是绘 制图表。功率测量模块如下图所示。*Vabcn=n% > >r>->1-1labc=Gain2 SoopeZ图fl-15功率测量模块1.7. 单机无穷大系统模型组合前文中的各个模块，在Simulink仿真窗口中搭建图fl-1所示的单机无 穷大模型。采用经典发电机模块所搭仿真系统如下:采用三绕组发电机模块所搭仿真系统和采用经典发电机模块所搭仿真系统 相近，只是发电机模块不同。1

14、.8. 型的验证为验证所搭的仿真系统是否准确，将仿真过程中各物理量的稳态值和L2 节计算得出的稳态值相对照。0.949Z17.570表f-2仿真稳态结果对照表L2节计算出的数据经典发电机模型系统三绕组发电机模型系统发电机机端TT功功率R0.90.90.9发电机机端无功功率Q,0.30303发电机机端电压耳10.9670.959E的相角6(49.84"49.8149.7定子绕组电流40.9490.9480.959无穷大系统母线电压Eb0.9950.99650.9971对照可知，无论是发电机经典模型系统还是发电机三绕组模型系统，各物理量稳态值与1.2节计算出的系统稳态值相近，且在容许的误

15、差范围内，故可将。准确性。已知发电机机端电压幅值为耳=1.0,无穷大母线电压后b = 0.995Z0% 发电机有功出力用=0.9,无功出力。=0,3。设5。为发电机机端电压相角， Xx为发电机端口到无穷大母线之间的电抗之和，则根据下列公式：D EEb .Pt = -TT-sinSo可得到50 = sin-1 （也）=shr1 （） = 36°。uEtEB1X0.9953 .算例系统仿真在进行仿真前，要设置仿真时采用的步长算法。由于模型是带有发电机的“刚性系统”，所以选择。de.23tb算法。ode-23tb 算法是在龙格库塔的第一阶段使用梯形法，第二阶段用二阶的Backward Di

16、fferentiation Formulas 算法，比 ode23t 算法和 odel5s 算法精度高。3.1. 不同扰动量情况下机电振荡变化在发电机经典模型仿真系统中，在发电机机端施加三相短路故障，观察系 统在不同的扰动量（通过控制故隙持续时间x来实现）的条件下，系统机电振 荡响应曲线：发电机功角6曲线、转速3曲线、输出电磁功率2曲线的变化。 像图fl-17那样记录响应曲线前三个峰值点的数据，数据如表f-2所示。图件17扰动量为3个周波时的功角响应曲线表f3不同扰动量的条件下机电振荡响应曲线数据扰动 情况故障 时间 （周波）峰值时间/幅值uniin/ 6nlinSmiii/ 6min匕/鼠水

17、cumin/311mJ max/3 minJmin/3藤mPminniaxPndn11039/97.01(23.05)0.92/73.96(0.53)1.43/93.51(0.51)0.67/0.9968(0.006)1.16/1.003(0.49)1.68/0.9976(0.52)0.228/-038(0.632)0.87/0.2518(0.642)138/-0.134(0.51)220392/101.4(29.44)0.94/71.96(0.548)1.44/96.75(0.5)0.665/0.996(0.008)1.157/1.004(0.492)1.715/0.997(0.558)0.

18、228/-0.402(0.687)0.887/0.285(0.66)1394/-0.154(0.507)330.42/108.4(39.08)0.985/69.72(0.565)1.51/101.7(0.525)0.7/0.9949(0.009)1.21/1.004(0.51)1.78/0.996(0.57)0.23/-0.477(0.819)0.937/ 03415 (0.707)1.46/-0.179(0.523)440.411/113.7(4538)1.019/68.32(0.608)1.558/1053(0.539)0.749/0.994(0.011)1.257/1.005(0.508

19、)1.848/0.996(0.591)0.228/-0.487(0.85)0.987/0.363(0.759)1.495/0.187(0.508)550.457/119.71.09/67.631.G4/108.90.815/0.99371.323/1.0061.93/0.99550.211/-0.4931.053/0.3371.5G1/-0.184(52.07)(0.633)(0.55)(0.012)(0.508)(0.607)(0.83)(0.842)(0.508)0.624/1.55/2.17/1.27/1.8/2.52/0.211/0.6/1.03/66125.876.63114.20.

20、99461.0050.9965-0.4160.058-0.118(49.17)(0.926)(0.62)(0.010)(0.53)(0.72)(0.474)(0.389)(0.43)76.5失稔整理得到的数据，得到扰动大小和振荡周期频率之间的关系如下表所示:表f4不同扰动量的条件下机电振荡响应曲线数据故障持续时间（周波）初始幅值振藩周期/频率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分析表中数据可知:系统所加扰动量不同的情况下，随着扰动量的

21、不断加大，系统所遭受的冲击 越大，相应的故障响应曲线峰值越大，振荡周期越长，非线性动态电力系统的 初始运行点越远离系统的稳定运行点。当扰动量增大到一定程度时，有可能超 出系统的稳定运行域，导致系统失稳。通过仿真发现当扰动量为6.5个周期时， 系统失去稳定。功角b失稳转速3失稳输出电磁功率P©失稳图fl-16扰动量为6.5个周波时系统失稳3.2. 相同扰动条件下频率响应的变化研究机电振荡中发电机功角b、转速co、输出电磁功率丹的振荡周期的变 化有利于我们分析系统中对于振荡期间对系统稳定性起主要所用的模式。在三绕组发电机模型仿真系统中重复上述实验，扰动量x为6个周波时， 仿真结果如下：表

22、f5扰动量为6个周波时机电振荡响应曲线数据序号63峰值时间 /幅值半周期 时间半周期 基值峰值时间/幅值半周期 时间半周期 差值峰值时间 /幅值半周期 时间半周期 差值10.59/1341.13 /0.00590.551 /0.048321.42 /73.230.831.67 /-0.00510.540.894 /-0.1466034332.04 /11430.62-0.212366 /0.0040.6960.1561.404 /0.28680.510.16742.66 /80.730.6202.909 /-0.00340.543-0.1531.928 /-0.12780.5240.01453

23、.25 /107.40.59-0.033.552 /0.00270.6430.12.638 /0.17060.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.010.5805.2780.558-0.05743940.573-0.053/88.08/-0.0016/-0.057895.59/100.40.5805.88 /0.00130.6020.0444.9

24、9 /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.01212734/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.

25、030.580.018.756 /-0.00050.569-0.0187.882 /-0.01290.55-0.033159.06/96.10.57-0.019335 /0.00050.5790.018.465/0.03620.5830.033169.63792.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

26、.00310.202/0.00010.5870.0125191136/94.740.58011.628/0.00020.573-0.00210.748/0.02510.546-0.0412011.94/92.910.58012.191/-0.00020.563-0.0111307 /0.00370.5590.013对表中数据进行分析可知:在系统所遭受到的扰动量相同时，由于系统中存在着正阻尼，随着时间的 增加，使振荡能量不断的消耗，振荡幅值减小。同时，由于系统非线性因素的 影响，振荡周期也相应减小。当扰动量为6个周波的时间长度时，能够得到相同的结论。这对于我们采用轨迹辨识的分析方法(比如Pron

27、y方法)，获取系统线性化的 低阶模型，得到相应传递函数，并且据此进一步设计电力系统稳定器(PSS)提供 帮助。4 .无穷大节点为电压源的建模建立一个三相对称的无穷大电源，可以模拟电压的幅值、频率、相角层面的扰动。这个无穷大电源了系统的ABC三相输出电压分别为：Va = V0(H-AV)sin27Tf0(l 4- Af) - t -F A0 X IOOOa/24Vb = V0(l + AV) sin2irf0(l + Af) - t + A0 4- - it x 1000 V22LVc = V0(l + AV) sin27Tf0(l + AQ - t + A0 4- -it x 1000a/2其

28、中v0-一.单机无穷大电源节点的电压幅值。V一电压幅值扰动量AV = ge。sin(2irf - t + 0Xf频率扰动量Af =t sin 2nf , t + % ,A0一一相角扰动量A。= C3e° - sin(2M-t+03。为扰动的衰减因子，f为扰动频率，G C2、C3分别为AV、AL A8的扰 动初始幅值。8、02、%是扰动量的初始相角，单位为角。公式末尾乘以 1000四的目的是为了将标幺值换算成有名值。基于上述含扰动的无穷大电源压源模型，在Matlab平台搭建的子系统如图 fl-4所示。fl-4含扰动的无穷大电源子系统下面分别介绍各个模块的功能。输入量模块：左侧为参量输入

29、模块，生成含10个参量的u向量u=V0, Ci，C2, C3, 0 , 0 , 0 , f, o, t JLCJ即ul = Vo, u2 =u3 = C3u10 = t Fen模块Fen模块是函数表达式模块，可对输入量进行函数运算。Fcnl模块的函数为：ul*(l+u2*exp(u9)*sin(2*pi*u8*u10+u5*pi/180)*1000*sqrt(2)是对输入向量U的各个元素、u310进行函数计算的结果。即Fen 1模块输出的是也就是+ sin(2irft+e 同)XV0(14-AV) X 1000V2Fcn2模块的函数为：t + % 嬴2*pi*60*u10*(l+u3*exp(

30、u9)*sin(2*pi*u8*u10+u6*pi/180)2irf0 t 1 + C2ea - sin(2irf 也就是27if0(l 4- Af) - tFcn3模块中的函数为u4*exp(u9)*sin(2*pi*u8*u10+u7*pi/180)也就是AB oFcn2和Fcn3输出的结果在进行求和运算（在sum模块中实现）后，在Fcn4、Fcn5、Fcn6中加入A、B、C各相的初相角并计算相应的正弦值。Fcn4中进行的运算为sin（Fcn2中的参数+Fcn3中的参数）Fcn5中进行的运算为sin(Fcn2中的参数+Fcn3中的参数4EFcn6中进行的运算为sin(Fcn2中的参数+Fc

31、n3中的参数+|it) 也就是Fcn4 输出sin2irfo(l + Af) t + A0Fcn5 输出sin27rfo(l 4- Af) - t 4- A0 4-Fcn6 输出sin2Ttfo(l 4- Af) - t + A0 4- -it3Product!,、Product2、Products模块分别输出的就是最终反馈至Simulink主 系统的Va、Vb、%的值。Productl输出的结果二Fcn4输出的结果XFcnl输出的结果Product2输出的结果二Fcn5输出的结果X Fcnl输出的结果Products输出的结果二Fcn6输出的结果X Fcnl输出的结果它们的具体表达式为：%

32、= V0(l + AV) sin2irf0(l 4- Af) t + A04VB = Vo (1 + AV) sin27rf0(l 4- Af) - t + A0 4- -tt2% = V。(1 + AV) sin 2irf0 (1 + Af) t + A。+ -tt1.1.不同扰动量情况下机电振荡变化在该模型中，改变无穷大系统内扰动源的扰动频率，观察系统机电振荡响 应曲线：发电机功角b曲线、转速3曲线、输出电磁功率Pe曲线的变化。当CO.l 时，6、3、Pe稳态峰值如下图所示。表f-55、3、Pg峰值表扰动频率f00) ( X IO-*)Pe ( X 10-3 )6 max6 min5min

33、3 m检wminWm.WminRtrnwR minRmin0.1084.80680.34874.45730.8711-0.60551.4766一一0.2084.94179.73255.20852.055-1.16273.2177,0.4085.92578.6175730754.9109-3.28578.19666.9881-0.01186.99990.6086.7026773512935147.696-8.947616.643614.9-27342.20.7086.114578.26997.84467.0213-8.872815.894120.9-29.850.70.8085.828877.2

34、8718.54177.8114-11.70219.513426.5-59.886.30.9085.507574.88310.624510.428-16.29426.72248.1-99.214731.084.778871.75913.019814.807-21.59836.40574.6-140.42151.183.325665.174418.151224.854-3137256.226119.1-206332541.1293.02450.742842.281264.995-72.339137334279.8420.3700.11.149737345.260352.112782.209-90.

35、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.846245.51548331280.604-87.125167.729390.1-500.98911385.202251.100134.102161.893-65.183127.076

36、337.5-399.7737.21.575.985258.499917.485336.4837.31773.797231.2-269.5500.71.775.522865.038710.484124.544-24.83149375173.7-209.1382.81.975.961868.79587.16618.689-18.85637.545146-177.9323.92.176.313870.86485.44915.322-15.40730.729132.2-161293.22.376.498772.16274.33613.111-13.11826.229124.4-149.2273.62.

37、576.79472.95393.840111.696-12.01123.707120.7-143.9264.6相角扰动量的幅值(C3=0.02)时， 6 . A3、A Pe随扰动频率的变化如下图所示:7nnzuu180 -i An>1OU/ i on/ u100/ X wgn/ O V60 -40 一- 7V on -/u0 -IZV9Z86 IZP9Z 8rMm SZ6 ES o o o o o o oiVUUUUiiiiizzzi i i i i图件7A Pe变化曲线当0=0.05时,6、3、Pe稳态峰值如下图所示。表f5b、3、Pe峰值表扰动频率 f63 (X10-4)Pc ( X

38、10-3)66 min6 mNmintomin3 m秋-3minPg macP° minPg max-Pp min0.183.136480.94760.3630.26690.283.195380.73340.8488-0.48040.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.9383-16.41.083.209880.6453

39、3.2935-3.701313-18.71.183.138480.65223.431-4.013814.8-22.41.283.160380.27133.98544842218.1-33.81382.964479.5785.0664-6.079623.3-48.91.581.93177.25698.581-9.580948.6-86.71.780.77576.11629.9273-10.53264.4-100.81.980389176.54929.0465-9.374565.8-97.62.180.211977.0917.9844-8.157863.9-922.380.149377.53157

40、.0888-7.156462.6-88.42.580.085477.82826.3532-6.369461.1-84.11.1283.195880.61423.4898-4.179615.8-23.11.141.1683.088780.4823.6913-4.495116.2-291.1783.122380.44263.754-4.5791163-28.91.1883.073180.4013.8214-4.670816.8-30.1无扰动在不同PSS情况下，系统的稳定情况表:KL=30KI=40KH=160KL=30等幅振荡增幅振荡收敛慢KI=40增幅振荡等幅振荡收敛快KH=160收敛慢收敛快

41、等幅振荡PSS系统的Bode图如下:Lf bvdF 3nd图fl-11系统的Bode图KL二30、KI=40、KH=160时，系统的稳定性较好。只考虑KI、KH （也就是KL=O、KI=40. KH=160时），比同时考虑KL、KI、KH（也就是KL=30、KI=40> KH=160时）的稳定性要高。当减小模型中FL、Fl、FH参数时，曲线左移。反之，右移。参考文献1 PRABHA KUNDUR.电力系统稳定与控制.北京:中国电力出版社，2002刘取,电力系统稳定性及发电机励磁控制.北京:中国电力出版社，20073西安交通大学等.电力系统计算,北京:中国电力出版社，19784动态电力系统

42、的理论和分析.电力系统计算.北京:中国电力出版社，1978扰动量为6个周波时，6、3、已的响应曲线6响应曲线3响应曲线Pe响应曲线扰动量为6个周波时，仿真结果如下:序号63峰值时间/幅值半周期 时间半周期 差值峰值时间/幅值半周期 时间半周期 差值峰值时间 /幅值半周期 时间半周期 差值10.59/1341.13 /0.00590.551/0.048321.42 3.230.831.67/ 0.00510.540.894 /-0.1466034332.04 /11430.62-0.212366 /0.0040.6960.1561.404 /0.28680.510.16742.66 /80.73

43、0.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.17060.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.0574394 /

44、-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.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.01212734 /91.230.590.037.60A0.00070.563-0.036.735 /-0.02320.5870.002137.91 /97.110.57-

45、0.028.187 /0.00060.5870.0247332 /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.63792.520.5709.905 /-0.00040.57-0.0099.0405/- 0.00530.5755-0.00751710.20/95330.57010.483/0.000320.5780.0089.615 /0.

46、02940.5745-0.0011810.78/92.790.580.0111.055/0.00030.575-0.00310.202/0.00010.5870.0125191136/94.740.58011.628/0.00020.573-0.00210.748/0.02510.546-0.0412011.94/92.910.58012.191 /-0.00020.563-0.0111307 /0.00370.5590.013扰动量为4个周波时，仿真结果如下:序号63峰值时间/幅值半周期 时间半周期 差值峰值时间 /幅值半周期B悯半周期 片值峰值时间 /幅值半周期 时间半周期 差值10.40

47、84/110.70.699 /0.00560.9375/0.321120.9691/68.040.56071.191 /-0.00490.4921.4449 /-0.17310.507431.493/103.30.5239-0.03681.782 /0.00410.5910.0992.021/0.19910.57610.068742.037/73.010.5440.02012.274 /-0.00360.492-0.0992.5288 /-0.11690.5078-0.068352.5583 /98.70.5213-0.02272.832/0.00310.5580.0663.0716/0.14280.54280.03563.0971/76360.53880.01753329 /-0.00270.497-0.0613.5959 /-0.08470.5243-0.018573.6172/95350.5201-0.01873.878 /0.00230.5490.0524.1214 /0.10810.52550.001284.1467/78.740.52950.00954379 /-0.0020.501-0.04