An inertial parallel algorithm for a finite family of $G$-nonexpansive mappings applied to signal recovery

1.   Introduction

In ^[1], S. Omri and L. T. Rachdi define the Gauss-Weierstrass transform $\mathscr{W}_{\nu, t}$ associated with the Hankel transform as follows:

where $\mathscr{J}_\nu(\cdot)$ is the normalized Bessel function defined in (2.2). This integral transform, which generalizes the usual Weierstrass transform ^[2,3,4], is used to solve the heat equation problem:

where the Bessel differential operator is given by

The authors in ^[1] established practical, real inversion formulas for Hankel-type heat diffusion, building on the ideas of Saitoh, Matsuura, Fujiwara, and Yamada ^[2,3,4,5,6], and utilizing the theory of reproducing kernels ^[2,3,4].

Recently, many researchers have adapted and applied this same method to various types of Gauss-Weierstrass integral transforms associated with several kinds of differential and difference-differential operators. For instance, Soltani pioneered the exploration of $L^p$-Fourier multipliers for the Dunkl operator on the real line ^[7], extremal functions on Sobolev-Dunkl spaces ^[8], multiplier operators and extremal functions related to the dual Dunkl-Sonine operator ^[9], and extremal functions on Sturm-Liouville hypergroups ^[10]. More recently, the same authors examined Dunkl-Weinstein multiplier operators ^[11]. Additional research was conducted by Dziri and Kroumi ^[12], as well as by Ghobber and Mejjaoli ^[13]. For further work related to existing results on inverse problems, some important findings can be found in ^[14,15,16].

In this work, we consider the space-fractional diffusion equation associated with the Bessel operator, which is given by

where the parameters $\nu$ and $\alpha$ are restricted by the condition $\nu \geq -\frac{1}{2},$ $1 \leq \alpha \leq 2,$ and the space-fractional Bessel operator $(-\mathscr{B}_\nu)^{\alpha/2}$, which is defined pointwise by the principal value integral ^[17],

here, the normalization constant $c_{\nu, \alpha}$ is given by

In this work, we introduce the generalized Gauss-Weierstrass transform associated with the Bessel operator by

where $\mathscr{S}^{\alpha, \nu}(x, y, \, t)$ is the fractional heat kernel, which will be defined later. For $\alpha = 2$, this integral transform simplifies to the Gauss-Weierstrass transform defined in (1.1). Thus, it can be considered a one-parameter extension of the transform (1.1). The principal motivation for considering the generalized Gauss-Weierstrass integral transform is that for $\nu = \frac{n}{2} - 1$, it reduces to the ordinary Gauss-Weierstrass transform for radial functions on the Euclidean space $\mathbf{R}^n$. Since the Bessel operator coincides with the radial part of the Laplace operator $\Delta = \sum_{i = 1}^{m} \partial_i^2$, this transform provides a significant extension. For more details, the reader is referred to the paper ^[18]. By using the properties of the Fourier-Bessel transform $\mathscr{F}_\nu$ and its connection with the $*$-convolution product (see Section 2), we first show that the transform $\mathscr{W}_{\alpha, \nu, t}$ is a one-to-one bounded linear operator from a Sobolev space $\mathcal{H}^{s}_\nu$ into $L^2_\nu(0, \infty)$. By the same argument as the standard Gauss-Weierstrass transform, we can assume that the operator $\mathscr{W}_{\alpha, \nu, t}^{-1}$ is unbounded or that its range is not closed, which causes the ill-posed problem in solving the operator equation

Then, for stable reconstruction of $\phi$, some regularization techniques are necessary. The Tikhonov regularization techniques are widely applicable (e.g., Bakushinsky and Goncharsky ^[19], Baumeister ^[20], Tikhonov and Arsenin ^[21], Tikhonov et al. ^[22]). In our case, the Tikhonov regularization can be stated as follows: For given data $\psi \in L^2_\nu(0, \infty)$, we search for a minimizer of a functional given by

with a parameter $\gamma > 0$, which is called a regularizing parameter. We show that the above variational problem has a unique solution denoted by $\mathscr{R}_{\gamma, \psi}$ and called the regularized solution; it is also referred to as the extremal solution by Soltani ^[7]. The following theorem is the main result of the paper, which provides a real inversion of the generalized Gauss-Weierstrass transform.

Theorem 1.1. Let $s > \nu + 1$. For every $\phi \in \mathcal{H}^{s}_{\nu}$ and $\psi = \mathscr{W}_{\alpha, \nu, t}(\phi)$, we have:

