CN110727271A - Robot motion primitive determining method and device - Google Patents

Robot motion primitive determining method and device Download PDF

Info

Publication number
CN110727271A
CN110727271A CN201911046662.1A CN201911046662A CN110727271A CN 110727271 A CN110727271 A CN 110727271A CN 201911046662 A CN201911046662 A CN 201911046662A CN 110727271 A CN110727271 A CN 110727271A
Authority
CN
China
Prior art keywords
robot
motion
primitive
time
determining
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.)
Pending
Application number
CN201911046662.1A
Other languages
Chinese (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.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
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 University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN201911046662.1A priority Critical patent/CN110727271A/en
Publication of CN110727271A publication Critical patent/CN110727271A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0212Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
    • G05D1/0223Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory involving speed control of the vehicle

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

The invention provides a robot motion primitive determining method and device, which can avoid the problems that the robot is easy to fall into local minimum or difficult to track a path when carrying out path tracking control. The method comprises the following steps: determining an initial state vector of the robot; determining a translational motion speed interval and a yaw velocity interval which are reached by the robot in the time of generating the motion elements, and equally dividing the translational motion speed interval into m-1 parts and equally dividing the yaw velocity interval into N-1 parts according to the number N of the generated motion elements, wherein N is m x N; and (3) taking the translational motion speed corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to form a control vector, determining the state of the robot according to the determined initial state vector of the robot and the kinematic constraint of the robot, and connecting the states of the robot to obtain a motion element generated under the control vector. The invention relates to the technical field of robot path planning.

Description

Robot motion primitive determining method and device
Technical Field
The invention relates to the technical field of robot path planning, in particular to a method and a device for determining a robot motion primitive.
Background
Path planning is a basic premise for realizing autonomous navigation of a robot. The robot acquires perception information of the environment by means of a sensor in the movement process, and meanwhile, an environment map is created. In the path planning process, the robot motion primitives determine the connection relationship between nodes, and one state node can be converted to another state node through the motion primitives. And (3) searching a series of motion primitives transformed from an initial state to a target state according to the barrier information marked on the environment map by the path planning, namely planning to generate a collision-free safe path.
The motion primitives are used as a representation mode of a node connection relationship in path planning, and usually a robot is regarded as a particle and generated by adopting straight line connection, such as a 4-connection grid and an 8-connection grid, and the motion primitives generated by adopting the method do not consider the kinematic constraint conditions of the robot, so that the robot is easy to fall into local minimum or difficult to track a path when performing path tracking control.
The Chinese patent with the application number of 201610348356.3 discloses a mobile robot path planning method based on a dynamic motion element learning model, the method adopts a handle to control the motion of a robot, records the motion track of the robot as a sample of the dynamic motion element model, and obtains the parameters of the dynamic motion element model by establishing the dynamic motion element model and training by utilizing the track sample, and the motion element generation process is complex and time-consuming and is greatly influenced by the manual operation level.
Disclosure of Invention
The invention aims to solve the technical problem of providing a robot motion primitive determining method and device to solve the problems that in the prior art, the motion primitive is generated without considering the kinematic constraint condition of a robot, so that the robot is easy to fall into local minimum or difficult to track a path when performing path tracking control, and the motion primitive generating process is complex and time-consuming.
In order to solve the above technical problem, an embodiment of the present invention provides a method for determining a robot motion primitive, including:
determining an initial state vector of the robot;
determining a translational motion speed interval and a yaw velocity interval which are reached by the robot in the time of generating the motion elements, and equally dividing the translational motion speed interval into m-1 parts and equally dividing the yaw velocity interval into N-1 parts according to the number N of the generated motion elements, wherein N is m x N;
and taking the translational motion speed corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to form a control vector, determining the state of the robot according to the determined initial state vector of the robot and the kinematic constraint of the robot, and connecting the states of the robot to obtain a motion element generated under the control vector, wherein i is 1, 2.
Further, the starting state vector is zs=[xs,yss]TWherein z represents a state vector, x and y represent displacements of the robot along the x direction and the y direction respectively, theta represents a heading angle of the robot, subscript s represents a start, and superscript T represents a matrix transposition.
Further, the step of determining a translational motion speed interval and a yaw rate interval reached by the robot within the time of generating the motion primitives, and equally dividing the translational motion speed interval into m parts and equally dividing the yaw rate interval into N parts according to the number N of generated motion primitives comprises the following steps:
determining the translation motion speed interval [ v ] reached by the robot in the time t generated by the motion primitives-at,vs+at]And yaw angular velocity interval [ omega ]s-αt,ωs+αt]Wherein v represents the speed of the translational motion of the robot, ω represents the yaw rate of the robot, and a and α represent the acceleration at the translational motion of the robot and the yaw acceleration at the yaw motion, respectively;
according to the number N of the robot motion elements, equally dividing a translational motion speed interval into m-1 parts, equally dividing a yaw angle speed interval into N-1 parts, and satisfying the following conditions: n — m × N, then: translational motion speed corresponding to endpoint iYaw rate corresponding to endpoint j
Figure BDA0002254300380000022
Further, the step of taking the translational motion velocity corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to form a control vector, determining the state of the robot according to the determined initial state vector of the robot and the kinematic constraint of the robot, and connecting the states of the robot to obtain a motion primitive generated under the control vector includes:
within the time of motion primitive generation, with u ═ vij]TThe constant speed is used as a control vector, and the robot is determined at the sampling time t according to the determined initial state vector of the robot and the kinematic constraint of the robotfThe time-dependent state is:
wherein u represents a control vector, (x)f,yff) Indicating that the robot is at the sampling time tfState vector of (x)f,yff) Is a discrete point on the motion primitive, tfAny sampling moment in time t generated by a motion primitive is designated;
connecting state vectors at different sampling moments in the time of generating the motion primitive to obtain a control vector u ═ vij]TThe next generated motion primitive.
Further, the state vectors at different sampling time are connected in the time of generating the motion primitive, and a control vector u ═ v is obtainedij]TThe following generated motion primitives include:
in the time of generating the motion primitive, if the termination state vector is not at the grid center of the environment map, translating the termination state vector to the grid center by adopting a rounding method;
connecting the state vectors at different sampling times to obtain a control vector u ═ vij]TThe next generated motion primitive.
An embodiment of the present invention further provides a device for determining a robot motion primitive, including:
the determining module is used for determining an initial state vector of the robot;
the robot comprises an equally dividing module, a judging module and a judging module, wherein the equally dividing module is used for determining a translational motion speed interval and a yaw angular velocity interval which are reached by the robot in the time of generating motion primitives, and equally dividing the translational motion speed interval into m-1 parts and equally dividing the yaw angular velocity interval into N-1 parts according to the generated number N of the motion primitives, wherein N is m x N;
and the generating module is used for taking the translational motion speed corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to form a control vector, determining the state of the robot according to the determined initial state vector of the robot and the kinematic constraint of the robot, and connecting the states of the robot to obtain a motion element generated under the control vector, wherein i is 1,2, and.
Further, the starting state vector is zs=[xs,yss]TWherein z represents a state vector, x and y represent displacements of the robot along the x direction and the y direction respectively, theta represents a heading angle of the robot, subscript s represents a start, and superscript T represents a matrix transposition.
Further, the halving module comprises:
a first determining unit for determining the translational motion speed interval [ v ] reached by the robot in the time t of the generation of the motion primitives-at,vs+at]And yaw angular velocity interval [ omega ]s-αt,ωs+αt]Wherein v represents the speed of the translational motion of the robot, ω represents the yaw rate of the robot, and a and α represent the acceleration at the translational motion of the robot and the yaw acceleration at the yaw motion, respectively;
and the equally dividing unit is used for equally dividing the translational motion speed interval into m-1 parts and the yaw velocity interval into N-1 parts according to the generated number N of the robot motion elements, and meets the following requirements: n — m × N, then: translational motion speed corresponding to endpoint iThe transverse line corresponding to the end point jAngular velocity of oscillation
Figure BDA0002254300380000042
Further, the generating module includes:
a second determination unit for generating a motion primitive with u ═ v in timeij]TThe constant speed is used as a control vector, and the robot is determined at the sampling time t according to the determined initial state vector of the robot and the kinematic constraint of the robotfThe time-dependent state is:
Figure BDA0002254300380000043
wherein u represents a control vector, (x)f,yff) Indicating that the robot is at the sampling time tfState vector of (x)f,yff) Is a discrete point on the motion primitive, tfAny sampling moment in time t generated by a motion primitive is designated;
a generating unit for connecting the state vectors at different sampling time points within the time of generating the motion primitive to obtain a control vector u ═ vij]TThe next generated motion primitive.
Further, the generating unit is configured to, during the time of generating the motion primitive, if the termination state vector is not at the grid center of the environment map, translate the termination state vector to the grid center by using a rounding method; also for concatenating the state vectors at different sampling instants, resulting in a control vector u ═ vij]TThe next generated motion primitive.
The technical scheme of the invention has the following beneficial effects:
in the scheme, the initial state vector of the robot is determined; determining a translational motion speed interval and a yaw velocity interval which are reached by the robot in the time of generating the motion elements, and equally dividing the translational motion speed interval into m-1 parts and equally dividing the yaw velocity interval into N-1 parts according to the number N of the generated motion elements, wherein N is m x N; and (3) taking the translational motion speed corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to form a control vector, determining the state of the robot according to the determined initial state vector of the robot and the kinematic constraint of the robot, and connecting the states of the robot to obtain a motion element generated under the control vector. In this way, the reachable state of the robot in the time of generating the motion primitives is determined according to the kinematic constraints of the robot, so that the motion primitives under various different control vectors are generated, the problems that the robot is easy to fall into local minimum or difficult to track paths when the robot carries out path tracking control can be avoided, and the generation process of the motion primitives is simple.
Drawings
Fig. 1 is a schematic flow chart of a robot motion primitive determination method according to an embodiment of the present invention;
fig. 2 is a schematic diagram of a coordinate system of a robot according to an embodiment of the present invention;
FIG. 3 is a diagram of a control vector u according to an embodiment of the present invention9=[1,0.5]TA schematic diagram of motion elements generated under action;
FIG. 4 shows a starting state vector z of the robot according to the embodiment of the present inventions=[0,0,0]TA schematic diagram of a lower motion primitive;
fig. 5 is a schematic structural diagram of a robot motion primitive determination apparatus according to an embodiment of the present invention.
Detailed Description
In order to make the technical problems, technical solutions and advantages of the present invention more apparent, the following detailed description is given with reference to the accompanying drawings and specific embodiments.
The invention provides a method and a device for determining a motion element of a robot, aiming at the problems that the motion element is generated without considering the kinematic constraint condition of the robot, so that the robot is easy to fall into local minimum or difficult in path tracking when performing path tracking control, and the generation process of the motion element is complex and time-consuming.
Example one
As shown in fig. 1, a method for determining a robot motion primitive provided by an embodiment of the present invention includes:
s101, determining an initial state vector of the robot;
s102, determining a translational motion speed interval and a yaw velocity interval which are reached by the robot within the time of generating motion elements, and equally dividing the translational motion speed interval into m-1 parts and equally dividing the yaw velocity interval into N-1 parts according to the generated number N of the motion elements, wherein N is m × N;
and S103, taking the translational motion speed corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to form a control vector, determining the state of the robot according to the determined initial state vector of the robot and the kinematic constraint of the robot, and connecting the states of the robot to obtain a motion element generated under the control vector, wherein i is 1,2, and m, j is 1, 2.
The robot motion primitive determining method of the embodiment of the invention determines the initial state vector of the robot; determining a translational motion speed interval and a yaw velocity interval which are reached by the robot in the time of generating the motion elements, and equally dividing the translational motion speed interval into m-1 parts and equally dividing the yaw velocity interval into N-1 parts according to the number N of the generated motion elements, wherein N is m x N; and (3) taking the translational motion speed corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to form a control vector, determining the state of the robot according to the determined initial state vector of the robot and the kinematic constraint of the robot, and connecting the states of the robot to obtain a motion element generated under the control vector. In this way, the reachable state of the robot in the time of generating the motion primitives is determined according to the kinematic constraints of the robot, so that the motion primitives under various different control vectors are generated, the problems that the robot is easy to fall into local minimum or difficult to track paths when the robot carries out path tracking control can be avoided, and the generation process of the motion primitives is simple.
The robot motion primitive determining method provided by the embodiment of the invention can be used for generating motion primitives for various differentially driven robots, so that the robot path planning is realized.
In this embodiment, it is assumed that the robot works on a flat ground, and therefore, the robot can be operated on a machineThe motion of the person is simplified by considering only the translation of the robot in the x, y (also can be written: X, Y) direction and the yaw motion around the z-axis, ignoring the displacement in the vertical direction, the pitch angle and the roll angle, as shown in fig. 2, where v denotes the velocity of the translation motion of the robot and v denotes the velocity of the translation motion of the robotx、vyThe translational motion speeds of the robot along the x direction and the y direction are respectively, omega represents the yaw velocity of the robot, and then x and y should satisfy the kinematic constraint condition of the formula (1):
wherein x and y respectively represent the displacement of the robot along the x direction and the y direction,
Figure BDA0002254300380000062
respectively representing first derivatives of x and y to time, v representing the translational motion speed of the robot, theta representing the heading angle of the robot, and theta being equal to 0 and 2 pi;
if z is [ x, y, theta ]]TRepresents a state vector, u ═ v, ω]TThe control vector is expressed, omega represents the yaw velocity of the robot, and the state transition equation expression of the robot can be obtained according to the formula (1):
Figure BDA0002254300380000071
in this embodiment, assume that the initial state vector of the robot is zs=[xs,yss]TThe starting control vector is us=[vss]TWherein z represents a state vector, subscript s represents a start, and superscript T represents a matrix transpose; the acceleration of the robot during the translational motion and the yaw angular acceleration of the robot during the yaw motion are respectively represented by a and alpha, and the translational motion speed section which can be reached by the robot within the time (also called simulation time) t generated by the motion elements is [ v [ v ] ]s-at,vs+at](wherein v ∈ [ v [ ]s-at,vs+at]) And a yaw angular velocity interval of [ omega ]s-αt,ωs+αt](wherein ω ∈ [ ω ])s-αt,ωs+αt])。
In the embodiment, according to the number N of generated robot motion elements, a translational motion speed interval is equally divided into m-1 parts, a yaw angular velocity interval is equally divided into N-1 parts, and the following conditions are met: n — m × N, then: translational motion speed corresponding to endpoint i
Figure BDA0002254300380000072
Yaw rate corresponding to endpoint j
Figure BDA0002254300380000073
In the present embodiment, during the time of motion primitive generation, u ═ vij]TThe constant speed is used as a control vector, and the robot is determined at the sampling time t according to the determined initial state vector of the robot and the kinematic constraint of the robotfThe time-dependent state is:
Figure BDA0002254300380000074
wherein u represents a control vector, (x)f,yff) Indicating that the robot is at the sampling time tfState vector of time, tfRefers to any sampling moment in the simulation time t generated by the motion primitive;
in this embodiment, the sampling time t of the robot within the simulation time t of the motion primitive generation can be obtained according to equation (3)fSampling state of time (x)f,yff) These sampling states are discrete points on the motion primitive, which are connected to form different motion primitives, resulting in a control vector u ═ vij]TThe next generated motion primitive.
In an embodiment of the foregoing method for determining a robot motion primitive, further, the connecting the state vectors at different sampling time instants within the time of generating the motion primitive to obtain the motion primitive generated under the control vector includes:
in the time of generating the motion primitive, if the termination state vector is not at the grid center of the environment map, translating the termination state vector to the grid center by adopting a rounding method;
connecting the state vectors at different sampling times to obtain a control vector u ═ vij]TThe next generated motion primitive.
In order to better understand the robot motion primitive determination method provided by the embodiment of the invention, the robot starting state vector z is useds=[0,0,0]TAnd control vector u ═ vij]TGenerating robot motion primitives for example:
first, the acceleration a of the robot during translational motion is set to 0.5m/s2And the yaw angular acceleration of the robot during the yaw movement is alpha which is 0.25rad/s2. Within the time (simulation time) 2s of the generation of the motion primitive, the translational motion speed interval of the robot is [ -1m/s, 1m/s]The yaw angular velocity interval is [ -0.5rad/s, 0.5rad/s]。
Then, the translational velocity interval of the robot is equally divided into two parts to obtain 3 end points, namely-1 m/s, 0m/s and 1m/s respectively, and the yaw velocity interval is equally divided into two parts to obtain 3 end points, namely-0.5 rad/s, 0rad/s and 0.5rad/s respectively, so that the number of robot motion elements which can be generated is 9. Combining the translational motion speed corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to obtain control vectors u1=[-1,-0.5]T,u2=[-1,0]T,u3=[-1,0.5]T,u4=[0,-0.5]T,u5=[0,0]T,u6=[0,0.5]T,u7=[1,-0.5]T,u8=[1,0]TAnd u9=[1,0.5]T. In the simulation time 2s, the control vector is substituted into the formula (3), and the state vector of the robot can be calculated; specifically, the method comprises the following steps: sampling is carried out by taking 0.05s as a sampling period, and the motion primitives generated under the corresponding control vectors can be obtained by connecting the state vectors of the robot.
As shown in FIG. 3, FIG. 3 is an on-controlSystem vector u9=[1,0.5]TAnd (3) generating a motion primitive under action, wherein arrows in the figure 3 represent the heading angles of the robot state sampling points. For the starting state of the robot, the robot motion primitive is calculated by all the control vectors, and the motion primitive of the robot in the starting state can be obtained, as shown in fig. 4, and fig. 4 shows the vector z of the robot in the starting states=[0,0,0]TAnd the lower motion element comprises two pivot-turning motion elements and a static motion element.
Example two
The present invention also provides a specific embodiment of a robot motion primitive determining apparatus, which corresponds to the specific embodiment of the robot motion primitive determining method, and the robot motion primitive determining apparatus can achieve the object of the present invention by executing the flow steps in the specific embodiment of the method, so the explanation in the specific embodiment of the robot motion primitive determining method is also applicable to the specific embodiment of the robot motion primitive determining apparatus provided by the present invention, and will not be described again in the following specific embodiment of the present invention.
As shown in fig. 5, an embodiment of the present invention further provides a robot motion primitive determining apparatus, including:
the determining module 11 is used for determining an initial state vector of the robot;
an equally dividing module 12, configured to determine a translational motion speed interval and a yaw rate interval reached by the robot within the time of generating the motion primitives, and equally divide the translational motion speed interval into m-1 parts and the yaw rate interval into N-1 parts according to the generated number N of the motion primitives, where N is m × N;
and the generating module 13 is configured to obtain a control vector formed by the translational motion speed corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j, determine the state of the robot according to the determined starting state vector of the robot and the kinematic constraint of the robot, and connect the states of the robot to obtain a motion primitive generated under the control vector, where i is 1, 2.
The robot motion primitive determining device determines an initial state vector of a robot; determining a translational motion speed interval and a yaw velocity interval which are reached by the robot in the time of generating the motion elements, and equally dividing the translational motion speed interval into m-1 parts and equally dividing the yaw velocity interval into N-1 parts according to the number N of the generated motion elements, wherein N is m x N; and (3) taking the translational motion speed corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to form a control vector, determining the state of the robot according to the determined initial state vector of the robot and the kinematic constraint of the robot, and connecting the states of the robot to obtain a motion element generated under the control vector. In this way, the reachable state of the robot in the time of generating the motion primitives is determined according to the kinematic constraints of the robot, so that the motion primitives under various different control vectors are generated, the problems that the robot is easy to fall into local minimum or difficult to track paths when the robot carries out path tracking control can be avoided, and the generation process of the motion primitives is simple.
In an embodiment of the foregoing robot motion primitive determining apparatus, further, the initial state vector is zs=[xs,yss]TWherein z represents a state vector, x and y represent displacements of the robot along the x direction and the y direction respectively, theta represents a heading angle of the robot, subscript s represents a start, and superscript T represents a matrix transposition.
In an embodiment of the foregoing robot motion primitive determining apparatus, further, the equally dividing module includes:
a first determining unit for determining the translational motion speed interval [ v ] reached by the robot in the time t of the generation of the motion primitives-at,vs+at]And yaw angular velocity interval [ omega ]s-αt,ωs+αt]Wherein v represents the speed of the translational motion of the robot, ω represents the yaw rate of the robot, and a and α represent the acceleration at the translational motion of the robot and the yaw acceleration at the yaw motion, respectively;
an equally dividing unit for generating the number according to the robot motion primitiveN, equally dividing the translational motion speed interval into m-1 parts, equally dividing the yaw angle speed interval into N-1 parts, and satisfying the following conditions: n — m × N, then: translational motion speed corresponding to endpoint i
Figure BDA0002254300380000101
Yaw rate corresponding to endpoint j
Figure BDA0002254300380000102
In an embodiment of the foregoing robot motion primitive determining apparatus, further the generating module includes:
a second determination unit for generating a motion primitive with u ═ v in timeij]TThe constant speed is used as a control vector, and the robot is determined at the sampling time t according to the determined initial state vector of the robot and the kinematic constraint of the robotfThe time-dependent state is:
Figure BDA0002254300380000103
wherein u represents a control vector, (x)f,yff) Indicating that the robot is at the sampling time tfState vector of (x)f,yff) Is a discrete point on the motion primitive, tfAny sampling moment in time t generated by a motion primitive is designated;
a generating unit for connecting the state vectors at different sampling time points within the time of generating the motion primitive to obtain a control vector u ═ vij]TThe next generated motion primitive.
In an embodiment of the foregoing robot motion primitive determining apparatus, further, the generating unit is configured to, during the time of generating the motion primitive, if the termination state vector is not at the grid center of the environment map, translate the termination state vector to the grid center by a rounding method; also for concatenating the state vectors at different sampling instants, resulting in a control vector u ═ vij]TThe next generated motion primitive.
It is noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention as defined in the appended claims.

