基于改进混合粒子群优化的特定谐波消除脉宽调制策略-陈炜.pdf

《基于改进混合粒子群优化的特定谐波消除脉宽调制策略-陈炜.pdf》由会员分享，可在线阅读，更多相关《基于改进混合粒子群优化的特定谐波消除脉宽调制策略-陈炜.pdf（8页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第49卷 第5期 2016年5月 天津大学学报(自然科学与工程技术版) Journal of Tianjin University(Science and Technology) Vol.49 No.5May 2016 收稿日期：2014-09-30；修回日期：2015-01-19. 基金项目：国家重点基础研究发展计划(973 计划)资助项目(2013CB035600)；国家自然科学基金资助项目(51577134)；天津市应用基础研究计划重点资助项目(13JCZDJC34700). 作者简介：陈 炜（1977 ），男，博士，副教授 通讯作者：陈 炜，chen_. 网络出版时间：2015-09-

2、18. 网络出版地址：http:/ DOI:10.11784/tdxbz201409079 基于改进混合粒子群优化的特定谐波消除脉宽调制策略 陈 炜1，陈玉健1，耿 强2(1. 天津大学电气与自动化工程学院，天津 300072； 2. 天津工业大学电工电能新技术天津市重点实验室，天津 300387) 摘 要：针对传统迭代法求解特定谐波消除脉宽调制策略(SHEPWM)开关角方程组需要合适的初值和难收敛的问题，提出一种改进的混合粒子群优化(HPSO)算法对其进行求解该算法对粒子群(PSO)算法的权重系数和学习因子进行了改进，并且提出了一种温度系数线性递减的模拟退火(SA)算法与粒子群算法结合，有效

3、弥补了传统粒子群(CPSO)算法求解开关角收敛速度慢和精度低的缺点仿真分析表明，该算法消除了对初值的依赖，提升了算法寻优的能力，从而提高了求解速度与精度，并且通过实验验证了该算法的可行性 关键词：特定谐波消除；优化算法；模拟退火；混合粒子群优化 中图分类号：TM464 文献标志码：A 文章编号：0493-2137(2016)05-0498-08 Selective Harmonic Elimination Pulse Width Modulation Based on Improved Hybrid Particle Swarm Optimization Chen Wei1，Chen Yuji

4、an1，Geng Qiang2(1School of Electrical Engineering and Automation，Tianjin University，Tianjin 300072，China； 2Tianjin Key Laboratory of Advanced Electrical Engineering and Energy Technology， Tianjin Polytechnic University，Tianjin 300387，China) Abstract：The requirement of an initial value and the speed

5、of convergence are the main challenges for the traditional iteration method in selective harmonic elimination pulse width modulation(SHEPWM)As to these problems，an improved hybrid particle swarm optimization(HPSO)algorithm was proposed in this paperIn the proposed strategy，the inertia weight and lea

6、rning factors of PSO algorithm was adjustedMeanwhile，a new kind of simulated annealing(SA)algorithm，whose temperature coefficient decreases with the increase of the number of iterations，is put for-ward to combine with PSO，making a favorable improvement in terms of the convergence speed and accuracy

7、com-pared with the conventional PSO(CPSO)algorithmFinally，the viability and performance of the proposed strategyare shown through simulation and experimental results in a laboratory prototype. Keywords：selective harmonic elimination；optimization algorithm；simulated annealing(SA)；hybrid particle swar

8、m optimization(HPSO) 在机车牵引、风力发电等大功率系统中，受到开关损耗和散热的限制，功率器件的开关频率较低但是，传统的调制策略在开关频率与输出电压频率之比较小的情况下会产生大量的低次谐波，严重影响系统的运行性能1 特定谐波消除脉宽调制策略(selective harmonic elimination pulse width modulation，SHEPWM)在功率器件开关频率较低的情况下能够有效消除低次谐波，提高逆变器输出的波形质量和直流电压利用率，适合作为大功率逆变器系统的调制策略1但是，该策略的开关角求解方程是傅里叶分析产生的非线性超越方程组，求解困难，所以开关角的求

