CN102236338A - 四阶定值发生器的实现方法 - Google Patents

四阶定值发生器的实现方法 Download PDF

Info

Publication number
CN102236338A
CN102236338A CN2010101536874A CN201010153687A CN102236338A CN 102236338 A CN102236338 A CN 102236338A CN 2010101536874 A CN2010101536874 A CN 2010101536874A CN 201010153687 A CN201010153687 A CN 201010153687A CN 102236338 A CN102236338 A CN 102236338A
Authority
CN
China
Prior art keywords
time period
acceleration
displacement
constraint
speed
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.)
Granted
Application number
CN2010101536874A
Other languages
English (en)
Other versions
CN102236338B (zh
Inventor
张志钢
王强
石勋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Micro Electronics Equipment Co Ltd
Original Assignee
Shanghai Micro Electronics Equipment Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Micro Electronics Equipment Co Ltd filed Critical Shanghai Micro Electronics Equipment Co Ltd
Priority to CN2010101536874A priority Critical patent/CN102236338B/zh
Publication of CN102236338A publication Critical patent/CN102236338A/zh
Application granted granted Critical
Publication of CN102236338B publication Critical patent/CN102236338B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)

Abstract

一种四阶定值发生器的实现方法,用以产生满足约束的四阶定值轨迹,其包括:定义djerk时间段Td、jerk时间段Tj、加速度时间段Ta、速度时间段Tv,其中jerk为加速度的微分,djerk为jerk的微分,该过程包括:令Td=td,其中td为常djerk时间段;令Tj=td+tj,其中tj为常jerk时间段;令Ta=td+tj+ta,其中ta为常加速度时间段;令Tv=td+tj+ta+tv,其中tv为常速度时间段;根据期望的位移以及djerk、jerk、加速度和速度的最大值约束求出Td、Tj、Ta、Tv的值;得到满足约束的四阶定值轨迹。以上方法采用了不同于现有技术的时间段定义方式,进而采用非递推计算过程来实现四阶定值发生器,其相对于现有技术,实现过程简洁,计算量小。

Description