Claims (10)

1. A robot motion primitive determination method, comprising:
determining an initial state vector of the robot;
determining a translational motion speed interval and a yaw velocity interval which are reached by the robot in the time of generating the motion elements, and equally dividing the translational motion speed interval into m-1 parts and equally dividing the yaw velocity interval into N-1 parts according to the number N of the generated motion elements, wherein N is m x N;
and taking the translational motion speed corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to form a control vector, determining the state of the robot according to the determined initial state vector of the robot and the kinematic constraint of the robot, and connecting the states of the robot to obtain a motion element generated under the control vector, wherein i is 1, 2.
2. A robot motion primitive determination method according to claim 1, characterized in that the starting state vector is zs=[xs,yss]TWherein z represents a state vector, x and y represent displacements of the robot along the x direction and the y direction respectively, theta represents a heading angle of the robot, subscript s represents a start, and superscript T represents a matrix transposition.
3. The robot motion primitive determining method according to claim 2, wherein the determining the translational motion velocity section and the yaw velocity section reached by the robot within the time of the motion primitive generation, and the equally dividing the translational motion velocity section into m parts and the equally dividing the yaw velocity section into N parts according to the number N of the motion primitive generation comprises:
determining the translation motion speed interval [ v ] reached by the robot in the time t generated by the motion primitives-at,vs+at]And yaw angular velocity interval [ omega ]s-αt,ωs+αt]Wherein v represents the speed of the translational motion of the robot, ω represents the yaw rate of the robot, and a and α represent the acceleration at the translational motion of the robot and the yaw acceleration at the yaw motion, respectively;
according to the number N of the robot motion elements, equally dividing a translational motion speed interval into m-1 parts, equally dividing a yaw angle speed interval into N-1 parts, and satisfying the following conditions: n — m × N, then: translational motion speed corresponding to endpoint i
Figure FDA0002254300370000011
Yaw rate corresponding to endpoint j
Figure FDA0002254300370000012
4. The method for determining robot motion primitives according to claim 3, wherein the step of taking the translational motion velocity corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to form a control vector, determining the state of the robot according to the determined initial state vector of the robot and the kinematic constraint of the robot, and connecting the states of the robot to obtain the motion primitives generated under the control vector comprises the steps of:
within the time of motion primitive generation, with u ═ vij]TThe constant speed is used as a control vector, and the robot is determined at the sampling time t according to the determined initial state vector of the robot and the kinematic constraint of the robotfThe time-dependent state is:
Figure FDA0002254300370000021
wherein u represents a control vector, (x)f,yff) Indicating that the robot is at the sampling time tfState vector of (x)f,yff) Is a discrete point on the motion primitive, tfAny sampling moment in time t generated by a motion primitive is designated;
connecting state vectors at different sampling moments in the time of generating the motion primitive to obtain a control vector u ═ vij]TThe next generated motion primitive.
5. The robot motion primitive determining method according to claim 4, wherein the state vectors at different sampling time are connected in the time of motion primitive generation to obtain a control vector u ═ vij]TThe following generated motion primitives include:
in the time of generating the motion primitive, if the termination state vector is not at the grid center of the environment map, translating the termination state vector to the grid center by adopting a rounding method;
connecting the state vectors at different sampling times to obtain a control vector u ═ vij]TThe next generated motion primitive.
6. A robot motion primitive determination apparatus, comprising:
the determining module is used for determining an initial state vector of the robot;
the robot comprises an equally dividing module, a judging module and a judging module, wherein the equally dividing module is used for determining a translational motion speed interval and a yaw angular velocity interval which are reached by the robot in the time of generating motion primitives, and equally dividing the translational motion speed interval into m-1 parts and equally dividing the yaw angular velocity interval into N-1 parts according to the generated number N of the motion primitives, wherein N is m x N;
and the generating module is used for taking the translational motion speed corresponding to the endpoint i and the yaw velocity corresponding to the endpoint j to form a control vector, determining the state of the robot according to the determined initial state vector of the robot and the kinematic constraint of the robot, and connecting the states of the robot to obtain a motion element generated under the control vector, wherein i is 1,2, and.
7. Robot motion primitive determination device according to claim 6, characterized in that the starting state vector is zs=[xs,yss]TWherein z represents a state vector, x and y represent displacements of the robot along the x direction and the y direction respectively, theta represents a heading angle of the robot, subscript s represents a start, and superscript T represents a matrix transposition.
8. The robot motion primitive determination device of claim 7, wherein the aliquoting module comprises:
a first determining unit for determining the translational motion speed interval [ v ] reached by the robot in the time t of the generation of the motion primitives-at,vs+at]And yaw angular velocity interval [ omega ]s-αt,ωs+αt]Wherein v represents the speed of the translational motion of the robot, ω represents the yaw rate of the robot, and a and α represent the acceleration at the translational motion of the robot and the yaw acceleration at the yaw motion, respectively;
and the equally dividing unit is used for equally dividing the translational motion speed interval into m-1 parts and the yaw velocity interval into N-1 parts according to the generated number N of the robot motion elements, and meets the following requirements: n — m × N, then: translational motion speed corresponding to endpoint i
Figure FDA0002254300370000031
Yaw rate corresponding to endpoint j
Figure FDA0002254300370000032
9. The robot motion primitive determination device of claim 8, wherein the generation module comprises:
a second determination unit for generating a motion primitive with u ═ v in timeij]TThe constant speed is used as a control vector, and the robot is determined at the sampling time t according to the determined initial state vector of the robot and the kinematic constraint of the robotfThe time-dependent state is:
Figure FDA0002254300370000033
wherein u represents a control vector, (x)f,yff) Indicating that the robot is at the sampling time tfState vector of (x)f,yff) Is a discrete point on the motion primitive, tfAny sampling moment in time t generated by a motion primitive is designated;
a generating unit for connecting the state vectors at different sampling time points within the time of generating the motion primitive to obtain a control vector u ═ vij]TThe next generated motion primitive.
10. The robot motion primitive determination device according to claim 9, wherein the generation unit is configured to translate the termination state vector to the grid center by a rounding method if the termination state vector is not at the grid center of the environment map during the time of motion primitive generation; also for concatenating the state vectors at different sampling instants, resulting in a control vector u ═ vij]TThe next generated motion primitive.
CN201911046662.1A 2019-10-30 2019-10-30 Robot motion primitive determining method and device Pending CN110727271A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911046662.1A CN110727271A (en) 2019-10-30 2019-10-30 Robot motion primitive determining method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911046662.1A CN110727271A (en) 2019-10-30 2019-10-30 Robot motion primitive determining method and device