9、解成为学者研究 2016年5月 陈 炜等：基于改进混合粒子群优化的特定谐波消除脉宽调制策略 499 SHEPWM的重点和难点 为了求解方程，应用较多的是牛顿迭代法，该方法在初g984合适的情况下能够有效求解开关角2-3g2976g40884g909用g1358g863g4778法g2090迭代法g6408行g2951g6408，g1816g2369了迭代g4778法g2949g2964的能力g2976g40885提出了g1778g773 Walsh 方程的求解方法，g6923低了g4778法g2090初g984的g5941求g700g6423方法求解g3024g6543g6985g5941合适

10、的初g984，g2290g708g6947g4427求解开关角的g1816多，g4778法g1836g3133性g1816g1194，g1194大了求解的难g2318g6447g6398研究发g4119，开关角的求解g6855g7076g1338以g6336g1232成为g3107g853g1232g6855g7076g3150解g1092，g6407有利g773g2891g5233求解初g984的g3147g4999，g2409到g3101g853求解方法6-7g2976g40888和g2976g40889分g1150g6611用了g3720合g6475传g4778法和传统g4833g200

11、4g5061g853g1232(conventional particle swarm optimi-zation，CPSO)算法g2090电压g1194g3134g2483谐波g4249变率和电压g2483谐波g4249变率g3107小g3024的开关角g6408行g853g1232求解，g1106小了输出电压的g2483谐波g4249变率g2976g408810用g5709g5061g4778法g2090开关角方程g6408行求解，有效提高了求解g4852g2318g2976g408811提出了g5763g5061g4778法用g773多电g2288逆变器开关角求解g815g700g642

12、3研究g1338以发g4119，g5061g853g1232g4778法适合g773求解 SHEPWM 开关角方程 在g5061g853g1232g4778法中，g4833g2004g5061g853g1232g4778法g1665g1060g2471g2566直g5945、g2051g4119g4772g1260g2290g708g1061有g2401高的g2671行效率成为研究热点，但g4833g2004g5061g4778法g3035g6941g1049g2124g6537g3107g853解12-13为提高g4778法的性能，g2976g408814g6611用g6947迭代次g297

13、2线性g6443g1106的g2564性g3134重，g909g2409g4778法g1169g3119g1061有较高的g1050g2124g2876g4873能力，g1360g3119g1061有较高的g2124g6537g2876g4873能力，g1194g2443了g2949g2964g6451g2318g7755在g3428g1778g4520g700，g2976g408815g6611用g3380g4857调g2974g3134重法g1198g2454调g5355g1050g2124和g2124g6537g2876g4873能力g2976g408816g2097g3380g2737

14、g6429g3873(simulated an-nealing，SA)g4778法引g1049到了g4833g2004g5061g4778法中，g909g2876g4873g6398程g1061有g3343率g6265变的能力，g1816g2369了g4778法g2092g853的能力 g3124g2976提出g694g4631g2951g6408的g3720合g4833g2004g5061g853g1232(hybrid particle swarm optimization，HPSO)g4778法，g2290g2097g1060应用g773 SHEPWM 开关角的求解该g4778法g209

15、0g4833g2004g5061g4778法的g3134重系g2972和学g759g1665g2004g6408行g2951g6408，g2290g708g2097g3741g2318系g2972g6947迭代次g2972g6443g1106的 SA g4778法与 PSO g4778法g4934合g3150提g1251g4778法g2092g853的能力g3124g2976g2090g2951g6408g4778法的性能和该g4778法求解的g4934g3174g6408行分析，g6447g6398g2051g7181g7181g6012求解g4934g3174的有效性 1 SHEPWM开关

16、角求解方程 g2090g773g719电g2288g699g4400逆变器，g2374g6611用g1316g3164性 PWM 调制方g2351g3024，产生的输出g4400电压有g3427电g2288开始和负电g2288开始g719g4631形g2351，分g1150称为 a 类和 b 类，g1060中 a 类波形在 1/4 周g3119内开关角的个g2972 N 为偶g2972，b 类为奇g2972图 1 为 a g4400g1316g3164性调制 PWM 输出波形，a 类以N2为例，b类以N3为例图1所示1、2、3为开关器件切换g3024刻的g4400位角，g4772称为开关角 （