Moreover, the set $\{\mathscr{R}_{\gamma, \psi}\}_{\gamma > 0}$ converges uniformly to $\phi$ as $\gamma \to 0^+$.

Our paper is organized as follows:

● Section 2 serves as an introductory section that provides an overview of fundamental concepts. Topics covered include the Fourier-Bessel transform, generalized translation, generalized convolution, fractional Bessel operator, and the space-fractional Bessel diffusion equation, setting the stage for understanding subsequent content.

● Section 3 is devoted to introducing the generalized Gauss-Weierstrass transform and establishing its principal properties.

● Section 4 states the main results of the paper and provides their proofs.

2.   Preliminaries

Before revealing our main results, it is essential to establish the groundwork by introducing key notations and collecting pertinent facts about the Bessel operator. This section serves as a primer, elucidating the significance of the Fourier-Bessel transform and the Delsarte translation, which will be pivotal for the subsequent analysis.

2.1.   Fourier-Bessel transform

The normalized Bessel function is defined as follows:

where $\Gamma(\cdot)$ is the Gamma function ^[23] and $J_\nu(\cdot)$ is the Bessel function of the first kind, see [23, (10.16.9)]. Then

The normalized Bessel function arises as the unique solution to the eigenvalue problem related to the Bessel equation. More precisely, the functions defined as $x\rightarrow \mathscr{J}_\nu(\lambda x)$ stand as the unique solution to the eigenvalue problem [23, (10.13.5)]

The function $\mathscr{J}_\nu(\cdot)$ is an entire analytic function with even symmetry. Notably, there are straightforward special cases that hold:

We introduce the following notation:

● $L^p_\nu(0, \infty)$ $(1\leq p)$ represents the Lebesgue space associated with the measure

The norm $\|\phi\|_{L^p_\nu(0, \infty)}$ is the conventional norm given by

● $S_*(\mathbb{R})$ signifies the space of even functions on $\mathbb{R}$ that are infinitely differentiable and decrease rapidly, along with all their derivatives.

For $\nu\geq -1/2$, the Fourier-Bessel transform $\mathscr{F}_\nu \phi$ of $\phi\in L^1_{\nu}(0, \infty)$ is defined as:

This integral transform can be extended to establish an isometry of $L^2_\nu(0, \infty)$. For any function $\phi$ belonging to $L^1_{\nu}(0, \infty)\cap L^2_{\nu}(0, \infty)$, the following relationships hold [24, Prop. 5.Ⅲ.2]

Furthermore, its inverse is expressed as:

Moving forward, our focus shifts to the exploration of the generalized translation operator linked to the Bessel operator. This operator is symbolized as $\tau_\nu^x$ and operates on functions belonging to $L^1_\nu(0, \, \infty)$ according to the following expression [25, §3.4.1]:

With the help of this translation operator, one defines the convolution of $\phi\in L^1_\nu(0, \infty)$ and $\psi\in L^p_\nu(0, \infty)$ for $p\in [1, \, \infty)$ as the element $f*_\nu g$ of $L^p_\nu(0, \, \infty)$ given by

The following properties are obvious.

● $\mathscr{F}_\nu (\tau^x_\nu \phi)(t) = \mathscr{J}_\nu(xt)\mathscr{F}_\nu \phi(t),$

● $\mathscr{F}_\nu(\phi*_\nu \psi)(x) = \mathscr{F}_\nu \phi(x) \mathscr{F}_\nu \psi(x).$

2.2.   Space-fractional diffusion equation

In this section, we consider the space-fractional Bessel diffusion equation ^[17,18]

where the parameters $\nu$ and $\alpha$ are restricted by the condition $\nu \geq -\frac{1}{2},$ $1 \leq \alpha \leq 2,$ and the space-fractional Bessel operator $(-\mathscr{B}_\nu)^{\alpha/2}$, which is defined pointwise by the principal value integral ^[17],

here, the normalization constant $c_{\nu, \alpha}$ is given by

Moreover, the Fourier-Bessel transform of the fractional Bessel operator is given by ^[17]

Let us denote the Fourier-Bessel transform of a function $u(x, t)$ with respect to $x$ as $\hat{u}(\xi, t)$, where $\xi\geq 0$. Applying the Fourier-Bessel transform to both sides of the equation in (1.3), we obtain:

Then

Therefore,

where

Using the following scaling rules for the Fourier-Bessel transform:

we obtain the following scaling property of the kernel $\mathscr{G}^{\alpha, \nu}(x, t)$

Consequently by introducing the similarity variable $x/t^\alpha,$ we can write

where

Particular cases of the density $\mathscr{K}^{\alpha, \nu}$ are the following ^[26]:

