CN111914981B - 基于粒子群-蚁群并行交叉算法的改进pi模型辨识方法 - Google Patents
基于粒子群-蚁群并行交叉算法的改进pi模型辨识方法 Download PDFInfo
- Publication number
- CN111914981B CN111914981B CN202010481274.2A CN202010481274A CN111914981B CN 111914981 B CN111914981 B CN 111914981B CN 202010481274 A CN202010481274 A CN 202010481274A CN 111914981 B CN111914981 B CN 111914981B
- Authority
- CN
- China
- Prior art keywords
- model
- algorithm
- improved
- ant
- node
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000002245 particle Substances 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005457 optimization Methods 0.000 title claims description 11
- 239000003016 pheromone Substances 0.000 claims description 22
- 230000006870 function Effects 0.000 claims description 20
- 238000013507 mapping Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 7
- 239000013598 vector Substances 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 241000257303 Hymenoptera Species 0.000 claims description 5
- 230000004927 fusion Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 230000008020 evaporation Effects 0.000 claims description 3
- 238000001704 evaporation Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 4
- 230000006872 improvement Effects 0.000 abstract description 2
- 230000001737 promoting effect Effects 0.000 abstract description 2
- 238000006073 displacement reaction Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 230000001174 ascending effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000002431 foraging effect Effects 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
- G05B11/42—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Software Systems (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Computational Mathematics (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Automation & Control Theory (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Feedback Control In General (AREA)
Abstract
一种基于粒子群‑蚁群并行交叉算法的改进PI模型辨识方法,属于系统辨识技术领域。本发明的目的是利用具有变化斜率的改进Play算子和死区算子实现了对传统PI模型的改进,使建立的模型具有描述非对称迟滞特性能力的基于粒子群‑蚁群并行交叉算法的改进PI模型辨识方法。本发明步骤是:得到改进的PI模型,设计粒子群‑蚁群并行交叉算法辨识改进PI模型的参数,搭建压电微定位平台采集辨识模型所需的数据,依据步骤二所述的粒子群‑蚁群并行交叉算法辨识最终的模型参数并给出建模结果。本发明对于推进研究消除压电微定位平台迟滞非线性的方法,推广压电微定位平台使用,具有很大的研究意义。
Description
技术领域
本发明属于系统辨识技术领域。
背景技术
在现代高精密仪器及超精密加工技术中压电微定位平台是核心驱动部件,具有体积小、质量轻、分辨率高和不受环境温湿度变化影响等优点,但其输入输出映射表现的固有非对称迟滞非线性严重制约了压电微定位平台的更广泛应用。因此建立压电微定位平台精确模型,有利于加快开展降低其迟滞特定影响的方法研究,目前用于描述迟滞特性的主要模型按照建模原理区分主要有基于材料物理机理建立的物理机理模型还有基于材料表征的迟滞特性现象所建立的唯象类模型。由于物理类模型依靠对特定材料内部机电关系的熟知和严格的推导,故该模型建立极其繁琐且很难推广;唯象模型是指根据迟滞现象所建立,可以不考虑材料内部物理特性,使得建模难度大大降低。其中PI模型是唯象类模型,由于传统的PI模型所描述的迟滞曲线为上升与下降部分关于算子中心对称且函数表达形式为两边凸起的环状,而通过实验所测得的压电微定位平台的迟滞环非中心对称,这就造成传统PI模型表征压电微定位平台能力不足,从而导致模型精度不高。
发明内容
本发明的目的是利用具有变化斜率的改进Play算子和死区算子实现了对传统PI模型的改进,使建立的模型具有描述非对称迟滞特性能力的基于粒子群-蚁群并行交叉算法的改进PI模型辨识方法。
本发明步骤是:
步骤一:根据传统PI模型和死区算子得到改进的PI模型:先依据传统单边play算子提出映射关系为曲线形式的改进单边play算子,然后串联死区算子得到非对称的改进单边play算子,进而加权叠加得到改进PI模型;
传统PI模型的离散形式为Play算子的加权和形式:
式中,k=1,2,...,K,K为系统输入数据个数,n=1,2,...N,N代表PI模型中Play算子总数目,wn表示第n个Play算子的权值系数,Fn[u](k)为第n个单边Play算子,传统单边Play算子形式为:
式中:
rn表示第n个算子的阈值,0≤r1<r2<…<rN≤rmax,为了实现较少算子情况下精确表述迟滞曲线,将传统PI模型进行改进;
首先提出改进Play算子,其映射关系为曲线形式,表达式如下:
其中wb(k)表示算子在变化过程中的斜率,wb(k)=Fn′[u](k-1)e-Δu(k),Δu(k)=u(k)-u(k-1);然后采用单边死区算子与改进Play算子串联结合,得到基于死区算子的非对称改进PI模型,其中单边死区算子形式如下:
最终得到改进的PI模型表达式如下:
步骤二:设计粒子群-蚁群并行交叉算法辨识改进PI模型的参数:首先分别给出粒子群与蚁群算法辨识模型参数过程,为了综合发挥两种算法的优点随后提出粒子群-蚁群并行交叉算法用于辨识模型参数;
2.1粒子群优化
用一个矩阵表示,个体的运动速度表示为选出其中距离目标最近的个体坐标作为最优的个体坐标,将代表当前种群中每个粒子最佳位置,gbest=[p1,p2,…,pD]是全局最佳位置,c1代表学习因子,c2代表社会系数,个体的坐标和速度更新限定在和内取随机值,G为最大迭代次数,粒子群优化算法步骤如下:
(1)粒子的初始化:
初始化每个粒子的位置和速度:
(2)计算每个粒子的适应度函数,并根据适应度函数选取第g次迭代的每个粒子最佳位置和全局最佳位置:
选取适应度函数:
gbest(g)=min(pbestm(g)) (10)
(3)每个粒子的位置和速度的更新公式为:
Vm(g+1)=Vm(g)+c1rand(pbestm(g)-Xm(g))+c2rand(gbest(g)-Xm(g)) (11)
Xm(g+1)=Xm(g)+Vm(g+1) (12)
(4)根据适应度函数选取第g次迭代的每个粒子最佳位置和全局最佳位置,如果迭代次数达到设定的值,此时将gbest为最优解;否则,令迭代次数加一返回步骤(2);
2.2蚁群算法
(1)参数初始化
蚁群算法在用于模型参数辨识时,首先确定待辨识参数个数为I个,那么蚁群需要经过I条相互平行且等距离间隔的线段上的节点,每条线段均等分为(J-1)分,即每条线段各有J个节点,节点值记为j=1,2,...,J,设定蚂蚁总数目为L,每只蚂蚁l经过的节点定义为一个矩阵其中h=1,2,...,H迭代次数,H为最大迭代次数,其中待辨识参数的初始取值由如下公式给出:
(2)计算第l只蚂蚁的适应度函数,
最佳适应度
flmin(h)=min(fl(h)),l=1,2,...,L(15)
选取适应度函数值最小的蚂蚁经过的节点为当前迭代周期的最佳路径矩阵Zgood(h)=[z1(h),z2(h),...,zI(h)];
(3)计算第l只蚂蚁由第i-1条线段上任意节点向下一条平行线上节点j移动的概率
第l只蚂蚁由第i-1条线段上任意节点向下一条线段上概率最大的节点移动;
(4)更新蚂蚁l经过的路径节点
ji为第h迭代中,第l只蚂蚁由第i-1条线段上任意节点向下一条线段上概率最大的节点值;
(5)更新信息素
表示第l只蚂蚁由第i-1条线段上任意节点向下一条线段上移动所路径上释放的信息素浓度增量,Δτ(zi(h),j)表示所有蚂蚁在第i-1条线段上所有节点向下一条线段上节点j移动连接路径上释放的信息素浓度之和;
(6)判断是否达到最大迭代次数H,若否,则当前迭代次数加一,返回步骤(2),若是,输出最优解;
2.3粒子群-蚁群并行交叉算法
(1)产生初始种群:在初始范围[a,b]内划分N个小区间,每个小区间随机产生两个数,记Q1,Q2,...,QN作为粒子群算法粒子初始位置,P1,P2,...,PN作为蚁群初始路径节点对应数值信息;
(2)依据粒子群算法与蚁群算法的执行步骤并行执行两种算法,执行10代以后,利用公式(8)和(15)分别计算两种算法得到的新种群对应的适应度值;
(3)对两种算法得到粒子最佳位置和蚁群最佳路径节点对应的数值,将此时粒子群算法的位置数据计为S1,S2,...,SN,蚁群算法的路径节点数据记为R1,R2,...,RN;
(4)将R1,R2,…,Rξ与SN-ξ+1,SN-ξ+2,…,SN进行数据融合交换,令S1,S2,..,Sξ,R1,R2,...,Rξ作为粒子群的新位置数据,SN-ξ+1,SN-ξ+2,...,SN,Rξ+1,Rξ+2,...,RN作为蚁群的新路径节点数据,ξ取1到N/2之间的随机整数;
(5)判断是否满足算法终止条件算法最大迭代次数70,若满足,则算法终止;否则,粒子群-蚁群并行交叉算法的初始迭代周期T加一,转到步骤(2)继续执行;
步骤三:搭建压电微定位平台采集辨识模型所需的数据,依据步骤二所述的粒子群-蚁群并行交叉算法辨识最终的模型参数并给出建模结果。
本发明开展对传统PI模型进行改进,使其具有描述非对称迟滞环的能力,并设计系统辨识方案得到能够精确表征压电微定位平台输入输出映射关系的模型,对于推进研究消除压电微定位平台迟滞非线性的方法,推广压电微定位平台使用,具有很大的研究意义。
附图说明
图1为本发明方法的基于粒子群-蚁群并行交叉算法流程图;
图2为本发明方法的压电微定位平台实验装置图;
图3为本发明方法的压电微定位平台实测迟滞特性曲线和改进PI模型曲线对比图;
图4为本发明方法的压电微定位平台实际输出与模型输出误差曲线。
具体实施方式
粒子群和蚁群算法都是智能优化算法,能够有效辨识系统模型参数的能力,本发明结合两种算法,提出了一种该方法基于粒子群-蚁群并行交叉算法,便于发变化挥两种算法各自的优点,又能使结果避免陷入局部最优,得到更高精度的迟滞模型。
步骤一:根据传统PI模型和死区算子得到改进的PI模型:先依据传统单边play算子提出映射关系为曲线形式的改进单边play算子,然后串联死区算子得到非对称的改进单边play算子,进而加权叠加得到改进PI模型。
步骤二:设计粒子群-蚁群并行交叉算法辨识改进PI模型的参数:首先分别给出粒子群与蚁群算法辨识模型参数过程,为了综合发挥两种算法的优点随后提出粒子群-蚁群并行交叉算法用于辨识模型参数。
步骤三:搭建压电微定位平台采集辨识模型所需的数据,依据步骤2所述的粒子群-蚁群并行交叉算法辨识最终的模型参数并给出建模结果。
以下对本发明做进一步详细说明:
步骤一:设计改进的PI模型:
传统PI模型的离散形式为Play算子的加权和形式:
式中,u为压电微定位平台实际输入,y为PI模型的输出,k=1,2,...,K,K为采集压电微定位平台系统输入数据的总个数,取K=4001,n=1,2,...N,N代表PI模型中Play算子总数目,其中N=8,wn表示第n个Play算子的权值系数,是待辨识参数,Fn[u](k)为第n个单边Play算子,传统单边Play算子形式为:
式中:
rn表示第n个算子的阈值,0≤r1<r2<…<rN≤rmax,传统单边Play算子的输入输出映射关系曲线由直线构成且上升与下降部分关于算子中心对称,故为提高描述迟滞曲线精度就需要较多算子。为了实现较少算子情况下精确表述迟滞曲线,将传统PI模型进行改进。
首先提出改进Play算子,其映射关系为曲线形式,表达式如下:
其中wb(k)表示算子在变化过程中的斜率,wb(k)=Fn′[u](k-1)e-Δu(k),Δu(k)=u(k)-u(k-1)。
然后采用单边死区算子与改进Play算子串联结合,得到基于死区算子的非对称改进PI模型,其中单边死区算子形式如下:
最终得到改进的PI模型表达式如下:
步骤二:设计粒子群-蚁群并行交叉算法:
2.1本发明首先给出粒子群算法的具体步骤
粒子群优化算法是群体智能优化算法的一种,对于一个有M个个体的种群,其中每个个体都可以向D个方向自由运动,在算法中每个个体的坐标都有记录,用一个矩阵表示,个体的运动速度表示为选出其中距离目标最近的个体坐标作为最优的个体坐标,将代表当前种群中每个粒子最佳位置,gbest=[p1,p2,…,pD]是全局最佳位置,c1代表学习因子,表示对自身的继承能力;c2代表社会系数,表示对群体的学习能力;个体的坐标和速度更新限定在和内取随机值,G为最大迭代次数。
粒子群优化算法步骤如下:
(1)粒子的初始化:设定最大迭代次数G=10,粒子数M=50,待辨识参数个数D=16,学习因子和社会系数设定为c1=c2=2.732,每个粒子的初始位置和初始速度设为:
(2)计算每个粒子的适应度函数,并根据适应度函数选取第g次迭代的每个粒子最佳位置和全局最佳位置:
选取适应度函数:
其中yd(k)是压电陶瓷微定位平台的实际输出,y(k)是改进PI模型输出。
gbest(g)=min(pbestm(g))(10)
(3)每个粒子的位置和速度的更新公式为:
Vm(g+1)=Vm(g)+c1rand(pbestm(g)-Xm(g))+c2rand(gbest(g)-Xm(g))(11)
Xm(g+1)=Xm(g)+Vm(g+1)(12)
(4)根据适应度函数选取第g次迭代的每个粒子最佳位置和全局最佳位置,如果迭代次数达到设定的最大迭代次数,将gbest为最优解。否则,令迭代次数加一返回步骤(2)。
2.2本发明给出蚁群算法的具体步骤
蚁群算法是由意大利学者Dorigo、Maniezzo等人根据蚂蚁觅食规则受到启发于20世纪90年代提出,利用多次行进过程中根据释放的信息素能找到距离食物最近的路线。
蚁群算法的用于辨识模型参数的流程如下:
(1)参数初始化
蚁群算法在用于模型参数辨识时,首先确定待辨识参数个数为I个,那么蚁群需要经过I条相互平行且等距离间隔的线段上的节点,每条线段均等分为(J-1)分,即每条线段各有J个节点,节点值记为j=1,2,...,J。设定蚂蚁总数目为L,每只蚂蚁l经过的节点定义为一个矩阵其中h=1,2,...,H迭代次数,H为最大迭代次数。
改进PI模型中确定待辨识参数个数为I=16个,那么蚁群需要经过16条相互平行且等距离间隔的线段上的节点,每条线段均等分为(J-1)分,即每条线段各有J=100个节点,节点值记为j=1,2,...,J。设定蚂蚁总数目为L=50,每只蚂蚁l经过的节点定义为一个矩阵其中h=1,2,...,H迭代次数,H=10为最大迭代次数。其中待辨识参数的初始取值由如下公式给出:
(2)计算第l只蚂蚁的适应度函数,
最佳适应度
flmin(h)=min(fl(h)),l=1,2,...,L(15)
选取适应度函数值最小的蚂蚁经过的节点为当前迭代周期的最佳路径矩阵Zgood(h)=[z1(h),z2(h),...,zI(h)]。
第l只蚂蚁由第i-1条线段上任意节点向下一条线段上概率最大的节点移动。
(4)更新蚂蚁l经过的路径节点
ji为第h迭代中,第l只蚂蚁由第i-1条线段上任意节点向下一条线段上概率最大的节点值。
(5)更新信息素
表示第l只蚂蚁由第i-1条线段上任意节点向下一条线段上移动所路径上释放的信息素浓度增量,Δτ(zi(h),j)表示所有蚂蚁在第i-1条线段上所有节点向下一条线段上节点j移动连接路径上释放的信息素浓度之和。
(6)判断是否达到最大迭代次数H,若否,则当前迭代次数加1,返回步骤(2),若是,输出最佳路径矩阵为最优解。
2.3粒子群-蚁群并行交叉算法
本发明将粒子群算法与蚁群算法融合,设计粒子群-蚁群并行交叉算法便于发挥两种算法各自的优点又能相互补充,使结果避免陷入局部最优,又能使辨识精度更高,粒子群-蚁群算法交叉融合辨识运行过程:
(1)产生初始种群:在初始范围[a,b]=[0,1]内划分N=16个小区间,每个小区间随机产生两个数,记Q1,Q2,...,QN作为粒子群算法粒子初始位置,P1,P2,...,PN作为蚁群初始路径节点对应数值信息,粒子群-蚁群并行交叉算法的初始迭代周期T=1。
(2)依据粒子群算法与蚁群算法的执行步骤并行执行两种算法,各执行10代以后,利用公式(8)和(15)分别计算两种算法得到的新种群对应的适应度值。
(3)对两种算法得到粒子最佳位置和蚁群最佳路径节点对应的数值,将此时粒子群算法的位置数据计为S1,S2,...,SN,蚁群算法的路径节点数据记为R1,R2,...,RN;
(4)将R1,R2,…,Rξ与SN-ξ+1,SN-ξ+2,…,SN进行数据融合交换,令
S1,S2,..,Sξ,R1,R2,...,Rξ作为粒子群的新位置数据,SN-ξ+1,SN-ξ+2,...,SN,Rξ+1,Rξ+2,...,RN作为蚁群的新路径节点数据,ξ取1到N/2之间的随机整数;
(5)判断是否满足算法终止条件算法最大迭代次数70,若满足,则算法终止;否则,粒子群-蚁群并行交叉算法的初始迭代周期T加一,转到步骤(2)继续执行。
步骤三:搭建用于采集实际数据的压电微定位平台,如图2所示,采集辨识模型所需的输出输出数据,验证基于粒子群-蚁群并行交叉算法的改进PI模型辨识方法的有效性,得到最终实验结果。
本发明工作过程:
为了验证所建立改进PI模型的描述能力,我们搭建如图2所示的实验平台。整个控制台主要包括台式计算机、苏州博实机器人有限公司生产的精密定位控制器(PPC-2CR0150)和压电微定位平台(MPT-2MRL102A),研华科技公司的数据采集卡(PCI-1716l)。
首先在计算机运行MATLAB/SIMULINK搭建的实验环境中编写数据采集程序,通过外部模式即RTW工作环境编译控制程序发出实时输出信号;通过数据采集卡完成将计算机输出的数字信号转化为模拟信号;然后经过精密定位控制器将信号放大后驱动压电微定位平台使其产生位移,然后通过平台内部的位移传感器将压电微定位平台的位移信号采集,通过精密定位控制器将所测电压按一定倍数缩小后传递到数据采集卡,通过采集卡的模数转换功能转化成计算机能够处理的数字信号,以此构成实验闭环回路,得到压电微定位平台实际的位移数据。
利用采集程序给定最大电压为90V,频率为1Hz,四个周期的正弦电压信号激励压电微定位平台,进而采集到实际的输入输出数据,采用本发明所提出的基于粒子群-蚁群并行交叉算法对改进PI模型进行参数辨识,实测压电微定位平台迟滞特性曲线和模型曲线如图3所示,建模误差曲线如图4所示,由图可以得出最大建模误差的绝对值小于0.63μm,所建立的模型能够精确描述压电微定位平台的非对称迟滞特性。通过基于粒子群-蚁群并行交叉算法辨识可以得到改进PI模型的待辨识参数如表1所示:
表1基于粒子群-蚁群并行交叉算法辨识的改进PI模型参数
以下是本发明所涉及的符号:
Claims (1)
1.一种基于粒子群-蚁群并行交叉算法的改进PI模型辨识方法,其特征在于:其步骤是:
步骤一:根据传统PI模型和死区算子得到改进的PI模型:先依据传统单边play算子提出映射关系为曲线形式的改进单边play算子,然后串联死区算子得到非对称的改进单边play算子,进而加权叠加得到改进PI模型;
传统PI模型的离散形式为Play算子的加权和形式:
式中,k=1,2,...,K,K为系统输入数据个数,n=1,2,...N,N代表PI模型中Play算子总数目,wn表示第n个Play算子的权值系数,Fn[u](k)为第n个单边Play算子,u为压电微定位平台实际输入,传统单边Play算子形式为:
式中:
rn表示第n个算子的阈值,0≤r1<r2<…<rN≤rmax,将传统PI模型进行改进;
首先提出改进Play算子,其映射关系为曲线形式,表达式如下:
其中wb(k)表示算子在变化过程中的斜率,wb(k)=Fn′[u](k-1)e-Δu(k),Δu(k)=u(k)-u(k-1);然后采用单边死区算子与改进Play算子串联结合,得到基于死区算子的非对称改进PI模型,其中单边死区算子形式如下:
最终得到改进的PI模型表达式如下:
步骤二:设计粒子群-蚁群并行交叉算法辨识改进PI模型的参数:首先分别给出粒子群与蚁群算法辨识模型参数过程,为了综合发挥两种算法的优点随后提出粒子群-蚁群并行交叉算法用于辨识模型参数;
2.1粒子群优化
一个有M个个体的种群,其中每个个体都可以向D个方向自由运动,表示这个种群中第m个个体在D个方向上的坐标向量,其中和分别表示第m个个体在第1个方向,第2个方向和第D个方向上的坐标,m=1,2,...,M,个体的运动速度表示为选出其中距离目标最近的个体坐标作为最优的个体坐标,将代表当前种群中每个粒子最佳位置,gbest=[p1,p2,…,pD]是全局最佳位置,c1代表学习因子,c2代表社会系数,个体的坐标和速度更新限定在和内取随机值,G为最大迭代次数,粒子群优化算法步骤如下:
(1)粒子的初始化:
初始化每个粒子的位置和速度:
(2)计算每个粒子的适应度函数,并根据适应度函数选取第g次迭代的每个粒子最佳位置和全局最佳位置:
选取适应度函数:
gbest(g)=min(pbestm(g)) (10)
(3)每个粒子的位置和速度的更新公式为:
Vm(g+1)=Vm(g)+c1rand(pbestm(g)-Xm(g))+c2rand(gbest(g)-Xm(g)) (11)
Xm(g+1)=Xm(g)+Vm(g+1) (12)
(4)根据适应度函数选取第g次迭代的每个粒子最佳位置和全局最佳位置,如果迭代次数达到设定的值,此时将gbest为最优解;否则,令迭代次数加一返回步骤(2);
2.2蚁群算法
(1)参数初始化
蚁群算法在用于模型参数辨识时,首先确定待辨识参数个数为I个,那么蚁群需要经过I条相互平行且等距离间隔的线段上的节点,每条线段均等分为(J-1)分,即每条线段各有J个节点,节点值记为j=1,2,...,J,设定蚂蚁总数目为L,每只蚂蚁l经过的节点定义为一个矩阵其中h=1,2,...,H迭代次数,H为最大迭代次数,其中待辨识参数的初始取值由如下公式给出:
(2)计算第l只蚂蚁的适应度函数,
最佳适应度
flmin(h)=min(fl(h)),l=1,2,...,L (15)
选取适应度函数值最小的蚂蚁经过的节点为当前迭代周期的最佳路径矩阵Zgood(h)=[z1(h),z2(h),...,zI(h)];
(3)计算第l只蚂蚁由第i-1条线段上任意节点向下一条平行线上节点j移动的概率
第l只蚂蚁由第i-1条线段上任意节点向下一条线段上概率最大的节点移动;
(4)更新蚂蚁l经过的路径节点
ji为第h迭代中,第l只蚂蚁由第i-1条线段上任意节点向下一条线段上概率最大的节点值;
(5)更新信息素
表示第l只蚂蚁由第i-1条线段上任意节点向下一条线段上移动所路径上释放的信息素浓度增量,Δτ(zi(h),j)表示所有蚂蚁在第i-1条线段上所有节点向下一条线段上节点j移动连接路径上释放的信息素浓度之和;
(6)判断是否达到最大迭代次数H,若否,则当前迭代次数加一,返回步骤(2),若是,输出最优解;
2.3粒子群-蚁群并行交叉算法
(1)产生初始种群:在初始范围[a,b]内划分N个小区间,每个小区间随机产生两个数,记Q1,Q2,...,QN作为粒子群算法粒子初始位置,P1,P2,...,PN作为蚁群初始路径节点对应数值信息;
(2)依据粒子群算法与蚁群算法的执行步骤并行执行两种算法,执行10代以后,利用公式(8)和(15)分别计算两种算法得到的新种群对应的适应度值;
(3)对两种算法得到粒子最佳位置和蚁群最佳路径节点对应的数值,将此时粒子群算法的位置数据计为S1,S2,...,SN,蚁群算法的路径节点数据记为R1,R2,...,RN;
(4)将R1,R2,…,Rξ与SN-ξ+1,SN-ξ+2,…,SN进行数据融合交换,令S1,S2,..,Sξ,R1,R2,...,Rξ作为粒子群的新位置数据,SN-ξ+1,SN-ξ+2,...,SN,Rξ+1,Rξ+2,...,RN作为蚁群的新路径节点数据,ξ取1到N/2之间的随机整数;
(5)判断是否满足算法终止条件算法最大迭代次数70,若满足,则算法终止;否则,粒子群-蚁群并行交叉算法的初始迭代周期T加一,转到步骤(2)继续执行;
步骤三:搭建压电微定位平台采集辨识模型所需的数据,依据步骤二所述的粒子群-蚁群并行交叉算法辨识最终的模型参数并给出建模结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010481274.2A CN111914981B (zh) | 2020-05-31 | 2020-05-31 | 基于粒子群-蚁群并行交叉算法的改进pi模型辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010481274.2A CN111914981B (zh) | 2020-05-31 | 2020-05-31 | 基于粒子群-蚁群并行交叉算法的改进pi模型辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111914981A CN111914981A (zh) | 2020-11-10 |
CN111914981B true CN111914981B (zh) | 2022-11-08 |
Family
ID=73237885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010481274.2A Active CN111914981B (zh) | 2020-05-31 | 2020-05-31 | 基于粒子群-蚁群并行交叉算法的改进pi模型辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111914981B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112883632B (zh) * | 2021-01-13 | 2022-12-02 | 合肥工业大学 | 一种基于改进蚁群算法的锂电池等效电路模型参数辨识方法 |
CN112835295B (zh) * | 2021-01-22 | 2022-07-19 | 武汉工程大学 | 基于pi模型的压电陶瓷执行器参数辨识和复合控制方法 |
CN114077196B (zh) * | 2021-11-16 | 2023-06-06 | 哈尔滨工业大学(深圳) | 一种基于改进Prandtl-Ishlinskii模型的压电驱动器复合控制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105389623A (zh) * | 2015-10-22 | 2016-03-09 | 东南大学 | 一种采用改进蚁群算法的热工过程模型参数辨识方法 |
CN106125574A (zh) * | 2016-07-22 | 2016-11-16 | 吉林大学 | 基于dpi模型的压电陶瓷微定位平台建模方法 |
CN106485339A (zh) * | 2015-08-31 | 2017-03-08 | 中车大连电力牵引研发中心有限公司 | 一种电力系统的负荷特性确定方法及系统 |
CN107146411A (zh) * | 2017-06-06 | 2017-09-08 | 成都信息工程大学 | 基于方势阱模型的量子行为粒子群优化方法 |
CN108898257A (zh) * | 2018-07-05 | 2018-11-27 | 上海海事大学 | 一种基于粒子群优化的灰色预测方法 |
CN109116724A (zh) * | 2018-08-22 | 2019-01-01 | 上海电力学院 | 一种基于粒子群改进细菌觅食算法的负荷频率控制方法 |
CN110530384A (zh) * | 2019-07-17 | 2019-12-03 | 江苏大学 | 一种采用混合势场蚁群算法的车辆路径规划方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9037526B2 (en) * | 2012-12-13 | 2015-05-19 | Korea Institute Of Industrial Technology | Apparatus and method for selecting motion signifying artificial feeling |
-
2020
- 2020-05-31 CN CN202010481274.2A patent/CN111914981B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106485339A (zh) * | 2015-08-31 | 2017-03-08 | 中车大连电力牵引研发中心有限公司 | 一种电力系统的负荷特性确定方法及系统 |
CN105389623A (zh) * | 2015-10-22 | 2016-03-09 | 东南大学 | 一种采用改进蚁群算法的热工过程模型参数辨识方法 |
CN106125574A (zh) * | 2016-07-22 | 2016-11-16 | 吉林大学 | 基于dpi模型的压电陶瓷微定位平台建模方法 |
CN107146411A (zh) * | 2017-06-06 | 2017-09-08 | 成都信息工程大学 | 基于方势阱模型的量子行为粒子群优化方法 |
CN108898257A (zh) * | 2018-07-05 | 2018-11-27 | 上海海事大学 | 一种基于粒子群优化的灰色预测方法 |
CN109116724A (zh) * | 2018-08-22 | 2019-01-01 | 上海电力学院 | 一种基于粒子群改进细菌觅食算法的负荷频率控制方法 |
CN110530384A (zh) * | 2019-07-17 | 2019-12-03 | 江苏大学 | 一种采用混合势场蚁群算法的车辆路径规划方法 |
Non-Patent Citations (3)
Title |
---|
A Novel Hybrid Ant Colony-Particle Swarm Optimization Techniques Based Tuning STATCOM for Grid Code Compliance;Omar Makram Kamel 等;《IEEE Access》;20200227;第8卷;41566-41587 * |
压电微定位平台动态滑模自抗扰控制方法研究;杨立平;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20200815(第08期);C042-110 * |
超磁致径向微进给机构结构优化与迟滞建模及控制方法研究;刘孝亮;《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅰ辑》;20190115(第12期);B022-82 * |
Also Published As
Publication number | Publication date |
---|---|
CN111914981A (zh) | 2020-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111914981B (zh) | 基于粒子群-蚁群并行交叉算法的改进pi模型辨识方法 | |
Wang et al. | A new hysteresis modeling and optimization for piezoelectric actuators based on asymmetric Prandtl-Ishlinskii model | |
Yang et al. | Parameter identification of the generalized Prandtl–Ishlinskii model for piezoelectric actuators using modified particle swarm optimization | |
CN109918708B (zh) | 一种基于异质集成学习的材料性能预测模型构建方法 | |
Cai et al. | Different-level redundancy-resolution and its equivalent relationship analysis for robot manipulators using gradient-descent and Zhang's neural-dynamic methods | |
CN111931411A (zh) | 压电驱动微定位平台Duhem动态迟滞建模方法 | |
CN107919813A (zh) | 基于模糊神经网络的超声电机转速控制 | |
CN111368400B (zh) | 基于pso算法的压电微驱动变频定位平台建模辨识方法 | |
CN110245430A (zh) | 改进Bouc-Wen模型迟滞建模方法 | |
CN104991997A (zh) | 自适应差分进化算法优化的广义率相关p-i迟滞建模方法 | |
Lin et al. | Evolutionary algorithm based feedforward control for contouring of a biaxial piezo-actuated stage | |
Liu et al. | Neural network control system of cooperative robot based on genetic algorithms | |
Han et al. | Tracking control of ball and plate system using a improved PSO on-line training PID neural network | |
Chu et al. | Classifier-based approximator for friction compensation in high accelerated positioning system | |
CN113204918B (zh) | 水泥窑生产能耗的建模方法 | |
CN111554355B (zh) | 一种基于非冯诺依曼架构的分子动力学计算方法 | |
CN111753466B (zh) | 一种三极磁轴承转子径向位移软测量建模方法 | |
CN117252085A (zh) | 基于神经网络与粒子群的斜拉桥索力容差优化方法及系统 | |
CN115167150B (zh) | 具有未知系统动态的批次过程二维离轨策略交错q学习最优跟踪控制方法 | |
CN112363398B (zh) | 一种控制输入受限下的桥吊系统有限时间滑模控制系统及方法 | |
Li et al. | Parameter identification based on PSO algorithm for piezoelectric actuating system with rate-dependent Prandtl-Ishlinskii hysteresis modeling method | |
Sang et al. | Double Inverted Pendulum control based on three-loop PID and improved BP Neural network | |
CN101122779A (zh) | 一种利用蚁群算法辨识高精度伺服系统摩擦参数的方法 | |
Cai et al. | Application of improved wavelet neural network in MBR flux prediction | |
Son et al. | Black-box modeling of nonlinear system using evolutionary neural NARX model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |