CN107980108A - 机器人运动轨迹规划方法及相关装置 - Google Patents

机器人运动轨迹规划方法及相关装置 Download PDF

Info

Publication number
CN107980108A
CN107980108A CN201780002229.0A CN201780002229A CN107980108A CN 107980108 A CN107980108 A CN 107980108A CN 201780002229 A CN201780002229 A CN 201780002229A CN 107980108 A CN107980108 A CN 107980108A
Authority
CN
China
Prior art keywords
point
movement locus
section
sections
bezier curve
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
CN201780002229.0A
Other languages
English (en)
Other versions
CN107980108B (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.)
Shenzhen A&E Intelligent Technology Institute Co Ltd
Original Assignee
Shenzhen A&E Intelligent Technology Institute 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 Shenzhen A&E Intelligent Technology Institute Co Ltd filed Critical Shenzhen A&E Intelligent Technology Institute Co Ltd
Publication of CN107980108A publication Critical patent/CN107980108A/zh
Application granted granted Critical
Publication of CN107980108B publication Critical patent/CN107980108B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02—Control of position or course in two dimensions
    • G05D1/021—Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0212—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02—Control of position or course in two dimensions

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Manipulator (AREA)

Abstract

本方案涉及运动控制领域，尤其涉及机器人运动轨迹规划技术。在一种机器人运动轨迹规划方法中，根据第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点，n不小于4(401)；根据上述特征点构建n次Bezier曲线，作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段(402)。由于采用n次Bezier曲线作为两CP运动轨迹段之间的平滑过渡轨迹段，该平滑过渡轨迹段在其起始点与第一CP运动轨迹段相切且曲率相同，同时在其终止点与第二CP运动轨迹段相切且曲率相同，从而可实现连续CP运动的平滑转接，同时，速度不必降为0，并且，加速度不跳变。

Description

