CN107643752B - Omnidirectional mobile robot path planning algorithm based on pedestrian trajectory prediction - Google Patents

Omnidirectional mobile robot path planning algorithm based on pedestrian trajectory prediction Download PDF

Info

Publication number
CN107643752B
CN107643752B CN201710322289.2A CN201710322289A CN107643752B CN 107643752 B CN107643752 B CN 107643752B CN 201710322289 A CN201710322289 A CN 201710322289A CN 107643752 B CN107643752 B CN 107643752B
Authority
CN
China
Prior art keywords
robot
pedestrian
path planning
point
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.)
Active
Application number
CN201710322289.2A
Other languages
Chinese (zh)
Other versions
CN107643752A (en
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.)
Qingyan Huayu Intelligent Robot Tianjin Co ltd
Original Assignee
Qingyan Huayu Intelligent Robot Tianjin 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 Qingyan Huayu Intelligent Robot Tianjin Co ltd filed Critical Qingyan Huayu Intelligent Robot Tianjin Co ltd
Priority to CN201710322289.2A priority Critical patent/CN107643752B/en
Publication of CN107643752A publication Critical patent/CN107643752A/en
Application granted granted Critical
Publication of CN107643752B publication Critical patent/CN107643752B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Manipulator (AREA)

Abstract

The invention provides an omnidirectional mobile robot path planning algorithm based on pedestrian track prediction, which comprises the following steps: s1, performing kinematic modeling by using a robot chassis and Mecanum wheels installed on the chassis; s2, planning a path by the robot according to the current coordinate and the target position coordinate; s3, the robot predicts the behavior track of the pedestrian in the visual field through the laser radar module, and calculates the possible range of the pedestrian within a certain time; and S4, planning the pedestrian avoiding path by the robot according to the possible range of the pedestrian. The omnidirectional mobile robot path planning algorithm based on pedestrian track prediction can enable the robot to move more flexibly; and tracking the pedestrians in the environment in real time through the laser radar module, and prejudging the positions of the pedestrians at the next moment according to the tracking result of the pedestrians.

Description

