
Citation: | LI Cuiran, YANG Qian, XIE Jianli, LYU Anqi. Throughput Performance Analysis and Optimization of Energy Harvesting Wireless Sensor Network Based on Simultaneous Wireless Information and Power Transfer[J]. Journal of Southwest Jiaotong University, 2024, 59(5): 1014-1022. doi: 10.3969/j.issn.0258-2724.20220625 |
In view of two-hop multi-relay transmission in energy harvesting wireless sensor networks (WSN), a wireless radio frequency power beacon (PB)-assisted energy harvesting relay model based on simultaneous wireless information and power transfer (SWIPT) was constructed. Under the condition that the relay node has the characteristics of capturing the source node, loop self interference, and PB signal energy, the outage probability and throughput of the destination node using two different receiving strategies, namely, selection combining (SC) and maximal ratio combining (MRC), were derived. Then, under multiple constraints such as ensuring communication quality of service (QoS), a relay selection algorithm was proposed to jointly optimize the time switching factor and power splitting factor with the goal of maximizing the throughput. Simulation and numerical results show that the PB transmit power, time switching factor, number of antennas, and power splitting factor significantly affect the system’s outage probability and throughput. When PB transmit power is 6 dBW, and the number of antennas is 3, compared with the random relay selection algorithm and the max-min relay selection algorithm, the system’s throughput gains under the SC strategy are 0.29 and 0.15 bit/(s·Hz), respectively, and those under the MRC strategy are 0.32 and 0.16 bit/(s·Hz), respectively.
高速磁浮具有高速快捷、运能强大、安全可靠、适应性好、绿色环保等优势,适用于中长途运输,能促进城际间同城化发展,也适用于大型城市群内部快速通勤. 相对高速轮轨铁路,高速磁浮列车摆脱了车轮与轨道黏滞系数的限制,能更经济地达到更高的速度[1-2].
长定子直线同步电机(LSM)具有高速度、高效率等优点,被应用于常导高速磁浮列车中. 同时,实现列车的牵引与悬浮. 悬浮系统通过动态调节励磁电流保持悬浮气隙的恒定[3-6],但是,励磁电流的变化会改变气隙磁场,从而引起推力波动[7-8]. 此外,高速磁浮列车在换步过程中会产生电流波动,也会引起推力波动[9-10]. 另一方面,直线电机的齿槽效应[11]和端部效应[12]等也会带来推力波动问题. 推力波动直接作用在磁浮列车上会产生振动,降低列车的乘坐舒适性. 磁悬浮车辆作为机械结构有机械固有频率,当推力波动频率与机械频率接近时,激发列车结构共振,影响列车安全运行[13]. 因此,对高速磁浮直线同步电机进行精确建模可为后续推力波动抑制方法设计与仿真打下基础.
传统的直线电机数学模型一般为集中参数模型[14]. 为考虑磁饱和与交叉饱和现象,可以将电感表示为电流的多项式函数,进一步提升模型精度[15]. 然而,集中参数模型基于线性化假设,无法准确体现推力波动现象,该建模方法具有一定局限性. 为反映电磁场中的电机谐波和非线性问题,研究人员提出一些利用磁场分析的先进半解析模型. 基于有限元分析结果和适当的基函数,场重构方法可以准确描述非饱和条件下磁场中的谐波特性[16]. 磁链重构的电机建模方法可以准确预测磁饱和情况下的转矩脉动[17]. 由于电机的磁链与电压、磁力和机械的动态方程均可通过磁共能进行描述[18],因此,磁共能的电机解析数学模型在简洁性与准确性上具有明显优势. 磁共能重构的方法被应用于开关磁阻电机[19]、直流励磁磁通切换电机[20]和旋转永磁同步电机[21],该方法可准确描述电机转矩脉动与磁饱和现象. 文献[22]将基于磁共能重构的分布参数建模方法应用于高速磁浮直线同步电机,验证分布参数建模方法可以精确表征磁场谐波及铁芯磁饱和对牵引控制的影响,但并未考虑高速磁浮实际运行中悬浮系统对牵引控制的影响.
集中参数推导的传统高速磁浮直线同步电机数学模型忽略悬浮系统调节、直线电机谐波和非线性特性的影响,导致模型精度不高,本文针对该问题,提出计及悬浮系统影响的高速磁浮直线同步电机分布参数建模方法. 通过仿真和硬件在环实验,验证本文所提建模方法的正确性.
常导高速磁浮列车车厢底部装有多个电磁铁模块,末节车厢每侧包含7.5个标准电磁铁模块,而中间车厢为8.0个. 电磁铁模块不仅提供悬浮力,同时与轨道上的定子铁芯和线圈构成电励磁长定子直线同步电机.
长定子直线同步电机的定子为分段结构,牵引变电站对不同的定子段分段供电. 因此,随着列车位置发生变化,直线电机的漏感、电阻以及线路阻抗参数等都会发生变化. 在建立长定子直线同步电机的数学模型时,需要将其分为3部分考虑:列车覆盖定子段部分、列车未覆盖部分定子段、轨旁馈电电缆部分.
首先,建立覆盖部分的d、q轴电压和磁链方程:
{udis=risid+pψdis−ψqisπv/τ,uqis=risiq+pψqis+ψdisπv/τ,ψdis=Ldisid+ψm,ψqis=Lqisiq, | (1) |
式中:udis和uqis分别为覆盖部分的d、q轴电压,id和iq分别为定子d、q轴电流,ψdis和ψqis分别为d、q轴磁链,ψm为动子励磁磁链,Ldis、Lqis分别为列车覆盖部分直线同步电机的d、q轴电感,p为微分算子,ris为覆盖部分的电阻,τ为悬浮电磁铁极距,v为列车运行速度.
磁极未覆盖部分的定子绕组电压和磁链方程为
{udos=Rosid+pψdos−ψqosπv/τ,uqos=Rosiq+pψqos+ψdosπv/τ,ψdos=(Los+Lon)id,ψqos=(Los+Lon)iq,Ros=ros+ron, | (2) |
式中:udos和uqos分别为列车磁极未覆盖部分d、q轴的定子绕组电压,ψdos和ψqos分别为列车磁极未覆盖部分d、q轴的定子绕组磁链,Los为磁极未覆盖部分定子绕组总漏感,Lon为馈电电缆的等效电感,Ros为磁极未覆盖部分总电阻,ros为磁极未覆盖部分绕组上的电阻,ron为馈电电缆等效电阻.
高速磁浮列车运动学方程为
{mdvdt=Fx−Fz,S=∫vdt, | (3) |
式中:m为列车质量,Fx为推力,Fz为阻力[23],S为行驶距离,t为时间.
本文利用有限元软件建立电磁铁模型,二维模型如图1,参数如表1所示.
参数名称 | 数值 | 参数名称 | 数值 | |
额定励磁电流/A | 20 | 动子极距/mm | 266.5 | |
定子绕组匝数/匝 | 1 | 铁芯厚度/mm | 185 | |
励磁绕组匝数/匝 | 270 | 定子槽距/mm | 86 | |
定子极距/mm | 258 | 额定气隙/mm | 10 |
电磁铁模块的磁共能计算式为
Wc(Is,If,β,x)=∫ρcoe(Is,If,β,x)dV, | (4) |
式中:Wc(·)为磁共能,Is为定子电流幅值,If为励磁电流幅值,β为定子电流相角,x为电磁铁位置,ρcoe为磁共能密度,V 为 电磁铁模块耦合段的体积.
考虑到磁浮列车采用多点悬浮控制,悬浮间隙可以稳定于额定气隙附近,但是励磁电流往往波动较大. 因此,本文重点考虑励磁电流波动对模型的影响.
不同工况下磁共能分布如图2所示. 由图2可知:磁共能在定子电流相角上的周期为2π,在动子位置角上的基波周期为π/3.
根据磁共能在定子电流相角和动子位置角2个维度上的周期性,对所有磁共能数值解进行二维傅里叶级数展开,如式(5).
Wc=N1∑k1=−N1N2∑k2=−N2Wck1k2(Is,If)ej6k1πxτ+jk2β, | (5) |
式中:$ {\text{W}}_{\text{c}{{k}}_{\text{1}}{{k}}_{\text{2}}} $(·)为x的第k1次,β的第k2次傅里叶系数,同时也是定子电流和励磁电流的多项式;N1和N2分别为磁共能在x和β维度上的最高阶次.
考虑到磁共能主要为基波和低次谐波,将N1和N2分别设定为磁共能在x、β维度上的最高阶次.
式(5)中的磁共能可以用矢量形式来表达,首先定义傅里叶基向量为
{V(x)=[e−j6N1πxτe−j6(N1−1)πxτ⋯ej6(N1−1)πxτej6N1πxτ]T,U(β)=[e−jN2βe−j(N2−1)β⋯ej(N2−1)βejN2β]. | (6) |
设计式(7)所示的二元多项式对式(5)中的傅里叶系数进行拟合.
Wck1k2(Is,If)=N3∑I=0N3−J∑J=0W(I)(J)ck1k2IIsIJf, | (7) |
式中:N3为定子电流和励磁电流的最高阶次,I和J分别为定子电流和励磁电流的阶次.
化简后的磁共能解析式为
Wc=U(β)Wck1k2(Is,If)V(x). | (8) |
对式(8)所示的磁共能模型求d、q轴电流的偏导数,即可得到定子d、q轴系下磁链的表达式为
{ψd(Is,If,β,x)=∂Wc(Is,If,β,x)∂id,ψq(Is,If,β,x)=∂Wc(Is,If,β,x)∂iq. | (9) |
d、q轴下的电压可表示为
[uduq]=[idiq]ris + [dψddtdψqdt]+πτ˙x[−ψqψd]. | (10) |
则电磁铁模块的推力为
Fx=32∂Wc(Is,If,β,x)∂x. | (11) |
高速磁浮直线同步电机建模方法如图3所示. 图中:k为当前时刻,k + 1为下一时刻;ifSM(k)为每个电磁铁模块的励磁电流;Fxn为第n个电磁铁模块的推力大小;Fxl和Fxr分别为左侧和右侧的电磁铁模块推力;FxL、FxR为高速磁浮列车左、右侧直线电机的推力,通过计算所有电磁铁模块的推力之和就可以得到左、右两侧直线电机的推力;NSM为高速磁浮列车电磁铁模块数量. 首先,根据列车编组配置,计算NSM;其次,计算逆变器输出的d、q轴电压,并根据式(2)计算出未覆盖部分的电压;随后,得到覆盖部分的电压,并均分到每个电磁铁模块;根据式(10),可以由电压方程计算得到d、q轴电流;接着,计算出每个电磁铁模块的推力;再次,根据列车编组数以及当前速度,计算列车的阻力;最后,左右两侧直线电机推力相加可以得到总的推力,并根据式(3)最终计算得到列车位置和速度.
通过MATLAB软件可以将式(8)计算得到的磁共能与有限元仿真获得的磁共能进行对比. 为验证磁共能重构模型在一个仿真周期内总的精度,本文采用均方根误差表示磁共能模型和有限元仿真结果之间的偏离程度. 图4是所有电流工况点下的磁共能误差分布情况,除了在边界点由于多项式拟合的边界效应导致磁共能误差略大外,其他工况下磁共能误差均低于0.1 J,说明本文拟合的磁共能模型精度较高.
在实验室中搭建基于RT-LAB的高速磁浮直线同步电机的硬件在环仿真平台,如图5所示. 图中:i为电流. 牵引系统在RT-LAB仿真机中,控制器芯片为TI的TMS320F28379D. 控制器和RT-LAB之间的通信通过一块接口板实现,接收RT-LAB中的速度和电流等模拟量,向RT-LAB发送控制指令.
本文提出的方法可以通过配置电磁铁模块的数量以模拟不同列车编组数的牵引工况. 图6显示了当列车编组设置为3、5、8时加速过程中的推力. 在3种编组情况下,电磁铁模块数量分别为23、39、63个. 3种情况下的平均推力分别为175、190、210 kN. 可以发现:高速磁浮列车在相同的加速度下,随着列车数量的增加会导致车身重量的增加,因此,为保证相同的加速度,推力也会随之增加;推力波动会随着推力幅值的增加而增加,因此,随着编组数量的增加,推力波动会进一步加剧.
图7是传统集中参数模型和所提出的建模方法下的推力对比结果. 可以发现:在3个阶段下,2种建模方法的平均推力是大致相等的,分别为190、38、−135 kN;通过计算得到的基于集中参数模型的推力波动分别为22.10、21.00、26.40 kN,所提方法的推力波动分别为37.30、28.45、28.20 kN,基于磁共能重构方法的推力波动在3个阶段均大于传统方法. 这是由于有限元计算得到的磁共能在一个基波周期内沿动子位置以6倍基波频率波动,并且关于电流矢量相角也是周期性的,波动的幅值受定子电流和励磁电流大小影响. 同时,传统集中参数表示的磁共能难以处理齿槽效应、端部效应和铁芯饱和造成的变化,而这些变化特征可以利用二维傅里叶向量和多项式系数矩阵来表示,如式(6)和式(7)所示. 通过建立磁共能解析模型,所提方法可反映出电机电磁场中的电机谐波和非线性问题.
为验证悬浮系统对牵引控制的影响,在所有电磁铁模块的额定励磁电流值上叠加一个谐波电流来模拟不同的励磁电流. 设置3组不同幅值和频率的励磁电流谐波:幅值为0.5 A,频率为5 Hz(谐波1);幅值为1.0 A,频率为10 Hz(谐波2);幅值为2.0 A,频率为20 Hz(谐波3). 不同阶段时的推力对比如图8所示. 由图可知,谐波3对应的推力波动最大,谐波1对应的推力波动最小. 表2列出3种励磁电流谐波时在加速、匀速和减速阶段时的推力波动和主要频率的占比. 可以发现:随着励磁电流波动幅值的增加,推力波动幅值也随之增加;推力的波动频率也与励磁电流的谐波频率相近.
阶段 | 类型 | 推力谐波占比/% | 推力波动/kN | ||
f =5 Hz | f =10 Hz | f =20 Hz | |||
加速 | 谐波 1 | 1.43 | 0.21 | 0.10 | 25.60 |
谐波 2 | 0.54 | 2.48 | 0.30 | 32.40 | |
谐波 3 | 0.29 | 0.16 | 5.14 | 39.30 | |
匀速 | 谐波 1 | 4.48 | 2.04 | 1.27 | 20.69 |
谐波 2 | 1.84 | 9.66 | 1.43 | 26.45 | |
谐波 3 | 2.10 | 0.63 | 21.75 | 31.46 | |
减速 | 谐波 1 | 3.86 | 0.10 | 1.01 | 22.00 |
谐波 2 | 0.59 | 7.98 | 0.37 | 33.40 | |
谐波 3 | 0.48 | 0.63 | 14.17 | 48.50 |
综合上述:当高速磁浮列车处于匀速运行时仅需要克服阻力,推力的幅值不大,但推力波动所占比例较大;当高速磁浮列车在加速和减速时,此时定子电流较大,推力波动幅值也比较大,但是,推力波动所占比例不大;不管是加速、匀速和减速工况,悬浮系统的调节都加剧了推力波动,并且在推力中引入了励磁电流同频的谐波,这将带来共振的危害. 针对悬浮系统导致的推力波动问题,一方面可以进一步优化悬浮系统的控制,降低悬浮系统对牵引控制的影响;另一方面,则需要改进高速磁浮的牵引控制,减小悬浮系统调节导致的推力波动问题.
针对集中参数模型难以表征高速磁浮悬浮系统的调节,电机谐波和非线性特性对牵引控制的影响,本文提出一种考虑悬浮系统影响的分布参数模型,通过对模型研究分析,得到以下结论:
1) 所提出方法可以通过配置高速磁浮电磁铁模块的数量反映不同列车编组时的牵引特性.
2) 与传统的集中参数模型相比,所提方法通过有限元获取磁共能数值解来建立直线同步电机数学模型,考虑电机的空间谐波以及铁芯的饱和效应,模型精度更高.
3) 高速磁浮直线同步电机推力波动幅值随着励磁电流谐波幅度的增大而增大,并且推力的波动频率与励磁电流的谐波频率相关.
本文提出的分布参数建模方法可以准确地表征悬浮系统对牵引控制的影响,是实现推力波动抑制的基础,有利于提高高速磁浮牵引控制性能.
[1] |
SUDEVALAYAM S, KULKARNI P. Energy harvesting sensor nodes: survey and implications[J]. IEEE Communications Surveys & Tutorials, 2011, 13(3): 443-461.
|
[2] |
SISINNI E, SAIFULLAH A, HAN S, et al. Industrial Internet of Things: challenges, opportunities, and directions[J]. IEEE Transactions on Industrial Informatics, 2018, 14(11): 4724-4734. doi: 10.1109/TII.2018.2852491
|
[3] |
王茜竹,胡洪瑞,徐勇军,等. 基于能量收集的UAV-D2D网络资源分配算法[J]. 电子与信息学报,2022,44(3): 976-986.
WANG Qianzhu, HU Hongrui, XU Yongjun, et al. Resource allocation for UAV-assisted D2D communications with energy harvesting[J]. Journal of Electronics & Information Technology, 2022, 44(3): 976-986.
|
[4] |
刘向丽,刘冬妮,李海娇,等. 能量协作下中继系统性能推导及仿真分析[J]. 西安电子科技大学学报,2019,46(5): 91-97.
LIU Xiangli, LIU Dongni, LI Haijiao, et al. Derivation of the performance of the relay system and simulation analysis of the system under energy cooperation[J]. Journal of Xidian University, 2019, 46(5): 91-97.
|
[5] |
HUANG K B, LAU V K N. Enabling wireless power transfer in cellular networks: architecture, modeling and deployment[J]. IEEE Transactions on Wireless Communications, 2014, 13(2): 902-912. doi: 10.1109/TWC.2013.122313.130727
|
[6] |
JIA X H, ZHANG C Z, KANG J M, et al. Joint beamforming design and time allocation for wireless powered asymmetric two-way multirelay network[J]. IEEE Transactions on Vehicular Technology, 2018, 67(10): 9641-9655. doi: 10.1109/TVT.2018.2860001
|
[7] |
LIU H W, KIM K J, KWAK K S, et al. Power splitting-based SWIPT with decode-and-forward full-duplex relaying[J]. IEEE Transactions on Wireless Communications, 2016, 15(11): 7561-7577. doi: 10.1109/TWC.2016.2604801
|
[8] |
ZHANG R, HO C K. MIMO broadcasting for simultaneous wireless information and power transfer[J]. IEEE Transactions on Wireless Communications, 2013, 12(5): 1989-2001. doi: 10.1109/TWC.2013.031813.120224
|
[9] |
施安妮,李陶深,王哲,等. 无线携能通信系统中的吞吐量最优化全双工中继选择策略[J]. 小型微型计算机系统,2021,42(9): 1906-1912.
SHI Anni, LI Taoshen, WANG Zhe, et al. Optimization scheme for the SWIPT-NOMA opportunity cooperative system[J]. Journal of Chinese Computer Systems, 2021, 42(9): 1906-1912.
|
[10] |
TRAN TIN P, VAN-DUC P, NGUYEN T N, et al. Performance analysis in the decode-and-forward full-duplex relaying network with SWIPT[J]. The Scientific World Journal, 2021, 2021: 1-8.
|
[11] |
TIN P T, NGUYEN T N, TRAN D H, et al. Performance enhancement for full-duplex relaying with time-switching-based SWIPT in wireless sensors networks[J]. Sensors, 2021, 21(11): 3847.1-3847.16.
|
[12] |
WANG D X, ZHANG R Q, CHENG X, et al. Relay selection in full-duplex energy-harvesting two-way relay networks[J]. IEEE Transactions on Green Communications and Networking, 2017, 1(2): 182-191. doi: 10.1109/TGCN.2017.2686325
|
[13] |
JING Y D, JAFARKHANI H. Single and multiple relay selection schemes and their achievable diversity orders[J]. IEEE Transactions on Wireless Communications, 2009, 8(3): 1414-1423. doi: 10.1109/TWC.2008.080109
|
[14] |
李陶深,施安妮,王哲,等. 基于SWIPT的吞吐量最优化NOMA全双工中继选择策略[J]. 通信学报,2021,42(5): 87-97.
LI Taoshen, SHI Anni, WANG Zhe, et al. Optimal relay selection for full duplex SWIPT-NOMA systems with maximal throughput[J]. Journal on Communications, 2021, 42(5): 87-97.
|
[15] |
WANG Z H, CHEN Z Y, XIA B, et al. Cognitive relay networks with energy harvesting and information transfer: design, analysis, and optimization[J]. IEEE Transactions on Wireless Communications, 2016, 15(4): 2562-2576. doi: 10.1109/TWC.2015.2504581
|
[16] |
罗轶,孔静恬,董健,等. 能量收集认知多跳中继网络中断性能分析及优化[J]. 电子与信息学报,2021,43(10): 2920-2927.
LUO Yi, KONG Jingtian, DONG Jian, et al. Outage performance analysis and optimization of energy harvesting cognitive multihop relay networks[J]. Journal of Electronics & Information Technology, 2021, 43(10): 2920-2927.
|
[17] |
HU Z W, YUAN C W, ZHU F C, et al. Weighted sum transmit power minimization for full-duplex system with SWIPT and self-energy recycling[J]. IEEE Access, 2016, 4: 4874-4881. doi: 10.1109/ACCESS.2016.2593914
|
[18] |
NGUYEN T N, MINH T H Q, TRAN P T, et al. Adaptive energy harvesting relaying protocol for two-way half-duplex system network over rician fading channels[J]. Wireless Communications and Mobile Computing, 2018, 2018: 1-10.
|
[19] |
GRADSHTEYN I S, RYZHIK I M, JEFFREY A, et al. Table of integrals, series, and products[M]. 7th ed. San Diego: Academic Press, 2007: 366-369.
|
[20] |
CHONG E K P, ŻAK S H. An introduction to optimization[M]. Hoboken: Wiley, 2008.
|
[21] |
WANG D X, ZHANG R Q, CHENG X, et al. Full-duplex energy-harvesting relay networks: capacity-maximizing relay selection[J]. Journal of Communications and Information Networks, 2018, 3(3): 79-85. doi: 10.1007/s41650-018-0027-0
|
[22] |
高云波,程璇,李翠然,等. T2T和T2G混合网络中的功率分配算法[J]. 西南交通大学学报,2023,58(5): 1126-1134,1179.
GAO Yunbo, CHENG Xuan, LI Cuiran, et al. Power allocation algorithm in T2T and T2G hybrid network[J]. Journal of Southwest Jiaotong University, 2023, 58(5): 1126-1134,1179.
|
[23] |
TRAN TIN P, THE HUNG D, NGUYEN T, et al. Secrecy performance enhancement for underlay cognitive radio networks employing cooperative multi-hop transmission with and without presence of hardware impairments[J]. Entropy, 2019, 21(2): 217.1-217.16.
|
参数名称 | 数值 | 参数名称 | 数值 | |
额定励磁电流/A | 20 | 动子极距/mm | 266.5 | |
定子绕组匝数/匝 | 1 | 铁芯厚度/mm | 185 | |
励磁绕组匝数/匝 | 270 | 定子槽距/mm | 86 | |
定子极距/mm | 258 | 额定气隙/mm | 10 |
阶段 | 类型 | 推力谐波占比/% | 推力波动/kN | ||
f =5 Hz | f =10 Hz | f =20 Hz | |||
加速 | 谐波 1 | 1.43 | 0.21 | 0.10 | 25.60 |
谐波 2 | 0.54 | 2.48 | 0.30 | 32.40 | |
谐波 3 | 0.29 | 0.16 | 5.14 | 39.30 | |
匀速 | 谐波 1 | 4.48 | 2.04 | 1.27 | 20.69 |
谐波 2 | 1.84 | 9.66 | 1.43 | 26.45 | |
谐波 3 | 2.10 | 0.63 | 21.75 | 31.46 | |
减速 | 谐波 1 | 3.86 | 0.10 | 1.01 | 22.00 |
谐波 2 | 0.59 | 7.98 | 0.37 | 33.40 | |
谐波 3 | 0.48 | 0.63 | 14.17 | 48.50 |
参数名称 | 数值 | 参数名称 | 数值 | |
额定励磁电流/A | 20 | 动子极距/mm | 266.5 | |
定子绕组匝数/匝 | 1 | 铁芯厚度/mm | 185 | |
励磁绕组匝数/匝 | 270 | 定子槽距/mm | 86 | |
定子极距/mm | 258 | 额定气隙/mm | 10 |
阶段 | 类型 | 推力谐波占比/% | 推力波动/kN | ||
f =5 Hz | f =10 Hz | f =20 Hz | |||
加速 | 谐波 1 | 1.43 | 0.21 | 0.10 | 25.60 |
谐波 2 | 0.54 | 2.48 | 0.30 | 32.40 | |
谐波 3 | 0.29 | 0.16 | 5.14 | 39.30 | |
匀速 | 谐波 1 | 4.48 | 2.04 | 1.27 | 20.69 |
谐波 2 | 1.84 | 9.66 | 1.43 | 26.45 | |
谐波 3 | 2.10 | 0.63 | 21.75 | 31.46 | |
减速 | 谐波 1 | 3.86 | 0.10 | 1.01 | 22.00 |
谐波 2 | 0.59 | 7.98 | 0.37 | 33.40 | |
谐波 3 | 0.48 | 0.63 | 14.17 | 48.50 |