机器人运动轨迹规划方法及相关装置
技术领域
本发明涉及运动控制领域，更具体地说，涉及机器人运动轨迹规划技术。
背景技术
在对机器人运动控制中，被控部件(例如操作臂)的轨迹运动方式可分为CP运动(连续运行，包括直线运动和圆弧运动)和PTP运动(点到点运动)。
在对机器人进行运动轨迹规划时，可能需要多段CP运动连续进行。以两段CP运动连续进行为例，请参见图1a，CP运动轨迹段1(也可称为轨迹段1)和CP运动轨迹段2(也可称为轨迹段2)可分别为直线段，此外，还可圆弧段与圆弧段(图1b)、圆弧段与直线段(图1c和图1d)连续进行。
在单段CP运动的规划中，其初速度和末速度都为0，如果在多段CP连续运动的过程中令每段CP运动的初始速度和末速度都降为0，会大大降低工作效率，并且频繁的加减速也会对电机和减速器的寿命有影响。因此，运动轨迹规则的趋势是不令多段CP连续运动中的各段CP运动的速度降为零。
这就面临另一问题：仍请参见图1a，若连续的直线段不在同一直线上，就会有尖角。尖角的顶点为O(O点是两轨迹段的交点)。同理，请参见图1b，若连续的两圆弧段切线不重合，也会有尖角。同理，若直线段与的圆弧段的切线不重合(请参见图1c和图1d)，也会有尖角。速度不为0通过尖角的话，会产生振动。
为此，可在连续的两CP运动轨迹段之间设计平滑过渡轨迹段。例如，请参见图2，可采用矢量叠加法得到平滑过渡轨迹段。这种方式通过空间向量加法，令图2中轨迹段1从拐出点A到交点O的插补点，与轨迹段2从拐入点B到交点O的插补点重合。叠加后的插补点会组成空间的一段弧形(也即图2中点A与点B之间的弧形)，使轨迹平滑。
也即，最终规划出的运行轨迹是，由C点运行至A点，之后拐出轨迹段1，进入点A和点B之间的平滑过渡轨迹段，再由点B拐入轨迹段2，最后到达D。
当然，对于超过两段CP运动的连续运动，可分别对两两连续的CP运动轨迹段进行如此设计。
矢量叠加法虽然对前后相连的CP运动进行了平滑，但在实际应用中，如图1a-1d所示，有直线段与直线段转接、直线段与圆弧段转接、圆弧段与圆弧段转接。矢量叠加法虽可用圆弧替代尖角，但其仅在直线段与直线段转接的场景下肯定适用，在直线段与圆弧段转接、圆弧段与圆弧段的转接中可能无法找到与两条轨迹都相切的圆弧。这是因为空间中相交的直线段与直线段是肯定共面的，但是直线段与圆弧段、圆弧段与圆弧段都不一定共面，在不共面的情况下就无法找出公切的圆弧转接。
因此，现在需要机器人运动轨迹规划技术解决方案，以令无论连续的运动轨迹段是否共面均能进行平滑转接，并在拐点处(拐入点和拐出点)速度连续(也即速度不必降为0)。
发明内容
有鉴于此，本发明实施例的目的在于提供机器人运动轨迹规划方法及相关装置，以令连续的CP运动平滑转接，同时，速度不必降为0。
本申请发明人在研究过程中发现：
Bezier曲线由线段和节点组成的数学曲线。其整体形状并不局限于圆弧，可实现非共面的轨迹平滑，从而可应用在机器人运动轨迹规划方面。
可令Bezier曲线作为连续的运动轨迹段间的平滑过渡轨迹段，Bezier曲线的起始点与一运动轨迹段(可称为lc轨迹段)的拐出点重合，Bezier曲线的终止点与另一运动轨迹段(可称为ln轨迹段)的拐入点重合。
发明人发现，欲实现在拐点处(拐出点和拐入点)速度连续，则Bezier曲线在其起始点需与lc轨迹相切，同时在其终止点需与ln轨迹相切。
n次(n次)Bezier曲线公式为：
其中bi,n(t)称为Bernstein基函数，t＝0对应Beizer曲线的起始点，t＝1对应Beizer曲线的终止点。
P0-Pn这n个特征点或控制点在平面或在三维空间定义了n次Bezier曲线，其组成的多边形称为特征多边形或控制多边形。其中，P0为起始点，Pn为终止点，P1-Pn-1为中间点。例如请参见图3，P0、P1、P2、P3四个点定义了三次Bezier曲线，三次Bezier曲线起始于P0走向P1，并从P2方向到P3。需要说明的是，一般三次Bezier曲线不经过P1、P2，这两个点只提供方向。
根据Bezier基函数的导函数性质可求取n次Bezier曲线始末点的切矢量为则：
当t＝0时，B′(0)＝n(P1-P0)，B′(0)即为n次Bezier曲线在起始点的切矢量；
当t＝1时，B′(1)＝n(Pn-Pn-1)，B′(1)即为n次Bezier曲线在终止点的切矢量。
观察上述切矢量可知，n次Bezier曲线在起始点和终止点处的切线方向(即切矢量)和特征多边型的第一条边及最后一条边走向一致。以图3所示三次Bezier曲线为例，B′(0)＝3(P1-P0)，B′(1)＝3(P3-P2)，而图3的第一条边是P1P0，最后一条边是P3P2。
为令n次Bezier曲线在始末点的切矢方向分别与lc轨迹段和ln轨迹段切矢方向相同，需要最少四个特征点或控制点，所以可选用三次Bezier曲线作为平滑过渡轨迹段。
然而，本申请发明人在进一步的研究过程中发现：
使用三次Bezier曲线作为平滑过渡轨迹段虽然能令连续的CP运动平滑转接，同时，速度不必降为0，但仅能保证在拐点处速度连续，而在拐点处的加速度则不一定连续。随着机器人科技的发展，对于机器人轨迹规划的要求越来越高，如果加速度不连续，在拐点处跳变的加速度会造成电机震动，对电机和减速器造成损伤。
本申请发明人通过分析加速度跳变的原因，发现加速度跳变的原因是因为突增的向心加速度。
加速度连续的充要条件为在速度连续的基础上曲率连续。
某一点的向心加速度有公式其中v表示该点的速度，ρ表示该点的曲率。在v连续的情况下，保证曲率连续就能保证轨迹加速度连续。
曲线的曲率公式有：其中B′(t)表示Bezier曲线的一阶导函数，B″(t)表示Bezier曲线的二阶导函数。通过曲率公式可以得出：确保Bezier曲线二阶可导即可确保Bezier曲线上任一点加速度连续。所以只要次数大于1的Bezier曲线，其轨迹一定曲率连续。
另一方面，为保证拐点处不跳变，要求拐点处速度和加速度均连续，这需要Bezier曲线与lc轨迹段在起始点(即拐出点)相切且曲率连续，同时，Bezier曲线与ln轨迹段在终止点(即拐入点)相切且曲率连续。
根据Bezier基函数的导函数性质可求取n次Bezier曲线始末点的二阶导矢：
求得：当t＝0时，B″(0)＝n(n-1)(P2-2P1+P0)
当t＝1时，B″(1)＝n(n-1)(Pn-2Pn-1+Pn-2)
上述公式表明，二阶导矢只与相邻的3个控制点有关。
为了二阶导矢方向与lc轨迹段和ln轨迹段的二阶导矢方向一致，需要最少五个控制点(或特征点)，所以可选用n次(n不小于4)Bezier曲线，四次Bezier曲线有公式：
B(t)＝P0(1-t)4+4P1t(1-t)3+6P2t2(1-t)2+4P3t3(1-t)+P4t4，t∈[0，1]
基于上述研究发现，本发明实施例提供如下技术方案：
一方面，本申请的实施例提供机器人运动轨迹规划方法，至少用于基于n次(n不小于4)Bezier曲线实现连续的两CP运动轨迹段之间的平滑转接，所述连续的两CP运动轨迹段包括第一CP运动轨迹段和第二CP运动轨迹段；该方法包括：根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点；所述特征点包括P0、P1、P2、P3，其中，P0为起始点、P3为终止点，P1、P2为中间点；根据所述特征点构建n次(n不小于4)Bezier曲线，作为所述第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段，所述平滑过渡轨迹段的起始点为所述第一CP运动轨迹段的拐出点，所述平滑过渡轨迹段的终止点为所述第二CP运动轨迹段的拐入点；其中，所述平滑过渡轨迹段与所述第一CP运动轨迹段在所述起始点上的切矢量方向相同且曲率相同；所述平滑过渡轨迹段与所述第二CP运动轨迹段在所述终止点上的切矢量方向相同且曲率相同。在本实施例中，采用n次(n不小于4)Bezier曲线作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段，由于平滑过渡轨迹段在始末点的切矢方向分别与所述第一CP运动轨迹段和第二CP运动轨迹段的切矢方向相同，因此，该平滑过渡轨迹段在其起始点与第一CP运动轨迹段相切，同时在其终止点与第二CP运动轨迹段相切，从而可实现在拐点处(拐出点和拐入点)速度连续，也即可实现连续的CP运动平滑转接，同时，速度不必降为0。此外，该平滑过渡轨迹段在其起始点的曲率与第一CP运动轨迹段在拐出点的曲率相同(起始点与拐出点重合，或也可认为起始点即拐出点)，在其终止点的曲率与第二CP运动轨迹段在拐入点的曲率相同(终止点与拐入点重合，或也可认为终止点即拐入点)，从而可实现在拐点处(拐出点和拐入点)加速度连续，不跳变。
在一个可能的设计中，在所述确定n次(n不小于4)Bezier曲线的特征点之前，还包括：规划连续的CP运动轨迹段。
在一个可能的设计中，所述n＝4，所述特征点包括P0、P1、P2、P3和P4，其中，P0为起始点、P4为终止点，P1、P2和P3为第一至第三中间点；所述第一CP运动轨迹段和第二CP运动轨迹段的交点表示为O点；所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点包括：确定所述交点O为所述第二中间点P2。
在一个可能的设计中，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点还包括：若所述第一CP运动轨迹段为直线段，选择所述交点O和起始点P0之间的线段上的一点作为所述第一中间点P1；若所述第一CP运动轨迹段为圆弧段，确定第一中间点P1位于起始点P0的切线上，并且，线段P1P0的长度为其中，r1表示所述第一CP运动轨迹段的半径，α1表示所述起始点处的圆心角。
在一个可能的设计中，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点还包括：若所述第二CP运动轨迹段为直线段，选择所述交点O和终止点P4之间的线段上的一点作为所述第三中间点P3；若所述第二CP运动轨迹段为圆弧段，确定所述第三中间点P3位于终止点P4的切线上，并且，线段P2P4的长度为其中，r2表示所述第二CP运动轨迹段的半径，α2表示终止点处的圆心角。
在一个可能的设计中，所述n＝5，所述特征点包括P0、P1、P2、P3、P4和P5，其中，P0为起始点、P5为终止点，P1至P4为第一至第四中间点；所述第一CP运动轨迹段和第二CP运动轨迹段的交点表示为O点。
在一个可能的设计中，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点包括：若所述第一CP运动轨迹段为圆弧段，确定所述第一中间点P1位于所述起始点P0的切线上，并且，线段P1P0的长度为其中，r1表示所述第一CP运动轨迹段的半径，α1表示所述起始点处的圆心角；确定所述第二中间点P2为，以所述起始点P0和交点O为端点的圆弧段上的三等分点。
在一个可能的设计中，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点包括：若所述第一CP运动轨迹段为直线段，选择以所述起始点P0和交点O为端点的直线段上的三等分点分别作为第一中间点P1和第二中间点P2。
在一个可能的设计中，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点还包括：若所述第二CP运动轨迹段为圆弧段，确定所述第三中间点P3位于终止点的切线上，并且，线段P3P5的长度为其中，r2表示所述第二CP运动轨迹段的半径，α2表示所述终止点处的圆心角；确定所述第四中间点P4端点为，以所述终止点和交点O为端点的圆弧段上的三等分点。
在一个可能的设计中，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点还包括：若所述第二CP运动轨迹段为直线段，选择以所述终止点和交点O为端点的直线段上的三等分点分别作为所述第三中间点P3和第四中间点P4。
另一方面，本发明实施例提供了一种机器人运动轨迹规划装置，该装置具有实现上述方法实际中机器人运动轨迹规划装置行为的功能。所述功能可以通过硬件实现，也可以通过硬件执行相应的软件实现。所述硬件或软件包括一个或多个与上述功能相对应的模块。
在一个可能的设计中，机器人运动轨迹规划装置的结构包括：处理器和存储器，所述处理器通过运行存储在所述存储器内的软件程序、调用存储在所述存储器内的数据，执行上述方法。
再一方面，本发明实施例提供了一种机器人，其包括上述的机器人运动轨迹规划装置，以及按照上述机器人运动轨迹规划装置规划出的连续的两CP运动轨迹段及平滑过渡轨迹段运行的被控器件。该控制器件例如可以是操作臂。
再一方面，本发明实施例提供了一种计算机存储介质，用于储存为上述机器人运动轨迹规划装置所用的计算机软件指令，其包含用于执行上述方面所设计的程序。
再一方面，本发明实施例提供了一种计算机存储介质，用于储存为上述机器人所用的计算机软件指令，其包含用于执行上述方面所设计的程序。
相较于现有技术，本实施例采用n次(n不小于4)Bezier曲线作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段，由于平滑过渡轨迹段在始末点的切矢方向分别与所述第一CP运动轨迹段和第二CP运动轨迹段的切矢方向相同，因此，该平滑过渡轨迹段在其起始点与第一CP运动轨迹段相切，同时在其终止点与第二CP运动轨迹段相切，从而可实现在拐点处(拐出点和拐入点)速度连续，也即可实现连续的CP运动平滑转接，同时，速度不必降为0。此外，该平滑过渡轨迹段在其起始点的曲率与第一CP运动轨迹段在拐出点的曲率相同(起始点与拐出点重合，或也可认为起始点即拐出点)，在其终止点的曲率与第二CP运动轨迹段在拐入点的曲率相同(终止点与拐入点重合，或也可认为终止点即拐入点)，从而可实现在拐点处(拐出点和拐入点)加速度连续，不跳变。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1a-1d为本发明实施例提供的连续的运动轨迹段的示意图；
图2为本发明实施例提供的采用矢量叠加法得到平滑过渡轨迹段示意图；
图3为本发明实施例提供的n次(n不小于4)Bezier曲线示意图；
图4-5，6a、7a、8a、9a、10a、11a、12a为本发明实施例提供的机器人运动轨迹规划方法示例性流程图；
图6b、7b、8b、9b、10b、11b、12b为本发明实施例提供的两连续运动轨迹段利用平滑过渡轨迹段过渡示意图；
图6c、7c、9c、10c、11c、12c为本发明实施例提供的曲率图；
图8c为本发明实施例提供的几何原理图；
图13为本发明实施例提供的机器人运动轨迹规划装置结构示意图；
图14为本发明实施例提供的机器人运动轨迹规划装置通用计算机架构示意图。
具体实施方式
本发明可能使用的技术名词、简写或缩写如下：
TCP点：工具坐标中心点；
直线运动：机器人TCP点在笛卡尔空间内沿直线轨迹运动；
圆弧运动：机器人CP点在笛卡尔空间内沿圆弧轨迹运动；
CP：连续运动，笛卡尔空间下的轨迹运动方式，包括直线运动和圆弧运动；
n次Bezier曲线：n次Bezier曲线公式为：称为Bernstein基函数，t＝0对应Beizer曲线的起始点，t＝1对应Beizer曲线的终止点。
P0-Pn这n个特征点或控制点在平面或在三维空间定义了n次Bezier曲线，其组成的多边形称为特征多边形或控制多边形。其中，P0为起始点，Pn为终止点，P1-Pn-1为中间点。例如请参见图3，P0、P1、P2、P3四个点定义了三次Bezier曲线，三次Bezier曲线起始于P0走向P1，并从P2方向到P3。需要说明的是，一般三次Bezier曲线不经过P1、P2，这两个点只提供方向。
本发明实施例提供了机器人运动轨迹规划方法及相关装置，以令连续的CP运动平滑转接，同时，速度不必降为0并且加速度不跳变。本发明实施例所涉及到的机器人运动轨迹规划装置应用于机器人中，其具体可为机器人中的控制器/处理器。
本发明的一个实施例提供一种实现机器人运动轨迹规划的方法，和基于这个方法的机器人运动轨迹规划装置。该装置通过执行上述方法,可实现连续的两CP运动轨迹段(第一CP运动轨迹段和第二CP运动轨迹段)之间的平滑转接，请参见图4，其具体操作包括：
在401部分：根据第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线(n不小于4)的特征点；
以n＝4为例，特征点可包括P0、P1、P2、P3和P4，其中，P0为起始点、P4为终止点，P1、P2、P3为中间点。
而若n＝5，则特征点可包括P0、P1、P2、P3、P4和P5，其中，P0为起始点、P5为终止点，其他点为中间点。
在402部分：根据上述特征点构建上述n次Bezier曲线(n不小于4)，作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段。
在确定出特征点后，根据特征点构建n次Bezier曲线(n不小于4)可采用现有的构建方式，在此不作赘述。
需要注意的是，上述平滑过渡轨迹段的起始点与第一CP运动轨迹段的拐出点相重合(或可认为平滑过渡轨迹段的起始点即为第一CP运动轨迹段的拐出点)，平滑过渡轨迹段的终止点与所述第二CP运动轨迹段的拐入点相重合(或可认为平滑过渡轨迹段的终止点为第二CP运动轨迹段的拐入点)。这样可保证路径连续。
此外，平滑过渡轨迹段在起始点P0上的切矢量方向、与第一CP运动轨迹段在拐出点上的切矢量方向相同，这样，平滑过渡轨迹段在其起始点与第一CP运动轨迹段相切，并且平滑过渡轨迹段在起始点P0上的曲率与第一CP运动轨迹段在拐出点上的曲率相同，这样，可保证在拐出点的速度及加速度均连续。
同理，平滑过渡轨迹段在其终止点上的切矢量方向和第二CP运动轨迹段在拐入点上的切矢量方向的方向相同。这样，平滑过渡轨迹段在其终止点与第二CP运动轨迹段相切。并且平滑过渡轨迹段在终止点上的曲率与第二CP运动轨迹段在拐入点上的曲率相同。这样，可保证在拐入点的速度及加速度均连续。
在一个示例中，拐出点和拐入点可根据用户的配置来确定。具体的，用户可指示拐出点与交点之间的距离(交点是第一CP运动轨迹段和第二CP运动轨迹段的交点)，拐入点与交点之间的距离与拐出点与交点之间的距离相等，这样，就可确定出拐入点和拐出点。
需要说明的是，对于连续的多段CP运动轨迹段，可针对每两个连续的CP运动轨迹段执行上述401部分和402部分。例如，对于连续的三段CP运动轨迹段1-3，可先规划出CP运动轨迹段1-2之间的平滑过渡轨迹段，再规划出CP运动轨迹段2-3之间的平滑过渡轨迹段。
下面结合附图5，对本发明的实施例提供的方案进行说明。
图5为本发明实施例提供的机器人运动轨迹规划方法的另一示例性流程图。
在500部分：规划连续的CP运动轨迹段。
连续的CP运动轨迹段也是需要预先规划的。500部分是后续501部分和502部分的基础。
在一个示例中，可一次性规划出所需的所有CP运动轨迹段。例如，一共需要连续的4个CP运动轨迹段以完成机器人的某个操作，则可一次性规划出4个连续的CP运动轨迹段。
在另一个示例中，也可分批次规划出所需的所有CP运动轨迹段。仍以需要连续的CP运动轨迹段1-4以完成机器人的某个操作为例，可先规划CP运动轨迹段1-2，然后，执行后续的501部分和502部分，以得到CP运动轨迹段1-2之间的平滑过渡轨迹段。之后，再规划出CP运动轨迹段3，再执行后续的501部分和502部分，以得到CP运动轨迹段2-3之间的平滑过渡轨迹段，以此类推。
在又一个示例中，也可采用前瞻X段的方式进行规划，X一般取3。以需要连续的CP运动轨迹段1-8以完成机器人的某个操作为例，可先规划CP运动轨迹段1-2，然后，执行后续的501部分和502部分，以得到CP运动轨迹段1-2之间的平滑过渡轨迹段。之后，再规划出CP运动轨迹段3，再执行后续的501部分和502部分，以得到CP运动轨迹段2-3之间的平滑过渡轨迹段。然后，等CP运动轨迹段1运行完后，再规划出CP运动轨迹段4，然后，执行后续的501部分和502部分，以得到CP运动轨迹段3-4之间的平滑过渡轨迹段。再等CP运动轨迹段2运行完后，规划出CP运动轨迹段5，然后，执行后续的501部分和502部分，以得到CP运动轨迹段4-5之间的平滑过渡轨迹段，以此类推。
在501部分：根据第一CP运动轨迹段和第二CP运动轨迹段确定n次(n不小于4)Bezier曲线的特征点。
501部分与401部分相类似，在此不作赘述。
在502部分：根据上述特征点构建n次(n不小于4)Bezier曲线，作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段。
502部分与402部分相类似，在此不作赘述。
前已述及，n＝4时，有五个特征点，n＝5时，有六个特征点。在n＝4的情况下，中间点的确定，与n＝5时，中间点的确定并不相同。并且，在实际应用中，第一CP运动轨迹段和第二CP运动轨迹段可分别为：直线段与直线段、直线段与圆弧段、圆弧段与直线段、圆弧段与圆弧段四种情况。下面将分情况，对本发明的实施例做进一步说明。
图6a为n＝4且第一CP运动轨迹段和第二CP运动轨迹段分别为直线段时，本发明实施例提供的又一种机器人运动轨迹规划方法示例性流程图。图6b显示的是第一CP运动轨迹段(lc)和第二CP运动轨迹段(ln)以及平滑过渡轨迹段。lc轨迹段和ln轨迹段相交于交点O。
在601部分：确定lc轨迹段的拐出点为平滑过渡轨迹段的起始点P0，确定ln轨迹段的拐入点为平滑过渡轨迹段的终止点P4。
在602部分：确定交点O为(第二)中间点P2。
前已证明，n次(4次)Bezier曲线在起始点的曲率与相邻的3个特征点(P0、P1、P2)相关，在终止点处的曲率也与相邻的3个特征点(P2、P3和P4)有关，则P2点同时需要与P0、P1、P3和P4共面，因此，选择交点O作为P2。
在603部分：选择线段P0O(也即交点O和起始点P0之间的线段)上的一点作为(第一)中间点P1；线段P0O是交点O和起始点P0之间的线段。
在一个示例中，中间点P1到交点O的距离，等于其到起始点P0的距离。因此，在求解P1坐标时，可使P0P1＝P1O，进而获得P1点坐标。P0P1＝P1O可以使特征点分布均匀，特征点均匀分布则可以使Beizer曲线的参数均匀，轨迹曲率变化平稳。
在604部分：选择线段OP4上的一点作为(第三)中间点P3。
线段OP4是交点O和终止点P4之间的线段。
在一个示例中，中间点P3到交点O的距离，等于中间点P3到终止点P4的距离。因此，在求解P3坐标时，可使P3P4＝P3O，从而获得P3点坐标。
在605部分：根据上述特征点构建四次Bezier曲线，作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段(图6b虚线部分即为平滑过渡轨迹段)。
605部分与前述的402和502部分相类似，在此不作赘述。
请参见图6c，通过图6a所示的方法，在始末点曲率皆为0，进而实现在拐点处加速度连续。
图7a为n＝4且第一CP运动轨迹段为直线段、第二CP运动轨迹段为圆弧段时，本发明实施例提供的又一种机器人运动轨迹规划方法示例性流程图。图7b显示的是第一CP运动轨迹段(lc)和第二CP运动轨迹段(ln)以及平滑过渡轨迹段。
在701部分：确定lc轨迹段的拐出点为平滑过渡轨迹段的起始点P0，确定ln轨迹段的拐入点为平滑过渡轨迹段的终止点P4。
701部分与601部分相类似，在此不作赘述。
在702部分：确定交点O为(第二)中间点P2。
702部分与602部分相类似，在此不作赘述。
在703部分：选择线段P0O(也即交点O和起始点P0之间的线段)上的一点作为中间点P1。
703部分与603部分相类似，在此不作赘述。
在704部分：确定中间点P3位于终止点P4的切线上，并且，线段P3P4的长度为(可参见图7b)。
其中，r2表示ln轨迹段的半径，α2表示终止点P4处的圆心角。
公式的推导过程请参见图8a所示实施例。
在705部分：根据上述特征点构建四次Bezier曲线，作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段(图7b虚线部分即为平滑过渡轨迹段)。
705部分与前述的402、502和605部分相类似，在此不作赘述。
请参见图7c，通过图7a所示的方法，在起始点曲率为0，在终止点曲率为半径(50)的倒数0.02，进而实现在拐点处加速度连续。
图8a为n＝4且第一CP运动轨迹段为圆弧段、第二CP运动轨迹段为直线段时，本发明实施例提供的又一种机器人运动轨迹规划方法示例性流程图。图8b显示的是第一CP运动轨迹段(lc)和第二CP运动轨迹段(ln)以及平滑过渡轨迹段。
在801部分，确定lc轨迹段的拐出点为平滑过渡轨迹段的起始点P0，确定ln轨迹段的拐入点为平滑过渡轨迹段的终止点P4。
801部分与601、701部分相类似，在此不作赘述。
在802部分：确定交点O为(第二)中间点P2。
802部分与602、702部分相类似，在此不作赘述。
在803部分：确定中间点P1位于起始点P0的切线上，并且，线段P1P0的长度为
其中，r1表示lc轨迹段的半径，α1表示起始点的圆心角。
下面讲述公式和的由来。
请参见图8c，有两条连续CP运动轨迹段，前一段为圆弧段，后一段为直线段，点A为交点，点C为拐出点，圆弧段半径为r，拐出点处的圆心角为α。根。
B点为P1点或P3点，B点位于圆弧在点C处的切线上，并满足n次Bezier曲线在起始点处(即C点)曲率等于圆弧曲率。已知圆弧段的曲率为半径的倒数，则列等式：
其中：
n为Bezier曲线次数；
∠θ为向量AB与向量BC的夹角。
建立直角三角形ABD，得|AB|·sinθ＝AD＝r-r·cosα。
化简等式两边得到
在804部分：选择线段OP4上的一点作为中间点P3。线段OP4是交点O和终止点之间的线段。
804部分与604部分相类似，在此不作赘述。
在805部分：根据上述特征点构建四次Bezier曲线，作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段(图8b虚线部分即为平滑过渡轨迹段)。
805部分与前述的402、502、605和705部分相类似，在此不作赘述。
图9a为n＝4且第一CP运动轨迹段和第二CP运动轨迹段均为圆弧段时，本发明实施例提供的又一种机器人运动轨迹规划方法示例性流程图。图9b显示的是第一CP运动轨迹段(lc)和第二CP运动轨迹段(ln)以及平滑过渡轨迹段。
在901部分，确定lc轨迹段的拐出点为平滑过渡轨迹段的起始点P0，确定ln轨迹段的拐入点为平滑过渡轨迹段的终止点P4。
901部分与601、701、801部分相类似，在此不作赘述。
在902部分：确定交点O为(第二)中间点P2。
902部分与602、702和802部分相类似，在此不作赘述。
在903部分：确定中间点P1位于起始点P0的切线上，并且，线段P1P0的长度为
903部分与803部分相类似，在此不作赘述。
在904部分：确定中间点P3位于终止点P4的切线上，并且，线段P2P4的长度为(可参见图9b)。
其中，r2表示ln轨迹段的半径，α2表示终止点P4处的圆心角。
904部分与704部分相类似，在此不作赘述。
在905部分：根据上述特征点构建四次Bezier曲线，作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段(图9b虚线部分即为平滑过渡轨迹段)。
905部分与前述的402、502、605、705、805部分相类似，在此不作赘述。
请参见图9c，通过图9a所示的方法，在始末点曲率皆为半径(50)的倒数0.02，进而实现在拐点处加速度连续。
考虑到四次Bezier曲线的五个控制点中，P2点选择两轨迹交点，自由度不够。为了使平滑轨迹可调，增加更多的选择，可以把四次Bezier曲线升阶为五次Bezier曲线。
下面将介绍基于五次Bezier曲线的运动轨迹规划方法。五次Bezier曲线的公式为：
B(t)＝P0(1-t)5+5P1t(1-t)4+10P2t2(1-t)3+10P3t3(1-t)2...+5P4t4(1-t)+P5t5，t∈[0，1]
图10a为n＝5且第一CP运动轨迹段和第二CP运动轨迹段分别为圆弧段时，本发明实施例提供的又一种机器人运动轨迹规划方法示例性流程图。图10b显示的是第一CP运动轨迹段(lc)和第二CP运动轨迹段(ln)以及平滑过渡轨迹段，lc轨迹段和ln轨迹段相交于交点O。特征点包括P0、P1、P2、P3、P4和P5，其中，P0为起始点、P5为终止点，P1至P4为第一至第四中间点。
在1001部分：确定lc轨迹段的拐出点为平滑过渡轨迹段的起始点P0，确定ln轨迹段的拐入点为平滑过渡轨迹段的终止点P5。
1001部分与601、701、801、901部分相类似，在此不作赘述。
在1002部分：确定(第一)中间点P1位于起始点P0的切线上，并且，线段P1P0的长度为
其中，r1表示lc轨迹段的半径，α1表示起始点的圆心角。
1002部分与803部分相类似，在此不作赘述。
在1003部分：确定(第二)中间点P2为，以起始点P0和交点O为端点的圆弧段(该圆弧段属于lc轨迹段的一部分)上的三等分点。
由于三等分点有两个，在一个示例中，可选取靠近交点的三等分点作为中间点P2。
在1004部分：确定第四中间点P4位于终止点P5的切线上，并且，线段P4P5的长度为
其中，r2表示ln轨迹段的半径，α2表示终止点处的圆心角。
1004部分与904部分相类似，在此不作赘述。
在1005部分：确定(第三)中间点P3端点为，以终止点P5和交点O为端点的圆弧段(该圆弧段属于ln轨迹段的一部分)上的三等分点。
由于三等分点有两个，在一个示例中，可选取靠近交点的三等分点作为中间点P3。
在1006部分：根据上述特征点构建五次Bezier曲线，作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段(图10b虚线部分即为平滑过渡轨迹段)。
请参见图10c，通过图10a所示的方法，在始末点曲率皆为半径(50)的倒数0.02，进而实现在拐点处加速度连续。
图11a为n＝5且第一CP运动轨迹段为直线段，第二CP运动轨迹段为圆弧段时，本发明实施例提供的又一种机器人运动轨迹规划方法示例性流程图。图11b显示的是第一CP运动轨迹段(lc)和第二CP运动轨迹段(ln)以及平滑过渡轨迹段，lc轨迹段和ln轨迹段相交于交点O。
在1101部分：确定lc轨迹段的拐出点为平滑过渡轨迹段的起始点P0，确定ln轨迹段的拐入点为平滑过渡轨迹段的终止点P5。
1101部分与601、701、801、901和1001部分相类似，在此不作赘述。
在1102部分：选择以起始点P0和交点O为端点的直线段上的三等分点分别作为第一中间点P1和第二中间点P2。
由于直线曲率为0，所以P0、P1、P2共线。
由于三等分点有两个，在一个示例中，可选取靠近交点的三等分点作为第二中间点P2，选择靠近起始点P0的三等分点作为第一中间点P1。
在1103部分：确定第四中间点P4位于终止点P5的切线上，并且，线段P4P5的长度为
其中，r2表示ln轨迹段的半径，α2表示终止点处的圆心角。
1103部分与904、1004部分相类似，在此不作赘述。
在1104部分：确定第三中间点P3端点为，以终止点和交点O为端点的圆弧段上的三等分点。
由于三等分点有两个，在一个示例中，可选取靠近交点的三等分点作为中间点P3。
1104部分与1005部分相类似，在此不作赘述。
在1105部分：根据上述特征点构建五次Bezier曲线，作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段(图11b虚线部分即为平滑过渡轨迹段)。
请参见图11c，通过图11a所示的方法，在起始点曲率为0，在终止点曲率为半径(50)的倒数0.02，进而实现在拐点处加速度连续。
图12a为n＝5且第一CP运动轨迹段为直线段，第二CP运动轨迹段为圆弧段时，本发明实施例提供的又一种机器人运动轨迹规划方法示例性流程图。
图12b显示的是第一CP运动轨迹段(lc)和第二CP运动轨迹段(ln)以及平滑过渡轨迹段，lc轨迹段和ln轨迹段相交于交点O。
在1201部分：确定lc轨迹段的拐出点为平滑过渡轨迹段的起始点P0，确定ln轨迹段的拐入点为平滑过渡轨迹段的终止点P5。
1201部分与601、701、801、901、1001和1101部分相类似，在此不作赘述。
在1202部分：选择以起始点P0和交点O为端点的直线段上的三等分点分别作为(第一)中间点P1和第二中间点P2。
1202部分与1102部分相类似，在此不作赘述。
在1203部分：选择以终止点P5和交点O为端点的直线段上的三等分点分别作为第三中间点P3和第四中间点P4。
由于直线曲率为0，所以P3、P4、P5共线。
由于三等分点有两个，在一个示例中，可选取靠近交点的三等分点作为第三中间点P3，选择靠近终止点P5的三等分点作为第四中间点P4。
在1204部分：根据上述特征点构建五次Bezier曲线，作为第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段(图12b虚线部分即为平滑过渡轨迹段)。
请参见图12c，通过图12a所示的方法，在始末点曲率皆为0，进而实现在拐点处加速度连续。
至于n＝5且第一CP运动轨迹段为圆弧段和第二CP运动轨迹段为直线段时，平滑过渡轨迹段的起始点P0、终止点P5的确定方式可参见601、701、801、901、1001和1201部分，第一中间点P1和第二中间点P2的确定方式可参见1002部分和1003部分，第三中间点P3和第四中间点P4的确定方式可参见1203部分，在此不作赘述。
图13示出了上述实施例中所涉及的机器人运动轨迹规划装置的一种可能的结构示意图，包括：
特征点确定单元131，用于根据第一CP运动轨迹段和第二CP运动轨迹段确定n次(n不小于4)Bezier曲线的特征点。
平滑过渡轨迹构建单元132，用于根据所述特征点构建n次(n不小于4)Bezier曲线，作为所述第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段，所述平滑过渡轨迹段的起始点为所述第一CP运动轨迹段的拐出点，所述平滑过渡轨迹段的终止点为所述第二CP运动轨迹段的拐入点；其中，所述平滑过渡轨迹段与所述第一CP运动轨迹段在所述起始点上的切矢量方向相同且曲率相同；所述平滑过渡轨迹段与所述第二CP运动轨迹段在所述终止点上的切矢量方向相同且曲率相同。
此外，还可包括规划单元133，用于在所述特征点确定单元确定n次(n不小于4)Bezier曲线的特征点之前，规划连续的CP运动轨迹段。
其中，特征点确定单元131可用于执行图4所示的401部分、图5所示的501部分、图6a所示的601-604部分、图7a所示的701-704部分、图8a所示的801-804部分，图9a所示的901-904部分，图10a所示的1001-1005部分，图11a所示的1101-1104部分，以及，图12a所示的1201-1203部分。
平滑过渡轨迹构建单元132可用于执行图4所示的402部分、图5所示的502部分、图6a所示的605部分、图7a所示的705部分、图8a所示的805部分，图9a所示的905部分，图10a所示的1006部分，图11a所示的1105部分，以及，图12a所示的1204部分。
规划单元133可用于执行图5所示的500部分。
图14示出了上述实施例中所涉及的机器人的一种可能的结构示意图，包括：
总线、控制器/处理器1、存储器2、通信接口3、输入设备4、输出设备5和被控器件6。处理器1、存储器2、通信接口3、输入设备4、输出设备5和被控器件6可通过总线相互连接。其中：
总线可包括一通路，在计算机***各个部件之间传送信息。
控制器/处理器1(机器人运动轨迹规划装置)可以是通用处理器，例如通用中央处理器(CPU)、网络处理器(NetworkProcessor，简称NP)、微处理器等，也可以是特定应用集成电路(application-specific integrated circuit，ASIC)，或一个或多个用于控制本发明方案程序执行的集成电路。还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。控制器/处理器1也可以是实现计算功能的组合，例如包含一个或多个微处理器组合，DSP和微处理器的组合等等。
存储器2中保存有执行本发明技术方案的程序，还可以保存有操作***和其他应用程序。具体地，程序可以包括程序代码，程序代码包括计算机操作指令。更具体的，存储器2可以是只读存储器(read-only memory，ROM)、可存储静态信息和指令的其他类型的静态存储设备、随机存取存储器(randomaccess memory，RAM)、可存储信息和指令的其他类型的动态存储设备、磁盘存储器等等。
输入设备4可包括接收用户输入的数据和信息的装置，例如键盘、鼠标、摄像头、扫描仪、光笔、语音输入装置、触摸屏等。
输出设备5可包括允许输出信息给用户的装置，例如显示屏、打印机、扬声器等。
通信接口3可包括使用任何收发器一类的装置，以便与其他设备或通信网络通信，如以太网，无线接入网(RAN)，无线局域网(WLAN)等。
控制器/处理器1可用于执行图4、5、6a、7a、8a、9a、10a、11a、12a中涉及机器人运动轨迹规划装置的处理过程和/或用于本申请所描述的技术的其他过程。控制器/处理器也可用于实现前述特征点确定单元131、平滑过渡轨迹构建单元132和规划单元133的功能。
被控器件6可用于按照控制器/处理器1规划出的连续的两CP运动轨迹段及平滑过渡轨迹段运行。
可以理解的是，图14仅仅示出了机器人的简化设计。在实际应用中，机器人可以包含任意数量的发射器，接收器，处理器，控制器，存储器，通信接口等，而所有可以实现本发明的机器人都在本发明的保护范围之内。
上述所有实施例所公开的方法和装置可运用于机械臂控制中，实现机械臂在笛卡尔空间下两条连续运行轨迹中的平滑转接。当然也可以拓展到任何需要规划路径的设备中，如小车、飞行器等等。
结合本发明公开内容所描述的方法或者算法的步骤可以硬件的方式来实现，也可以是由处理器执行软件指令的方式来实现。软件指令可以由相应的软件模块组成，软件模块可以被存放于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、移动硬盘、CD-ROM或者本领域熟知的任何其它形式的存储介质中。一种示例性的存储介质耦合至处理器，从而使处理器能够从该存储介质读取信息，且可向该存储介质写入信息。当然，存储介质也可以是处理器的组成部分。处理器和存储介质可以位于ASIC中。
本领域技术人员应该可以意识到，在上述一个或多个示例中，本发明所描述的功能可以用硬件、软件、固件或它们的任意组合来实现。当使用软件实现时，可以将这些功能存储在计算机可读介质中或者作为计算机可读介质上的一个或多个指令或代码进行传输。计算机可读介质包括计算机存储介质和通信介质，其中通信介质包括便于从一个地方向另一个地方传送计算机程序的任何介质。存储介质可以是通用或专用计算机能够存取的任何可用介质。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的技术方案的基础之上，所做的任何修改、等同替换、改进等，均应包括在本发明的保护范围之内。

