Next Article in Journal
On General Reduced Second Zagreb Index of Graphs
Next Article in Special Issue
Improved Hille-Type Oscillation Criteria for Second-Order Quasilinear Dynamic Equations
Previous Article in Journal
Numerical Computation of Hybrid Carbon Nanotubes Flow over a Stretching/Shrinking Vertical Cylinder in Presence of Thermal Radiation and Hydromagnetic
Previous Article in Special Issue
Numerical Solutions of Fractional-Order Electrical RLC Circuit Equations via Three Numerical Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Time-Optimal Control Problem of Two Non-Synchronous Oscillators

Institute of Control Sciences of RAS, 117997 Moscow, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(19), 3552; https://doi.org/10.3390/math10193552
Submission received: 22 August 2022 / Revised: 23 September 2022 / Accepted: 26 September 2022 / Published: 29 September 2022
(This article belongs to the Special Issue Mathematical Modeling and Simulation of Oscillatory Phenomena)

Abstract

:
The time-optimal control problem for a system consisting of two non-synchronous oscillators is considered. Each oscillator is controlled with a shared limited scalar control. The objective of the control is to accelerate the oscillatory system to a given specific position, where the first oscillator must have non-zero phase coordinates, but the second one must remain motionless at the terminal moment. For an arbitrary number of unknown switching moments that determine the optimal relay control, the necessary extremum conditions in the form of nonlinear matrix equalities are proposed. The study of the necessary/sufficient conditions of the extremum made it possible to describe the reachability set in the phase space of the first oscillator, to find an analytical form of the curve corresponding to the two-switching control class, which also separates the reachability set of the three switching-control class. The corresponding theorems are proved and the dependence of the criteria on control constraints is shown. Matrix conditions for different classes of control switchings are found. All of the obtained analytical results are numerically validated and illustrated with mathematical modeling.

1. Introduction

Problems surrounding the lack of control resources (when the dimensions of a control vector are less than the dimensions of the state space of a physically-controlled system) have widespread applications in practice. Oscillatory systems, such as electric circuits [1], mechanical systems [2,3], quantum oscillators [4,5], systems involving physical nature [6,7], and controlled with one external force, are instances of such systems. Using the oscillatory system, it is possible to investigate periodic solutions and the stability of the model of two coupled oscillators in the optics of chiral molecules [6], as well as consider the behaviors of biologically complex mechanisms [7]. Increasing the number of oscillators in these models to sufficiently large N makes it possible to study (analytically and numerically) the macro-parameters and properties of similar systems [8,9].
Thus, the Timoshenko–Ehrenfest beam motion control, which is modeled with a system of N oscillators that perform forced oscillations under the actions of scalar and limited control forces, was considered in work [2]. The time-optimal problem for a platform with attached oscillators was investigated in [3]. The maximum variation of the oscillation energy of the oscillator system for a given time was examined in [10]. Due to the difficulties in finding an analytical solution for the control problems, it is of interest to asymptotically synthesize the optimal control of a system of an arbitrary number of linear oscillators (with general limited control), as obtained in [11].
Optimal control problems were investigated using the Pontryagin maximum principle, which was formulated during the seminars on the theory of oscillations and automatic control by L.S. Pontryagin and M.A. Aizerman in the 1950s, where the problem of optimal performance was considered [12]. For systems considered in [1,5], it is often necessary for one of the subsystems to get into the required position as quickly as possible, while the others should remain at rest at the terminal moment. Thus, for a single oscillator, V.G. Boltyansky obtained the synthesis of optimal control in the time-optimal problem [13]. The optimal control problems for the systems of many oscillators with various criteria have already been studied by F.L. Chernousko in [14], where a significant complexity of obtaining an analytical solution was noted, and the problem involving two oscillators was solved for a very special case of their frequency ratio.
In this paper, we review the time-optimal control problem for two non-synchronous oscillators controlled by limited scalar impacts with initial zero-phase coordinates but different terminal ones, which are features of the problem statement that considerably complicate it.
One of the key aspects of such control problems is the study of the accessibility and controllability of the oscillatory system [15]. It can be shown with the methods of geometric control theory [16], namely: the Sussman–Djurdjevic and the Poincare theorems, as well as with the results of the classical optimal control theory, the La Salle–Conti theorem [17].
This article is structured as follows. In Section 2, the optimal control problem is formulated. In Section 3, the optimal control law is presented using the Pontryagin maximum principle. Section 4 examines the necessary and sufficient optimality conditions for different classes of control switchings. Section 5 discusses the continuity of the problem criterion on the control parameter. Section 6 illustrates the obtained results, in particular, the mathematical modeling of the set of reachability of the first oscillator, the criterion of the problem, and the time intervals of the switching controls. In the conclusion, we suggest directions for further research.

2. Optimal Control Problem Statement

A system consisting of two non-synchronous oscillators, i.e., having different natural frequencies ( ω 1 ω 2 ) , was considered. The control u ( t ) impacts the momenta p 1 , p 2 of both oscillators, but not their coordinates q 1 , q 2 . The motion of this linear control system is described with the following system of equations:
q ˙ 1 ( t ) = p 1 ( t ) , p ˙ 1 ( t ) = ω 1 2 q 1 ( t ) + u ( t ) , q ˙ 2 ( t ) = p 2 ( t ) , p ˙ 2 ( t ) = ω 2 2 q 2 ( t ) + u ( t ) , x = ( q 1 , p 1 , q 2 , p 2 ) T R 4 = M .
The system (1) can be rewritten in matrix form x ˙ ( t ) = A x ( t ) + B u ( t ) , where
A = 0 1 0 0 ω 1 2 0 0 0 0 0 0 1 0 0 ω 2 2 0 , B = 0 1 0 1 .
The number of non-zero components of the control vector B is less than the dimensionality of the system; moreover, the control is limited by a value ε , which is a parameter of the problem:
u ( t ) [ ε , ε ] = U .
At the initial moment, the phase states of both oscillators coincide with the origins of the phase coordinate planes, i.e., their coordinates and momenta are equal to zero. At the terminal moment, the phase coordinates of the second oscillator become zero again, while the first oscillator moves to the state ( q 1 T 0 , p 1 T 0 ) , which are also parameters of the problem. The explicit forms of the initial and terminal conditions are
x ( 0 ) = x 0 = ( 0 , 0 , 0 , 0 ) T ,
x ( T 0 ) = x T 0 = ( q 1 T 0 , p 1 T 0 , 0 , 0 ) T ,
where T 0 is a terminal moment and T defines a transpose operation.
The total time of the system’s movement from the initial state to the terminal one is chosen as the criterion of the problem and needs to be minimized:
T 0 [ u ] = 0 T 0 d t min u ( · ) U .
The proofs of global controllability and strong accessibility of the system (1) with the condition (2) are given in [15,18].

3. Solution of the Optimal Control Problem