● The density $\mathscr{K}_{2, \nu}(x)$, where $\nu \geq -\frac{1}{2}$, corresponds to the Gaussian density kernel:

● The density $\mathscr{K}_{1, \nu}$, where $\nu \geq -\frac{1}{2}$, corresponds to the Poisson density:

More generally, for $1 < \alpha < 2$, we have [17, Proposition 4.1]:

Definition 2.1. For $\nu \geq -\frac{1}{2}$ and $0 < \alpha < 2$, the generalized heat kernel is defined as:

Lemma 2.1. The heat kernel $\mathscr{G}^{\alpha, \nu}(x, y, t)$ possesses the following properties:

$i)$ $\mathscr{G}^{\alpha, \nu}(x, y, t) = \tau^y_\nu\mathscr{G}^{\alpha, \nu}(x, t).$

$ii)$ $\mathscr{G}^{\alpha, \nu}(x, y, t) = \mathscr{G}^{\alpha, \nu}(y, x, t).$

$iii)$ $\mathscr{G}^{\alpha, \nu}(x, y, t) > 0.$

$iv)$ $\| \mathscr{G}^{\alpha, \nu}(x, y, t)\|_{L^1_\nu(0, \, \infty)} = 1.$

$v)$ $\mathscr{G}^{\alpha, \nu}(x, y, t) = t^{-2(\nu+1)/\alpha} \mathscr{G}^{\alpha, \nu}(xt^{-1/\alpha}, yt^{-1/\alpha}, 1).$

$vi)$ $\mathscr{F}_\nu(\mathscr{G}^{\alpha, \nu}(\cdot, y, t))(\xi) = e^{-t\xi^\alpha} \mathscr{J}_\nu(y\xi).$

Proof. The proofs for properties $i)$, $ii)$, $v)$, and $vi)$ are straightforward. Property $iii)$ follows from [18, Theorem 7] and the positivity of the generalized translation operator $\tau^y_\nu$. The proof for property $iv)$ is derived by setting $\xi = 0$ in the formula from property $vi$. □

3.   The generalized Weierstrass transform

Definition 3.1. The generalized Weierstrass transform associated with the fractional Bessel operator, denoted as $\mathscr{W}_{\alpha, \nu, t}$, is defined on $L^2(d\sigma_\nu)$ by the following expression:

where $\mathscr{S}^{s, \nu}(x, y, \, t)$ is the generalized heat defined in (2.17).

Let $s > \nu + 1$. We define the space $\mathcal{H}^{(s)}_{\nu}$ as follows ^[1]:

This space is equipped with an inner product defined by:

and a norm:

The space $\mathcal{H}^{(s)}_{\nu}$ features a reproducing kernel, which is defined by:

Additionally, this space satisfies the following inclusions:

For further details concerning the space $\mathcal{H}^{(s)}_{\nu}$, readers are referred to the paper ^[1].

Theorem 3.1. $i)$ Let $\phi \in C_0(\mathbb{R}) \cap L^2_\nu(0, \infty)$. For $t > 0$ and $x \in [0, \infty)$, the function $\mathscr{W}_{\alpha, \mu, t}\phi(x)$ solves the following heat equation.

with the initial condition

$ii)$ The integral transform $\mathscr{W}_{\alpha, \nu, t}$, for $t > 0$, is a one-to-one bounded linear operator from $\mathcal{H}^{(s)}_{\nu}$ into $L^2_\nu(0, \infty)$, and we have:

Proof. The claim $i)$ follows from [17, Theorem 4.5]. From Lemma 2.1, for all $\phi \in L^2_\nu(0, \infty)$, we have:

This inequality shows that the transform $\mathscr{W}_{\alpha, \nu, t}$ is indeed bounded. To complete the proof of $ii),$ it remains to show that this transform is one-to-one. Let $\phi \in \mathcal{H}^{(s)}_\nu$ such that $\mathscr{W}_{\alpha, \nu, t}\phi = 0$. Then

from the injectivity of the Fourier-Bessel transform, we get $\phi = 0$. This shows that $\mathscr{W}_{\alpha, \nu, t}$ is one-to-one. □

We denote by $\mathcal{H}_{\nu, \alpha, \gamma}^{(s)}$, the space $\mathcal{H}_{\nu}^{(s)}$ equipped with the inner product

and the norm

Then, we have the following main result:

Theorem 3.2. Let $\xi$, $t > 0$ and $s > \nu + 1$. Then the Hilbert space $\mathcal{H}_{\nu, \alpha, \gamma}^{(s)}$ admits the following reproducing kernel:

that is

