Next Article in Journal
Maximum Correntropy with Variable Center Unscented Kalman Filter for Robust Power System State Estimation
Next Article in Special Issue
Some Further Results on the Fractional Cumulative Entropy
Previous Article in Journal
New, Spherical Solutions of Non-Relativistic, Dissipative Hydrodynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Positive Definite Kernels of Integral Operators Corresponding to the Boundary Value Problems for Fractional Differential Equations

by
Mukhamed Aleroev
1 and
Temirkhan Aleroev
2,*
1
Department of Higher Mathematics, Moscow Automobile and Road Construction State Technical University (STU-MADI), Leningradsky Ave. 64, 125319 Moscow, Russia
2
Department of Applied Mathematics, The National Research Moscow State University of Civil Engineering (NRU MGSU), Yaroslavskoe Highway 26, 129337 Moscow, Russia
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(4), 515; https://doi.org/10.3390/e24040515
Submission received: 23 February 2022 / Revised: 31 March 2022 / Accepted: 31 March 2022 / Published: 6 April 2022

Abstract

:
In the spectral analysis of operators associated with Sturm–Liouville-type boundary value problems for fractional differential equations, the problem of positive definiteness or the problem of Hermitian nonnegativity of the corresponding kernels plays an important role. The present paper is mainly devoted to this problem. It should be noted that the operators under study are non-self-adjoint, their spectral structure is not well investigated. In this paper we use various methods to prove the Hermitian non-negativity of the studied kernels; in particular, a study of matrices that approximate the Green’s function of the boundary value problem for a differential equation of fractional order is carried out. Using the well-known Livshits theorem, it is shown that the system of eigenfunctions of considered operator is complete in the space L 2 ( 0 , 1 ) . Generally speaking, it should be noted that this very important problem turned out to be very difficult.

1. Introduction

In the theory of oscillations, the general physical process of reciprocity (when the source and the observer can exchange places) finds its mathematical expression in the self-adjointness of the corresponding boundary value problems [1,2]. When we talk about processes corresponding to structures with fractal geometry, the corresponding processes are described by non-self-adjoint operators, and here, first of all, we are talking about integral operators with a persymmetric kernel [1,2]. The spectral structure of such operators has hardly been studied. The present manuscript is devoted to this question.

2. Main Results

