一种基于粒子群算法pid分离的参数整定方法
【专利摘要】本发明公开了一种基于粒子群算法PID分离的参数整定方法，包括以下步骤：初始化设定粒子群的粒子个数，惯性系数，加速系数，粒子最大迭代次数，初始速度，对应每个粒子的空间位置分量和PID初始值；计算每个粒子对应的目标函数值；根据目标函数值更新个体最优的粒子位置、全局最优粒子位置、每个粒子的位置分量，再根据PID映射公式，更新对应每个粒子此次迭代的PID的值；重复迭代，直到迭代结束；输出全局最优PID参数值和对应的PID值。本发明在常规的粒子群优化算法中加入了PID映射的算法，可以将粒子空间位置分量映射到三维的PID参数空间中，这样便于算法筛选出最优的粒子位置，提高PID参数整定的效率。
【专利说明】
一种基于粒子群算法PID分离的参数整定方法
技术领域
[0001]本发明涉及一种交流电机的PID参数整定方法，特别涉及一种基于粒子群算法PID 分离的参数整定方法。
【背景技术】
[0002] PID(比例（proportion)、积分（integral)、导数(derivative))控制器作为最早实 用化的控制器已有近百年历史。PID控制器简单易懂，使用中不需精确的系统模型等先决条 件，因而成为应用最为广泛的控制器，广泛应用于工业控制，比如电机速度控制系统等。但 传统的PID值通常使用手动调节，耗费较多的时间且往往不容易得到最优的PID参数。因此 国内外研究人员对PID的整定放大进行了大量的研究，但交流电机的速度控制器PID参数整 定仍然是一个待完善和待解决的问题。
【发明内容】
[0003] 为了解决上述技术问题，本发明提供一种方便、快捷的基于粒子群算法PID分离的 参数整定方法。
[0004] 本发明解决上述问题的技术方案是：一种基于粒子群算法PID分离的参数整定方 法，包括以下步骤：
[0005] 步骤一:初始化设定粒子群的粒子个数，惯性系数，加速系数，粒子最大迭代次数， 初始速度，对应每个粒子的空间分量和PID初始值；
[0006] 步骤二:计算每个粒子在交流电机中对应的目标函数值；
[0007] 步骤三:根据目标函数值，更新个体最优的粒子位置，更新全局最优粒子位置;更 新每个粒子的位置分量，再根据PID映射公式，更新对应每个粒子此次迭代的PID的值；
[0008] 步骤四：判断当前迭代次数是否和最大迭代次数相等，若是，则跳出迭代，执行步 骤五;若否，则当前迭代次数加1，跳到步骤二；
[0009] 步骤五:输出全局最优PID参数值和对应的PID值。
[0010] 上述基于粒子群算法PID分离的参数整定方法，所述步骤二中，目标函数为：
[0011] f = /t I S-Sref I dt
[0012]其中，S为电机实测速度，Sref为预期速度，t为时间。
[0013] 上述基于粒子群算法PID分离的参数整定方法，所述步骤三中，更新个体最优的粒 子位置具体步骤为:将此次迭代个体最优值与历史个体最优值进行比较，如果此次迭代个 体最优值大于或者等于历史个体最优值，则不更新个体最优的粒子位置;如果此次迭代个 体最优值小于历史个体最优值，则将此次迭代个体最优值作为当前个体最优值，更新个体 最优值。
[0014] 上述基于粒子群算法PID分离的参数整定方法，所述步骤三中，更新全局最优粒子 位置具体步骤为:判断找到的此次迭代的全局最优粒子位置是否小于历史全局最优粒子位 置;如果是，则将此次迭代的全局最优粒子位置作为当前全局最优粒子位置;否，则不更新， 继续沿用历史全局最优粒子位置。
[0015] 上述基于粒子群算法PID分离的参数整定方法，所述步骤三中，每个粒子的位置分 量更新公式为：
[0016] r*+1 = inf + α{}\(Phestf - xf ) + a2r2(Gbest - x*)
[0017] 4a=4+vhi
[0018] 其中，为第k个粒子第i+1次迭代的速度，:为第k个粒子第i次迭代的速度;xf为 第k个粒子第i次迭代的位置分量，为第k个粒子第i+Ι次迭代的位置分量;〇和^为(0，1) 之间的随机数，为第k个粒子第i次迭代时找到的个体最优值，Gbest为当前找到的全 局最优值。
[0019] 上述基于粒子群算法PID分离的参数整定方法，所述步骤三中，PID映射公式为：
[0023]其中，靖丨，心和球!为对应粒子xL的PID参数，即为第i + Ι次迭代的P，I和D; ，/f和￡)f分别为对应于第K个粒子第i次迭代的P，I和D参数;τΡ，τ^Ρτ?)分别为P，I和D分 离系数;λΡ，λ^Ρλυ分别为映射的方向系数，具体公式如下：
[0027] 其中，；if，/f和Z)f为对应第k个粒子在第i次迭代时的Ρ，Ι和D参数。
[0028] 本发明的有益效果在于:本发明在常规的粒子群优化算法中加入了 PID映射的算 法，可以将粒子空间位置分量映射到三维的PID参数空间中。其中，方向系数λ决定着P，I和D 映射的方向，保证Ρ，Ι和D的值都为正值;分离系数τΡ，τ4Ρτ?)决定着Ρ，Ι和D映射的步长，保证 PID系数在解空间里快速地增加。在本方法中，方向系数和分离系数可以更好地对应空间粒 子的PID参数进行控制，这样便于算法筛选出最优的粒子位置，可以快速方便地整定出交流 电机速度控制器的PID的值，节省了PID整定的时间，提高PID参数整定的效率，保证系统的 正常运行。
【附图说明】
[0029] 图1为本发明的流程图。
[0030] 图2为本发明的原理图。
【具体实施方式】
[0031]下面结合附图和实施例对本发明作进一步的说明。
[0032]如图1所示，本发明步骤如下：
[0033] 步骤一:初始化设定粒子群的粒子个数η，惯性系数w，加速系数ai和a2,粒子最大迭 代次数m，初始速度v，每个粒子的空间分量位置和与每个粒子相对应的PID初始值；
[0034] 步骤二:将粒子群中每个粒子对应的PID值代入交流电机的速度控制器中，如图2 所示。计算并得到对应的目标函数值f，f = /t I S-Srrf I dt,其中，S为电机实测速度，Srrf为预 期速度，t为时间。
[0035] 步骤三:根据目标函数值，更新个体最优的粒子位置，更新全局最优粒子位置;更 新每个粒子的位置分量，再根据PID映射公式，更新对应每个粒子此次迭代的PID的值。具体 步骤为：
[0036] 3-1)根据目标函数，找到目标函数值里的最小值，即为此次迭代个体最优的粒子 位置，再根据找到的所有粒子的此次迭代的个体最优值的最小值，即为此次迭代的全局最 优粒子位置。
[0037] 3-2)更新个体最优的粒子位置，具体步骤为:将此次迭代个体最优值与历史个体 最优值进行比较，如果此次迭代个体最优值大于或者等于历史个体最优值，则不更新个体 最优的粒子位置;如果此次迭代个体最优值小于历史个体最优值，则将此次迭代个体最优 值作为当前个体最优值，更新个体最优值。
[0038] 3-3)更新全局最优粒子位置，具体步骤为:判断步骤3-1)中找到的此次迭代的全 局最优粒子位置是否小于历史全局最优粒子位置;如果是，则将此次迭代的全局最优粒子 位置作为当前全局最优粒子位置;否，则不更新，继续沿用历史全局最优粒子位置。
[0039] 3-4)更新每个粒子的速度与分量位置，每个粒子的位置分量更新公式为：
[0040] = irvf +〇[/[( Pbest^ - x*) + a2r2 (Gbest - )
[0041 ] 4：n=xf+vf+l
[0042]其中，ν^Η为第k个粒子第i + 1次迭代的速度，vf.为第k个粒子第i次迭代的速度；^ 为第k个粒子第i次迭代的位置分量，为第k个粒子第i+1次迭代的位置分量;〇和^为(0， 1)之间的随机数，为第k个粒子第i次迭代时找到的个体最优值，Gbest为当前找到的 全局最优值。
[0043] 3-5)根据PID映射公式，更新对应每个粒子此次迭代的PID的值，PID映射公式为：
[0047]其中，硪i，祀i和At为对应粒子的PID参数，即为第i + 1次迭代的P，I和D; if，<和分别为对应于第K个粒子第i次迭代的P，I和D参数;τΡ，τ4Ρτ?)分别为P，I和D分 离系数;λ Ρ，λ^Ρλυ分别为映射的方向系数，具体公式如下：
[0051 ] 其中，；/f , /f和D/为对应第k个粒子在第i次迭代时的Ρ，I和D参数。此处的方向 系数λ起着决定映射的PID随着优化迭代方向的作用。分离系数τ起着映射的PID参数随着优 化迭代的变化步长的大小。
[0052]步骤四：判断当前迭代次数是否和最大迭代次数m相等，若是，则跳出迭代，执行步 骤五;若否，则当前迭代次数加1，跳到步骤二。
[0053]步骤五:输出全局最优PID参数值和对应的PID值。
【主权项】
1. 一种基于粒子群算法PID分离的参数整定方法，包括w下步骤： 步骤一:初始化设定粒子群的粒子个数，惯性系数，加速系数，粒子最大迭代次数，初始 速度，对应每个粒子的空间分量和PID初始值； 步骤二:计算每个粒子在交流电机中对应的目标函数值； 步骤Ξ:根据目标函数值，更新个体最优的粒子位置，更新全局最优粒子位置;更新每 个粒子的位置分量，再根据PID映射公式，更新对应每个粒子此次迭代的PID的值； 步骤四：判断当前迭代次数是否和最大迭代次数相等，若是，则跳出迭代，执行步骤五； 若否，则当前迭代次数加1，跳到步骤二； 步骤五:输出全局最优PID参数值和对应的PID值。2. 根据权利要求1所述的基于粒子群算法PID分离的参数整定方法，其特征在于:所述 步骤二中，目标函数为： f = /t I S-Sref I dt 其中，S为电机实测速度，Sref为预期速度，t为时间。3. 根据权利要求1所述的基于粒子群算法PID分离的参数整定方法，其特征在于:所述 步骤Ξ中，更新个体最优的粒子位置具体步骤为:将此次迭代个体最优值与历史个体最优 值进行比较，如果此次迭代个体最优值大于或者等于历史个体最优值，则不更新个体最优 的粒子位置;如果此次迭代个体最优值小于历史个体最优值，则将此次迭代个体最优值作 为当前个体最优值，更新个体最优值。4. 根据权利要求3所述的基于粒子群算法PID分离的参数整定方法，其特征在于:所述 步骤Ξ中，更新全局最优粒子位置具体步骤为:判断找到此次迭代的全局最优粒子位置是 否小于历史全局最优粒子位置;如果是，则将此次迭代的全局最优粒子位置作为当前全局 最优粒子位置;否，则不更新，继续沿用历史全局最优粒子位置。5. 根据权利要求4所述的基于粒子群算法PID分离的参数整定方法，其特征在于:所述 步骤Ξ中，每个粒子的位置分量更新公式为：其中，戊1为第k个粒子第i+1次迭代的速度，vf为第k个粒子第i次迭代的速度;坏为第k个 粒子第i次迭代的位置分量，為1为第k个粒子第i + 1次迭代的位置分量;ri和η为(0，1)之间 的随机数，掛故《为第k个粒子第i次迭代时找到的个体最优值，Gbest为当前找到的全局最 优值。6. 根据权利要求5所述的基于粒子群算法PID分离的参数整定方法，其特征在于:所述 步骤Ξ中，PID映射公式为：其中，韓，，複1和辕1:为对应粒子记1的PID参数，即为第i+1次迭代的P，I和D;Pik，婷和 蜡分别为对应于第K个粒子第i次迭代的P，I和D参数;τρ，τι和TD分别为P，I和D分离系数;λρ， λι和λ〇分别为映射的方向系数，具体公式如下：其中，；Pik，皆和Df为对应第k个粒子在第i次迭代时的Ρ，巧阳参数。
【文档编号】G05D13/62GK106094910SQ201610704401
【公开日】2016年11月9日
【申请日】2016年8月22日
【发明人】于文新, 王俊年, 隋永波, 王振恒
【申请人】湖南科技大学