$(i)$For all $x \in [0, \infty)$, the function $y \mapsto \mathscr{K}_{\nu, s, \alpha, \gamma}(x, y)$ belongs to $\mathcal{H}_{\nu, \alpha, \gamma}^{(s)}$.

$(ii)$ For all $\phi \in \mathcal{H}_{\nu, \alpha, \gamma}^{(s)}$ and any $y \in [0, \infty)$.

Proof. For all $x \in [0, \infty)$, consider the function

which belongs to both $L^1_\nu(0, \infty)$ and $L^2_\nu(0, \infty)$. Then, by the Plancherel theorem for the Fourier-Bessel transform, the function

is well-defined. Following this,

is a member of $L^2_\nu(0, \infty)$. This demonstrates that for all $y \geq 0$, the function $\mathscr{K}_{\nu, s, \alpha, \gamma}(\cdot, y)$ belongs to $\mathcal{H}^{(s)}_{\nu}$. This establishes part (ⅰ) of the theorem.

Let $\phi \in \mathcal{H}^{(s)}_{\nu}$ and $y \in [0, \infty)$. By Eq (3.3), we have

From the relation

the action of $\mathscr{W}_{\alpha, \nu, t}$ on the kernel $\mathscr{K}_{\nu, s, \alpha, \gamma}(\cdot, y)$ is then:

Therefore

Combining Eqs (3.4) and (3.5), we get

This confirms the reproducing property (ⅱ). □

4.   Tikhonov regularized method

We now consider the variational functional associated with the generalized Weierstrass integral transform $\mathscr{W}_{\alpha, \nu, t},$ defined as

For $\gamma > 0,$ the functional $J_\gamma$ is strictly convex and $J_\gamma(\phi)\geq \frac{\gamma}{2}\|\phi\|_{\nu}$. Hence, $J_\gamma$ has a unique minimizer, which can be characterized by the first-order condition