In the present paper we consider the following operator
A ρ u = 1 Γ ( ρ 1 ) 0 x ( x t ) 1 ρ 1 u ( t ) d t 0 1 ( x x t ) 1 ρ 1 u ( t ) d t ,   1 / 2 < ρ < 1
Corresponding to the following problem [3]
D 1 / ρ u = 1 Γ ( ρ 1 ) d d x 0 x u ( t ) d t ( x t ) 1 / ρ 1 = λ u ,
u ( 0 ) = u ( 1 ) = 0 .
By the Livshits theorem [2,4], the system of eigenfunctions of this operator is complete in L 2 (0, 1); precisely, we have
Theorem 1.
(Livshits): if K ( x , y ) , ( a x , y b ) is a bounded kernel, and the “real part” 1 2 ( K + K * ) of it is non-negative kernel, then the inequality holds
j = 1 R e ( 1 / λ j ) a b R e K ( t , t ) d t ,
where λ j is the characteristic numbers of kernel K. The system of main functions of the kernel K is complete in the domain of values of the integral operator K f if and only if, when there is an equal sign in the inequality above.
Here, in this theorem characteristic numbers are eigenvalues. It is known that the operator A is called positive definite (the definition of positivity of the operator and its properties can be found in [5]) if ( A u , u ) > 0 , ( u 0 ) . However, it is very difficult to verify this condition directly. Therefore, we will use the matrix approximation of the operator A ρ [1]. As in [1], we denote the corresponding matrix by T n 1 ( μ ) , μ = 1 ρ 1
T n 1 ( μ ) = ( 1 n ) μ ( n 1 n ) μ ( 1 n ) μ ( n 2 n ) μ ( 1 n ) μ ( 1 n ) μ ( 2 n ) μ ( n 1 n ) μ ( 1 n ) μ ( 2 n ) μ ( n 2 n ) μ ( 2 n ) μ ( 1 n ) μ ( n 1 n ) μ ( n 1 n ) μ ( n 1 n ) μ ( n 1 n ) μ ( n 2 n ) μ ( n 3 n ) μ ( n 1 n ) μ ( 1 n ) μ .
The matrix T n 1 ( μ ) has many useful properties. In particular, this matrix is positive, persymmetric, indecomposable, etc. It is known [1] that one of the necessary conditions for the positive definiteness of a matrix is the positivity of all its lead main minors. The fact that these minors are positive was shown in [1,6]. We give a detailed proof of one theorem from which the above follows.
Theorem 2.
The minors
A i 1 i 2 i r j 1 j 2 j r
of the matrix
T n 1 ( μ ) = ( 1 n ) μ ( n 1 n ) μ ( 1 n ) μ ( n 2 n ) μ ( 1 n ) μ ( 1 n ) μ ( 2 n ) μ ( n 1 n ) μ ( 1 n ) μ ( 2 n ) μ ( n 2 n ) μ ( 2 n ) μ ( 1 n ) μ ( n 1 n ) μ ( n 1 n ) μ ( n 1 n ) μ ( n 1 n ) μ ( n 2 n ) μ ( n 3 n ) μ ( n 1 n ) μ ( 1 n ) μ .
for i k j k , 1 k r , are positive. Moreover, they are equal to
( n μ ) r 1 ( n j r ) μ ( i 1 ) μ b r ( r 1 ) b ( r 1 ) ( r 2 ) b 21 ,
where
b k i = ( i k j i ) μ , i k > j i 0 , i k < j i
Proof. 
Let us consider the minor
M r = A i 1 i 2 i r j 1 j 2 j r .
For i k j k , 1 k r , we to overwrite M r as follows
M r = i 1 μ i 2 μ i 3 μ i r 1 μ ( ( n j 1 ) μ ( n j 2 ) m u ( n j r ) μ ) n μ 0 0 0 0 b 21 0 0 0 b r 1 b r 2 b r 3 b r ( n 1 )
for
b k i = ( i k j i ) μ , i k > j i 0 , i k < j i
To calculate the determinant, M r we consider
det ( M r λ I ) = ( 1 ) r det ( λ I M r ) = ( 1 ) r det ( A ˜ x y T ) = ( 1 ) r λ r ( 1 y T A ˜ 1 x ) .
Here,
A ˜ = n μ λ n μ 0 0 0 b 21 λ n μ 0 0 b r 1 b r 2 b r 3 λ n μ ,
x = i 1 μ i 2 μ i 3 μ i r 1 μ , y T = ( ( n j 1 ) μ ( n j 2 ) μ ( n j r ) μ ) .
It is clear that
x r 1 = ( n μ λ ) r 1 b r ( r 1 ) b ( r 1 ) ( r 2 ) b 21 x 1 = = ( 1 ) r 1 ( n μ λ ) r 1 b r ( r 1 ) b ( r 1 ) ( r 2 ) b 21 +
So,
det ( M r λ I ) = ( 1 ) r λ r ( 1 y T A ˜ 1 x ) = ( 1 ) r λ r ( 1 ( n j r ) μ i 1 r x r 1 + )
from this follows
det ( M r ) = ( n μ ) r 1 ( n j r ) μ ( i 1 ) μ b r ( r 1 ) b ( r 1 ) ( r 2 ) b 21
that proves Theorem 2.  □
To prove that matrix T n 1 ( μ ) is positive defined, we have
T R = 1 2 ( T n 1 ( μ ) + T n 1 * ( μ ) ) .
It is obvious that the matrix T R , in addition to everything else, is also bisymmetric (symmetric with respect to both the main and secondary diagonals). Using the high-level mathematical package MATLAB, the eigenvalues of the matrix T R were considered for various values of μ and the dimension of the matrix N. It was shown that all eigenvalues of the matrix T R , for any N 3000 and μ > 0 , are positive; that is, the above calculations confirm the hypothesis that the matrix T R is positive definite. This became the basis for us to assume that the matrix T ( n 1 ) ( μ ) under study is positive definite. Naturally, the operator A ρ corresponding to the matrix T ( n 1 ) ( μ ) , will also be positive definite [7,8].
We give a strong proof of the positive definiteness of the matrix T R ( μ ) . First, let us write out the matrices T 6 ( 1 / 2 ) , T 6 * ( 1 / 2 ) , T R ( 1 / 2 ) using the MATLAB package
T 6 ( 1 / 2 ) = 2.4495 2.2361 2 1.7321 1.4142 1 0.8184 3.1623 2.8284 2.4495 2 1.4142 0.5010 1.2272 3.4641 3 2.4495 1.7321 0.3164 0.7305 1.3542 3.4641 2.8284 2 0.1857 0.4174 0.7305 1.2272 3.1623 2.2361 0.0839 0.1857 0.3164 0.5010 0.8184 2.4495
T 6 * ( 1 / 2 ) = 2.4495 0.8184 0.5010 0.3164 0.1857 0.0839 2.2361 3.1623 1.2272 0.7305 0.4174 0.1857 2 2.8284 3.4641 1.3542 0.7305 0.3164 1.7321 2.4495 3 3.4641 1.2272 0.5010 1.4142 2 2.4495 2.8284 3.1623 0.8184 1 1.4142 1.7321 2 2.2361 2.4495
T R ( 1 / 2 ) = 2.4495 1.5272 1.2505 1.0242 0.8 0.5420 1.5272 3.1623 2.0278 1.5900 1.2087 0.8 1.2505 2.0278 3.4641 2.1771 1.59 1.0242 1.0242 1.5900 2.1771 3.4641 2.0278 1.2505 0.8 1.2087 1.59 2.0278 3.1623 1.5272 0.5420 0.8 1.0242 1.5272 1.5272 2.4495
A simple analysis of these matrices shows that the elements from the main diagonal (including the diagonal itself) increase in rows and columns from the edges to the main diagonal. That is, the following statements hold:
Lemma 1.
For any fixed i 0 j , the relations
a i 0 , j a i 0 , j + 1 , i 0 j ;
a i 0 , j < a i 0 , j + 1 , i 0 > j .
Proof. 
We write the formula for the general element of the matrix
a i j = ( N i i j ) μ θ ( i , j ) ( N i N j ) μ ,
where
θ ( i , j ) = 0 , j i 1 , j < i .
Obviously, the elements under the main diagonal are calculated as follows
a i j = ( N i i j ) μ ( N i N j ) μ , i > j ,
and the elements under the main diagonal are
a i j = ( N i i j ) μ , i > j .
From these formulas, it follows that the elements located above the main diagonal decrease. To consider the elements under the main diagonal, we introduce the generating function
φ ( x ) = ( N i i x ) μ θ ( i , j ) ( N i N x ) μ , μ ( 0 , 1 ) , x [ 1 , N ] .
Obviously, the derivative of this function is positive on the segment x [ 1 , N ] , which means that the function φ ( x ) increases on the segment x [ 1 , N ] . That is the prove.    □
Lemma 2.
For any fixed j 0 < i , the relations
a i , j 0 a i + 1 , j 0 , i j 0 ;
a i , j 0 < a i + 1 , j 0 , i < j 0
hold.
Proof. 
The proof of Lemma 2 is similar to the proof of Lemma 1.    □
Lemma 3.
The statements of Lemmas 1 and 2 are valid for the matrices T n T ( μ ) ( T n T ( μ ) is the transposed matrix).
Lemma 4.
The statements of Lemmas 1,2,3 are also valid for the matrices T R ( μ ) = T n ( μ ) + T n T ( μ ) 2 .
Using these lemmas, we prove the following theorem
Theorem 3.
The matrix T R ( μ ) = T n ( μ ) + T n T ( μ ) 2 is a positive defined for μ [ 0 , 1 ] .
Proof. 
It is obvious that all the main lead minors of the matrix T R ( 0 ) are non-negative. In the same way, all main lead minors of the matrix T R ( 1 ) are positive.
Let us show for μ ( 0 , 1 ) that all main lead minors of the matrix are T R ( μ ) 0 . To do this, it is enough to prove that all the rows (columns) of the leading main lead minors of the matrix are linearly independent. In proving this statement, without loss of generality, for definiteness, we consider rows with numbers k and k + 1 . Then, it suffices to note that, by Lemma 4, a k , 1 a k + 1 , 1 < 1 and a k , k + 1 a k , k > 1 , which proves the linear independence of these rows.
Let us introduce the following function
d e t T R ( μ ) = Δ ( μ ) , μ [ 0 , 1 ] .
It is known that Δ ( 0 ) 0 and Δ ( 1 ) 0 .
From the Theorem 3 it follows that the operator A ρ is positive definite for 1 / 2 < ρ < .
The proof of the positive definiteness of the operator A ρ for 1 / 2 < ρ < can also be carried out as follows. Let us define
1 Γ ( 1 ρ ) 0 x ( x t ) 1 ρ 1 u ( t ) d t 0 1 x 1 ρ 1 ( 1 t ) 1 ρ 1 u ( t ) d t = v ( x ) .
We act on both sides of this equation by the operator D 1 ρ , where D 1 ρ is the fractional differentiation operator in the Riemann–Liouville sense, then u ( x ) = D 1 ρ v [5]. Then we may show that the form ( A ρ u , u ) > 0 .
In reality,
( A ρ u , u ) = ( v , D 1 ρ v ) = 1 Γ ( 1 ρ ) d d x 0 x f ( t ) ( x t ) 1 ρ 1 d t , f ( x ) = 1 Γ ( 1 ρ ) d d x 0 x f ( t ) ( x t ) 1 ρ 1 d t , f ( x ) = ( J 1 ρ f , f ) ,
where
( J 1 ρ f ) ( t ) = 1 Γ ( 1 ρ ) 0 t ( t s ) 1 1 ρ f ( s ) d s
-is the operator of fractional integration in the Riemann–Liouville sense of order 1 ρ .    □
Since [1] ( J 1 ρ f , f ) > 0 , for 1 ρ < 1 , and taking into account that operator A ρ is kernel [5,9], we prove the following theorem
Theorem 4.
The system of eigenfunctions and associated functions A ρ for 1 2 < 1 / ρ < 1 is complete in L 2 ( 0 , 1 ) .
Corollary 1.
Since the operator A ρ is positive definite, then all matrices T ( n 1 ) ( μ ) for n > N are positive definite.
Remark 1.
The matrix T ( n 1 ) ( μ ) may be presented as
T ( n 1 ) ( μ ) = B * B
where B is the triangular matrix.

