基于单纯形法的PID参数寻优_pid参数的经验整定法
基于单纯形法的PID参数寻优由刀豆文库小编整理,希望给你工作、学习、生活带来方便,猜你可能喜欢“pid参数的经验整定法”。
基于单纯形法的PID参数寻优
1.引言
由于PID 控制具有直观、易于实现、稳态精度高以及鲁棒性和稳定性较好等优点,并且人们对其原理和物理意义等比较熟悉,已经建立了比较完善的理论体系,尤其在工业现场控制过程中应用尤为普遍,特别适用于对象动态特性未完全掌握、得不到精确地数学模型、难以用控制理论来分析和综合的场合。
常规的PID 参数,通常做法是相应于某一点求取相对最优的调节参数后, 应用于整个运行过程,是不太合理的。此外, 参数的整定十分麻烦、复杂, 需要进行被控对象特性和过渡过程特性的测试和计算, 或者需要借助于积累的调试经验, 才能获得比较满意的整定效果,为了减少PI D 调节器参数整定的麻烦, 克服因被控对象特性变化或扰动作用造成的系统性能的降低。本文介绍采用数字PI D 调节器的参数自寻最优控制。所谓自寻最优控制是利用计算机的快速运算和强大的逻辑判断能力, 按照选定的寻优方法, 不断探测、不断调整, 自动寻找最优的数字P ID 调节参数, 使得系统的性能处于最优状态。
参数的寻优方法有很多,如单纯形法、最速下降法、拟牛顿法等,都是对目标函数进行计算,来得到所需参数的目的,计算量很大.但是.由于单纯形加速法具有控制参数收敛快、计算工作量小、简单实用等特点.在线数字PID 参数自寻最优控制中比较普遍使用单纯形加速法。
2.数字PID 控制算法
PID控制器是一种线性控制器,它根据给定值与实际输出值的偏差构成控制偏差,将偏差的比例P、积分I和微分D通过线性组合构成控制量,从而对被控对象进行控制。PID控制器的控制规律为:
DsUsEskP(11TDs)(2.1)TIs其中:kP是比例系数;TI是积分时间常数;TD是微分时间常数。采用后项差分变换法,将上式离散化,得到
DzKpT1KpTDU(z)1kp(1z)1E(z)TI1zT11K(1z)D11z(2.2)KpKI
其中,KlKpT/Tl为积分系数;KDKpTD/T为微分系数。比例控制的作用是通过加大比例系数来增加系统的动态响应速度;积分控制的作用是消除系统的稳态误差;微分控制作用是改善系统的动态性能。
由式(2.2)有:
uk[kpKl1KD(1z1)]e(k)(2.3)11z由上式计算出来的u(k)是控制量的绝对大小,如果控制量是阀门,则它反映了阀门的开度,因此
称为全位置式PID控制器。由式(2.3),有
uk-1[kpKl由式(2.3)与式(2.4),1KD(1z1)]z1e(k)(2.4)11z112ukk1zkk12zzPIDek(2.5)
进而得到
uku(k1)[kp(1z1)KIKD(12z1z2)]e(k)(2.6)式(2.5)计算得到的是控制量的增量,因此称为增量式或者速度式PID控制器,式(2.6)称为递推位置式PID控制器。
位置式PID和增量式PID本质相同,只是形式不同而已,对系统的控制作用完全相同,但是增量式PID要求执行机构具有记忆功能。实际应用过程中常使用增量式PID或者递推位置式PID,因为全位置式PID中含有积分项,需要存储过去全部偏差,计算编程不方便。
3.对象离散化数学模型
连续传递函数与离散传递函数的转换, 根据应用场合的不同, 通常有双线性变换法, 零极点匹配法、冲激响应不变法和零阶保持法等几种不同的转换方法。在计算机控制系统中, 离散信号uk后面通常接有保持器, 而最常用的是零阶保持器、对于调节对象前面含有零阶保持器的连续传递函数的离散化, 采用零阶保持法可以获得准确的离散等效。故本文采用零阶保持法求取调节对象连续传递函数的等效离散传递函数:
GsGz1z1Z(3.1)s本系统离散化后得到的传递函数为:
Y(z)0.048z10.472z2Gz(3.2)12U(z)11.9048z0.9048z得到uk与ek的关系,如下:
ukuk1b0ekb1ek1b2ek2(3.3)式中b0为kpkIkD,b1为-kP-2kDb2为kD。
得到yk与uk的关系,如下:
yk1.9048yk10.9048yk20.048uk10.0472uk2(3.4)
4.单纯形算法介绍和PID 参数的单纯形寻优流程
4.1单纯形算法介绍
在n维空间中的单纯形是一种多胞形,它具有n+1个不在同一超平面的顶点.若各个棱长彼此相等,则称为正规单纯形.所谓单纯形就是一定的空间中的最简单图形。N维的单纯.就是N+1个顶点组成的图形,如二维空间,单纯形是三角形。
设二元函JX1,X2构成二维空间,有不在一条直线上的三个点,XH,XGXL构成了一个单纯形。由三个顶点计算出相应的函数值JH,JG,JL。若,对于求极小值问题来说,JH最差,JG次之,JL最好。可以想象函数的变化趋势:一般情况下.JG次之,JL最好。可以想象函数的变化趋势:一般情况下.好点在差点对称位置的可能性比较大,因此将XG,XL的中点XF与XH连接.并在XFXG的射线方向上取XR,使XHXFXFXR.见图表一.以XR作为计算点,计算其函数JR。(1)若JRJG
说明步长太大,以致XR并不比XH好多少.因此需要压缩步长,可在XR与XH问另选新点XS。(2)若JRJG
说明情况有好转,而且还可以加大步长。可以XHXR的延长线上取一新点XE。若JEJR.取XE作为新点XS。若JEJR.取XR作为新点XS。总之.是可以得到一个新点XS。若JSJG,说明情况确有改善.可舍弃原求的XH点,而以XG,XL,XS三点构成一新的单纯形(XGXLXS),称作单纯形扩张,然后,重复上述步骤。
若JSJG.说明XS代替XH改善不大.可把原来的单纯形(XGXLXS),按照一定的比例缩小,例如边长都缩小一半,构成新的单纯形(XFXLXM).称作单纯形收缩。然后重复以前的步骤,直到满足给定的收敛条件。
图4.1 单纯形法示意图
4.2 PID 参数的单纯形寻优
采用单纯形法时,计算机必须根据一定的性能指标来确定控制性能的好坏。这个性能指标则通过目标函数来体现。本题目采用二次性能指标为
J=[e2(k)u2(k)]k=0(4.1)PID 的参数整定需同时确定3 个参数单纯形法不用梯度这个信息,而是利用对许多参数点的目标函数值的比较来确定寻优方向,适合于变量不太多的场合。从它们的大小关系看出函数变化趋势,为函数的下降方向提供参考。在搜索过程中主要采用延伸、收缩、压缩和扬弃的方法对各点进行处理。
单纯形寻优的初始化步骤如下:
1)输入初始点X0,单纯形边长l,反射系数r,延伸因子e,收缩因子c,收敛误差和最大搜索次数Maxl.
2)令k=0。按XiXihEi。方法构建4个点以形成初始单纯形.
具体流程图如图4.2。
5.实验仿真
5.1仿真模型
图5.1系统模型
5.2结果分析
(1)变化对参数影响
0.5,KP1.5003,KI2.7066,KP1.2067 0,KP3.9503,KI0.4672,KP0
可见,可参数随着的增大,参数会有变化,比例度变大,积分时间变小,微分时间常数为0,即可满足要求。
原因主要强调控制效果,较少计较控制代价,即主要考虑控制精度。比例度大一些,可提高整个系统的灵敏度,也可以相应减小余差。积分时间越小,控制作用越强,有助消除余差。微分控制作用的特点是:动作迅速,具有超前调节功能,可有效改善被控对象有较大时间滞后的控制品质;但是它不能消除余差。
因此,应该根据实际情况,合理选择加权系数。既要考虑控制效果,又要考虑控制代价。(2)控制效果分析
选用加权系数为0.5分析控制效果。
图5.1 0.5时,控制器输出
图5.20.5,系统输出图 图5.1表明控制的输出经过一段时间的波动后,趋于常值,因此控制器输出序列是稳定的。由图5.2可知,加入PID控制器之后系统输出超调量小,过渡时间较短,稳态误差小,得到了较满意的效果。
由以上可知,本文设计的控制器满足了闭环系统稳定性和控制器输出序列稳定性的要求。经过单纯形法参数寻优的PID控制器设计是令人满意的.
6.结论
由以上可知,本文设计的控制器满足了闭环系统稳定性和控制器输出序列稳定性的要求。经过单纯形法参数寻优的PID控制器设计是令人满意的.
但是还存在一些不足:
(1)计算延时时间过长,不能满足实时性。(2)容易陷入局部最优。
(3)对于PID参数给定的是经验值,实际对于不同的问题,这需要做大量的试验。(4)本文采用性能指标为未直接考虑调节时间,所以快速性不好。(5)未对扰动影响进行分析。
因此,优化算法有待进一步的研究和改进。