Next Article in Journal
Generalizations of Hermite–Hadamard Type Integral Inequalities for Convex Functions
Next Article in Special Issue
Robust Stability of Time-Varying Markov Jump Linear Systems with Respect to a Class of Structured, Stochastic, Nonlinear Parametric Uncertainties
Previous Article in Journal
Mathematical Analysis of a Fractional COVID-19 Model Applied to Wuhan, Spain and Portugal
Previous Article in Special Issue
Nash Equilibrium Sequence in a Singular Two-Person Linear-Quadratic Differential Game
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two Inverse Problems Solution by Feedback Tracking Control

by
Vladimir Turetsky
Department of Mathematics, Ort Braude College of Engineering, 51 Snunit Str., P.O. Box 78, Karmiel 2161002, Israel
Axioms 2021, 10(3), 137; https://doi.org/10.3390/axioms10030137
Submission received: 10 May 2021 / Revised: 8 June 2021 / Accepted: 21 June 2021 / Published: 28 June 2021
(This article belongs to the Special Issue Advances in Analysis and Control of Systems with Uncertainties)

Abstract

:
Two inverse ill-posed problems are considered. The first problem is an input restoration of a linear system. The second one is a restoration of time-dependent coefficients of a linear ordinary differential equation. Both problems are reformulated as auxiliary optimal control problems with regularizing cost functional. For the coefficients restoration problem, two control models are proposed. In the first model, the control coefficients are approximated by the output and the estimates of its derivatives. This model yields an approximating linear-quadratic optimal control problem having a known explicit solution. The derivatives are also obtained as auxiliary linear-quadratic tracking controls. The second control model is accurate and leads to a bilinear-quadratic optimal control problem. The latter is tackled in two ways: by an iterative procedure and by a feedback linearization. Simulation results show that a bilinear model provides more accurate coefficients estimates.

1. Introduction

The problems of restoration of a system input and coefficients arise in various applications such as metrology [1,2,3,4], image processing [5,6], spectroscopy [7], geophysics [8,9], process control [10], engineering [11,12], medicine [13,14], and many others. The identification of the virus spread models, i.e., the restoration of the coefficients of the corresponding differential equations, became crucial in the period of coronavirus pandemic (see, e.g., [15,16]) for constructing reliable forecasts of the pandemic evolution. These problems belong to a wide class of ill-posed (inverse) problems [17,18,19]. The general framework for the solution of these problems is given by the Tikhonov regularization [17], but its specific implementation varies widely depending on the system description, assumptions on the noise etc.
In this paper, we deal with the input and time-varying coefficients restoration problems for a system, modeled by a linear differential equation, based on noisy output measurements. Mathematically, they are formulated as a first-kind operator equations with compact linear operators mapping an input (coefficients) to the output. In practical formulations, it is assumed that both output and input (coefficient) signals belong to Hilbert spaces. Even if the inverse operator exists, it is unbounded [20], which makes an exact solution unstable w.r.t. output measurement errors. In the Tikhonov regularization method, the solution of the original equation is approximated by the solution of a minimization problem for a so-called smoothing functional, constructed as a sum of a squared discrepancy and some stabilizing functional (stabilizer) with a small coefficient (regularization parameter). The main result of Tikhonov [17] is that by a proper choice of a stabilizer and a regularization parameter, the approximate solution tends to the exact one while the measurement noise tends to zero.
Ill-posed problems in general and the problems of input and coefficient restoration, in particular, are tackled in the literature in numerous forms and by numerous approaches which mainly represent specific adaptations of Tikhonov’s regularization: for linear and non-linear systems [21], in the time [22] and in the frequency [23] domains, in statistical [3,24,25] and deterministic [19] formulations, etc. However, to the best knowledge of the author, a challenging problem of restoration of time-dependent coefficients of a high-order linear differential equation was not formulated previously.
In this paper, we adopt an approach developed in the Ekaterinburg (Sverdlovsk) School on Control Theory [22,26,27,28], where an ill-posed problem is reformulated as an optimal control problem with a Tikhonov-type tracking cost functional. Thus, an inverse problem is solved as a direct control problem, and the control realization approximates the input or the coefficients to be restored. Such a solution enjoys both the regularizing properties of the Tikhonov’s functional minimizer and the stability properties of a feedback control.
In this paper, we mainly concentrate on the problem of restoration of time-dependent coefficients of a linear differential equation. However, the input restoration problem is also formulated, and its solution, based on the auxiliary linear-quadratic tracking problem, is presented. It is important both for methodological purposes (the two problems are solved in similar lines), and from the practical point of view (the stable differentiation, which is the particular case of the input restoration, is implemented for solving the coefficient restoration problem). The novel approach, proposed in this paper, implies reformulating the problem of restoration of time-dependent coefficients of a linear differential equation of order n as an input restoration problem for an auxiliary bilinear system of order n 2 . Then, for a bilinear system, an optimal control problem with a Tikhonov-type cost function is posed, where unknown coefficients play the role of control. For a first-order differential equation, such approach was presented in [29].
Thus, an optimal control version of the coefficient restoration problem yields a bilinear-quadratic formulation. This control problem is treated both by a linear-quadratic approximation and in an exact formulation. The approximate approach needs the output derivatives which are approximated by optimal controls in another auxiliary linear-quadratic problem as it was carried out, e.g., in [30]. The exact, bilinear-quadratic, problem is solved by two approaches. The first one comprises an iterative solution of the corresponding Bellman equation based on [31]. In the second one, a feedback linearization [32] is implemented similarly to [33].