3. Discussion

Operators generated (induced) by a differential expression of a fractional order and boundary conditions of the Sturm–Liouville type are non-self-adjoint and their spectral structure is almost not studied. The methods proposed by the authors are fundamentally new. They allow study of the completeness of systems of eigenfunctions and associated functions of these operators.

4. Conclusions

Thus, our spectral analysis of the operators generated by boundary value problems for fractional differential equations and boundary conditions of the Sturm–Liouville type, using matrix calculus, shows that the spectral structure of these operators can be studied by the matrices we studied above.

Author Contributions

All the main results of this paper were provided by M.A. Discussion and editing were provided by T.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Aleroev, M.T.; Aleroev, T.S.; Kirane, M.; Tang, Y. On one class of persymmetric matrices generated by boundary value problems for differential equations of fractional order. Appl. Math. Comput. 2015, 268, 151–163. [Google Scholar] [CrossRef]
  2. Aleroev, M.; Aleroeva, H.; Aleroev, T. Proof of the completeness of the system of eigenfunctions for one boundary-value problem for the fractional differential equation. AIMS Math. 2019, 4, 714–720. [Google Scholar] [CrossRef]
  3. Aleroev, T.S.; Aleroeva, H.T. On a class of non-selfadjoint operators, corresponding to differential equations of fractional order. Izv. Vyss. Uchebnykh Zaved. Mat. 2014, 10, 3–12. [Google Scholar]
  4. Livshits, M.S. On the spectral decomposition of linear non-self-adjoint operators. Mat. Sb. 1954, 34, 145–199. [Google Scholar]
  5. Aleroev, T.S.; Aleroeva, H.T. Problems of Sturm-Liouville Type for Differential Equations with Fractional Derivatives; Kochubei, A., Luchko, Y., Eds.; Fractional Differential Equation; Walter de Gruyter GmbH: Berlin, Germany, 2019; pp. 21–46. [Google Scholar]
  6. Bellman, R. Introduction to Matrix Analysis; SIAM: Philadelphia, PA, USA, 1997; 431p. [Google Scholar]
  7. Gantmacher, F.R. The Theory of Matrices; Chelsea Pub: New York, NY, USA, 1959; Volumes 1–2. [Google Scholar]
  8. Faddeev, D.K.; Faddeeva, V.N. Computational Methods of Linear Algebra, 2nd ed.; Fizmatgiz: Moscow, Russia, 1963. (In Russian) [Google Scholar]
  9. Aleroev, T.S.; Aleroeva, H.T. On one class of nonself-adjoint operators associated with differential equations of fractional order. Ukr. Math. Bull. 2015, 12, 293–310. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Aleroev, M.; Aleroev, T. On Positive Definite Kernels of Integral Operators Corresponding to the Boundary Value Problems for Fractional Differential Equations. Entropy 2022, 24, 515. https://doi.org/10.3390/e24040515

AMA Style

Aleroev M, Aleroev T. On Positive Definite Kernels of Integral Operators Corresponding to the Boundary Value Problems for Fractional Differential Equations. Entropy. 2022; 24(4):515. https://doi.org/10.3390/e24040515

Chicago/Turabian Style

Aleroev, Mukhamed, and Temirkhan Aleroev. 2022. "On Positive Definite Kernels of Integral Operators Corresponding to the Boundary Value Problems for Fractional Differential Equations" Entropy 24, no. 4: 515. https://doi.org/10.3390/e24040515

APA Style

Aleroev, M., & Aleroev, T. (2022). On Positive Definite Kernels of Integral Operators Corresponding to the Boundary Value Problems for Fractional Differential Equations. Entropy, 24(4), 515. https://doi.org/10.3390/e24040515

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