2024年2月28日发(作者:2017园区初三数学试卷)

ELECTRICDRIVE2021Vol.51No.6电气传动2021年第51卷第6期双三相永磁同步电机模型预测电流控制研究宋文祥,任航(上海大学机电工程与自动化学院,上海200444)摘要:六相逆变器为双三相永磁同步电机提供了丰富的电压矢量资源,能够使预测电流控制变得更加精准,但更多的电压矢量会带来算法计算量过大的问题,同时双三相电机的谐波电流会使电机的损耗增加,需要对其进行抑制。提出了一种改进的模型预测电流控制方法。利用最外围大矢量与次外围中矢量在z1z2子平面的。根据定子磁链所在扇区的位置确定出更小范围内的8个预测电压矢量,从而减少了系统运算量。同时以d,q轴电流误差项作为价值函数,消除z1z2子平面的电流误差项,如此可避免权重系数的整定。通过实验对所研究方法进行了验证,结果表明所提MPCC方法可以有效地降低谐波电流,并且具有良好的控制性能。关键词:双三相永磁同步电机;模型预测电流控制;谐波电流抑制;预测电压矢量中图分类号:TM28文献标识码:ADOI:10.19457/22559方向相反的特性,在一个控制周期内将两个矢量按照相应比例结合并作用于电机,可实现抑制谐波电流的目ResearchonModelPredictiveCurrentControlofDualThreePhasePermanentMagnetSynchronousMotorSONGWenxiang,RENHang(SchoolofMechatronicEngineeringandAutomation,ShanghaiUniversity,Shanghai200444,China)Abstract:Sixphaseinverterprovidesabundantvoltagevectorresourcesfordualthree-phasepermanentmagnetsynchronousmotor,r,ametime,theharmoniccurrentofdualthree-phasemotorwillincreasethemotorloss,ore,animprovedmodelpredictivecurrentcontrol(MPCC)ingtothecharacteristicthatthelargestvectoroftheoutermostregionisoppositetothemiddlevectorofthesubperipheryinthez1z2sub-plane,thetwovectorswerecombinedaccordingtoacertainpropormore,thestatorfluxlinkagepositionwasobserv,qaxiscurrenterrortermwasusedasthevaluefunctiontoeliminatethecurrenterrortermofz1z2sub-plane,ultsshowthatMPCCcands:dualthreephasepermanentmagnetsynchronousmotor(DTP-PMSM);modelpredictivecurrentcontrol(MPCC);harmoniccurrentsuppression;predictivevoltagevector随着电力电子技术、微控制器技术和电机控制理论的发展,以及工业应用场合的需求,多相电机及驱动系统以其低压大功率输出、高可靠性、低转矩脉动的特点吸引了越来越多的学者研究[1-3]。其中,双三相永磁同步电机驱动系统是当前研究热点之一。同时,模型预测电流控制技术具有结构简单、响应速度快、灵活度高的特点而受到人们广泛关注[4-7]。作者简介:宋文祥(1973—),男,博士,教授,Email:**************.cncontrol,MPCC)是一种通过衡量不同的电压矢量对电机状态产生的影响,筛选出下一时刻最优的MPCC取代了传统的比例积分(PI)控制器和脉宽题,并且具有更快的转矩响应能力[8]。MPCC的研模型预测电流控制(modelpredictivecurrent工作电压矢量的控制策略。相比于矢量控制,调制器,避免了控制器参数整定和复杂计算的问究广泛应用于三相电机。对于三相电机而言,3

电气传动2021年第51卷第6期MPCC电机离散数学模型,只需针对静止坐标系或者旋转坐标系下的通过预测模型计算下一时刻的电流值,利用价值函数在线寻优得出最优电压矢量[9-11]。而对于六相电机,根据矢量空间解耦理论[12],解耦变换后有三个不同的子平面,需要同时考虑两个子平面,即基波子平面和谐波子平面(零序子平面为0)。然而,只有基波子平面负责产生转矩,谐波子平面不参与机电能量转换,但该平面很小的电阻和漏感能够造成很大的谐波电流[13]。同时,六相电压源型逆变器产生的电压矢量在不同的子平面有不同的幅值和方向,这就大大地增加了设计预测控制系统的难度。目前,对多相电机模型预测控制的研究并不多见。文献[14-15]对六相感应电机MPCC进行了研究,考虑了电压矢量个数和不同的价值函数对模型预测的影响,并讨论了实现该控制算法所需要的时间,证明MPCC策略在多相电机应用方面的有效性。文献[16]针对五相永磁同步电机提出一种虚拟电压矢量MPCC控制,充分利用电压矢量的特点,将虚拟矢量作为预测矢量,从而达到降低谐波电流的目的。文献[14-16]虽都实现了多相电机模型预测控制,且降低了谐波电流,但价值函数过于复杂,需要对权重系数整定,整定过程根据仿真实验结果反复调整,无疑会增加实验的复杂性。文献[17]运用多目标方式取代权重系数,通过建立两个价值函数,利用一个平均评价准则来选取使磁链、转矩误差都较小的电压矢量。文献[18]建立了一种自适应机制,运用两个价值函数,使转矩和磁链误差限制在初始界限内,从而消除了权重系数。文献[17-18]都是以三相电机为研究对象,六相电机比三相电机更为复杂,解耦后多一个谐波子平面,因此关于六相电机权重系数的问题需进一步研究。文献[19]针对双三相永磁同步电机提出一种模型预测转矩控制策略,通过分析电压矢量在基波子平面和谐波子平面的方向和幅值,建立抑制谐波电流的开关表,新的开关表有效地抑制了谐波电流并降低了系统的计算量。文献[20]针对T型三电平双三相永磁同步电机驱动系统的共模电压,利用系统空间零共模电压矢量,显著抑制了系统共模电压,并很好地控制了电机电流。phase为了抑制双三相永磁同步电机(dualthree统计算量。提出了一种改进的PMSM,DTP-PMSM)的谐波电流,DTP-PMSM并降低系模型4宋文祥,等:双三相永磁同步电机模型预测电流控制研究预测电流控制方法。通过在一个控制周期内选择作用两个电压矢量,即一个外围大矢量和一个次外围矢量,调整两个矢量的作用时间使z1z2子平面的平均电压幅值为零,从而实现谐波电流的抑制。对定子磁链进行观测,确定磁链所在扇区,将预测电压矢量个数由12个减少为8个,降低了系统的计算量。以αβ子平面的电流跟踪误差作为价值函数,消除了权重系数,选择出最优电压矢量。最后通过实验研究验证了该方案的有效性。1双三相永磁同步电机数学模型图1为双三相永磁同步电机结构示意图。图1双三相永磁同步电机结构Fig.1Structuremagnetsynchronousofdoublethreemotorphasepermanent不同于六相对称永磁同步电机,DTP-PMSM的两套三相绕组采用互差30°的方式放置,并且两套绕组的中性点相互隔离,所以又称六相不对space称电机。根分量映decompation据射到相互,正VSD矢量交)空间解耦变换(vector的,αβ将双三相电机的所有,z面[12],其中,αβ子平面为基波子平面,1z2及o1o2只有该平3个子平面分量参与机电能量转换,z平面,该平面分量不对气隙磁通和转矩做贡献,1z2子平面为谐波子o1o2子平面为零序子平面,由于中性点相互隔VSD离,零序谐波分量不能流动,变换矩阵为该平面分量为零。éêê1ê1-12-223êê1-0ùúú1êêêê0-1úúúúúê123Τs=3êêêê1-êê--23121232ê2-0úúúúúúúúêêêê01-23ê231232321-1úúúúëúúú(1û1)

宋文祥,等:双三相永磁同步电机模型预测电流控制研究然后利用旋转变换式,将静止坐标系转换到同步旋转坐标系。由于只有αβ子平面参与机电能量转换,所以无需考虑z旋转变换式如下:1z2子平面、o1o2子平面。éêêë-sincosθθsincosθθ0ùTs/r=式中:θ为转子位置角度:0I为04维单位矩阵。I0úú(2)4û同步旋转坐标系下αβ4子平面的数学方程为éêudùéRëuú=ês0ùúéêidùú+péêΨdùú+ωeéê-Ψqùú(3)qûëéêΨ0RsûëiqûdùéLequdëΨú0ëΨqûëΨdûùúéêidùúqû=êë+éΨfù其中0Lequêqûëiqûë0úû(4)Lequdequi=3Ld+LlLq式中:u=3Lq+Lld,uq,d,iq,Ψd,Ψq为αβ子平面的定子电压、电流、磁链;p为微分算子;R角频率;Ψf为永磁体磁链;Ls为定子电阻;ωequd,Leque为电q分别为d,q轴电感;Ld,Lq分别为d,q轴主自感;L自感。l为定子绕组漏z1z2子平面的数学方程为éêuz1ùëuú=éêRs0ùúéêiz1ùú+péêΨz1ùúé0RsûëiëΨ(5)z2ûëz2ûz2ûêΨz1ùéLëΨú=êl0ùúéêiz1ùú(6)z2ûë式中:ui0Llûëiz2ûz1,uz2,z1,iz2,Ψz1,Ψz2分别为z1z电流、磁链。2子平面的定子电压、电磁转矩方程为Te=3pn(iqΨd-idΨq(7)式中:T模型预测电流控制e为电磁转矩;p)n为极对数。2.12预测模型首先,基于αβ子平面的双三相永磁同步电机数学模型,将磁链方程(4)代入电压方程(3),通过离散化处理得预测电流模型:ìïTïid(k+1)=idï(k)+Lequsd[ud(k)-ïíRsid(k)+ωLequqiqïïïiq(k+1)=iqï(k)+T(k)]Lequsq[uq(k)-îRsiq(k)-ωLequdid(k)-ωΨf]式中:k)(k)(8)id(,iq,ud(k),uq电流、电压值;+1),i(k)为k时刻的d,q轴id(kq(k+1)为k+1时刻的d,电气传动2021年第51卷第6期q轴电流值。同理可得z式所示:1z2子平面的预测电流模型,如下ìïïiz1í(k+1)=iTz1(k)+Lsl[uz1(k)-Rsiz1(k)]ïïîi+1)=iTz2(kz2(k)+Ls[u(9)z式中:iul2(k)-Rsiz2(k)]z1(k),iz2(k),z1z(k),uz11z时刻的2子平面电流、电压值;k2(+k)1分别为)(k+k时刻的iz1(,iz21)为k+选取电压矢量z1z2子平面电流值。2.2图2为六相逆变器电压矢量图,从图中可以看出,六相电压源逆变器共有49个有效矢量,可根据αβ子平面矢量的幅值,将其分为4组矢量和一个零矢量,每组矢量12个,分别为最外围大矢量|umax|=0.644Udc、次外围中矢量|umidl|=0.471Udc、次0.173内围中小矢量U|umids|=0.333Udc、最内围小矢量|umin|=dc。图2六相逆变器电压矢量Fig.2Voltagevectorsofsix-phaseinverter5

电气传动2021年第51卷第6期电流值,MPCC选取最接近给定电流值的电压矢量作为通过计算电压矢量作用于电机产生的最优电压矢量。将49个矢量全部代入预测方程进行计算,这种方式存在计算量大、效率低、谐波电流大等缺点。为了减小计算量、降低谐波电流,一般采用最外围12个大矢量作为预测电压矢量。然而最外围大矢量在谐波子平面具有一定的电压,会引起较大的谐波电流。本文将最外围大矢量和次外围中矢量相结合作为预测矢量,以实现最优电压矢量的选取和抑制谐波电流的目的。如图2所示,最外围大矢量与次外围中矢量在αβ子平面方向相同,而在z如uu1z2子平面则方向相反,周期内作用两个矢量,64,46。根据这个特性,可选择在一个控制即一个大矢量和一个中矢量,通过调整它们的作用时间,使谐波子平面的平均电压幅值为零,如下式所示:ìíTs|uαβ|=μ|umax|+(TsîT-μ)|umidl|s|u(10)z1z2|=μ|umin|-(Ts由两个子平面建立两个等式,第一个等式为-μ)|umidl|αβ子平面的等式。依据伏秒积平衡原理,在一个控制周期Ts二个等式为μ内,),等效的电压矢量幅值为大矢量作用时间为μ,中矢量作用时间为(Ts-z|uαβ|。同理,第为|u1z2子平面的等式,等效的电压幅值z1z2令|。|uz1z2|=0,解式{(μ10=)0.731可得:Ts图3为预测电压矢量,|uαβ|=0.596一个大矢量和一个中U(11)dc矢量组成一组作为预测电压矢量,共12组。图3预测电压矢量Fig.3Predictedvoltagevector由式(11)可知,当选择一组电压矢量作用时,如U2(u64,u46),让u64作用0.73Ts,u46作用6宋文祥,等:双三相永磁同步电机模型预测电流控制研究0.270.596TsU,即可得到一个等效在αβ子平面幅值为dc、在z1z2子平面幅值为0的电压矢量。2.3减少预测电压矢量个数由前文分析,MPCC选取12个预测电压矢量代入预测模型进行迭代计算,相比49个电压矢量,预测矢量已减少许多,但对于系统总的运算量依然较大,需要较久的计算时间。为降低运算量,本文通过观测定子磁链位置,将预测电压矢量由12个减少到8个,实现减少计算时间的目的。图4为预测电压矢量优化选择图。图4预测电压矢量优化选择根据定子磁链所在扇区选择预测电压矢量,Fig.4Optimalselectionofpredictivedvoltagevector若定子磁链位于扇区Ⅰ,为避免电压矢量选取变化过大,降低电流波动,U6,U7与定子磁链方向相反,不作为预测电压矢量。同时希望电流具有快速的动态响应,预测矢量不考虑最接近扇区Ⅰ的两个矢量U1,U12,因此,仅需选择8个电压矢量,即U2,U3,U4,U5,U8,U9,U10,U11作为预测电压矢量,从而减少预测电压矢量个数。其它扇区的预测电压矢量可依次推理得出,在此不再赘述。定子磁链利用VSD变换的定子电流iα,iβ以及转子的位置进行观测,将其带入下两式即可得到定子磁链所在位置é(Lequdequé(LequθΨequ。(LequequêΨêëΨúê+2Lq)+d-2Lq)cos(2θ)d-2Lq)sin(2θ)ùúαùβû=êêë(Lequd-Lequq)(Lequd+Lequq)(Lequd-Lequúú·é2sin(2θ)2-2q)cos(2θ)úûêiαùéëiúû+êëcossinθθùúûΨβfθ(12)Ψ=arctan(Ψβ/Ψα)(13)2.4价值函数对于双三相永磁同步电机,MPCC一般以d,q

宋文祥,等:双三相永磁同步电机模型预测电流控制研究轴电流跟踪差值与谐波电流值之和作为价值函数,如下式所示:g1=|i*dλ[|i-*id(k+1)|+|i∗q-iq(k+1)|+(14z1式中:i*d,i∗-iz1(k+1)|+|i∗)z2-iz2(k+1)|]q为d,q轴电流给定值;i*z1,i∗z2为zz;λ为权重系数,一般按照同等权重的1,流给定值02轴电原则设计,在实际应用中需优化调整,以获得良好的控制效果。由于文章采用两个电压矢量抑制谐波电流,对于价值函数而言,无需添加谐波电流作为评价指标,只需通过d,q轴电流误差跟踪来选取最优的一组电压矢量作为价值函数,双三相永磁同步电机模型预测电流控制系i如下式所示:g2=|i*d-d(k+1)|+|i∗q-iq(k+1)|(15)统框图如图5所示。图5双三相永磁电机模型预测电流控制系统框图Fig.5BlockcontroldiagramsystemofformodelDTP-PMSMpredictivecurrent控制方法如下:旋转1变)将换k,时刻采样的电流值通过得到ikikikikVSDik变换和α,β,d,q,代入ikα,β,θ到式12)得到定子磁链所在扇区,进而确定8个预测电压矢量。i预测模型式2)把预测电压矢量和(8)中,对k+1时刻电流d,q轴电流iikd,kq代入到kd+1,ikq+1进行计算预测。d,q轴电流跟踪误差最小的电压矢量,3)利用价值函数式(15)在线寻优,将其作为筛选出使最优电压矢量。34实验研究)将最优的电压矢量作用于电机。实验研究所用的双三相永磁同步电机控制系统平台配置,如图6所示。其中,控制系统采用的微控制器为电气传动2021年第51卷第6期TMS320F2812kW,双三相永磁同步电机功率为5.5下:,额定功率对拖一台直流发电机作为负载,5.5kW,额定电压38012.0V电机参数如,mH额定电流感4.642.0A,定子电阻mH,永磁体磁链1.52Ω,0.585d轴电感Wb,极对数为,q3轴电,额定转速1500r/min,系统采样频率5kHz,直流母线电压Udc为540V。图6实验系统配置Fig.6Experimentalsystemconfiguration实验研究中首先对只使用最外围大矢量的模型预测电流控制算法进行了验证,MPCC此处简称为1。然后进一步地采用最外围大矢量和次外围中矢量相结合的模型预测电流控制算法对电机MPCC进行控制,简称MPCC2。为了使MPCC1与MPCC2在同一个采样控制的周期下进行比较,1也使用8个预测电压矢量进行控制。图7为转速100r/min下MPCC7a为整体电流波形,负载为1控制的电流波形。图20N·m,加减载前后相电流波形未有明显变化。分别观察空载与加载时电流运行放大波形,如图7b所示,电机空载运行时,A,U相电流并不为零,并且有效值较大为9A,同时谐波子平面的电流iz1,iz2有效值也在9A左右,其原因在于双三相电机的谐波子平面阻抗很小,即便较小的电压也将会引起很大的谐波电流,iz1,iz2幅值高,导致通过双三相电机的相电流过大。从图7c带载运行结果看出,相电流幅值较空载时未有明显增加,且电流波形严重畸变,原因在于谐波电流的幅值很大,增加的基波电流不足以使相电流产生明显变化。图8给出了转速300r/min下MPCC运行结果。观察发现,在突加60%额定负载,1加减载负载转矩为20N·m,转速略有波动,能够保持给定转速运行,但存在较大的转矩波动,同时相电流在加减载前后未有明显变化,电流有效值约10A。从图8b稳态结果看到,相电流波形严重畸变,主要原因如前所述,谐波电流含量过大,使得相电流不能反映基波电流的变化,转矩波动10N·m左右,电机运行效果较差。7(

电气传动2021年第51卷第6期图7转速100r/min的相电流与谐波电流波形(MPCCFig.7Waveformscurrentat100ofphase1)r/mincurrent(MPCCandharmonic1)图8转速300r/min的负载转矩突增与突减运行结果(MPCCOperation1)Fig.8andreductionresultsat300ofsuddenr/min(loadMPCCincrease1)8宋文祥,等:双三相永磁同步电机模型预测电流控制研究图9为MPCC2运行的加减载工作波形。图9转速100r/min的相电流与谐波电流波形(MPCCFig.92)Waveformcurrentatof100phaser/mincurrent(MPCCandharmonic2)图9a为整体电流波形,负载为20N·m,可以看出,9b电流的幅值在加减载前后有明显变化。图电流有效值约为电机空载情况下的相电流与谐波电流,1A,与未加谐波抑制算法的谐波谐波电流相比,60电流波形正常,%额定负载运行结果,有显著的抑制效果。图9c为电机带体现了MPCC相电流有效值达到4.0A,控制效果。2在基波子平面的良好图10为不同转速情况下双三相电机MPCC的稳态运行实验结果,负载转矩20N·m。由图2102可以看出,不同转速情况下,转矩波动在上下带载运行。N·m,同时转速保持平稳,在不同转速情况下稳定图11给出了转速900r/min时,电机突加与突减负载动态性能实验结果。由图11可以看出,电机转矩在负载突变时能够快速响应,对动态中的转速变化响应迅速,电机运行状态良好。

宋文祥,等:双三相永磁同步电机模型预测电流控制研究图10双三相电机稳态运行结果(MPPCFig.10Steadystateoperation2)(MPPCresultsofdualthreephasemotor2)4结论本文研究了DTP-PMSM驱动系统的模型预测电流控制策略,采用两个电压矢量相结合的方法,根据定子磁链位置选取预测电压矢量,无需将每个电压矢量代入预测模型计算,最终通过实验研究对模型预测电流控制算法进行了验证,得到如下结论:电压大矢量作用于电机,1)传统的模型预测电流控制仅采用最外围虽然谐波子平面的平均电压值小,但该平面具有很小的阻抗,能够引起较大的谐波电流,使得相电流幅值过大,不能正常运行。波子平面方向相反的特性,2)利用最外围大矢量与次外围中矢量在谐每个采样周期选择两个电压矢量作用于电机,使得谐波子平面的平均电压幅值为零,可以有效抑制双三相电机的谐波电流。3)根据定子磁链所在扇区选择预测电压矢电气传动2021年第51卷第6期图11电机负载突变的动态性能波形Fig.11withDynamicsuddenperformanceloadchangewaveformsofmotor量,可将12个预测电压矢量减少至8个,预测矢量减少后的控制算法运行稳定,具有良好的稳态和动态性能。实验结果验证了所研究控制方案的有效性和可行性。参考文献[1]王永兴,制[J].中国电机工程学报,温旭辉,赵峰.多相永磁同步电机多维优化矢量控2015,35(10):2534-2543.[2]电工技术学报,刘自程,李永东,2017郑泽东,32(.多相电机控制驱动技术研究综述24):17-29.[J].[3]陶涛,述[J].赵文祥,中国电机工程学报,程明,等.多相电机容错控制及其关键技术综2019,39(2):316-326,629.[4]ZhangpleYongchang,BaiYuning,YangHaitao.(drives[J].-vector-8):sedTransactionsmodelpredictiveAuniversalmulti⁃onPowercontrolElectronicsofinduction,2018motor,33[5]齐昕,电流控制吴文昊,[J].中国电机工程学报,吴琳,等.基于时间辅助信息的感应电机预测2019,39(16):4927-4934,[6]4995.王治国,郑泽东,李永东,等.基于增量观测器的MPCC一拍9