17、a）a类(N2) （b）b类(N3) 图1 双极性调制PWM输出波形 Fig.1 Output waveforms of bipolar modulation PWM g2090逆变器输出的g4400电压波形g6408行傅里叶级g2972展开，g1060表达g2351为 a1() sin( ) cos( )nnnut A nt B nt=+(1)g2351中：n为谐波次g2972g7755为g1778波角频率 由g773输出电压存在 1/4 周g3119g2090称性和半波g2090称性，所以 Bn0，Ang1338以用 1/4 周g3119中的开关角 1，2，N表示，即 dc121 2 (

18、1) cos( )0NlnllnAU nnB= + =(2)g2351中：l 为开关角编号g7755N 为 1/4 周g3119内开关角的个g2972g7755Udc为逆变器输g1049侧的直流电压g7755“”用g773 a类波形的计g4778g7755“”用g773b类波形的计g4778 定义调制比为 1dc2UmU= (3)g2351中U1为参考电压 在g699g4400逆变器中，由g773波形存在g699g4400g2090称、半波g2090称和1/4周g3119g2090称，在输出的线电压中3及3的倍g2972次谐波和偶g2972次谐波为零令g1778波幅g984等g773参考电压，

19、指定次g2972的谐波幅g984为零，即g1338g2409到SHEPWM的开关角g2318求解方程组为 1141 2 ( 1) cos 12 (1)cos( )0NlllNlllmq=+ =+ =(4)g2351中 q为g5941消除谐波的次g2972，q5，7，3,g1(g为偶g2972)求解N个开关角g6985g5941N个方程，g2351(4)中只g6985取 N1 个 q g984即g1338，即g3428g3024能消除 N1 个不g1358 500 天津大学学报(自然科学与工程技术版) 第49卷 第5期 次g2972的谐波 g2351(4)为多变量的非线性超越方程组，g2290g

20、708存在多组解，g1060求解比较困难为求解g3428方程组，g1338g2097g1060g6336换成为g3107g853g1232g6855g7076，即 1151122 2115412 (1)cos12 (1)cos(5)12 (1)cos( )(, )( )NlllNlllNlqllN qCCCqFCmCC= = + =+ =+ = +vertellipsismidhorizellipsismidhorizellipsis(5)g2351中012N/2 求解 F(1，N)的g3107小g984，该g3107小g984所g2090应的开关角就是g2351(4)所表示的开关角方程组的解

21、g3124g2976g2097 F(1，N)定义为g4778法的适应g2318函g2972，g6611用g2951g6408的g3720合g4833g2004g5061g853g1232g4778法g6408行g3107g853g1232求解 2 基于改进混合粒子群算法的开关角求解 2.1 自适应调节的权重系数和学习因子 在 CPSO 中，每个g4833g2004g6543有位置、g6451g2318和由目标函g2972g1092定的适应g2318g984，g2290根据“自身经g7181”和“g4631g5061经g7181”在g2876g4873空间中向g3101好的位置移g1198，g2

22、876g4873g1050g2124g3107g853解CPSO中的各系g2972均为常g2972，g1060g4778法g3101新方程为 -1 111 B122 B()()ikk kii ikiwcrcr=+ +vv PXGX(6)1kkkiii=+XX v (7)g2351中：k 为g2374g1169迭代次g2972g7755i 为g4833g2004编号(i1，2，s，s 为g5061体规g3380)g7755kiX 为第 i 个g4833g2004的位置，T12, , kkk kiii iNxx x=XT12, , kk kii iN ，即第 i 组开关角的g984g7755kiv

23、 为第 i 个g4833g2004的g6451g2318，即第 i 组开关角的g6451g2318，T12, , kkk kiii iNvv v=v g7755BiP 为第 i 个g4833g2004的g3107g853点，即第 i 组开关角自身迭代的g3107g853g984，1Bb,iip=P 2Tbb, iiNpp g7755BG 即 s 组开关角中的g3107g853g984，为g1050g2124g3107g853点，12TBbb b, , Ngg g=G g7755w为g2564性g3134重g7755c1和c2为学g759g1665g2004g7755r1和r2为(0，1)之间的

24、g6947机g2972 在 CPSO 中，由g773g2564性g3134重w和学g759g1665g2004c1、c2g6543是常g2972，g4778法不能根据g1061体情况g1198g2454调g2974，g2092g853的性能下g6923，所以g6985g5941g4778法的参g2972 w、c1、c2g6947迭代次g2972的g1816g1194g1198g2454调g2974 在求解 SHEPWM 方程组g3024，根据方程组的特点，g2374开关角的适应g2318g98412(, , )kk kii iNF midhorizellipsis 较大g3024，g2564性