Publications (1)

Publication Number Publication Date
CN110727271A true CN110727271A (en) 2020-01-24

Family

ID=69222536

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911046662.1A Pending CN110727271A (en) 2019-10-30 2019-10-30 Robot motion primitive determining method and device

Country Status (1)

Country Link
CN (1) CN110727271A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112180726A (en) * 2020-09-29 2021-01-05 北京航空航天大学 Spacecraft relative motion trajectory planning method based on meta-learning
CN115495882A (en) * 2022-08-22 2022-12-20 北京科技大学 Method and device for constructing robot motion primitive library under non-flat terrain

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106774389A (en) * 2016-12-09 2017-05-31 武汉科技大学 A kind of four rotor wing unmanned aerial vehicles electricity tower method for inspecting based on motor learning
CN107357293A (en) * 2017-07-31 2017-11-17 上海应用技术大学 Method for planning path for mobile robot and system
CN107894766A (en) * 2017-10-18 2018-04-10 哈尔滨工业大学深圳研究生院 Paths planning method for interior decoration robot
CN108919794A (en) * 2018-06-01 2018-11-30 广州视源电子科技股份有限公司 A kind of global path planning method and its device of two-wheel differential type mobile robot
CN109799822A (en) * 2019-01-30 2019-05-24 中国石油大学(华东) Mobile robot global smooth paths planing method
CN109795502A (en) * 2018-09-27 2019-05-24 吉林大学 Intelligent electric automobile path trace model predictive control method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106774389A (en) * 2016-12-09 2017-05-31 武汉科技大学 A kind of four rotor wing unmanned aerial vehicles electricity tower method for inspecting based on motor learning
CN107357293A (en) * 2017-07-31 2017-11-17 上海应用技术大学 Method for planning path for mobile robot and system
CN107894766A (en) * 2017-10-18 2018-04-10 哈尔滨工业大学深圳研究生院 Paths planning method for interior decoration robot
CN108919794A (en) * 2018-06-01 2018-11-30 广州视源电子科技股份有限公司 A kind of global path planning method and its device of two-wheel differential type mobile robot
CN109795502A (en) * 2018-09-27 2019-05-24 吉林大学 Intelligent electric automobile path trace model predictive control method
CN109799822A (en) * 2019-01-30 2019-05-24 中国石油大学(华东) Mobile robot global smooth paths planing method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
RADHE SHYAM SHARMA等: "DMP Based Trajectory Tracking for a Nonholonomic Mobile Robot With Automatic Goal Adaptation and Obstacle Avoidance", 《2019 INTERNATIONAL CONFERENCE ON ROBOTICS AND AUTOMATION》 *
姜明浩等: "基于动态运动基元的移动机器人路径规划", 《高技术通讯》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112180726A (en) * 2020-09-29 2021-01-05 北京航空航天大学 Spacecraft relative motion trajectory planning method based on meta-learning
CN115495882A (en) * 2022-08-22 2022-12-20 北京科技大学 Method and device for constructing robot motion primitive library under non-flat terrain
CN115495882B (en) * 2022-08-22 2024-02-27 北京科技大学 Method and device for constructing robot motion primitive library under uneven terrain