where $J_\gamma^{'}(\phi)$ is the Fréchet differential of $J_\gamma.$ We denote by $\mathscr{R}_{\gamma, \psi}$ the regularized solution of the Eq (4.2), that is,

The following theorem is our second main result.

Theorem 4.1. For $\nu\geq -1/2,$ $\gamma > 0,$ and $\psi \in L^2_\nu(0, \, \infty).$ Then there is a unique function $\mathscr{R}_{\gamma, \psi}\in \mathcal{H}^{(s)}_{\nu},$ where the infimum of the functional $J_\gamma$, defined by

is attained. Furthermore, the regularized function $\mathscr{R}_{\gamma, \psi}$ is given by

where

Proof. Observe that

where $\Delta \phi$ denotes the increment.

Since,

Taking into account Eq (4.6) and using the Plancheral formula for the Fourier-Bessel transform to get

Hence, the Fréchet differential of $J_\gamma$ can be written as

By the Perseval formula for the Fourier-Bessel transform, it follows that the regularized solution $\mathcal{R}_{\gamma, \psi}(\xi)$ is given by

Therefore

It is easy to see that

By the inversion formula for the Fourier-Bessel transform, we have

We have

where

□ In the following theorem, we will provide an error estimate for the inversion formula.

Theorem 4.2. Let $s > \nu + 1$. For all $\psi_1, \psi_2 \in L^2_\nu(0, \infty)$, the following inequality holds:

Proof. Consider any $\psi_1, \psi_2 \in L^2_\nu(0, \infty)$. The squared norm of the difference between the operators $\mathscr{R}_{\gamma, \psi_1}$ and $\mathscr{R}_{\gamma, \psi_2}$ in the space $\mathcal{H}^{(s)}_{\nu}$ are given by:

Using the formula for the Fourier-Bessel transform of the extremal function $\mathscr{R}_{\gamma, \psi_i}$:

we can express the integral as:

By using the inequality

we can further estimate:

This completes the proof. □

Proposition 4.1. Let $s > \nu + 1$, $\gamma > 0$, and $\psi \in L^2_{\nu}(0, \infty)$. We have the following estimate:

where

Proof. From (4.5) and applying the Cauchy-Schwarz inequality, we have:

Integrating over $[0, \infty)$ with respect to the measure $\sigma_\nu(dx)$, we obtain:

However,

it follows that

Therefore,

We complete the proof by using the relation (3.2) and the fact that:

and

□

Theorem 4.3. Let $s > \nu + 1$. For every $\phi \in \mathcal{H}^{(s)}_{\nu}$ and $\psi = \mathscr{W}_{\alpha, \nu, t}(\phi)$, we have:

Moreover, the set $\{\mathscr{R}_{\gamma, \psi}\}_{\gamma > 0}$ converges uniformly to $\phi$ as $\gamma \to 0^+$.

Proof. Let $\phi \in \mathcal{H}^{(s)}_{\nu}$ and $\psi = \mathscr{W}_{\alpha, \nu, t}(\phi)$. Utilizing the formula given in Eq (4.6), the Fourier-Bessel transform of the extremal function $\mathscr{R}_{\gamma, \psi}$ takes the form:

We can express the norm of the difference between $\mathscr{R}_{\gamma, \psi}$ and $\phi$ in $\mathcal{H}^{(s)}_{\nu}$ as:

Using the dominated convergence theorem and observing that

and given that $\phi \in \mathcal{H}^{(s)}_{\nu}$, we deduce that

The function $\mathscr{F}_\nu(\phi)$ belongs to both $L^1_\nu(0, \, \infty)$ and $L^2_\nu(0, \, \infty)$. Applying the inversion formula for the Fourier-Bessel transform, we compute the deviation of the function $\phi$ under the operator $\mathscr{R}_{\gamma, \psi}$ as follows:

Thus, for all $\xi \in [0, +\infty[$, the magnitude of the deviation is bounded by:

Employing the dominated convergence theorem and noting that:

we deduce that the supremum over all $\xi \geq 0$ of the deviation approaches zero as $\gamma$ tends towards zero

This completes the proof of convergence in the $\mathcal{H}^{(s)}_\nu$ norm and uniform convergence as $\gamma \to 0^+$. □

Example 4.1. As an illustrative example, consider the fractional heat equation on $(0, \infty) \times (0, \infty)$

with the initial condition

To apply the Tikhonov regularization method to this fractional heat equation, we consider the integral operator $\mathscr{W}_{1, -1/2, t} \colon \mathcal{H}^{(1)}_{-1/2} \to L^2(0, \infty)$ defined by

where the Sobolev space $\mathcal{H}^{(1)}_{-1/2}$ is realized by the reproducing kernel Hilbert space $\mathscr{K}_1(x, y)$ given by

where $x, y \geq 0$.

We now consider the following best approximation problem, that is, the Tikhonov functional. For any $\psi \in L^2(0, \infty)$ and $\gamma > 0$, we aim to solve

Then for the RKHS $\mathcal{H}_{-1/2, 1, \gamma}^{(s)}$ consisting of all the members of $\mathcal{H}^{(1)}_{-1/2}$ with the norm

the reproducing kernel $\mathscr{N}_{-1/2, 1, 1, \gamma}(x, y)$ can be calculated directly using the Fourier integrals as follows:

Hence, the unique member of $\mathcal{H}^{(1)}_{-1/2}$ with the minimum $\mathcal{H}^{(1)}_{-1/2}$ norm—the function $\mathscr{R}_{\gamma, \psi}$ which attains the infimum—is given by

For $\phi \in \mathcal{H}^{(1)}_{-1/2}$ and for $\psi = \mathscr{W}_{1, -1/2, t}\phi$, we have the formula

uniformly on $(0, \infty)$.

5.   Conclusions

In this paper, we have explored the generalized Gauss-Weierstrass integral transform associated with the Bessel operator, emphasizing its application to space-fractional diffusion equations. This extension is significant because the Bessel operator coincides with the radial part of the Laplace operator, thereby broadening the scope of the classical transform.

We utilized the properties of the Fourier-Bessel transform and its connection with the $*$-convolution product to demonstrate that the transform $\mathscr{W}_{\alpha, \nu, t}$ is a one-to-one bounded linear operator from the Sobolev space $\mathcal{H}^{s}_\nu$ into $L^2_\nu(0, \infty)$. Given the ill-posed nature of the inverse problem, we applied Tikhonov regularization techniques to achieve stable reconstruction of functions. Our main theorem established the convergence of the regularized solution $\mathscr{R}_{\gamma, \psi}$ to the original function $\phi$ as the regularization parameter $\gamma$ approaches zero.

The implications of our findings extend beyond theoretical interest, offering potential applications in solving inverse problems associated with fractional diffusion equations. Future research could further investigate the numerical implementation of the regularized inversion process and explore other types of fractional differential operators within this framework.

Overall, our work provides a robust foundation for the generalized Gauss-Weierstrass transform's utility in fractional calculus and opens avenues for further investigation into its applications in various mathematical contexts.

Use of AI tools declaration

The author declares he/she has not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

The author would like to extend their sincere appreciation to the Researchers Supporting Project number (RSPD2024R974), King Saud University, Riyadh, Saudi Arabia.

Conflict of interest

The author declares that there is no conflict of interest regarding the publication of this paper.