25、g3134重系g2972w 应该较大，以保g6012g2876g4873的g1050g2124性g7755而g2374开关角的适应g2318g984较小g3024，g2564性g3134重系g2972w 应该较小，以保g6012在该g4833g2004处g2876g4873的g4852确性g1665g3428，g2564性g3134重的选取应根据开关角的适应g2318g984大小g6408行自适应选择 为了g1198g2454地调g5355g4778法的g2124g6537和g1050g2124g2876g4873的能力，定义g2564性g3134重系g2972为 max min minmin

26、 avgavg minmax avg()()()iiiwwFFwFFFFw+=(8)g2351中：wmin和 wmax分g1150为g2564性g3134重的g3107小g984和g3107大g984g7755Fi为第 i 个g4833g2004的适应g2318g984g7755Favg为g2374g1169所有g4833g2004适应g2318的g2288均g984g7755Fmin为g2374g1169所有g4833g2004适应g2318的g3107小g984 为了g6408g694步提高g4778法性能，g2090学g759g1665g2004也g6408行g2951g6408，g64

27、47g6398g909g719个g1665g2004分g1150变g1232，g909g4778法在g853g1232的初始阶段g1061有较大的自我学g759能力和较小的社会学g759能力，g1194g2369g1050g2124g2876g4873能力g7755而在g1360g3119则刚好g4400反，有利g773g2949g2964到g1050g2124g3107g853解 为g3428，定义g719个学g759g1665g2004分g1150为 11 10110maxcccc kk=+ (9)21 20220maxcccc kk=+ (10)g2351中：c10和 c20分g1150

28、为 c1和 c2的迭代初g984g7755c11和 c21分g1150为 c1和 c2的迭代终g984g7755k 和 kmax分g1150表示g2374g1169迭代次g2972和g3107大迭代次g2972 2.2 改进的模拟退火混合粒子群算法 CPSO g4778法g3124身不g1061备g4833g2004g6265变的能力，g909g4833g2004容g3035在g2124g6537g3107g853解处停滞g3124g2976g2097g3380g2737g6429g3873g2471g2566引g1049到g4833g2004g5061g4778法中，g909g2409g47

29、78法在g2974个g2876g4873g6398程中g1061有g3343率g6265变的能力传统的g6429g3873方g2351g3741g2318系g2972设置为常g2972，g3124g2976提出了g3741g2318系g2972g6947迭代次g2972的g1816g1194而线性g6443g1106的g6429g3873策略，保g6012迭代开始g3024，g3741g2318系g2972较大，保g6012g2876g4873的g1050g2124性g7755而迭代g1360g3119g3741g2318系g2972较小，缩短g2949g2964g3024间 假设g3380g

30、2737g6429g3873g4778法的初始g3741g2318系g2972和g3107终g3741g2318系g2972分g1150为t0和t1，则第k次迭代g3024的g3741g2318g1338表示为 01max11max()( )kkttk kTtTk=+ (11)g2374g1169g3741g2318下突变的g3343率g1338以表示为 /ekFTP= (12)g2351中GP_BestFF F= ，GB()FF= G 表示g1050g2124g3107g853g4833g2004的适应g2318g984，P_BestF 表示 s 个个体g3107g853g4833g2004

31、的g3107小适 2016年5月 陈 炜等：基于改进混合粒子群优化的特定谐波消除脉宽调制策略 501 应g2318g984，12P_Best P P Pmin , , , sFFFF=，PB()iiFF= P 表示第i个个体g3107g853g4833g2004的适应g2318g984 g1360g2976中g2097g2951g6408的自适应g3134重系g2972和学g759g1665g2004的g3380g2737g6429g3873g4833g2004g5061g4778法g4772称为g2951g6408的g3720合g4833g2004g5061g4778法图2为g3124g29

32、76提出的g2951g6408g3720合g4833g2004g5061g4778法的流程 迭代开始g3024，g6985赋g984的参g2972有g3107大迭代次g2972kmax，g4631g5061规g3380 s，g3134重系g2972 wmin和 wmax，学g759g1665g2004c10、c11、c20、c21，g3741g2318系g2972t0和t1为了g909g6429g3873g4778法g815初始g3741g2318逐渐下g6923，g3741g2318系g2972应在(0，1)范围内，设置t01，g909g4778法g815初始g3741g2318开始g692

