Generation of Achievable Three-Dimensional Trajectories for Autonomous Wheeled Vehicles via Tracking Differentiators

: Planning an achievable trajectory for a mobile robot usually consists of two steps: (i) ﬁnding a path in the form of a sequence of discrete waypoints and (ii) transforming this sequence into a continuous and smooth curve. To solve the second problem, this paper proposes algorithms for automatic dynamic smoothing of the primary path using a tracking differentiator with sigmoid corrective actions. Algorithms for setting the gains of the differentiator are developed, considering a set of design constraints on velocity, acceleration, and jerk for various mobile robots. When tracking a non-smooth primary path, the output variables of the differentiator generate smooth trajectories implemented by a mechanical plant. It is shown that the tracking differentiator with a different number of blocks also generates derivatives of the smoothed trajectory of any required order, taking into account the given constraints. Unlike standard analytical methods of polynomial smoothing, the proposed algorithm has a low computational load. It is easily implemented in real time on the on-board computer. In addition, simple methods for modeling a safety corridor are proposed, taking into account the dimensions of the vehicle when planning a polygon with stationary obstacles. Conﬁrming results of numerical simulation of the developed algorithms are presented.


Introduction
Nowadays, mobile robots are used in a wide variety of processes, and their applications are constantly expanding.They are remotely controlled or have different degrees of autonomy.In this paper, nonholonomic wheeled robots are the plant of research.Currently, two-wheeled unmanned vehicles are very popular due to their efficient design.They have two actuators with three degrees of mobility, so they cannot realize three independent movements in longitudinal and transverse directions, as well as rotation around the vertical axis [1,2].
With remote control, the route of the mobile robot is not known in advance.The operator generates real-time control signals that determine the reference direction and the velocity of the robot.In range conditions or in an environment with fixed obstacles, the autonomous robot's movement can be designed in advance, and the reference trajectory can be stored on the on-board computer.In both cases, the basic goal of automatic robot control is to track set-point inputs designed either online or offline.
In path stabilization problems, the robot's velocity is not controlled and is constant [3].In trajectory control problems, it is necessary to determine not only the path but also the reference velocity and acceleration of motion at different sections of the curvilinear path [4].For a complete design of an automatic motion control system for a mobile robot, information about the coordinates of the reference trajectory and their higher-order derivatives is generally required.A closed-loop tracking and software control system will give good performance when the reference trajectory is smooth and achievable.If the reference trajectory has special points where the design constraints on vehicle velocity, acceleration, or jerk are violated, then spikes in control torques should be expected at these points.If such points occur regularly, this will lead to premature wear of the actuator.When tracking a non-smooth trajectory, the robot may stray from the route and have an accident.Thus, in trajectory control problems, the most important problem is to design smooth trajectories with continuous curvature that satisfy the design constraints of a particular wheeled robot and are achievable [5].
Planning an achievable trajectory for a mobile robot usually consists of two steps: (i) finding a path in the form of a sequence of discrete waypoints and (ii) converting this sequence into a continuous and smooth curve considering the movement time.Currently, many algorithms have been developed with different approaches that effectively solve the first problem depending on the reference conditions [6].These include graph search algorithms (A* [7], etc.), sampling-based algorithms (RRT [8], etc.), reaction-based algorithms (DWA [9], etc.), and their various modifications and combinations [10].
The subject of this research is the second problem.Interpolating curve algorithms are currently used to solve it [11].They allow one to obtain an analytical description of the trajectory and, consequently, its derivatives of any required order.To obtain smooth curves, they use fragments of circles, polynomial splines [12], clothoids [13], generalized Cornu spirals [14], etc.A typical modern solution is B-splines [15].The problem is reduced to the computation of polynomial coefficients of a stepped function of time (polynomial) that connects several neighboring points considering a specified time.For example, a smooth trajectory fragment can be obtained using a cubic B-spline and five reference points.To generate a trajectory subject to the design constraints of the mobile robot, additional algorithmization or increasing the degree of the polynomial is required.For example, a polynomial curve of degree seven is used to constrain the first derivative in [16].In addition, a smooth join of the separate fragments of the trajectory needs to be ensured.The implementation of analytical interpolation methods at the trajectory planning stage does not cause difficulties.However, they are not suitable for smoothing complex routes with a large number of reference points and fragments of different shapes in real time due to the large computational load and computation time [17].In addition, the complex analytical trajectory can take up much memory when stored on the on-board computer.
To solve this problem, we developed a universal dynamic smoothing method that automatically generates a smooth curve that satisfies the design constraints of a particular robot for any number of reference points.At the same time, complicating the shape of the reference curve does not lead to a more complex algorithm and increased computational operations.Only nominal values of vehicle velocity and acceleration, etc., are required for its presetting.The algorithm can be implemented both offline at the motion planning stage and online on an on-board computer with any primitive software.The main difference from the above-mentioned interpolating curve algorithms is that within the framework of the dynamic smoothing method, we obtain a smooth trajectory and its derivatives of any required order in the form of continuous vector signals without analytical description.It is possible to store sequences of waypoints in the on-board computer memory (it takes very little space) for different operational scenarios and activate their processing in real time.
The complexity of the above advantages is provided via a dynamic tracking differentiator.In traditional applications, this computational algorithm is used to process the sensor signal of a state variable of the control plant and recover its first derivative.This allows no installation of velocity or acceleration sensors (they usually produce noisy signals [18]).
The tracking differentiator is an autonomous system of second-order differential equations having a canonical input-output form [19][20][21][22][23]. The input corrective influence is formed as a nonlinear function of the differentiator variables and the external sensor signal to stabilize the tracking error (the uncertainty between the differentiator output variable and the sensor signal) and to filter out noise disturbances if they are present in the measurements.As a result, the state variables of the differentiator reproduce the useful measurement signal and its derivative, which are used in the feedback loop.In the deterministic case, discontinuous correction [21,22] is used, and in the general case of noisy measurements, continuous analogs of the sign function (e.g., hyperbolic tangent [20]), more complex fuzzy logic methods [23], and others are used.
In this paper, we propose to apply a tracking differentiator with continuous correction to track the initial non-smooth trajectory of a wheeled robot.First, we need to smooth the external signal subject to design constraints.Second, we recover its derivatives of any desired order.Increasing the dynamic order of the differentiator solves the second problem.To solve the first problem, we developed a special algorithm for the design of corrective influences with the generation of S-shaped sigmoidal local links.Sigma functions (modifications of the hyperbolic tangent) are smooth and bounded [24].The gains that depend on the design constraints of the robot define the domain of their values.In a closed-loop system, the output variable of the differentiator tracks the external signal with some accuracy, which is a measure of the smoothing of the external signal at angular points.The other variables of the differentiator (derivatives of the output variable) "track" the introduced sigmoid functions and, therefore, do not move out of the specified ranges.Thus, the tracking differentiator automatically generates a smooth trajectory and its derivatives in real time with specified constraints.These signals are used in the tracking system of the control plant as reference achievable influences.
In our previous paper [25], this method was tested to generate achievable 4D trajectories for an aircraft-type unmanned aerial vehicle (UAV) center of mass and showed excellent performance.In this paper, we extend this method to a set of motion-planning problems for a wheeled robot.The contributions of this work and its scientific novelty are described below: 1.
A proof of the criterion of boundedness of the solution of an elementary dynamical system with sigmoidal feedback under the influence of an external uncontrolled disturbance is provided.The radii of the convergence domain of the solution are formalized (Section 2.2), which are further used to adjust the gains of the tracking differentiator.