四阶定值发生器的实现方法
技术领域
本发明涉及伺服控制技术领域,特别是涉及一种四阶定值发生器的实现方法。
背景技术
为满足某种目的,产生运动和对物体运动进行控制是人类最重要的活动之一。所谓伺服控制就是指对物体运动的有效控制,即对物体运动的速度、位置、加速度进行控制。这种控制正在变得越来越普遍。
目前,随着机电制造的发展,对于高精度、高速度的伺服控制的要求越来越高。但是当前制约伺服控制性能的一个主要问题在于运动轨迹的加速度和减速度的突变会激发系统的谐振,增加建立时间(settling time)。针对这个问题,许多学者提出了各种定值发生器(set point generator),通过定值发生器来产生更加平滑的加速度轨迹,来避免激发对象的高频谐振模态。
例如,2003年2月公开的论文“Trajectory planning and feedforward design forelectromechanical motion systems”揭露了一种三阶定值发生器的实现方法;2004年8月公开的公开号为US2004/0166425的美国专利申请揭露了一种三阶定值发生器。这种三阶定值发生器虽然可以使得加速度的曲线变得平滑,但是加速度的微分(jerk)曲线仍然是不连续的。而jerk的突变也会导致较大的建立时间,降低生产效率。因此,很多学者开始致力于研究更高阶的定值发生器,如四阶定值发生器。
例如,2003年2月公开的论文“Traj ectory planning and feedforward design forelectromechanical motion systems”揭露了一种四阶定值发生器的实现方法;2005年6月公开的公开号为US2005/0128460的美国专利申请揭露了一种四阶定值发生器。其中提出虽然阶次升高会导致整个轨迹的执行时间加长,但是由于jerk曲线更加平滑,减小了建立时间,从而提高了生产效率。
四阶定值发生器的目的是要根据期望的位移得到时间最短且满足jerk的微分(djerk)、jerk、加速度a、速度v等最大值约束的运动轨迹。设
Figure GSA00000094645800022
Figure GSA00000094645800023
分别为djerk、jerk、加速度和速度的最大值约束。标准的四阶轨迹如图1所示,其中曲线101为djerk曲线,呈方波形状;对djerk曲线101进行积分,可得曲线102所示的jerk曲线,该曲线呈梯形;对jerk曲线102进行积分,可得曲线103所示的加速度曲线,此时的加速度曲线103完全平滑,可以极大的减少由于加速度突变所产生的振荡,减小建立时间;继续对加速度曲线103进行积分就可得曲线104所示的速度曲线,对速度曲线104再进行积分可得位移曲线105。
下面简要回顾一下以上四阶定值发生器的实现方法:
首先,按照如图2所示的方式定义时间段,其中td为djerk保持最大值的时间段,称为常djerk时间段;tj为jerk等于
Figure GSA00000094645800026
的时间段,称为常jerk时间段;ta为加速度等于
Figure GSA00000094645800027
的时间段,称为常加速度时间段;tv是速度等于的时间段,称为常速度时间段。然而,根据位移以及对djerk、jerk、加速度和速度的最大值约束求出四个时间段td、tj、ta和tv的值,然后就可以得到djerk的曲线。而后,再通过多次积分,可得速度位移等曲线。
其中,根据位移以及对djerk、jerk、加速度和速度的最大值约束求出四个时间段td、tj、ta和tv的值的过程的具体流程如图3所示,包括如下步骤:
步骤S1:首先,设立一个标记变量flag,并令flag=0,根据位移x计算td
步骤S2:根据td计算轨迹的最大速度vmax
步骤S3:比较vmax
Figure GSA00000094645800029
的大小,当时,执行步骤S4,否则执行步骤S5;
步骤S4:根据
Figure GSA000000946458000211
计算td,并令flag=1;
步骤S5:根据td计算轨迹的最大加速度amax
步骤S6:比较amax的大小,当
Figure GSA000000946458000213
时,执行步骤S7,否则执行步骤S8;
步骤S7:根据
Figure GSA000000946458000214
计算td,并令flag=2;
步骤S8:根据td计算轨迹的最大jerk,jmax
步骤S9:比较jmax的大小,当
Figure GSA000000946458000216
时,执行步骤S10,否则执行步骤S11;
步骤S10:根据
Figure GSA00000094645800031
计算td,并令flag=3;
步骤S11:判断flag的值:当flag=3时,执行步骤S12;当flag=2时,执行步骤S20;当flag=1时,执行步骤S24;当flag=0时,代表整个运动轨迹的djerk、jerk、a、v都满足其最大值约束,则此时可以结束。
步骤S12:根据x计算tj,并令flag=0;
步骤S13:根据td和tj计算轨迹的最大速度vmax
步骤S14:比较vmax
Figure GSA00000094645800032
的大小,当
Figure GSA00000094645800033
时,执行步骤S15,否则执行步骤S16;
步骤S15:根据
Figure GSA00000094645800034
计算tj,并令flag=1;
步骤S16:根据td和tj计算轨迹的最大加速度amax
步骤S17:比较amax的大小,当
Figure GSA00000094645800036
时,执行步骤S18,否则执行步骤S19;
步骤S18:根据
Figure GSA00000094645800037
计算tj,并令flag=2;
步骤S19:判断flag的值,当flag=2时,执行步骤S20;当flag=1时,执行步骤S24;当flag=0时,则结束。
步骤S20:根据x计算ta
步骤S21:根据td、tj和ta计算轨迹的最大速度vmax
步骤S22:比较vmax的大小,当
Figure GSA00000094645800039
时,则结束,否则执行步骤S23;
步骤S23:根据计算ta
步骤S24:根据x计算tv
虽然以上方法实现了四阶定值发生器,但是其采用的是一种递推式的实现方法,实现过程较为繁琐,包括一些重复的计算。而且其对于计算过程中得到的每个时间段,不论是中间量还是最终的值都要进行离散化,同时还要修改最大的djerk,如此制约着系统控制性能的提高。
另外,前面所得到的时间段都是连续的时间段,但实际上所得到的轨迹都需要进行离散化来实现。以上四阶定值发生器的实现方法所对应的离散化方法是在设计轨迹的过程中不断的进行离散化,即设计轨迹与离散化是同时进行的。例如图3中所示的步骤S1至S8是一个求取td的过程,其中最多可能要计算四次td,而每次计算完td都要对其进行离散化,使其成为采样周期ts的倍数。离散化公式如下:
t d ‾ ′ = ceil ( t d ‾ t s ) × t s
其中
Figure GSA00000094645800042
由于时间段经过了重新计算,因此可能导致规划的轨迹超出了约束的范围,因此离散化的同时还要对进行重新计算。对于其他的时间段也是相同的方法,因此以上四阶定值发生器的离散化方法计算量很大。
可见,现有四阶定值发生器的实现过程复杂,计算量大,且对于系统控制性能的提高有一定的影响。
发明内容
本发明的目的在于提供一种四阶定值发生器的实现方法,以解决现有四阶定值发生器实现过程复杂,计算量大等技术问题。
为解决以上技术问题,本发明提供一种四阶定值发生器的实现方法,用以产生满足约束的四阶定值轨迹,其包括:定义djerk时间段Td、jerk时间段Tj、加速度时间段Ta、速度时间段Tv,其中jerk为加速度的微分,djerk为jerk的微分,该过程包括:令Td=td,其中td为常djerk时间段;令Tj=td+tj,其中tj为常jerk时间段;令Ta=td+tj+ta,其中ta为常加速度时间段;令Tv=td+tj+ta+tv,其中tv为常速度时间段;根据期望的位移
Figure GSA00000094645800045
以及djerk、jerk、加速度和速度的最大值约束
Figure GSA00000094645800046
求出Td、Tj、Ta、Tv的值;得到满足约束的四阶定值轨迹。
进一步的,所述根据期望的位移
Figure GSA00000094645800047
以及djerk、jerk、加速度和速度的最大值约束
Figure GSA00000094645800048
求出Td、Tj、Ta、Tv的值的过程包括:根据求出td;根据和td求出最大位移xmax、速度vmax和加速度amax;分析轨迹是否超出期望的位移
Figure GSA000000946458000411
以及速度、加速度的最大值约束
Figure GSA000000946458000412
的限制,当超出限制时,利用期望的位移
Figure GSA000000946458000413
或相应的最大值约束调整td;根据调整的后的td,计算Td、Tj、Ta、Tv的值。
进一步的,所述分析轨迹是否超出期望的位移
Figure GSA000000946458000414
以及速度、加速度的最大值约束
Figure GSA000000946458000415
的限制,当超出限制时,利用期望的位移
Figure GSA000000946458000416
或相应的最大值约束调整td的过程包括:比较以上最大位移xmax、速度vmax、加速度amax与期望的位移速度、加速度的最大值约束
Figure GSA000000946458000418
以判断位移、速度、加速度是否超出限制;当加速度、速度、位移都超出限制时,分别根据
Figure GSA00000094645800051
Figure GSA00000094645800052
来求取三个时间段tda、tdv、tdx,取其中的最小值作为td;当仅加速度、速度超出限制时,分别根据
Figure GSA00000094645800053
Figure GSA00000094645800054
来求取两个时间段tda、tdv,取其中的最小值作为td;当仅加速度、位移超出约束时,分别根据
Figure GSA00000094645800055
Figure GSA00000094645800056
来求取两个时间段tda、tdv,取其中的最小值作为td;当仅速度、位移超出约束时,分别根据
Figure GSA00000094645800057
Figure GSA00000094645800058
来求取两个时间段tdv、tdx,取其中的最小值作为td;当仅加速度超出约束时,根据
Figure GSA00000094645800059
求取一个时间段tda,且令td=tda;当仅速度超出约束时,根据来求取一个时间段tdv,且令td=tdv;当仅位移超出约束时,根据
Figure GSA000000946458000511
来求取一个时间段tdx,令td=tdx
进一步的,所述常jerk时间段tj的求取方法包括:分别求出对应于加速度的最大值约束
Figure GSA000000946458000512
期望的位移
Figure GSA000000946458000513
以及速度的最大值约束
Figure GSA000000946458000514
的三个时间段tja、tjx和tjv,取其中最小时间段作为tj
进一步的,所述常加速度时间段ta的求取方法包括:分别求出对应于期望的位移以及速度的最大值约束
Figure GSA000000946458000516
的两个时间段tax和tav,取其中最小时间段作为ta
进一步的,所述常速度时间段tv的求取方法包括:求出对应于期望的位移
Figure GSA000000946458000517
的两个时间段tvx,作为常速度时间段tv
进一步的,所述四阶定值发生器的实现方法,还包括:在已经得到各个时间段以后,对每个时间段进行离散化;利用离散化后的时间段,得到离散化后的轨迹的位移;调整djerk的最大值约束
Figure GSA000000946458000518
来保证离散化后的轨迹的位移等于期望的位移并且保证jerk、加速度和速度等不超过其最大值约束
Figure GSA000000946458000520
可见,在以上方法中,采用了不同于现有技术的时间段定义方式,进而采用非递推计算过程来实现四阶定值发生器,其相对于现有技术,实现过程简洁,计算量小。把以上方法与现有方法用Matlab实现,经实验验证,使用本发明提供的方法可以使求取四阶轨迹的最长时间比现有方法缩短约30%。因此,本发明提供的方法能够在更短的时间内计算出四阶轨迹,有利于系统控制性能的提高。
附图说明
图1为一种标准的四阶轨迹示意图;
图2为现有技术的四阶定值发生器实现方法中时间段定义方式的示意图;
图3为现有技术的四阶定值发生器实现方法的流程示意图;
图4为本发明一实施例所提供的四阶定值发生器实现方法的流程示意图;
图5为本发明一实施例所提供的时间段定义方式的示意图;
图6为本发明另一实施例所提供的时间段定义方式的示意图;
图7为本发明一实施例所提供的各个时间段的求取流程示意图;
图8为本发明一实施例所提供四阶定值发生器实现方法的流程示意图;
图9为本发明一实施例所提供常加速度求取方法的流程示意图。
具体实施方式
为让本发明的上述特征和优点能更明显易懂,下文特举示例性实施例,并配合附图,作详细说明如下。
在背景技术中已经提及,四阶定值发生器的目的是要根据期望的位移
Figure GSA00000094645800061
得到时间最短且满足jerk(加速度的微分)的微分(djerk)、jerk、加速度a、速度v等最大值约束的运动轨迹。设
Figure GSA00000094645800062
Figure GSA00000094645800063
分别为djerk、jerk、加速度和速度的最大值约束。
由于现有方法中采取了递推式的实现方法,使得实现过程较为繁琐,包括了一些重复的计算。针对这个问题,本发明的发明人经过研究与实验,提出了一个非递推式的实现方法,其流程更加简单,计算量更小。下面结合图4对该方法概述如下:
S100:首先,定义时间段,在此步骤中,采用了与现有技术不同的定义方式。具体,请参考图5与图6两个实施例中对于时间段的定义:定义djerk时间段Td为常djerk时间段td,即令Td=td;定义jerk时间段Tj为常djerk时间段td加上常jerk时间段tj,即令Tj=td+tj。以此类推,定义加速度时间段Ta为常djerk时间段td、常jerk时间段tj与常加速度时间段ta的和,即令Ta=td+tj+ta;定义速度时间段Tv为常djerk时间段td、常jerk时间段tj、常加速度时间段ta与常速度时间段tv的和,即令Tv=td+tj+ta+tv
S200:而后,根据期望的位移
Figure GSA00000094645800064
以及djerk、jerk、加速度和速度的最大值约束求出Td、Tj、Ta、Tv的值。
S300:进而,得到满足约束的四阶定值轨迹。
由于采用了步骤S100所述的时间段定义方式,步骤S200便可以实现一种非递推式的计算过程,具体如图7所示,包括如下步骤:
S201:根据
Figure GSA00000094645800071
求出td
S202:根据
Figure GSA00000094645800072
和td求出最大位移xmax、速度vmax和加速度amax
S203:分析轨迹是否超出期望的位移
Figure GSA00000094645800073
以及速度、加速度的最大值约束
Figure GSA00000094645800074
的限制,当超出限制时,利用期望的位移
Figure GSA00000094645800075
或相应的最大值约束调整td
S204:根据调整的后的td,计算Td、Tj、Ta、Tv的值。
为使以上步骤更为明显易懂,下面以图6的实施例为例,并结合图8,引入以下公式,来详细描述以上各个步骤:
步骤S201:根据
Figure GSA00000094645800076
求出td,相应公式如下:
t d = j ‾ / d ‾ - - - ( 1 )
此时,Td=td,Tj=td,Ta=2td,Tv=4td,jerk的轨迹如图6所示。
步骤S202:根据
Figure GSA00000094645800078
和td求出最大位移xmax、速度vmax和加速度amax,相应公式如下:
v max = 2 d ‾ t d 3 - - - ( 2 )
a max = d ‾ t d 2 - - - ( 3 )
x max = 8 d ‾ t d 4 - - - ( 4 )
由于图6中的轨迹并没有考虑加速度、速度和位移的约束,因此下面需要分析该轨迹是否满足约束。可以分多种情况:当加速度、速度、位移都超出限制时,分别根据
Figure GSA000000946458000713
来求取三个时间段tda、tdv、tdx,取其中的最小值作为td;当仅加速度、速度超出限制时,分别根据
Figure GSA000000946458000714
来求取两个时间段tda、tdv取其中的最小值作为td;当仅加速度、位移超出约束时,分别根据
Figure GSA000000946458000715
Figure GSA000000946458000716
来求取两个时间段tda、tdx,取其中的最小值作为td;当仅速度、位移超出约束时,分别根据
Figure GSA000000946458000717
来求取两个时间段tdv、tdx,取其中的最小值作为td;当仅加速度超出约束时,根据
Figure GSA000000946458000719
求取一个时间段tda,且令td=tda;当仅速度超出约束时,根据
Figure GSA000000946458000720
来求取一个时间段tdv,且令td=tdv;当仅位移超出约束时,根据
Figure GSA000000946458000721
来求取一个时间段tdx,令td=tdx。对以上各种情况详细描述如下:
S203:分析轨迹是否超出期望的位移
Figure GSA00000094645800081
以及速度、加速度的最大值约束
Figure GSA00000094645800082
Figure GSA00000094645800083
的限制,当超出限制时,利用期望的位移
Figure GSA00000094645800084
或相应的最大值约束调整td
1. a max > a ‾
说明该轨迹的最大加速度已经超出了限制,那么需要对td进行修改。以加速度的最大值约束
Figure GSA00000094645800086
为标准来调整时间段td,公式如下:
t da = a ‾ d ‾ - - - ( 5 )
下面再看位移和速度是否满足要求:
1.1 x max > x ‾ , v max > v ‾
说明位移、速度和加速度都超出了限制。再根据
Figure GSA000000946458000811
调整时间段td,公式如下:
t dv = v ‾ 2 d ‾ 3 - - - ( 6 )
t dx = x ‾ 8 d ‾ 4 - - - ( 7 )
比较tda、tdv和tdx三个量,哪个最小,就以哪个值作为最终的td,这样才能满足所有的要求:
如果tda最小,td=tda,此时加速度已经达到最大,但位移还不满足要求,因此后续需要计算常加速度段ta,以计算得到Ta
如果tdv最小,td=tdv,此时速度已经达到最大,但位移还不满足要求,因此后续需要计算常速度段tv,以计算得到Tv
如果tdx最小,td=tdx,满足所有的约束,而且位移也满足要求,因此结束。
1.2 x max < x &OverBar; , v max > v &OverBar;
说明速度和加速度超出了限制。根据按式(6)计算tdv,然后判断tda和tdv大小:
如果tda<tdv,则td=tda,此时加速度已经达到最大,但位移还不满足要求,因此后续需要计算常加速度段ta,以计算得到Ta
如果tda>tdv,则td=tdv,此时速度已经达到最大,但位移还不满足要求,因此后续需要计算常速度段tv,以计算得到Tv
1.3 x max > x &OverBar; , v max < v &OverBar;
说明位移和加速度都超出限制。根据
Figure GSA00000094645800093
按式(7)计算tdx,然后判断tda和tdx两个时间段的大小:
如果tda<tdx,则td=tda,此时加速度已经达到最大,但位移还不满足要求,因此后续需要计算常加速度段ta,以计算得到Ta
如果tda>tdx,则td=tdx,满足所有约束,且位移也达到要求,因此结束。
1.4 x max < x &OverBar; , v max < v &OverBar;
说明位移和速度都没有超出限制,但是加速度超出限制。因此设td=tda。此时加速度已经达到最大,但位移还不满足要求,因此后续需要计算常加速度段ta,以计算得到Ta
2. a max &le; a &OverBar;
2.1 x max < x &OverBar; , v max > v &OverBar;
说明速度超出了限制。根据按式(6)来求tdv,此时速度已经达到最大,但位移还不满足要求,因此后续需要计算常速度段tv,以计算得到Tv
2.2 x max > x &OverBar; , v max > v &OverBar;
说明位移超出了限制。首先按式(6)根据来求tdv,然后按式(7)根据位移来求tdx,然后判断tdv和tdx两个时间段的大小:
如果tdx>tdv,说明到达最大速度的时间小于到达位移的时间,因此设td=tdv,此时速度已经达到最大,但位移还不满足要求,因此后续需要计算常速度段tv,以计算得到Tv
如果tdx<tdv,说明以位移来计算tdx并没有超出速度的限制,因此td=tdx,结束。
2.3 x max > x &OverBar; , v max < v &OverBar;
说明位移超出限制,速度和加速度没有超出限制。因此轨迹形状不变,只需要按式(7)根据位移重新计算td即可,结束。
2.4 x max < x &OverBar; , v max < v &OverBar;
说明所有的参数都没有超出限制。那么后续需要计算常jerk段tj,以计算得到Tj
下面详细描述常jerk时间段tj、常加速度时间段ta和常速度时间段tv的求取方法。
常jerk时间段tj的求取方法包括:分别求出对应于加速度的最大值约束
Figure GSA00000094645800101
期望的位移
Figure GSA00000094645800102
以及速度的最大值约束
Figure GSA00000094645800103
的三个时间段tja、tjx和tjv,取其中最小时间段作为tj。具体如图8与以下公式:
首先,根据最大加速度
Figure GSA00000094645800104
来求tja为:
t ja = a &OverBar; j &OverBar; - - - ( 8 )
然后根据来求tjx
p = - 1 9 t d 2 q = - 1 27 t d 3 - x &OverBar; 4 d &OverBar; t d D = p 3 + q 2 r = - q + D 3 - - - ( 9 )
根据上式可得tjx为:
t jx = r - p r - 5 3 t d + t d - - - ( 10 )
根据
Figure GSA00000094645800109
来求tjv
t jv 2 + t d &OverBar; t jv - v &OverBar; j &OverBar; = 0 - - - ( 11 )
比较三个值的大小,可得:
如果tja最小,设此时加速度已经达到最大,但位移还不满足要求,因此后续需要计算常加速度段ta,以计算得到Ta
如果tjv最小,设tj=tjv,此时速度已经达到最大,但位移还不满足要求,因此后续需要计算常速度段tv,以计算得到Tv
如果tjx最小,则tj=tjx,结束。
常加速度时间段ta的求取方法如图9所示,包括:分别求出对应于期望的位移
Figure GSA00000094645800111
以及速度的最大值约束
Figure GSA00000094645800112
的两个时间段tax和tav,取其中最小时间段作为ta,具体如下:
首先,根据
Figure GSA00000094645800113
求出tav为:
t av = v &OverBar; a &OverBar; - - - ( 12 )
然后根据
Figure GSA00000094645800115
求出tax
t ax 2 + ( t d + t j ) t ax - x &OverBar; a &OverBar; = 0 - - - ( 13 )
而后比较tax和tav的大小:
tax>tav,说明虽然已经达到最大速度,但是位移仍然不满足要求,因此ta=tav,因此后续需要计算常速度段tv,以计算得到Tv
tax<tav,说明满足位移要求的同时没有超出速度的限制,设ta=tax,结束。
常速度时间段tv的求取方法包括:求出对应于期望的位移
Figure GSA00000094645800117
的两个时间段tvx,作为常速度时间段tv,具体如以下公式:
根据位移计算tv
t v = x &OverBar; v &OverBar; - - - ( 14 )
下面比较一下以上实施例所提供的方法与现有技术在设计四阶轨迹时的效率。用[td,tj,ta,tv]来表示根据期望的位移和约束条件,按照四阶轨迹规划方法所求出的常djerk时间段、常jerk时间段、常加速度时间段和常速度时间段。由于期望的位移及约束条件的不同,常jerk时间段、常加速度时间段和常速度时间段都可能不存在,即对应的时间为0。例如,[td,0,ta,0]就表示所得到的轨迹有常djerk时间段和常加速度时间段,而没有常jerk时间段和常速度时间段的情况。但是,当4个参数都需要,即td、tj、ta和tv都不为0的时候,所需要求的参数最多,即计算量最大,是整个方法的最坏情况。针对最坏情况,比较现有技术的方法和以上实施例提出的方法,要求出所有的参数,现有技术需要经过10次两两比较与16次计算,而以上实施例提供的方法仅需要经过6次两两比较与10次的计算。把两种方法用Matlab实现,经实验验证,使用以上实施例提供的方法可以使求取四阶轨迹的最长时间比现有技术的方法缩短约30%。因此,以上实施例提供的方法能够在更短的时间内计算出四阶轨迹,有利于系统控制性能的提高。
前面所得到的时间段都是连续的时间段,但实际上所得到的轨迹都需要进行离散化来实现。
在背景技术中已经提及现有技术的离散化方法是在设计轨迹的过程中不断的进行离散化,即设计轨迹与离散化是同时进行的。例如现有技术中求td的流程如图3所示,整个流程中最多可能要计算四次td,而每次计算完td都要对其进行离散化,使其成为采样周期ts的倍数:
t d &OverBar; &prime; = ceil ( t d &OverBar; t s ) &times; t s - - - ( 15 )
其中由于时间段经过了重新计算,因此可能导致规划的轨迹超出了约束的范围,因此离散化的同时还要对
Figure GSA00000094645800124
进行重新计算。对于其他的时间段也是相同的方法,因此现有的离散化方法计算量很大。
本发明以上实施例所提出的方法也是连续时间段的计算方法,为此本发明同时也提出了对应的离散化方法。在已经得到各个时间段以后,再对每个时间段进行离散化,计算量更小。具体如下:
首先在连续情况下可知:
t d = j &OverBar; d &OverBar; , t j = a &OverBar; j &OverBar; , t a = v &OverBar; a &OverBar; , t v = x v &OverBar; - - - ( 16 )
通过上式可得:
j &OverBar; = t d d &OverBar; - - - ( 17 )
a &OverBar; = t d t j d &OverBar; - - - ( 18 )
v &OverBar; = t d t j t a d &OverBar; - - - ( 19 )
x = t d t j t a t v d &OverBar; - - - ( 20 )
在离散的情况下,要对以上的时间段都按式(15)的方法进行离散化,使得每个时间段成为采样周期ts的整数倍。设离散化以后的轨迹的位移是:
x &prime; = t d &prime; t j &prime; t a &prime; t v &prime; d &OverBar; - - - ( 21 )
其中td′、tj′、ta′和tv′是经过离散化的时间段,由于离散化的时间段大于连续情况下的时间段,因此离散化后的位移x′肯定大于期望的位移x。针对这个问题,可以通过调整
Figure GSA00000094645800133
的大小来保证位移等于期望值,并且保证jerk、加速度和速度等不超过其最大值约束。
根据期望的位移和实际的位移可得:
f = x x &prime; = t d t j t a t v t d &prime; t j &prime; t a &prime; t v &prime; - - - ( 22 )
Figure GSA00000094645800135
那么利用调整后的
Figure GSA00000094645800136
来计算的轨迹的位移如下:
x &prime; = t d &prime; t j &prime; t a &prime; t v &prime; d &OverBar; &prime; = t d &prime; t j &prime; t a &prime; t v &prime; t d t j t a t v t d &prime; t j &prime; t a &prime; t v &prime; d &OverBar; = x - - - ( 23 )
轨迹的最大的jerk为:
j ^ max = t d &prime; d &OverBar; &prime; = t d &prime; t d t j t a t v t d &prime; t j &prime; t a &prime; t v &prime; d &OverBar; = j &OverBar; t j t a t v t j &prime; t a &prime; t v &prime; - - - ( 24 )
由于离散化后的时间都大于连续情况时的时间,因此
Figure GSA00000094645800139
同理可得
a ^ max < a &OverBar; , v ^ max < v &OverBar; .
这样的话,只需要在连续的情况下计算好轨迹的时间段以后,通过计算一个系数f,就可以实现轨迹的离散化处理,与现有技术提出的方法相比要简单得多。
综上所述,在以上方法中,采用了不同于现有技术的时间段定义方式,进而采用非递推计算过程来实现四阶定值发生器,其相对于现有技术,实现过程简洁,计算量小。把以上方法与现有方法用Matlab实现,经实验验证,使用本发明提供的方法可以使求取四阶轨迹的最长时间比现有方法缩短约30%。因此,本发明提供的方法能够在更短的时间内计算出四阶轨迹,有利于系统控制性能的提高。
以上仅为举例,并非用以限定本发明,本发明的保护范围应当以权利要求书所涵盖的范围为准。

