CC BY
28
ISSN 1992-6502 (Print)_TB&C/ftVH/U/IC ^¡jf/^Q^^lj _ISSN 2225-2789 (Online)
Vol. 18, no. 5 (66), pp. 31-41, 2014 http://journal.ugatu.ac.ru
UDC 004.65
Classes of trajectories of persecution
1 2 A. V. Zakharov , M. R. Bogdanov
1 ralzax@yandex.ru, 2 bogdanov_marat@mail.ru
1,2 Bashkir State Pedagogical University named after M. Akmullah, Russia 1 Bashkir State University, Russia
Submitted 2014, June 10
Annotation. Classification of laws of persecution motion is offered in the paper. Researching of persecution modes in number of specific cases was performed. Various modes of motion control using GPS/GLONASS are considered.
Key words: management of dynamic motion; persecution trajectories; GPS; GLONASS.
ABSTRACT
Precision of positioning in GPS/GLONASS systems is getting better extending areas of application of global positioning systems data. Improving of positioning accuracy based on enhanced mathematical algorithms shown in our work [1]. It is getting possible to control of motion of carrier. L.S. Pontryagin and coworkers developed effective in practical applications theory of control of dynamic motion [2]. Dynamic motions defined by the system of equation are considering.
= f1 (x 1, ...,xn,u1,. . .,ur,t),i = 1 . . .n, (1)
where x1, ..., xn are the phase coordinates and velocities of the managed object, u1, ..., ur are functions managing turn wheel, brakes and accelerators forming the closed limited set U in the r-dimen-sional space of object control, t is the time [3, 4].
Dynamic movement control problem is to selection the motion control signals ul = u'(t) for which the system (1) gives an unique solution X = xl(t), corresponding to the initial data xl(t0) = (x0)i, i = 1, ... , n, which satisfy to our predefined condition of motion. The integral feature usually determines optimality condition of motion
J = it2fn+1 (x 1.....Xn,u 1.....ur,t) dt , (2)
which has the property of transitivity for autonomous systems.
In the case where the optimality criterion is a requirement of minimal motion time between two given points, or closed region of the phase space of generalized coordinates and velocities of the managed moving body, the problem of managing be-
comes a problem of performance according to L. S. Pontryagin.
Performance problems are simplest because op-timality function in this case becomes equal to one. These ones are well known and there are number of solutions for specific classes of movement managing problems [2]. L. S. Pontryagin and coworkers found partially continuous managing signals u(t) located on the border of the limited set of values U during all the time are satisfy to solution of performance problems. Therefore, to fastest movement from one point of trajectory to another driver should speed up at one point and have to maximum push on the brake at another point. In many practical situations, that a strategy is a good choice, but every driver knows the driving rule: accelerates smoothly and brake smoothly. It satisfies to another optimality condition of motion. Each optimality condition defines the own class of optimal controls and optimal trajectories. In the area of intersection of classes of managed movement trajectories simultaneously satisfies two or more of the conditions of optimality.
The simplest class of managed motion is a straight linear motion with constant velocity. Managing signals are constant and equal to zero [5].
We choice an optimal motion with minimal integral curvature as a generalized straight linear motion
; = £ 1+tg ( t) u ( t) d t (3)
or more common class - the motion with minimal torsion. Trajectories of performance with partially constant acceleration are given by quadratic splines. They cannot always approximate trajectory of per-
secution on a finite time interval, which have a smooth managing.
The problem of motion managing for a minimum curvature and torsion we call problem of persecution.
Communication between GLONASS/GPS receiver with satellite system is carried out at an interval At. However, persecution time may be both larger and smaller At.
Determining of the persecution time is a very important problem both in theory of managing and in theory of persecution. If the time between communication sessions is less than the persecution time At < x, then each session gives final trajectory persecution correction, time of persecution and managing signals [6].
Anyway, final time of persecution satisfies to own optimal law of persecution x(i), usually different from the straight linear motion.
We proffer classification of the available laws of persecution based on geometric and algebraic approaches.
We can solve the generalized two point Cauchy problem for each class of trajectories and find control commands in each class of trajectories for the dynamic laws of persecution based on the selected class of the persecution.
Finally, the computer subsystem initially selects a class of trajectory then finds laws of dynamic control and solves the persecution problem based on the relative positions and velocities of interceptor and target.
In the simplest case, when the target does not deviate from the predetermined path, the command of persecution control is equal to zero. If the target is retarding linearly, the persecution command for increasing of linear tangent velocity is run. There is a finite set of situations that require more complex trajectories of persecution, which are analyzed below.
THE GEOMETRICAL APPROACH
1. CRITERIA OF OPTIMALITY AND EQUATION OF PERSECUTION
Its known geometrical characteristics of trajectory of movement includes unit vectors for accom-
dv
panying orthoreper: — = x - unit vector of velocity
d2r
tangent to the trajectory — = o • n, n - unit vector of normal to the trajectory, c - curvature of the tra-
[civ
—l = o^ , b - unit vector of binormal of the trajectory, s - length of the trajectory's arc, - a natural parameter; o2 = — o _ 2 •
[dr d2r~l (d3r\
— x —I • (— 1 - second curvature or twisting of
the trajectory [3]. A straight line has zero curvature o = 0 and as a result zero twisting o2 =0. They forms the first geometrical class of trajectories widely using in dynamic games. The next class of trajectories are the zero torsion trajectories with smallest curvature.
1.1. Plane motion with minimal integral curvature. Let us require functional (3) on the optimal trajectories of persecution should be smallest without using of dynamic equations of managed object (1). Using the time t as a parameter we have following:
x(t)y(t)z(t) - x(t)z(t)y(t) - x(t)y(t)z(t) + +x(t)z(t)y(t) + x(t)y'(t)z(t) — x(t)z(t)y(t)=0,
(1.1.1)
which means the linear dependence of the velocity vectors, acceleration and acceleration of acceleration:
"x(0" xit) xit)'
• m = a(t) ■ m №
m m m
(1.1.2)
Equation (1.1.2) represents a general type of linear homogeneous differential equation of the second order relatively components of interceptor's vector of velocity. This one has two fundamental solutions relatively vector of velocity leading to the plane motion. Indeed, the movement will be flat if the twisting of trajectory
_-> rdr d2r~l (d3r\ ,, ,
*2 = — * (113)
will be equal to zero in every point [3]. However, mixed multiplication of vectors in the right part of statement (1.1.3) becomes zero due to linear dependency (1.1.2), as the determinant of a matrix with linearly dependent rows.
In the particular case of ai(t) = 0 the order of equation is decreasing, trajectories becomes straight linear and under a(t) = 0 motion becomes inertial. Simplified variants of plane motion are discussed above. In the general case of plane motions we can consider a(t) = 1. We can consider a plane passing through the radius-vector connecting interceptor at initial time of persecution and the target at final time of persecution and vector of velocity of interceptor at initial time of persecution as a plane of persecution. Plane of persecution may change when adjusting the persecution.
Arbitrary coefficients of equation connected to the control signals and have to be defined by methods of reversal problems of dynamics. Solutions founded of conjugated equations (1.1.1) or (1.1.2) with some control signals a1(t), a(t), b(t) we substi-
tute into equation (1), from which we find the control signals u(t) controlled dynamical system [7].
1.2. Movement with minimal integral torsion. Planar motion of persecution is quite capable to solve the first problem of persecution and is not able to solve the second problem of persecution if the velocity vector at the time of completion of the persecution has a component orthogonal to the plane of the persecution.
In this case, a wider class of trajectories, such as converging spiral or helix - trajectory with torsion that goes out of the plane of persecution is needed.
To do this, we can choose the limitation in form of the generalization of statement (1.1.1) or the generalization of equation (1.1.2). We may require the minimal of functional
s=¡r^^m2 )ds.
(1.2.1)
c2 here is a second curvature of trajectory defined by formula (1.1.1) [3].
Using the orthogonality properties of the unit vectors of the accompanying benchmark, we can cause the integral (1.2.2) to the form
S = f (b(s)dn(s) + Ads)
(1.2.2)
Variation of functional (1.2.3) leads to the complex nonlinear equation interesting in some cases
/d2f d3r\ /_ J_ [is2 X ds3) V &2. d_ ( td3f(s) - ds\ V
ds3
X
dr(s) ds
) ( ff(s)2)
+
a(s)2 '
d2 (fdr d2r\ /_ J_\N - ds2 I \ls X ds2 ) I _ a2)
d J 4 d2f\ ,df ^ d2f. d3f. 22df -0 ds vcr3 ds2) ds ds2 ds3 ds
(1.2.3)
However, following a geometrical logic, as the conjugated dynamic equations we choose the generalized equation (1.1.3) of planar motion
r(t) = a (t) r(t) + b (t) r (t) + c (t) r (t) , (1.2.4)
which have three fundamental solutions. Two fundamental solutions gives a planar curative motion and the third fundamental solution may give a torsion of trajectory [8].
2. CLASSES OF PLANAR TRAJECTORIES OF PERSECUTION
2.1. Planar circular motion. The purpose of this and subsequent subsections is the solution of targeting problem at a particular class of trajecto-
ries. In case of constant speed planar motion the requirement of minimal integral curvature leads to the constant curvature round motion. This class of persecution is the second one after straight-line motion.
Algorithm is following. For a given initial velocity vector {Vx, Vy} find the line perpendicular to the initial velocity vector and passing through the initial point {X0, Y0}. For given initial {X0, Y0} and final {X1, Y1} points of trajectory of persecution define coordinates of the center of curvature of the trajectory {x, y} and the radius of curvature.
R =
HXI-ZXqX^XI+YI-ZY^O+Y2) (V}+V2x)
I (-VyXo+VyX^Y^x+YoVx)2
X =
_ -VyXI +VyX2 +Y2Vy-2Y1VyY0-2Y1VxX0+Y2Vy+2Y0VxX0 ~ l^-VyXo+VyX^Vx+YoVx) '
y =
_ 2Y0VyX0-2Y0VyX1-Y2 Vx+VxX2-2VxX0X1 +VxX2+VxY2
2(-VyX0+VyX1-Y1Vx+Y0Vx)
q) = ar ctan (--) .
(2.1.1)
Formula (2.1.1) solves the Cauchy problem and we can definitely determine the law of prosecution:
ut
x (t) = x + R c o s (— + q) ) ,
ut
y (0 = y + R co s (— + q ) , z(t) = 0.
2 R \X V 1
The time of persecution t = — a rc s i n ( x ° ) .
Here X - vector from initial to final point of persecution, V0 - vector of initial velocity.
This class may occur circular trajectories with a constant centripetal acceleration, but with variable velocity and curvature when it is necessary to catch up with the target or close to it with a minimum speed.
One of variant is R = K-u2 with constant K.
The second variant is the twisting spiral.
Initial spiral's slope is defined by vector of velocity at the start of persecution. Indeed let us suppose
a(t) = -2A,, b(t) = - A,2-ra2 and find the solution
x(t) = a-exp(-Àt)-cos(rat+À), y(t) = a-exp(-Àt)-sin(rat+À), z(t) = 0.
(2.1.2)
Parameter of the low of persecution are defined from initial coordinates and velocities of the target and interceptor.
Note for simplicity we use a relativity principle and we suppose initial point is stationary at the origin. Parameters of the law of persecution are calculated via relative initial coordinates and velocities
a = JXo+yo ,
x =
00 = —
-vyoyo-XoVX0
xo + Vi ' vx0y0-x0vy0
(2.1.3)
cp = ar ctan (—) .
2.2. Planar trajectories with available management. It has known in the problems of performance control signal take their boundary values from the area of management with an abrupt change [2]. Available control we understand as finite functions belong to the compact set of management U [2]. However, in contrast to the results of the theory of performance, our control signals are differentiable functions on the whole of set U.
Equation (1.1.3) with substitution V(t)=a(t)-Z(t) allows transition to equation
Z ( t)=f ( t)-Z ( t) ,
(2.2.1)
where f (t) = (b (t) —a ( t) j and function
}a(t)dt
o (t) = e 2 , which doesn't consist the first derivative. Among all kinds of functional dependencies the managing signal a(t) in form
a(t)
_ vN
v^Af ak ^ = 1 t-tk
(2.2.2)
is of particular interest, as it leads to expression
= YlLi (t - tky
ak
(2.2.3)
which represent the limited function if all ak are positive values. The numbers tk are the moments of time of shifting retarding to accelerating when zeroing relative velocity occur. These moments in time are the moments of finishing of persecution t. To complain this condition the coefficients of acceleration ak should be positive and possibly integer values. The equation (2.2.1) have a solution:
Z ( t) = C1-eM ( (2.2.4)
where function ^(t) is found as a solution of equation
j-J (t) + f (t) 2=f (t) .
(2.2.5)
Now we can limit the force of equation (2.2.1). Note the scale factor "force" Kf(t) is transferred to
the ^(t) and t: t1 = tK, ^ = that don't appear in
the integral (2.2.4) .
2.2.1 The first case. Let f (t) = 0. The cause by b (t) + —— a (t) = 0. For limb control signals and monotone subject that most closely approaching to Pontrjagin's "relay" type control, put
a(t)2 . ,.. n*(n+2) f . n*(n+2) —--a(t)=———, b(t) = — -
( t -t) 2
In this case, we find that
( t -t) 2
Z ( t) = c1( t + c 2 ) , V(t) = Ci (t-T)n (t + C2), xk(t) =
(Clk(t - T)n+1((n + 1) * t + x + c2k)
+ xok.
((n + l)(n + 2)) Or by beginning coordinates and velocity
XfcCO = (Oifc - *ok) (l + (n + 1) -V0 kt )(t-l ) (n+V + Xi k. (2.3.1)
The time of persecution have evolved as
x = —(-(n + l)Vx0 +
ax
+V ( (n + 1 ) Vx0) 2 + ax (n + 2 ) (n+ 1 ) (xx — x0) ).
(2.3.2)
2.2.2 The second case. Let b(t) = 0 . At that cause we have
Vk(t) =
2 J(n+1)2__ 2k
Unlike the first case, the relative velocity does not vanish while reaching the time of persecution t but it becomes the maximum physical and technical conditions for all positive integers even n, except for n = 0. When n = 0, we have uniformly accelerated motion with acceleration.
Ctyyi
■vk(t) =
Cifc(l
: /
2 C2... Vk (T) = C-f ,
Vfc ( 0 )=^— , Vk (t) = Vk (0 ) + ak^t .
For arbitrary n > 0 motion, as the velocity increases rapidly in the final stage of persecution. This situation can be compared with the slow persecution before the targeting in a straight line and then shooting a laser beam. This one has conditional practical application.
xk(t) =
\c2k)
= C
lk
: +
n+i+kn+iy-i
+ ■
/t-r^-J^+D2-!
\c2k)
\
3 .
- + n ■
2
+
(n + iy-\)
+Xo k
2.2.3. The third case. The control function b(t) such that the controlling force f (t), defined on a finite time interval t is given by Legendre polynomials
Pn(x) =
— • —((x2
2n-n\ dxn yy
I)"),
where = - . Those complex movements can be useful for the persecution of managed aircraft when an effort against the persecution takes place. In simple cases may be useful more simple management techniques, which are defined by special cases of solutions of equation (1.1.3).
2.3. The differentiable case. Let us b(t) =
=' a (t) , in equation (3), Then equation (1.1.3) reduces to
U (t) = a (t)U (t) + ñ ,
where U(t) is a vector of velocity and Q are integration constants. Let us denote
u (t) = e~$a (-) d then general type solution of the equation becames
x(t)
m m
nx Cx
ny •J u(t)dt+ Cy
az vcz\
u(t)
(2.4.1)
Existing of two vectors reduces this solution to the plane case. Solution (2.4.1) in the case of a constant a(t), shows that vector of velocity shifts from one direction to another one exponentially under persecution. Without losing of significance, we can assign
a (t) = —A, x(t)=^+Cxe
y(0=^+Cye-z (t) = 0 .
-At
■At
Parameters of the law of motion are express by coordinates and velocities of the initial an final points of persecution and time of persecution:
x ( t) = Vxlt + ^^e-Xt + Xo ,
y ( t) = Vnt + V-^e-^ + yo , 2 = VXoVyi-VXiVyo
XQ
Yo
-Vx1Y1+Vx1Y0+X1Vy1-X0Vy1 '
-Vx^yoXo + Vx\Yx - Vx, X, Vy,
+
+
Vx0 Vyx - Vx-^Vyo
Vx1X0Vy1 -vx^vxo +Vx1Y0Vx0 +X1Vx0Vy1
+
Vx0Vy1-Vx1Vy0 VyoX1Vy1-Vy0X0Vy1+Vy1Vx1Y1Vx1Y1
Vx0 — Vx-L Vy0 -Vy1Vx1Y0-X1Vy2+X0Vy2+Y0Vx0Vy1-Vy0Vx1Y1
+
T =
Vx0Vyx-VxxVy0
(2.4.2)
Time of persecution is estimated by formula:
VyoX-t-VyoXo+Vx^-Vx-CYo +
-X^y^XoVy^VXo+YoVXo
Vx0Vy1-Vx1Vy0 +
Vx0Vy1-Vx1Vy0
We have more general case of solution of the equation (1.1.3). Let us transform the equation (3) to the form
^(f(0 - a(t)f(t) + (d(t) - 6(0)r(0) =
= r ( t)± (a(t)-b(t)) .
and require should be constant.
Then we have equation
r (t) — a (t) r (t) + ar (t) = p ,
and it's solution under
a (t) = a • t ,
r(t) =
.,fa-a 3 at2\ , „ , T,fa-a 3 at2\ , „
- is a more general type of planar motion. There M and U are special functions.
3. PERSECUTION WITH TORSION
3.1. Screwed spiral. The most important in the theory of the prosecution motion with torsion is a movement along a twisting of spiral.
x 0
y (t) = R e xp ( — |rc) s in (o x) + y0, z(t) = azt2 + Vz0t + z0 . (3.1.1)
Parameters of the law of motion are found from initial and final points of the phase space as following. Initially we have to find radii of curvature R in the iteration process, R initial value is 10. Then we have to fined Z parameter as approximately solution of equation
- (x-l -*0) -tan(Z) +y0 -y1
R =
tan (Z)
(3.1.2)
Then we substitute Z found into the equation for correction of A:
A= exp
2ZVxc
vy0
A
(— tan (Z) x1 + tan(Z) x0 + y0 — yO2'
A = 2y1y0tan(Z)2 - y^tan(Z)2 + tan(Z)2x^ -—2tan(Z)2x1x0 + tan(Z)2Xo —
-yl + 2y1y0 - yl - y2tan(Z)2.
Then we do correction
2000
and hav-
ing a new value of Z we calculate a new correction to (3.1.3) formula
VZ0-VZ1
az = —■
2t
Control signals are calculated by formulas:
2t[i3 +3[i2+2t(ù2 h-Cl)2
a(t) = -6(t) =
t\i2+2\i+t(ù2
(n4+2n2ûJ2+ûJ4)t t[i2+2[i+t(ù2 '
(3.1.4)
c(t) =
y14+2[12(Ù2+(Ù4
t\i2 +2\i+ta>2
Where the functions a(t) and b(t) tend to a finite limit, and c(t) - tends to zero. Omitting cumbersome calculations show dependence of the rate of persecution time (Fig. 1).
We can see in this class of trajectories control signals are finite and vary only in the start of persecution until by torsion do not go on the plane prosecution (Fig. 2 and 3).
Fig 1. Dependence of velocity of persecution on the time
Fig. 2. Dependence of the torsion and curvature of the trajectory on the time of persecution: red line is torsion, green line is curvature
Fig. 3. Dependence of controlling forces on time
4. TRAJECTORIES OF PERSECUTION UNDER GLONASS CONTROL
4.1. Slowly moving objects. We divide the problem into two classes of objects: slowly moving bodies and fast moving bodies to select of trajectories of persecution. In case of slowly moving bodies deviations between real positioning and components of velocity vectors and their tabular values are small within GLONASS sessions receiving time At. Ships, cars, people, animals, birds etc. are slowly moving bodies. In case of slowly moving bodies, we can choice the quadratic law of motion - motion with constant acceleration. In the next session of communication, the acceleration may be abruptly slightly changed depending on value and direction of deviations. Let us assign
x = a0 + att + a2t2,y = 60 + btt + b2t2
(4.1.1)
For which equation we have
Z = C0 +1C2t2 + Cit.
The law of motion (4.1.1) has nine arbitrary constants, which we have to determine based on initial values. Let us suppose coordinates, velocities and accelerations qx, qy, qz, VX, Vy, Vz,
ax, ay az were calculated at the moment of previous session t1.
Rx ( t) = qx + Vxi+-ax t2,
Ry ( 0 = Ry + V yx +\ay t2, 1 ?
Rz(t) = V +-azt2.
(4.1.2)
While in this interval of time, the body should move under the law
Qx(t) = Qx + Ux-t + -wx-t2 Qy ( t) = Qy + Uy-t + jwy t2, Qz ( t) = Qz + Uz-t + \wz t2.
(4.1.3)
According to condition of persecution under t=0, a0=qx, b0=qy, z0=qz. At the moment x according to condition of persecution (1) and (2) we have equations:
% + 2a2x = Ux + wxx, x0 + % x + a2x2 =
1
= Qx + UXT + — wxx
from which we find
_ -2x0+2Qx+Uxr al — ~ > a2
-2x0-wxt2 +2QX
2t2
. (4.1.4)
Analogous,
h =
-2 y0 + 2 Qy + Uyt
b7 =
c[l] = c[2] = -
T }
—2y0 — wyr2 + 2 Qy
2T2 )
-2z0 + 2 Qz + UZ t
t
-2z0-wzt2+2Qz
(4.1.5)
Thus, using formulas (4.1.1) with coefficients (4.1.4) and (4.1.5) we can solve the problem of slowly moving bodies persecution. In this case, GPS/GLONASS receiver should accept current values of position, velocity and acceleration.
4.2. Fast moving objects. In the case of fast moving objects defining of quadratic on time law of motion is not enough. We should take into consideration variations of acceleration approximated by polynomials of third and fourth power. However in case of free of torsion planar trajectories we have relatively complicated third power polynomials analytical expressions making solving of targeting problem too hard. Indeed, in case of free of torsion planar motion with minimal curvature we found
x ( t) = a0 + axt + a2t2 + a3t 3 ,
y (t) = bQ + bit + b2t2 + b3t3 ,
z(t) = z0 + qect - r)y(t) +
+C2e(t - t) f eF^dt
(4.1.5)
where
F = —2W4 -
-2
-2
— 6 agôf + 6a3b3a2b2 — + 9(2301% + 6 agôf — 6a3b3a2b2
M1
W. -
W-, -
1
' -9a\bl - lQal^a^ + ^ +18a303a102 — 18a3a1ô|a2 + + 12a3a2b2b3 + 18a|ô1ô2 + 18 alb^a-t — 12a3a2btb3 — 12a|a1ô| — 12 a3a2b2 — \ —9 a3alb3 + 12a3a1Ô203a2
M3W3
/
-6—(-6a36|ai6i + 3albl ~
-4albl - Aasb^b^ --4a1ô2a2ô32 + 3a3b3b^ + +180301«! - 18alèfèi + +12a3b2b3 + 12ô|a1a2 — -8 bla22b2 + 4ô32a22ô1 -—12 b2a3a2b3 + 4a1ô|a3ô3 + +1 2a3b3a2b2 - 1 2a3bla±b2) M2. (4.1.6)
W values are a functional dependencies, M values are numerical expressions of a[£], b[k] constants.
M1 = 3 a\bl + 3 60,3(1^2 + 12 a16|a2 — — 12a3b3a1b2 — 6a3b1b3a2 — 4a1b2b3a2 +
+4a3a16| - 8a2b2b3 + Aajb^ +
+12a3b3a2b1 — 24a3b3a2b2 + 12a3b2 +
+3 a 2b2 — 12 a'2b1b2 + 8 a 3a2b2 — 4 a3b2a2 bx ,
-1
M2 = V3ô!ô3 - b2
9b2
M3 = ^¡Зa^x3 — a^ ,
Wt = arctan (M2(b2 + 3 b3t + 3 b3)) ,
W2 = 1 n (at + 2a2t + 3 a3t2) ,
W 3 = arctan (M3(a2 + 3a3 t)) ,
W.4 = 1 n (b 1 + 2b2t + 3 b3t2) ,
W5 = 1 n (b! + 2 b2t + 3 b3t2 - 2b2 - 6 b3t) .
Integral parameter in formula (4.1.5) is specific for motion with variable acceleration. M1 value at denominator of formula becomes equal to zero under zero values a3 and b3, which are responsible for variable acceleration. Solutions of the integral are elementary functions, thus we can find twelve constants ak, bk using algebraic sums. The simplest maneuver of acceleration allow perform transition from one predefined trajectory to another one within a finite interval of time. In this case, both initial positions and velocities and final ones are coincide. Using these parameters, we can find these twelve
constants of the law of motion with above-mentioned formulas. The problem of management of fast moving objects using GPS/GLONASS is that positioning is performing based on quasi range from receiver to constellation within cone of visibility limited by angle of 5-10 degrees above horizon. Lower angle rays cause multipath interference, which decreases accuracy of positioning.
However in case of fast moving objects the acquiring of satellite signals will occur within the sequence of points Mj(x(ti), y(ti), z(t,)) not in the one point of trajectory. We find the coordinates and velocities related to the last point along the path defined by the law of persecution {x(t), y(t), z(t)}. We correct pseudo ranges of previous points with cosine theorem. Thus, we eliminate dilution of precision of positioning of fast moving objects, which may exceed hundreds of meters. In our previous work, we show more precise method of positioning using analytical approach [1]. The approach allow find more accurate trajectories of persecution between communication sessions. Thus, the problem becomes interconnected and under certain conditions, to be special studies it is becomes convergent to the exact limit.
5. ALGEBRAIC CLASSIFICATION OF TRAJECTORIES OF PERSECUTION
Geometrical classification of trajectories of persecution is based on equations describing trajectories of persecution as lines that bend elastic thread with appropriate boundary conditions. There is alternative. We can consider surfaces emerging under deformation of the elastic film and the congruence of geodesic lines on curved surfaces such as the classes of trajectories persecution. In addition, as the world lines on these surfaces, which define the torsion of the Lorentz force. It is known that the equations of Einstein's gravity describe four-dimensional space-time continuum as pulsating elastic film in the space of higher dimension. Therefore, we consider the equation of motion in spaces of general relativity in the form of the Ham-ilton-Jacobi equation, as a model equation. It is known algebraic classification of these equations based on the algebra of first integrals of motion. This method can be generalized for other Hamilto-nian dynamical systems, specifically for the Pontryagin's Hamilton in modern theory of management of dynamic motion.
5.1. Classification of solutions of Hamilton-Jacobi equation
Let us consider the Hamilton-Jacobi equation of test charged particle in general relativity theory
as a model equation. The arbitrary gravity field in conjunction with electromagnetic field are able to form various combinations of curvature and torsion of tested charged particles trajectories. It is known the equation of Hamilton-Jacobi in general theory of relativity may be presented as following [4].
vik(x) (¿-ks—qA 0 (£ts—qA 0=m 2,
(5.1.1)
where xi are generalized coordinates, gik(x) are metric coefficients expressed as coordinate functions, S
is function of action, p [k] = S(x,x [0 ] ) is gen-
"xk
eralized momentum, qAk(x) are charge and vector -potential of the electromagnetic field. It is known also first integrals oa=oa(x, p) of Hamilton's system form the Lie algebra relatively Jacobi brackets.
First integral commutator can be expressed linearly from integrals of algebra with structural constants of algebra which are anti symmetric on lower indexes and satisfy to the Jacobi identity.
{°a> °b} — £r=i Ca,b°r
(5.1.2)
{ac, { a a, ab}} + {a a, { ab, oc}} + {ab, { a a, ac} } = 0 .
(5.1.3)
Hamiltonian
H(P, x) = gik[-£-S- qAk) (4-S- qA,
dxu
dx¡
satisfy to Jacobi brackets because it is integral of motion also.
{H (P,X), O a (p,x) } = 0 .
(5.1.4)
Algorithm of composing of classes of trajectories is following.
1. The system of homogeneous quadratic equations relatively unknown components of algebra (5.1.3) structural constants containing n linear first integrals of motion is solving.
oa = Ja(x)lpl, a = 1...n, i = 1.4, n = 1.10 . (5.1.5)
Based on the found coefficients structural constants excepting isomorphic cases are composing.
2. For each algebra defined by structural constants crb the system (5.1.2) is solving and first integrals oa = Ja(x) lpi or Kiling's vectors Ja(x)l, a = 1...n are fining.
3. System of equations (2.1.4), which decomposes for linear integrals (2.1.5) into a system of
equations for the components of the Killing metric gik(x) is solving:
D J fe (x) + D „] i(x) = О ,
(5.1.6)
where Di denotes the covariant derivative with respect to the desired metric and to the Kiling equation for the vector component of electromagnetic field:
Ja(xy^(Aj(x)g(xym)-
—Aj (x) g (x) jl^(Ja (x) m) = 0 (5.1.7)
Solving systems of linear partial differential equations (5.1.6) and (5.1.7) we restore the power functions of the each algebraic class with more or less certainty, depending on the number of integrals of algebra n, and with them, and Hamiltonian.
4. As a result, availability of composing of algebraic classes of trajectories is getting possible. There are two approaches. The first one is a direct method. Coordinates on the deformed film, we transfer on locally planar space. Because analytical expression of the law of motion of a certain algebraic class is desirable. To find the law of motion we compose an equation of geodetic line - trajectory without torsion or we compose an equation of world line in present of electromagnetic field - trajectory with torsion not transferring the body from surface of elastic space. In the second case, we solve the problem of embedding a deformed film in the planar space of higher dimension and define a class of trajectories in planar space.
5. Thus, it is possible to solve the problem of persecution in this algebraic class of trajectories. To do this, we solve the Cauchy problem of the transformed equation (5.1.8) at the start and final data of trajectory of persecution. A number of authors solved that hard problem. Levi-Chivita found structural constants of a three membered Lie algebra [5]. A. Z. Petrov and V. R. Caigorodov at Kazan University, Russia found all remaining algebraic structures with real structural constants until tenth order. These authors found canonic forms of Kiling vectors respecting to linear first integrals and satisfying those metrics [6]. One of author of this researching found all electromagnetic fields on Kiling-Petrov-Kaigorodov's vectors [7]. For commutative algebras of first integrals, the problem was solved for classical system by Jarov-Jarovoi and for general case of general relativity was solved by Jarov-Jarovoi's method [8].
Using Maple environment for computer aided mathematical simulation the solution of the problem is possible in general form during a short time
not only in canonical Petrov-Kaigorodoc-Zakharov's form. Authors have a number of cases proving the method. Let us consider a case of filth order algebra.
A metric tensor
9ij =
-10 0 0 0 -1 -x2-x0 0 -x -10 L 0 0 0 1J
satisfy algebra of integrals
Ol = P 1— УР2 , а2=Р2 , аз = Р 3 , 04=P4 ,
ö"5 = "УР 2 + +
7 7
у -ДГ
P 3 .
Geodetic lines equations
x"'(t) - x(t)y'(t)2 - y'(tV(t) = 0,
z"(0 + 2
have solutions:
y"(t) + x(t)x'(t)y'(t) + y'(t)z'(t) = o, x(t)2 1
+ -)x'(t)y'(t)-x(t)x'(t)z'(t) = 0, Z'(t) = О
1. Rectilinear uniform motion along the axis ox or axis oz;
2. Uniform motion along the axis oy and uniformly accelerated motion along the axis oz;
3. Uniform motion along the axis ox and oy, uniformly accelerated along the axis oz.
4. The variant that met in Part 1:
y (t) = — -fc — +
_L £5 2есзесз
x (t) = S^mmy(-)d t + xo ,
z ( t) = j — y(-)+x(-)*(-)y(-)dt + zo . v y J x(t) u
5.2 Algebraic classification of trajectories of persecution by Pontryagin
It is known the Pontryagin's theory of management motion is reduced to the Hamilton form [2].
H(xp,x,u) = Zn+1ipaf(x,u)a . (5.2.1) дх1 _ д dVi _ а , .
= , ~-г = — -г" . (5 22)
Moreover, generalized momentum is a gradient of some function of action, which leads to the form of the Hamiltonian system of Hamilton - Ja-cobi equations [2]. Let us transfer the general relativity algorithm of algebraic classification of Ham-ilton-Jacobi equations to the Hamiltonian system of
40
mathematical modeling and decision-making
controlled motion (5.2.1), (5.2.2). In case of planar motion, we have four-dimensional space of variables xi - two generalized coordinates and two generalized velocities. For this case, structural constants of Lie algebra and respecting first integrals are transferring without changing from researching of Petrov and Caigorodov [4]. Kiling type equations for the Hamiltonian (5.2.1) of autonomous system with r-membered Lie algebra of first integrals
°a=Ja(x)l, a = 1...r, i = 1...H+1.
can be presented as d
71+1
dxi
jfl(x, u)
7 = 1
n+1
7=1
i,j = 1 ... (n + 1), a = 1 ... r.
(5.2.3)
Killing vectors belonging to the Lie algebra ensures consistency of the system (5.2.3) and the existence of solutions for the required power functions , which becomes functions of coordinates of extended phase space [9]. This situation is typical for autonomous dynamical systems, when the control signals are not functions of time, and the functions of the coordinates, and they have entered the control forces implicitly [2]. We have to find the law of motion of managed object and control signals based on found power functions - right sides of equations (1). For movements in three-dimensional configuration space, we have the six-dimensional phase space and the seven-dimensional extended phase space. As far as authors know algebraic structure these spaces are not currently constructed.
Let us consider for instance the planar motion in xoy plane. Let assign x=xi, y=x2, x'(t) = x3, y'(0 = X4■
We will looking for power forces in form of
fl = X3>f2 = X4'/3 = X4>
0 (xs),h = - X30 (xs),fs = 0 (xs) .
Let us suppose that first killing integral
Jl = x3> J2 = x4, /3 = w ' x4' J4 = —w ' x3> J5 = Fs(x3,xs) is known. Thus equation (5.2.3) gives equation . And equations of persecution
becomes equations of performance with trajectories
x(t) = -
C3 COs(ûJt)-C4 sin(ûJt)-C2ûJ
Vx(t) = C3 Sin (<Dt) + C4COS (COt), Vy(t) = C3 cos(ct) — C4sin (ct), r = t + Cs.
CONCLUSIONS
1. Trajectories of persecution satisfying to the condition of minimal curvature and torsion may differ significantly from piecewise time-quadratic or time-cubic laws of performance.
That may be a decisive factor for the solving of problems of persecution. Fig. 4 show above mentioned.
y(0 =
C3 sin(ûJt)+C4 COsCûJiO + CiûJ
b
Fig. 4. The trajectory of persecution as a line with minimal curve (a) and the trajectory of persecution as a cubic relationship (b), t2 = 6,2033
Initial data in cases of A and B are the same but the final conditions are differ significantly. Curves are coincide at the initial stage and are disperse outside curve of minimum curvature decomposition in power series of Taylor (Fig. 5).
2. We do not need finding the quadratic laws of persecution motion at every small region of persecution because computational load will increase
significantly and both calculation errors and iteration errors will increase extensively. We can use a certain class of laws of motion at the finite trajectory of persecution within a finite time instead. Using this approach, we should do decomposition of trajectories into classes. A geometrical approach containing the set of simplest geometry trajectories could not give us comprehensive and invariant classification of trajectories of persecution. Algebraic classification is better for solving the problem of classification of trajectories.
3. Uniformly accelerated law of motion is using in the theory of performance. Complex trajectory is approximating by splines with appropriate triggering of control signals. In the theory of persecution a set of predefined trajectories is using. Each class of trajectories have own analytical expression and constants as equation parameters. The problem of targeting is solving separately for each class of equations. These solutions include time of persecution and numerical values of trajectory parameters allowing transition from initial point of phase space to the final one. The law of motion with certain numerical parameters allow find control signals as continuous functions of time. Algebraic classification of trajectories of motion is very valuable because it is general and invariant.
REFERENCES
1. Захаров А. В., Богданов М. Р., Мухаметъянов А. Р.
Аналитический метод решения задачи позиционирования // Информационные технологии. 2012. № 8. С. 63-65. [ A. V. Zakharov, M. R. Bogdanov, A. R. Mukhametjanov, "Analytical method for solving of problem of positioning," (in Russian), Informazionnye technologii, no. 8, pp. 63-65, 2012. ]
2. Понтрягин Л. С., Болтянский В. Г., Гамкрелид-зе Р. В., Мищенко Е. Ф. Математическая теория оптимальных процессов. 4-е изд. М: Наука, Гл. ред. физ.-мат. лит., 1983. 392 с. [ L. S. Pontryagin, V. G. Boltyansky, R. V. Gam-krelidze, E. F. Mishenko, The Mathematical Theory of Optimal Processes, (in Russian). Moscow: Nauka, 1983. ]
3. Балонин Н. А. Новый курс теории управления движением. СПб.: СПбУ, 2000. 160 с. [ N. A. Balomin, The new course of theory of control of motion, (in Russian). St. Petersburg: Publishing of Saint Petersburg University, 2000. ]
4. Смирнов В. И. Курс высшей математики. Т. 2. М.: Наука, 1972. С. 396-402. [ V. I. Smirnov, The course of high math, (in Russian), vol. 2, pp. 396-402. Moscow: Nauka, 1972. ]
5. Красовский Н. Н. Управление динамической системой. Задача о минимуме гарантированного результата. М.: Наука. Гл. ред. физ.-мат. лит., 1985. 520 с. [ N. N. Krasovsky, Control of dynamic system. The problem of minimal of guaranteed result, (in Russian). Moscow: Nauka, 1985. ]
6. Петров А. З. Новые методы в ОТО. М.: Наука, 1966. [ A. Z. Petrov, Novel methods in the general theory of relativity, (in Russian). Moscow: Nauka, 1966. ]
7. Захаров А. В., Сингатуллин Р. С. Структура обобщенного импульса пробной заряженной частицы и обрат-
ная задача в ОТО // Известия вузов. Физика. 1981. № 7. С. 60-64. [ A. V. Zakharov, R. S. Singatullin, "The structure of generalized momentum of test charged particle and reversal problem in the general theory of relativity," (in Russian), Izvestia vuzov, no. 7, pp. 60-64, 1981. ]
8. Захаров А. В. Решение обратной задачи динамики пробной заряженной частицы в ОТО по алгебрам первых интегралов: автореф. дис. ... канд. физ.-мат. наук. Уфа: БГПУ, 1998. [ A. V. Zakharov, Solving of reversal problem of reversal problem of dynamics of test charged particle in general theory of relativity based on first integrals algebras, (in Russian). Ufa: BSPU, 1998. ]
9. Понтрягин Л. С. Непрерывные группы преобразований. М. Наука. Гл. ред. физ.-мат. лит., 1982. 402 с. [ L. S. Pontryagin, Continuous groups of transformation, (in Russian). Moscow: Nauka, 1982. ]
ABOUT AUTHORS:
ZAKHAROV, Aleksander Vasilievich, Assistant Professor, Institute of Professional Education and Information Technologies. BOGDANOV, Marat Robertovich, Assistant Professor, Institute of Professional Education and Information.
МЕТАДАННЫЕ
Название: Классификация траекторий преследования.
1 2 Авторы: А. В. Захаров , М. Р. Богданов
Организации:
Башкирский государственный педагогический университет им. М. Акмуллы (БГПУ), Россия. Башкирский государственный университет (БГУ), Россия.
Email: 1 ralzax@yandex.ru, 2 bogdanov_marat@mail.ru. Язык: английский.
Источник: Вестник УГАТУ. 2014. Т. 18, № 5 (66). С. 31-41.
ISSN 2225-2789 (Online), ISSN 1992-6502 (Print). Аннотация: Предложена классификация законов движения преследования. Проведено исследование режимов преследования в ряде конкретных случаев. Рассмотрены варианты управления движением при использовании систем ГЛОНАСС и GPS. Ключевые слова: управление динамическим движением; траектории преследования; GPS; ГЛОНАСС.
Об авторах:
ЗАХАРОВ Александр Васильевич, доц. каф. информационных и полиграфических систем и технологий БГПУ, доц. инж. фак-та БГУ. Дипл. физик (БГПУ, 1973). Канд. физ.-мат. наук (там же, 1984). Иссл. в обл. алгебраической классификации решений обратной задачи динамики в ОТО.
БОГДАНОВ Марат Робертович, доц. каф. прикладной информатики. Дипл. преп. биол. и хим. (БГУ, 1991). Канд. биол. наук (Всерос. ин-т защ. растений, 1997). Иссл. в обл. цифровой обработки сигналов, геогр. инф. систем, биоинформатики.