Similar Documents

Publication Publication Date Title
CN108656112B (en) Mechanical arm zero-force control experiment system for direct teaching
CN103792885B (en) A kind of numerical controlled bending of pipe machining simulation method and device
CN107263466B (en) Base undisturbed control method of space robot based on quadratic programming problem
CN111427370A (en) Sparse pose adjustment-based Gmapping mapping method for mobile robot
CN107024863A (en) A kind of UUV Trajectory Tracking Control methods for avoiding differential from exploding
Wang et al. Precise discrete-time steering control for robotic fish based on data-assisted technique and super-twisting-like algorithm
Li et al. Localization and navigation for indoor mobile robot based on ROS
CN103878770A (en) Space robot visual delay error compensation method based on speed estimation
CN110561421B (en) Mechanical arm indirect dragging demonstration method and device
CN104076688A (en) Master-slave type cooperative control method for autonomous underwater vehicles
CN110727271A (en) Robot motion primitive determining method and device
CN110561419B (en) Arm-shaped line constraint flexible robot track planning method and device
CN110116409A (en) A kind of four-way remote operating bilateral control method based on disturbance observer
CN106896817A (en) A kind of many AUV formation control methods based on viscous damping mode
CN112155940A (en) Rehabilitation motion control method, system, equipment and medium based on rehabilitation robot
CN116429116A (en) Robot positioning method and equipment
CN114347018A (en) Mechanical arm disturbance compensation method based on wavelet neural network
CN114211492B (en) Optimal track planning method of multi-degree-of-freedom mechanical arm based on model
CN115122325A (en) Robust visual servo control method for anthropomorphic manipulator with view field constraint
Zhang et al. Research on obstacle avoidance algorithm of multi-uav consistent formation based on improved dynamic window approach
Huang et al. Immersive virtual simulation system design for the guidance, navigation and control of unmanned surface vehicles
CN115828463B (en) Heading machine pose estimation method and system based on digital twin six-dimensional model
CN103170976B (en) Two-freedom-degree robot motion control method
Gu et al. Dexterous obstacle-avoidance motion control of Rope Driven Snake Manipulator based on the bionic path following
CN114839990A (en) Cluster robot experiment platform

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200124