Claims (7)

1.一种四阶定值发生器的实现方法,用以产生满足约束的四阶定值轨迹,其特征是,包括:
定义djerk时间段Td、jerk时间段Tj、加速度时间段Ta、速度时间段Tv,其中jerk为加速度的微分,djerk为jerk的微分,该过程包括:
令Td=td,其中td为常djerk时间段;
令Tj=td+tj,其中tj为常jerk时间段;
令Ta=td+tj+ta,其中ta为常加速度时间段;
令Tv=td+tj+ta+tv,其中tv为常速度时间段;
根据期望的位移
Figure FSA00000094645700011
以及djerk、jerk、加速度和速度的最大值约束
Figure FSA00000094645700012
Figure FSA00000094645700013
求出Td、Tj、Ta、Tv的值;
得到满足约束的四阶定值轨迹。
2.根据权利要求1所述的四阶定值发生器的实现方法,其特征是,所述根据期望的位移
Figure FSA00000094645700014
以及djerk、jerk、加速度和速度的最大值约束
Figure FSA00000094645700015
求出Td、Tj、Ta、Tv的值的过程包括:
根据
Figure FSA00000094645700016
求出td
根据
Figure FSA00000094645700017
和td求出最大位移xmax、速度vmax和加速度amax
分析轨迹是否超出期望的位移
Figure FSA00000094645700018
以及速度、加速度的最大值约束
Figure FSA00000094645700019
的限制,当超出限制时,利用期望的位移
Figure FSA000000946457000110
或相应的最大值约束调整td
根据调整的后的td,计算Td、Tj、Ta、Tv的值。
3.根据权利要求2所述的四阶定值发生器的实现方法,其特征是,所述分析轨迹是否超出期望的位移
Figure FSA000000946457000111
以及速度、加速度的最大值约束
Figure FSA000000946457000112
的限制,当超出限制时,利用期望的位移
Figure FSA000000946457000113
或相应的最大值约束调整td的过程包括:
比较以上最大位移xmax、速度vmax、加速度amax与期望的位移
Figure FSA000000946457000114
速度、加速度的最大值约束
Figure FSA000000946457000115
以判断位移、速度、加速度是否超出限制;
当加速度、速度、位移都超出限制时,分别根据
Figure FSA000000946457000116
Figure FSA000000946457000117
来求取三个时间段tda、tdv、tdx,取其中的最小值作为td
当仅加速度、速度超出限制时,分别根据
Figure FSA000000946457000118
来求取两个时间段tda、tdv,取其中的最小值作为td
当仅加速度、位移超出约束时,分别根据
Figure FSA00000094645700021
Figure FSA00000094645700022
来求取两个时间段tda、tdv,取其中的最小值作为td
当仅速度、位移超出约束时,分别根据
Figure FSA00000094645700023
Figure FSA00000094645700024
来求取两个时间段tdv、tdx,取其中的最小值作为td
当仅加速度超出约束时,根据
Figure FSA00000094645700025
求取一个时间段tda,且令td=tda
当仅速度超出约束时,根据
Figure FSA00000094645700026
来求取一个时间段tdv,且令td=tdv
当仅位移超出约束时,根据
Figure FSA00000094645700027
来求取一个时间段tdx,令td=tdx
4.根据权利要求1所述的四阶定值发生器的实现方法,其特征是,所述常jerk时间段tj的求取方法包括:
分别求出对应于加速度的最大值约束期望的位移
Figure FSA00000094645700029
以及速度的最大值约束
Figure FSA000000946457000210
的三个时间段tja、tjx和tjv,取其中最小时间段作为tj
5.根据权利要求1所述的四阶定值发生器的实现方法,其特征是,所述常加速度时间段ta的求取方法包括:
分别求出对应于期望的位移
Figure FSA000000946457000211
以及速度的最大值约束
Figure FSA000000946457000212
的两个时间段tax和tav,取其中最小时间段作为ta
6.根据权利要求1所述的四阶定值发生器的实现方法,其特征是,所述常速度时间段tv的求取方法包括:
求出对应于期望的位移
Figure FSA000000946457000213
的两个时间段tvx,作为常速度时间段tv
7.根据权利要求1所述的四阶定值发生器的实现方法,其特征是,还包括:
在已经得到各个时间段以后,对每个时间段进行离散化;
利用离散化后的时间段,得到离散化后的轨迹的位移;
调整djerk的最大值约束
Figure FSA000000946457000214
来保证离散化后的轨迹的位移等于期望的位移
Figure FSA000000946457000215
,并且保证jerk、加速度和速度等不超过其最大值约束
Figure FSA000000946457000216
CN2010101536874A 2010-04-22 2010-04-22 四阶定值发生器的实现方法 Active CN102236338B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101536874A CN102236338B (zh) 2010-04-22 2010-04-22 四阶定值发生器的实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101536874A CN102236338B (zh) 2010-04-22 2010-04-22 四阶定值发生器的实现方法

