Next Article in Journal
Asymmetric Fractal Characteristics and Market Efficiency Analysis of Style Stock Indices
Next Article in Special Issue
Axisymmetric Fractional Diffusion with Mass Absorption in a Circle under Time-Harmonic Impact
Previous Article in Journal
Domain Adaptation with Data Uncertainty Measure Based on Evidence Theory
Previous Article in Special Issue
Heterogeneous Diffusion and Nonlinear Advection in a One-Dimensional Fisher-KPP Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On an Approximate Solution of the Cauchy Problem for Systems of Equations of Elliptic Type of the First Order

by
Davron Aslonqulovich Juraev
1,2,†,
Ali Shokri
3,† and
Daniela Marian
4,*,†
1
Department of Natural Science Disciplines, Higher Military Aviation School of the Republic of Uzbekistan, Karshi 180117, Uzbekistan
2
Department of Mathematics, Anand International College of Engineering, Jaipur 303012, Rajasthan, India
3
Department of Mathematics, Faculty of Basic Sciences, University of Maragheh, Maragheh 83111-55181, Iran
4
Department of Mathematics, Technical University of Cluj-Napoca, 28 Memorandumului Street, 400114 Cluj-Napoca, Romania
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Entropy 2022, 24(7), 968; https://doi.org/10.3390/e24070968
Submission received: 24 May 2022 / Revised: 8 July 2022 / Accepted: 9 July 2022 / Published: 13 July 2022
(This article belongs to the Special Issue Dynamical Systems, Differential Equations and Applications)

Abstract

:
In this paper, on the basis of the Carleman matrix, we explicitly construct a regularized solution of the Cauchy problem for the matrix factorization of Helmholtz’s equation in an unbounded two-dimensional domain. The focus of this paper is on regularization formulas for solutions to the Cauchy problem. The question of the existence of a solution to the problem is not considered—it is assumed a priori. At the same time, it should be noted that any regularization formula leads to an approximate solution of the Cauchy problem for all data, even if there is no solution in the usual classical sense. Moreover, for explicit regularization formulas, one can indicate in what sense the approximate solution turns out to be optimal.

1. Introduction

Most of actively developing modern area of scientific knowledge is the theory of correctly and incorrectly posed problems, most of which have practical value and require decision making in uncertain or contradictory conditions. The development and justification of methods for solving such a complex problems as ill-posed ones is intensely investigated of the present time. The results regarding ill-posed problems are a scientific research apparatus for many scientific areas, such as differentiation of approximately given functions, solving inverse boundary value problems, solving problems of linear programming and control systems, solving systems of linear equations, degenerate or ill-conditioned, etc.
The concept of a “well-posed problem” was first introduced by the French mathematician J. Hadamard in 1923 when he considered for partial differential equations of mathematical physics the extension of boundary value problems. The concept of correctness of problems was the basis for the classification of boundary value problems. In this case, the correctness of the problem statement was ensured by the fulfillment of two conditions: the existence of a solution and its uniqueness. The requirement of stability of the solution was subsequently attached to the first two by other mathematicians already during a more in-depth study of this class of problems. Problems in which any of the three conditions for the correct formulation of the problem (stability, existence or uniqueness) is not fulfilled belong to the class of ill-posed problems. The need to solve unstable problems like the one above requires a more precise definition of the solution to the problem (example Hadamard, see, for instance [1], p. 39).
We will say that the problem is correctly posed according to Tikhonov (See [2]) if:
(1)
the solution of the problem exists in some class;
(2)
the solution is unique in this class;
(3)
the solution of the problem depends continuously on the input data.
The Cauchy problem for systems of elliptic equations with constant coefficients belongs to the family of ill-posed problems: the solution of the problem is unique, but unstable. For more details on this subject can be consulted [2,3,4,5,6,7,8,9,10]. The paper studies the construction of exact and approximate solutions to the ill-posed Cauchy problem for matrix factorizations of the Helmholtz equation. Such problems naturally arise in mathematical physics and in various fields of natural science (for example, in electro-geological exploration, in cardiology, in electrodynamics, etc.). In general, the theory of ill-posed problems for elliptic systems of equations has been sufficiently developed thanks to the works of A.N. Tikhonov, V.K. Ivanov, M.M. Lavrent’ev, N.N. Tarkhanov and others famous mathematicians. Among them, the most important for applications are the so-called conditionally well-posed problems, characterized by stability in the presence of additional information about the nature of the problem data. One of the most effective ways to study such problems is to construct regularizing operators. For example, this can be the Carleman-type formulas (as in complex analysis) or iterative processes (the Kozlov-Maz’ya-Fomin algorithm, etc.) [10]. Boundary problems, as well as numerical solutions of some problems, are considered in works [11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32].
We construct, in this paper, an explicit Carleman matrix, regarding the Cauchy problem for Helmholtz’s equation, based on works [7,8,9,10]. Using this, a regularized solution of the Cauchy problem for the matrix factorization of the Helmholtz equation is given. Some formulas of Carleman type for certain equations and systems of elliptic type are given in [7,8,9,10,33,34,35,36,37,38,39]. In work [33] it was considered the Cauchy problem for the Helmholtz equation in an arbitrary bounded plane domain with Cauchy data, known only on the region boundary. In [40], the Cauchy problem for the Helmholtz equation in a bounded domain was considered. In the present study, we have constructed an approximate solution of the Cauchy problem for matrix factorizations of the Helmholtz equation in a two-dimensional unbounded domain.
In many well-posed problems it is not easy to compute the values of the function on the whole boundary. Thus, one of the important problems in the theory of differential equations is the reconstructing of the solution of systems of equations of first order elliptic type, factorizing the Helmholtz operator (see, for instance [34,35,36,37,38,39]).
The Cauchy problem for elliptic equations was investigated in [6,7,40] and subsequently developed in [9,10,33,35,36,37,38,39].
Next we establish the notations used in the paper.
Let x = ( x 1 , x 2 ) R 2 , y = ( y 1 , y 2 ) R 2 . We consider in R 2 an unbounded domain, simply-connected, Ω R 2 . We suppose that its border Ω is piece wise smooth and is composed of the plane T: y 2 = 0 and a smooth curve Σ lying in the half-space y 2 > 0 , that is, Ω = Σ T .
Let:
r = y x , α = y 1 x 1 , z = i a 2 + α 2 + y 2 , a 0 ,
x = x 1 , x 2 T , x ξ T , ξ T = ξ 1 ξ 2 transposed vector ξ ,
V ( x ) = ( V 1 ( x ) , , V n ( x ) ) T , v 0 = ( 1 , , 1 ) R n , n = 2 m , m = 2 ,
E ( w ) = w 1 0 0 0 w 2 0 0 0 0 w n diagonal matrix , w = ( w 1 , , w n ) R n .
We consider a ( n × n ) —dimensional matrix D ( ξ T ) such that
D * ( ξ T ) D ( ξ T ) = E ( ( ξ 2 + λ 2 ) v 0 ) ,
where D * ( ξ T ) is the Hermitian conjugate matrix D ( ξ T , ) λ R and the elements of D ( ξ T ) are linear functions with constant coefficients of the complex plane.
We also consider the system of differential equations:
D x V ( x ) = 0 , x Ω ,
D x being the matrix of first-order differential operators.
Let A Ω = { V : Ω ¯ R n V is continuous on Ω ¯ = Ω Ω and V satisfies the system (1)}.