To investigate and solve problems (1)–(5), the Pontryagin maximum principle (PMP) is used. The maximum principles in the canonical coordinates consist of the following statements:
  • The Hamiltonian for the optimal control problem is expressed as
    h u = λ 0 + ξ 1 p 1 + ξ 2 p 2 + η 1 ( ω 1 2 q 1 + u ) + η 2 ( ω 2 2 q 2 + u ) , λ 0 > 0 ,
    where ξ 1 , ξ 2 , η 1 , η 2 are adjoint variables and λ 0 is constant.
  • A Hamiltonian system, including motion equations and an adjoint system of equations:
    q ˙ 1 ( t ) = p 1 ( t ) , p ˙ 1 ( t ) = ω 1 2 q 1 ( t ) + u ( t ) , q ˙ 2 ( t ) = p 2 ( t ) , p ˙ 2 ( t ) = ω 2 2 q 2 ( t ) + u ( t ) , ξ ˙ 1 ( t ) = h u q 1 = ω 1 2 η 1 ( t ) , η ˙ 1 ( t ) = h u p 1 = ξ 1 ( t ) , ξ ˙ 2 ( t ) = h u q 2 = ω 2 2 η 2 ( t ) , η ˙ 2 ( t ) = h u p 2 = ξ 2 ( t ) .
  • Maximum condition
    h u = max u ( · ) U h u = λ 0 + ξ 1 p 1 + ξ 2 p 2 η 1 ω 1 2 q 1 η 2 ω 2 2 q 2 + u ( η 1 + η 2 ) .
Transversality conditions are absent because the optimal problem statement involves fixed endpoints.
The Hamiltonian (8) is linear with respect to u and, thus, the optimal control is given by the bang–bang control: u = ε for η 1 + η 2 < 0 and u = ε for η 1 + η 2 > 0 , which can be written as
u ( t ) = ε sign ( S F ( t ) ) ,
S F ( t ) = η 1 ( t ) + η 2 ( t ) ,
where S F ( t ) is a switching function. The adjoint system can be rewritten in the following form
η ¨ 1 ( t ) = ω 1 2 η 1 ( t ) , η ¨ 2 ( t ) = ω 2 2 η 2 ( t ) .
From the system (11) η 1 ( t ) and η 2 ( t ) can be found:
η 1 ( t ) = C 1 cos ( ω 1 t ) + C 2 sin ( ω 1 t ) , η 2 ( t ) = C 3 cos ( ω 2 t ) + C 4 sin ( ω 2 t ) ,
where C = C 1 , C 2 , C 3 , C 4 is a constant vector. Combining the expressions (9) and (12), an explicit form of the optimal control law can be obtained
u ( t ) = ε sign ( C 1 cos ( ω 1 t ) + C 2 sin ( ω 1 t ) + C 3 cos ( ω 2 t ) + C 4 sin ( ω 2 t ) ) .
This function is uniquely defined and cannot be equal to zero on the whole interval, perhaps with the exception of the isolated points, which lead to the absence of special control modes.
The relay form of the optimal control obtained in this Section allows explicitly writing solutions of the motion Equation (1) taking into account how many times the control (13) will change the sign and obtain the necessary extremum conditions for any number of such switchings.

4. The Necessary Extremum Conditions