Omnidirectional mobile robot path planning algorithm based on pedestrian trajectory prediction
Technical Field
The invention belongs to the field of robot path planning, and particularly relates to an omnidirectional mobile robot path planning algorithm based on pedestrian trajectory prediction.
Background
Since the advent of the traditional industrial robots, robotics has made significant progress in as little as 60 years. Nowadays, robots are no longer limited to industrial applications, and the research focus and development direction thereof are more oriented towards daily life applications, such as catering, welcoming, entertainment, companions, and the like. Among many forms of mobile robots, omni-directional wheeled mobile robots are also receiving attention for their flexible kinematic performance. With the rapid development of modern sensing technology, computer technology and artificial intelligence technology, path planning of wheeled mobile robots and obstacle avoidance technology in the motion process also become research hotspots of robot technology.
The path planning problem can be divided into global path planning where the environmental information is completely known and local path planning where the environmental information is partially or even completely unknown. When the environmental information is partially or completely unknown, the robot senses the external environmental information through a sensing device mounted on the robot, and then reaches the destination in an approximate optimal path without collision. In the traditional local path planning, the robot is usually only enabled to avoid external obstacles in real time in the walking process, and the movement of external moving obstacles such as pedestrians is not pre-judged.
Disclosure of Invention
In view of the above, the present invention aims to provide a pedestrian trajectory prediction-based path planning algorithm for an omnidirectional mobile robot, so as to solve the problem that the current mobile robot has a single function and cannot make a motion prediction.
In order to achieve the purpose, the technical scheme of the invention is realized as follows:
an omnidirectional mobile robot path planning algorithm based on pedestrian track prediction comprises the following steps:
s1, performing kinematic modeling by using a robot chassis and Mecanum wheels installed on the chassis;
s2, planning a path by the robot according to the current coordinate and the target position coordinate;
s3, the robot predicts the behavior track of the pedestrian in the visual field through the laser radar module, and calculates the possible range of the pedestrian within a certain time;
and S4, planning the pedestrian avoiding path by the robot according to the possible range of the pedestrian.
Further, in step S1, the robot chassis is equipped with three mecanum wheels, and the mecanum wheels are left-handed mecanum wheels.
Further, in step S1, the robot kinematics modeling method specifically includes:
s101, establishing a robot body coordinate system on a robot chassis, and establishing a wheel coordinate system on wheels of the robot;
s102, setting the X axis of the robot body coordinate system to each wheel coordinate system XiIncluded angle of axis being gammaiThe included angle between the rolling shaft of each wheel and the wheel shaft is alpha, and the connecting line between the central point O and the central point of each wheel is LiX-axis to L in robot body coordinate systemiIs thetai
S103, setting the rotation angular velocities of the three wheel shafts to
Figure BDA0001290166480000021
The central movement speed of the robot is [ v ]x,vy,w]TThe radius of the outer wheel is R, and the radius of the roller is R;
s104, the inverse kinematics formula can be expressed as:
Figure BDA0001290166480000022
wherein,
Figure BDA0001290166480000023
the angle between the rolling axis of each wheel and the axle
Figure BDA0001290166480000024
Further, in step S2, the robot path planning implements path planning according to a path planning algorithm, the robot path planning includes global path planning and local path planning, the global path planning performs path planning according to a global grid map, current coordinates, and target position coordinates, an optimal global path from a start point to a destination point is correctly obtained, and after the optimal global path is determined, local path planning is started.
Furthermore, the local path planning adopts a DWA algorithm, the DWA algorithm collects multiple groups of speed data in a speed space, simulates tracks of the robot at the speeds within a certain time, then performs scoring evaluation on the tracks according to a scoring function, selects an optimal track, and generates a robot speed control command.
Further, the specific method for local path planning is as follows:
s201, because the robot moves in all directions, the position relation between a starting point and an end point needs to be analyzed first, a space quadrant of a speed space is determined, and the current moving speed of the robot is set to be (v)cx,vcy,wc):
S202, setting the robot to have kinematic speed limit, and setting an x-direction speed limit, a y-direction speed limit and a w-rotation speed limit as follows respectively: v. ofx_min、vx_max、vy_min、vy_max、wmin、wmax
Vs={vx∈[vxmin,vxmax],vy∈[vymin,vymax],w∈[wmin,wmax]}
S203, influenced by the performance and the motion stability of the motor, the robot has the maximum acceleration in the motion, so that a cubic area V exists in the forward simulation period (delta t) of the track of the robotdThe speed in the cubic region is the speed that the robot can actually reach:
Figure BDA0001290166480000031
s204, the robot must stop before colliding with the obstacle due to safety limitation, so that an anti-collision speed area V can be obtaineda
Figure BDA0001290166480000032
Wherein docc(vx,vyW) is the velocity (v)x,vyAnd w) the closest distance to the obstacle on the simulated trajectory.
The search space V of the DWA algorithm improved by the steps S202, S203 and S204 can be obtainedr=Vs∩Vd∩Va
S205, after obtaining the velocity space, in order to optimize the local path, each velocity (v) in Vr needs to be scored according to a certain scoring standardx,vyW) scoring;
the scoring function is as follows:
score=α*dpath+β*dgoal+γ*docc
wherein, alpha, beta and gamma are respectively a path fitting weight, a distance target point position weight and a distance weight with an obstacle; doccIs velocity (v)x,vyW) maximum obstacle distance value on trajectory, dgoalIs velocity (v)x,vyW) distance of the trajectory end point to the global target point, dpathIs the degree of conformance of the trajectory to the global path.
Further, in step S3, the pedestrian trajectory prediction is implemented by a behavior trajectory prediction algorithm, when the system finds a person, the system outputs a relative coordinate between the pedestrian and the lidar module, and the relative coordinate is located in the robot body coordinate system, which is specifically implemented as follows:
let the coordinates of the current and the previous three pedestrians in the global coordinate system be pk=(xk,yk),pk-1=(xk-1,yk-1),pk-2=(xk-2,yk-2),pk-3=(xk-3,yk-3) The time when the system acquires the coordinates is t4,t3,t2,t1. From these four sets of data, the average velocity of the robot over this period of time can be calculated:
Figure BDA0001290166480000041
Figure BDA0001290166480000042
let the predicted time be
Figure BDA0001290166480000043
Then
Figure BDA0001290166480000044
Position coordinates p at which a person may be after a timepre_k=(xpre_k,ypre_k) Can be expressed as:
Figure BDA0001290166480000045
Figure BDA0001290166480000046
with ppre_kAs the center of a circle, drawing a circle with radius r, and within the circle, predicting time
Figure BDA0001290166480000056
The range S in which the person is likely to bek
When a person is lost or retraced, the loop needs to be restarted and all variables initialized.
Further, in step S4, the pedestrian-avoiding path planning is implemented by an algorithm, which specifically includes:
s301, according to a local path planning algorithm, the robot plans an optimal track of future delta t time, and sets a last point p of the trackiThe coordinate is (x)i,yi). If the pedestrian tracking result exists, the coordinate point p of the pedestrian is assumed to be predicted at the momentpreThe coordinate is (x)pre,ypre)。
S302, if
Figure BDA0001290166480000051
The robot path planning point is not in the radiation circle of the pedestrian prediction point, the robot is relatively safe, and the robot can plan to walk according to the optimal track;
s303, if
Figure BDA0001290166480000052
The robot path planning point is located inside or at the edge of the radiation circle, and the next moment may touch the pedestrian, and then the path planning for fast escaping from the radiation circle is performed.
Further, when the robot is located in or at the edge of the radiation ring of the pedestrian prediction point, the moving linear speed of the robot escapes from the radiation ring at the fastest speed along the tangential direction of the circle, and the specific method is as follows:
let the coordinate of the last point of the predicted point of the local path of the robot be pi=(xi,yi) Pedestrian prediction point ppreThe coordinate is (x)pre,ypre) At this time, the distance between two points
Figure BDA0001290166480000053
Let the maximum linear velocity of robot motion be vmaxMinimum linear velocity of vminAnd the distance between the two points and the linear velocity have a linear relation:
therefore, the escape linear velocity of the robot at this time should be:
Figure BDA0001290166480000054
the speed direction is divided into two parts which are respectively v1And v2According to a rotation matrix
Figure BDA0001290166480000055
Unit vectors for two speed directions are obtained:
Figure BDA0001290166480000061
according to the formula A and the formula B, the speeds in two tangential directions are obtained
Figure BDA0001290166480000062
And
Figure BDA0001290166480000063
according to a local path planning algorithm, respectively
Figure BDA0001290166480000064
And
Figure BDA0001290166480000065
and taking positive and negative ranges of 30 degrees as the center, and taking the positive and negative ranges as the planning speed to plan the local path of the robot.
Compared with the prior art, the omni-directional mobile robot path planning algorithm based on pedestrian track prediction has the following advantages:
the omnidirectional mobile robot path planning algorithm based on pedestrian track prediction can enable the robot to move more flexibly; tracking the pedestrians in the environment in real time through a laser radar module, and pre-judging the positions of the pedestrians at the next moment according to the tracking result of the pedestrians; and according to the pedestrian track pre-judgment result, adopting a speed combination near the direction of the predicted point connecting line perpendicular line to re-plan the path so as to achieve the purpose of avoiding pedestrians in advance.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate an embodiment of the invention and, together with the description, serve to explain the invention and not to limit the invention. In the drawings:
fig. 1 is a schematic diagram illustrating a kinematics modeling principle of an omnidirectional mobile robot according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of the improved DWA algorithm according to an embodiment of the present invention;
fig. 3 is a schematic diagram illustrating a principle of avoiding a pedestrian prediction point according to an embodiment of the present invention.
Detailed Description
It should be noted that the embodiments and features of the embodiments may be combined with each other without conflict.
In the description of the present invention, it is to be understood that the terms "center", "longitudinal", "lateral", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", and the like, indicate orientations or positional relationships based on those shown in the drawings, and are used only for convenience in describing the present invention and for simplicity in description, and do not indicate or imply that the referenced devices or elements must have a particular orientation, be constructed and operated in a particular orientation, and thus, are not to be construed as limiting the present invention. Furthermore, the terms "first", "second", etc. are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first," "second," etc. may explicitly or implicitly include one or more of that feature. In the description of the present invention, "a plurality" means two or more unless otherwise specified.
In the description of the present invention, it should be noted that, unless otherwise explicitly specified or limited, the terms "mounted," "connected," and "connected" are to be construed broadly, e.g., as meaning either a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meaning of the above terms in the present invention can be understood by those of ordinary skill in the art through specific situations.
The present invention will be described in detail below with reference to the embodiments with reference to the attached drawings.
The omnidirectional mobile robot path planning algorithm based on pedestrian track prediction comprises the following steps:
s1, performing kinematic modeling by using a robot chassis and Mecanum wheels installed on the chassis;
s2, planning a path by the robot according to the current coordinate and the target position coordinate;
s3, the robot predicts the behavior track of the pedestrian in the visual field through the laser radar module, and calculates the possible range of the pedestrian within a certain time;
and S4, planning the pedestrian avoiding path by the robot according to the possible range of the pedestrian.
Wherein the robot comprises three left-handed mecanum wheels mounted to a bottom of the robot.
As shown in fig. 1, the robot kinematics modeling method specifically includes the following steps:
s101, establishing a robot body coordinate system on a robot chassis, and establishing a wheel coordinate system on wheels of the robot;
s102, setting the X axis of the robot body coordinate system to each wheel coordinate system XiIncluded angle of axis being gammaiThe included angle between the rolling shaft of each wheel and the wheel shaft is alpha, and the connecting line between the central point O and the central point of each wheel is LiX-axis to L in robot body coordinate systemiIs thetai
S103, setting the rotation angular velocities of the three wheel shafts to
Figure BDA0001290166480000081
The central movement speed of the robot is [ v ]x,vy,w]TThe radius of the outer wheel is R, and the radius of the roller is R;
s104, the inverse kinematics formula can be expressed as:
Figure BDA0001290166480000082
wherein,
Figure BDA0001290166480000083
the angle between the rolling axis of each wheel and the axle
Figure BDA0001290166480000084
When Mecanum wheel and groundIntersecting roller and LiWhen vertical, the matrix a is always invertible, i.e. the positive kinematics always has a solution.
The robot path planning algorithm comprises global path planning and local path planning, wherein the global path planning carries out path planning according to a global grid map, current coordinates and target position coordinates, an optimal global path from a starting point to a terminal point is correctly obtained, and the local path planning is started after the optimal global path is determined.
The local path planning adopts a DWA algorithm, the DWA algorithm collects multiple groups of speed data in a speed space, simulates tracks of the robot at the speeds within a certain time, then performs scoring evaluation on the tracks according to a scoring function, selects an optimal track, and generates a robot speed control command.
As shown in fig. 2, the specific method of the local path planning is as follows:
a conventional DWA (Dynamic Window Approach) algorithm is generally applied to a two-wheeled differential mobile robot, and only forward velocity vx and rotation angular velocity w need to be considered, so that a velocity Window is a rectangle in a plane. However, the movement characteristic of the omnidirectional mobile robot is omnidirectional movement in a plane, and in order to ensure the flexibility of the movement, the speed of transverse movement needs to be planned at the same time, that is, the window pane is not a plane rectangle any more, but a space cuboid. Therefore, the traditional dynamic window lattice method no longer meets the flexibility requirement of the omnidirectional mobile robot, and needs to be improved.
S201, because the robot moves in all directions, the position relation between a starting point and an end point needs to be analyzed first, a space quadrant of a speed space is determined, and the current moving speed of the robot is set to be (v)cx,vcy,wc):
S202, setting the robot to have kinematic speed limit, and setting an x-direction speed limit, a y-direction speed limit and a w-rotation speed limit as follows respectively: v. ofx_min、vx_max、vy_min、vy_max、wmin、wmax
Vs={vx∈[vxmin,vxmax],vy∈[vymin,vymax],w∈[wmin,wmax]}
S203, influenced by the performance and the motion stability of the motor, the robot has the maximum acceleration in the motion, so that a cubic area V exists in the forward simulation period (delta t) of the track of the robotdThe speed in the cubic region is the speed that the robot can actually reach:
Figure BDA0001290166480000091
s204, the robot must stop before colliding with the obstacle due to safety limitation, so that an anti-collision speed area V can be obtaineda
Figure BDA0001290166480000101
Wherein docc(vx,vyW) is the velocity (v)x,vyAnd w) the closest distance to the obstacle on the simulated trajectory.
The search space V of the DWA algorithm improved by the steps S202, S203 and S204 can be obtainedr=Vs∩Vd∩Va
After obtaining the velocity space, it is necessary to perform the scoring for each velocity (v) in Vr based on a certain scoring criterion in order to optimize the local path (avoid the obstacle, fit the predetermined optimal global path as much as possible, and minimize the position error from the final point)x,vyW) scoring;
the scoring function is as follows:
score=α*dpath+β*dgoal+γ*docc
wherein, alpha, beta and gamma are respectively a path fitting weight, a distance target point position weight and a distance weight with an obstacle; doccIs velocity (v)x,vyW) maximum obstacle distance value on trajectory, dgoalIs velocity (v)x,vyW) distance of the trajectory end point to the global target point, dpathIs the degree of conformance of the trajectory to the global path.
When the system finds a person, the pedestrian track prediction algorithm outputs a relative coordinate between the pedestrian and the laser radar module, and the relative coordinate value is located in the robot body coordinate system, and the specific method comprises the following steps:
let the coordinates of the current and the previous three pedestrians in the global coordinate system be pk=(xk,yk),pk-1=(xk-1,yk-1),pk-2=(xk-2,yk-2),pk-3=(xk-3,yk-3) The time when the system acquires the coordinates is t4,t3,t2,t1. From these four sets of data, the average velocity of the robot over this period of time can be calculated:
Figure BDA0001290166480000102
Figure BDA0001290166480000103
let the predicted time be
Figure BDA0001290166480000104
Then
Figure BDA0001290166480000105
Position coordinates p at which a person may be after a timepre_k=(xpre_k,ypre_k) Can be expressed as:
Figure BDA0001290166480000116
Figure BDA0001290166480000117
with ppre_kAs the center of a circle, drawing a circle with radius r, and within the circle, predicting time
Figure BDA0001290166480000118
The range S in which the person is likely to bek
When a person is lost or retraced, the loop needs to be restarted and all variables initialized.
The pedestrian-avoiding path planning algorithm specifically comprises the following steps:
s301, according to a local path planning algorithm, the robot plans an optimal track of future delta t time, and sets a last point p of the trackiThe coordinate is (x)i,yi). If the pedestrian tracking result exists, the coordinate point p of the pedestrian is assumed to be predicted at the momentpreThe coordinate is (x)pre,ypre)。
S302, if
Figure BDA0001290166480000111
The robot path planning point is not in the radiation circle of the pedestrian prediction point, the robot is relatively safe, and the robot can plan to walk according to the optimal track;
s303, if
Figure BDA0001290166480000112
The robot path planning point is located inside or at the edge of the radiation circle, and the next moment may touch the pedestrian, and then the path planning for fast escaping from the radiation circle is performed.
As shown in fig. 3, when the robot is located in or at the edge of the radiation circle of the pedestrian prediction point, the moving linear speed of the robot can fast escape from the radiation circle along the tangential direction of the circle, and the specific method is as follows:
let the coordinate of the last point of the predicted point of the local path of the robot be pi=(xi,yi) Pedestrian prediction point ppreThe coordinate is (x)pre,ypre) At this time, the distance between two points
Figure BDA0001290166480000113
Let the maximum linear velocity of robot motion be vmaxMinimum linear velocity of vminAnd the distance between the two points and the linear velocity have a linear relation:
therefore, the escape linear velocity of the robot at this time should be:
Figure BDA0001290166480000114
the speed direction is divided into two parts which are respectively v1And v2According to a rotation matrix
Figure BDA0001290166480000115
Unit vectors for two velocity directions can be obtained:
Figure BDA0001290166480000121
based on the formulas A and B, two tangential speeds can be obtained
Figure BDA0001290166480000122
And
Figure BDA0001290166480000123
according to a local path planning algorithm, respectively
Figure BDA0001290166480000124
And
Figure BDA0001290166480000125
and taking positive and negative ranges of 30 degrees as the center, and taking the positive and negative ranges as the planning speed to plan the local path of the robot.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims (6)

