216 2011,47(25) ComputerEngineeringandApplications计算机工程与应用 @工程与应用@ 新型的PID控制器参数整定方法 徐志成 ,王树青 XU Zhicheng , NG Shuqing 1.常州机电职业技术学院,江苏常州213164 2,浙江大学工业控制技术国家重点实验室,先进控制研究所,杭州310027 1.Changzhou Institute of Mechatronic Technology,Changzhou,Jiangsu 213164,China 2.National Lab of ndusItrial Control Technology,Institute of Advanced Process Control,Zhejiang University,Hangzhou 310027,China XU Zhicheng,WANG Shuqing.Novel PID controller parameters tuning method.Computer Engineering and Applications, 2011,47(25):216—219. Abstract:There exists the disadvantages such as premamrity in particle swarm optimization because of the decrease of swarm diversity.n order to solIve this problem,a new and improved version of particle swarm optimization algorithm is pro- posed combining the global best and local best mode/,termed GLBest-PSO algorithm.This algorithm incorporates global-local best inertia weight with global-local best acceleration coeficifent.The velocity equation of the GLBest-PSO algorithm is sim- pliied afnd the performance of the algorithm is improved.The abiliy of tthe GLBest-PSO algorithm is tested with a set of bench mark problems and the results of the simulation show the validity and better optimization performance.The algorithm is proposed to design the parameter optimization of PID controller.The simulation results show that the optimal PID control- ler based on the proposed method has a satisfying performance and is superior to the conventional PID controller based on the conventional tuning methods. Key words:particle swarm optimization;global-local:PID controller;parameters tuning ‘ 摘一要:粒子群优化算法在进化中随种群多样性降低易出现早熟收敛等问题。针对这一问题,结合全局.局部最优模型,提出了 种改进的粒子群优化算法,称为全局.局部参数最优的粒子群优化算法。算法利用全局一局部最优惯性权重及全局一局部最优加 速度常数,算法的速度更新方程被简化,性能得到改善。利用一组bench mark问题对该算法进行测试,仿真结果表明了算法的有 效性和高效性。将该算法应用到对传统PID控制器的参数优化当中,仿真结果表明方法可以获得满意的控制效果,各项控制性 能指标优于传统方法整定得到的PID控制器。 关键词:粒子群优化;全局.局部;PID控制器;参数整定 DOI:10.3778 ̄.issn.1002.8331.2011.25.057 文章编号:1002.8331(2011)25.0216.04 文献标识码:A 中图分类号:TP18 l 引言 粒子群优化算法(Particle Swarm Optimization,PSO)“ 本原则Ⅲ,而“粒子”则是一个折衷的选择,它一方面将群体中 的成员描述为没有质量和体积的个体,另一方面赋予其速度 和加速度。 由Kennedy和Eberhart等人于1995年提出,鉴于其发展历史较 短,理论基础及应用推广仍存在着一些问题,本文拟针对提高 基本PSO算法收敛速度及解决易陷入局部最优解的特点,提 出了一种改进的PSO算法,即全局-局部参数最优粒子群优化 算法(Globa1.Local Best PSO,GLBest.PSO),以期对基本PSO 算法做出改进。 PSO算法的原理为 :D维空间中存在m个粒子,每个粒 子坐标为x =( , ,…, ),同时每个粒子具有速度vi=(vfl’ v …,v 。该粒子的历史最好位置P = f1’P …,p ,记 P 。 ;记群 中所 代的第J维速度和位置。 t+1pijvg() =vij(O+r ̄q(一 Pg=(gf1,ga,…,go), 记g 。对第t代的第f个粒子,PsO算法根据下式计算第f+1 )+r2 2 一 ) (1) 2基于全局.局部参数最优的粒子群优化算法 2.1标准PSO算法 PSO算法中,“群(swarm)”来源于Millonas在开发应用于 人工生命(artiifcial life)的模型时所提出的群体智能的5个基 ,十D: (f)+V (f+1) .,、,、式中, 、,:为【0,1】的随机数;c。和c:为加速度常数(acceler- 基金项目:国家自然科学基金(the National Natural Science Foundation of China under Grant No.60421002)。 作者简介:徐志成(198O一),硕士研究生,主要研究方向:智能优化和鲁棒控制;王树W(1939--),教授,博士生导师。E—mail:zexu_ezmec@163.eom 收稿日期:2010.05—14;修回日期:2010—09—06 徐志成,王树青:新型的PID控制器参数整定方法 ation constants),此外,粒子速度v 被最大速度 w口 ,所限制。 文献【2】研究表明:c =C =2.0时算法的寻优性能更佳,因 此,大多数学者默认情况下一般取C =c,=2.0,但一些研究表 明 :表征粒子“认知”能力的C 值大于表征粒子“社会”能力的 后来,Shi Y等人又增设了惯性权重(inertia weight)因子 将式(1)变为: (,+1) wry(t)+rjcj 一 (t))+r2c2(glj一 ( ) f,, 、 C,时,效果更加,但必须满足表达式(5): c1+C2 4 (5) (f+1)= (D+vi, +1) 、一 式(2)通常被认为是基本PSO算法。 Clerc M引入了限定因子 ,构成PSO算法的限定因子 法 ,即 V (H 1) X(Vo(f)+ l — (f))+ 2 一Xi)O))) 基于这种思想,文献【5]对c.和C,的值提出了一种线性变 化加速常数法,即随着进化次数的增加,减少粒子的“认知”能 力,增加粒子的“社会”能力,c,和C,随着时间线性变化。在搜 … f,(f+1): , )+VqO+1) 南 4, 一般选择 为4.1,则 为0.729 8,这相当于W为定值0.729 8, c =C,=1.494 45的基本PSO算法。 2.2算法参数分析及改进 在基本PSO算法中包含3个权重因子:惯性权重W,加速 度常数C 和C:。惯性权重W使粒子保持运动惯性,使其具有 扩展搜索空间的趋势,有能力探索新的区域。引入陨性权重W 可消除对 的需要,因为它们的作用都是维护全局和局部 搜索能力的平衡。这样,当 增加时,可通过减小W来达到 平衡搜索。而w的减小可使得所需的进化次数变小。从这个 意义上看,可以将1,一 .d固定为只对W进行调节。式(2)中,如果令w=0,则速度只取决于粒子当前位置和 其历史最好位置Ph 和gh。 速度本身没有记忆性。假设一 个粒子位于全局最好位置,它将保持静止。而其他粒子则飞 向它本身最好位置P 和全局最好位置g 。 的加权中心。在 这种条件下,粒子群将收缩到当前的全局最好位置,此时,PSO 算法更像一个局部算法。 对全局搜索,通常好的方法是在前期具有较高的探索能 力以得到合适的种子,而在后期有较高的开发能力以加快收 敛速度。基于这种思想,Shi Y采用线性递减权重(Linearly Decreasing Weight,LDW)策略 ,即 w:(w1一w,)(—iter_m ̄—x—-iter—一)+w, (4) ‘Her㈨ ‘ 其中,W,和W,为惯性权重W的初值和终值;iter是算法进化 过程中的进化代数;f 一是最大进化代数,可将w的范围设 为由1.4到0,由0.9到O.4,由0.95到0.2等,目前大多数学者在 使用基本及改进的PSO算法时,对于惯性权重W的取值一般 在上述范围内。 加速度常数c。和C:表示将粒子推向p 位置和g 。 位 置的统计加速项的权重。低的值允许粒子在被拉回之前可以 在目标区域外徘徊,而高的值则导致粒子突然冲向或越过目 标区域。 式(2)中,如果令C,=O,则粒子没有认知能力,也就是“只 有社会”的模型。在粒子的相互作用下,有能力到达新的搜索 空间。它的收敛速度比标准版本更快,但对复杂问题,则比标 准版本更容易陷入局部最优点。 式(2)中,如果令C,=O,则粒子之间没有社会信息共享,也 就是“只有认知”的模型。因为个体问没有交互,一个规模为 的群体等价于m个粒子各自独立运行,因而得到解的概率 非常小。 索初期,c,的值较大,c,的值较小,以达到允许粒子在更大范 围内搜索,而代替原来单纯地向自身最优值的移动;随着进化 代数的增加,C,线性减小,而C,线性增加,使粒子向全局最优 解移动,基于该思想的PSO算法称为时变加速常数粒子群 (Time.Varying Acceleration Coefifcient PSO,TVAC.PSO)算 法。在该算法中,C.和C,的取值可以用下面两式表示。 C1:(c 一c r)(iter ̄_m_-iter.)+c。, (6) ’ max ,口 — ,口 c2:(c2f—c2r)( )+c2r (7) ~,max 其中,c 。和C, 是加速度常数的初值;C ,和C,,是加速度常数 的终值。大量的实验数据表明 :随着进化次数的增加,C,的 取值从2,5到0,5线性变化;c,的取值从O.5到2.5线性变化,优 化效果最佳。 研究表明:若PSO算法过早收敛,粒子速度将下降至0,粒 子群将趋于当前的极值P 。 和g ,而它们往往为局部极值, 尚未达到全局最优。因此对PSO算法的改进不能着眼于收敛 性,而应调节算法的搜索范围,以及全局和局部搜索能力。鉴 于贤洼权值对粒子速度的影响以及加速度常数C 和c,是决定 粒子“认知”和“社会”能力的关键参数,本文中,借鉴文献【2,5】 中的思想,惯性权重W,加速常数c.和C,既不取恒值,也不随 进化次数的增加而线性变化,而是表示成局部最优和全局最 优的适应度函数,可以用下面两个表达式来表示: (8) c=(1十 gb est,..J (9) 式中,W 为每个进化代数对应的惯性权重,gbest 为该进化代 数的全局最优位置;pbest.为该进化代数对应的所有粒子历 史最优位置的平均值。表达式(8)称为全局.局部平均最优惯 性权重(Globa1.Average Local Best IW,GLBest IW);表达 式(9)称为全局.局部最优加速度常数(Globa1.Local Best.AC, GLBest.AC),则相应的速度更新表达式为: vg(t+1)=WV (f)十crl(pbesti+gbestf一2xf,O)) (10) (,+1)= ,(,)+v O+1) (11) 表达式(8)~(11)所表示的PSO算法称为全局.局部最优 粒子群算法(Globa1.Local Best PSO,GLBest.PSO)。 由式(8)~(11)可以看出,当全局最优值等于局部最优值 时,全局.局部平均最优惯性权重的值达到最小,这实际上使粒 子在全局最优值附近搜索,并且迅速地向最优值收敛。同样, 当全局最优值等于局部最优值时,全局.局部最优加速度常数 等于2,并且在整个搜索过程中,其值始终位于2附近。该两个 参数帮助GLBest.PSO算法获得更佳的寻优性能。 Computer Engineering册 4 疗c口砌 计算机工程与应用 2-3算法性能测试 Rosenbroek函数是很难极小化的病态单峰二次函数,一般 形式的Rosenbrock函数表示为: n—l minf(x)=∑[100(x…一 2) + 一1) ] j=l ∈卜10,10】) 图2反馈系统结构示意图 Rosenbrock函数具有全局极小点 =(1,1,…,1),其函数 式中, 、 、 分别表示PID控制器的比例、积分、微分参 数。经过拉氏变换后,PID控制器传递函数可描述为: c( = = + + (13) 值为. =0,其全局极小点位于一个狭长、抛物线形状的极 为平坦的函数曲面山谷中,图1为2维该函数的几何特性。一 般来说,要想到收敛Rosenbrock函数全局极小点十分困难,尤 其是高维l青况下。 式中, = 佤, = 。图2中G 表示工业对象或过 程的传递函数。典型的工业过程通常可以近似为一阶纯滞后 ×104 二 4 图1 2维Rosenbrock函数的几何特性 利用上述3种算法对该函数取10维进行测试,种群粒子 数为30,LDW-PSO算法中,c,=C,=2,w从0.9到0.4随着进化 代数增加而线性减小;TVAC—PSO算法中,c,和C,分别按式(6) 和式(7)在[2.5 0.5]和[0.5 2.5】区间线性变化,W从0.9到0.4 随着进化次数的增加而线性减小。GLBest.PSO算法中,w和 C的初值在[0.1 1】和[1 2】中随机产生。3种算法的终止条件均 为进化代数达到1 000,对3种算法分别进行50次实验,表1给 出了3种算法所达到精度的平均值。 表l给定进化代数时各方法达到精度比较表 表l数据表明,在给定的进化代数下,本文算法所达到的 精度最高,搜索性能强缘于它的收敛性好,能迅速有效地靠近 函数唯一的最优值,与其相比,TVAC.PSO算法的优化性能次 之,LDW-PSO算法最差。 3基于GLBest.PSO算法的PID控制器参数优化 3.1问题描述 PID控制器由于结构简单、容易实现以及鲁棒性好等特 点,在目前工业控制回路中仍然占据着重要的地位。PID控制 器性能的好坏与对其参数的整定有很大关系,对PID控制器参 数的优化整定一直是控制领域研究的重要问题,随着一些智 能优化算法的提出,出现了利用遗传算法、粒子群优化算法等 对PID控制器参数进行优化整定的方法一。 ,并且取得了很好 的效果,在此研究基础上,现利用本文提出的GLBest.PSO算 法对PID控制器的参数优化进行尝试。 考虑如图2所示的反馈控制系统,图2中 )表示PID控 制器的传递函数,假设印)表示系统误差, (f)表示控制器输 出,则 “(f)=K,te ̄0+ (r)+ ] (12) 对象或者二阶滞后对象,其传递函数表示为: G = (14) 』 十l G :— 。 s +As+B (15) 这里f表示过程滞后时间, , , 为过程参数。 PID参娄嫩问题涉 U寻拢一{ 数 =L ,K, 】, 使得控制系统在快速响应的同时又满足超调量小、稳定时间 短等性能目标。 3.2实例仿真 飞锯nu是对连续生产的无限长工件进行定长切割的锯床, 由于具有高效率、高精度的特点,目前自动化生产线上飞锯的 使用比较广泛。整个飞锯系统由滑块刀架、位置伺服系统、自 动卸料机构以及可编程序控制器PLC控制的辅助设备组成。 这里采用PID控制器来实现对位置伺服结构控制,飞锯交流调 速系统的传递函数可近似为nu: G = e (16) 式中,对象参数为 =1,f=0.95, =1.72。 采用传统PID控制器对该系统进行控制,分别利用Z-N 法、LDW-PSO算法、TVAC.PSO算法和GLBest.PSO算法来优 化整定得到PID控制器的3个参数以及对应的目标优化值(算 法参数的设定与本文2.3节算法性能测试中所选参数相同,以 ISE值作为目标函数值,如表2所示),单位阶跃响应曲线如 图3所示。 表2 PID控制器参数整定结果 整定方法 尼 ISE值 Z-N 4.564 2 2.452 1 1.821 3 9.636 7 LDW-PSO 1.477 5 0.812 2 0.319 9 6.664 7 TV C.PSO 1.943 2 0.453 3 0.123 5 3.421 3 GLBest—PSO 0.245 7 0.412 1 0.153 8 2.654 5 …『I= 一Z-N 一,\ ^ r—— . 、 Dw_】 So 一 \ TV_AC.P ;o I \GI esh PSO 0 1O 2O 30 40 50 60 70 80 90 l0o time/s 图3四种PID参数整定控制方法比较图 徐志成,王树青:新型的PID控制器参数整定方法 2011,47(25) 219 可以看出,本文提出的GLBest.PSO算法整定得到的PID convergence in a multidimensional complex space[J].IEEE Traas・ 控制器具有快速、无超调的响应特性,相对于传统的PID控制 actions on Evolmionary Computation,2002,6(1):58—73. 器参数整定方法z.N算法、LDW-PSO算法、TVAC.PSO算法, [4】Shi Y,Eberhart R C.Fuzzy adaptive particle swarm optimiza- 性能指标均有所提高,优化目标函数值变小。 tion[C]//Proc of the IEEE Congress on Evolutionary Computa- tion.Seoul。Korea:IEEE Press,2001:101—106. 4结论 【5]Ratnaweera A,Halgamuge S K,Watson C.Self-organizing hierar- chical particle swarlTl optimizer with time—varying acceleration 作为群集智能的代表性方法之一,粒子群优化给大量非 coefifcients[J].IEEE Trans on Evolutionary Computation.IS.1.】: 线性不可微和多峰值复杂优化问题的求解提供了一种新的思 IEEE Press,2004,8(3):240—255. 路和解决方法。鉴于PSO算法的理论基础以及应用推广仍存 [6]Chen Debao,Zhao Chunxia.Particle swalTn optimization、vith 在着一些问题,本文拟针对提高基本PSO算法收敛速度及解 adaptive population size and its application[J].Applied SoR 决易陷入局部最优解的特点,提出了一种改进的粒子群优化 Computing,2009,9:39—48. GLBest.PSO算法,利用全局.局部最优惯性权重及最优加速度 [7】Anghinolfi D,Paolucci M.A new discrete particle swarln optimi— 常数,提高了算法的收敛速度和寻优性能。将方法应用于对 zation approach for the single—machine total weighted tardiness 传统PID控制器的参数优化中,仿真实例得到了令人满意的结 scheduling problem with sequence—dependent setup times[J].Euro・ pean Journal of Operational Research,2009,1 93:73-85. 果,表明了该方法的有效性和合理性。 [8]Arumugam M S,Rao M V C,Tan A W C.A novel and effec— tire particle swanTl optimization like algorithm witll extrapola— 参考文献: tion technique[J].Applied Soft Computing,2009,9:308—320. [1]Kennedy J,Eberhart R.Particle swarm optimization[C]//Proc IEEE [9]刘楠楠.改进的多目标遗传算法及其在PID优化设计中的应用[J】. Int Conf on Neural Networks.Perth:IEEE Press,1 995. 应用科学学报,2010,28(1):83—89. [2]Shi Yuhui,Eberhart R.A modiifed particle swarln optimizer[C]// [1O]吴文进.基于粒子群优化算法的P1D液位控制[J】.合肥工业大学学 Proc IEEE Int Conf on Evolutionary Computation.Anchorage: 报:自然科学版,2009,32(1):1674—1677. IEEE Press,1997:303—308. [11】李先银,胡乾斌,李光斌 最优PID控制算法在飞锯位置伺服系统 [3】Clerc M,Kennedy J.The particle swarm explosion,stability,and 中的应用fJ].电气传动,2001,5(1):41—43. (上接181页) Iar Gabor filters[C]//16th IEEE IntemationaI Conference on Pat- [5]焦峰,崔国勤,高文,等.基于局部特征分析的人脸识别方法[J]l计 tern Recognition,2002:901—904. 算机辅助设计与图形学学报,2003,15(1):53—58. [8]Porter R,Canagarajah N.Robust rotation-invarinat texture classiif- [6】Ayinde 0,Yang Yee—hong.Face recognition approach based on cati0n:wavelet,Gabor filter and GMRF based schemes[J].1EE Proceedings Vision,Image and Singal Processing,1997,144: rank correlation of Gabor-filtered images[J].Pattern Recogniiton, 18O一188. 2002,35(6):1275 1289. [9]Hyvafinen A,Karhunen J,Oja E.Independent component an ̄ysisBVq. 【7]Zhang J,Tan T,Ma L.Invariant texture segmentation via circu— New York:John Wiley,2001. (上接204页) 参考文献: 分类器与最近邻分类器的分类识别率相同(而在非最优参数 [1]李坚,刘一星.木材涂饰与视觉物理量【M].哈尔滨:东北林业大学 组合下, 近邻分类器的分类识别率高于最近邻分类器的分 出版社,1998:325.540. 类识别率),集成神经网络的总体分类识别率达到94.0%, 丘 [2]白雪冰,王克奇,王辉.基于灰度共生矩阵的木材纹理分类方法的 邻分类器的总体识另IJ/ ̄91.0%,获得了较高了分类识别率。 研究[J].哈尔滨工业大学学报,2005,37(12):1667—1670. [3】王克奇,石岭.基于高斯一马尔可夫随机场的木材表面纹理分析【J]. 5结论 林业科技,2005,30(6):4648. 5阶GMRF随机场参数对自然纹理有较强的描述能力,用 [4]熊四昌.基于计算机视觉的刀具磨损状态监测技术的研究[D】.杭 州:浙江大学,2003. 模拟退火算法除去了参数的冗余,选择最优参数。分别设计 [5】Woods J W.Two—dimentional discrete markovina fields[J].IEEE 了 近邻分类器、最近邻分类器和集成神经网络分类器,结果 Transactions on Information Theory,1972,18:232-240. 表明集成神经网络分类器识别率最高达到了94.0%,分类精度 [6】边肇祺,张学工.模式识别[M】.2版.北京:清华大学出版社,2000. 基本能够满足木材加工企业要求,从而提高了木材加工工业 [7]王克奇,王辉,白雪冰.基于模拟退火算法和最近邻分类器识别率 的自动化水平。 的特征选择方法[J].自动化技术与应用,2007,26(I):27—29.