2.
It is shown that by designing a single-block tracking differentiator with the minimum possible dynamic order, it is possible to smooth the external signal considering the design constraints on the speed and acceleration of the mobile robot.Design methods and a setting algorithm for an arbitrary-order tracking differentiator with sigmoidal local links are developed.(Section 2.3.3).

3.
We propose computationally simple algorithms for modeling a safety corridor considering the dimensions of the wheeled platform (Section 2.3.4).
The remainder of this paper is organized as follows.The following section presents the theoretical methods and their rationale.In Section 2.1, the basic dimensions and structure of the tracking differentiator for the generation of achievable 3D trajectories considering constraints on the velocity, acceleration, and jerk of the wheeled robot based on the transformation of the mathematical model of the motion of the center of mass of the wheeled platform into the canonical Brunovsky form are determined.Section 2.2 describes the properties of the sigma function and sigmoid feedback.Section 2.3 presents principles for designing and tuning tracking differentiators of various dimensions.Reasons for designing tracking differentiators with different block numbers are listed.Some aspects of designing paths and polygons are considered.In Section 3, numerical simulation results of the developed algorithms are presented and analyzed.In particular, the results of smoothing a reference sequence of waypoints with different sets of design constraints are demonstrated.

Equivalent Transformation of the Equations of Motion of a Wheeled Robot
The plant for trajectory planning is an unmanned wheeled platform with two driving wheels.Consider a dynamical system that describes the motion of the center of mass of this vehicle in a stationary Cartesian coordinate system Oxy [4,26]: where the first three equations describe kinematics, and x, y are the coordinates of the platform center of mass (output variables); θ is the angle between the axis Ox and the center line of the platform, which coincides with the direction of the velocity vector; θ determines the orientation of the platform concerning the stationary coordinate system.The variables in the last two equations of the system (1) are linear velocity v and path curvature κ; κv is the angular velocity of rotation of the vehicle from its center of mass in motion mode, v = 0; u v , u κ are the control torques generated by the actuators.The dynamics of the actuators are not considered in this model.
To generate motion trajectories with specified velocity and acceleration constraints, it is more appropriate to represent the model of the control plant in canonical form concerning the output variables and their derivatives up to and including second order.In [26], it is shown that the system (1) is differentially flat concerning flat output y 11 = x, y 12 = y.By extending the state space (namely by introducing an auxiliary dynamic variable ξ(t)), System (1) can be represented in the canonical Brunowski form with two inputs and two outputs: ... where .
The canonical system (2) is equivalent to the extended system (1) of the form .
A specific mobile robot, depending on its design and the actuators used, has physical constraints on the acceptable velocity, acceleration, and jerk, i.e., on the first, second, and third derivatives of the output variables of System (1).Let us introduce the corresponding constraints for the canonical variables of System (2) and new controls Here and below * is the l ∞ -norm of the vector, namely y T , y (i) y 1j (v 4 ), j = 1, 2. Let us note that, if necessary, constraints on higher-order derivatives can be introduced in a similar way.For example, y 1 (t) ≤ Y 4 limits the rate of control change.In general, the order of the derivatives to be constrained is equal to the relative degree of the full model of the control plant.