1. The omnidirectional mobile robot path planning algorithm based on pedestrian track prediction is characterized by comprising the following steps:
s1, performing kinematic modeling by using a robot chassis and Mecanum wheels installed on the chassis;
s2, planning a path by the robot according to the current coordinate and the target position coordinate;
s3, the robot predicts the behavior track of the pedestrian in the visual field through the laser radar module, and calculates the possible range of the pedestrian within a certain time;
s4, planning a pedestrian avoiding path by the robot according to the possible range of the pedestrian;
in step S3, the pedestrian trajectory prediction is implemented by a behavior trajectory prediction algorithm, and when the system finds a person, the system outputs a relative coordinate between the pedestrian and the lidar module, and the relative coordinate is located in the robot body coordinate system, and the specific method is as follows:
let the coordinates of the current and the previous three pedestrians in the global coordinate system be pk=(xk,yk),pk-1=(xk-1,yk-1),pk-2=(xk-2,yk-2),pk-3=(xk-3,yk-3) The time when the system acquires the coordinates is t4,t3,t2,t1(ii) a From these four sets of data, the average velocity of the robot over this period of time can be calculated:
Figure FDA0002743843950000011
Figure FDA0002743843950000012
let the predicted time be
Figure FDA0002743843950000015
Then
Figure FDA0002743843950000016
Position coordinates p at which a person may be after a timepre_k=(xpre_k,ypre_k) Can be expressed as:
Figure FDA0002743843950000013
Figure FDA0002743843950000014
with ppre_kAs the center of a circle, drawing a circle with radius r, and within the circle, predicting time
Figure FDA0002743843950000017
The range S in which the person is likely to bek
When a person is lost or tracked again, the circulation needs to be restarted, and all variables are initialized;
in step S4, the pedestrian avoidance path planning is implemented by an algorithm, which specifically includes:
s301, according to a local path planning algorithm, the robot plans an optimal track of future delta t time, and sets a last point p of the trackiThe coordinate is (x)i,yi) (ii) a If the pedestrian tracking result exists, the coordinate point p of the pedestrian is assumed to be predicted at the momentpreThe coordinate is (x)pre,ypre);
S302, if
Figure FDA0002743843950000021
The robot path planning point is not in the radiation circle of the pedestrian prediction point, the robot is relatively safe, and the robot can plan to walk according to the optimal track;
s303, if
Figure FDA0002743843950000022
Machine for testingThe robot path planning point is positioned in or at the edge of the radiation ring, and the next moment is likely to touch the pedestrian, so that path planning for rapidly escaping from the radiation ring is performed;
when the robot is positioned in or at the edge of a radiation ring of a pedestrian prediction point, the moving linear speed of the robot escapes from the radiation ring fastest along the tangential direction of a circle, and the specific method is as follows:
let the coordinate of the last point of the predicted point of the local path of the robot be pi=(xi,yi) Pedestrian prediction point ppreThe coordinate is (x)pre,ypre) At this time, the distance between two points
Figure FDA0002743843950000023
Let the maximum linear velocity of robot motion be vmaxMinimum linear velocity of vminAnd the distance between the two points and the linear velocity have a linear relation:
therefore, the escape linear velocity of the robot at this time should be:
Figure FDA0002743843950000024
the speed direction is divided into two parts which are respectively v1And v2According to a rotation matrix
Figure FDA0002743843950000025
Unit vectors for two speed directions are obtained:
Figure FDA0002743843950000026
according to the formula A and the formula B, the speeds in two tangential directions are obtained
Figure FDA0002743843950000027
And
Figure FDA0002743843950000028
according to a local path planning algorithm, respectively
Figure FDA0002743843950000029
And
Figure FDA00027438439500000210
and taking positive and negative ranges of 30 degrees as the center, and taking the positive and negative ranges as the planning speed to plan the local path of the robot.
2. The pedestrian trajectory prediction-based omnidirectional mobile robot path planning algorithm of claim 1, wherein: in step S1, the robot chassis is provided with three mecanum wheels, which are left-handed mecanum wheels.
3. The pedestrian trajectory prediction-based omnidirectional mobile robot path planning algorithm according to claim 1 or 2, wherein: in step S1, the robot kinematics modeling method specifically includes:
s101, establishing a robot body coordinate system on a robot chassis, and establishing a wheel coordinate system on wheels of the robot;
s102, setting the X axis of the robot body coordinate system to each wheel coordinate system XiIncluded angle of axis being gammaiThe included angle between the rolling shaft of each wheel and the wheel shaft is alpha, and the connecting line between the central point O and the central point of each wheel is LiX-axis to L in robot body coordinate systemiIs thetai
S103, setting the rotation angular velocities of the three wheel shafts to
Figure FDA0002743843950000031
The central movement speed of the robot is [ v ]x,vy,w]ΤThe radius of the outer wheel is R, and the radius of the roller is R;
s104, the inverse kinematics formula can be expressed as:
Figure FDA0002743843950000032
wherein,
Figure FDA0002743843950000033
the angle between the rolling axis of each wheel and the axle
Figure FDA0002743843950000034
4. The pedestrian trajectory prediction-based omnidirectional mobile robot path planning algorithm of claim 2, wherein: in step S2, the robot path planning implements path planning according to a path planning algorithm, the robot path planning includes global path planning and local path planning, the global path planning performs path planning according to a global grid map, current coordinates, and target position coordinates, an optimal global path from a start point to an end point is correctly obtained, and after the optimal global path is determined, local path planning is started.
5. The pedestrian trajectory prediction-based omnidirectional mobile robot path planning algorithm of claim 4, wherein: the local path planning adopts a DWA algorithm, the DWA algorithm collects multiple groups of speed data in a speed space, simulates tracks of the robot at the speeds within a certain time, then carries out scoring evaluation on the tracks according to a scoring function, selects an optimal track and generates a robot speed control command.
6. The pedestrian trajectory prediction-based omnidirectional mobile robot path planning algorithm of claim 5, wherein: the specific method for local path planning is as follows:
s201, because the robot moves in all directions, the position relation between the starting point and the end point needs to be analyzed first, and the space image of the speed space is determinedLimit, set the current moving speed of the robot as (v)cx,vcy,wc):
S202, setting the robot to have kinematic speed limit, and setting an x-direction speed limit, a y-direction speed limit and a w-rotation speed limit as follows respectively: v. ofx_min、vx_max、vy_min、vy_max、wmin、wmax
Vs={vx∈[vxmin,vxmax],vy∈[vymin,vymax],w∈[wmin,wmax]}
S203, influenced by the performance and the motion stability of the motor, the robot has the maximum acceleration in the motion, so that a cubic area V exists in the forward simulation period (delta t) of the track of the robotdThe speed in the cubic region is the speed that the robot can actually reach:
Figure FDA0002743843950000041
s204, the robot must stop before colliding with the obstacle due to safety limitation, so that an anti-collision speed area V can be obtaineda
Figure FDA0002743843950000051
Wherein d isocc(vx,vyW) is the velocity (v)x,vyW) the closest distance to the obstacle on the simulated trajectory;
the search space V of the DWA algorithm improved by the steps S202, S203 and S204 can be obtainedr=Vs∩Vd∩Va
S205, after obtaining the velocity space, in order to optimize the local path, each velocity (v) in Vr needs to be scored according to a certain scoring standardx,vyW) scoring;
the scoring function is as follows:
score=α*dpath+β*dgoal+γ*docc
wherein, alpha, beta and gamma are respectively a path fitting weight, a distance target point position weight and a distance weight with an obstacle; doccIs velocity (v)x,vyW) maximum obstacle distance value on trajectory, dgoalIs velocity (v)x,vyW) distance of the trajectory end point to the global target point, dpathIs the degree of conformance of the trajectory to the global path.
CN201710322289.2A 2017-05-09 2017-05-09 Omnidirectional mobile robot path planning algorithm based on pedestrian trajectory prediction Active CN107643752B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710322289.2A CN107643752B (en) 2017-05-09 2017-05-09 Omnidirectional mobile robot path planning algorithm based on pedestrian trajectory prediction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710322289.2A CN107643752B (en) 2017-05-09 2017-05-09 Omnidirectional mobile robot path planning algorithm based on pedestrian trajectory prediction