For the relay control (13), two values are defined: the moment of switching t i , i = 1 , K 1 ¯ , K N , when the control changes signs and the duration of control intervals τ n , n = 1 , K ¯ . Here, a further designation 1 , K ¯ , which means “ 1 , 2 , , K ”, is used for convenience. The optimal control u ( t ) with K 1 switchings and K control intervals is illustrated in Figure 1:
Definition 1.
The term “ ( K 1 ) -switching control class” defines an optimal solution where the optimal control has exactly K 1 switchings.
The solution to the system of motion Equation (1) with boundary conditions (3) and (4) has a general form
q 1 ( T 0 ) = 1 ω 1 0 T 0 sin ( ω 1 ( T 0 τ ) ) u ( τ ) d τ = q 1 T 0 , p 1 ( T 0 ) = 0 T 0 cos ( ω 1 ( T 0 τ ) ) u ( τ ) d τ = p 1 T 0 , q 2 ( T 0 ) = 1 ω 2 0 T 0 sin ( ω 2 ( T 0 τ ) ) u ( τ ) d τ = 0 , p 2 ( T 0 ) = 0 T 0 cos ( ω 2 ( T 0 τ ) ) u ( τ ) d τ = 0 .
The structure of optimal control (13) allows to write a solution of the system that takes into account not only the number of switchings K 1 , but also the sign of the control at the first interval.
2 j = 1 K ( 1 ) j + 1 cos ω 1 i = j K τ i cos ω 1 i = 1 K τ i = ( 1 ) K 1 + ( 1 ) k + 1 q 1 T 0 ω 1 2 ε , 2 j = 1 K ( 1 ) j + 1 sin ω 1 i = j K τ i sin ω 1 i = 1 K τ i = ( 1 ) k p 1 T 0 ω 1 ε , 2 j = 1 K ( 1 ) j + 1 cos ω 2 i = j K τ i cos ω 2 i = 1 K τ i = ( 1 ) K 1 , 2 j = 1 K ( 1 ) j + 1 sin ω 2 i = j K τ i sin ω 2 i = 1 K τ i = 0 ,
where an additional parameter k equal to 0 or 1 corresponds to the initial control ε and ε , respectively.
Remark 1.
A system (15) consisting of four equations can be used to study the three-switching control class, where four control intervals need to be determined.
The two-switching control class occurs as degeneration of other classes. For example, when the outermost interval is zeroed in on a three-switching control class, or the internal interval is zeroed in on a four-switching control class.
Lemma 1.
In the two-switching control class, the functional dependencies for the duration of the control intervals τ 1 , τ 2 , τ 3 are valid
cos ( ω 2 τ 1 ) = cos ( ω 2 τ 3 ) ,
cos ( ω 2 τ 2 ) = 2 cos 2 ( ω 2 τ 1 ) 4 cos ( ω 2 τ 1 ) + 3 5 4 cos ( ω 2 τ 1 ) i f τ 2 < 2 π ω 2 .
The proof of Lemma 1 is given in Appendix A.
To investigate control classes with more than three switchings, the necessary extremum conditions are introduced.
Theorem 1
(Necessary Extremum Conditions). Any solution to the problem (1)–(5) in the class of piecewise constant controls (13) satisfies the system (15) and K 4 equalities:
det cos ( ω 1 t i ) cos ( ω 1 t i + 1 ) cos ( ω 1 t i + 2 ) cos ( ω 1 t i + 3 ) sin ( ω 1 t i ) sin ( ω 1 t i + 1 ) sin ( ω 1 t i + 2 ) sin ( ω 1 t i + 3 ) cos ( ω 2 t i ) cos ( ω 2 t i + 1 ) cos ( ω 2 t i + 2 ) cos ( ω 2 t i + 3 ) sin ( ω 2 t i ) sin ( ω 2 t i + 1 ) sin ( ω 2 t i + 2 ) sin ( ω 2 t i + 3 ) = 0 , i = 1 , K 4 ¯ .
Proof. 
The switching function (10) that determines the optimal control (13) becomes zero at moments t i , i = 1 , K 1 ¯ , when control switchings occur:
S F ( t i ) = C 1 cos ( ω 1 t i ) + C 2 sin ( ω 1 t i ) + C 3 cos ( ω 2 t i ) + C 4 sin ( ω 2 t i ) = 0 , i = 1 , K 1 ¯
For a more compact notation, Equation (17) can be rewritten in matrix form
( C , Ω i ) = 0 , i = 1 , K 1 ¯ ,
where
C = C 1 , C 2 , C 3 , C 4 , Ω i = cos ( ω 1 t i ) sin ( ω 1 t i ) cos ( ω 2 t i ) sin ( ω 2 t i ) .
Equation (18) form the system consisting of K 1 equations
cos ( ω 1 t 1 ) cos ( ω 1 t 2 ) cos ( ω 1 t K 2 ) cos ( ω 1 t K 1 ) sin ( ω 1 t 1 ) sin ( ω 1 t 2 ) sin ( ω 1 t K 2 ) sin ( ω 1 t K 1 ) cos ( ω 2 t 1 ) cos ( ω 2 t 2 ) cos ( ω 2 t K 2 ) cos ( ω 2 t K 1 ) sin ( ω 2 t 1 ) sin ( ω 2 t 2 ) sin ( ω 2 t K 2 ) sin ( ω 2 t K 1 ) T C 1 C 2 C 3 C 4 = 0 .
The condition of the existence of the non-zero vector C in (19) is equivalent to K 4 equalities
det Ω i , Ω i + 1 , Ω i + 2 , Ω i + 3 = 0 , i = 1 , K 4 ¯ .
In other words, (20) are the non-degeneracy conditions of vector C, defining the optimal control (13). □
Remark 2.
In the case when K 1 = 4 , the condition of the existence of non-zero vector C has the form
det cos ( ω 1 t 1 ) cos ( ω 1 t 2 ) cos ( ω 1 t 3 ) cos ( ω 1 t 4 ) sin ( ω 1 t 1 ) sin ( ω 1 t 2 ) sin ( ω 1 t 3 ) sin ( ω 1 t 4 ) cos ( ω 2 t 1 ) cos ( ω 2 t 2 ) cos ( ω 2 t 3 ) cos ( ω 2 t 4 ) sin ( ω 2 t 1 ) sin ( ω 2 t 2 ) sin ( ω 2 t 3 ) sin ( ω 2 t 4 ) = 0 ,
which is equivalent to the following equation
sin ( w 2 τ 4 ) sin ( w 1 τ 2 ) sin ( w 2 ( τ 3 + τ 4 ) ) sin ( w 1 ( τ 2 + τ 3 ) ) sin ( w 2 ( τ 2 + τ 3 ) ) sin ( w 1 ( τ 3 + τ 4 ) ) + + sin ( w 2 τ 3 ) sin ( w 1 ( τ 2 + τ 3 + τ 4 ) ) + sin ( w 2 ( τ 2 + τ 3 + τ 4 ) ) sin ( w 1 τ 3 ) + sin ( w 2 τ 2 ) sin ( w 1 τ 4 ) = 0 .
If K 1 > 4 , there are C K 1 4 (number of permutations) of the non-zero conditions for the vector C
det cos ( ω 1 t m 1 ) cos ( ω 1 t m 2 ) cos ( ω 1 t m 3 ) cos ( ω 1 t m 4 ) sin ( ω 1 t m 1 ) sin ( ω 1 t m 2 ) sin ( ω 1 t m 3 ) sin ( ω 1 t m 4 ) cos ( ω 2 t m 1 ) cos ( ω 2 t m 2 ) cos ( ω 2 t m 3 ) cos ( ω 2 t m 4 ) sin ( ω 2 t m 1 ) sin ( ω 2 t m 2 ) sin ( ω 2 t m 3 ) sin ( ω 2 t m 4 ) = 0 , m 1 m 2 m 3 m 4 ; m 1 , m 2 , m 3 , m 4 = 1 , K 1 ¯ .
Then, in order to find τ i , it is required to use (15) together with K 4 equations from the set (23), so that all the moments of t i are contained there. A set that satisfies these conditions is (20), for example.
Remark 3.
Theorem 1 allows finding all of the required switching moments t i , i = 1 , K 1 ¯ that define vector C as the kernel of the corresponding linear mapping in (19).
This implies another necessary condition of the extremum, which is formulated in the form of the following lemma.
Lemma 2.
The optimal control u ( t ) with vector C obtained from (19) contains exactly K 1 zeros and, according to Definition 1, corresponds to the ( K 1 ) -switching control class.
In the next section, a theorem considering the continuity of the problem criterion from the parameter ε is formulated and proved.

5. Continuity of T 0 ( ε )

