CN113970905B - 用于高精度运动平台控制的任意阶s型曲线运动规划方法 - Google Patents
用于高精度运动平台控制的任意阶s型曲线运动规划方法 Download PDFInfo
- Publication number
- CN113970905B CN113970905B CN202111244799.5A CN202111244799A CN113970905B CN 113970905 B CN113970905 B CN 113970905B CN 202111244799 A CN202111244799 A CN 202111244799A CN 113970905 B CN113970905 B CN 113970905B
- Authority
- CN
- China
- Prior art keywords
- order
- curve
- motion planning
- shaped
- matrix
- 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
Classifications
-
- 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
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/416—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by control of velocity, acceleration or deceleration
-
- 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
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/34—Director, elements to supervisory
- G05B2219/34169—Coarse interpolator, path calculator delivers position, speed, acceleration blocks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种用于高精度运动平台控制的任意阶S型曲线运动规划方法,包括:对高精度运动平台控制所必须的规划指令,即任意阶S型运动规划曲线进行定义;针对所定义的任意阶S型运动规划曲线,构造该运动规划曲线的求导矩阵;根据高精度运动平台的规划距离、最大速度、最大加速度、速度连续、加速度连续等要求,建立等效的方程组,通过求解方程组得到任意阶次的S型曲线的表达式;利用该表达式进行高精度运动平台高阶S型曲线的运动规划。本发明所提出的方法可以为高精度运动平台的闭环控制提供任意阶次的S型曲线规划指令,便于编程实现,可以快速生成任意阶次的S型运动规划曲线。
Description
技术领域
本发明涉及高精度控制及制造领域,具体涉及一种用于高精度运动平台控制的任意阶S型曲线运动规划方法。
背景技术
在现代化制造领域中,高精度运动平台,如高精度直线平台、高精度旋转平台、SCARA机器人等得到广泛应用。高精度运动平台的控制策略一般采用闭环控制,闭环控制中常需要用到规划指令信号与反馈信号,尤其在进行对精度和速度有较高要求的加工任务时,如何在保持定位精度的前提下尽可能地缩短定位时间对于提高设备的生产效率和改善加工质量至关重要,解决这类问题的关键技术之一就是高速高精度运动规划。
在高速高精度运动规划中,高阶运动规划曲线对执行部件的速度平滑和性能的充分发挥有着极大的提升。然而受限于现有控制器或控制卡的功能,目前在工程领域中应用较多的运动规划曲线多为低阶的运动规划曲线,对于高阶运动规划曲线的推导构造鲜有研究涉及。
高精度运动平台的控制多需要结合控制器或者控制卡实现,目前应用较为广泛的有固高运动控制卡、PMAC控制器、dSPACE控制器等。目前,控制器或控制卡中使用最为广泛的是二阶运动规划曲线和三阶运动规划曲线,即梯形速度规划和三阶S型速度规划。通常而言,采用梯形速度规划可以极大地缩短运动规划的时间,但是梯形速度规划曲线存在加速度突变,容易激发较大的残余振动,进而导致运动系统需要较多的振动衰减时间来满足定位精度。所以一般在半导体制造等精密工程领域中,使用较多的为三阶S型速度规划,然而三阶S型速度规划曲线中,仍存在急动度不连续的情况。在一些场合中,高阶非对称曲线的规划方法可以得到相比于低阶运动规划曲线更好的效果。然而,目前工程中常用的控制器或控制卡并不具备高阶S型曲线运动规划的功能,大多情况下需要使用者自己写出高阶S型运动规划曲线的表达式,然后根据采样时间逐步计算规划指令,最后发送给控制器或控制卡实现闭环控制。
目前尚没有一个统一的方法来设计用于高精度运动平台控制的任意阶次S型运动规划曲线,这对于工程应用的编程实现极为不便。
发明内容
本发明的目的是提供一种用于高精度运动平台控制的任意阶S型曲线运动规划方法,用以解决现有技术中存在的问题。
为了实现上述任务,本发明采用以下技术方案:
一种用于高精度运动平台控制的任意阶S型曲线运动规划方法,包括以下步骤:
步骤1,对高精度运动平台控制所必须的规划指令,即任意阶S型运动规划曲线进行定义;
步骤2,针对所定义的任意阶S型运动规划曲线,构造该运动规划曲线的求导矩阵;
步骤3,根据高精度运动平台的规划距离、最大速度、最大加速度、速度连续、加速度连续等要求,建立等效的方程组,通过求解方程组得到任意阶次的S型曲线的表达式;利用该表达式进行高精度运动平台高阶S型曲线的运动规划。
进一步地,所述任意阶S型运动规划曲线进行定义,包括:
n阶S型运动规划曲线的位移表达式是2n-1段且每段均为n次多项式的表达式,Si(t)表示n阶S型运动规划曲线的第i段曲线的位移表达式:
Si(t)=ai,0+ai,1t+ai,2t2+ai,3t3+...+ai,ntn,i=1,2,3,...,(2n-1) (1)
其中,ai,0,ai,1,ai,2,...,ai,n表示第i段曲线位移表达式的n+1个系数,t表示时间。
进一步地,所述构造该运动规划曲线的求导矩阵,包括:
n阶S型运动规划曲线第i段位移曲线表达式Si(t)及其一系列导数的表达式如下:
其中,上标(n)表示n阶导数;用列向量x表示第i段曲线未知系数向量,如式(5)所示:
x=[ai,0,ai,1,ai,2,...,ai,n]T (5)
在进行高精度运动平台的控制时,所用规划指令中通常包含n阶S型运动规划曲线的导数,将每阶导数表达式中的各项系数、求导系数、时间t及其幂级数分开,形成未知系数向量x、时间矩阵Dn、以及求导系数矩阵Cn,其中Dn和Cn其中如式(6)和式(7)所示:
将式(4)拆分成由x、Dn、Cn经过特定运算得到的形式,如式(8)所示:
Dn(.*)Cn·x(8)
运算(.*)表示将Dn和Cn中相同位置索引的元素相乘;式(8)即为式(4)的矩阵形式,其中Dn(.*)Cn为n阶S型运动规划曲线的求导矩阵,记为Gn:
Gn=Dn(.*)Cn (9)
进一步地,所述根据高精度运动平台的规划距离、最大速度、最大加速度、速度连续、加速度连续等要求,建立等效的方程组,包括:
高精度运动平台控制所需的n阶S型运动规划曲线的位移表达式是2n-1段且每段均为n次多项式的表达式;每一位移分段表达式包含nc个未知数需要确定,则对于n阶S型运动规划曲线,需要确定的未知系数个数记为sn;
其中:
sn=ns×nc
通过n阶S型运动规划曲线的首尾时间节点处的约束条件、n阶S型运动规划曲线的0阶导至(n-1)阶导在第2个时间节点至第np-1个时间节点处的连续条件、以及n阶S型运动规划曲线的2阶导至n阶导的特殊约束条件来构造sn个方程。
进一步地,对于n阶S型运动规划曲线的首尾时间节点处的约束条件,表示如下:
将式(1)代入式(12)可得:
末时间节点处的约束条件可以用式(14)表示:
将式(1)代入式(14)得:
将式(13)和式(15)写成A1X=b1的形式,其中X为未知系数列向量:
A1为(2n×sn)的常数矩阵,b1为(2n×1)的常数列向量。利用式(9)的求导矩阵,可以较为方便的写出A1,对于第一个时间节点的约束条件,矩阵A1的部分元素为:
A1(j,1:n+1)表示矩阵A1第j行,第1列至第n+1列的元素,下同;G(j,:)|t=0表示时间t=0时的求导矩阵第j行的元素。对于末时间节点的约束条件:
表示/>时的求导矩阵第j行的元素;同时,b1部分元素为:
b1(n+1,1)=Q (18)
其中,时间节点向量
除了式(16)、(17)、(18)所指定的值,A1和b1的其他元素均为0。
进一步地,n阶S型运动规划曲线的0阶导至(n-1)阶导在第2个时间节点至第np-1个时间节点处的连续条件,表示为:
将式(1)代入式(19)得:
将式(19)写成矩阵A2X=b2的形式,A2为n×(np-2)行,sn列的常数矩阵,b2为[n×(np-2)]×1的常数列向量;利用式(9)的求导矩阵,可以较为方便的写出A2:
表示时间t=Tj的求导矩阵第k行的元素;除了式(21)所指定的元素外,矩阵A2其余元素均为0,b2的元素均为0。
进一步地,所述阶S型运动规划曲线的2阶导至n阶导的特殊约束条件,表示为:
n阶S型运动规划曲线的2阶导至n阶导存在恒等于0的曲线段,且2阶导至n阶导恒等于0的曲线段索引与曲线的阶数和导数有关;通过归纳总结可知,n阶S型运动规划曲线的g阶导数中恒等于0的曲线段索引为:
向量Vn,g为由n阶S型运动规划曲线的g阶导数中恒等于0的曲线段索引组成的行向量,Ln,g为向量Vn,g的长度;
将以上约束条件同样写成矩阵形式A3X=b3,常数向量b3的元素均为0,常数矩阵A3的元素由0和1组成,当n阶S型运动规划曲线的第i段曲线段的n阶导恒等于0时,则该曲线段的系数ai,n等于0;(n-1)阶导恒等于0时,系数ai,n-1和ai,n等于0;(n-2)阶导恒等于0时,系数ai,n-2、ai,n-1、ai,n等于0,以此类推;因此,若想让向量X中的某些参数等于0,则让矩阵A3中对应位置的元素值为1即可;经过归纳总结,可以得到A3:
其中,Bg为n阶S型运动规划曲线第g阶导中恒等于0的方程组构成的常数矩阵,Vn,g(k)表示向量Vn,g的第k个元素。
进一步地,通过求解方程组得到任意阶次的S型曲线的表达式,包括:
求得的A1、A2、A3以及b1、b2、b3,可以构造矩阵A和向量b:
通过下式即可求得未知系数向量X,则n阶S型运动规划曲线便构造完毕:
X=A-1b (29)
进一步地,所述利用该表达式进行高精度运动平台高阶S型曲线的运动规划,包括:
在求得系数向量X之后,便可以求得n阶S型运动规划曲线的表达式,根据式(9)可以求得其一系列导数,然后根据采样时间dt计算每个控制周期内发送给控制器的规划指令,从而完成高精度运动平台高阶S型曲线的运动规划。
一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,处理器执行计算机程序时实现用于高精度运动平台控制的任意阶S型曲线运动规划方法的步骤。
一种计算机可读存储介质,计算机可读存储介质存储有计算机程序,计算机程序被处理器执行时实现用于高精度运动平台控制的任意阶S型曲线运动规划方法的步骤。
与现有技术相比,本发明具有以下技术特点:
1.本发明所提出的方法可以为高精度运动平台的闭环控制提供任意阶次的S型曲线规划指令,且对任意阶次S型运动规划曲线的构造提供了一种通用的构造方法。
2.本发明所提出的构造n阶S型运动规划曲线的方法便于编程实现,可以快速生成任意阶次的S型运动规划曲线。
附图说明
图1为高精度运动平台闭环控制图;图中S(t)为n阶S型曲线运动规划指令;sc为位移指令;vc为速度指令;ac为加速度指令;sd为位移反馈信号;vd为速度反馈信号;ad为加速度反馈信号。
图2为二阶S型运动规划曲线图,图中T1、T2、T3为时间节点。
图3为三阶S型运动规划曲线图,图中T1、T2、T3、T4、T5、T6、T7为时间节点。
图4为四阶S型运动规划曲线图,图中T1、T2、T3、T4、T5、T6、T7、T8、T9、T10、T11、T12、T13、T14、T15为时间节点。
图5为五阶S型运动规划曲线图。
具体实施方式
本发明针对dSPACE等控制器可以自定义规划曲线的特点,从任意阶次S型运动规划曲线的特性和工程实践性出发,对任意阶次的S型运动规划曲线进行定义,对运动规划曲线中需要满足的连续性要求和特殊约束要求采用理论推导与归纳的方式进行总结,从而获得构造曲线的所必须的各项系数。采用本发明所提出的方法,仅需少量的代码即可快速生成用于高精度运动平台控制的任意阶次的S型运动规划曲线。
参见附图,本发明提出的一种用于高精度运动平台控制的任意阶S型曲线运动规划方法,包括以下步骤:
步骤1,对高精度运动平台控制所必须的规划指令,即任意阶S型运动规划曲线进行定义。
高精度运动平台的闭环控制框图大致可以用图1进行描述,其中规划指令的设计是闭环控制中重要的一环。目前工程中常用的控制卡或控制器中一般只具备二阶S型曲线运动规划和三阶S型运动规划的功能,而对于更高阶的S型曲线运动规划的功能比较欠缺。通过控制器或控制卡进行高精度运动平台的高阶S型曲线运动规划,通常的做法是先生成任意阶次的S型运动规划曲线,该曲线是关于时间t的表达式,然后根据采样时间dt计算出每一个控制周期里的规划指令值。本发明所要解决的问题便是图1中的规划指令S(t)的生成问题。
在高精度运动平台的控制中除了用到位置规划指令,经常也会用到速度规划指令和加速度规划指令等作为指令前馈,这可以大大提高平台的定位精度和运动的响应速度。速度规划指令和加速度规划指令是位置规划指令的一阶导数和二阶导数,因此在进行运动规划时,最重要的是求得位置规划指令的表达式。二阶运动规划曲线如图2所示,其位移表达式是3段的分段表达式,每段表达式是关于时间t的2次多项式;三阶运动规划曲线如图3所示,其位移表达式是7段的分段表达式,每段表达式是关于时间t的3次多项式。因此任意阶S型曲线运动规划的位移表达式可以用分段的n次多项式进行定义。Si(t)表示n阶S型运动规划曲线的第i段曲线的位移表达式,如式(1)所示。
Si(t)=ai,0+ai,1t+ai,2t2+ai,3t3+...+ai,ntn,i=1,2,3,...,(2n-1) (1)
其中,ai,0,ai,1,ai,2,...,ai,n表示第i段曲线位移表达式的n+1个系数,当该段曲线的系数确定时,则该段曲线也确定。将n阶S型运动规划曲线的段数记为ns;每段曲线的未知系数的个数记为nc;总的时间节点个数记为np。可以知道:
时间节点包含运动时间起点与终点,以及第i段和第i+1段曲线交点处所对应的时间,通常用向量T表示时间节点向量,如式(3)所示。
步骤2,针对步骤1中定义的任意阶S型运动规划曲线,构造该运动规划曲线的求导矩阵。
在进行高精度运动平台的控制时,通常需要用到位移信息S(t)、速度信息V(t)、加速度信息A(t)、加加速度信息J(t)等,即位移曲线S(t)及其一系列导数等。n阶S型运动规划曲线第i段位移曲线表达式Si(t)及其一系列导数的表达式如式(4)所示:
用列向量x表示第i段曲线未知系数向量,如式(5)所示。
x=[ai,0,ai,1,ai,2,...,ai,n]T (5)
在进行高精度运动平台的控制时,所用规划指令中通常包含n阶S型运动规划曲线的导数,将每阶导数表达式中的各项系数、求导系数、时间t及其幂级数分开,形成未知系数向量x、时间矩阵Dn、以及求导系数矩阵Cn,其中Dn和Cn其中如式(6)和式(7)所示:
Dn,Cn均为(n+1)×(n+1)的矩阵。为方便后续编程实现,将式(4)拆分成由x、Dn、Cn经过特定运算得到的形式,如式(8)所示:
Dn(.*)Cn·x (8)
运算(.*)表示将Dn和Cn中相同位置索引的元素相乘。式(8)即为式(4)的矩阵形式,其中Dn(.*)Cn为n阶S型运动规划曲线的求导矩阵,记为Gn:
Gn=Dn(.*)Cn (9)
步骤3,最后根据高精度运动平台的规划距离、最大速度、最大加速度、速度连续、加速度连续等要求,建立等效的方程组,通过求解方程组得到任意阶次的S型曲线的表达式。
从步骤1可知,高精度运动平台控制所需的n阶S型运动规划曲线的位移表达式是2n-1段且每段均为n次多项式的表达式。每一位移分段表达式包含nc个未知数需要确定,则对于n阶S型运动规划曲线,需要确定的未知系数个数sn为:
sn=ns×nc (10)
用列向量X表示n阶S型运动规划曲线未知系数向量,如式(11)所示:
因此如需确定n阶S型运动规划曲线,则需确定向量X,X向量中包含sn个未知数,即需要构造sn个方程。本发明通过n阶S型运动规划曲线的首尾时间节点处的约束条件、n阶S型运动规划曲线的0阶导至(n-1)阶导在第2个时间节点至第np-1个时间节点处的连续条件、以及n阶S型运动规划曲线的2阶导至n阶导的特殊约束条件来构造sn个方程。
1)n阶S型运动规划曲线的首尾时间节点处的约束条件,即在第一个时间节点处,高精度运动平台的位移、速度、加速度、加加速度等均等于0,即n阶S型运动规划曲线的0阶导至n-1导均等于0;在末时间节点处,高精度运动平台的规划距离等于Q,而其速度、加速度、加加速度等均等于0,即n阶S型运动规划曲线的0阶导等于规划的距离Q,1阶导至n-1阶导均等于0。第一个时间节点的约束条件可以用式(12)表示:
将式(1)代入式(12)可得:
末时间节点处的约束条件可以用式(14)表示:
将式(1)代入式(14)得:
将式(13)和式(15)写成A1X=b1的形式,其中X为未知系数列向量,A1为(2n×sn)的常数矩阵,b1为(2n×1)的常数列向量。利用式(9)的求导矩阵,可以较为方便的写出A1,对于第一个时间节点的约束条件,矩阵A1的部分元素为:
A1(j,1:n+1)表示矩阵A1第j行,第1列至第n+1列的元素,下同;G(j,:)|t=0表示时间t=0时的求导矩阵第j行的元素。对于末时间节点的约束条件:
表示/>时的求导矩阵第j行的元素。同时,b1部分元素为:
b1(n+1,1)=Q (18)
除了式(16)、(17)、(18)所指定的值,A1和b1的其他元素均为0。
2)高精度运动平台在高速运动过程中,如若要保持较高的精度,除了其位移曲线需保证在运动时间内连续外,其速度、加速度等曲线也需保证连续性。如图2所示的二阶S型曲线运动规划,其加速度段在时间节点处不连续,容易导致高精度运动平台在运动过程中产生振动,从而丢失定位精度。高精度运动平台要保证较高的定位精度和加减速的平滑性,需保证其规划指令n阶S型运动规划曲线的0阶导至(n-1)阶导在第2个时间节点至第np-1个时间节点处的连续条件,可以用式(19)表示:
将式(1)代入式(19)得:
将式(19)写成矩阵A2X=b2的形式,A2为n×(np-2)行,sn列的常数矩阵,b2为[n×(np-2)]×1的常数列向量。利用式(9)的求导矩阵,可以较为方便的写出A2:
表示时间t=Tj的求导矩阵第k行的元素。除了式(21)所指定的元素外,矩阵A2其余元素均为0,b2的元素均为0。
3)为了充分发挥高精度运动平台电机的性能,规划指令n阶S型运动规划曲线的2阶导至n阶导需要满足特殊约束条件,即n阶S型运动规划曲线的2阶导至n阶导在某些曲线段恒等于0,而这些曲线段的索引存在一定的规律。通过二阶运动规划曲线、三阶运动规划曲线、四阶运动规划曲线中2阶导至n阶导恒等于0的曲线段索引的规律进行归纳总结,从而推广至n阶运动规划曲线的情况。
二阶运动规划曲线如图1所示,可以看出,加速度曲线的第2段曲线恒等于0,即:
三阶运动规划曲线如图2所示,可以看出,加速度曲线的第4段曲线恒等于0,加加速度曲线的第2、4、6段曲线即/>恒等于0,即:
四阶运动规划曲线如图3所示,同样可得,加速度曲线的第8段曲线恒等于0,加加速度曲线的第4、8、12段曲线恒等于0,四阶导的第2、4、6、8、10、12、14段曲线恒等于0,即:
从式(22)、(23)、(24)可知,2阶导至n阶导恒等于0的曲线段索引与曲线的阶数和导数有关。通过归纳总结可知,n阶S型运动规划曲线的g阶导数中恒等于0的曲线段索引为:
向量Vn,g为由n阶S型运动规划曲线的g阶导数中恒等于0的曲线段索引组成的行向量,Ln,g为向量Vn,g的长度。
将以上约束条件同样写成矩阵形式A3X=b3,常数向量b3的元素均为0,常数矩阵A3的元素由0和1组成。当n阶S型运动规划曲线的第i段曲线段的n阶导恒等于0时,则该曲线段的系数ai,n等于0;(n-1)阶导恒等于0时,系数ai,n-1和ai,n等于0;(n-2)阶导恒等于0时,系数ai,n-2、ai,n-1、ai,n等于0,以此类推。因此,若想让向量X中的某些参数等于0,则让矩阵A3中对应位置的元素值为1即可。经过归纳总结,可以得到A3:
其中,Vn,g(k)表示向量Vn,g的第k个元素;Bg为n阶S型运动规划曲线第g阶导中恒等于0的方程组构成的常数矩阵。
综合上述三种情况求得的A1、A2、A3以及b1、b2、b3,可以构造矩阵A和向量b:
通过下式即可求得未知系数向量X,则n阶S型运动规划曲线便构造完毕。
X=A-1b (29)
在求得系数向量X之后,便可以求得n阶S型运动规划曲线的表达式,根据式(9)可以求得其一系列导数,然后根据采样时间dt计算每个控制周期内发送给控制器的规划指令,从而完成高精度运动平台高阶S型曲线的运动规划。
实施例:构造用于高精度直线运动平台控制的五阶S型运动规划曲线
通过本发明的方法生成任意阶的S型运动规划曲线时,需要给定规划的距离Q和时间节点向量T。假设高精度直线平台运动的距离为0.15m,运动规划曲线的阶次为5阶,则根据式(2)可得:
令时间节点向量T为:
通过式(5)(7)计算时间矩阵D5和求导系数矩阵C5
然后通过式(16)(17)(21)即可计算出A1和A2。通过式(25)计算V5,5、V5,4、V5,3、V5,2以及L5,5、L5,4、L5,3、L5,2:
通过式(26)即可计算出B2、B3、B4、B5,然后通过式(27)即可得出A3。最后通过式(28)(29)即可计算所有的未知系数。图4为构造出的五阶S型运动规划曲线。
以上实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围,均应包含在本申请的保护范围之内。
Claims (3)
1.一种用于高精度运动平台控制的任意阶S型曲线运动规划方法,其特征在于,包括以下步骤:
步骤1,对高精度运动平台控制所必须的规划指令,即任意阶S型运动规划曲线进行定义;
步骤2,针对所定义的任意阶S型运动规划曲线,构造该运动规划曲线的求导矩阵;
步骤3,根据高精度运动平台的规划距离、最大速度、最大加速度、速度连续、加速度连续要求,建立等效的方程组,通过求解方程组得到任意阶次的S型曲线的表达式;利用该表达式进行高精度运动平台高阶S型曲线的运动规划;
对所述任意阶S型运动规划曲线进行定义,包括:
n阶S型运动规划曲线的位移表达式是2n-1段且每段均为n次多项式的表达式,Si(t)表示n阶S型运动规划曲线的第i段曲线的位移表达式:
Si(t)=ai,0+ai,1t+ai,2t2+ai,3t3+...+ai,ntn,i=1,2,3,...,2n-1 (1)
其中,ai,0,ai,1,ai,2,...,ai,n表示第i段曲线位移表达式的n+1个系数,t表示时间;
所述构造该运动规划曲线的求导矩阵,包括:
n阶S型运动规划曲线第i段位移曲线表达式Si(t)及其一系列导数的表达式如下:
其中,上标n表示n阶导数;用列向量x表示第i段曲线未知系数向量,如式(5)所示:
x=[ai,0,ai,1,ai,2,...,ai,n]T (5)
在进行高精度运动平台的控制时,所用规划指令中包含n阶S型运动规划曲线的导数,将每阶导数表达式中的各项系数、求导系数、时间t及其幂级数分开,形成未知系数向量x、时间矩阵Dn、以及求导系数矩阵Cn,其中Dn和Cn如式(6)和式(7)所示:
将式(4)拆分成由x、Dn、Cn经过特定运算得到的形式,如式(8)所示:
Dn(.*)Cn·x(8)
运算(.*)表示将Dn和Cn中相同位置索引的元素相乘;式(8)即为式(4)的矩阵形式,其中Dn(.*)Cn为n阶S型运动规划曲线的求导矩阵,记为Gn:
Gn=Dn(.*)Cn (9)
所述根据高精度运动平台的规划距离、最大速度、最大加速度、速度连续、加速度连续要求,建立等效的方程组,包括:
高精度运动平台控制所需的n阶S型运动规划曲线的位移表达式是2n-1段且每段均为n次多项式的表达式;将n阶S型运动规划曲线的段数记为ns,总的时间节点个数记为np,每一位移分段表达式包含nc个未知系数需要确定,则对于n阶S型运动规划曲线,需要确定的未知系数个数记为sn;
其中:
sn=ns×nc
通过n阶S型运动规划曲线的首尾时间节点处的约束条件、n阶S型运动规划曲线的0阶导至n-1阶导在第2个时间节点至第np-1个时间节点处的连续条件、以及n阶S型运动规划曲线的2阶导至n阶导的特殊约束条件来构造sn个方程;
对于n阶S型运动规划曲线的首尾时间节点处的约束条件,表示如下:
将式(1)代入式(12)可得:
末时间节点处的约束条件可以用式(14)表示:
将式(1)代入式(14)得:
将式(13)和式(15)写成A1X=b1的形式,其中X为未知系数列向量:
A1为2n×sn的常数矩阵,b1为2n×1的常数列向量;利用式(9)的求导矩阵,写出A1,对于第一个时间节点的约束条件,矩阵A1的部分元素为:
A1(j,1:n+1)表示矩阵A1第j行,第1列至第n+1列的元素,下同;G(j,:)|t=0表示时间t=0时的求导矩阵第j行的元素;对于末时间节点的约束条件:
表示/>时的求导矩阵第j行的元素;同时,b1部分元素为:
b1(n+1,1)=Q (18)
其中,时间节点向量
除了式(16)、(17)、(18)所指定的值,A1和b1的其他元素均为0;
n阶S型运动规划曲线的0阶导至n-1阶导在第2个时间节点至第np-1个时间节点处的连续条件,表示为:
将式(1)代入式(19)得:
将式(19)写成矩阵A2X=b2的形式,A2为n×(np-2)行,sn列的常数矩阵,b2为[n×(np-2)]×1的常数列向量;利用式(9)的求导矩阵,写出A2:
表示时间t=Tj的求导矩阵第k行的元素;除了式(21)所指定的元素外,矩阵A2其余元素均为0,b2的元素均为0;
所述n阶S型运动规划曲线的2阶导至n阶导的特殊约束条件,表示为:
n阶S型运动规划曲线的2阶导至n阶导存在恒等于0的曲线段,且2阶导至n阶导恒等于0的曲线段索引与曲线的阶数和导数有关;通过归纳总结可知,n阶S型运动规划曲线的g阶导数中恒等于0的曲线段索引为:
向量Vn,g为由n阶S型运动规划曲线的g阶导数中恒等于0的曲线段索引组成的行向量,Ln,g为向量Vn,g的长度;
将以上约束条件同样写成矩阵形式A3X=b3,常数向量b3的元素均为0,常数矩阵A3的元素由0和1组成,当n阶S型运动规划曲线的第i段曲线段的n阶导恒等于0时,则该曲线段的系数ai,n等于0;n-1阶导恒等于0时,系数ai,n-1和ai,n等于0;n-2阶导恒等于0时,系数ai,n-2、ai,n-1、ai,n等于0,以此类推;因此,若想让向量X中的某些参数等于0,则让矩阵A3中对应位置的元素值为1即可;经过归纳总结,可以得到A3:
其中,Bg为n阶S型运动规划曲线第g阶导中恒等于0的方程组构成的常数矩阵,Vn,g(k)表示向量Vn,g的第k个元素;
通过求解方程组得到任意阶次的S型曲线的表达式,包括:
求得的A1、A2、A3以及b1、b2、b3,可以构造矩阵A和向量b:
通过下式即可求得未知系数向量X,则n阶S型运动规划曲线便构造完毕:
X=A-1b (29)。
2.根据权利要求1所述的用于高精度运动平台控制的任意阶S型曲线运动规划方法,其特征在于,所述利用该表达式进行高精度运动平台高阶S型曲线的运动规划,包括:
在求得系数向量X之后,便可以求得n阶S型运动规划曲线的表达式,根据式(9)可以求得其一系列导数,然后根据采样时间dt计算每个控制周期内发送给控制器的规划指令,从而完成高精度运动平台高阶S型曲线的运动规划。
3.一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,处理器执行计算机程序时实现根据权利要求1-2中任一权利要求所述的用于高精度运动平台控制的任意阶S型曲线运动规划方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111244799.5A CN113970905B (zh) | 2021-10-26 | 2021-10-26 | 用于高精度运动平台控制的任意阶s型曲线运动规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111244799.5A CN113970905B (zh) | 2021-10-26 | 2021-10-26 | 用于高精度运动平台控制的任意阶s型曲线运动规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113970905A CN113970905A (zh) | 2022-01-25 |
CN113970905B true CN113970905B (zh) | 2023-07-18 |
Family
ID=79588514
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111244799.5A Active CN113970905B (zh) | 2021-10-26 | 2021-10-26 | 用于高精度运动平台控制的任意阶s型曲线运动规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113970905B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104008250A (zh) * | 2014-06-10 | 2014-08-27 | 广东工业大学 | 基于主频能量时域最优分布的非对称变加速度规划方法 |
CN105182906A (zh) * | 2015-09-24 | 2015-12-23 | 哈尔滨工业大学 | 基于高阶s型运动轨迹的位置与速度控制方法 |
CN106393106A (zh) * | 2016-09-13 | 2017-02-15 | 东南大学 | 参数自适应密化的机器人nurbs曲线运动插补方法 |
CN111596543A (zh) * | 2020-06-10 | 2020-08-28 | 嘉兴市像景智能装备有限公司 | 一种快速规划6阶s型运动曲线的方法 |
CN112506143A (zh) * | 2020-12-07 | 2021-03-16 | 大连大学 | 基于s型曲线的高质量加工五次多项式速度规划方法 |
CN113253678A (zh) * | 2021-04-29 | 2021-08-13 | 中国科学院微电子研究所 | 精密运动平台的运动轨迹规划方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7479751B2 (en) * | 2007-01-29 | 2009-01-20 | Rockwell Automation Technologies, Inc. | Elimination of unintended velocity reversals in s-curve velocity profiles |
-
2021
- 2021-10-26 CN CN202111244799.5A patent/CN113970905B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104008250A (zh) * | 2014-06-10 | 2014-08-27 | 广东工业大学 | 基于主频能量时域最优分布的非对称变加速度规划方法 |
CN105182906A (zh) * | 2015-09-24 | 2015-12-23 | 哈尔滨工业大学 | 基于高阶s型运动轨迹的位置与速度控制方法 |
CN106393106A (zh) * | 2016-09-13 | 2017-02-15 | 东南大学 | 参数自适应密化的机器人nurbs曲线运动插补方法 |
CN111596543A (zh) * | 2020-06-10 | 2020-08-28 | 嘉兴市像景智能装备有限公司 | 一种快速规划6阶s型运动曲线的方法 |
CN112506143A (zh) * | 2020-12-07 | 2021-03-16 | 大连大学 | 基于s型曲线的高质量加工五次多项式速度规划方法 |
CN113253678A (zh) * | 2021-04-29 | 2021-08-13 | 中国科学院微电子研究所 | 精密运动平台的运动轨迹规划方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113970905A (zh) | 2022-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10747193B2 (en) | Machine learning apparatus, servo control apparatus, servo control system, and machine learning method | |
Bristow et al. | A survey of iterative learning control | |
Chen et al. | Composite nonlinear feedback control for linear systems with input saturation: theory and an application | |
US11144017B2 (en) | Quadratic program solver for MPC using variable ordering | |
US10877442B2 (en) | Machine learning device, control device, and machine learning method | |
JP2017102617A (ja) | 補正装置、補正装置の制御方法、情報処理プログラム、および記録媒体 | |
CN108983615B (zh) | 基于反双曲正弦吸引律的离散双周期重复控制器 | |
CN113031439B (zh) | 一种双运动台精密协同控制系统及方法 | |
WO2022160894A1 (zh) | 基于迭代学习的数控机床快速误差补偿控制系统及方法 | |
Abdul-Adheem et al. | From PID to nonlinear state error feedback controller | |
Chang et al. | Analysis and design of two types of digital repetitive control systems | |
CN113253678B (zh) | 精密运动平台的运动轨迹规划方法 | |
CN113970905B (zh) | 用于高精度运动平台控制的任意阶s型曲线运动规划方法 | |
Kim et al. | Asymptotic state tracking in a class of nonlinear systems via learning-based inversion | |
CN113219907A (zh) | 一种精密运动台扰动力辨识系统与方法 | |
Coelingh | Design support for motion control systems: a mechatronic approach | |
CN1367920A (zh) | 用于在传动装置饱和时维持伺服稳定性的方法和装置 | |
CN1110729C (zh) | 基于对象特征模型描述的黄金分割智能控制方法 | |
Chen et al. | Nonlinear contour tracking of a voice coil motors-driven dual-axis positioning stage using fuzzy fractional PID control with variable orders | |
CN111522225A (zh) | 基于改进的粒子群算法对调阻机pid控制系统参数优化方法 | |
CN113156894B (zh) | 芯片封装控制的速度规划方法及芯片封装控制装置 | |
CN115167111A (zh) | 基于改进型分数阶pid的伺服电机控制方法及系统 | |
Li et al. | High precision position control of voice coil motor based on single neuron PID | |
CN112757299B (zh) | 质心轨迹的获取方法、装置、机器人及存储介质 | |
CN112975986A (zh) | 基于径向基函数的机械臂点对点轨迹规划方法及装置 |
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 |