Publications (2)

Publication Number Publication Date
CN107643752A CN107643752A (en) 2018-01-30
CN107643752B true CN107643752B (en) 2020-12-08

Family

ID=61110250

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710322289.2A Active CN107643752B (en) 2017-05-09 2017-05-09 Omnidirectional mobile robot path planning algorithm based on pedestrian trajectory prediction

Country Status (1)

Country Link
CN (1) CN107643752B (en)

Families Citing this family (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108445503B (en) * 2018-03-12 2021-09-14 吉林大学 Unmanned path planning algorithm based on fusion of laser radar and high-precision map
CN108681243B (en) * 2018-05-03 2020-12-18 南京理工大学 Robot trajectory tracking method
CN108961313B (en) * 2018-06-29 2021-06-29 大连民族大学 Overlooking pedestrian risk quantification method of two-dimensional world coordinate system
CN108983787B (en) * 2018-08-09 2021-09-10 北京智行者科技有限公司 Road driving method
CN109062227A (en) * 2018-09-25 2018-12-21 电子科技大学 A kind of more curve local paths planning methods based on dynamic window
CN111103875B (en) * 2018-10-26 2021-12-03 科沃斯机器人股份有限公司 Method, apparatus and storage medium for avoiding
CN109460029A (en) * 2018-11-29 2019-03-12 华南农业大学 Livestock and poultry cultivation place inspection mobile platform and its control method
CN109572857B (en) * 2018-12-26 2023-06-02 石家庄铁道大学 Mecanum wheel intelligent storage AGV and path planning method thereof
US11986964B2 (en) 2018-12-27 2024-05-21 Honda Motor Co., Ltd. Path determination device, robot, and path determination method
CN109668484B (en) * 2019-01-18 2023-05-02 北京瀚科科技集团有限公司 Target aircraft maneuvering flight control method and system for interaction of target aircraft and attack aircraft
CN109827586B (en) * 2019-02-20 2022-03-22 百度在线网络技术(北京)有限公司 Method, device and equipment for planning speed of automatic driving vehicle
CN109901575A (en) * 2019-02-20 2019-06-18 百度在线网络技术(北京)有限公司 Vehicle routing plan adjustment method, device, equipment and computer-readable medium
IL265495B (en) * 2019-03-19 2022-09-01 Rober Ohrenstein Traffic allowance method
CN110488811B (en) * 2019-07-22 2023-04-07 上海有个机器人有限公司 Method for predicting pedestrian track by robot based on social network model
CN112445211B (en) * 2019-08-15 2023-10-27 纳恩博(北京)科技有限公司 Robot control method and device, storage medium and electronic device
CN110531698B (en) * 2019-08-16 2022-03-11 达闼机器人有限公司 Mechanical arm path planning method and device and mechanical arm
CN111123947B (en) * 2019-12-31 2023-07-07 达闼机器人股份有限公司 Robot travel control method and device, electronic equipment, medium and robot
CN111288995B (en) * 2020-03-12 2022-05-13 深圳市人工智能与机器人研究院 Route planning method and route planning device of mobile robot and terminal equipment
CN111391828B (en) * 2020-03-13 2021-06-15 南京航空航天大学 Method for planning vehicle collision avoidance track under pedestrian crossing working condition
CN111429791B (en) * 2020-04-09 2022-11-18 浙江大华技术股份有限公司 Identity determination method, identity determination device, storage medium and electronic device
CN111308473A (en) * 2020-04-14 2020-06-19 安徽卫盾安全设备科技有限公司 Active three-dimensional through-wall radar equipment
CN111986224B (en) * 2020-08-05 2024-01-05 七海行(深圳)科技有限公司 Target behavior prediction tracking method and device
CN112035921B (en) * 2020-08-25 2021-06-29 中船文化科技(北京)有限公司 Shelter moving line planning method and device, electronic equipment and storage medium
CN113359715A (en) * 2021-05-26 2021-09-07 上海大学 Mecanum wheel-based mobile robot local trajectory planning method and system
CN113296519A (en) * 2021-05-26 2021-08-24 上海大学 Mecanum wheel-based mobile robot motion planning method and system
CN113515122B (en) * 2021-06-24 2021-12-28 深圳市普渡科技有限公司 Robot, robot welcome path planning method and device and readable storage medium
CN113344303B (en) * 2021-07-19 2023-05-23 安徽工程大学 Time window dynamic obstacle avoidance method for optimizing energy consumption of multi-mobile robot under three-dimensional terrain
CN113721633B (en) * 2021-09-09 2023-10-17 南京工业大学 Mobile robot path planning method based on pedestrian track prediction
CN113848893A (en) * 2021-09-14 2021-12-28 武汉联一合立技术有限公司 Robot navigation method, device, equipment and storage medium
CN114313878B (en) * 2021-11-19 2024-04-19 江苏科技大学 Kinematic modeling and path planning method for material transmission platform
CN113822253B (en) * 2021-11-24 2022-02-18 天津大学 Man-machine cooperation method and system
CN116013074A (en) * 2023-01-05 2023-04-25 北京清丰智行科技有限公司 Intelligent travel system based on car Lu Yun cooperation in park
CN115826588B (en) * 2023-02-14 2023-05-12 北京国安广传网络科技有限公司 Planning module for correcting path of health management robot based on LBS data
CN117091608B (en) * 2023-10-17 2023-12-29 山东怀思人工智能科技有限公司 Path planning method, equipment and medium for educational robot

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101359229A (en) * 2008-08-18 2009-02-04 浙江大学 Barrier-avoiding method for mobile robot based on moving estimation of barrier
CN103558856A (en) * 2013-11-21 2014-02-05 东南大学 Service mobile robot navigation method in dynamic environment
CN105955273A (en) * 2016-05-25 2016-09-21 速感科技(北京)有限公司 Indoor robot navigation system and method
CN106406320A (en) * 2016-11-29 2017-02-15 重庆重智机器人研究院有限公司 Robot path planning method and robot planning route

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101359229A (en) * 2008-08-18 2009-02-04 浙江大学 Barrier-avoiding method for mobile robot based on moving estimation of barrier
CN103558856A (en) * 2013-11-21 2014-02-05 东南大学 Service mobile robot navigation method in dynamic environment
CN105955273A (en) * 2016-05-25 2016-09-21 速感科技(北京)有限公司 Indoor robot navigation system and method
CN106406320A (en) * 2016-11-29 2017-02-15 重庆重智机器人研究院有限公司 Robot path planning method and robot planning route

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Bi-steerable robot navigation using a modified Dynamic Window Approach;Karima Rebai等;《2009 6th International Symposium on Mechatronics and its Applications》;IEEE;20090714;第1-6页 *
Human-centered X–Y–T space path planning for mobile robot in dynamic environments;Ippei Nishitani等;《Robotics and Autonomous Systems》;ELSEVIER;20150430;第66卷;第18-26页 *
基于麦克纳姆轮全向车设计与性能研究;朱艳杰;《中国优秀硕士学位论文全文数据库 信息科技辑》;中国学术期刊(光盘版)电子杂志社;20170115(第01期);第26-30、56页 *

Also Published As

Publication number Publication date
CN107643752A (en) 2018-01-30

Similar Documents

Publication Publication Date Title
CN107643752B (en) Omnidirectional mobile robot path planning algorithm based on pedestrian trajectory prediction
CN106406338B (en) Autonomous navigation device and method of omnidirectional mobile robot based on laser range finder
Zhu et al. A convex optimization approach to smooth trajectories for motion planning with car-like robots
Jian et al. Dynamic control barrier function-based model predictive control to safety-critical obstacle-avoidance of mobile robot
CN112183133B (en) Aruco code guidance-based mobile robot autonomous charging method
KR102210788B1 (en) Automated guided vehicle, system with a computer and an automated guided vehicle, method for operating an automated guided vehicle
WO2023116667A1 (en) Charging device, and method for controlling mechanical arm to charge
CN108124472B (en) Method and device for avoiding barrier and aircraft
KR101278708B1 (en) Underwater robot and Method for controlling the same
CN112762928B (en) ODOM and DM landmark combined mobile robot containing laser SLAM and navigation method
CN106774295A (en) A kind of distributed guided robot recharging system
US11712802B2 (en) Construction constrained motion primitives from robot maps
JPWO2017221859A1 (en) Mobile body, mobile body control system, mobile body control method, interface device, and program
CN113110521A (en) Mobile robot path planning control method, control device thereof and storage medium
Asadi et al. An integrated aerial and ground vehicle (UAV-UGV) system for automated data collection for indoor construction sites
CN117055556A (en) Multi-robot formation path planning method and device, electronic equipment and storage medium
CN116520879A (en) Unmanned aerial vehicle tracking method based on vision
JP2020064029A (en) Mobile body controller
Nasti et al. Obstacle avoidance during robot navigation in dynamic environment using fuzzy controller
CN112585616A (en) Method for predicting at least one future speed vector and/or future posture of a pedestrian
Zhu et al. A nonlinear model predictive control based control method to quadrotor landing on moving platform
US20220339786A1 (en) Image-based trajectory planning method and movement control method and mobile machine using the same
Dorj et al. Multi-sensor fusion based effective obstacle avoidance and path-following technology
Xiaodong et al. Obstacle avoidance for outdoor flight of a quadrotor based on computer vision
US20230040374A1 (en) Control device, control method, and storage medium

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