To find sufficient conditions for the continuity of a function T 0 ( ε ) , the implicit function theorem is applied. Let us write out the Jacobian of a system consisting of equations of motion (15) and additional conditions (16) obtained by Theorem 1 on the necessary extremum conditions.
The equations of the joint system are denoted by f 1 , f 2 , , f K , where f 1 , f 2 , f 3 , f 4 relate to the equations of motion (15), and f 5 , , f K are responsible for the non-degeneracy conditions (20). The Jacobian J K 1 corresponds to the ( K 1 ) -switching control class and consists of partial derivatives of K functions with K variables.
J K 1 = det f 1 τ 1 f 1 τ 2 f 1 τ K f 2 τ 1 f 2 τ 2 f 2 τ K f K τ 1 f K τ 2 f K τ K =
= det 2 ω 1 sin ω 1 i = 1 K τ i 2 ω 1 j = 1 2 1 j + 1 sin ω 1 i = j K τ i 2 ω 1 j = 1 K 1 j + 1 sin ω 1 i = j K τ i 2 ω 1 cos ω 1 i = 1 K τ i 2 ω 1 j = 1 2 1 j + 1 cos ω 1 i = j K τ i 2 ω 1 j = 1 K 1 j + 1 cos ω 1 i = j K τ i 2 ω 2 sin ω 2 i = 1 K τ i 2 ω 2 j = 1 2 1 j + 1 sin ω 2 i = j K τ i 2 ω 2 j = 1 K 1 j + 1 sin ω 2 i = j K τ i 2 ω 2 cos ω 2 i = 1 K τ i 2 ω 2 j = 1 2 1 j + 1 cos ω 2 i = j K τ i 2 ω 2 j = 1 K 1 j + 1 cos ω 2 i = j K τ i f 5 τ 1 f 5 τ 2 f 5 τ K f K τ 1 f K τ 2 f K τ K +
+ ω 1 sin ω 1 i = 1 K τ i ω 1 sin ω 1 i = 1 K τ i ω 1 sin ω 1 i = 1 K τ i ω 1 cos ω 1 i = 1 K τ i ω 1 cos ω 1 i = 1 K τ i ω 1 cos ω 1 i = 1 K τ i ω 2 sin ω 2 i = 1 K τ i ω 2 sin ω 2 i = 1 K τ i ω 2 sin ω 2 i = 1 K τ i ω 2 cos ω 2 i = 1 K τ i ω 2 cos ω 2 i = 1 K τ i ω 2 cos ω 2 i = 1 K τ i 0 0 0 0 0 0 .
Let us transform the above determinant, leading the Jacobian to a simpler form. The determinant of the system will not change if the previous column is subtracted from each column, except for the first.
J K 1 = det ω 1 sin ω 1 i = 1 K τ i 2 ω 1 sin ω 1 i = 2 K τ i 2 ω 1 ( 1 ) K + 1 sin ω 1 τ K ω 1 cos ω 1 i = 1 K τ i 2 ω 1 cos ω 1 i = 2 K τ i 2 ω 1 ( 1 ) K + 1 cos ω 1 τ K ω 2 sin ω 2 i = 1 K τ i 2 ω 2 sin ω 2 i = 2 K τ i 2 ω 2 ( 1 ) K + 1 sin ω 2 τ K ω 2 cos ω 2 i = 1 K τ i 2 ω 2 cos ω 2 i = 2 K τ i 2 ω 2 ( 1 ) K + 1 cos ω 2 τ K f 5 τ 1 f 5 τ 2 f 5 τ 1 f 5 τ K f 5 τ K 1 f K τ 1 f K τ 2 f K τ 1 f K τ K f K τ K 1 .
After multiplying the first and third rows of the resulting matrix by −1 and then replacing the sign for each even column, the transformed Jacobian has the following form
J K 1 = ( 1 ) [ K 2 ] 8 ω 1 2 ω 2 2 det sin ω 1 i = 1 K τ i sin ω 1 i = 2 K τ i sin ω 1 τ K cos ω 1 i = 1 K τ i cos ω 1 i = 2 K τ i cos ω 1 τ K sin ω 2 i = 1 K τ i sin ω 2 i = 2 K τ i sin ω 2 τ K cos ω 2 i = 1 K τ i cos ω 2 i = 2 K τ i cos ω 2 τ K 2 f 5 τ 1 f 5 τ 2 + f 5 τ 1 ( 1 ) K + 1 f 5 τ K f 5 τ K 1 2 f K τ 1 f K τ 2 + f K τ 1 ( 1 ) K + 1 f K τ K f K τ K 1 .
The Jacobian, according to Equation (25), for the three-switching control class is written as follows
J 3 = 8 ω 1 2 ω 2 2 det sin ω 1 ( τ 1 + τ 2 + τ 3 + τ 4 ) sin ω 1 ( τ 2 + τ 3 + τ 4 ) sin ω 1 ( τ 3 + τ 4 ) sin ω 1 τ 4 cos ω 1 ( τ 1 + τ 2 + τ 3 + τ 4 ) cos ω 1 ( τ 2 + τ 3 + τ 4 ) cos ω 1 ( τ 3 + τ 4 ) cos ω 1 τ 4 sin ω 2 ( τ 1 + τ 2 + τ 3 + τ 4 ) sin ω 2 ( τ 2 + τ 3 + τ 4 ) sin ω 2 ( τ 3 + τ 4 ) sin ω 2 τ 4 cos ω 2 ( τ 1 + τ 2 + τ 3 + τ 4 ) cos ω 2 ( τ 2 + τ 3 + τ 4 ) cos ω 2 ( τ 3 + τ 4 ) cos ω 2 τ 4 .
The expression (26) has an explicit form
J 3 8 ω 1 2 ω 2 2 = sin ( w 2 τ 3 ) sin ( w 1 τ 1 ) sin ( w 2 ( τ 2 + τ 3 ) ) sin ( w 1 ( τ 1 + τ 2 ) ) sin ( w 2 ( τ 1 + τ 2 ) ) sin ( w 1 ( τ 2 + τ 3 ) ) + sin ( w 2 τ 2 ) sin ( w 1 ( τ 1 + τ 2 + τ 3 ) ) + + sin ( w 2 ( τ 1 + τ 2 + τ 3 ) ) sin ( w 1 τ 2 ) + sin ( w 2 τ 1 ) sin ( w 1 τ 3 ) .
Expressions (21) and (26) can be compared, or rather their explicit forms (22) and (27), respectively. It turns out that with an index shift by one at τ i , the functional parts of these expressions are equal. Expression (26) includes four control intervals in the three-switching class and Expression (21) includes five intervals in the four-switching class, where the first control interval becomes zero at the boundary value ε , and that is how the transition into the three-switching class occurs.
This result is derived from several interrelated observations. The first four lines of the Jacobian (25), ordered by cos, sin, and ω 1 , ω 2 , which obviously do not change the sign of the determinants, can be written as
S 4 , K = cos ω 1 T 0 cos ω 1 ( T 0 t 1 ) cos ω 1 ( T 0 t K 1 ) sin ω 1 T 0 sin ω 1 ( T 0 t 1 ) sin ω 1 ( T 0 t K 1 ) cos ω 2 T 0 cos ω 2 ( T 0 t 1 ) cos ω 2 ( T 0 t K 1 ) sin ω 2 T 0 sin ω 2 ( T 0 t 1 ) sin ω 2 ( T 0 t K 1 ) = ( Γ 0 , Γ 1 , , Γ K 1 )
and then converted to the form
S 4 , K = cos ω 1 T 0 sin ω 1 T 0 sin ω 1 T 0 cos ω 1 T 0 1 cos ω 1 t 1 cos ω 1 t K 1 0 sin ω 1 t 1 sin ω 1 t K 1 cos ω 2 T 0 sin ω 2 T 0 sin ω 2 T 0 cos ω 2 T 0 1 cos ω 2 t 1 cos ω 2 t K 1 0 sin ω 2 t 1 sin ω 2 t K 1 .
The matrix S 4 , K includes the matrix
S K = 1 cos ω 1 t 1 cos ω 1 t K 1 0 sin ω 1 t 1 sin ω 1 t K 1 1 cos ω 2 t 1 cos ω 2 t K 1 0 sin ω 2 t 1 sin ω 2 t K 1 = ( Ω 0 , Ω 1 , , Ω K 1 ) ,
where the notation Ω 0 = ( 1 , 0 , 1 , 0 ) T is introduced. For the three-switching control class K = 4 , the equality is fulfilled
det S 4 , 4 = det S 4 .
The last equality follows from the next lemma.
Lemma 3.
If a square matrix A of size N × N is composed of M blocks of matrices A i of size M i × N , such that i = 1 M M i = N , and a square matrix B of size N × N is made up of M blocks of the product of matrices B i A i , where B i are square matrices of size M i × M i , then
det B = det A i = 1 M det B i .
Proof. 
Matrix B can be written as the product of matrices B = B ˜ A , where matrix A has the form
A = A 1 A 2 A M ,
and the auxiliary matrix B ˜ containing B i , i = 1 , M ¯ has the form
B ˜ = B 1 0 0 0 B 2 0 0 0 B M
Due to the multiplicativity of the determinant and the properties of the block–diagonal matrix, the determinant of matrix B has the form
det B = det ( B ˜ A ) = det B ˜ det A = det A i = 1 M det B i .
According to Lemma 3, the corresponding matrices for the three-switching control class are equal
B = S 4 , 4 , A = S 4 , B i = cos ω i T 0 sin ω i T 0 sin ω i T 0 cos ω i T 0 , i = 1 , 2 ¯ .
Moreover, det B i = 1 , i = 1 , 2 ¯ , so the equality (31) is established.
One can see that for the three-switching control class, the condition det S 4 = 0 is possible only if the necessary conditions of the extremum (23) and the equality t 1 = 0 for for the four-switching control class are fulfilled simultaneously.
The Jacobian J K 1 (25) can be written via the determinants of the square matrices of size 4 × 4 , composed of Ω i , i = 0 , , K 1 of the matrix S K .
J K 1 = 0 j 0 < j 1 < j 2 < j 3 K 1 α j 0 , j 1 , j 2 , j 3 det ( Γ j 0 , Γ j 1 , Γ j 2 , Γ j 3 ) = = 0 j 0 < j 1 < j 2 < j 3 K 1 α j 0 , j 1 , j 2 , j 3 det ( Ω j 0 , Ω j 1 , Ω j 2 , Ω j 3 ) .
The expression (35) is a decomposition of (25) by the last K 4 rows, while the first four rows are divided into C K 4 corresponding minors.
For further narration, it is necessary to write down the implicit function theorem [19] for the case of the following system of equations.
F 1 ( x 1 , , x m , y 1 , , y n ) = 0 , F n ( x 1 , , x m , y 1 , , y n ) = 0 ,
which is solved with respect to the variables y 1 , , y n . The following system of functional relations (locally equivalent to system (36)) is studied.
y 1 = g 1 ( x 1 , , x m ) , y n = g n ( x 1 , , x m ) .
The following designations are introduced:
x = ( x 1 , , x m ) , y = ( y 1 , , y n ) .
System (36) is denoted with F ( x , y ) = 0 and the mapping (37) – with y = g ( x ) .
g ( x ) = g 1 x 1 g 1 x m g n x 1 g n x m ( x ) , F x ( x , y ) = F 1 x 1 F 1 x m F n x 1 F n x m ( x , y ) , F y ( x , y ) = F 1 y 1 F 1 y n F n y 1 F n y n ( x , y ) ,
A square matrix F y ( x , y ) is reversible only when its determinant is non-zero.
Theorem 2
(The implicit function theorem as in [19]). If a mapping F : U R n , defined in the neighborhood U of a point ( x 0 , y 0 ) R m + n is such that the following requirements are met
  • F C ( p ) ( U ; R n ) , p 1 ,
  • F ( x 0 , y 0 ) = 0 ,
  • F y ( x 0 , y 0 ) − reversible matrix,