Publications (2)

Publication Number Publication Date
CN102236338A true CN102236338A (zh) 2011-11-09
CN102236338B CN102236338B (zh) 2013-11-13

Family

ID=44887078

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101536874A Active CN102236338B (zh) 2010-04-22 2010-04-22 四阶定值发生器的实现方法

Country Status (1)

Country Link
CN (1) CN102236338B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106898034A (zh) * 2017-02-28 2017-06-27 北京润科通用技术有限公司 一种目标运动轨迹模拟方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5004968A (en) * 1987-12-21 1991-04-02 Fanuc Ltd. Method for acceleration and deceleration control of servomotors
JP2001092533A (ja) * 1999-09-21 2001-04-06 Yaskawa Electric Corp 速度データの動的最適化方法
JP2007272597A (ja) * 2006-03-31 2007-10-18 Yaskawa Electric Corp ロボットの加減速パターン生成方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5004968A (en) * 1987-12-21 1991-04-02 Fanuc Ltd. Method for acceleration and deceleration control of servomotors
JP2001092533A (ja) * 1999-09-21 2001-04-06 Yaskawa Electric Corp 速度データの動的最適化方法
JP2007272597A (ja) * 2006-03-31 2007-10-18 Yaskawa Electric Corp ロボットの加減速パターン生成方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106898034A (zh) * 2017-02-28 2017-06-27 北京润科通用技术有限公司 一种目标运动轨迹模拟方法和装置
CN106898034B (zh) * 2017-02-28 2020-04-03 北京润科通用技术有限公司 一种目标运动轨迹模拟方法和装置