33、3g3741，设置t1110-3g909g3741g2318系g2972线性g6443g1106，提g1251g6923g3741g6398程的g6451g2318设初始g3741g2318T01赋g984给g4833g2004g3164限g6451g2318vmax0.05, 求取初g984的步骤为：首先在(0，/2)之间g6947机选取开关角0iX xi1，xi2，xiNg7755在(vmax，vmax)之间g6947机选取初始g6451g23180iv vi1，vi2，viN，i1，2，s然g1360，根据g2351(5)计g4778出各组初始开关角的适应g2318g984 FiF(0i

34、X )，g2290g708根据定义计g4778初始个体g3107g853g4833g2004为00Bii=PXg7755初始个体g3107g853g4833g2004的适应g2318g984为0PB()iiFF= P g7755初始个体g3107g853g4833g2004的适应g2318g3107小g984为120000P_Best P P Pmin , , , sFFFF=g7755初始个体g3107g853g4833g2004中适应g2318g984g3107小的g4833g2004为0Bestarcmin=X1200PP, ,FF0PsF . g3107g1360，求解初始g1050g

35、2124g3107g853g4833g2004的适应g2318g984为 FG0P_BestF g7755初始g1050g2124g3107g853g4833g2004为GB0BestX 在g2951g6408的HPSOg4778法的迭代g6398程中，g2090g773第i个g4833g2004，如g3174g2374g1169迭代次g2972下的g984kiX g853g773个体g3107g853g9841BikP g3024，就g909kiX 替代1BikP ，即个体g3107g853的适应g2318g9841PPmin ( ), iikkiFFF= X g7755个体g3107g85

36、3Barcmin ( ),ikkiF=PX 1P)ikF，表示kiX 和1BikP 中适应g2318较小的那个g3428g3024，11P_Best P P Pmin , , , skkkkFFFF=，表示第k次迭代g3024s个个体g3107g853g4833g2004的g3107小适应g2318g984，1Best Parcmin ,kkF=X 1PP, skkFF 表示第k次迭代g3024s个个体g3107g853g4833g2004中适应g2318g984g3107小的g4833g2004 在g3101新g1050g2124g3107g853 GBg3024，g2374计g4778g2

37、409到的 F0g3024，以g3343率1接受GBP_BestkFBestkX ，FGP_BestkF g7755而g2374F0 g3024，以g3343率 P 接受GBBestkX ，FGP_BestkF g7755g3428g3024g1050g2124g3107g853g4833g2004以g3343率 P 发生g6265变，g909g4778法g1061有g694定接受较差解的能力，提高了g4778法g1050g2124g2876g4873的能力 g6611用g2951g6408的g6429g3873g4778法不但能够提高g2951g6408g4778法的g2671行g6451g

38、2318，g1358g3024g6947g4427迭代次g2972g1816g1194，P 逐渐g1106小，在迭代初g3119 P g984较大，保g6012了在g2876g4873初g3119g4833g2004以较大g3343率接受差解，有利g773g6265出g2124g6537g3107g853g7755到迭代g1360g3119，Pg984较小，g4778法接受较差解的g3343率较小，g4833g2004能在某个g2124g6537空间g6408行细致g2876g4873 图2 改进混合粒子群算法优化流程 Fig.2 Flow chart of the improved HPSO

39、 algorithm 2.3 改进的HPSO算法的求解结果 g3124g2976g6611用g2951g6408的g3720合g4833g2004g5061g4778法分g1150g2090 1/4 周g3119含有 2、3、5、7、9 个开关角的情况g6408行求解，g1060中既包含 a 类波形(N2)，也包含 b 类波形(N3，5，7，9)每g4631情况下求g2409的解各g719组，g4934g3174如图 3 所 502 天津大学学报(自然科学与工程技术版) 第49卷 第5期 （a）N2 （b）N3 （c）N5 （d）N7 （e）N9 图3 改进混合粒子群优化算法求解结果 Fig.