2. Input Restoration Problem

2.1. Ill-Posed Problem Formulation

Consider a linear time-invariant system given by a strictly proper transfer function
G ( s ) = b 1 s n 1 + b 2 s n 2 . . . + b n 1 s + b n s n + a 1 s n 1 + . . . + a n 1 s + a n .
The problem is a stable (w.r.t. measurement noise) restoration of an unknown scalar input u ( t ) based on noisy measurements
y ˜ ( t ) = y ( t ) + η ( t )
of the scalar output y ( t ) where
| η ( t ) | δ , t 0 ,
for some δ > 0 . Mathematically, this means to solve the operator (integral) equation
G u ( · ) = y ˜ ( · ) ,
with the compact operator G : L 2 [ 0 , t f ] L 2 [ 0 , t f ] , t f is a given time moment. An operator G can be, for example, an integral operator of optical instrument in the image restoration problem [5], or a convolution operator of a measuring device in the metrological signal restoration problem [1].
As it is acceptable in the literature on the ill-posed problems (see, e.g., [17] and references therein), it is assumed that for η ( t ) 0 , there exists an exact solution
u 0 ( · ) = G 1 y ( · ) .
It is well known [17] that the Equation (4) represents an ill-posed problem, i.e., the exact solution in the presence of noise u δ ( · ) = G 1 y ˜ ( · ) is unstable w.r.t. data errors. The problem is to construct an approximate solution u δ ( · ) such that
lim δ 0 | | u δ ( · ) u 0 ( · ) | | = 0 ,
where
| | u ( · ) | | = 0 t f u 2 ( t ) d t 1 / 2 .
Remark 1.
In Tikhonov’s regularization method [17], the approximate (regularized) stable solution of the Equation (4) is given by a minimizer u α ( · ) of a smoothing functional which can be chosen as
M α ( u ( · ) ) = | | G u ( · ) y ˜ ( · ) | | 2 + α | | u ( · ) | | 2 ,
where α > 0 is a regularization parameter. The minimizer u α ( · ) satisfies the integral Euler-Lagrange equation
G * G u ( · ) + α u ( · ) = G * y ˜ ( · ) ,
where G * is an operator adjoint to G . Explicitly,
u α ( · ) = G * G + α I 1 y ˜ ( · ) ,
where I is an identity operator. Moreover, [17] α = α ( δ ) exists, such that u δ ( · ) = u α ( δ ) ( · ) satisfies (6).

2.2. Optimal Control Problem Reformulation

The system determined by the transfer function (1) can be represented as a differential equation in an observable canonical form ([34], Ch. 2):
x ˙ = A x + b u ,
where
A = a 1 1 0 0 a 2 0 1 0 a n 1 0 0 1 a n 0 0 0 , b = b 1 b 2 b n , x 1 = y .
For the system (11), consider the linear-quadratic tracking problem with the cost functional
J α = 0 t f D x ( t ) y ˜ ( t ) 2 d t + α 0 t f u 2 ( t ) d t ,
where
D = [ 1 , 0 , 0 , . . . , 0 ] ,
α > 0 is a penalty coefficient for the control effort expenditure. The control problem is minimized (13) by a feedback control u ( t , x ) . Since G u ( · ) = y ( · ) = D x ( · ) = x 1 ( · ) , the cost (13) coincides with the Tikhonov’s functional (8). Thus, the problem of minimization of a Tikhonov’s smoothing functional can be solved by a time realization of the optimal strategy in the linear-quadratic tracking problem, (11), (13).
The problem (11), (13) is a particular case of a linear-quadratic regulator problem, and its solution is well-known [35]—namely, the optimal feedback strategy is
u α * ( t , x ) = 1 2 α b T Φ T ( t f , t ) 2 R α ( t ) Φ ( t f , t ) x + r α ( t ) ,
where Φ ( t , τ ) is the transition matrix of x ˙ = A x . The matrix function R α ( t ) satisfies the Riccati differential equation
R ˙ = R Q α ( t ) R S ( t ) , R ( t f ) = 0 ,
where
Q α ( t ) = 1 α Φ ( t f , t ) b b T Φ T ( t f , t )
S ( t ) = Φ T ( t , t f ) D T D Φ ( t , t f ) .
The vector function r α ( t ) satisfies the linear differential equation
r ˙ = R α ( t ) Q α ( t ) r + 2 Φ T ( t , t f ) D T y ˜ ( t ) , r ( t f ) = 0 .
Remark 2.
Both functions u α ( t ) = u α * ( t , x ( t ) ) , based on the feedback strategy (15), and u α ( · ) given by (10), solve the input restoration problem. However, a feedback solution is simpler from the computational point of view and it is more stable numerically.

2.3. Stable Differentiation Problem