Claims (23)

1.一种机器人运动轨迹规划方法，其特征在于，至少用于实现连续的两CP运动轨迹段之间的平滑转接，所述连续的两CP运动轨迹段包括第一CP运动轨迹段和第二CP运动轨迹段；
所述方法包括：
根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点，n为不小于4的正整数；
根据所述特征点构建n次Bezier曲线，作为所述第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段，所述平滑过渡轨迹段的起始点为所述第一CP运动轨迹段的拐出点，所述平滑过渡轨迹段的终止点为所述第二CP运动轨迹段的拐入点；
其中，所述平滑过渡轨迹段与所述第一CP运动轨迹段在所述起始点上的切矢量方向相同且曲率相同；所述平滑过渡轨迹段与所述第二CP运动轨迹段在所述终止点上的切矢量方向相同且曲率相同。
2.如权利要求1所述的方法，其特征在于，在所述确定n次Bezier曲线的特征点之前，还包括：规划连续的CP运动轨迹段。
3.如权利要求1或2所述方法，其特征在于，
所述n＝4，所述特征点包括起始点P0、终止点P4，以及，第一至第三中间点P1、P2、P3；所述第一CP运动轨迹段和第二CP运动轨迹段的交点表示为O点；
所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点包括：
确定所述交点O为所述第二中间点P2。
4.如权利要求3所述的方法，其特征在于，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点还包括：
若所述第一CP运动轨迹段为直线段，选择所述交点O和起始点P0之间的线段上的一点作为所述第一中间点P1；
若所述第一CP运动轨迹段为圆弧段，确定第一中间点P1位于起始点P0的切线上，并且，线段P1P0的长度为其中，r1表示所述第一CP运动轨迹段的半径，α1表示所述起始点处的圆心角。
5.如权利要求3或4所述的方法，其特征在于，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点还包括：
若所述第二CP运动轨迹段为直线段，选择所述交点O和终止点P4之间的线段上的一点作为所述第三中间点P3；
若所述第二CP运动轨迹段为圆弧段，确定所述第三中间点P3位于终止点P4的切线上，并且，线段P3P4的长度为其中，r2表示所述第二CP运动轨迹段的半径，α2表示终止点处的圆心角。
6.如权利要求1或2所述方法，其特征在于，
所述n＝5，所述特征点包括起始点P0、终止点P5以及第一至第四中间点P1、P2、P3、P4；所述第一CP运动轨迹段和第二CP运动轨迹段的交点表示为O点。
7.如权利要求6所述方法，其特征在于，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点包括：
若所述第一CP运动轨迹段为圆弧段，确定所述第一中间点P1位于所述起始点P0的切线上，并且，线段P1P0的长度为其中，r1表示所述第一CP运动轨迹段的半径，α1表示所述起始点处的圆心角；
确定所述第二中间点P2为，以所述起始点P0和交点O为端点的圆弧段上的三等分点。
8.如权利要求6所述方法，其特征在于，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点包括：
若所述第一CP运动轨迹段为直线段，选择以所述起始点P0和交点O为端点的直线段上的三等分点分别作为第一中间点P1和第二中间点P2。
9.如权利要求7或8所述方法，其特征在于，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点还包括：
若所述第二CP运动轨迹段为圆弧段，确定所述第三中间点P3为，以所述终止点和交点O为端点的圆弧段上的三等分点；
确定所述第四中间点P4位于终止点P5的切线上，并且，线段P4P5的长度为其中，r2表示所述第二CP运动轨迹段的半径，α2表示所述终止点处的圆心角。
10.如权利要求7或8所述方法，其特征在于，所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点还包括：
若所述第二CP运动轨迹段为直线段，选择以所述终止点和交点O为端点的直线段上的三等分点分别作为所述第三中间点P3和第四中间点P4。
11.一种机器人运动轨迹规划装置，其特征在于，至少用于实现连续的两CP运动轨迹段之间的平滑转接，所述连续的两CP运动轨迹段包括第一CP运动轨迹段和第二CP运动轨迹段；所述装置包括：
特征点确定单元，用于根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点，n为不小于4的正整数；
平滑过渡轨迹构建单元，用于根据所述特征点构建n次Bezier曲线，作为所述第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段，所述平滑过渡轨迹段的起始点为所述第一CP运动轨迹段的拐出点，所述平滑过渡轨迹段的终止点为所述第二CP运动轨迹段的拐入点；
其中，所述平滑过渡轨迹段与所述第一CP运动轨迹段在所述起始点上的切矢量方向相同且曲率相同；所述平滑过渡轨迹段与所述第二CP运动轨迹段在所述终止点上的切矢量方向相同且曲率相同。
12.如权利要求11所述的装置，其特征在于，还包括：
规划单元，用于在所述特征点确定单元确定n次Bezier曲线的特征点之前，规划连续的CP运动轨迹段。
13.如权利要求11或12所述装置，其特征在于，所述n＝4，所述特征点包括起始点P0、终止点P4，以及，第一至第三中间点P1、P2、P3；所述第一CP运动轨迹段和第二CP运动轨迹段的交点表示为O点；
在所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点的方面，所述平滑过渡轨迹构建单元用于：确定所述交点O为所述第二中间点P2。
14.如权利要求13所述的装置，其特征在于，在所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点的方面，所述平滑过渡轨迹构建单元还用于：
若所述第一CP运动轨迹段为直线段，选择所述交点O和起始点P0之间的线段上的一点作为所述第一中间点P1；
若所述第一CP运动轨迹段为圆弧段，确定第一中间点P1位于起始点P0的切线上，并且，线段P1P0的长度为其中，r1表示所述第一CP运动轨迹段的半径，α1表示所述起始点处的圆心角。
15.如权利要求13或14所述的装置，其特征在于，在所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点的方面，所述平滑过渡轨迹构建单元还用于：
若所述第二CP运动轨迹段为直线段，选择所述交点O和终止点P4之间的线段上的一点作为所述第三中间点P3；
若所述第二CP运动轨迹段为圆弧段，确定所述第三中间点P3位于终止点P4的切线上，并且，线段P3P4的长度为其中，r2表示所述第二CP运动轨迹段的半径，α2表示终止点处的圆心角。
16.如权利要求11或12所述的装置，其特征在于，所述n＝5，所述特征点包括起始点P0、终止点P5以及第一至第四中间点P1、P2、P3、P4；所述第一CP运动轨迹段和第二CP运动轨迹段的交点表示为O点。
17.如权利要求16所述的装置，其特征在于，在所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点的方面，所述平滑过渡轨迹构建单元用于：
若所述第一CP运动轨迹段为圆弧段，确定所述第一中间点P1位于所述起始点P0的切线上，并且，线段P1P0的长度为其中，r1表示所述第一CP运动轨迹段的半径，α1表示所述起始点处的圆心角；
确定所述第二中间点P2为，以所述起始点P0和交点O为端点的圆弧段上的三等分点。
18.如权利要求16所述的装置，其特征在于，在所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点的方面，所述平滑过渡轨迹构建单元用于：若所述第一CP运动轨迹段为直线段，选择以所述起始点P0和交点O为端点的直线段上的三等分点分别作为第一中间点P1和第二中间点P2。
19.如权利要求17或18所述的装置，其特征在于，在所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点的方面，所述平滑过渡轨迹构建单元还用于：
若所述第二CP运动轨迹段为圆弧段，确定所述第三中间点P3为，以所述终止点和交点O为端点的圆弧段上的三等分点；
确定所述第四中间点P4位于终止点P5的切线上，并且，线段P4P5的长度为其中，r2表示所述第二CP运动轨迹段的半径，α2表示所述终止点处的圆心角。
20.如权利要求17或18所述的装置，其特征在于，在所述根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点的方面，所述平滑过渡轨迹构建单元还用于：若所述第二CP运动轨迹段为直线段，选择以所述终止点和交点O为端点的直线段上的三等分点分别作为所述第三中间点P3和第四中间点P4。
21.一种机器人运动轨迹规划装置，其特征在于，至少用于实现连续的两CP运动轨迹段之间的平滑转接，所述连续的两CP运动轨迹段包括第一CP运动轨迹段和第二CP运动轨迹段；
所述装置包括：处理器和存储器，所述处理器通过运行存储在所述存储器内的软件程序、调用存储在所述存储器内的数据，至少执行如下步骤：
根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点，n为不小于4的正整数；
根据所述特征点构建n次Bezier曲线，作为所述第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段，所述平滑过渡轨迹段的起始点为所述第一CP运动轨迹段的拐出点，所述平滑过渡轨迹段的终止点为所述第二CP运动轨迹段的拐入点；
其中，所述平滑过渡轨迹段与所述第一CP运动轨迹段在所述起始点上的切矢量方向相同且曲率相同；所述平滑过渡轨迹段与所述第二CP运动轨迹段在所述终止点上的切矢量方向相同且曲率相同。
22.一种机器人，其特征在于，包括机器人运动轨迹规划装置和被控器件，其中：
所述机器人运动轨迹规划装置用于：根据所述第一CP运动轨迹段和第二CP运动轨迹段确定n次Bezier曲线的特征点，n为不小于4的正整数；根据所述特征点构建n次Bezier曲线，作为所述第一CP运动轨迹段和第二CP运动轨迹段之间的平滑过渡轨迹段，所述平滑过渡轨迹段的起始点为所述第一CP运动轨迹段的拐出点，所述平滑过渡轨迹段的终止点为所述第二CP运动轨迹段的拐入点；其中，所述平滑过渡轨迹段与所述第一CP运动轨迹段在所述起始点上的切矢量方向相同且曲率相同；所述平滑过渡轨迹段与所述第二CP运动轨迹段在所述终止点上的切矢量方向相同且曲率相同；
所述被控器件用于：按照所述机器人运动轨迹规划装置规划出的连续的两CP运动轨迹段及平滑过渡轨迹段运行。
23.如权利要求22所述的机器人，其特征在于，所述被控器件具体为：操作臂。
CN201780002229.0A 2017-01-04 2017-01-04 机器人运动轨迹规划方法及相关装置 Active CN107980108B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2017/070133 WO2018126355A1 (zh) 2017-01-04 2017-01-04 机器人运动轨迹规划方法及相关装置