Also Published As

Publication number Publication date
CN102236338B (zh) 2013-11-13

Similar Documents

Publication Publication Date Title
US7919940B2 (en) System and method for jerk limited trajectory planning for a path planner
Ngoduy Instability of cooperative adaptive cruise control traffic flow: A macroscopic approach
WO2017073184A1 (ja) 車速制御装置
CN102207729B (zh) 基于理想误差动态的半周期重复控制器
CN103076759A (zh) 一种获取数控系统的s型曲线的方法及其装置
CN104036331A (zh) 一种基于改进粒子群算法的电力系统动态经济调度方法
CN111341104B (zh) 一种交通流跟驰模型的速度时滞反馈控制方法
Corcau et al. A fuzzy logic controller as a power system stabilizer
De Wit Best-effort highway traffic congestion control via variable speed limits
Azmi et al. Optimization of the PID-PD parameters of the overhead crane control system by using PSO algorithm
CN114077226A (zh) S型曲线速度规划方法、控制终端及计算机可读存储介质
Sakhdari et al. Robust tube-based mpc for automotive adaptive cruise control design
Ji et al. Bifurcation control of solid angle car-following model through a time-delay feedback method
CN102236338B (zh) 四阶定值发生器的实现方法
JP2019077346A (ja) 車間距離制御装置
Dahmane et al. Stochastic mpc for optimal energy management strategy of hybrid vehicle performing acc with stop&go maneuvers
Wang et al. Driver support and cooperative systems control design: Framework and preliminary results
CN111130104B (zh) 一种电力系统混沌振荡的无源性滑模控制方法
Ichchou et al. A reduced energy supply strategy in active vibration control
CN111798060A (zh) 一种基于机组爬坡率估计的功率指令优化分配方法
Liu et al. A Centralized Cooperative Lane-Changing Strategy for Merging into CAV Platoon
CN113968263B (zh) 一种列车运行曲线优化方法、装置及电子设备
Zhang et al. Feasibility Analysis and Proof of Model Predictive Control
Lu et al. A Comprehensive Comparison of Three Desired Acceleration Calculation Methods for Mass Production of Adaptive Cruise System
CN113805485B (zh) 一种暖启动c/gmres方法、系统、设备及介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 1525 Zhang Dong Road, Shanghai, No. 201203

Patentee after: Shanghai microelectronics equipment (Group) Limited by Share Ltd

Address before: 1525 Zhang Dong Road, Shanghai, No. 201203

Patentee before: Shanghai Micro Electronics Equipment Co., Ltd.