Properties of the Sigma Function and Sigmoid Feedback
In further constructions, we will use as corrective actions in the tracking differentiator a modification of the hyperbolic tangent th(x) = 1 − 2/(1 + exp(2x)), namely the sigma function σ(x) = −th(−x/2), with two gains p, l = const > 0 in the form The sigma function ( 5) belongs to the class of logistic functions with the following features: they are defined on the entire numerical axis, smooth, odd, bounded, and have an S-shape.Unlike other logistic functions, such as rational and root sigmoid and arctangent, etc., the sigma function is easy to compute since the corresponding Maclaurin series converges on the entire numerical axis, and its derivative has a recursive form: In a small neighborhood of zero, the sigma function is equivalent to the linear function pσ(lx) ∼ x→0 0.5plx, and as the gain l increases, it is closer to the sign function pσ(lx) ∼ l→+∞ psign(x).The above features allow us to constrain the sigma function from below to a piecewise linear function where the points σ(±l∆) separate the sigma function into a conditionally constant and a conditionally linear function.Let us denote Let us consider the features of sigmoidal feedback in the problem of stabilizing a state variable under the influence of an external disturbance.We will use the example of a first-order system .
where x ∈ R is the state variable, u ∈ R is the sigmoid control, and does not decay with time, then in the closed-loop system (9), it is possible to ensure stabilization of the state variable with some accuracy if certain conditions are fulfilled.

Lemma 1.
In System (9), let the choice interval gain p > 0 be bounded from above: p ≤ p.Then, for any initial condition x(0), he solution of the system x(t), t ≥ 0 will be bounded if and only if F < p.