2. Statement of the Cauchy Problem

Let f C ( Σ , R n ) . We formulate the following Cauchy problem for the system (1):
Let V ( y ) A ( Ω ) such that
V ( y ) Σ = f ( y ) , y Σ .
We specify that V ( y ) is defined on Ω , knowing f ( y ) , y Σ .
If V ( y ) A ( Ω ) , then
V ( x ) = Ω L ( y , x ; λ ) V ( y ) d s y , x Ω ,
L ( y , x ; λ ) = E φ 2 ( λ r ) v 0 D * x D ( t T ) ,
where t = ( t 1 , t 2 ) means the unit exterior normal at a point y Ω and φ 2 ( λ r ) represents the fundamental solution of the Helmholtz equation in R 2 , that is
φ 2 ( λ r ) = i 4 H 0 ( 1 ) ( λ r ) ,
H 0 ( 1 ) ( λ r ) being the the Hankel function of the first kind [41].
An entire function K ( z ) is introduced, taking real values for real part of z, ( z = a + i b , a , b R ) and such that:
K ( z ) 0 , sup b 1 b p K ( p ) ( z ) = B ( a , p ) < , < a < , p { 0 , 1 , 2 } .
Let
Ψ ( y , x ; λ ) = 1 2 π K ( x 2 ) 0 Im K ( z ) z x 2 a I 0 ( λ a ) a 2 + α 2 d a , for y x ,
where I 0 ( λ a ) = J 0 ( i λ a ) is the zero order Bessel function of the first kind [4].
We remark that (3) holds if we consider
Ψ ( y , x ; λ ) = φ 2 ( λ r ) + g ( y , x ; λ ) ,
instead φ 2 ( λ r ) , g ( y , x ) being the regular solution of the Helmholtz equation with respect to the variable y, including the case y = x .
Hence (3) becomes:
V ( x ) = Ω L ( y , x ; λ ) V ( y ) d s y , x Ω ,
L ( y , x ; λ ) = E Ψ ( y , x ; λ ) v 0 D * x D ( t T ) .
Formula (8) can be generalized for the case when Ω is unbounded.
Suppose Ω lies inside a strip of the smallest width defined by:
0 < y 2 < h , h = π ρ , ρ > 0 ,
and Ω extends to infinity.
So next we consider an unbounded domain Ω R 2 finitely connected, having a piecewise smooth boundary Ω ( Ω —extends to infinity).
Let Ω R be the part of Ω situated inside a circle centered at zero, having radius R:
Ω R = y : y Ω , y < R , Ω R = Ω \ Ω R , R > 0 .
Theorem 1.
Consider V ( y ) A ( Ω ) . If x Ω , x fixed, we have
lim R Ω R L ( y , x ; λ ) V ( y ) d s y = 0 ,
then the Formula (8) is true.
Proof. 
For x Ω ( x < R ) , x fixed, using (8) into account, we get
Ω L ( y , x ; λ ) V ( y ) d s y = Ω R L ( y , x ; λ ) V ( y ) d s y + Ω R L ( y , x ; λ ) V ( y ) d s y = V ( x ) + Ω R L ( y , x ; λ ) V ( y ) d s y , x Ω R .
Taking into account condition (9), for R , we obtain (8).
We also suppose
Ω exp d 0 ρ 0 y 1 d s y < , 0 < ρ 0 < ρ ,
for some d 0 > 0 , and
V ( y ) exp exp ρ 2 y 1 , ρ 2 < ρ , y Ω .
In (6) we put
K ( z ) = exp d i ρ 1 z h 2 d 1 i ρ 0 d h 2 , K ( x 2 ) = exp d cos ρ 1 x 2 h 2 + d 1 cos i ρ 0 x 2 h 2 , 0 < ρ 1 < ρ , 0 < x 2 < h ,
where
d = 2 c e ρ 1 x 1 , d 1 > d 0 cos ρ 0 h 2 , c 0 , d > 0 .
Hence (8) holds.
Consider x Ω be fixed and y . In the following we estimate the function Ψ ( y , x ; λ ) and also its derivatives Ψ ( y , x ; λ ) y j , j { 1 , 2 } . For the estimation Ψ ( y , x ; λ ) y j we use equalities
2 π K ( x 2 ) Ψ ( y , x ; λ ) y 1 = ( y 1 x 1 ) Re K ( z 0 ) sign ( y 1 x 1 ) ( y 2 x 2 ) Im K ( z 0 ) r 2 ( y 1 x 1 ) λ 0 a 2 + α 2 Re K ( w ) ( y 2 x 2 ) Im K ( z ) a 2 + r 2 · I 1 ( λ a ) d a a 2 + α 2 , y x , z 0 = i y 1 x 1 + y 2 , I 1 ( λ a ) = I 0 ( λ a )
and
2 π K ( x 2 ) Ψ ( y , x ; λ ) y 2 = ( y 2 x 2 ) Re K ( z 0 ) ( y 1 x 1 ) Im K ( z 0 ) r 2 λ 0 ( y 2 x 2 ) Re K ( z ) a 2 + α 2 Im K ( z ) a 2 + r 2 I 1 ( λ a ) d u , y 1 x 1 ,
which are obtained from (6).
Really,
exp d i ρ 1 z h 2 d 1 i ρ 0 z h 2 = exp Re d i ρ 1 z h 2 d 1 i ρ 0 z h 2 = exp d ρ 1 a 2 + α 2 cos ρ 1 y 2 h 2 d 1 ρ 0 a 2 + α 2 cos ρ 0 y 2 h 2 .
As
π 2 ρ 1 ρ · π 2 ρ 1 ρ · π 2 < π 2 , π 2 ρ 1 ρ · π 2 ρ 0 y 2 h 2 ρ 1 ρ · π 2 < π 2 .
Consequently,
cos ρ y 2 h 2 > 0 , cos ρ 0 y 2 h 2 cos h ρ 0 2 > δ 0 > 0 ,
It does not vanish in the region Ω and
Ψ ( y , x ; λ ) = O exp ε ρ 1 y 1 , ε > 0 , y , y Ω Ω , Ψ ( y , x ; λ ) y j = O exp ε ρ 1 y 1 , ε > 0 , y , y Ω Ω , j { 1 , 2 } .
We now choose ρ 1 with the condition ρ 2 < ρ 1 < ρ . Then condition (10) is fulfilled and the integral formula (8) is true. Theorem 1 is proved. □
Condition (12) may be weakened.
Consider
A ρ ( Ω ) = V ( y ) A ( Ω ) , V ( y ) exp O exp ρ y 1 , y , y Ω .
Theorem 2.
If V ( y ) A ρ ( Ω ) satisfies the growth condition
V ( y ) C exp c cos ρ 1 y 2 h 2 exp ρ 1 y 1 , c 0 , 0 < ρ 1 < ρ , y Ω ,
C constant, then (8) is true.
Proof. 
Divide Ω by a line y 2 = h 2 into the following two domains:
Ω 1 = y : 0 < y 2 < h 2   and   Ω 2 = y : h 2 < y 2 < h
.
Consider first the domain Ω 1 . We put K 1 ( z ) in (6),
K 1 ( z ) = K ( z ) exp δ i τ z h 2 δ 1 i ρ z h 2 , ρ < τ < 2 ρ , δ > 0 , δ 1 > o ,
K ( z ) being given by (12). With this notation, (10) is true.
Really,
exp i τ z h 4 δ 1 i ρ z h 4 = exp δ τ a 2 + α 2 cos τ y 2 h 4 = exp δ τ a 2 + α 2 exp δ exp τ y 1 ,
π 2 τ π 4 τ y 2 h 4 τ π 2 < h 2   and   cos τ y 2 h 4 cos τ h 4 δ 0 > 0
.
Let denote by Ψ + ( y , x ; λ ) the corresponding function Ψ ( y , x ; λ ) . As
cos τ y 2 h 4 δ 0 , y Ω 1 Ω 1 ,
then for fixed x Ω 1 , y Ω 1 Ω 1 ,
Ψ + ( y , x ; λ ) = O exp ( δ 0 exp τ y 1 , y , ρ < τ < 2 ρ , Ψ + ( y , x ; λ ) y j = O exp ( δ 0 exp τ y 1 , y , ρ < τ < 2 ρ , j { 1 , 2 } .
Suppose that V ( y ) A ρ ( Ω 1 ) satisfies:
V ( y ) C exp exp 2 ρ ε y 1 , ε > 0 , y Ω 1 .
We consider τ such that 2 ρ ε < τ < 2 ρ in (17).
Hence (17) is satisfied for the region Ω 1 , so
V ( x ) = Ω 1 L ( y , x ; λ ) V ( y ) d s y , x Ω 1 ,
L ( y , x ; λ ) = E Ψ + ( y , x ; λ ) v 0 D * x D ( t T ) .
If V ( y ) A ρ ( Ω 2 ) satisfies the growth condition (16) in Ω 2 , and 2 ρ ε < τ < 2 ρ , then
V ( x ) = Ω 2 L ( y , x ; λ ) V ( y ) d s y , x Ω 2 ,
L ( y , x ; λ ) = E Ψ ( y , x ; λ ) v 0 D * x D ( t T ) .
Here Ψ ( y , x ; λ ) it is defined by the formula (6), in which K ( z ) it is replaced by the function K 2 ( z ) :
K 2 ( z ) = K ( z ) exp δ i τ z h 1 δ 1 i ρ z h 2
where
h 1 = h 2 + h 4 , h 2 < y 2 < h , h 2 < x 2 < h 1 , δ > 0 , δ 1 > 0 .
In the formulas obtained with this formula, the integrals (according to (11)) converge uniformly for δ 0 , when V ( y ) A ρ ( Ω ) . In these formulas we put δ = 0 , hence
V ( x ) = Ω L ( y , x ; λ ) V ( y ) d s y , x Ω , x 2 h 2 ,
L ( y , x ; λ ) = E Ψ ˜ ( y , x ; λ ) v 0 D * x D ( t T ) .
(integrals over the cross section y 2 = h 2 are mutually destroyed)
Ψ ˜ ( y , x ; λ ) = ( Ψ + ( y , x ; λ ) ) δ = 0 = ( Ψ ( y , x ; λ ) ) δ = 0 .
Ψ ˜ ( y , x ; λ ) is obtained here by (6), K ( z ) being given by (17), where δ = 0 is considered. Using now the continuation principle, (22) holds, x Ω . Under condition (18) and (22) holds, δ 1 0 . Considering δ 1 = 0 , Theorem 2 is proved. □
Choosing
K ( z ) = 1 z x 2 + 2 h exp ( σ z ) , K ( x 2 ) = 1 2 h exp ( σ x 2 ) , 0 < x 2 < h , h = π ρ ,
in (6), we get
Φ σ ( y , x ) = e σ x 2 π ( h ) 1 0 Im exp ( σ z ) ( z x 2 + 2 h ) ( z x 2 ) a I 0 ( λ a ) a 2 + α 2 d a .
Hence (8) becomes:
V ( x ) = Ω L σ ( y , x ; λ ) V ( y ) d s y , x Ω ,
L σ ( y , x ; λ ) = E Ψ σ ( y , x ; λ ) v 0 D * x D ( t T ) .

3. Regularized Solution of the Cauchy Problem

Theorem 3.
Let V ( y ) A ρ ( Ω ) satisfying
V ( y ) M , y T .
If
V σ ( x ) = Σ L σ ( y , x ; λ ) V ( y ) d s y , x Ω ,
then:
V ( x ) V σ ( x ) K ρ ( λ , x ) σ M e σ x 2 , x Ω ,
V ( x ) x j V σ ( x ) x j K ρ ( λ , x ) σ M e σ x 2 , σ > 1 , x Ω , j { 1 , 2 } ,
where K ρ ( λ , x ) are bounded functions on compact subsets of the domain Ω .
Proof. 
We prove first (28). Using (25) and (27), we have
V ( x ) = Σ L σ ( y , x ; λ ) U ( y ) d s y + T L σ ( y , x ; λ ) V ( y ) d s y = L σ ( x ) + T L σ ( y , x ; λ ) V ( y ) d s y , x Ω .
Using now (26), we obtain
V ( x ) V σ ( x ) T L σ ( y , x ; λ ) V ( y ) d s y T L σ ( y , x ; λ ) V ( y ) d s y M T L σ ( y , x ; λ ) d s y , x Ω .
We estimate now T Ψ σ ( y , x ; λ ) d s y and T Ψ σ ( y , x ; λ ) y j d s y , j { 1 , 2 } .
Using (24), we have
Ψ σ ( y , x ) = e σ ( y 2 x 2 ) π ( h ) 1 0 β + β 1 cos σ α 1 α 1 2 + β 1 2 α 1 2 + β 2 + α 1 2 + β 1 β α 1 2 + β 1 2 α 1 2 + β 2 sin σ α 1 α 1 a I 0 ( λ a ) d a ,
where
α 1 2 = a 2 + α 2 , β = y 2 x 2 , β 1 = y 2 x 2 + 2 h .
Given (31) and the inequality
I 0 ( λ a ) 2 λ π a ,
we have
T Ψ σ ( y , x ; λ ) d s y K ρ ( λ , x ) σ e σ x 2 , σ > 1 , x Ω .
Using now
Ψ σ ( y , x ; λ ) y j = Ψ σ ( y , x ; λ ) s s y j = 2 ( y j x j ) Ψ σ ( y , x ; λ ) s , s = α 2 , j { 1 , 2 } ,
according to (31) and (32) we get
T Ψ σ ( y , x ; λ ) y 1 d s y K ρ ( λ , x ) σ e σ x 2 , σ > 1 , x Ω ,
According to (31) and (32), we have
T Ψ σ ( y , x ; λ ) y 2 d s y K ρ ( λ , x ) σ e σ x 2 , σ > 1 , x Ω ,
Using the inequalities (33), (35), (36) and (30), we get the estimate (28).
We prove now (29). From (25) and (27) we get:
V ( x ) x j = Σ L σ ( y , x ; λ ) x j V ( y ) d s y + T L σ ( y , x ; λ ) x j V ( y ) d s y , V σ ( x ) x j = Σ L σ ( y , x ; λ ) x j V ( y ) d s y , x Ω , j { 1 , 2 } .
According to (37) and (26), we have
V ( x ) x j σ V ( x ) x j T L σ ( y , x ; λ ) x j V ( y ) d s y T L σ ( y , x ; λ ) x j V ( y ) d s y M T L σ ( y , x ; λ ) x j d s y , x Ω , j { 1 , 2 } .
We estimate now T Ψ σ ( y , x ; λ ) x 1 d s y and T Ψ σ ( y , x ; λ ) x 2 d s y on the part T of the plane y 2 = 0 .
We use
Ψ σ ( y , x ; λ ) x 1 = Ψ σ ( y , x ; λ ) s s x 1 = 2 ( y 1 x 1 ) Ψ σ ( y , x ; λ ) s , s = α 2 ,
for the estimation of the first integral.
From (31) and (32) and (39), we have
T Ψ σ ( y , x ; λ ) x 1 d s y K ρ ( λ , x ) σ e σ x 2 , σ > 1 , x Ω .
According to (31) and (32), we have
T Ψ σ ( y , x ; λ ) x 2 d s y K ρ ( λ , x ) σ e σ x 2 , σ > 1 , x Ω .
From inequalities (40), (41) and (38), we get (29). □
Corollary 1.
We have
lim σ V σ ( x ) = V ( x ) , lim σ V σ ( x ) x j = V ( x ) x j , j { 1 , 2 } , x Ω .
Let
Ω ¯ ε = ( x 1 , x 2 ) Ω , q > x 2 ε , q = max T ψ ( x 1 ) , 0 < ε < q ,
ψ ( x 1 ) being a curve and Ω ¯ ε Ω a compact set.
Corollary 2.
If x Ω ¯ ε , then the families of functions V σ ( x ) and V σ ( x ) x j converge uniformly for σ , i.e.,
V σ ( x ) V ( x ) , V σ ( x ) x j V ( x ) x j , j { 1 , 2 } .
We specify that the set E ε = Ω \ Ω ¯ ε is as a layer boundary for this problem.
Consider now the boundary of the domain Ω being composed of a hyper plane y 2 = 0 and a smooth curve Σ extending to infinity and lying in the strip
0 < y 2 < h , h = π ρ , ρ > 0 .
We consider Σ given
y 2 = ψ ( y 1 ) , < y 1 < ,
where ψ ( y 1 ) satisfies the condition
ψ ( y 1 ) P < , P = c o n s t .
We consider
q = max T ψ ( y 1 ) , l = max T 1 + ψ 2 ( y 1 ) .
Theorem 4.
If V ( y ) A ρ ( Ω ) satisfies (26), and on a smooth curve Σ satisfies
V ( y ) δ , 0 < δ < 1 ,
then
V ( x ) K ρ ( λ , x ) σ M 1 x 2 q δ x 2 q , σ > 1 , x Ω .
V ( x ) x j K ρ ( λ , x ) σ M 1 x 2 q δ x 2 q , σ > 1 , x Ω , j { 1 , 2 } .
Proof. 
We prove first (43). From (25), we obtain
V ( x ) = Σ L σ ( y , x ; λ ) V ( y ) d s y + T L σ ( y , x ; λ ) ) V ( y ) d s y , x Ω ,
and hence
V ( x ) Σ L σ ( y , x ; λ ) V ( y ) d s y + T L σ ( y , x ; λ ) V ( y ) d s y , x Ω .
From (42), we have
Σ L σ ( y , x ; λ ) V ( y ) d s y Σ L σ ( y , x ; λ ) V ( y ) d s y δ Σ L σ ( y , x ; λ ) d s y , x Ω .
We estimate now Σ Ψ σ ( y , x ; λ ) d s y , Σ Ψ σ ( y , x ; λ ) y 1 d s y and Σ Ψ σ ( y , x ; λ ) y 2 d s y on Σ .
Given equality (31) and (32), we have
Σ Ψ σ ( y , x ; λ ) d s y K ρ ( λ , x ) σ e σ ( q x 2 ) , σ > 1 , x Ω .
Using now (31), (32) and (34), we get
Σ Ψ σ ( y , x ; λ ) y 1 d s y K ρ ( λ , x ) σ e σ ( q x 2 ) , σ > 1 , x Ω .
From (31) and (32), we have
Σ Ψ σ ( y , x ; λ ) y 2 d s y K ρ ( λ , x ) σ e σ ( q x 2 ) , σ > 1 , x Ω .
From (48)–(50) and applying (49), we get
Σ L σ ( y , x ; λ ) V ( y ) d s y K ρ ( λ , x ) σ δ e σ ( q x 2 ) , σ > 1 , x Ω .
We know that
T L σ ( y , x ; λ ) V ( y ) d s y K ρ ( λ , x ) σ M e σ x 2 , σ > 1 , x Ω .
According to (51), (52) and (46), we obtain
V ( x ) K ρ ( λ , x ) σ 2 ( δ e σ q + M ) e σ x 2 , σ > 1 , x Ω .
Considering
σ = 1 q ln M δ ,
we get (43).
We prove now (44). From (25) we get:
V ( x ) x j = Σ L σ ( y , x ; λ ) x j V ( y ) d s y + T L σ ( y , x ; λ ) x j V ( y ) d s y = V σ ( x ) x j + T L σ ( y , x ; λ ) x j V ( y ) d s y , x Ω , j { 1 , 2 } ,
where
V σ ( x ) x j = Σ L σ ( y , x ; λ ) x j V ( y ) d s y .
We get
V ( x ) x j Σ L σ ( y , x ; λ ) x j V ( y ) d s y + T L σ ( y , x ; λ ) x j V ( y ) d s y V σ ( x ) x j + T L σ ( y , x ; λ ) x j V ( y ) d s y , x Ω , j { 1 , 2 } .
From (42), we have:
Σ L σ ( y , x ; λ ) x j V ( y ) d s y Σ L σ ( y , x ; λ ) x j V ( y ) d s y δ Σ L σ ( y , x ; λ ) x j d s y , x Ω , j { 1 , 2 } .
Now we deal with Σ Ψ σ ( y , x ; λ ) x 1 d s y , and Σ Ψ σ ( y , x ; λ ) x 2 d s y on Σ .
From (31), (32) and (39), we have
Σ Ψ σ ( y , x ; λ ) x 1 d s y K ρ ( λ , x ) σ e σ ( q x 2 ) , σ > 1 , x Ω ,
From (31) and (32), it follows:
Σ Ψ σ ( y , x ; λ ) x 2 d s y K ρ ( λ , x ) σ e σ ( q x 2 ) , σ > 1 , x Ω ,
From (59) and (60), bearing in mind (58), we have
Σ L σ ( y , x ; λ ) x j V ( y ) d s y K ρ ( λ , x ) σ δ e σ ( q x 2 ) , σ > 1 , x Ω , j { 1 , 2 } .
We known that
T L σ ( y , x ; λ ) x j V ( y ) d s y K ρ ( λ , x ) σ M e σ x 2 , σ > 1 , x Ω , j { 1 , 2 } .
According to (61), (62) and (57), we obtain
V ( x ) x j K ρ ( λ , x ) σ 2 ( δ e σ q + M ) e σ x 2 , σ > 1 , x Ω , j { 1 , 2 } .
Considering σ as in (54) we obtain (44). □
Assume that V ( y ) A ( Ω ) and instead of V ( y ) on Σ its continuous approximations f δ ( y ) are given, with error 0 < δ < 1 . We have
max Σ V ( y ) f δ ( y ) δ .
We put
V σ ( δ ) ( x ) = Σ N σ ( y , x ; λ ) f δ ( y ) d s y , x Ω .
Theorem 5.
If V ( y ) A ( Ω ) satisfies (26) on the plane y 2 = 0 , then
V ( x ) V σ ( δ ) ( x ) K ρ ( λ , x ) σ M 1 x 2 q δ x 2 q , σ > 1 , x Ω ,
V ( x ) x j V σ ( δ ) ( x ) x j K ρ ( λ , x ) σ M 1 x 2 q δ x 2 q , σ > 1 , x Ω . j { 1 , 2 } .
Proof. 
From (25) and (65), we get
V ( x ) V σ ( δ ) ( x ) = Ω L σ ( y , x ; λ ) L ( y ) d s y Σ L σ ( y , x ; λ ) f δ ( y ) d s y = Σ L σ ( y , x ; λ ) V ( y ) d s y + T L σ ( y , x ; λ ) V ( y ) d s y Σ L σ ( y , x ; λ ) f δ ( y ) d s y = Σ L σ ( y , x ; λ ) V ( y ) f δ ( y ) d s y + T L σ ( y , x ; λ ) L ( y ) d s y .
and
V ( x ) x j V σ ( δ ) ( x ) x j = Ω L σ ( y , x ; λ ) x j V ( y ) d s y Σ L σ ( y , x ; λ ) x j f δ ( y ) d s y = Σ L σ ( y , x ; λ ) x j V ( y ) d s y + T L σ ( y , x ; λ ) x j V ( y ) d s y Σ L σ ( y , x ; λ ) x j f δ ( y ) d s y = Σ L σ ( y , x ; λ ) x j V ( y ) f δ ( y ) d s y + T L σ ( y , x ; λ ) x j V ( y ) d s y , j { 1 , 2 } .
From (26) and (64), we obtain:
V ( x ) V σ ( δ ) ( x ) = Σ L σ ( y , x ; λ ) V ( y ) f δ ( y ) d s y + T L σ ( y , x ; λ ) V ( y ) d s y Σ L σ ( y , x ; λ ) V ( y ) f δ ( y ) d s y + T L σ ( y , x ; λ ) V ( y ) d s y δ Σ L σ ( y , x ; λ ) d s y + M T L σ ( y , x ; λ ) d s y .
and
V ( x ) x j V σ ( δ ) ( x ) x j = Σ L σ ( y , x ; λ ) x j V ( y ) f δ ( y ) d s y + T L σ ( y , x ; λ ) x j V ( y ) d s y Σ L σ ( y , x ; λ ) x j U ( y ) f δ ( y ) d s y + T L σ ( y , x ; λ ) x j V ( y ) d s y δ Σ L σ ( y , x ; λ ) x j d s y + M T L σ ( y , x ; λ ) x j d s y , j { 1 , 2 } .
Analog as in Theorems 3 and 4, we can prove that
V ( x ) V σ ( δ ) ( x ) K ρ ( λ , x ) σ 2 ( δ e σ q + M ) e σ x 2 ,
V ( x ) x j V σ ( δ ) ( x ) x j K ρ ( λ , x ) σ 2 ( δ e σ q + M ) e σ x 2 , j { 1 , 2 } .
Considering σ as in (54), we get (66) and (67). □
Corollary 3.
We have
lim δ 0 V σ ( δ ) ( x ) = V ( x ) , lim δ 0 V σ ( δ ) ( x ) x j = V ( x ) x j , j { 1 , 2 } , x Ω .
Corollary 4.
If x Ω ¯ ε , then the families of functions V σ ( δ ) ( x ) and V σ ( δ ) ( x ) x j are convergent uniformly, for δ 0 , i.e.,
V σ ( δ ) ( x ) V ( x ) , V σ ( δ ) ( x ) x j V ( x ) x j , j { 1 , 2 } .
The following example illustrates the possibility of incorrect formulation of the classical Cauchy problem for system (1).
Example 1.
Prove that the Cauchy problem for the following systems of linear partial differential equations is ill-posed:
x 1 V 1 x 2 V 2 = 0 , x 2 V 1 + x 1 V 2 = 0 , x 1 V 3 + x 2 V 4 = 0 , x 2 V 3 + x 1 V 4 = 0 .
Solutions to this system will be sought in the form
V 1 = U 1 e i ( λ x 1 + μ x 2 ) , V 2 = U 2 e i ( λ x 1 + μ x 2 ) , V 3 = U 3 e i ( λ x 1 + μ x 2 ) , V 4 = U 4 e i ( λ x 1 + μ x 2 ) .
Substituting these into the system, we obtain
λ 2 + μ 2 = 0 , U 1 = λ μ U 2 , λ 2 + μ 2 = 0 , U 3 = λ μ U 4 .
We choose the following μ = n , λ = i n . Then
V 1 n = U 1 n e n x 1 i n x 2 , V 2 n = i U 1 n e n x 1 i n x 2 , V 3 n = U 3 n e i ( λ x 1 + μ x 2 ) , V 4 n = i U 3 n e n x 1 i n x 2 . .
Separating the real part, we find the solutions
V 1 n = U 1 n e n x 1 cos n x 2 , V 2 n = U 1 n e n x 1 sin n x 2 , V 3 n = U 3 n e n x 1 cos n x 2 , V 4 n = U 3 n e n x 1 sin n x 2 .
The constants U 1 n and U 3 n are given by the formula U 1 n = U 3 n = e n .
Hence
V 1 n = e n e n x 1 cos n x 2 , V 2 n = e n e n x 1 sin n x 2 , V 3 n = e n e n x 1 cos n x 2 , V 4 n = e n e n x 1 sin n x 2 .
The solutions ( V 1 n , V 2 n ) , ( V 3 n , V 4 n ) satisfy at x 1 = 0 the following initial data:
V 1 n ( 0 , x 2 ) = φ 1 n ( x ) = e n cos n x 2 , V 2 n ( 0 , x 2 ) = φ 2 n ( x ) = e n sin n x 2 , V 3 n ( 0 , x 2 ) = φ 3 n ( x ) = e n cos n x 2 , V 4 n ( 0 , x 2 ) = φ 4 n ( x ) = e n sin n x 2 .
At n , these initial data tend to zero. Moreover, their derivatives φ 1 n ( k ) ( x ) , φ 2 n ( k ) ( x ) , φ 3 n ( k ) ( x ) , φ 4 n ( k ) ( x ) of orders k = 1 , 2 , , p tend to zero as n (here, p is an arbitrary fixed natural number). Indeed,
φ 1 n ( x ) = ± n k e n cos n x 2 φ 2 n ( x ) = ± n k e n sin n x 2 , if k is even,
φ 1 n ( x ) = ± n k e n sin n x 2 φ 2 n ( x ) = ± n k e n cos n x 2 , if k is odd,
φ 3 n ( x ) = ± n k e n cos n x 2 φ 4 n ( x ) = ± n k e n sin n x 2 , if k is even,
φ 3 n ( x ) = ± n k e n sin n x 2 φ 4 n ( x ) = ± n k e n cos n x 2 , if k is odd.
On the other hand, V 1 n ( x 1 , x 2 ) , V 2 n ( x 1 , x 2 ) , V 3 n ( x 1 , x 2 ) , V 4 n ( x 1 , x 2 ) is unbounded for any x 1 .
We see that no matter what norm we choose to estimate the value of the initial data, we will not be able to assert that the smallness of this norm implies the smallness of the solution (the solution is estimated here by the maximum of its modulus). As admissible norms for the initial data, we here admit the following norms:
φ 1 ( x ) p = max 0 k p sup x 2 φ 1 ( k ) ( x ) , φ 2 ( x ) p = max 0 k p sup x 2 φ 2 ( k ) ( x ) , φ 3 ( x ) p = max 0 k p sup x 2 φ 3 ( k ) ( x ) , φ 4 ( x ) p = max 0 k p sup x 2 φ 4 ( k ) ( x ) .
That is, there is no continuous dependence on the initial data and, therefore, the problem is set incorrectly. Thus, this problem does not have stability properties and, therefore, is ill-posed. We have seen that the solution of the Cauchy problem for this system is unstable. If we narrow the class of solutions under consideration to a compact set, then the problem becomes conditionally well-posed. To estimate the conditional stability, we can apply the results of the above theorems.
Example 2.
Let a system of partial differential equations of first order of the form
V 1 x 1 V 2 x 2 + i V 4 = 0 , V 1 x 2 + V 2 x 1 + i V 3 = 0 , V 3 x 1 + V 4 x 2 i V 2 = 0 , V 3 x 2 + V 4 x 1 + i V 1 = 0 .
Check that the following relation holds:
D * ( ξ T ) D ( ξ T ) = E ( ( ξ 2 + λ 2 ) v 0 ) , v 0 = ( 1 , , 1 ) R n .
Assuming x 1 ξ 1 and x 2 ξ 2 , compose the following matrices
D ( ξ T ) = ξ 1 ξ 2 0 i ξ 2 ξ 1 i 0 0 i ξ 1 ξ 2 i 0 ξ 2 ξ 1 , D * ( ξ T ) = ξ 1 ξ 2 0 i ξ 2 ξ 1 i 0 0 i ξ 1 ξ 2 i 0 ξ 2 ξ 1 .
The relation (68) is easily checked.

4. Conclusions

We have explicitly determined a regularized solution of the Cauchy problem for the matrix factorization Helmholtz’s equation in an unbounded two-dimensional domain. We specify that the approximate values of V ( x ) and V ( x ) x j , x Ω , j { 1 , 2 } must be determined, for solving applicable problems.
We have built a vector-functions family V ( x , f δ ) = V σ ( δ ) ( x ) and V ( x , f δ ) x j = V σ ( δ ) ( x ) x j , ( j { 1 , 2 } ) depending on σ (which is a parameter) and we have proved that for certain choices of σ = σ ( δ ) , δ 0 , and under certain conditions, the family V σ ( δ ) ( x ) and V σ ( δ ) ( x ) x j converges to V ( x ) and respectively to V ( x ) x j , x Ω . Hence, V σ ( δ ) ( x ) and V σ ( δ ) ( x ) x j determine the regularization of the solution of problems (1) and (2).

Author Contributions

Conceptualisation, D.A.J.; methodology, A.S. and D.M.; formal analysis, D.A.J., A.S. and D.M.; writing—original draft preparation, D.A.J., A.S. and D.M. All authors read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the editor and reviewers in advance for helpful comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hadamard, J. The Cauchy Problem for Linear Partial Differential Equations of Hyperbolic Type; Nauka: Moscow, Russia, 1978. [Google Scholar]
  2. Tikhonov, A.N. On the solution of ill-posed problems and the method of regularization. Dokl. Akad. Nauk. SSSR 1963, 151, 501–504. [Google Scholar]
  3. Aizenberg, L.A. Carleman’s Formulas in Complex Analysis; Nauka: Novosibirsk, Russia, 1990. [Google Scholar]
  4. Bers, A.; John, F.; Shekhter, M. Partial Differential Equations; Mir: Moscow, Russia, 1966. [Google Scholar]
  5. Carleman, T. Les Fonctions Quasi Analytiques; Gautier-Villars et Cie: Paris, France, 1926. [Google Scholar]
  6. Goluzin, G.M.; Krylov, V.M. The generalized Carleman formula and its application to the analytic continuation of functions. Sb. Math. 1933, 40, 144–149. [Google Scholar]
  7. Lavrent’ev, M.M. On the Cauchy problem for second-order linear elliptic equations. Rep. USSR Acad. Sci. 1957, 112, 195–197. [Google Scholar]
  8. Lavrent’ev, M.M. On Some Ill-Posed Problems of Mathematical Physics; Nauka: Novosibirsk, Russia, 1962. [Google Scholar]
  9. Tarkhanov, N.N. A criterion for the solvability of the ill-posed Cauchy problem for elliptic systems. Dokl. Math. 1990, 40, 341–345. [Google Scholar]
  10. Tarkhanov, N.N. The Cauchy Problem for Solutions of Elliptic Equations; Akademie-Verlag: Berlin, Germany, 1995; Volume 7. [Google Scholar]
  11. Pankov, P.S.; Zheentaeva, Z.K.; Shirinov, T. Asymptotic reduction of solution space dimension for dynamical systems. TWMS J. Pure Appl. Math. 2021, 12, 243–253. [Google Scholar]
  12. Berdawood, K.; Nachaoui, A.; Saeed, R.; Nachaoui, M.; Aboud, F. An efficient D-N alternating algorithm for solving an inverse problem for Helmholtz equation. Discret. Contin. Dyn. Syst.-S 2022, 15, 57–78. [Google Scholar] [CrossRef]
  13. Bulnes, J. An unusual quantum entanglement consistent with Schrödinger’s equation. Glob. Stoch. Anal. 2022, 9, 78–87. [Google Scholar]
  14. Bulnes, J. Solving the heat equation by solving an integro-differential equation. Glob. Stoch. Anal. 2022, 9, 89–97. [Google Scholar]
  15. Corcino, B.C.; Corcino, R.B.; Damgo, B.A.A.; Cañete, J.A.A. Integral representation and explicit formula at rational arguments for Apostol–Tangent polynomials. Symmetry 2022, 14, 35. [Google Scholar] [CrossRef]
  16. Giang, N.H.; Nguyen, T.-T.; Tay, C.C.; Phuong, L.A.; Dang, T.-T. Towards predictive Vietnamese human resource migration by machine learning: A case study in northeast Asian countries. Axioms 2022, 11, 151. [Google Scholar] [CrossRef]
  17. Fayziyev, Y.; Buvaev, Q.; Juraev, D.A.; Nuralieva, N.; Sadullaeva, S. The inverse problem for determining the source function in the equation with the Riemann-Liouville fractional derivative. Glob. Stoch. Anal. 2022, 9, 43–52. [Google Scholar]
  18. Ibrahimov, V.R.; Mehdiyeva, G.Y.; Yue, X.G.; Kaabar, M.K.A.; Noeiaghdam, S.; Juraev, D.A. Novel symmetric numerical methods for solving symmetric mathematical problems. Int. J. Circuits Syst. Signal Process. 2021, 15, 1545–1557. [Google Scholar] [CrossRef]
  19. Ramazanova, A.T. Necessary conditions for the existence of a saddle point in one optimal control problem for systems of hyperbolic equations. Eur. J. Pure Appl. Math. 2021, 14, 1402–1414. [Google Scholar] [CrossRef]
  20. Shokri, A.; Saadat, H. P-stability, TF and VSDPL technique in Obrechkoff methods for the numerical solution of the Schrödinger equation. Bull. Iran. Math. Soc. 2016, 42, 687–706. [Google Scholar]
  21. Shokri, A.; Tahmourasi, M. A new two-step Obrechkoff method with vanished phase-lag and some of its derivatives for the numerical solution of radial Schrödinger equation and related IVPs with oscillating solutions. Iran. J. Math. Chem. 2017, 8, 137–159. [Google Scholar] [CrossRef]
  22. Shokri, A. The Symmetric P-Stable Hybrid Obrenchkoff Methods for the numerical solution of second Order IVPS. TWMS J. Pure Appl. Math. 2012, 5, 28–35. [Google Scholar]
  23. Shokri, A. An explicit trigonometrically fitted ten-step method with phase-lag of order infinity for the numerical solution of the radial Schrödinger equation. J. Appl. Comput. Math. 2015, 14, 63–74. [Google Scholar]
  24. Shokri, A.; Shokri, A.A. The hybrid Obrechkoff BDF methods for the numerical solution of first order initial value problems. Acta Univ. Apulensis Math. Inform. 2014, 38, 23–33. [Google Scholar]
  25. Marian, D.; Ciplea, S.A.; Lungu, N. Ulam-Hyers stability of Darboux-Ionescu problem. Carpathian J. Math. 2021, 37, 211216. [Google Scholar] [CrossRef]
  26. Marian, D.; Ciplea, S.A.; Lungu, N. Hyers-Ulam Stability of Euler’s Equation in the Calculus of Variations. Mathematics 2021, 9, 3320. [Google Scholar] [CrossRef]
  27. Marian, D.; Ciplea, S.A.; Lungu, N. On the Ulam-Hyers Stability of Biharmonic Equation. U.P.B. Sci. Bull. Ser. A 2020, 8, 141–148. [Google Scholar]
  28. Marian, D. Semi-Hyers–Ulam–Rassias Stability of the Convection Partial Differential Equation via Laplace Transform. Mathematics 2021, 9, 2980. [Google Scholar] [CrossRef]
  29. Marian, D. Laplace Transform and Semi-Hyers–Ulam–Rassias Stability of Some Delay Differential Equations. Mathematics 2021, 9, 3260. [Google Scholar] [CrossRef]
  30. Musaev, H.K. The Cauchy problem for degenerate parabolic convolution equation. TWMS J. Pure Appl. Math. 2021, 12, 278–288. [Google Scholar] [CrossRef]
  31. Grzegorzewski, P.; Ladek, K.G. On some dispersion measures for fuzzy data and their robustness. TWMS J. Pure Appl. Math. 2021, 12, 88–103. [Google Scholar]
  32. Adiguzel, R.S.; Aksoy, U.; Karapinar, E.; Erhan, I.M. On the solutions of fractional differential equations via Geraghty type hybrid contractions. Appl. Comput. Math. 2021, 20, 313–333. [Google Scholar]
  33. Arbuzov, E.V.; Bukhgeim, A.L. Carleman’s formula for the system of equations of electrodynamics on the plane. Sib. Electron. Math. Rep. 2008, 5, 432–488. [Google Scholar] [CrossRef]
  34. Juraev, D.A. The Cauchy problem for matrix factorizations of the Helmholtz equation in an unbounded domain. Sib. Electron. Math. Rep. 2017, 14, 752–764. [Google Scholar] [CrossRef]
  35. Juraev, D.A. On the Cauchy problem for matrix factorizations of the Helmholtz equation in an unbounded domain in R 2 . Sib. Electron. Math. Rep. 2018, 15, 1865–1877. [Google Scholar] [CrossRef] [Green Version]
  36. Juraev, D.A. Solution of the ill-posed Cauchy problem for matrix factorizations of the Helmholtz equation on the plane. Glob. Stoch. Anal. 2021, 8, 1–17. [Google Scholar]
  37. Juraev, D.A.; Gasimov, Y.S. On the regularization Cauchy problem for matrix factorizations of the Helmholtz equation in a multidimensional bounded domain. Azerbaijan J. Math. 2022, 12, 142–161. [Google Scholar]
  38. Juraev, D.A. On the solution of the Cauchy problem for matrix factorizations of the Helmholtz equation in a multidimensional spatial domain. Glob. Stoch. Anal. 2022, 9, 1–17. [Google Scholar]
  39. Zhuraev, D.A. Cauchy problem for matrix factorizations of the Helmholtz equation. Ukr. Math. J. 2018, 69, 1583–1592. [Google Scholar] [CrossRef]
  40. Yarmukhamedov, S. On the extension of the solution of the Helmholtz equation. Rep. Russ. Acad. Sci. 1997, 357, 320–323. [Google Scholar]
  41. Kythe, P.K. Fundamental Solutions for Differential Operators and Applications; Birkhauser: Boston, MA, USA, 1996. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Juraev, D.A.; Shokri, A.; Marian, D. On an Approximate Solution of the Cauchy Problem for Systems of Equations of Elliptic Type of the First Order. Entropy 2022, 24, 968. https://doi.org/10.3390/e24070968

AMA Style

Juraev DA, Shokri A, Marian D. On an Approximate Solution of the Cauchy Problem for Systems of Equations of Elliptic Type of the First Order. Entropy. 2022; 24(7):968. https://doi.org/10.3390/e24070968

Chicago/Turabian Style

Juraev, Davron Aslonqulovich, Ali Shokri, and Daniela Marian. 2022. "On an Approximate Solution of the Cauchy Problem for Systems of Equations of Elliptic Type of the First Order" Entropy 24, no. 7: 968. https://doi.org/10.3390/e24070968

APA Style

Juraev, D. A., Shokri, A., & Marian, D. (2022). On an Approximate Solution of the Cauchy Problem for Systems of Equations of Elliptic Type of the First Order. Entropy, 24(7), 968. https://doi.org/10.3390/e24070968

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