Publications (2)

Publication Number Publication Date
CN107980108A true CN107980108A (zh) 2018-05-01
CN107980108B CN107980108B (zh) 2021-07-09

Family

ID=62006173

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201780002229.0A Active CN107980108B (zh) 2017-01-04 2017-01-04 机器人运动轨迹规划方法及相关装置

Country Status (2)

Country Link
CN (1) CN107980108B (zh)
WO (1) WO2018126355A1 (zh)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109773791A (zh) * 2019-01-31 2019-05-21 北京华航唯实机器人科技股份有限公司 路径生成方法及装置
CN109884979A (zh) * 2019-01-31 2019-06-14 北京华航唯实机器人科技股份有限公司 设备运动点的处理方法及装置
CN110370267A (zh) * 2018-09-10 2019-10-25 北京京东尚科信息技术有限公司 用于生成模型的方法和装置
CN110488758A (zh) * 2019-08-09 2019-11-22 南京埃斯顿自动化股份有限公司 一种基于PLCopen规范的轨迹过渡方法
CN110730934A (zh) * 2018-08-01 2020-01-24 深圳市大疆创新科技有限公司 轨迹切换的方法和装置
CN110749332A (zh) * 2019-10-18 2020-02-04 广州文远知行科技有限公司 Rs曲线的曲率优化方法、装置、计算机设备及存储介质
CN111127590A (zh) * 2019-12-26 2020-05-08 新奥数能科技有限公司 一种二阶贝塞尔曲线绘制方法及装置
WO2020093253A1 (zh) * 2018-11-06 2020-05-14 深圳配天智能技术研究院有限公司 机器人的运动控制方法、控制***和存储装置
WO2020143357A1 (zh) * 2019-01-08 2020-07-16 苏州极目机器人科技有限公司 飞行器作业路径规划方法以及控制装置和控制设备
CN111670093A (zh) * 2018-11-06 2020-09-15 深圳配天智能技术研究院有限公司 机器人的运动控制方法、控制***和存储装置
CN112731933A (zh) * 2020-12-24 2021-04-30 江苏新冠亿科技有限公司 一种自主规划路径的agv跟踪控制方法、设备及存储介质
CN112861238A (zh) * 2021-03-02 2021-05-28 广联达科技股份有限公司 生成道路轮廓线的方法、装置、设备及可读存储介质
CN113741336A (zh) * 2021-09-10 2021-12-03 上海景吾智能科技有限公司 基于实时运动的笛卡尔空间轨迹规划方法和***
WO2021244545A1 (zh) * 2020-06-05 2021-12-09 深圳市道通智能航空技术股份有限公司 一种无人机制导方法、无人机及存储介质
CN114055459A (zh) * 2020-08-06 2022-02-18 库卡机器人(广东)有限公司 轨迹规划的方法、装置、电子设备及存储介质
CN115509240A (zh) * 2022-11-21 2022-12-23 合肥井松智能科技股份有限公司 一种基于可通行区域的路径规划方法
WO2023016194A1 (zh) * 2021-08-11 2023-02-16 灵动科技（北京）有限公司 用于移动机器人的路径规划方法及程序产品
CN116149260A (zh) * 2023-01-31 2023-05-23 苏州浩智工业控制技术有限公司 3c金属加工中圆弧与圆弧平滑过渡的方法及***

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110502010B (zh) * 2019-08-15 2021-06-04 同济大学 一种基于贝塞尔曲线的移动机器人室内自主导航控制方法
TWI716127B (zh) * 2019-10-01 2021-01-11 東元電機股份有限公司 機器人及其路徑插值規劃命令產生系統
CN112621739B (zh) * 2019-10-08 2022-03-15 东元电机股份有限公司 机器人及其路径插值规划命令产生***
CN113467498B (zh) * 2021-07-14 2022-07-01 西北工业大学 一种基于Bezier-凸优化的运载火箭上升段轨迹规划方法
CN114234968B (zh) * 2021-12-17 2023-12-05 江西洪都航空工业集团有限责任公司 一种基于a星算法的移动机器人自主导航方法
CN114310921B (zh) * 2022-03-16 2022-06-10 珞石（北京）科技有限公司 一种最小曲率的直线过渡路径生成方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102207736A (zh) * 2010-03-31 2011-10-05 中国科学院自动化研究所 基于贝塞尔曲线的机器人路径规划方法及装置
EP2378383A2 (de) * 2010-03-25 2011-10-19 KUKA Laboratories GmbH Verfahren zum Betreiben eines holonomen/omnidirektionalen Flurförderfahrzeugs
KR20120098152A (ko) * 2011-02-28 2012-09-05 한국과학기술연구원 이동 로봇의 경로 생성 시스템
CN103528585A (zh) * 2013-09-26 2014-01-22 中北大学 一种不等距分割可通行区域的路径规划方法
CN103907068A (zh) * 2011-09-02 2014-07-02 布鲁克斯自动化公司 用于机器人传送装置的时间最佳轨迹
CN104808688A (zh) * 2015-04-28 2015-07-29 武汉大学 一种无人机曲率连续可调路径规划方法
CN105435997A (zh) * 2015-12-14 2016-03-30 江苏大学 基于Bézier曲线的喷涂机器人路径规划方法
CN105573315A (zh) * 2015-12-01 2016-05-11 珞石（北京）科技有限公司 用于工业机器人的笛卡尔空间轨迹的几何平滑方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2378383A2 (de) * 2010-03-25 2011-10-19 KUKA Laboratories GmbH Verfahren zum Betreiben eines holonomen/omnidirektionalen Flurförderfahrzeugs
CN102207736A (zh) * 2010-03-31 2011-10-05 中国科学院自动化研究所 基于贝塞尔曲线的机器人路径规划方法及装置
KR20120098152A (ko) * 2011-02-28 2012-09-05 한국과학기술연구원 이동 로봇의 경로 생성 시스템
CN103907068A (zh) * 2011-09-02 2014-07-02 布鲁克斯自动化公司 用于机器人传送装置的时间最佳轨迹
CN103528585A (zh) * 2013-09-26 2014-01-22 中北大学 一种不等距分割可通行区域的路径规划方法
CN104808688A (zh) * 2015-04-28 2015-07-29 武汉大学 一种无人机曲率连续可调路径规划方法
CN105573315A (zh) * 2015-12-01 2016-05-11 珞石（北京）科技有限公司 用于工业机器人的笛卡尔空间轨迹的几何平滑方法
CN105435997A (zh) * 2015-12-14 2016-03-30 江苏大学 基于Bézier曲线的喷涂机器人路径规划方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JILE JIAO: "Bezier Curve Based Path Planning for A Mobile Manipulator in Unknown Environments", 《INTERNATIONAL CONFERENCE ON ROBOTICS》 *
余伶俐等: "基于贝塞尔曲线的机器人非时间轨迹跟踪方法", 《仪器仪表学报》 *
陈成等: "基于四阶贝塞尔曲线的无人车可行轨迹规划", 《自动化学报》 *

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110730934A (zh) * 2018-08-01 2020-01-24 深圳市大疆创新科技有限公司 轨迹切换的方法和装置
WO2020024134A1 (zh) * 2018-08-01 2020-02-06 深圳市大疆创新科技有限公司 轨迹切换的方法和装置
CN110370267A (zh) * 2018-09-10 2019-10-25 北京京东尚科信息技术有限公司 用于生成模型的方法和装置
CN111684380A (zh) * 2018-11-06 2020-09-18 深圳配天智能技术研究院有限公司 机器人的运动控制方法、控制***和存储装置
WO2020093253A1 (zh) * 2018-11-06 2020-05-14 深圳配天智能技术研究院有限公司 机器人的运动控制方法、控制***和存储装置
CN111670093B (zh) * 2018-11-06 2023-12-22 深圳配天机器人技术有限公司 机器人的运动控制方法、控制***和存储装置
CN111670093A (zh) * 2018-11-06 2020-09-15 深圳配天智能技术研究院有限公司 机器人的运动控制方法、控制***和存储装置
WO2020143357A1 (zh) * 2019-01-08 2020-07-16 苏州极目机器人科技有限公司 飞行器作业路径规划方法以及控制装置和控制设备
CN109884979B (zh) * 2019-01-31 2020-02-21 北京华航唯实机器人科技股份有限公司 设备运动点的处理方法及装置
CN109884979A (zh) * 2019-01-31 2019-06-14 北京华航唯实机器人科技股份有限公司 设备运动点的处理方法及装置
CN109773791A (zh) * 2019-01-31 2019-05-21 北京华航唯实机器人科技股份有限公司 路径生成方法及装置
CN110488758A (zh) * 2019-08-09 2019-11-22 南京埃斯顿自动化股份有限公司 一种基于PLCopen规范的轨迹过渡方法
CN110749332B (zh) * 2019-10-18 2021-04-20 广州文远知行科技有限公司 Rs曲线的曲率优化方法、装置、计算机设备及存储介质
CN110749332A (zh) * 2019-10-18 2020-02-04 广州文远知行科技有限公司 Rs曲线的曲率优化方法、装置、计算机设备及存储介质
CN111127590A (zh) * 2019-12-26 2020-05-08 新奥数能科技有限公司 一种二阶贝塞尔曲线绘制方法及装置
CN111127590B (zh) * 2019-12-26 2023-06-20 新奥数能科技有限公司 一种二阶贝塞尔曲线绘制方法及装置
WO2021244545A1 (zh) * 2020-06-05 2021-12-09 深圳市道通智能航空技术股份有限公司 一种无人机制导方法、无人机及存储介质
CN114055459A (zh) * 2020-08-06 2022-02-18 库卡机器人(广东)有限公司 轨迹规划的方法、装置、电子设备及存储介质
CN112731933A (zh) * 2020-12-24 2021-04-30 江苏新冠亿科技有限公司 一种自主规划路径的agv跟踪控制方法、设备及存储介质
CN112861238A (zh) * 2021-03-02 2021-05-28 广联达科技股份有限公司 生成道路轮廓线的方法、装置、设备及可读存储介质
CN112861238B (zh) * 2021-03-02 2024-03-22 广联达科技股份有限公司 生成道路轮廓线的方法、装置、设备及可读存储介质
WO2023016194A1 (zh) * 2021-08-11 2023-02-16 灵动科技（北京）有限公司 用于移动机器人的路径规划方法及程序产品
CN113741336A (zh) * 2021-09-10 2021-12-03 上海景吾智能科技有限公司 基于实时运动的笛卡尔空间轨迹规划方法和***
CN115509240A (zh) * 2022-11-21 2022-12-23 合肥井松智能科技股份有限公司 一种基于可通行区域的路径规划方法
CN116149260A (zh) * 2023-01-31 2023-05-23 苏州浩智工业控制技术有限公司 3c金属加工中圆弧与圆弧平滑过渡的方法及***
CN116149260B (zh) * 2023-01-31 2024-04-05 苏州浩智工业控制技术有限公司 3c金属加工中圆弧与圆弧平滑过渡的方法及***