Proof of Lemma 1.
To analyze the steady state of System ( 9), let us use the second Lyapunov method and estimate the derivative of the Lyapunov function V = x 2 /2.Outside the neighborhood of |x| ≤ ∆ = c/l, the sigma function from below is bounded by the constant (7), so considering (8), we have .
V > 0, and for all p < F/σ(c), System (9) is not stable, and its solution is unbounded lim t→+∞ x(t) = ∞.If F = pσ(c), then .V = 0, and System ( 9) is neutral, is on the stability boundary, and its solution can be considered constant x(t) ≈ x(0).If p ≤ p ≤ F ⇔ F/p ≥ 1, then by virtue of 0 < σ(c) < 1 ( 5), (8) in both cases, we have a strict inequality p < F/σ(c), i.e., the stability of System ( 9) cannot be ensured.Thus, it has been proved that it is necessary to fulfill the condition F < p.
V < 0 and for all p > F/σ(c), by virtue of ( 7), the state variable converges to the following neighborhood: i.e., the solution of System ( 9) is bounded.If p > F ⇔ F/p < 1, then there exists an acceptable value of c * : F/p < σ(c * ) < 1, at which there is a non-empty interval for choosing a gain p * : F/σ(c * ) < p * ≤ p, that ensures (10).Thus, the sufficiency of the fulfillment of condition F < p is proved.
For fixed values of c * , p * chosen from the above intervals, consider different variants of the initial conditions of System (9).If |x(0)| ≤ ∆, then inequality ( 10) is correct at t ≥ 0. If |x(0)| > ∆, the state variable reaches the specified neighborhood in a finite time t * > 0, and inequality (10) Thus, if F < p, then for any initial conditions, the solution of System ( 8) is bounded.
Note that in the neighborhood of |x| ≤ ∆ (10), the sigma function is bounded from below by the oblique line (7), and for the derivative of the Lyapunov function, considering (8), let us take the following estimate: .
to which the state variable converges.The radius of the domain (11) is smaller than the primary estimate (10).It follows from ( 11) that the radius of the convergence domain ( 11) is directly proportional to the parameter c > 0 and inversely proportional to the gains p, l > 0. The control u in System ( 9) is bounded by the amplitude p : |u(t)| < p, t ≥ 0, and the control rate .u also depends on the second gain l.The maximum values of p, l in practical problems are usually constrained for physical reasons.As a criterion for choosing the parameter c > 0, we take the minimum of the basic estimate of the modulus of the control rate [24].Considering (6), the following estimates .
are correct.If ∆, F are fixed, the reference condition becomes c/σ 2 (c) → min, c > 0 and is fulfilled at c ≈ 2.2.In the following research, this value, as well as l∆ = 2.2, σ(2.2) ≈ 0.8; 1/σ(2.2) ≈ 1.25, will be used.At that, in the indicated intervals of ( 7 In the next subsection, the discussed sigmoidal feedback will be used to design the synthesis of the tracking differentiator.Under the made assumptions, we can introduce the notation t 1 := 0, t 2 := T, with the velocity vector v, which coincides with the center line of the platform, directed tangentially to the reference trajectory.Under the assumptions made, we can introduce the notations t 1 := 0, t 2 := T, with y 1 (0) ≈ χ 1 (0),

Design of
= 0, and the velocity vector → v (0), which coincides with the center line of the platform, directed tangentially to the reference trajectory.
Let us highlight the main problems that need to be comprehensively considered when defining a trajectory for a single-wheeled robot: 1.
The path must be implementable by a mechanical plant, i.e., it must be sufficiently smooth and have continuous curvature, minimum requirements In addition, the reference motion velocity and the acceleration and jerk must not exceed the robot's constraint values (4), i.e., 2.
The path must be safe and not lead to collisions with fixed and dynamic obstacles.It is therefore necessary to consider the dimensions of the vehicle, the configuration of the polygon, etc.

3.
A reference trajectory must satisfy various criteria, which, depending on the robot's mission and work scenario, are formulated as various terminal and optimization problems.For example, to reach the end point of the route in minimum time, to create the shortest route by avoiding reference points, to perform a task with minimum energy consumption, to solve pursuit or evasion problems, etc.
In the following subsubsections of this section, the first problem is discussed in various formulations, and the second problem in terms of vehicle dimensions in planning a polygon with stationary obstacles.The third problem is beyond the framework of this research and is not considered here.

Designing a Three-Block Tracking Differentiator
Let the reference influences χ 1 (t) = (χ 11 (t), χ 12 (t)) T enter the information and control the on-board system of the robot in real time in the form of a deterministic vector signal.Their analytical description is unknown, and hence, the current values of reference signal derivatives are also unknown.It is assumed that the reference trajectory satisfies the safety requirements, but its realizable conditions (15) and ( 16) are partially fulfilled.Namely, the time functions χ 11 (t), χ 12 (t) are either continuous but non-smooth or piecewise continuous.On the continuity intervals, the derivatives of these functions do not exceed the allowed values (16).The number of special points (junctions or points of the discontinuity of the first kind) is bounded.At these points, there exist bounded left and right first derivatives.Within the framework of the method used, an additional condition of type ( 14) is introduced: .
In real time, the problems of smoothing the vector signal t and restoring its first and second derivatives are posed considering the constraints of the robot under consideration (4).To solve these problems, we propose using a dynamic model (tracking differentiator), which has a canonical form similar to (2): .
where y 1 (t), respectively.The analogs of the controls u 1 , u 2 of the canonical system (2) are the corrective sigmoidal influences w ∈ R 2 , by means of which it is necessary to ensure that the output variables of the differentiator z 1 (t) track the external non-smooth signals χ 1 (t) with some accuracy.Tracking accuracy depends on the constraints of a particular robot (4), which need to be ensured in the closed-loop system (18) by choosing the gains of corrective influences: Under these conditions, the differentiator variables (18) generate the achievable trajectory in signal form and its first and second derivatives.These variables are new reference influences and are used to generate control influences u 1 , u 2 in the tracking system of the mobile robot.
To design a tracking differentiator (18) with the indicated features by means of nongenerated replacement of variables, we proceed to a new coordinate frame, which is formed by the tracking errors and sigmoidal local connections: where e i = (e i1 , e i2 ) T , i = 1, 2, 3; Let us also introduce a sigmoidal corrective influence where p 3 , l 3 = const > 0, σ(l 3 e 3 ) = (σ(l 3 e 31 ), σ(l 3 e 32 )) T , and note the closed-loop system  .
The tracking problem is transformed into the problem of stabilization of the virtual system (22), considering the constraints (19).Let us first formulate sufficient conditions of the boundedness of solutions of the system (22) without considering (19).To simplify the setting, let us establish the following initial conditions in the tracking differentiator (18): z 1 (0) = χ 1 (0), z i (0) = → 0 , i = 2, 3, which by virtue of σ(0) = 0, will ensure zero initial conditions for the state variables in the virtual system ( 22): As we see, the structure of Equation ( 22) is similar to the structure of the closed-loop system (9).We can treat all terms on the right side of each equation, except for the sigmoidal ones, as external bounded disturbances.While analyzing the system (22), let us take as a base the proof of lemma 1.By virtue of (24), the variables of the system (22) are in domains similar to (10) from the beginning at t = 0: Sufficient conditions under which inequalities (25) will be correct at t ∈ [0; T], similarly to (14), have the form Then, for the derivatives of .e 1 , .e 2 , and, consequently, for Equation (23), similarly to (12), we obtain the estimations .
By virtue of a priori assumptions, namely the lower inequalities (4), considering ( 14), we assume the system (29) to be joint.It has infinitely many solutions.
Under the assumption of non-smoothness of the external signal χ 1 (t), the accuracy of stabilization of the variables e 2 (t), e 3 (t) (25) is not important.The tracking error at special points is bounded e 1 (t) ≤ 2.2/l 1 .In the steady-state mode, we have an estimate similar to (11): To minimize the steady-state tracking error (30), an iterative procedure for setting the gains of sigmoid functions based on the system of inequalities ( 29) is required.For example, the following variant is proposed: where β 1,2 are very small positive constants that are introduced so that the left inequalities in the system (29) remain strict.This procedure allows us to accept the maximum possible value of p 1 .However, in general, this does not provide an absolute minimum of the tracking error estimate (30).If necessary, this figure can be improved by decreasing p 1 and increasing l 1 .
As a result, the tracking differentiator (18) with the corrective influence (21) will be realized as a closed-loop system .
The dynamic smoothing algorithm (32) along with the gains setting (31) has two universal features.First, it can be used for all wheeled robots whose dynamic model is flat and represented in the form of (2).The peculiarities and constraints of a particular vehicle will be considered when the corresponding constants are substituted into the Equation (31).The values of the corresponding gains p i , l i are calculated at the preparatory stage.They do not depend on the features of the external signal χ 1 (t) and remain unchanged.
At different nominal values of (4), the signals z i (t) will automatically approximate the external signals and their derivatives in the neighborhood of special points in different ways.In the steady-state mode, in the admissible cases, their forms will be almost identical.The results of the corresponding numerical simulations are presented in Section 3.
Second, by changing the dimensionality of the tracking differentiator blocks, it is possible to automatically generate smoothed trajectories with defined constraints in phase spaces of different dimensionalities.For example, the system (32) at ∀z i ∈ R generates an achievable trajectory for a single-channel control plant (e.g., a single-link manipulator), and at ∀z i ∈ R 3 it generates a spatial trajectory for a UAV [25].
From a computational viewpoint, the dynamic smoothing algorithm (32) is a sigmoid computation with three integration operations.Their implementation is not difficult using any software.The counting time of the algorithm is negligible and does not lead to lag in real-time operation when vector signals z χ 1 (t), and w(t) ≈ ... χ 1 (t) are input to the tracking system of the mobile robot as needed.

Reasons for Designing Tracking Differentiators with Different Block Numbers
As was shown in the previous subsection, it is sufficient to apply a dynamic model consisting of three blocks (32) to smooth external signals and to recover their derivatives of the first, second, and even third orders.However, the number of blocks of the tracking differentiator can be both reduced and increased depending on the conditions of the problem to be solved.The main factors influencing the dimensionality of the differentiator are as follows: 1.
Methods used to design control in a tracking system of a mobile robot; 2.
The considered dynamic order of the control plant; 3.
For example, to solve the problem of tracking by the output variables x(t), y(t) of System (1) of reference signals under the influence of parametric and external disturbances, we can apply a decomposition method for the design of sigmoid feedback, similar to the method we applied in the previous section for the design of a tracking differentiator [24].In this method, the first derivatives of the reference influences are treated as external disturbances (as in System ( 22)), and only the reference influences are used to form the static feedback.Another approach is to use mixed-variable observers in tracking systems [27].In this case, for feedback design, it is not necessary to obtain separate estimates of the derivatives of the reference influences.Thus, in both approaches, we only need to know the reference signals.If these signals do not contain parasitic disturbances, a singleblock tracking differentiator with the minimum possible dynamic order can be used to smooth them: . where The inequalities for the choice of gains p 1 , l 1 = const > 0 considering the constraints on the velocity and acceleration of the mobile robot, similar to ( 29), (31), are of the form The shapes of the smoothed curves obtained using the tracking differentiators (32) and (34) under the same constraints (4) will be almost identical.
The signals of the derivatives of the reference influences are usually needed in program control systems [28] as well as when using the feedback linearization method [3].The order of the required derivatives depends on the number of dynamic links to be considered and is equal to the relative degree of the model used in the control plant.
Let derivatives up to and including n-th order be required for the design of the tracking system of a wheeled robot.The tracking differentiator, which smooths the external deterministic signal χ 1 (t) and recovers its n derivatives, consists of (n + 1) blocks and has a common order 2(n + 1): .
where z i ∈ R 2 , i = 1, n + 1 are differentiator variables, the vector variable z 1 (t) defines the position of the center point of the wheeled platform, and the variable z i (t) is its To set the tracking differentiator (35), the constraints of a particular mobile robot on the higher derivatives are used.To fulfill the constraints on the state variables of the differentiator, its corrective influences, and their rate of change (i.e., constraints on the robot's control influences and their rate of change), nominal values (n + 2) are required Equations for selection gains of tracking differentiator (35) at which the constraint of its variables is ensured z are similar to the system of dual inequalities (29) and have the following form: The iterative procedure for setting the gains of sigmoid functions based on the system of inequalities (36) minimizing the steady-state tracking error (30) is similar to (31): where β i , i = 1, n are very small positive constants.Finally, the third factor influencing the choice of the dynamic order of the tracking differentiator is related to the need to filter the external signal if it additively contains parasitic noise χ 1 (t) = χ 1 + η(t), for example, white noise η(t) with zero mathematical expectation and bounded variance.
Note that the noisy signal enters the input of the tracking differentiator (35) and passes through a chain of integrators.Thus, there is natural filtering of the input signal, which improves with increasing number of integrators.Analytical analysis and simulation results have shown that extending the dynamic order of the tracking differentiator is an alternative to installing low-pass prefilters [29].Two or three integration operations are often sufficient to obtain good performance.For example, if the recovery of the derivatives of the external signal is not a problem, a two or three (32) block tracking differentiator should be used instead of a single-block differentiator (33) for smoothing and filtering the external signal.In turn, instead of a three-block differentiator (32), a four-or five-block tracking differentiator must be used in order to obtain signals of satisfactory quality up to and including the second derivative.
Note that, as in the Kalman filter setting, when selecting gains of sigmoidal functions, we have to compromise between tracking accuracy and filtering features of the tracking differentiator [23].Therefore, another algorithm is needed instead of the procedure for setting gains p i , l i (37).However, a detailed analysis of this problem is beyond the framework of this research.

Some Aspects of Designing Paths and Polygons
The use of tracking differentiators can greatly simplify the computational aspect of planning achievable trajectories offline on a polygon with stationary obstacles.
Let a wheeled robot have a primary operating scenario in the form of a sequence of 3D points defined in a fixed Cartesian coordinate system Oxy, considering time Connecting all neighboring pairs of point (x i , y i , t i ) and (x i+1 , y i+1 , t i+1 ) by segments we obtain a primitive non-smooth 3D trajectory and, accordingly, the primary reference influences for both output variables x(t), y(t) of the system (1): In standard approaches, special analytical approximations are used to smooth the angles of the composite trajectory (39) [11].In this case, the problem of ensuring the specified constraints on the velocity and acceleration of the robot requires additional algorithmization.If the singular points are sufficiently numerous, analytical smoothing methods provide a large computational load, which makes it difficult to use them in real time.In addition, composite analytical reference influences take up much memory when stored on the on-board computer.
The application of dynamic smoothing algorithms for primitive trajectories (39) completely eliminates these problems.Only a sequence of points (38) is stored in the on-board computer memory, and smoothing of the corresponding primitive trajectories (39) is performed in real time using a tracking differentiator (35), (37) of the required dimensionality.
Note that this method is not intended for solving terminal problems and does not provide an exact fit to the reference points.However, it can be used for offline simulation on a compressed time scale.We can quickly correct the smoothed trajectory obtained after the first run of the algorithm and achieve the desired result by correcting the coordinates of particular points (38).At the same time, visual results will be obtained quickly and without cumbersome analytical calculations.For such a process, it is sufficient to use a single-block tracking differentiator (33).
Let us consider some aspects of planning a safe route for an oversized wheeled vehicle on a polygon with stationary obstacles.The opposite task is to arrange objects on the polygon (factory workshop, warehouse) so that the robot can safely fulfill the operating scenario.Therefore, it is necessary to consider the dimensions of the vehicle when planning the polygon.The developed dynamic smoothing algorithm can be applied to graphically represent the position change not only for the center of mass of the wheeled platform but also for its corner points.In this problem, it is better (but not necessarily) to use a two-block tracking differentiator .
Its parameters are determined from the procedure (37) at n = 1 and have the form: Let us present the corresponding computations for a symmetrical rectangular wheeled platform whose center of mass is located at the intersection of diagonals in the middle of the center line.Let us introduce the following notations: ρ > 0 is the distance (with a small margin) from the center of mass to each corner point of the platform; α < π/2 is the value of the angle between the center line of the platform and its diagonals; (z 11 (t); z 12 (t)) is the current coordinates of the center of mass of the platform; and (z 11,j (t); z 12,j (t)) is the current coordinates of the corner points of the platform, j = 1, 4. Recall that in System (1), θ is the angle between the Ox axis and the center line of the platform, which coincides with the direction of the velocity vector, θ = arctg( .y 12 / .y 11 ) (3).First, the waypoints of the mobile robot (38) are planned on the polygon.The primary reference influences (39) are entered into the input of the tracking differentiator (41) with the output variables z 11 (t), z 12 (t), which simulate the position of the platform's center of mass in the stationary Cartesian coordinate system.The formulas for computing the current coordinates of the platform corner points are as follows and are correct at 0 ≤ θ ≤ π/2.Technically, at other directions of the velocity vector, we should change the signs ± in Equation (42) and use additional logic.Let us not complicate the algorithm and track the direction of motion of the platform.When changing the quadrant in which the angle θ is located, the angular points on the right side z 12 (t), z 13 (t) will represent the path of the points on the left side z 11 (t), z 14 (t) and conversely.Thus, curves (z 11 (t), z 12 (t)) and (42) shown on the same plot will indicate the dimensional footprint of the wheeled platform, which is required.
Such representations are a visual and convenient tool for polygon design, object placement, and trajectory planning.
The developed algorithms complex is a convenient tool for mobile robot motion planning and polygon design.Figure 1 shows a flowchart of the robot motion planning process on a polygon with stationary obstacles.21 22 and are correct at .Technically, at other directions of the velocity vector, we should change the signs ± in Equation (42) and use additional logic.Let us not complicate the algorithm and track the direction of motion of the platform.When changing the quadrant in which the angle θ is located, the angular points on the right side ) ( ), ( 13 and (42) shown on the same plot will indicate the dimensional footprint of the wheeled platform, which is required.
Such representations are a visual and convenient tool for polygon design, object placement, and trajectory planning.
The developed algorithms complex is a convenient tool for mobile robot motion planning and polygon design.Figure 1 shows a flowchart of the robot motion planning process on a polygon with stationary obstacles.This process includes the following: 1.
Designing a base set of 3D points (38) for a specific workspace considering obstacles, robot velocity and turning radius constraints, route length, etc.; 2.
Smoothing of the primitive trajectory using a tracking differentiator.This process is a numerical solution of differential Equation (32), i.e., calculation of the sigmoid and integration operations; 4.
Appendix A presents a structural diagram created by simulating the algorithms developed in the MATLAB Simulink software environment.The main interface of Simulink is a graphical block diagraming tool and a customizable set of block libraries.Therefore, the structure of the graphical block diagram in the MATLAB Simulink software environment (Figure A1) is similar to that of the flowchart shown in Figure 1.At the same time, Figure A1 clearly demonstrates the computational process and the computational operations performed.
Note that in the first stage, machine learning methods can be effectively used to generate primitive trajectories considering various criteria [30].The result is the solution of discrete optimization problems in the form of a set of discrete waypoints.As mentioned, designing a smooth and achievable trajectory on their basis requires additional tools [16].
In the next subsection, the results of numerical simulations of this and other developed algorithms will be presented.

Results
To verify the developed algorithms, simulation modeling in the MATLAB Simulink environment was carried out.We used the Euler method with a constant step of 10 −2 for numerical integration.The results of two experiments using a three-block tracking differentiator (32) are presented below.The first experiment demonstrates the smoothing processes of primitive trajectories for the robot's center of mass under different sets of constraints (4).In the second experiment, for one of the smoothed trajectories, a motion trajectory is simulated considering the vehicle dimensions.
In the first experiment, the design constraints of two civilian wheeled robots moving in an urban environment were used to set the gains of the tracking differentiator.Based on the inequalities (29), the following gains of the tracking differentiator (32) were adopted for the first robot with a velocity of up to 2.3 m/s: For the second robot with a developing velocity of up to 1.9 m/s, the following gains of the tracking differentiator (32) were adopted: The robots' motion was planned along a rhombus defined in a fixed Cartesian coordinate system.Thirteen waypoints (x i , y i ) were randomly marked on the sides of the rhombus.Based on them, we generated two sets of 3D waypoint sequences (x i , y i , t i ), i = 1, 13 (38), which differ only in the values of time t i , by which the average velocity of motion on each path segment is defined and limited.In the first sequence, the maximum velocity was limited by 2.5 m/s and in the second by 2 m/s.We purposely used values slightly higher than the above-mentioned robot velocities.
In future research, it is planned to develop the algorithm for considering the dimensions of the control plants under consideration, considering the obstacles present in the course of movement along the developed route.However, at the moment, the dimensions of the control plants are already considered on the basis of the obtained data using the tracking differentiator (32) with the computationally simple algorithm (42).This is also an effective tool for trajectory planning, considering the objects and obstacles present on the polygon.
vector variables of the differentiator, which are analogs of the robot position y 1 (t) = (y 11 = x, y 12 = y) T , its velocity .y 1 (t), and its acceleration ..

Figure 1 .
Figure 1.Flowchart of the trajectory planning process.Figure 1. Flowchart of the trajectory planning process.

Figure 1 .
Figure 1.Flowchart of the trajectory planning process.Figure 1. Flowchart of the trajectory planning process.
Tracking Differentiators 2.3.1.Problem Definition Let χ 1 (t) = (χ 11 (t), χ 12 (t)) T be a vector of reference influences for the output variables of System (2) y 1 (t) = (y 11 = x, y 12 = y) T , which determines the reference trajectory of the wheeled robot at a time interval t ∈ [t 1 ; t 2 ]; namely, it specifies in the stationary Cartesian coordinate system Oxy the track path (χ 11 , χ 12 ) ∈ R 2 , as well as the motion velocity .χ 1 (t) = ( .χ11 (t), .χ12 (t))T , acceleration .. χ 1 (t) = ( .. χ 11 (t), .. χ 12 (t)) T , etc.The problems of starting, stopping the vehicle, and reaching a reference point of a reference trajectory with a reference orientation are not considered in this paper.It is assumed that the interval [t 1 ; t 2 ] corresponds to a continuous motion of the robot without stops, and the time of staying on the route t 2 − t 1 is acceptable.