电气传动2021年第51卷第6期滞后补偿方法研究[J].电气传动,2020,50(6):119-123.[7]徐艳平,型预测电流控制王极兵,王建渊,[J].电气传动,等.考虑预测误差的改进双矢量模2018,48(9):62-66.[8]ZhangtipleIEEE-vectorYongchangTransactions-based,modelXuDonglinonIndustrialpredictive,eralized,PMSM2018,65drives[J].mul⁃9356-9366.(12):[9]张永昌,机模型预测控制杨海涛,[J].魏香龙电工技术学报,.基于快速矢量选择的永磁同步电2016,31(6):66-73.[10]林宏民,机模型预测直接转矩控制吴晓新,乐胜康,等.[J].基于三电平优化矢量的异步电电机与控制学报,2018,22(8):65-74.[11]WangdictiveZhiqiang,YuAnbo,LiXinmin,lcontrolbasedonfuzzyalgorithmAfornovelPMSM[J].currentIEEEpre⁃2019,7(of2)Emerging:ectedTopicsinPowerElectronics,[12]ZhaothreeYifan,ectorPWMcontrolofdual(tion[J].-phase5):ductionTransactionsmachineonIndustryusingvectorApplicationsspacedecomposi⁃,1995,31[13]周长攀,制谐波电流抑制研究苏健勇,杨贵杰[J]..双三相永磁同步电机直接转矩控电机与控制学报,2015,19(9):[14]oofpredictiveF,ArahalcurrentM,GregorR,ofconceptstudythree-phaseACmachines[J].controlforIEEEVSITransactions-drivenasymmetricalonIndustrialdual10宋文祥,等:双三相永磁同步电机模型预测电流控制研究[15]ElectronicsBarrerodictiveF,,Arahal2009,M56(,Gregor6):1937-1954.R,-stepmodulationpre⁃phaseinductioncurrentcontrolmachine[J].methodIEEEfortheTransactionsasymmetricalondualIndustrialthree-[16]ElectronicsXuemodelCheng,,2009Song,56Wensheng(6):dictivecurrentcontrol,-phaseFinitepermanentcontrol-set[17]IETsynchronousRojasElectricandfluxCcontrolA,RodriguezPowermachineApplicationsbasedwithoutJweighting,Villarroel,2017onvirtualvoltagevectors[J].-mag⁃factors[J].F,,11et(al.5):onIndustrialElectronics,2013,60(2):ansactions[18]WangdictiveFengxiangtorquecontrol,XieforHaotianinduction,ChenmachinesQing,elpre⁃(factors[J].2):ansactionsonPowerElectronics,weighting2020,35[19]LuoforaYixiaoIEEEdual,TransactionsthreeLiuChunhua.-phaseAsimplifiedmodelpredictivecontrolonPMSMIndustrialwithElectronicsreducedharmonic,2018,65currents9079-9089.[J].(11):[20]徐质闲,驱动零共王政,模电王学庆,压模型等预.测T型三电平双三相永磁同步电机控制[J].中国电机工程学报,2020,40(13):4301-4310.收稿日期:2020-10-19修改稿日期:2020-10-28


更多推荐

矢量,电流,预测