Also Published As

Publication number Publication date
WO2018126355A1 (zh) 2018-07-12
CN107980108B (zh) 2021-07-09

Similar Documents

Publication Publication Date Title
CN107980108A (zh) 机器人运动轨迹规划方法及相关装置
CN107980109A (zh) 机器人运动轨迹规划方法及相关装置
Kallmann Shortest Paths with Arbitrary Clearance from Navigation Meshes.
Ezair et al. Planning high order trajectories with general initial and final conditions and asymmetric bounds
Isaiah et al. Motion planning algorithms for the Dubins travelling salesperson problem
CN107952243B (zh) 路径确定方法及装置
WO2022161315A1 (zh) 机器人路径规划方法、运行方法、机器人以及介质
CN110231041B (zh) 一种车道切换的导航方法和装置
CN103731377B (zh) 处理报文的方法与设备
CN113033925B (zh) 用于控制自动驾驶车辆行驶、装置、电子设备和介质
CN111539345A (zh) 用于确定变道动作的方法、装置、设备及可读存储介质
WO2019217043A1 (en) Drag gesture animation
CN114514536A (zh) 分布式***中的神经网络训练
CN110017840A (zh) 导航数据的处理方法和装置、电子设备、计算机可读介质
US11392262B2 (en) Method, system, and non-transitory computer readable record medium for providing chatroom in 3D form
CN105184372B (zh) 知识网络构建方法及装置
WO2020034887A1 (zh) 智能体行进路径的确定方法及装置
JP2011203785A (ja) ポリゴン分割装置およびポリゴン分割方法
CN113890513A (zh) 量子硬件的控制方法及装置
CN111151003A (zh) 虚拟对象的控制方法和装置
JP3975986B2 (ja) 画像表示装置
CN112578796B (zh) 一种基于曲率约束的引导线生成方法及装置
JP2016049924A (ja) プログラム及び運転曲線作成装置
WO2020230249A1 (ja) ウォークスルー表示装置、ウォークスルー表示方法、及びウォークスルー表示プログラム
CN113534794A (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