40、3 Results of the improved HPSO algorithm 示，g2051线为第1组，虚线为第2组 g2090g773g3124g2976求解的 SHEPWM g853g1232g6855g7076，g3720和g4833g2004g5061g4778法求g2409的g719组开关角均是g3107g853解，均能消除指定次谐波，达到消谐目的消除指定次g2972的谐波之g1360，还存在的主g5941谐波含量是g3107低的g719次不为零的谐波. 为了研究调制策略的电压谐波特性，g1338以g909用电压谐波g4249变g1665g2972(harmonic distor

41、tion factor，HDF)6g3150g6408行评g840，定义g3428g1665g2972为 2231 3512232 341HDFNNNNUUNUUUNU+=+为偶数为奇数(13)g2351中 U3N1、U3N5、U3N2、U3N4分g1150为 3N1、3N5、3N2、3N4次谐波幅g984 g2097在各个调制比 m 下的开关角代g1049到谐波计g4778表达g2351(2)中即g1338g2409到消谐之g1360g1179g887的g3107低g719次谐波的含量，利用g2351(13)g1338以计g4778出不g1358 N g3024，在求g2409的g719组开

42、关角作用下输出电压的HDF，图4为HDFg6947m变g1232的g3099线 图4 两组解的HDF对比 Fig.4 Comparison of HDF between two sets of solutions g815图 4 中g719组开关角作用下输出电压的电压谐波的 HDF g3150g4409，第 2 组开关角作用下的特性g5941好g773第1组开关角，g2290g708g719组开关角作用下的特性差g6255g6947g4427开关角g2972目的g1816多而逐渐g1106小 3 改进混合粒子群算法分析 为了g7181g6012g3124g2976所提g4778法的性能，以 m

43、0.6 g3024求解g719个开关角为例，g2090传统g4833g2004g5061(CPSO)、g2564性g3134重线性变g1232g4833g2004g5061(LPSO)、传统g3380g2737g6429g3873g4833g2004g5061(SAPSO)g4778法与g2951g6408的g3720合g4833g2004g5061(HPSO)g4778法的求解g6398程g6408行了比较，各g4833g2004g5061g4778法参g2972如表 1 所示，比较g4934g3174如图5所示 g815图 5 g1338以g4409出，CPSO g6941g1049了g2

44、124g6537g3107g853解g7755LPSO g2949g2964g3024迭代次g2972较g2101，g2949g2964g3024适应g2318g984较小， 2016年5月 陈 炜等：基于改进混合粒子群优化的特定谐波消除脉宽调制策略 503 g6060g3033求解角g2318g4852g2318比 CPSO 高g7755SAPSO 求解的g4852g2318高g773 LPSO 和 CPSO，但是g1060g2949g2964g3024迭代次g2972较多g7755而g3124g2976g6611用的 HPSO 则在求解到g3107g853解的g1358g3024，保g60

45、12了g4778法以g3107g2101的迭代次g2972g2949g2964 表1 不同粒子群优化算法参数 Tab.1 Parameters of different PSO algorithms 算 法 kmaxs c1c2w CPSO 200 30 2 2 0.8 LPSO 200 30 2 2 wmin0.2 wmax0.9 SAPSO 200 30 c102.5 c110.5 c200.5 c212.5 wmin0.2 wmax0.9 HPSO 200 30 c102.5 c110.5 c200.5 c212.5 wmin0.2 wmax0.9 图5 不同粒子群算法求解过程对比 Fi

46、g.5 Comparison of solving process of different PSO algorithms 为了g6408g694步g7181g6012g4778法的g853越性，在g4400g1358g3149件下g2090牛顿迭代法和各g4631g4833g2004g5061g4778法g6408行性能g2090比分析以 m0.6 g3024求解g719个开关角为例，在g6947机初g984的情况下g2090迭代法和各g4631g4833g2004g5061g4778法均运行 20 次，g2090g2288均g1260次迭代耗g3024、g2949g2964g3024g22

47、88均迭代次g2972、g2949g2964g3024g2288均g2483体迭代耗g3024以及g2288均适应g2318g984g6408行g2090比分析，g4934g3174如表2所示 表2 不同算法性能对比 Tab.2 Comparison of performance of different algorithms 算 法 初值 范围 单次耗 时/s 迭代 次数 总体耗 时/ms 适应 度值 迭代法 (0，/2) 不收敛 不收敛 不收敛 不收敛CPSO (0，/2) 096.5 49 4.73 110-3LPSO (0，/2) 098.5 39 3.84 110-4SAPSO (0，/2) 104.1 72 7.50 110-7HPSO (0，/2) 108.0 30 3.24 110-9g815表 2 g1338以g4409出，g6611用牛顿迭代法g3024，赋与g4833g2004g5061g4778法g1358g3239(0，/2)之间的g6947机初g984，g4778法不g2949g2964，而g4833