1. Introduction
In [
1], an algorithm was constructed for estimating the parameters of a system of first-order ordinary differential equations by a large number of inaccurate observations in the vicinity of one selected point. We will further call such points reference points. Estimates of the solution of the system and its derivative at the reference point are constructed by the method of linear regression analysis. According to them, estimates of the system parameters are determined by the method of moments, and their consistency is proved. However, the development of this topic requires the selection of several reference points for estimating the parameters of a linear system of ordinary differential equations. Such a problem arises, for example, if the number of unknown parameters is greater than the number of equations in the system and if the order of the system is higher than the first.
The estimation of the parameters of a system of linear differential equations with constant coefficients based on inaccurate observations of its solution is of particular interest. The choice of reference points forming a finite arithmetic progression allows us, in this case, to construct ratios for evaluating the elements of the coefficient matrix, including the inverse matrix, the square root of the matrix, etc. [
2,
3,
4]. At the same time, it turns out that the Vandermonde determinant plays the main role in the circulation of the matrices considered in the paper. The reversibility of the Vandermonde matrix requires the realness and difference of the eigenvalues of the coefficients matrix in the system of equations.
However, consideration of a system of first-order differential equations is insufficient in the study of oscillatory dynamical systems, described by systems of the second order. Meanwhile, the extension of the problem of estimating the parameters of systems by inaccurate data to the class of oscillatory systems is of interest. For example, it is an inverse problem of oscillations theory [
5,
6] when it is necessary to restore the parameters of a distributed system by one or another full-scale data. This formulation of the question is closely related to the models of mechanical systems in the problems of mechatronics and robotics, which have received great development and dissemination in recent years [
7,
8,
9,
10,
11]. This task is also closely related to the problems of technical systems reliability, containing similar mechanical components. If we talk about reliability, then the issues related to the possibility of resonance [
12,
13,
14,
15,
16] are of particular interest. Therefore, it is desirable to obtain more accurate estimates of the parameters of such differential equations systems, especially frequencies, which play an important role in the analysis of resonant phenomena.
Direct reduction in general oscillatory systems to first-order systems (in the case of linear systems with constant coefficients) leads to matrices whose spectrum can contain multiples of eigenvalues. The reversal of the matrices arising, in this case, requires additional research. In the case of conservative mechanical systems, this leads to a matrix inversion algorithm based on the Vandermonde determinant.
In this paper, we consider two problems of estimating the parameters of systems of first and second-order differential equations using inaccurate observations from analytical and computational points of view. The solution to both problems is reduced to matrix calculations, in which the reversibility of the corresponding matrices plays an important role. Moreover, an important element of solving these problems is the choice of reference points in the vicinity of which numerous measurements are made. These points form an arithmetic progression, which makes it possible to use the Vandermonde determinant to estimate the parameters of the models under consideration. For models of mechanical systems, this problem is solved by switching from a system of second-order differential equations to a system of first-order equations and estimating the parameters of this system by inaccurate observations. Analytical calculations are supplemented by computational experiments confirming the possibility of using the proposed methods.
The main result of the work is the selection of reference points for estimating parameters, an analytical study of the constructed estimates, and a computational experiment to determine the errors of the estimates obtained. These results can be used to analyze vibrations in mechatronics and robotics systems and to determine the reliability of these systems in terms of their protection from resonance. Estimates of matrices containing solutions of systems at reference points and matrices containing derivatives of solutions at reference points were constructed. Estimates of the coefficient matrices of the systems under consideration were based on them. At the same time, the conditions for the reversibility of the matrices were established using the Vandermonde formula. The constructed estimates of the coefficient matrices are consistent. This is confirmed in the course of computational experiments. The method of parameter estimation proposed in the article allows us not only to build sufficiently accurate estimates of parameters but also to control this accuracy by choosing the necessary number of observations and reference points in the vicinity of which observations are carried out, i.e., ensuring that these estimates are consistent (converging in probability with an increase in the number of observations).
It is worth mentioning that this paper is devoted to the estimation of the frequencies and amplitudes of linear oscillations in a conservative system. This task is complex and requires both the study of the properties of the differential equations themselves, describing a conservative system, and statistical estimates of the oscillation parameters. A large number of papers have been devoted to the study of the properties of such differential equations. Among them, one should point to the classical monograph [
2], in which linear Lagrange equations, describing a conservative system, are solved by methods of matrix theory. As for the estimates of the parameters of differential equations of a fairly general form, we should point to the monographs [
17,
18,
19] and articles developing this direction (see, for example, [
20,
21,
22,
23]). They give estimates of the parameters of differential equations using least squares error minimization between the response of the model and the actual response of the system. In some cases, the asymptotic normality of the obtained estimates is established.
In this paper, the question of adapting such estimates to differential equations, describing conservative systems, is raised. To do this, we had to move from differential equations in Lagrange variables, studied in [
2], to differential equations in Hamilton variables. This made it possible to obtain solutions to the equations in the form of exponential matrices and to establish the conditions for the reversibility of the matrices used to construct statistical estimates. As for the properties of statistical estimates, by constructing consistent estimates of the values of the observed functions and their derivatives, using linear regression analysis for a large number of observations in the vicinity of selected (reference) points, it is possible to construct estimates of matrices, included in linear differential equations, describing conservative systems, and with their help to estimate the frequencies and amplitudes of oscillations.
2. Systems of Linear Differential Equations of the First Order
Consider a system of linear differential equations of the form
Suppose that the time points
are given (let us call these moments reference points) and then the equalities are fulfilled
Let inaccurate observations be made in the vicinity of each reference point. It is required to construct consistent estimates (converging in probability at to the estimated parameter) of all elements of the matrix according to the constructed consistent estimates of matrix elements
Let us rewrite the system of Equalities (
2) in matrix form
Then, the matrix
A can be reconstructed using the matrix relation
by matrices
, if the matrix
Y is revisable. Let us find out the conditions under which the matrix
Y is revisable.
Suppose that the matrix
A has
m of various real eigenvalues
and their corresponding eigenvectors-columns
forming a basis in
m-dimensional space. Let us rewrite the system of differential Equation (
1) in the basis
assuming
where
Here, the multipliers
are components of the column vector
defining the initial conditions for a system of differential Equation (
5).
Let us define a set of reference points
for some
by the relations
Let us define the square matrix
by the equality
Then, the matrix
Y for
satisfies the following relations:
From the Formula (
6) and the formula for calculating the Vandermond determinant (see [
24], for example), we obtain
Since the eigenvalues are pairwise different and , when performing the relations , the matrix Y is reversible.
Suppose that at time points
inaccurate observations are obtained. We construct estimates
for matrix elements
Here,
are independent identically distributed random variables with zero mean and finite variance. For
in [
1], convergence in probability (and hence the consistency of estimates) is proved
Denote the matrices
Each element of the matrix
(matrix element
) is a consistent estimate of the corresponding matrix element
(matrix element
Y). The inverse matrix
coincides with the transposed matrix of algebraic complements corresponding to the elements of the matrix
Y divided by the determinant of the matrix
Y (see, for example, [
24]). Therefore, each element of the matrix
is a consistent estimate of the corresponding element of the matrix
It follows that each element of the matrix
is a consistent estimate of the corresponding element of the matrix
Computational Experiment
A computational experiment was conducted for the Cauchy problem
The solution of this system of equations has the form
We assumed that by observing the process described by this system of equations, inaccurate observations were obtained at time points
,
100,000,
Here,
are independent random variables distributed uniformly on the segment
Using Formula (
9), the matrix
A is evaluated by the following matrices
Next, a computational experiment was conducted for the Cauchy problem
The solution of this system of equations has the form
We assumed that by observing the process described by this system of equations, inaccurate observations were obtained at time points
,
100,000,
where
are independent random variables distributed uniformly on the segment
As a result
It should be noted that for the proposed version of the parameter estimation, it is necessary to establish the conditions for the reversibility of the matrices involved in this assessment and correctly select the step of the arithmetic progression between the reference points and the step between neighboring observation points in the vicinity of the reference points.
3. Equations of Oscillations of a Conservative System
Consider a mechanical conservative system described by a vector of generalized coordinates
with kinetic energy
and potential energy
. We assume that the matrices
,
are matrices of positive definite quadratic forms and, therefore, reversible. Then, the Lagrange equation for the system under consideration has the form
or
Suppose that the eigenvalues of the matrix K are positive and distinct. It should be noted that the requirement of the uniqueness of the eigenvalues of the matrix K can be relaxed.
Let us put
, using the transition from Lagrangian variables
q to Hamiltonian variables
we obtain
Then, the system (
10) will take the form (
1). Its solution is described by a matrix exponent:
Let us write down this solution with respect to the basis of the eigenvectors of the matrix A. To do this, we denote by the eigenvalues of the matrix . The numbers are the roots of the equation and satisfy the equalities Let us put and .
Denote
R and
T transition matrices as the basis of the eigenvectors of the matrices
K and
, respectively. We can consider
R to be real since such are the eigenvalues of the matrix
K. By definition of the transition matrix,
and
. A non-trivial solution to these equations is
. Indeed, let us check the equality
and, therefore,
Elementary calculations check the equality
Obviously, all the introduced functions from matrices are diagonal. It is not difficult to verify that using the equality
and Formula (
12), it is possible to rewrite Formula (
11) as
and, then,
Let us now investigate the reversibility of the matrix
,
Using (
13), it is easy to check that
where
is a column of units. Thus, the reversibility of
Y is determined by the reversibility of the matrices
and
. The matrix
is reversible if and only if all its diagonal elements are not zero, i.e., when all components of the vector
are non-zero (vector
with the components
and
, subjected to the specified condition, is called the general position vector). In turn, the matrix
is a Vandermonde matrix, the determinant of which is
is not zero if and only if all exponents
are different. If
where
then all the exponents will be different.
By analogy with Formulas (
8) and (
9), it is now possible to construct estimates of the matrix
A and estimates of the matrix
K based on them. In turn, the estimates of the matrix
K allow us to construct estimates of the frequencies
and coefficients for
Computational Experiment for a System of Two Coupled Pendulums
As an example, consider a system of two identical connected linear pendulums (
Figure 1).
In this case, the kinetic and potential energies are given by the expressions
Here,
m—mass of each pendulum,
l—length,
and
—angles of deviation from the vertical axis of the first and second pendulums, respectively,
—stiffness of the spring connecting the pendulums,
h—the distance from the pendulum suspension point to the spring attachment point. Then
Solving the problem of eigenvalues and eigenvectors for the matrix
K, we find
Substituting these expressions into (
13), we obtain
The conditions for the reversibility of the observation matrix are formulated as follows. The vector of the initial data of the general position is an arbitrary vector of the form
with non-zero pairs
and
. For example, such is the vector
. Meanwhile, the vector
is not a vector of the general position. When choosing a general position vector as the initial data vector, reversibility obviously takes place under the condition
. To conduct a computational experiment, we choose
then
Estimates of
(see formula (
8)) can be used to construct estimates of matrix elements
and so estimates of matrix elements
Therefore, it is possible to estimate the frequencies
To perform the calculations, we additionally put
then
We assume that by observing the process, described by this system of equations, inaccurate observations are obtained at time points
100,000,
where
are independent random variables distributed uniformly on the segment
The segment
characterizes the spread of random observation errors and their variation. As a result, we have
Thus, we obtained estimates
of frequencies
It is interesting to investigate the oscillation amplitudes (multipliers at
) of the components of the vector
To do this, using Formula (
14), we obtain the equality
in which the matrix
F is representable as
Let us now construct an estimate of
of the matrix
replacing the values
,
in it by their estimates
Then the elements of the matrix
determine estimates of oscillation amplitudes with frequencies
For
, these estimates converge in probability to the corresponding elements of the matrix
In our case
4. Discussion
When analyzing a system of differential equations, describing oscillations in a conservative mechanical system, a transition is used from (Lagrangian) variables included in a system of second-order differential equations to (Hamiltonian) variables included in a system of first-order differential equations. As a result of such a replacement, it becomes possible to use exponents from diagonal matrices. This circumstance makes it possible to establish conditions for the reversibility of the matrix used in estimating the parameters of a conservative system. During the computational experiment, it was found that the step of the arithmetic progression, which determines the set of reference points (in the vicinity of which observations are carried out), significantly affects the accuracy of the estimates obtained. In turn, the estimation of the parameters of a conservative mechanical system is mainly needed to estimate the oscillation frequencies. These frequencies allow us to determine how close this mechanical system is to resonance, and, therefore, what is its stability and reliability.
5. Conclusions
The paper selects a set of reference points in the form of an arithmetic sequence for estimating the parameters of a system of ordinary differential equations, based on inaccurate observations. This makes it possible to obtain conditions for the reversibility of the matrix involved in estimating the parameters of a system of first-order differential equations. The estimation of the parameters of a mechanical conservative system, based on inaccurate observations, is carried out by switching to a system of first-order differential equations. This transition is made by replacing Lagrangian variables describing the analyzed conservative system with Hamiltonian variables. As a result, the reversibility of the matrix used in parameter estimation is established using the Vandermonde formula. A rather voluminous computational experiment has been carried out showing how, with an increase in the number of observations in the vicinity of reference points, the accuracy of parameter estimates of the analyzed systems increases. The influence of the arithmetic progression step, determining the set of reference points on the accuracy of parameter estimates, is investigated. The problem of estimating the characteristics of oscillations in a conservative system, considered in this paper, is multi-parametric. Therefore, to solve it, it is necessary to combine the methods of the theory of ordinary differential equations with constant coefficients and methods for constructing consistent estimates of coefficients based on the definition of reference points for planning an experiment. This circumstance makes it possible to construct modern measuring systems based on strict mathematical methods and to increase the accuracy of determining the parameters of the conservative system based on inaccurate observations. It should be noted that the proposed method of parameter estimation can be applied to systems consisting of a large number of elements since it is based on classical calculations of such systems, using Lagrange equations and matrix methods. However, in order to construct and study the quality of the estimates obtained, it is necessary to switch from Lagrange variables q (generalized coordinates) to generalized Hamilton variables, .