then there exists a ( m + n ) -dimensional interval I = I x m × I y n U where
I x m = { x R m | | x x 0 | < α } , I y n = { y R n | | y y 0 | < β } ,
and a mapping g C ( p ) ( I x m ; I y n ) , such that for any point ( x , y ) I x m × I y n the following condition holds
F ( x , y ) = 0 y = g ( x ) ,
and
g ( x ) = [ F y ( x , g ( x ) ) ] 1 [ F x ( x , g ( x ) ) ]
Now, it is possible to prove the continuity, not only of τ 1 ( ε ) , , τ K ( ε ) but also of the criterion T 0 ( ε ) by the implicit function theorem applied to our problem.
Lemma 4.
For continuity of the functions T 0 ( ε ) , τ 1 ( ε ) , , τ K ( ε ) for K 1 -switching control class, the following condition is required
J K 1 0 .
Proof. 
The implicit function theorem 2 is applied for the considered problem. Let n = K , then
F 1 ( x 1 , , x m , y 1 , , y n ) = f 1 ( ε , τ 1 , , τ K ) = 0 , F n ( x 1 , , x m , y 1 , , y n ) = f K ( ε , τ 1 , , τ K ) = 0 ,
x = ε , y = ( τ 1 , , τ K ) ,
where, as before, f 1 , f 2 , , f K refer to the equations of the systems (15) and (20). For the sake of the reader’s comprehension, the arguments of the functions are omitted.
F y = f 1 τ 1 f 1 τ K f K τ 1 f K τ K , g = g 1 ε g K ε , F x = f 1 ε f K ε = ( 1 ) k + 1 q 1 T 0 ω 1 2 ε 2 ( 1 ) k p 1 T 0 ω 1 ε 2 0 0 .
The first condition of Theorem 2 is satisfied if the parameter ε is non-zero. For the solution of the optimal control problem, taking into account the term of Theorem 1, the second condition of Theorem 2 is fulfilled. The reversibility of F y requires the condition det ( F y ) 0 , which is equivalent to J K 1 0 . All conditions of the implicit function theorem are satisfied, so there is a ( K + 1 ) -dimensional interval I = I ε 1 × I τ 1 , , τ K K and continuous mapping f on this interval.
τ 1 = g 1 ( ε ) , τ K = g K ( ε ) .
Consequently, the durations of the control intervals τ 1 , , τ K are continuous functions of ε . Since T 0 = τ 1 + + τ K , then T 0 ( ε ) is a continuous function of ε too. □
Lemma 5
(Sufficient conditions of continuity of T 0 ( ε ) ). T 0 ( ε ) is continuous if the non-zero conditions of (23) are met in the ( K 1 ) -switching class.
Proof. 
We investigate the Jacobian J K 1 , which we will rewrite using the formula (35) in a form
J K 1 = 0 j 0 < j 1 < j 2 < j 3 K 1 α j 0 , j 1 , j 2 , j 3 det ( Ω j 0 , Ω j 1 , Ω j 2 , Ω j 3 ) = = j 0 , j 1 , j 2 , j 3 I 1 α j 0 , j 1 , j 2 , j 3 det ( Ω j 0 , Ω j 1 , Ω j 2 , Ω j 3 ) + j 0 , j 1 , j 2 , j 3 I 2 α j 0 , j 1 , j 2 , j 3 det ( Ω j 0 , Ω j 1 , Ω j 2 , Ω j 3 ) ,
where the total sum is divided into two parts. The set of indexes I 1 always contains a null time element that corresponds to Ω 0 . I 2 corresponds to all other sets without a null time element. The fulfillment of the conditions (23) leads to zeroing of the second sum in (41), due to the fact that each summand is zero. Then, it holds
J K 1 = j 0 , j 1 , j 2 , j 3 I 1 α j 0 , j 1 , j 2 , j 3 det ( Ω j 0 , Ω j 1 , Ω j 2 , Ω j 3 ) .
Due to the linear dependence of any set consisting of four different Ω j , j = 1 , K 1 ¯ by the condition (23), any vector Ω j except Ω 0 can be written as a linear combination of Ω 1 , Ω 2 , Ω 3 . Then (42) can be expressed as follows
J K 1 = j β j det ( Ω 0 , Ω 1 , Ω 2 , Ω 3 ) ,
where β j combines both the coefficients α j 0 , j 1 , j 2 , j 3 and the coefficients of expansion Ω j , j = 1 , K 1 ¯ by Ω 1 , Ω 2 , Ω 3 .
Further, two cases are possible, namely J K 1 0 , when T 0 ( ε ) is continuous by Lemma 4, or J K 1 = 0 . In this case, if det ( Ω 0 , Ω 1 , Ω 2 , Ω 3 ) = 0 , then Ω 0 , Ω 1 , Ω 2 , Ω 3 are linearly dependent. Moreover, if det ( Ω 0 , Ω 1 , Ω 2 , Ω 3 ) 0 , then due to J K 1 = 0 and the presence in the first column of only the first and third elements other than zero, it follows that the first and third rows match. Which again entails a linear dependence of Ω 0 , Ω 1 , Ω 2 , Ω 3 . This means that the necessary conditions of the extremum include Ω 0 corresponding to t 1 = 0 and increasing the number of control switchings by one, which changes the total number of control switchings from K 1 to K. Thus, the conditions of the lemma for the number of control switchings are violated, which means that J K 1 0 . By Lemma 4, the function T 0 ( ε ) in the ( K 1 ) -switching control class is continuous. □
Corollary 1.
In a given control switching class, T 0 ( ε ) is a right-continuous function.
Remark 4.
It can be concluded from entry (25) that due to the absence of τ 1 in the explicit form f 5 , , f K , the sum of I 2 in (41) always equals zero.
Corollary 2.
Fulfillment of C K 1 4 Equation (23) in the K 1 -switching control class under the condition t 1 = 0 leads to J K 2 = 0 .
Lemma 6.
T 0 ( ε ) is a continuous function for ε ( 0 , + ] .
Proof. 
Let two adjacent control classes consist of P and Q switchings, respectively. Indeed, due to Lemma 5 and Corollary 1 in a given class of control switchings, T 0 ( ε ) is continuous. Let us denote the value of the criterion in each of the classes T P P ( ε ) = T 0 ( ε ) with ε ( ε 1 , ε P ] and T Q Q ( ε ) = T 0 ( ε ) at ε ( ε P , ε Q ] , where the upper index indicates the number of control switchings, and the lower index indicates the optimal value of control switchings. Next, for T Q P ( ε ) , let us denote the value of the criterion in the P-switching control class with their optimal number—Q for ε ( ε P , ε Q ] , and vice versa T P Q ( ε ) , the value of the criterion in the Q-switching control class with their optimal number is P and ε ( ε 1 , ε P ] . Then obviously the inequalities are fulfilled, T P P ( ε ) T P Q ( ε ) and T Q Q ( ε ) T Q P ( ε ) . Again, let us use Lemma 5 and Corollary 1. So T P Q ( ε ) and T Q P ( ε ) are continuous on ε ( ε 1 , ε Q ] = ( ε 1 , ε P ] ( ε P , ε Q ] . Then T P Q ( ε P ) = T Q P ( ε P ) . Thus, taking into account the above inequalities, it follows that T P P ( ε P ) = lim ε ε P + 0 T Q Q ( ε ) . The statement of the lemma follows from the union of all intervals of control switchings. □
Continuity of T 0 ( ε ) means that T 0 is continuous in the phase plane of the first oscillator on terminal positions ( p 1 T 0 , q 1 T 0 ) in each direction q 1 T 0 = λ p 1 T 0 . Then, applying the equations of characteristics (7) according to [20], where T 0 p 1 , T 0 q 1 exist, we obtain
δ T 0 = T 0 p 1 p 1 = p 1 T 0 δ p 1 + T 0 q 1 q 1 = q 1 T 0 δ q 1 = ξ 1 ( T 0 ) δ p 1 + η 1 ( T 0 ) δ q 1 .
The adjoint variables ξ 1 ( t ) , η 1 ( t ) are bounded, including at t = T 0 , so the variation of the criterion δ T 0 will be quite small under small variations of the terminal positions δ p 1 , δ q 1 .