Let the transfer function be
G ( s ) = 1 s ,
i.e., the system is the integrator. The integral Equation (4) becomes
0 t u ( τ ) d τ = y ˜ ( t ) y ˜ ( 0 ) , t [ 0 , t f ] ,
and the exact solution without measurement noise η ( t ) is u 0 ( t ) = y ˙ ( t ) . Thus, in this case, the signal restoration problem becomes the problem of stable differentiation of a noisy signal. The system can also be described by the differential equation
x ˙ = u , x ( 0 ) = x 0 .
The cost functional (13) in the auxiliary tracking problem becomes
J α = 0 t f ( x ( t ) y ˜ ( t ) ) 2 d t + α 0 t f u 2 ( t ) d t .
In the linear-quadratic problem (22)–(23), A = 0 , Φ ( t , τ ) 1 , b = 1 , D = 1 , and the optimal control (15) is given by
u α * ( t , x ) = 1 2 α 2 R α ( t ) x + r α ( t ) .
The differential Equation (16) for R α ( t ) becomes
R ˙ = 1 α R 2 1 , R ( t f ) = 0 .
By simple algebra, the Equation (25) admits the explicit solution
R α ( t ) = 1 μ tanh μ ( t f t ) ,
where
μ 1 α .
The Equation (19) for r α ( t ) becomes
r ˙ = μ tanh μ ( t f t ) + 2 y ˜ ( t ) , r ( t f ) = 0 ,
and finds the solution
r α ( t ) = 2 t t f cosh μ ( t f τ ) cosh μ ( t f t ) .
By substituting (26) and (29) into the optimal control (15), by integrating the differential Equation (22) for u = u α * ( t , x ) and by assuming that the function y ( t ) is twice continuously differentiable, the differentiation error is derived explicitly as
Δ α ( t ) u α * ( t , x ( t ) ) y ˙ ( y ) = μ sinh μ ( t f t ) cosh μ ( t f t ) Δ y 0 cosh ( μ t ) cosh ( μ t f ) y ˙ ( t f ) +
cosh ( μ t ) cosh ( μ t f ) 0 t f cosh μ ( t f τ ) w ( τ ) d τ 0 t cosh μ ( t τ ) w ( τ ) d τ ,
where
Δ y 0 = y ( 0 ) x 0 ,
w ( t ) = y ¨ ( t ) μ 2 η ( t ) .
By standard evaluations of the error function (30), the following theorem is established.
Theorem 1.
Let exist M 1 > 0 such that
| y ¨ ( t ) | M 1 , t [ 0 , t f ] .
Let the penalty coefficient α = α ( δ ) satisfy
lim δ 0 α ( δ ) = 0 ,
and exist δ 0 > 0 , M 2 > 0 , such that for δ δ 0 ,
δ α ( δ ) M 2 .
Then, for δ 0 , the realization u α ( t ) = u α * ( t , x ( t ) ) converges to y ˙ ( t ) uniformly on any interval
[ ε 1 , t f ε 2 ] [ 0 , t f ] , ε 1 , ε 2 > 0 .
In Figure 1, the results of the feedback differentiation are shown for y ( t ) = t 2 + t + 0.1 , t f = 2 . The realizations of the functions R α ( t ) and r α ( t ) were obtained by a backward numerical integration of differential Equations (25) and (28) with the integration step h = 2 N on the grid t k = k h , k = 0 , 1 , . . . , N , for N = 5000 . Then the Equation (22) was numerically integrated in a forward time. The measurement noise was realized at discrete time moments as random numbers η ( t k ) = η k uniformly distributed on [ δ , δ ] , k = 0 , 1 , . . . , N , for δ = 0.2 . The results are presented for x 0 = 0 and for three values of the penalty coefficient (regularization parameter) α = c δ ( c = 0.1 , 0.02 , 0.01 ) satisfying the conditions of Theorem 1. It is seen that at both ends of the interval [ 0 , t f ] , the error increases.
The error at the left-hand end of the interval [ 0 , t f ] is caused by the error in the initial value Δ y 0 0 . It can be proved that choosing x 0 = y ˜ ( 0 ) , for which | Δ y 0 | δ , guarantees the uniform convergence of u α ( t ) to y ˙ ( t ) for δ 0 on any interval [ 0 , t f ε ] [ 0 , t f ] , ε > 0 .
The error at the right-hand end of [ 0 , t f ] is caused by the nature of an optimal linear-quadratic control, satisfying u α ( t f ) = 0 independently of the function y ( t ) , the value of δ and the value of α . We see two possible approaches to circumvent this phenomenon. The first approach is to modify the cost functional by adding a terminal square term:
J ¯ α = | x ( t f ) x ¯ | 2 + 0 t f x ( t ) y ˜ ( t ) 2 d t + α 0 t f u 2 ( t ) d t ,
where x ¯ is an additional control parameter. The optimal control in the linear-quadratic problem (22), (36) given by
u ¯ α ( t , x ) = 1 2 α 2 R ¯ α ( t ) x + r ¯ α ( t ) ,
where R ¯ α ( t ) and r ¯ α ( t ) satisfy
R ˙ = 1 α R 2 1 , R ( t f ) = 1 ,
and
r ˙ = 1 α R α ( t ) r + + 2 y ˜ ( t ) , r ( t f ) = 2 x ¯ .
The second approach comprises extending the control interval from [ 0 , t f ] to 0 , t f + d for some d > 0 , by extrapolation of the noisy output y ˜ ( t ) for t [ t f , t f + d . In such a setting, u ¯ α ( t f + d ) = 0 and the right-hand error layer shifts to the right.
In Figure 2, the result of above mentioned improvements is presented for the same differentiation problem as in Figure 1 for α = 0.01 δ . In this case, the initial condition was x 0 = y ˜ ( 0 ) . The control interval was extended to [ 0 , 2.2 ] ( d = 0.2 ) by the linear extrapolation based on least squares fitting of y ˜ ( t ) at the interval [ 1.9 , 2 ] . The terminal target point was x ¯ = y ˜ ( 2.1 ) . It is seen that the differentiation is more accurate at both ends of the original interval.

3. Linear System Coefficients Restoration Problem

3.1. Problem Statement

Consider a linear homogeneous differential equation of order n
y ( n ) + u 1 ( t ) y ( n 1 ) + u 2 ( t ) y ( n 2 ) + . . . + u n 1 ( t ) y + u n ( t ) y = 0 , t [ 0 , t f ] ,
where u i ( t ) , i = 1 , . . . , n are the time-dependent coefficients. Let
y ( t ) = ( y 1 ( t ) , y 2 ( t ) , . . . , y n ( t ) ) T ,
be the vector of linearly independent solutions y i ( t ) , i = 1 , . . . , n , of (40). Substituting y i ( t ) , i = 1 , . . . , n , into (40) yields the linear algebraic system
H 0 ( t ) u ( t ) = y ( n ) ( t ) ,
w.r.t. to the coefficients vector
u ( t ) = ( u 1 ( t ) , u 2 ( t ) , . . . , u n ( t ) ) T ,
where
H 0 ( t ) = y 1 ( t ) y 2 ( t ) . . . y n ( t ) y 1 y 2 ( t ) . . . y n ( t ) y 1 ( n 1 ) ( t ) y 2 ( n 1 ) ( t ) . . . y n ( n 1 ) ( t ) T .
Note that det H 0 ( t ) is equal to the system Wronskian, i.e., det H 0 ( t ) 0 , and the system (42) admits a unique solution
u 0 ( t ) = H 0 1 ( t ) y ( n ) ( t ) .
However, the problem is to restore the coefficients (43) based on noisy output measurements
y ˜ ( t ) = y ( t ) + η ( t ) , t [ 0 , t f ] ,
where the error vector satisfies
| η ( t ) | δ , t [ 0 , t f ] ,
| y | denotes the euclidian norm of the vector y R n . The solution (45) is unstable w.r.t. the measurement noise, even if the exact derivatives of y ( t ) are available. Thus, the problem of coefficient restoration is ill-posed and requires regularization.
For the sake of illustration, consider the first-order equation
y ˙ + u ( t ) y = 0 ,
yielding the exact solution
u 0 ( t ) = y ˙ ( t ) y ( t ) .
Assume that the output y ( t ) is distorted by a high-frequency noise η ( t ) = δ sin ( ω t ) and an exact differentiation is available. Then, the exact solution becomes
u δ ( t ) = y ˜ ˙ ( t ) y ˜ ( t ) = y ˙ ( t ) + δ ω cos ( ω t ) y ( t ) + δ sin ( ω t ) ,
which, depending on the values δ and ω , can yield | | u ( · ) u 0 ( · ) | | 1 . In Figure 3, the functions (50) and (49) are shown for u 0 ( t ) = t ( y ( t ) = exp ( t 2 / 2 ) ), δ = 0.001 , ω = 1000 . It can be seen that the restoration error is very large.

3.2. Optimal Control Problem Reformulation

Let us introduce a state vector
x = ( x 1 , x 2 , . . . , x n 2 ) T R n 2 ,
where
x ( k 1 ) n + m = y k ( m 1 ) , k , m = 1 , . . . , n .
Then, for η ( t ) 0 , the system (42) can be rewritten as a bilinear system of differential equations
x ˙ = A x + B ( x ) u .
In the Equation (53), A is a block-diagonal n 2 × n 2 matrix
A = A 0 O O O O A 0 O O O O O A 0 ,
where A 0 is the n × n matrix
A 0 = 0 1 0 . . . 0 0 0 1 . . . 0 0 0 0 . . . 1 0 0 0 . . . 0 . ,
O denotes the n × n zero matrix. The n 2 × n matrix B ( x ) is
B ( x ) = B 1 ( x ) B 2 ( x ) B n ( x ) ,
where B k ( x ) are n × n matrices
B k ( x ) = 0 0 . . . 0 0 0 . . . 0 0 0 . . . 0 x ( k 1 ) n + n x ( k 1 ) n + n 1 . . . x ( k 1 ) n + 1 , k = 1 , . . . , n .
The output (41) is given by
y = F ( u ) = ( x 1 , x n + 1 , . . . , x ( n 1 ) n + 1 ) T = D x ,
where the n × n 2 matrix D is
D = 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 .
For example, the second-order differential equation
y ¨ + u 1 ( t ) y ˙ + u 2 ( t ) y = 0 ,
having two linearly independent solutions y 1 ( t ) and y 2 ( t ) , is represented as a bilinear system of four differential equations for the state vector x = ( x 1 , x 2 , x 3 , x 4 ) T = ( y 1 , y ˙ 1 , y 2 , y ˙ 2 ) T :
x ˙ 1 = x 2 x ˙ 2 = x 2 u 1 x 1 u 2 , x ˙ 3 = x 4 x ˙ 4 = x 4 u 1 x 3 u 2 ,
and the matrices A and B ( x ) in (53) are given by (54) and (56) with
A 0 = 0 1 0 0 , B 1 ( x ) = 0 0 x 2 x 1 , B 2 ( x ) = 0 0 x 4 x 3 .
The output is y = ( x 1 , x 3 ) T .
Let us formulate for the bilinear system (53) the tracking problem with the cost functional
J α = 0 t f | D x ( t ) y ˜ ( t ) | 2 d t + α 0 t f | u ( t ) | 2 d t ,
to be minimized by a feedback strategy u ( t , x ) .
Remark 3.
The cost (63) is the Tikhonov’s regularizing functional. Thus, if u * ( t , x ) is an optimal feedback in the bilinear-quadratic tracking problem (53), (63), then its time realization u * ( t ) = u * ( t , x ( t ) ) solves the problem of the coefficients restoration. The value of a regularization parameter α = α ( δ ) can be estimated, for example, by employing the discrepancy method [17].

4. Solution of Bilinear-Quadratic Tracking Problem

In this section, several approaches for solving the bilinear-quadratic tracking problem (53), (63) are presented.

4.1. Linear-Quadratic Approximation

In this approach, the matrix B ( x ) given by (56), is approximated by the matrix
B ˜ ( t ) = B ˜ 1 ( t ) B ˜ 2 ( t ) B ˜ n ( t ) ,
where
B ˜ k ( t ) = 0 0 . . . 0 0 0 . . . 0 0 0 . . . 0 y ˜ k ( n 1 ) ( t ) y ˜ k ( n 2 ) ( t ) . . . y ˜ k ( t ) , k = 1 , . . . , n .
The approximation (64) and (65) yields a linear system
x ˙ = A x + B ˜ ( t ) u ,
approximating the bilinear system (53). Thus, the bilinear-quadratic tracking problem (53), (63) is replaced by an approximating linear-quadratic problem (66), (63).
As it is seen from (65), this approach requires to restore the derivatives of the noisy data y ˜ ( t ) up to the order n 1 . In this paper, we use the linear-quadratic differentiation presented in Section 2.3. The advantage of this approach is that the solution of the linear-quadratic tracking problem is well-known. The optimal strategy is similar to (15):
u α * ( t , x ) = 1 2 α B ˜ ( t ) T Φ T ( t f , t ) 2 R α ( t ) Φ ( t f , t ) x + r α ( t ) ,
where the transition matrix Φ ( t , τ ) now corresponds to the matrix A given by (54), the matrix function R α ( t ) satisfies the Riccati differential equation (16), the vector function r α ( t ) satisfies the linear differential equation (19) with replacing b by B ˜ ( t ) , and D by D .

4.2. Iterative Solution

In this section, a modification of the iterative algorithm, proposed in [31] for a bilinear-quadratic problem with a slightly different cost functional, is applied. Assume that the function x ( 0 ) ( t ) is given. Then, at each iteration step, the linear-quadratic problem is solved for the system
x ˙ ( k + 1 ) = A x ( k + 1 ) + B ( x ( k ) ) u , k = 0 , 1 , . . .
with the cost functional
J α k + 1 = 0 t f | D x k + 1 ( t ) y ˜ ( t ) | 2 d t + α 0 t f | u k + 1 ( t ) | 2 d t , k = 0 , 1 , . . . .
As in (67), the optimal control in the problem (68) and (69) is
u α k + 1 ( t , x k + 1 ) = 1 2 α B ( x k + 1 ) T Φ T ( t f , t ) 2 R α k + 1 ( t ) Φ ( t f , t ) x k + 1 + r α k + 1 ( t ) ,
where the matrix function R α k + 1 ( t ) satisfies the Riccati
R ˙ = R Q α ( k ) ( t ) R S ( k ) ( t ) , R ( t f ) = 0 ,
where
Q α ( k ) ( t ) = 1 α Φ ( t f , t ) B ( x ( k ) ) B T ( x ( k ) ) Φ T ( t f , t )
S ( k ) ( t ) = Φ T ( t , t f ) D T D Φ ( t , t f ) .
The vector function r α ( k + 1 ) ( t ) satisfies the linear differential equation
r ˙ = R α ( t ) Q α ( k ) ( t ) r + 2 Φ T ( t , t f ) D T y ˜ ( t ) , r ( t f ) = 0 .
The convergence of this iterative procedure by an appropriate choosing α = α ( δ ) was proved in [31].

4.3. Feedback Linearization

In this section, the bilinear-quadratic control problem is treated by using feedback linearization [32]. The relative degree of the system (53) w.r.t. to each control u i , i = 1 , . . . , n , is n, yielding the full relative degree n 2 [36] w.r.t. to the vector control u. This implies that the system is exactly feedback linearizabile. Namely, define the auxiliary controls
U k = x ( k 1 ) n + n u 1 x ( k 1 ) n + n 1 u n 1 . . . x ( k 1 ) n + 1 u 1 , k = 1 , . . . , n ,
or,
U = B ¯ ( x ) u ,
where U = ( U 1 , . . . , U n ) T ,
B ¯ ( x ) = x n x n 1 . . . x 1 x 2 n x 2 n 1 . . . x n + 1 0 0 . . . 0 x n 2 x n 2 1 . . . x n 2 n + 1 .
Then, by denoting z = x , the system (53) is rewritten as
z ˙ = A z + B ^ U ,
where
B ^ = O 1 0 . . . 0 0 O 0 1 . . . 0 0 O 0 0 . . . 0 1 .
Let us formulate the tracking problem for the linear system (78) with the cost functional
J ^ α = 0 t f | D z ( t ) y ˜ ( t ) | 2 d t + α 0 t f | U ( t ) | 2 d t .
The optimal strategy in the linear-quadratic problem (78), (80) is similar to (67):
U ^ α ( t , z ) = 1 2 α B ^ T Φ T ( t f , t ) 2 R ^ α ( t ) Φ ( t f , t ) z + r ^ α ( t ) ,
where the matrix function R ^ α ( t ) satisfies the Riccati differential Equation (16), the vector function r ^ α ( t ) satisfies the linear differential Equation (19) with replacing b by B ^ , and D by D . Due to (76), the optimal control in the bilinear-quadratic problem is
u ^ α ( t , x ) = B ¯ 1 ( x ) U ( t , x ) .
Remark 4.
The state variables x ( k 1 ) n + 1 , k = 1 , . . . , n , track the functions y k ( t ) , k + 1 , . . . , n . Then, the value of | det ( B ¯ ( x ) ) | should be close to the wronskian of the linear differential Equation (40), providing the invertibility of the matrix B ¯ ( x ) . However, due to the measurement noise in y ( t ) and the tracking inaccuracy, this matrix can become singular. In order to avoid inverting a singular matrix, in practice we use
u ^ α ( t , x ) = B ¯ ( x ) + k I n 1 U ( t , x ) ,
where I n is the unit matrix, the parameter k > 0 is small enough.

4.4. Numerical Example

Let us consider for example the differential equation
y ¨ 2 ( t 1 ) y ˙ + ( t 1 ) 2 y = 0 ,
which linearly independent solutions are y 1 = exp ( t ( t 2 ) / 2 ) cos t , y 2 = exp ( t ( t 2 ) / 2 ) sin t . In Figure 4 and Figure 5, the results of the coefficients restoration by the linear-quadratic approximation are presented for t f = 1 , δ = 0.001 , α = 10 7 . The derivatives y ˙ 1 , y ˙ 2 were restored by using the linear-quadratic differentiation described in Section 2.3.
The respective restoration results by using the iterative procedure (for 5 iterations) and feedback linearization are shown in Figure 6, Figure 7, Figure 8 and Figure 9, respectively.
Similarly to the differentiation, one can observe the error at the left-hand end ( t = 0 ) caused by the inaccuracy in the initial condition.
Let u * ( t ) = ( u 1 * ( t ) , u 2 * ( t ) ) T denote the optimal control time realization, approximating the coefficients u ( t ) . The restoration L 2 -errors
ε i = 0 t f u i ( t ) u i * ( t ) 2 d t 1 / 2 ,
are presented in Table 1. It is seen that the most accurate restoration of u 1 ( t ) and u 2 ( t ) was obtained by using the iterative algorithm and the feedback linearization, respectively, i.e. by the approaches tackling the original bilinear-quadratic problem.

5. Conclusions

Two inverse ill-posed problems ((i) input restoration of a linear system and (ii) coefficient restoration of a linear homogeneous differential equation) were considered. Both problems were tackled in a unified manner by reformulating as an optimal control problem with Tikhonov’s cost functional. In both cases, the solution is represented by a feedback control, which adds to the robustness of a proposed approach. The first problem was reformulated as a linear-quadratic tracking problem. A specific attention was paid to a practically important differentiation problem, which solution was further employed in the coefficient restoration. By choosing an initial condition and by modifying the cost functional, the errors at the ends of the time interval were partially circumvented.
The coefficient restoration problem was reformulated as a bilinear-quadratic tracking problem. Three approaches for its solution were proposed: (i) a linear-quadratic approximation, (ii) an iterative method, and (iii) a feedback linearization. In the numerical example, it is shown that the most accurate restoration was obtained by the latter two approaches, tackling an original bilinear-quadratic control problem. The future research should be concentrated in decreasing the errors occurring at the ends of the time interval.
The following abbreviations are used in this manuscript:

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Solopchenko, G.N. Inverse problems in measurement. Measurement 1987, 5, 10–19. [Google Scholar] [CrossRef]
  2. Solopchenko, G.N. Formal metrological components of measuring systems. Measurement 1994, 13, 1–12. [Google Scholar] [CrossRef]
  3. Aronov, P.M.; Lyashenko, E.A.; Ryashko, L.B. Method of optimum linear-estimation for determining dynamic characteristics of measuring devices. Meas. Tech. 1991, 34, 1091–1096. [Google Scholar] [CrossRef]
  4. Granovskii, V.A. Models and methods of dynamic measurements: Results presented by St. Petersburg metrologists. In Advanced Mathematical and Computational Tools in Metrology and Testing X; Pavese, F., Bremser, W., Chunovkina, A., Fischer, N., Forbes, A., Eds.; World Scientific Publishing Company: Singapore, 2015; Volume 86, pp. 29–37, Series on Advances in Mathematics for Applied Sciences. [Google Scholar]
  5. Frieden, B.R. Image enhancement and restoration. In Picture Processing and Digital Filtering; Huang, T.S., Ed.; Springer: New York, NY, USA, 1979. [Google Scholar]
  6. Bates, R.H.T.; McDonnell, J. Image Restoration and Reconstruction; The Oxford Engineering Science Series; Clarendon Press: Oxford, UK, 1986; Volume 16. [Google Scholar]
  7. Jansson, P. Deconvolution: With Applications in Spectroscopy; Academic Press: New York, NY, USA, 1984. [Google Scholar]
  8. Claerbout, J.F. Imaging the Earth’s Interior; Blackwell Scientific Publications: Oxford, UK, 1985. [Google Scholar]
  9. Prutkin, I. Gravitational and magnetic models of the core-mantle boundary and their correlation. J. Geodyn. 2008, 45, 146–153. [Google Scholar] [CrossRef] [Green Version]
  10. Beck, M.B. Applications of System Identification and Parameter Estimation in Water Quality Modeling; IIASA Working Paper WP-79-099; IIASA: Laxenburg, Austria, 1979. [Google Scholar]
  11. Isaksson, A.J. Some aspects of industrial system identification. IFAC Proc. Vol. 2013, 46, 153–159. [Google Scholar] [CrossRef] [Green Version]
  12. Worden, K.; Barthorpe, R.J.; Cross, E.J.; Dervilis, N.; Holmes, G.R.; Manson, G.; Rogers, T.J. On evolutionary system identification with applications to nonlinear benchmarks. Mech. Syst. Signal Process. 2018, 112, 194–232. [Google Scholar] [CrossRef]
  13. Bellman, R. Mathematical Methods in Medicine; Modern Applied Mathematics; World Scientific: Singapore, 1983; Volume 1. [Google Scholar]
  14. Schwilden, H.; Honerkamp, J.; Eister, C. Pharmacokinetic model identification and parameter estimation as an ill-posed problem. Eur. J. Clin. Pharmacol. 1993, 45, 545–550. [Google Scholar] [CrossRef] [PubMed]
  15. Marinov, T.T.; Marinova, R.S. Dynamics of COVID-19 using inverse problem for coefficient identification in SIR epidemic models. Chaos Solitons Fractals X 2020, 5, 1–15. [Google Scholar] [CrossRef]
  16. Korolev, I. Identification and estimation of the SEIRD epidemic model for COVID-19. J. Econom. 2021, 220, 63–85. [Google Scholar] [CrossRef] [PubMed]
  17. Tikhonov, A.N.; Arsenin, V.I. Solution of Ill-Posed Problems; John Wiley & Sons: New York, NY, USA, 1977. [Google Scholar]
  18. Bakushinsky, A.; Goncharsky, A. Ill-Posed Problems: Theory and Applications; Kluwer: Dordrecht, The Netherlands, 1994; Volume 301, Mathematics and Its Applications; [Google Scholar]
  19. Ivanov, V.; Vasin, V.; Tanana, V. Theory of Linear Ill-Posed Problems and Its Applications; Inverse and Ill-Posed Problems Series; Walter De Gruyter: Berlin, Germany, 2002; Volume 36. [Google Scholar]
  20. Balakrishnan, A.V. Applied Functional Analysis; Springer: New York, NY, USA, 1976. [Google Scholar]
  21. Wang, Y.; Yagola, A.M.; Yang, C. (Eds.) Optimization and Regularization for Computational Inverse Problems and Applications; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  22. Kryazhimskiy, A.V.; Osipov, Y.S. Inverse Problems for Ordinary Differential Equations: Dynamical Solutions; Taylor and Francis: Basel, Switzerland, 1995. [Google Scholar]
  23. Abidi, M.A.; Gribok, A.V.; Paik, J. Frequency-domain implementation of regularization. In Optimization Techniques in Computer Vision: Ill-Posed Problems and Regularization; Springer International Publishing: Cham, Switzerland, 2016; pp. 113–130. [Google Scholar]
  24. Ljung, L. System Identification (2nd Ed.): Theory for the User; Prentice Hall PTR: Upper Saddle River, NJ, USA, 1999. [Google Scholar]
  25. Perona, P.; Porporato, A.; Ridolfi, L. On the trajectory method for the reconstruction of differential equations from time series. Nonlinear Dyn. 2000, 23, 13–33. [Google Scholar] [CrossRef]
  26. Maksimov, V.I. Dynamical identification of unknown characteristics in systems of the second order. In Proceedings of the 2001 European Control Conference, Porto, Portugal, 4–7 September 2001; pp. 1261–1265. [Google Scholar]
  27. Subbotina, N.N.; Tokmantsev, T.B.; Krupennikov, E.A. On the solution of inverse problems of dynamics of linearly controlled systems by the negative discrepancy method. Proc. Steklov Inst. Math. 2015, 291, 253–262. [Google Scholar] [CrossRef]
  28. Subbotina, N.N.; Krupennikov, E.A. The method of characteristics in an identification problem. Proc. Steklov Inst. Math. 2017, 299, 205–216. [Google Scholar] [CrossRef]
  29. Turetsky, V.Y.; Lemsh, E.S. A feedback control in linear system identification. IFAC Proc. Vol. 1996, 29, 1–4. [Google Scholar] [CrossRef]
  30. Jha, B.; Turetsky, V.; Shima, T. Robust path tracking by a Dubins ground vehicle. IEEE Trans. Control Syst. Technol. 2019, 27, 2614–2621. [Google Scholar] [CrossRef]
  31. Hofer, E.; Tibker, B. An iterative method for the finite-time bilinear-quadratic control problem. J. Optim. Theory Appl. 1988, 57, 411–427. [Google Scholar] [CrossRef]
  32. Isidori, A. Nonlinear Control Systems, 2nd ed.; Springer-Verlag: London, UK, 1989. [Google Scholar]
  33. Turetsky, V. Robust trajectory tracking for a feedback linearizable nonlinear system. IFAC Pap. 2016, 49, 540–545. [Google Scholar] [CrossRef]
  34. Mellodge, P. A Practical Approach to Dynamical Systems for Engineers; Woodhead Publishing Series in Mechanical Engineering; Elsevier: Amsterdam, The Netherlands, 2016. [Google Scholar]
  35. Kwakernaak, H.; Sivan, R. Linear Optimal Control Systems; Wiley-Interscience: New York, NY, USA, 1972. [Google Scholar]
  36. Slotine, J.; Li, W. Applied Nonlinear Control; Prentice Hall: Englewood Cliffs, NJ, USA, 1991. [Google Scholar]
Figure 1. Feedback differentiation of y ( t ) = t 2 + t + 0.1 .
Figure 1. Feedback differentiation of y ( t ) = t 2 + t + 0.1 .
Axioms 10 00137 g001
Figure 2. Improved feedback differentiation of y ( t ) = t 2 + t + 0.1 .
Figure 2. Improved feedback differentiation of y ( t ) = t 2 + t + 0.1 .
Axioms 10 00137 g002
Figure 3. Exact solution is unstable w.r.t. measurement noise.
Figure 3. Exact solution is unstable w.r.t. measurement noise.
Axioms 10 00137 g003
Figure 4. Restoration of u 1 ( t ) by linear-quadratic approximation.
Figure 4. Restoration of u 1 ( t ) by linear-quadratic approximation.
Axioms 10 00137 g004
Figure 5. Restoration of u 2 ( t ) by linear-quadratic approximation.
Figure 5. Restoration of u 2 ( t ) by linear-quadratic approximation.
Axioms 10 00137 g005
Figure 6. Restoration of u 1 ( t ) by iterations.
Figure 6. Restoration of u 1 ( t ) by iterations.
Axioms 10 00137 g006
Figure 7. Restoration of u 2 ( t ) by iterations.
Figure 7. Restoration of u 2 ( t ) by iterations.
Axioms 10 00137 g007
Figure 8. Restoration of u 1 ( t ) by feedback linearization.
Figure 8. Restoration of u 1 ( t ) by feedback linearization.
Axioms 10 00137 g008
Figure 9. Restoration of u 2 ( t ) by feedback linearization.
Figure 9. Restoration of u 2 ( t ) by feedback linearization.
Axioms 10 00137 g009
Table 1. L 2 -errors of restoration.
Table 1. L 2 -errors of restoration.
u 1 ( t ) u 2 ( t )
LQ-approximation 0.07 0.08
Iterative solution 0.008 0.07
Feedback linearization 0.07 0.03
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Turetsky, V. Two Inverse Problems Solution by Feedback Tracking Control. Axioms 2021, 10, 137. https://doi.org/10.3390/axioms10030137

AMA Style

Turetsky V. Two Inverse Problems Solution by Feedback Tracking Control. Axioms. 2021; 10(3):137. https://doi.org/10.3390/axioms10030137

Chicago/Turabian Style

Turetsky, Vladimir. 2021. "Two Inverse Problems Solution by Feedback Tracking Control" Axioms 10, no. 3: 137. https://doi.org/10.3390/axioms10030137

APA Style

Turetsky, V. (2021). Two Inverse Problems Solution by Feedback Tracking Control. Axioms, 10(3), 137. https://doi.org/10.3390/axioms10030137

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