王建波1
(1.国网陕西省电力公司电力科学研究院摘
要
∗
刘毅力2
西安
梁继国2丁换换2
西安
710048)
710100)(2.西安工程大学电子信息学院
为了更好地解决电力系统稳定器(PSS)参数寻优的问题,将差分进化粒子群算法(DEPSO)应用到PSS参数优
化中。DEPSO算法通过PSO种群和DE种群之间的信息交流机制以寻求全局最优位置,从而使算法具有动态自适应性,能够较容易地跳出局部最优。文中首先阐述了在低频范围(0.1~2.5Hz)内用差分进化粒子群算法对PSS超前—滞后环节及增益参数进行寻优的方法。然后利用Matlab软件构建单机无穷大系统模型并进行仿真研究,结果表明利用该方法设计的PSS稳定性有了较大的提高,且效果优于基本PSO算法。采用该方法,基于机组实测无补偿相位的特性,对安徽电网某机组PSS参数进行了整定,验证了该方法的有效性。
关键词
电力系统稳定器;差分进化粒子群算法;单机无穷大;参数优化;参数整定
TP391
DOI:10.3969/j.issn.1672-9722.2018.01.013
中图分类号
OptimizationandSettingofthePSSParametersBasedonthe
DifferentialEvolutionaryParticleSwarmOptimizationAlgorithm
WANGJianbo1
LIUYili2
LIANGJiguo2
DINGHuanhuan2
710100)710048)
(1.ElectricPowerResearchInstituteofStateGridShaanxiElectricPowerCompany,Xian
(2.SchoolofElectronicandInformation,XianPolytechnicUniversity,Xian
Abstract
particleswarmoptimization(DEPSO)isappliedtoPSSparameteroptimization.DEPSOalgorithmcanseekglobaloptimallocation
Inordertobettersolvetheproblemsofpowersystemstabilizer(PSS)parameteroptimization,differentialevolution
throughinformationexchangemechanismbetweenPSOpopulationandDEpopulation,sothatthealgorithmhasdynamicadaptabili⁃byDEPSOalgorithminthelowfrequencyrange(0.1~2.5Hz).ThenconstructingasinglemachineinfinitesystemmodelusingMat⁃powergrid,basedontheunitmeasureduncompensatedphasecharacteristics,thealgorithmsvalidityistested.tion,settingparameter
ClassNumberKeyWords
tytoskiplocaloptimumeasily.First,thepaperexpatiatesontheoptimizationmethodofPSSlead-lagginglinkandgainparameterslabandisstudiedbysimulation,theresultsshowthatthePSSstabilitydesignedbyDEPSOalgorithmhasbeenimprovedgreatly,andtheeffectofDEPSOisbetterthanPSOalgorithms.Usingtheproposedmethod,adjustingPSSparametersofsomeunits,Anhui
PSS,differentialevolutionaryparticleswarmoptimization,singlemachineinfinitebus,Parameteroptimiza⁃TP391
1引言
随着电网互联规模的扩大,高放大倍数快速励
能力和电力系统的暂态稳定水平。但在这种条件下,系统受到干扰后,也可能由于阻尼不足或负阻尼引起发电机的转子角、转速,以及相关电气量发生近似等幅或增幅的频率较低的振荡。在大型电
磁技术的广泛采用,有效地提高了发电机电压调节
∗
收稿日期:2017年7月5日,修回日期:2017年8月26日基金项目:国家自然科学基金青年科学项目“硅橡胶材料表面水珠/水带运动规律及其对沿面放电特性影响机理”(编号:51707141);陕西省智能电网输变电设备状态监测技术创新团队(编号:2014KCT-16);西安市科技局技术转移促进工程(编号:CX12185⑤)资助。作者简介:王建波,男,硕士研究生,高级工程师,研究方向:电能质量、网源协调、新能源发电检测技术研究。刘毅力,男,副教授,硕士生导师,研究方向:电能质量、分布式光伏发电技术、电力系统自动化、电力电子与电力传动。梁继国,男,硕士研究生,研究方向:电能质量、分布式光伏发电技术、电力系统自动化、电力电子与电力传动。丁换换,女,硕士研究生,研究方向:电能质量、分布式光伏发电技术、电力系统自动化、电力电子与电力传动。
54王建波等:基于差分进化粒子群算法的PSS参数优化与整定第46卷
力系统中,低频振荡是一个很普遍的问题,它会引起联络线过流跳闸、系统与系统或机组与系统之间的失步而解列,严重威胁电力系统的稳定。目前易于实现且经济可靠的方法是采用电力系统稳定器(PowerSystemStabilizer,PSS)抑制低频振荡的阻尼力矩[1]。
统稳定性的改善具有重要作用。合理配置PSS作为一种附加的励磁控制装置对电力系
PSS的
参数可以取得理想的系统动态性能,提高系统的稳定性。目前,国内外学者在PSS参数优化方面做了一定的研究工作。采用多种优化算法如遗传算法、混沌算法、Swarm数优化,Optimization进化策略和基本粒子群优化起到了一定的优化效果
,PSO)算法等来实现对(Particle[2~5]
,但存在收敛性
PSS的参差、容易陷入局部极值等缺点。本文在经典PSS参数优化研究工作的基础上,提出采用差分进化粒子群Optimization优化(Differential优化。该算法作为一种全局智能寻优方法,,DEPSO)Evolution算法来实现对BasedParticlePSS参数进行Swarm不仅兼
顾了PSO算法收敛速度快,全局搜索能力强的特点,还拥有(DifferentialEvolution,DE)算法特有的记忆能力,使其可以动态跟踪当前的搜索情况,以调整其搜索策略,具有较强的全局收敛能力和鲁棒性。这种方法不同于以往只寻找机电振荡模式下阻尼比最小的PSS优化方法,而是可以根据最优控制原理综合考虑PSS与励磁系统的性能,优化的控制目标设为系统输出按最小误差跟踪给定值的能力,将PSS参数优化协调转化为带有不等式约束的优化问题。仿真和实验结果证明了该方法的有效性。
2
PSS数学模型及参数设计
2.1
励磁系统模型
文中发电机数学模型选用经典的六阶模型,励
磁系统数学模型选用IEEE-ST1型[1]。IEEE-ST1型励磁模型与PSS的结构匹配如图1所示。
图1
IEEE-ST1型励磁模型与PSS结构图
电压励磁系统可以表示为
E′Fd=T1[K(Vref-Vt+UPSS)-Efd]
(1)
A
式中,Vt和Vref分别代表机端电压和参考电压,
UPSS是稳定器信号,K和TA是瞬态增益以及励磁时间常数。
2.2PSS数学模型
典型的IEEEPSS1A控制器模型包括增益、隔
直环节以及两个超前滞后补偿环节,传递函数如式(2)所示:
U=K
sTω(1+sT1)(1+sT3)
1+sTDω
ω(1+sT2)(1+sT4)(2)
式中,U为发电机的PSS输出信号;Dω为发电机的转速偏差输入信号;Tω为隔直时间常数;T1、
T2、T3和T[6]4为补偿环节的时间常数。在本文研究中,T2和T4为与被控对象相关的时间常数;PSS增益K和时间常数T1、T3为待优化的参数。
2.3PSS优化的目标函数
为了增加系统阻尼,在整个低频振荡频段上,
引入两个基于目标函数的特征值,表达式如式(3)、式(4)所示:
J1=max{Re(λi)}(3)J2=min{ζi}
(4)
这里Re(λi)和ζi为第i个机电模式相关的特征值的实数部分和阻尼比。在优化的过程中,可以减小J1,使特征值向S平面左边移动,提高阻尼;也可以增加J2,即增加系统阻尼。所以,根据参数的边界特性,PSS设计转化为如下优化问题[7]:
ìï
OptimizeJïKPSSmin£KPSS£KíPSSmax
ïT(5)ï1min
£T1£T1maxîT
3min£T3£T3max
这里KPSS的取值范围为[0.01,50],T1、T3的典型取值范围是[0.06,1.0]。考虑到式(3)、式(4)中给出的目标函数,这里选用DEPSO算法来解决这种优化问题,并找到一组最适合的PSS参数用于实际数据的实验。
3
差分进化粒子群优化算法
3.1
PSO算法基本原理
tion粒子群优化算法PSO(ParticleSwarm的启发,)是美国的于1995Kenned年提出的一种生物进化算法。和Eberhar受鸟群觅食行为
Optimiza⁃
PSO算法采用速度—位置搜索模型。每个粒子代表解空间的一个候选解,解的优劣程度由适应函数决
定。每个粒子以一定的速度在解空间运动,并向自身历史最佳位置pbest和邻域历史最佳位置gbest聚集,实现对候选解的进化[8~9]。
2018年第1期计算机与数字工程55具体算法描述是,在D维解空间中,有s个粒子组成的粒子群,其中第i个粒子位置可以表示成D维向量,xˉi(n)=[xi1xi2xijxiD],j表示变
量xˉi的第j维分量;粒子飞行的速度为
vˉi(n)=[vi1vi2vijviD];
每一次迭代,粒子通过动态跟踪两个极值来更新其速度和位置。第一个是粒子从初始到当前迭代次数搜索产生的最优解:个体极值pi=(pi1pi2piD)。第二个是粒子种群目前的最优解:全局极值gi=(gi1gi2giD)。粒子根据公式(6~7)来更新其速度和位置:vi(t+1)=ω vi(t)+c1Rand() (pi-xi(t))
+c2Rand() (gi-xi(t))
(6)xi(t+1)=xi(t)+vi(t+1)
(7)
其中ω为惯性权值;Rand()在[0,1]范围内变化的随机数;一般学习因子取c1=c2=2;粒子数i=1,2,…,s。粒子在解空间内不断跟踪个体极值
与全局极值进行搜索,直到达到规定的迭代次数或满足规定的误差标准为止。粒子在每一维飞行的速度不能超过算法设定的最大速度vmax。设置较大的vmax可以保证粒子种群的全局搜索能力,
vmax较小则粒子种群的局部搜索能力加强[10]。
3.2
演化算法一样,DE差分进化算法基本原理
算法是由DEStorn是一种模拟生物进化的随机模等于1995年提出的,和其它
型,通过反复迭代,使得那些适应环境的个体被保存了下来。其基本思想是:从某一随机产生的初始群体开始,按照变异(Mutation)、交叉(Crossover)和选择(Selection)的操作规则不断迭代计算,并根据每一个体的适应度值,保留优良个体,淘汰劣质个体,引导搜索过程向最优解逼近[11]。3种操作描述如下:
变异个体
1)变异操作:对于个体x(t),根据式(8)生成
xi(t+1)=xr3
(t)+F((xr1
(t)-xr2
(t)
)(8)
其中,随机数r1r2r3Î{12N},F为缩放比例因子,用于控制差向量的影响大小。
体x(2t)交叉操作:)和变异结果进行二项分布杂交产生新种群
可以增加种群的多样性。由个
x′
i(t+1)=(x′′′
i1(t+1)xij(t+1)xid(t+1))。
具体操作如下:
x′+1)=ìíxij(t+1) if b(j)£CR or j= r(i)
ij(tîx(9)
ij
(t) if b(j)³CR or j¹ r(i)其中,jÎ{12d},0£b(j)£1是[0,1]之间的随机数,0£CR£1是变异概率,r(i)Î{12d}是随机选择指数。
子代个体优于父代个体时,3)选择操作:采取贪婪策略,即f(即只有当产生的
x′i(t+1))£f(xi(t))时才被保留,否则父代个体被保留至下一代。
3.3
区域。可是,PSO基于收敛速度比较快,DEPSO算法的PSS当全局历史最优值陷入局部最优区域也可以较好地探索求解
参数优化
时,群体会迅速地收敛到该区域,使得种群多样性消失,出现早熟停滞的现象。DE的变异算子有利于增加全局搜索能力,保证种群的多样性;交叉算子可以提高局部搜索能力,加快收敛速度;选择算子具有一定的记忆能力,能够保留优秀个体。为了充分结合两种算法的优势,将二者有机结合,按粒子搜索方式不同将整个粒子群体分为两个新的分群,分别命名为PSO分群(P群)和DE分群
D群)[12~16]。
设Pbest(t)、Pgbest(t)分别代表P群t时刻个体历史极值解和全局历史极值解,Dgbest(t)代表t时刻D群全局历史极值解,gbest(t)代表t时刻整群全局历史极值解。DEPSO混合优化算法流程为:
步骤1初始化每个分群群体规模M,最大迭代次数n,求解精度ε,最大惯性权重ωmax,最小惯性权重ωmin,控制因子k,加速因子c1和c2,缩放因子σ,变异概率CR。
步骤2评价P群的适应度fitness,用每个粒子对应当前系统的特征值,根据式(3)、式(4)求出目标函数J,即特征值最小,阻尼比最大。
步骤3对D群体中所有个体执行差分进化算法的变异、交叉、选择操作。其中,进行选择操作时,选择的标准是粒子目标函数J的适应值最优。
步骤4按照适应度fitness选择P群中的最佳个体Pgbest(t)。
步骤5按照适应度fitness选择D群中的最佳个体Dgbest(t)。
步骤6比较Pgbest(t)和Dgbest(t)优劣,选择最佳个体作为P群和D群下一代进化依据。
步骤7保留当前整个群体中最佳个体gbestt),如果满足精度要求或整个进化已达到最大迭代次数,则终止算法,输出处于最优位置的粒子,对该粒子进行解码生成2。
PSS模型参数;否则转到步骤((56王建波等:基于差分进化粒子群算法的PSS参数优化与整定第46卷
4
4.1
仿真及实验
Simulink建模与仿真
本文基于单机无穷大系统进行参数优化,使用
'\"
0.264,Xd\"=0.225,Xq\"=0.248,Tdo=8.86s,Tdo=
f=50Hz,H=3.7s,Xd=2.61,Xq=2.48,Xd'=
V=27KV,发电机主要参数:Pn=1110MVA,
Matlab的电力系统仿真模块集SimPowerSystemsBlockset搭建仿真模型如图3所示。模型具体参数包括:
Out1
In1
0.036s,Tq\"=0.085s,rs=0.001。
励磁调节器主要参数:KAVR=300,TA=0.03。输电线参数:Xl=0.3,R=0.026。
PSS
wrefPrefPmwePe0 gatedwHTGvrefvd Vfvqvstab dw(pu) Continuous 10.7516 PmVf_ mAB VaABC Volts>pu ac Va(pu) ABCsource 1.0 C SynchronousMachine b Three-phaseTransformer ABCABCABCSpeed(pu) ExcitationSystem Vf(pu)3-PhaseFault 图2Simulink仿真模型 RotorSpeedDeviationDω(pu)0.0150.0050.020.010 DEPSO WithoutPSS T4=0.02;Tω=10。需优化的参数为K、T1、T3。通 PSS采用典型IEEEPSS1A数学模型,其中T2= 过采用DEPSO计算的最终优化结果为:K=6.63,得到转速偏差随时间变化的仿真图。如图3所示有加PSS时的转速偏差Dω随时间变化的比较; RotorSpeedDeviationDω(pu)0.0150.0050.0100.02 WithoutPSS PSO T1=0.277s,T3=0.161s。将优化结果带入仿真模型是系统采用PSO算法优化参数后的PSS与系统没 -0.005-0.015-0.025 -0.01-0.02 0 1 2 3 456 Time(s) 7 8 9 10 图4基于DEPSO优化的发电机Dω-t仿真曲线 4.2实验与分析 针对PSS参数优化设计的目的和要求,开发了 -0.005-0.015-0.025-0.02-0.01 一套基于DEPSO算法的PSS参数优化整定软件,软件内嵌了可仿真分析PSS参数优化效果的励磁系 0 1 2 3 456 Time(s) 7 8 9 10 统仿真模型和优化PSS参数的DEPSO算法。利用该软件,一方面可以根据不同的发电机组模型在系统中填写电机模型等相关参数,通过DEPSO算法获得PSS的优化参数,进行4.1节所述的仿真实验;另一方面,软件可以根据现场数据进行真机PSS的参数优化。即采用频谱仪实际测得电厂发电机组励磁系统的频率特性和变送器的频率特性,软件提取频率特性数据后将这些数据分别存入自定义的向量中,画出励磁系统的相频特性曲线和变送器的 图3基于PSO优化的发电机Dω-t仿真曲线 图4是系统采用DEPSO算法优化参数后的PSS与系统没有加PSS时的转速偏差Dω随时间变化的比较。图3中采用了PSO优化参数的PSS明显能够较快地平息低频振荡。而由图4可见,采用DEPSO算法的优化结果优于采用PSO算法的优化结果,能更快更稳地消除低频振荡。 2018年第1期计算机与数字工程57相频特性曲线,根据频率特性求解PSS目标函数,利用DEPSO算法进行优化计算,获得优化的PSS参数,并根据“励磁系统+PSS”的相频特性判断参数的优劣。 因此,实验的目的是检验在低频范围内,励磁系统通过附加采用DEPSO算法优化参数后的PSS,能否将原本系统偏差较大的频率特性快速、稳定的补偿到规定的范围内,从而更快更稳地消除低频振荡。实验选用安徽铜陵某电厂5号发电机组作为对象,通过频谱仪采集励磁系统及变送器频率特性数据。软件读取全部数据的采样点有256个,低频率范围在0~3Hz之间的采样点有100个,等均匀提取其中31个采样数据作为软件计算PSS参数的原始参量,画出相频特性如图5(a)、(b)所示。利用差分进化粒子群优化产生的参数组建PSS模型,并画出PSS的频率特性曲线如图5(c)。励磁系统相位与变送器相位的差值再加上PSS补偿环节的相位即为经PSS矫正后的系统相位,其频率特性曲线如图5(d)所示。 -50位相-100-150-200 0.5 1 1.5频率 2 2.5 3 -40(a)励磁系统相频特性 位-60相-100-80 -120-1400.5 1 1.5频率 2 2.5 3 20(b)励磁系统-变送器的相频特性 位相0-20-400.51 1.5频率 2 2.53 -80(c)PSS相频特性 位-90 相-100-110-120-1300.51 1.5频率 2 2.53 (d)励磁系统-变送器+PSS的相频特性 图5实际试验数据频率特性 对于输入量为转角偏差的PSS模型而言,其指标要求补偿后的相位应该在-80°~-145°之间。将所有的计算参数放入数据表中进行观察。这里,为了结果更加直观明了,将低频范围内的31组数据等均匀地选取11组予以表示,如表1所示。 表1 实际试验数据相频特性列表 序频率励磁系统励磁系统号(Hz)(度)VT器(变送度)PSS(度)120.27470.5493-64.8677-7.5265-14.6586-38.3534PSS-27.6511-95.6947(度)+3-102.961940.8240-89.969451.0986-104.5775-101.951361.3733-125.4950-21.250671.6479-140.7847-27.3586-14.2560-3.8148-97.5830-91.1023-104.382281.9226-137.2029-32.84023.5622-91.803892.1973-145.4522-37.58228.5184-90.3695102.4719-149.5936-41.968911.6795-90.5241112.7466-153.9648-45.69643.0212-153.5771-49.027413.5276-157.8047-51.761314.4133-54.353214.588814.2366-87.2270在结果中,“励磁系统+PSS(度)”这一列的矫 -89.2149正数据均在-80°~-145°之间,满足指标要求。由此可见,PSS采用差分进化粒子群算法优化参数后的实际问题。 能够有效的抑制发电系统中的低频振荡,解决5结语 将粒子群优化算法与差分进化相结合的差分 进化粒子群优化(DEPSO)算法,可以解决基本PSO算法早熟停滞现象,具有更快的收敛速度和全局搜索能力。 本文首先利用DEPSO算法来实现对PSS参数 优化,即优化计算PSS的增益及超前—滞后环节的时间常数。然后在Simulink环境下搭建单机无穷大系统仿真模型,将DEPSO算法优化得到的参数与基本PSO算法优化获得的参数进行仿真对比,结果表明使用DEPSO算法优化得到的PSS参数相对于基本PSO算法可以更加有效地获得PSS参数优 化问题全局最优解,即采用DEPSO算法优化参数后的PSS能更快更稳地消除低频振荡。最后,开发了基于DEPSO算法的用于PSS参数整定的优化软件,通过采集电厂发变机组真实数据来对参数优化结果进行实验验证,结果表明该方法是有效的。 参考文献 1]刘取.电力系统稳定性及发电机励磁控制[M].北京:中 国电力出版社,2007. [58王建波等:基于差分进化粒子群算法的PSS参数优化与整定 第46卷 LIUtionQu.controlStability[M].ofBeijingthepower:ChinasystemElectricandgeneratorPowerexcitaPress,⁃[2]王明东,2003. 刘宪林,于继来.基于改进遗传算法的发电机 调速器侧模糊PSS优化设计[J].继电器,2006,34(13):25-28. WANGsignoffuzzyMingdongPSS,forLIUgovernorXianlinof,YUgeneratorJilai.Optimizationde⁃provedGA[J].Relay,2006,34(13):25-28based.onanim⁃[3]陈刚,简华阳,龚啸.自适应混沌粒子群算法在PSS设计 中的应用[J].电力系统及其自动化学报,2012,24(4):82-87. CHENerGang,JIANHuayang,GONGXiao.DesignofPow⁃mOptimizationSystemStabilizersUsingAdaptiveProceedingsChaosofParticleSwartheCSU-EP⁃SA,2012,24(4)Algorithm:82-87.[J].⁃[4]牛振勇,杜正春,方万良,等.基于进化策略的多机系统 22-27. PSS参数优化[J].中国电机工程学报,2004,24(2):NIUParameterZhenyong,DUZhengchun,FANGWanliang,etal.StabilizersOptimizationUsingEvolutionaryofMulti-Machine].ProceedingsPowerSystemtheCSEE,2004,24(2):22-27.Strategy[Jof[5]吴峰,陈维荣,李奇,等.基于粒子群优化算法的PSS参 数优化[J].电力系统保护与控制,2009,37(10):53-58. WUmizationFengof,CHENGpowersystemWeirongJ].Powerstabilizer,LIQi,SystemonetProtectionparticleal.Parameterswamoptiandopti⁃Con⁃trolmization,2009algorithm,37(10)[:53-58. ⁃[6]朱方,方思力.电力系统稳定器的原理及应用[M].北 京:中国电力出版社,1996. ZHUPowerFangPress,1996. System,FANGStabilitySili.[TheM].PrinciplesBeijing:ChinaandApplicationsElectricPowerof[7]M.A.Abido.OptimalDesignofPowerSystemStabilizers UsingonEnergyParticleConversionSwarm,OptimizationVol.17,NO.[3.J].2002.IEEETransactions[8]崔志华,曾建潮.微粒群优化算法[M].北京:科学出版 社,2011. CUIZhihua,ZENGJianchao.ParticleSwarmOptimization[M].Beijing:SciencePress,2011. [9]KennedyJ,EberhartRC.ParticleSwarmOptimization [A].In:Proc.ofthe1995IEEEInternationalConference[10]祁万春,onNeural房鑫炎Networks.[基于混合粒子群优化算法的C].1995:1942-1984. PSS参数 优化[J].继电器,2005,33(13):21-24. QIpowerWanchunmizationsystem,algorithmstabilizerFANGXinyan.[J].Relayusing,hybridParameters2005,33particlesoptimizationof(13):21-24.warmopti⁃[11]唐俊,胡志敏,王敏.差分进化粒子群算法求解发电资 源调度问题[J].计算机工程与应用,2012,48(8):239-242. TANGtionulingparticleJun,HUtions,problemswarmZhimin2012,48[(J8]).algorithm,WANG:239-242. ComputerforMin.EngineeringpowerDifferentialresourcesandApplicaschedevolu⁃⁃⁃[12]杨妍,陈如清,俞金寿.差分进化粒子群混合优化算法 的研究与应用[J].计算机工程与应用,2010,46(25):238-241. YANGentialoptimizationevolutionYan,CHENalgorithmparticleRuqingandswarm,YUitsapplicationoptimizationJinshou.Study[Jbasedon].Computerhybriddiffer⁃[13]荣雅君,Engineering贾艳,and刘琳,Applications等.基于粒子群优化算法的电力 ,2010,46(25):238-241.系统稳定器设计[J].继电器,2008,36(4):24-26.RONGsystemYajunstabilizer,JIAYan,LIULin,etal.Designofpowerlay,2008,36(4):on24-26. particleswarmoptimization[J].Re⁃[14]倪以信,陈寿孙,张宝霖.动态电力系统的理论和分析 [M].北京:清华大学出版社,2002. NIAnalysisYixin,ofCHENDynamicShousunPower,ZHANGSystem[Baolin.M].BeijingTheory:TsingandhuaUniversityPress,2002. ⁃[15]PrabhaKundur,JohnPaserba,VenkatAjjarapu,etc. [DefinitionJ].IEEEandTrans.ClassificationonPowerSystemsofPower,Vol.19System,pp.1390Stability,[16]2004K.W.Wang.. Multimachineeigenvaluesensitivitiesof powerSystemssystem,2000,15parameters(2):741-747.[J].IEEETrans.onPower 因篇幅问题不能全部显示,请点此查看更多更全内容