6. Modeling

To illustrate the aforesaid theoretical results, numerical modeling for various configurations of the (1)–(5) problem were presented. All numerical experiments were conducted via Python scripts, which were developed specifically for this problem. In particular, the SciPy library was used to find solutions to the systems of nonlinear equations, and Matplotlib was used to visualize the obtained results.
Firstly, in Section 6.1, the motion of the system to state ( 0 , 1 , 0 , 0 ) (with optimal control u ( t ) (13) for control classes with different switchings) was described. In these examples, the dependencies of the optimal number of control switchings, the criterion of the problem T 0 , and others on parameter ε , were investigated. In Section 6.2, a range of terminal states of the first oscillator was considered, but now with a fixed value of ε . In this case, all the values were considered to be functions of the parameters ( q 1 T 0 , p 1 T 0 ) .

6.1. Fixed Terminal Position of the First Oscillator

The problem (1)–(5) is considered with the following parameters chosen without loss of generality:
ω 1 = 1 , ω 2 = 1.4 , q 1 T 0 = 0 , p 1 T 0 = 1 , q 2 T 0 = 0 , p 2 T 0 = 0 .
Control classes with three, four, and five switchings ( K 1 = 3 , 4 , 5 ) are considered. For a range of values ε in each control class, the solutions that satisfied Theorem 1 and minimized the criterion (5) were found.
An important step is checking that the resulting system solutions form, according to Remark 3, the control law with the number of switchings that are initially set. In other words, Lemma 2 must be fulfilled. For example, an attempt to find a solution for a three-switching control class for ε < 0.6 results in switching function S F ( t ) containing more than three switchings that do not satisfy Lemma 2. Figure 2 shows such an example for the parameter ε = 0.599 . The points obtained by Theorem 1 are highlighted in red. An extra zero of function S F ( t ) is marked with a red cross. The constant vector C = 0.57 , 0.57 , 0.26 , 0.52 , calculated according to Remark 3, defines the switching function S F ( t ) , which is pictured with a blue line.
Lemma 2 is not satisfied in this case.
However, Lemma 2 will be fulfilled for ε = 0.76 , as illustrated in Figure 3. Moreover, the optimal control law u ( t ) defined by S F ( t ) is shown.
Vector C in this case: C = 0.14 , 0.72 , 0.48 , 0.46 . The phase portraits of the first and second oscillators satisfy the given boundary conditions (44) and are presented in Figure 4.
The dependence of the criterion T 0 on the constraint ε is shown in Figure 5. The blue, red, and black curves show solutions for three-, four-, and five-switching control classes, which are defined by the boundaries of each control class. Figure 5 shows that T 0 ( ε ) is a continuous function according to Lemma 6.
Figure 6 illustrates the dependencies of the control interval durations from ε . The implicit function theorem shows that the durations of control intervals are continuous on ε inside each class, which can be observed in the graph. In area I, which corresponds to the five-switching class, the duration of the first control interval becomes zero when the boundary of the selected class ( I II ) is reached. Region II corresponds to the four-switching class, where the second interval of I becomes the first interval, and it also becomes zero when it reaches the boundary ( II III ) . The third interval from region I becomes the second interval in II and the first interval in the three-switching region (III). If one of the internal intervals of the four-switching class is nullified, then the transition to the two-switching class takes place.
Figure 7 shows the dependencies of the Jacobian for different switching control classes on ε .
As one can see from the graph, when the class boundaries are reached, the Jacobians become equal to zero. Moreover, the Jacobian J 4 of the four-switching class is zero when passing to the two-switching class at one point. Lemma 1 is valid for the two-switching control class. Two-, three-, and four-switching control classes for different finite states of the first oscillator ( q 1 T 0 , p 1 T 0 ) and ε = 0.4 are presented next. As before, it was checked that the obtained solutions for each final state form, according to Remark 3, the control law with the appropriate number of switchings.

6.2. The Phase Plane of the First Oscillator for Different Control Classes

Different control classes are shown in Figure 8. The red and green areas correspond to the three- and four-switching control classes with initial controls ε , ε , respectively. The reflections of these regions relative to the origins lead to blue and orange regions for the three- and four-switching control classes on the initial intervals ε , ε , respectively. The three-switching control classes are separated from each other by parametric curves obtained by Lemma 1. This is due to the fact that nullifying the lateral control interval in the three-switching control class also leads to the two-switching control class (orange and red curves). Zeroing the internal control interval in the ’four-switchings’ results in curves (blue and green) corresponding to the two-switching control classes at different initial controls. The yellow and gray boundary curves were obtained by moving from the four-switching class to the three-switching class by zeroing out the last control interval, with the control at the initial interval coinciding with the initial controls in the regions on either side of the boundary. The white and purple curves correspond to the zeroing of the first control interval in the four-switching class, which means the transition to the three-switching. The simultaneous nullifying of the side interval of the control constancy in the four-switching class corresponds to the blue points where the curves obtained under the conditions τ 1 = 0 or τ 2 = 0 intersect.
Figure 9 shows the dependence T 0 ( p 1 T 0 , q 1 T 0 ) of the criterion on the phase coordinates of the first oscillator. Three- and four-switching control classes are highlighted with red and green points. The colors of all curves correspond to the description in Figure 8.

7. Conclusions

Necessary and sufficient extremum conditions are proposed for the time-optimal control problem of a system consisting of two non-synchronous oscillators, and controlled by limited scalar impact. In addition, lemmas and theorems on the continuity of the criterion by the parameter are proposed. These results can be further extended to systems of any number of oscillators. Thus, a unified mathematical apparatus for describing the set of reachability of the first oscillator for systems of many oscillators will be developed in future works. The construction of control switching curves of the first oscillator in the phase plane will allow us to develop an approach to solving the problem of the synthesis of time-optimal control for multiple oscillators.

Author Contributions

Conceptualization, A.G. and P.L.; methodology, A.G.; software, L.B. and P.L.; validation, L.B. and A.G.; formal analysis, L.B. and A.G.; investigation, L.B., A.G. and P.L.; writing—original draft preparation, L.B. and A.G.; writing—review and editing, L.B., A.G. and P.L.; visualization, L.B.; supervision, A.G.; project administration, A.G.; funding acquisition, A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a grant to support youth scientific schools of ICS RAS “Methods for trajectories optimization of controlled objects”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof 
(Proof of Lemma 1.).
From dynamic Equation (7) for any terminal position ( p 1 T 0 , q 1 T 0 ) of the first oscillator and K = 3 of the next equation system is valid
cos ( ω 1 ( τ 1 + τ 2 + τ 3 ) ) 2 cos ( ω 1 ( τ 2 + τ 3 ) ) + 2 cos ( ω 1 ( τ 3 ) ) = 1 ( 1 ) k + 1 q 1 T 0 ω 1 2 ε , sin ( ω 1 ( τ 1 + τ 2 + τ 3 ) ) 2 sin ( ω 1 ( τ 2 + τ 3 ) ) + 2 sin ( ω 1 ( τ 3 ) ) = ( 1 ) k p 1 T 0 ω 1 ε , cos ( ω 2 ( τ 1 + τ 2 + τ 3 ) ) 2 cos ( ω 2 ( τ 2 + τ 3 ) ) + 2 cos ( ω 2 ( τ 3 ) ) = 1 , sin ( ω 2 ( τ 1 + τ 2 + τ 3 ) ) 2 sin ( ω 2 ( τ 2 + τ 3 ) ) + 2 sin ( ω 2 ( τ 3 ) ) = 0 .
Let us consider separately the last two equations of the system (A1) in the follow- ing form
cos ( ω 2 ( τ 1 + τ 2 + τ 3 ) ) 2 cos ( ω 2 ( τ 2 + τ 3 ) ) = 1 2 cos ( ω 2 ( τ 3 ) ) , sin ( ω 2 ( τ 1 + τ 2 + τ 3 ) ) 2 sin ( ω 2 ( τ 2 + τ 3 ) ) = 2 sin ( ω 2 ( τ 3 ) ) .
Equations of the system (A2) are squared and summed
1 4 cos ( ω 2 ( τ 3 ) ) + 4 = 1 + 4 4 cos ( ω 2 ( τ 1 ) ) ,
cos ( ω 2 ( τ 3 ) ) = cos ( ω 2 ( τ 1 ) ) .
The first statement of Lemma 1 follows from (A4).
In the next step, the last two equations of the system (A1) are squared, summed, and then transformed. Moreover, it is noticed that τ 2 0 and τ 2 < 2 π ω 2 according to lemma condition.
2 + cos ( ω 2 ( τ 1 + τ 2 ) ) 2 cos ( ω 2 τ 2 ) = cos ( ω 2 τ 1 ) ,
2 + cos ( ω 2 τ 1 ) cos ( ω 2 τ 2 ) sin ( ω 2 τ 1 ) sin ( ω 2 τ 2 ) 2 cos ( ω 2 τ 2 ) = cos ( ω 2 τ 1 ) ,
cos ( ω 2 τ 2 ) ( cos ( ω 2 τ 1 ) 2 ) + 2 cos ( ω 2 τ 1 ) = sin ( ω 2 τ 1 ) sin ( ω 2 τ 2 ) ,
( cos ( ω 2 τ 1 ) 2 ) ( cos ( ω 2 τ 2 ) 1 ) = sin ( ω 2 τ 1 ) sin ( ω 2 τ 2 ) ,
( cos ( ω 2 τ 1 ) 2 ) 2 ( 1 cos ( ω 2 τ 2 ) ) 2 = ( 1 cos 2 ( ω 2 τ 1 ) ) ( 1 cos 2 ( ω 2 τ 2 ) ) ,
( cos ( ω 2 τ 1 ) 2 ) 2 ( 1 cos ( ω 2 τ 2 ) ) 2 = ( 1 cos 2 ( ω 2 τ 1 ) ) ( 1 cos ( ω 2 τ 2 ) ) ( 1 + cos ( ω 2 τ 2 ) ) ,
( cos ( ω 2 τ 1 ) 2 ) 2 ( 1 cos ( ω 2 τ 2 ) ) = ( 1 cos 2 ( ω 2 τ 1 ) ) ( 1 + cos ( ω 2 τ 2 ) ) .
From this follows the required equality. Lemma 1 is proved. □

References

  1. Mirzaei, A.; Darabi, H. Mutual pulling between two oscillators. IEEE J. Solid-State Circuits 2014, 49, 360–372. [Google Scholar] [CrossRef]
  2. Salobutina, E. Regimes of more and more frequent switchings in the optimal control problem of oscillations of n oscillators. J. Math. Sci. 2008, 151, 3603–3610. [Google Scholar] [CrossRef]
  3. Kayumov, O. Time-optimal movement of platform with oscillators. Mech. Solids 2021, 56, 1622–1637. [Google Scholar] [CrossRef]
  4. Ansel, Q.; Chepelianskii, A.D.; Lages, J. Enhancing quantum exchanges between two oscillators. arXiv 2022, arXiv:2207.11156v1. [Google Scholar]
  5. Zhao, Y.; Chen, G.H. Two oscillators in a dissipative bath. Phys. A Stat. Mech. Appl. 2003, 317, 13–40. [Google Scholar] [CrossRef]
  6. Li, J.; Chen, Y.; Zhao, S. Periodic solutions and stability analysis for two-coupled-oscillator structure in optics of chiral molecules. Mathematics 2022, 10, 1908. [Google Scholar] [CrossRef]
  7. Firippi, E.; Chaves, M. Period-control in a coupled system of two genetic oscillators for synthetic biology. IFAC-PapersOnLine 2019, 52, 70–75. [Google Scholar] [CrossRef]
  8. Galyaev, A.A. Scalar control of a group of free-running oscillators. Autom. Remote Control 2016, 77, 1511–1523. [Google Scholar] [CrossRef]
  9. Galyaev, A.A. Impact of a system of material points against an absolutely rigid obstacle: A model for its impulsive action. Autom. Remote Control 2006, 67, 856–867. [Google Scholar] [CrossRef]
  10. Andresen, B.; Salamon, P.; Hoffmann, K.H.; Tsirlin, A.M. Optimal processes for controllable oscillators. Autom. Remote Control 2018, 79, 2103–2113. [Google Scholar] [CrossRef]
  11. Fedorov, A.; Ovseevich, A. Asymptotic control theory for a system of linear oscillators. Mosc. Math. J. 2016, 16, 561–598. [Google Scholar] [CrossRef]
  12. Pesch, H.J.; Plail, M. The Maximum Principle of optimal control: A history of ingenious ideas and missed opportunities. Control Cybern. 2009, 38, 973–995. [Google Scholar]
  13. Boltyansky, V.G. Mathematical Methods of Optimal Control, 2nd ed.; Nauka: Moscow, Russia, 1969. [Google Scholar]
  14. Chernousko, F.L.; Akulenko, L.D.; Sokolov, B.N. Oscillation Control; Nauka: Moscow, Russia, 1980. [Google Scholar]
  15. Berlin, L.M.; Galyaev, A.A. Extremum conditions for constrained scalar control of two nonsynchronous oscillators in the time-optimal control problem. Dokl. Math. 2022, 106, 286–290. [Google Scholar] [CrossRef]
  16. Sachkov, Y.L. Introduction to Geometric Control. arXiv 2021, arXiv:1903.00211v2. [Google Scholar]
  17. La Salle, J.P. The Time Optimal Control Problem. In Contributions to the Theory of Nonlinear Oscillations; Princeton University Press: Princeton, NJ, USA, 1960; Volume 5, pp. 1–24. [Google Scholar]
  18. Berlin, L.M.; Galyaev, A.A.; Lysenko, P.V. About time-optimal control problem for system of two non-synchronous oscillators. In Proceedings of the 16th International Conference on Stability and Oscillations of Nonlinear Control Systems (Pyatnitskiy’s Conference), Moscow, Russian, 1–3 June 2022. [Google Scholar] [CrossRef]
  19. Zorich, V.A. Mathematical Analysis I, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2015; p. 616. [Google Scholar]
  20. Isaacs, R. Differential Games: A Mathematical Theory with Applications to Warfare and Pursuit, Control and Optimization; John Wiley and Sons: New York, NY, USA, 1965. [Google Scholar]
Figure 1. Control intervals for optimal control u ( t ) .
Figure 1. Control intervals for optimal control u ( t ) .
Mathematics 10 03552 g001
Figure 2. Switching function S F ( t ) for ε = 0.599 .
Figure 2. Switching function S F ( t ) for ε = 0.599 .
Mathematics 10 03552 g002
Figure 3. Switching function S F ( t ) and optimal control u ( t ) for ε = 0.76 .
Figure 3. Switching function S F ( t ) and optimal control u ( t ) for ε = 0.76 .
Mathematics 10 03552 g003
Figure 4. The phase portraits of the first and second oscillators.
Figure 4. The phase portraits of the first and second oscillators.
Mathematics 10 03552 g004
Figure 5. Dependence T 0 ( ε ) .
Figure 5. Dependence T 0 ( ε ) .
Mathematics 10 03552 g005
Figure 6. Dependencies τ i ( ε ) .
Figure 6. Dependencies τ i ( ε ) .
Mathematics 10 03552 g006
Figure 7. Dependence J ( ε ) .
Figure 7. Dependence J ( ε ) .
Mathematics 10 03552 g007
Figure 8. Terminal points of the phase plane of the first oscillator.
Figure 8. Terminal points of the phase plane of the first oscillator.
Mathematics 10 03552 g008
Figure 9. The dependence of the criterion T 0 on the phase coordinates of the first oscillator ( p 1 T 0 , q 1 T 0 ) .
Figure 9. The dependence of the criterion T 0 on the phase coordinates of the first oscillator ( p 1 T 0 , q 1 T 0 ) .
Mathematics 10 03552 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Berlin, L.; Galyaev, A.; Lysenko, P. Time-Optimal Control Problem of Two Non-Synchronous Oscillators. Mathematics 2022, 10, 3552. https://doi.org/10.3390/math10193552

AMA Style

Berlin L, Galyaev A, Lysenko P. Time-Optimal Control Problem of Two Non-Synchronous Oscillators. Mathematics. 2022; 10(19):3552. https://doi.org/10.3390/math10193552

Chicago/Turabian Style

Berlin, Leonid, Andrey Galyaev, and Pavel Lysenko. 2022. "Time-Optimal Control Problem of Two Non-Synchronous Oscillators" Mathematics 10, no. 19: 3552. https://doi.org/10.3390/math10193552

APA Style

Berlin, L., Galyaev, A., & Lysenko, P. (2022). Time-Optimal Control Problem of Two Non-Synchronous Oscillators. Mathematics, 10(19), 3552. https://doi.org/10.3390/math10193552

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop