The Truncation Regularization Method for Identifying the Initial Value on Non-Homogeneous Time-Fractional Diffusion-Wave Equations

Abstract

In the essay, we consider an initial value question for a mixed initial-boundary value of time-fractional diffusion-wave equations. This matter is an ill-posed problem; the solution relies discontinuously on the measured information. The truncation regularization technique is used for restoring the initial value functions. The convergence estimations are given in a priori regularization parameter choice regulations and a posteriori regularization parameter choice regulations. Numerical examples are given to demonstrate this is effective and practicable.

1. Summary

Recently, abnormal diffusion has attracted intensive attention. In [1,2], we learn it has been widely used in semiconductors, nuclear magnetic resonance, quantum optics, molecular spectroscopy, porous media, polymer, turbulence, solid surface diffusion, economic, financial research, etc. The motion of normal diffusion particles is Brownian motion; its mean square and the displacement is linear with the time of movement, which is essentially a Markovian localized motion. In [3], for anomalous diffusion, it is a non-Markovian nonlocalized motion. Therefore, anomalous diffusion must take into account the long-range dependence (space) and historical dependence (time) of the motion process, so the fractional differential equation becomes an important tool for describing this process. Average migration in anomalous diffusion is non-linearly dependent on migration time. When $0<\alpha <1$, the corresponding anomalous diffusion is called subdiffusion. This phenomenon exists in many systems. In [4], it exits in the transportation of charged particles in amorphous semiconductors. In [5], it exits in the porous seepage system. In [6], it exits in the vibration system in polymers. In [7], it exits in the nuclear magnetic resonance in disordered media. In [8], it exits in the transmission in fractal geometry. When $1<\alpha <2$, the corresponding anomalous diffusion is called superdiffusion. This phenomenon exists in many respects. In [9], it exits in the swirl in a special region. In [10], it exits in the collective sliding diffusion on solid surface. In [11], it exits in the stratified velocity field. In [12], it exits Richardson turbulent diffusion. In [13], it exits in quantum optics. In [14], it exits in single molecule spectroscopy. In [15], it exits in the micelle system, transportation in heterogeneous rocks, etc. Based on abnormal superdiffusion of particles, in [16,17], the fractional diffusion-wave equation can accurately characterize the behavior of molecular dynamics of viscoelastic materials, mechanical, electromagnetic and biological reactions, etc. In [18], the fractional diffusion-wave equation can accurately characterize the propagation of mechanical scattering in viscoelastic media with power law creep properties. At present, there are a lot of research findings on the positive matter concerning the fractional-diffusion formula, i.e., the initial value question and initial-boundary value question. However, when solving practical problems by using the fractional diffusion equation, some boundary or initial conditions are unknown, some coefficients in the equation are unknown or the source term is unknown and some additional conditions need to be used to identify these unknown terms; this is the inverse problem of the fractional diffusion equation. For the sake of controlling pollution sources, the mathematical models concerning these problems can be reduced to the inverse problems of source term identification, initial condition reconstruction and diffusion coefficient identification of the fractional diffusion equation. Compared with the positive problems, the inverse problems are mostly ill-posed and need to be solved by the regularization method. At present, there are a lot of study achievements concerning the inverse problem for the fractional order reaction (sub-diffusion) diffusion equation, but the research for the inverse problem of the fractional order reaction (superdiffusion) diffusion equation is still in the preliminary stage. In [19], Yamamoto only gave the existence and uniqueness of solutions of the time fractional wave equation. In [20,21], Siskova considered the problem of identifying time source terms for time fractional inhomogeneous and non-linear wave equations, and this kind of problem is solved by means of the regularization method, but only the existence, uniqueness and a priori estimation formula of the solution are given, the posterior error estimation is not given. In [22], Wei proposed the inverse initial value problems for homogeneous time-fractional wave-equations, and Tikhonov was used to give the regularization; estimates of convergence rate among the regular solution and the analytical solution were given, but the error estimates were saturated. In [23], Yang et al. used the Landweber regularization method to identify the initial value for a time fractional nonhomogeneous wave equation with one measurement, but the nonhomogeneous term only depended on the spatial variable x.

Reference [22] represented a time-fractional diffusion-wave model with homogeneous Dirichlet bound conditions. In the following article, the mixed initial-boundary value problem of the time-fractional diffusion-wave nonhomogeneous equation is studied:

$$\left\{\begin{array}{cc}{D}_t^{\alpha }u=u_{xx}(x,t)+F(x,t),\hfill & 0<x<\pi ,0<t<T,\hfill \\ u(0,t)=u(\pi ,t)=0,\hfill & 0\le t\le T,\hfill \\ u(x,0)=a(x),\hfill & 0\le x\le \pi ,\hfill \\ u_t(x,0)=b(x),\hfill & 0\le x\le \pi ,\hfill \end{array}\right.$$

(1)

when the initial value $a(x)$ and $b(x)$ are known, we can solve $u(x,t)$ by means of nonhomogeneous terms and boundary conditions; this is a typical direct problem. Now that the initial conditions of $a(x)$, $b(x)$ are unknown, we remark that $u(x,T)=g(x)$, and we use a pair of measurement data $(F(x,t),g(x))$ to invert the initial conditions of $a(x)$ and $b(x)$. Here, $D_t^{\alpha }$ is the Caputo time-fractional derivatives with order $\alpha (1<\alpha <2)$, which is defined as [19]:

$$D_t^{\alpha }u=\frac{{\partial }^{\alpha }u}{\partial t^{\alpha }}:=\frac{1}{\Gamma (2-\alpha )}{\int }_0^t\frac{u_{ss}(x,s)}{{(t-s)}^{\alpha -1}}ds,1<\alpha <2,$$

(2)

$\Gamma (x)$ is the function of Gamma.

In Equation (1), the source term $F(x,t)\in L^{\infty }((0,T];L^2(0,\pi ))$ and the final value $g(x)\in L^2(0,\pi )$. The measurement data are all error-bound; assuming the measurement data functions with error are marked as $F^{\delta }$, $g^{\delta }$, they satisfy:

$$\parallel F^{\delta }{-F\parallel }_{L^{\infty }((0,T];L^2(0,\pi ))}\le \delta ,{\parallel g^{\delta }-g\parallel }_{L^2(0,\pi )}\le \delta ,$$

(3)

among them, $\delta >0$ is the error level.

The above inversion initial value problem can be decomposed into two inverse problems.

BP1: When $b(x)$ is known, the unknown initial value function $a(x)$ is inverted using the measured data of $F^{\delta }$, $g^{\delta }$.

BP2: When $a(x)$ is known, the unknown initial value function $b(x)$ is inverted using the measured data of $F^{\delta }$, $g^{\delta }$.

The truncated regularization method is used for resolving the questions BP1 and BP2. In each problem, the exact solution and conditional stability results of the problem are given. A priori error estimates and a posterior error estimates among the regular solution and the analytic solution are also presented. Finally, different numerical instances are displayed to illustrate the feasibility of the regularization mean.

The structure of this essay is as follows: The second part gives auxiliary mathematical conclusions; the third part illustrates that the problem is ill posed and gives the exact solution and conditional stability result; the fourth part gives error estimates among the regularization solutions and the analytic solution; the fifth part gives numerical consequences; and the sixth part is a summary of the article.

2. Some Auxiliary Results

In order to discuss the problem of Equation (1), this section introduces the Mittag–Leffler function and some of its properties.

Definition 1

([24,25]). For any constant $\alpha >0$ and $\beta \in \mathbb{R}$, the Mittag–Leffler function is defined as:

$$E_{\alpha ,\beta }(z)=\sum _{k=0}^{\infty }\frac{z^k}{\Gamma (\alpha k+\beta )},z\in \mathbb{C}.$$

(4)

Lemma 1

([22]). For any $\lambda >0$, $\alpha >0$ and a positive constant $m\in \mathbb{N}$, there are:

$$\begin{array}{cc}\hfill & \frac{d^m}{d{t}^m}{E}_{\alpha ,1}(-\lambda t^{\alpha })=-\lambda t^{\alpha -m}{E}_{\alpha ,\alpha -m+1}(-\lambda t^{\alpha }),t>0,\hfill \end{array}$$

(5)

$$\begin{array}{cc}\hfill & \frac{d}{dt}(t{E}_{\alpha ,2}(-\lambda t^{\alpha }))=E_{\alpha ,1}(-\lambda t^{\alpha }),t\ge 0.\hfill \end{array}$$

(6)

Lemma 2

([23]). Let $0<\alpha <2$, $\beta \in \mathbb{R}$, assume μ satisfies $\frac{\pi \alpha }{2}<\mu <min\{\pi ,\pi \alpha \}$, there is a constant $C=C(\alpha ,\beta ,\mu )>0$ to make:

$$|E_{\alpha ,\beta }(z)|\le \frac{C}{1+|z|},\mu \le |arg(z)|\le \pi .$$

(7)

Lemma 3

([24]). When $1<\alpha <2$, $\beta \in \mathbb{R}$, $\eta >0$, there is:

$$E_{\alpha ,\beta }(-\eta )=\frac{1}{\Gamma (\beta -\alpha )\eta }+O(\frac{1}{{\eta }^2}),\eta \to \infty .$$

(8)

3. Exact Solution and Conditional Stability

In this section, we give the solution for Equation (1) and the properties of Mittag–Leffler functions.

Theorem 1.

Assume $a,b\in L^2(0,\pi )$, $F\in L^{\infty }((0,T];L^2(0,\pi ))$, then the weak solutions of Equation (1) are:

$$u(x,t)={\sum }_{n=1}^{\infty }\{(a,{\phi }_n)E_{\alpha ,1}(-n^2{t}^{\alpha })+(b,{\phi }_n)t{E}_{\alpha ,2}(-n^2{t}^{\alpha })+{\int }_0^t{E}_{\alpha ,\alpha }(-n^2{(t-s)}^{\alpha }){(t-s)}^{\alpha -1}{F}_n(s)ds\}{\phi }_n(x),$$

(9)

where, $u\in C([0,T];L^2(0,\pi ))\cap C((0,T];H^2(0,\pi )\cap H_0^1(0,\pi ))$, $(a,{\phi }_n)$ represents the Fourier coefficients of $a(x)$, $(b,{\phi }_n)$ represents the Fourier coefficients of $b(x)$, $F_n(s)=(F(x,t),{\phi }_n(x))$, ${\phi }_n(x)=\sqrt{\frac{2}{\pi }}sin(nx)$ is a group of standard orthogonal bases under $L^2(0,\pi )$, $n=1,2,\cdots$.

Define:

$$H^p(0,\pi )=\{\psi \in L^2(0,\pi );\sum _{n=1}^{\infty }{(1+n^2)}^p|(\psi ,{\phi }_n){|}^2<\infty \},$$

(10)

here $(·,·)$ is the inner product in the sense of $L^2(0,\pi )$, then the norm of $H^p(0,\pi )$ is:

$${\parallel \psi \parallel }_{H^p(0,\pi )}=(\sum _{n=1}^{\infty }{(1+n^2)}^p|(\psi ,{\phi }_n){{|}^2)}^{\frac{1}{2}}.$$

(11)

Using the measurement data $u(x,T)=g(x)$, we acquire:

$$\begin{array}{cc}\hfill {\displaystyle g(x)=u(x,T)}& =\sum _{n=1}^{\infty }\{(a,{\phi }_n)E_{\alpha ,1}(-n^2{T}^{\alpha })+(b,{\phi }_n)T{E}_{\alpha ,2}(-n^2{T}^{\alpha })\hfill \\ & +{\int }_0^T{E}_{\alpha ,\alpha }(-n^2{(T-s)}^{\alpha }){(T-s)}^{\alpha -1}{F}_n(s)ds\}{\phi }_n(x).\hfill \end{array}$$

(12)

Let:

$$g_1(x)=g(x)-\sum _{n=1}^{\infty }[(b,{\phi }_n)T{E}_{\alpha ,2}(-n^2{T}^{\alpha })+{\int }_0^T{E}_{\alpha ,\alpha }(-n^2{(T-s)}^{\alpha }){(T-s)}^{\alpha -1}{F}_n(s)ds]{\phi }_n(x),$$

$$g_2(x)=g(x)-\sum _{n=1}^{\infty }[(a,{\phi }_n)E_{\alpha ,1}(-n^2{T}^{\alpha })+{\int }_0^T{E}_{\alpha ,\alpha }(-n^2{(T-s)}^{\alpha }){(T-s)}^{\alpha -1}{F}_n(s)ds]{\phi }_n(x),$$

then:

$$\sum _{n=1}^{\infty }(a,{\phi }_n)E_{\alpha ,1}(-n^2{T}^{\alpha }){\phi }_n(x)=g_1(x),$$

(13)

$$\sum _{n=1}^{\infty }(b,{\phi }_n)T{E}_{\alpha ,2}(-n^2{T}^{\alpha }){\phi }_n(x)=g_2(x).$$

(14)

For the sake of solving this question, BP1 and BP2 transfer over to define the following integral equations:

$$\widetilde{K_1}a(\xi ):={\int }_{(0,\pi )}\widetilde{k_1}(x,\xi )a(\xi )d\xi =g_1(x),0\le x\le \pi ,$$

(15)

$$\widetilde{K_2}b(\xi ):={\int }_{(0,\pi )}\widetilde{k_2}(x,\xi )b(\xi )d\xi =g_2(x),0\le x\le \pi ,$$

(16)

combine Equations (13) and (14), then the kernel functions for the operators $\widetilde{K_1}$ and $\widetilde{K_2},$ respectively, are:

$$\widetilde{k_1}(x,\xi )=\sum _{n=1}^{\infty }{E}_{\alpha ,1}(-n^2{T}^{\alpha }){\phi }_n(x){\phi }_n(\xi ),$$

(17)

$$\widetilde{k_2}(x,\xi )=\sum _{n=1}^{\infty }T{E}_{\alpha ,2}(-n^2{T}^{\alpha }){\phi }_n(x){\phi }_n(\xi ).$$

(18)

Based on $\widetilde{k_1}(x,\xi )=\widetilde{k_1}(\xi ,x)$ and $\widetilde{k_2}(x,\xi )=\widetilde{k_2}(\xi ,x)$, we know $\widetilde{K_1}$, $\widetilde{K_2}$ are all self-adjoint operators. From [23], we know $\widetilde{K_1}$ and $\widetilde{K_2}$ are linear compact operators in $L^2(0,\pi )\to L^2(0,\pi )$, thus BP1 and BP2 are both ill-posed problems.

Let ${\widetilde{K_1}}^{*}$ be the adjoint operator of $\widetilde{K_1}$; ${\{{\phi }_n\}}_{n=1}^{\infty }$ is a group of standard orthogonal bases on $L^2(0,\pi )$, and we have:

$${\widetilde{K_1}}^{*}\widetilde{K_1}{\phi }_n(\xi )=E_{\alpha ,1}^2(-n^2{T}^{\alpha }){\phi }_n(\xi ),$$

thus, ${\sigma }_n^{(1)}=|E_{\alpha ,1}(-n^2{T}^{\alpha })|$ is the singular value of $\widetilde{K_1}$. Define:

$$Y_n^{(1)}(x)=\left\{\begin{array}{c}{\phi }_n(x),E_{\alpha ,1}(-n^2{T}^{\alpha })\ge 0,\hfill \\ -{\phi }_n(x),E_{\alpha ,1}(-n^2{T}^{\alpha })<0.\hfill \end{array}\right.$$

We can see $\{Y_n^{(1)}{\}}_{n=1}^{\infty }$ are orthonormal on $L^2(0,\pi )$; then there exist:

$$\widetilde{K_1}{\phi }_n(\xi )={\sigma }_n^{(1)}{Y}_n^{(1)}(x)=E_{\alpha ,1}(-n^2{T}^{\alpha }){\phi }_n(x),$$

$${\widetilde{K_1}}^{*}{Y}_n^{(1)}(x)={\sigma }_n^{(1)}{\phi }_n(\xi )=E_{\alpha ,1}(-n^2{T}^{\alpha })Y_n^{(1)}(\xi ).$$

Thus, the singular system for the operator $\widetilde{K_1}$ is $({\sigma }_n^{(1)};{\phi }_n,Y_n^{(1)})$.

The singular system for the operator $\widetilde{K_2}$ is $({\sigma }_n^{(2)};{\phi }_n,Y_n^{(2)})$. In it ${\sigma }_n^{(2)}=|T{E}_{\alpha ,2}(-n^2{T}^{\alpha })|$:

$$Y_n^{(2)}(x)=\left\{\begin{array}{c}{\phi }_n(x),E_{\alpha ,2}(-n^2{T}^{\alpha })\ge 0,\hfill \\ -{\phi }_n(x),E_{\alpha ,2}(-n^2{T}^{\alpha })<0,\hfill \end{array}\right.$$

Lemma 4

([22]). When $1<\alpha <2$, $T>0$ is a given constant, there are finite points to make $E_{\alpha ,1}(-n^2{T}^{\alpha })=0$ and $E_{\alpha ,2}(-n^2{T}^{\alpha })=0$, respectively. Note the point set to let $E_{\alpha ,1}(-n^2{T}^{\alpha })=0$ is $I_1=\{n_1,n_2,n_3,\cdots ,n_N\}$, the point set to let $E_{\alpha ,2}(-n^2{T}^{\alpha })=0$ is $I_2=\{m_1,m_2,m_3,\cdots ,m_M\}$.

Proof. 

According to Lemma 3, there exists $L_0>0$ to make:

$$E_{\alpha ,1}(-n^2{T}^{\alpha })\le \frac{1}{2\Gamma (1-\alpha )n^2{T}^{\alpha }}<0,n^2{T}^{\alpha }>L_0,$$

if $1<\alpha <2$, then we have $E_{\alpha ,1}(-n^2{T}^{\alpha })=0$ only under $n^2{T}^{\alpha }\le L_0$. Then for ${lim}_{n\to +\infty }{n}^2=+\infty$ there exits finite points n to satisfy $n^2{T}^{\alpha }\le L_0$. The proof of $E_{\alpha ,2}(-n^2{T}^{\alpha })$ is similar to the above process. □

Remark 1.

When $I_1=\varphi$, $I_2=\varphi$, the kernel functions for $\widetilde{K_1}$ are as in Equation (17), the kernel function for $\widetilde{K_2}$ is as in Equation (18). In this paper, this situation is regarded as a special case.

Next we talk about the general situation, which is: $I_1\ne \varphi$, $I_2\ne \varphi$. When $I_1\ne \varphi$, $I_2\ne \varphi$, the kernel function for $\widetilde{K_1}$ and $\widetilde{K_2}$ are, respectively, expressed as,

$$\widetilde{k_1}(x,\xi )=\sum _{n=1,n\notin I_1}^{\infty }{E}_{\alpha ,1}(-n^2{T}^{\alpha }){\phi }_n(x){\phi }_n(\xi ),$$

$$\widetilde{k_2}(x,\xi )=\sum _{n=1,n\notin I_2}^{\infty }T{E}_{\alpha ,2}(-n^2{T}^{\alpha }){\phi }_n(x){\phi }_n(\xi ).$$

Then, the kernel space for the operator $\widetilde{K_1}$ is, when $I_1\ne \varphi$:

$$N(\widetilde{K_1})=span\{{\phi }_n;n\in I_1\}.$$

when $I_1=\varphi$:

$$N(\widetilde{K_1})=\{0\}.$$

The kernel space for the operator $\widetilde{K_2}$ is, when $I_2\ne \varphi$:

$$N(\widetilde{K_2})=span\{{\phi }_n;n\in I_2\},$$

when $I_2=\varphi$:

$$N(\widetilde{K_2})=\{0\},$$

The ranges for $\widetilde{K_1}$ and $\widetilde{K_2}$ are, respectively, written as:

$$R(\widetilde{K_1})=\{g_1\in L^2(0,\pi )|(g_1,{\phi }_n)=0,n\in I_1;\sum _{n=1,n\notin I_1}^{\infty }(\frac{(g_1,{\phi }_n)}{E_{\alpha ,1}(-n^2{T}^{\alpha })}{)}^2<+\infty \},$$

$$R(\widetilde{K_2})=\{g_2\in L^2(0,\pi )|(g_2,{\phi }_n)=0,n\in I_2;\sum _{n=1,n\notin I_2}^{\infty }(\frac{(g_2,{\phi }_n)}{T{E}_{\alpha ,2}(-n^2{T}^{\alpha })}{)}^2<+\infty \}.$$

Lemma 5

([22]). When $1<\alpha <2$, there exist positive constants $\underline{C}$ and $\overline{C}$ which only rely on α, T and $n^2$ to make:

$$\frac{\underline{C}}{n^2}\le |E_{\alpha ,1}(-n^2{T}^{\alpha })|\le \frac{\overline{C}}{n^2},n\notin I_1,$$

(19)

$$\frac{\underline{C}}{n^2}\le |E_{\alpha ,2}(-n^2{T}^{\alpha })|\le \frac{\overline{C}}{n^2},n\notin I_2.$$

(20)

Theorem 2.

When $I_1\ne \varnothing$, $g_1\in R(\widetilde{K_1})$, the integral Equation (15) has infinite solutions, but in $L^2(0,\pi )$ it has only one optimal approximate solution:

$$a(x)=\sum _{n=1,n\notin I_1}^{\infty }\frac{(g_1,{\phi }_n)}{E_{\alpha ,1}(-n^2{T}^{\alpha })}{\phi }_n(x),$$

(21)

when $I_1=\varnothing$, $g_1\in R(\widetilde{K_1})$, the exact solution for BP1 is:

$$a(x)=\sum _{n=1}^{\infty }\frac{(g_1,{\phi }_n)}{E_{\alpha ,1}(-n^2{T}^{\alpha })}{\phi }_n(x).$$

(22)

Proof. 

Making use of Equation (15), it is easy to get the above conclusion. □

Theorem 3.

When $I_2\ne \varnothing$ and $g_2\in R(\widetilde{K_2})$, the integral Equation (16) has infinite solutions, but in $L^2(0,\pi )$ it has only one optimal approximate solution:

$$b(x)=\sum _{n=1,n\notin I_2}^{\infty }\frac{(g_2,{\phi }_n)}{T{E}_{\alpha ,2}(-n^2{T}^{\alpha })}{\phi }_n(x),$$

(23)

when $I_2=\varnothing$, $g_2\in R(\widetilde{K_2})$, the exact solution for BP2 is:

$$b(x)=\sum _{n=1}^{\infty }\frac{(g_2,{\phi }_n)}{T{E}_{\alpha ,2}(-n^2{T}^{\alpha })}{\phi }_n(x).$$

(24)

Proof. 

Making use of Equation (16), it is easy to get the above conclusion. □

Define an a priori boundary for any $a(x)\in H^p(0,\pi )\cap N{(\widetilde{K_1})}^{\perp }$ and we have:

$${\parallel a\parallel }_{H^p(0,\pi )}\le E_1,p>0,E_1>0,$$

(25)

and for any $b(x)\in H^p(0,\pi )\cap N{(\widetilde{K_2})}^{\perp }$, we have:

$${\parallel b\parallel }_{H^p(0,\pi )}\le E_2,p>0,E_2>0.$$

(26)

Theorem 4.

For any $a(x)\in H^p(0,\pi )\cap N{(\widetilde{K_1})}^{\perp }$ which satisfies Equation (25), we get:

$$\parallel a\parallel \le C_1{\parallel \widetilde{K_1}a\parallel }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}},$$

(27)

in it, $C_1={\underline{C}}^{-\frac{p}{p+2}}.$

Proof. 

For any $a(x)\in H^p(0,\pi )\cap N{(\widetilde{K_1})}^{\perp }$, there is:

$$a=\sum _{n=1,n\notin I_1}^{\infty }(a,{\phi }_n){\phi }_n.$$

By $\parallel \widetilde{K_1}{a\parallel }^2={\sum }_{n=1,n\notin I_1}^{\infty }{E}_{\alpha ,1}^2(-n^2{T}^{\alpha }){(a,{\phi }_n)}^2,$ Lemma 5 and the inequality of H$\ddot{o}$lder, we obtain:

$$\begin{array}{cc}\hfill & {\parallel a\parallel }^2=\sum _{n=1,n\notin I_1}^{\infty }{(a,{\phi }_n)}^2\le \sum _{n=1,n\notin I_1}^{\infty }\frac{(|E_{\alpha ,1}(-n^2{T}^{\alpha }){|(a,{\phi }_n))}^{\frac{4}{p+2}}}{E_{\alpha ,1}^2(-n^2{T}^{\alpha })}(|E_{\alpha ,1}(-n^2{T}^{\alpha }){|(a,{\phi }_n))}^{\frac{2p}{p+2}}\hfill \\ \hfill & \le {(\sum _{n=1,n\notin I_1}^{\infty }\frac{{(a,{\phi }_n)}^2}{|E_{\alpha ,1}^p(-n^2{T}^{\alpha })|})}^{\frac{2}{p+2}}\times (\sum _{n=1,n\notin I_1}^{\infty }|E_{\alpha ,1}(-n^2{T}^{\alpha }){{|}^2{(a,{\phi }_n)}^2)}^{\frac{p}{p+2}}\hfill \\ \hfill & \le {(\sum _{n=1,n\notin I_1}^{\infty }({(\frac{n^2}{\underline{C}})}^p{(a,{\phi }_n)}^2))}^{\frac{2}{p+2}}{\parallel \widetilde{K_1}a\parallel }^{\frac{2p}{p+2}}\hfill \\ \hfill & \le {\underline{C}}^{-\frac{2p}{p+2}}{E}_1^{\frac{4}{p+2}}{\parallel \widetilde{K_1}a\parallel }^{\frac{2p}{p+2}}.\hfill \end{array}$$

So:

$$\parallel a\parallel \le {\underline{C}}^{-\frac{p}{p+2}}{\parallel \widetilde{K_1}a\parallel }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}}.$$

 □

Theorem 5.

For any $b(x)\in H^p(0,\pi )\cap N{(\widetilde{K_2})}^{\perp }$ that satisfies Equation (26), we get:

$$\parallel b\parallel \le {(T\underline{C})}^{-\frac{p}{p+2}}{\parallel \widetilde{K_2}b\parallel }^{\frac{p}{p+2}}{E}_2^{\frac{2}{p+2}}.$$

(28)

Proof. 

The method of proof is consistent with Theorem 4. There is no further elaboration here. □

4. Regular Solutions and Error Estimations

Consulting [26,27,28], regular solutions of BP1 and BP2 are given by the truncated regularization method. Based on the singular systems of linear compact operators ${\tilde{K}}_1$ and ${\tilde{K}}_2$:

$$a_{N_1}^{\delta }(x)=\sum _{n=1,n\notin I_1}^{N_1}\frac{(g_1^{\delta },{\phi }_n)}{E_{\alpha ,1}(-n^2{T}^{\alpha })}{\phi }_n(x),$$

(29)

In this, $N_1$ is the regularization parameter,

$$g_1^{\delta }=g^{\delta }-\sum _{n=1}^{\infty }[(b,{\phi }_n)T{E}_{\alpha ,2}(-n^2{T}^{\alpha })+{\int }_0^T{E}_{\alpha ,\alpha }(-n^2{(T-s)}^{\alpha }){(T-s)}^{\alpha -1}{F}_n^{\delta }(s)ds]{\phi }_n.$$

$$b_{N_2}^{\delta }(x)=\sum _{n=1,n\notin I_2}^{N_2}\frac{(g_2^{\delta },{\phi }_n)}{T{E}_{\alpha ,2}(-n^2{T}^{\alpha })}{\phi }_n(x),$$

(30)

In this, $N_2$ is the regularization parameter:

$$g_2^{\delta }=g^{\delta }-\sum _{n=1}^{\infty }[(a,{\phi }_n)E_{\alpha ,1}(-n^2{T}^{\alpha })+{\int }_0^T{E}_{\alpha ,\alpha }(-n^2{(T-s)}^{\alpha }){(T-s)}^{\alpha -1}{F}_n^{\delta }(s)ds]{\phi }_n.$$

Next, the convergence estimate among the optimal approximate solution $a(x)$ and the regular solution $a_{N_1}^{\delta }(x)$ is given under the a priori regularization parameter selection rule.

Lemma 6

([29]). Suppose $F\in L^{\infty }((0,T];L^2(0,\pi ))$, $0<d<4$, there is a constant M to make:

$$\begin{array}{cc}\hfill & \parallel g(x)-\sum _{n=1}^{\infty }{\int }_0^T{E}_{\alpha ,\alpha }(-n^2{(T-s)}^{\alpha }){(T-s)}^{\alpha -1}{F}_n(s)ds{\phi }_n(x)\parallel ,\hfill \\ \hfill & \le \sqrt{2(\parallel g{\parallel }^2+M\parallel F{\parallel }_{L^{\infty }((0,T];L^2(0,\pi ))}^2)},\hfill \end{array}$$

(31)

In this, $F_n(s)=(F(x,t),{\phi }_n(x))$, the constant $M=\frac{1}{c^2}{\sum }_{n=1}^{\infty }\frac{1}{n^{\frac{4}{d}}}$ and c is a positive constant which only relies on n.

Define an orthogonal projection operator $Q_1:L^2(0,\pi )\to \overline{R(\widetilde{K_1})}$, making use of Equations (3) and (31), we have:

$$\parallel Q_1{g}_1^{\delta }-Q_1{g}_1\parallel \le \parallel g_1^{\delta }-g_1\parallel =\parallel g_1^{\delta }-g_1\parallel \le \sqrt{2(M+1)}\delta .$$

Define an orthogonal projection operator $Q_2:L^2(0,\pi )\to \overline{R(\widetilde{K_2})}$; making use of Equations (3) and (31), we have:

$$\parallel Q_2{g}_2^{\delta }-Q_2{g}_2\parallel \le \parallel g_2^{\delta }-g_2\parallel =\parallel g_2^{\delta }-g_2\parallel \le \sqrt{2(M+1)}\delta .$$

Theorem 6.

$a(x)$ is as in Equation (21); the regularization solution $a_{N_1}^{\delta }(x)$ is as in Equation 29); assuming $a(x)$ satisfies Equation (25) and the measurement error hypothesis Equation (3) is held, then the regular parameter is:

$$N_1=[{(\frac{E_1}{\delta })}^{\frac{1}{p+2}}],$$

In this, $[{(\frac{E_1}{\delta })}^{\frac{1}{p+2}}]$ is the biggest integer and is no more than ${(\frac{E_1}{\delta })}^{\frac{1}{p+2}}$, the estimate between $a(x)$ and $a_{N_1}^{\delta }(x)$ is:

$$\parallel a_{N_1}^{\delta }-a\parallel \le C_2{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}},$$

(32)

In this, $C_2=1+\frac{\sqrt{2(M+1)}}{\underline{C}}$.

Proof. 

By the trigonometric inequality, we have:

$$\parallel a_{N_1}^{\delta }-a\parallel \le \parallel a_{N_1}^{\delta }-a_{N_1}\parallel +\parallel a_{N_1}-a\parallel .$$

First, we prove that the first term is derived from Lemmas 5 and 6:

$$\begin{array}{cc}\hfill & \parallel a_{N_1}^{\delta }-a_{N_1}{\parallel }^2=\parallel \sum _{n=1,n\notin I_1}^{N_1}(\frac{{g_1}_n-{g_1^{\delta }}_n}{E_{\alpha ,1}(-n^2{T}^{\alpha })}){\phi }_n(x){\parallel }^2\hfill \\ \hfill & \le [\underset{1\le n\le N_1,n\notin I_1}{sup}(\frac{1}{E_{\alpha ,1}(-n^2{T}^{\alpha })}){]}^2\times \sum _{n=1,n\notin I_1}^{N_1}{({g_1}_n-{g_1^{\delta }}_n)}^2\hfill \\ \hfill & \le [\underset{1\le n\le N_1,n\notin I_1}{sup}(\frac{1}{\frac{\underline{C}}{n^2}}){]}^2(2(M+1)){\delta }^2\le \frac{2(M+1)}{{\underline{C}}^2}{[\underset{1\le n\le N_1,n\notin I_1}{sup}\left(n^2\right)]}^2{\delta }^2\hfill \\ \hfill & \le \frac{2(M+1)}{{\underline{C}}^2}{N}_1^4{\delta }^2.\hfill \end{array}$$

In this, ${g_1}_n=(g_1,{\phi }_n)$, ${g_1^{\delta }}_n=(g_1^{\delta },{\phi }_n)$. Then:

$$\parallel a_{N_1}^{\delta }-a_{N_1}\parallel \le \frac{\sqrt{2(M+1)}}{\underline{C}}{N}_1^2\delta .$$

(33)

According to Equation (25), we estimate the second term of trigonometric inequality:

$$\begin{array}{cc}\hfill & \parallel a_{N_1}{-a\parallel }^2=\parallel \sum _{n=N_1+1,n\notin I_1}^{\infty }(\frac{{g_1}_n}{E_{\alpha ,1}(-n^2{T}^{\alpha })}){\phi }_n(x){\parallel }^2\hfill \\ \hfill & =\sum _{n=N_1+1,n\notin I_1}^{\infty }{(a_n)}^2=\sum _{n=N_1+1,n\notin I_1}^{\infty }(a_n{(1+n^2)}^{\frac{p}{2}}{(1+n^2)}^{\frac{-p}{2}}{)}^2\hfill \\ \hfill & \le [\underset{n\ge N_1+1,n\notin I_1}{sup}({(1+n^2)}^{\frac{-p}{2}}){]}^2{E}_1^2\le N_1^{-2p}{E}_1^2.\hfill \end{array}$$

Thus:

$$\parallel a_{N_1}-a\parallel \le N_1^{-p}{E}_1.$$

(34)

Then:

$$\parallel a_{N_1}^{\delta }-a\parallel \le C_2{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}},$$

In this, $C_2=1+\frac{\sqrt{2(M+1)}}{\underline{C}}$.

Next, the rule of parameter selection for posterior error regularization is given by the Morozov principle. We define the following operator:

$${P_1}_m{g}_1^{\delta }=\sum _{n=1,n\notin I_1}^m{g_1}_n^{\delta }{\phi }_n(x),$$

(35)

and select the regularization parameter to satisfy:

$$\parallel (Q_1-{P_1}_{N_1})g_1^{\delta }\parallel \le {\tau }_1\delta \le \parallel (Q_1-{P_1}_{N_1-1})g_1^{\delta }\parallel ,$$

(36)

In this, ${\tau }_1>\sqrt{2(M+1)}$ is a constant, $\parallel Q_1{g}_1^{\delta }\parallel >{\tau }_1\delta >0$. □

Lemma 7.

Suppose $\rho (N_1):=\parallel (Q_1-{P_1}_{N_1})g_1^{\delta }\parallel$; it has the following properties:

(a) 

$\rho (N_1)$ is a successive function;

(b) 

$\underset{N_1\to 0}{lim}\rho (N_1)=\parallel Q_1{g}_1^{\delta }\parallel ;$

(c) 

$\underset{N_1\to \infty }{lim}\rho (N_1)=0;$

(d) 

$\rho (N_1)$ is a strictly decreasing function.

Lemma 8.

Suppose Equations (3) and (25) are held, then the regularization parameter $N_1$ is obtained by the posterior regularization parameter selection rule, and we have:

$$N_1\le {\overline{C}}^{\frac{1}{p+2}}{(\frac{E_1}{({\tau }_1-\sqrt{2(M+1)})\delta })}^{\frac{1}{p+2}}.$$

(37)

Proof. 

On the one hand, according to Lemma 5 and Equations (15) and (35), we have:

$$\begin{array}{cc}& {\displaystyle \parallel (Q_1-{P_1}_{N_1-1})g_1{\parallel }^2=\parallel \sum _{n=N_1,n\notin I_1}^{\infty }({g_1}_n){\phi }_n(x){\parallel }^2=\sum _{n=N_1,n\notin I_1}^{\infty }{(a_n{E}_{\alpha ,1}(-n^2{T}^{\alpha }))}^2}\hfill \\ & =\sum _{n=N_1,n\notin I_1}^{\infty }(a_n{E}_{\alpha ,1}(-n^2{T}^{\alpha }){(1+n^2)}^{\frac{p}{2}}{(1+n^2)}^{\frac{-p}{2}}{)}^2\hfill \\ \hfill & \le [\underset{n\ge N_1,n\notin I_1}{sup}({(1+n^2)}^{\frac{-p}{2}}(E_{\alpha ,1}(-n^2{T}^{\alpha }))){]}^2\times \sum _{n=N_1,n\notin I_1}^{\infty }(a_n{(1+n^2)}^{\frac{p}{2}}{)}^2\hfill \\ \hfill & \le [\underset{n\ge N_1,n\notin I_1}{sup}({(1+n^2)}^{\frac{-p}{2}}(\frac{\overline{C}}{n^2})){]}^2{E}_1^2\le {\overline{C}}^2\frac{1}{N_1^{2p+4}}{E}_1^2,\hfill \end{array}$$

then:

$$\parallel (Q_1-{P_1}_{N_1-1})g_1\parallel \le \overline{C}{N}_1^{-(p+2)}{E}_1.$$

(38)

On the other hand, according to Lemma 6 and Equation (36), we have:

$$\begin{array}{cc}\hfill {\displaystyle \parallel (Q_1-{P_1}_{N_1-1})g_1\parallel }& =\parallel ({P_1}_{N_1-1}-Q_1)g_1^{\delta }-(Q_1-{P_1}_{N_1-1})(g_1-g_1^{\delta })\parallel \hfill \\ & \ge \parallel ({P_1}_{N_1-1}-Q_1)g_1^{\delta }\parallel -\parallel (Q_1-{P_1}_{N_1-1})(g_1-g_1^{\delta })\parallel \hfill \\ & \ge {\tau }_1\delta -\sqrt{2(M+1)}\delta =({\tau }_1-\sqrt{2(M+1)})\delta ,\hfill \end{array}$$

then:

$$\parallel (Q_1-{P_1}_{N_1-1})g_1\parallel \ge ({\tau }_1-\sqrt{2(M+1)})\delta .$$

(39)

Combine Equations (38) and (39); then:

$$N_1\le {\overline{C}}^{\frac{1}{p+2}}{(\frac{E_1}{({\tau }_1-\sqrt{2(M+1)})\delta })}^{\frac{1}{p+2}}.$$

 □

Lemma 9.

Suppose $a_{N_1}^{\delta }(x)$ is a regular solution of BP1 corresponding to measurement data with errors and $a_{N_1}(x)$ is the regular solution corresponding to the measurement data without errors; then we have:

$$\parallel a_{N_1}^{\delta }-a_{N_1}\parallel \le C_3{({\tau }_1-\sqrt{2(M+1)})}^{-\frac{2}{p+2}}{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}},$$

(40)

In this, $C_3=(C_2-1){\overline{C}}^{\frac{2}{p+2}}$.

Proof. 

Applying Equation (33) and Lemma 8, we get:

$$\begin{array}{cc}& \parallel a_{N_1}^{\delta }-a_{N_1}\parallel \le N_1^2\delta \frac{\sqrt{2(M+1)}}{\underline{C}}\le ({\overline{C}}^{\frac{1}{p+2}}{(\frac{E_1}{({\tau }_1-\sqrt{2(M+1)})\delta })}^{\frac{1}{p+2}}{)}^2\delta \frac{\sqrt{2(M+1)}}{\underline{C}}\hfill \\ & =\frac{\sqrt{2(M+1)}}{\underline{C}}{\overline{C}}^{\frac{2}{p+2}}{({\tau }_1-\sqrt{2(M+1)})}^{-\frac{2}{p+2}}{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}}\hfill \\ & =(C_2-1){\overline{C}}^{\frac{2}{p+2}}{({\tau }_1-\sqrt{2(M+1)})}^{-\frac{2}{p+2}}{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}}\hfill \\ & =C_3{({\tau }_1-\sqrt{2(M+1)})}^{-\frac{2}{p+2}}{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}}.\hfill \end{array}$$

 □

Lemma 10.

$a(x)$ is as in Equation (21) and $a_{N_1}(x)$ is the regular solution corresponding to the measurement data without errors; suppose Equation (3) and the a priori boundary Equation (25) are held; then we have:

$$\parallel a_{N_1}-a\parallel \le C_1{({\tau }_1+\sqrt{2(M+1)})}^{\frac{p}{p+2}}{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}}.$$

(41)

Proof. 

Making use of Lemmas 5 and 6, Equation (36) and the H$\ddot{o}$lder inequality, we obtain:

$$\begin{array}{cc}\hfill {\displaystyle \parallel a_{N_1}{-a\parallel }^2}& =\parallel \sum _{n=1,n\notin I_1}^{N_1}(\frac{{g_1}_n}{E_{\alpha ,1}(-n^2{T}^{\alpha })}){\phi }_n(x)-\sum _{n=1,n\notin I_1}^{\infty }(\frac{{g_1}_n}{E_{\alpha ,1}(-n^2{T}^{\alpha })}){\phi }_n(x){\parallel }^2\hfill \\ & =\sum _{n=N_1+1,n\notin I_1}^{\infty }(\frac{{g_1}_n}{E_{\alpha ,1}(-n^2{T}^{\alpha })}{)}^2=\sum _{n=N_1+1,n\notin I_1}^{\infty }(\frac{{g_1}_n^{\frac{p}{p+2}}{g_1}_n^{\frac{2}{p+2}}}{E_{\alpha ,1}(-n^2{T}^{\alpha })}{)}^2\hfill \\ & \le (\sum _{n=N_1+1,n\notin I_1}^{\infty }({g_1}_n^{\frac{2p}{p+2}}{)}^{\frac{p+2}{p}}{)}^{\frac{p}{p+2}}(\sum _{n=N_1+1,n\notin I_1}^{\infty }(\frac{{g_1}_n^{\frac{4}{p+2}}}{E_{\alpha ,1}^2(-n^2{T}^{\alpha })}{)}^{\frac{p+2}{2}}{)}^{\frac{2}{p+2}}\hfill \\ & =(\sum _{n=N_1+1,n\notin I_1}^{\infty }{({g_1}_n)}^2{)}^{\frac{p}{p+2}}(\sum _{n=N_1+1,n\notin I_1}^{\infty }({g_1}_n^2\frac{1}{E_{\alpha ,1}^{p+2}(-n^2{T}^{\alpha })}){)}^{\frac{2}{p+2}}\hfill \\ & \le (\sum _{n=N_1+1,n\notin I_1}^{\infty }({g_1}_n-{g_1}_n^{\delta }+{g_1}_n^{\delta }){)}^{\frac{2p}{p+2}}(\sum _{n=N_1+1,n\notin I_1}^{\infty }(a_n^2\frac{1}{E_{\alpha ,1}^p(-n^2{T}^{\alpha })}){)}^{\frac{2}{p+2}}\hfill \\ & \le (\sqrt{2(M+1)}\delta +{\tau }_1\delta {)}^{\frac{2p}{p+2}}(\sum _{n=N_1+1,n\notin I_1}^{\infty }(a_n^2{(1+n^2)}^p\frac{1}{E_{\alpha ,1}^p(-n^2{T}^{\alpha }){(1+n^2)}^p}){)}^{\frac{2}{p+2}}\hfill \\ & \le (\sqrt{2(M+1)}\delta +{\tau }_1\delta {)}^{\frac{2p}{p+2}}{E}_1^{\frac{4}{p+2}}(\underset{n\ge N_1+1,n\notin I_1}{sup}\frac{1}{{(1+n^2)}^p{E}_{\alpha ,1}^p(-n^2{T}^{\alpha })}{)}^{\frac{2}{p+2}}\hfill \\ & \le (\sqrt{2(M+1)}\delta +{\tau }_1\delta {)}^{\frac{2p}{p+2}}{E}_1^{\frac{4}{p+2}}(\underset{n\ge N_1+1,n\notin I_1}{sup}\frac{1}{{(1+n^2)}^p{(\frac{\underline{C}}{n^2})}^p}{)}^{\frac{2}{p+2}}\hfill \\ & \le (\sqrt{2(M+1)}+{\tau }_1{)}^{\frac{2p}{p+2}}(\frac{1}{\underline{C}}{)}^{\frac{2p}{p+2}}{\delta }^{\frac{2p}{p+2}}{E}_1^{\frac{4}{p+2}}.\hfill \end{array}$$

Then:

$$\parallel a_{N_1}-a\parallel \le {(\frac{1}{\underline{C}})}^{\frac{p}{p+2}}{({\tau }_1+\sqrt{2(M+1)})}^{\frac{p}{p+2}}{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}}.$$

 □

Theorem 7.

Suppose Equation (3) and the a priori boundary Equation (25) are held and the regularization parameter is given by Equation (36); then the a posteriori estimation between $a(x)$ and $a_{N_1}^{\delta }(x)$ is:

$$\parallel a_{N_1}^{\delta }-a\parallel \le [C_3{({\tau }_1-\sqrt{2(M+1)})}^{-\frac{2}{p+2}}+C_1{({\tau }_1+\sqrt{2(M+1)})}^{\frac{p}{p+2}}]{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}}.$$

(42)

Proof. 

By aid of trigonometric inequality and Equations (40) and (41), we get:

$$\begin{array}{cc}\hfill {\displaystyle \parallel a_{N_1}^{\delta }-a\parallel }& =\parallel a_{N_1}^{\delta }-a_{N_1}+a_{N_1}-a\parallel \hfill \\ & \le \parallel a_{N_1}^{\delta }-a_{N_1}\parallel +\parallel a_{N_1}-a\parallel \hfill \\ & \le C_3{({\tau }_1-\sqrt{2(M+1)})}^{-\frac{2}{p+2}}{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}}+C_1{({\tau }_1+\sqrt{2(M+1)})}^{\frac{p}{p+2}}{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}}\hfill \\ & =[C_3{({\tau }_1-\sqrt{2(M+1)})}^{-\frac{2}{p+2}}+C_1{({\tau }_1+\sqrt{2(M+1)})}^{\frac{p}{p+2}}]{\delta }^{\frac{p}{p+2}}{E}_1^{\frac{2}{p+2}}.\hfill \end{array}$$

 □

The following theorems give a priori error estimates and a posteriori error estimates of problem BP2. The proof method is similar to problem BP1, so only the result is given and the proof process is not given.

Theorem 8.

$b(x)$ is as in Equation (23) and $b_{N_2}^{\delta }(x)$ is the regular solution corresponding to the measured data with errors; assuming Equations (3) and (26) are held, then the regularization parameter is:

$$N_2=[{(\frac{E_2}{\delta })}^{\frac{1}{p+2}}]$$

In this, $[{(\frac{E_2}{\delta })}^{\frac{1}{p+2}}]$ is the biggest integer and is no more than ${(\frac{E_2}{\delta })}^{\frac{1}{p+2}}$; the estimate is:

$$\parallel b_{N_2}^{\delta }-b\parallel \le C_4{\delta }^{\frac{p}{p+2}}{E}_2^{\frac{2}{p+2}},$$

(43)

In this, $C_4=1+\frac{\sqrt{2(M+1)}}{\underline{C}T}$.

A posteriori error estimates of $b(x)$ are given by the Morozov principle. We define the following operator:

$${P_2}_m{g}_2^{\delta }=\sum _{n=1,n\notin I_2}^m{g_2^{\delta }}_n{\phi }_n(x),$$

(44)

In this, ${g_2^{\delta }}_n=(g_2^{\delta },{\phi }_n)$; select regularization parameters to satisfy:

$$\parallel (Q_2-{P_2}_{N_2})g_2^{\delta }\parallel \le {\tau }_2\delta \le \parallel (Q_2-{P_2}_{N_2-1})g_2^{\delta }\parallel ,$$

(45)

In this, ${\tau }_2>\sqrt{2(M+1)}$ is a constant. Similar to Lemma 4.3, suppose $\parallel Q_2{g}_2^{\delta }\parallel >{\tau }_2\delta$.

Theorem 9.

$b(x)$ is as in Equation (23) and $b_{N_2}^{\delta }(x)$ is the regular solution corresponding to the measured data with errors; suppose Equations (3) and (26) are held, if the regularization parameter is given by Equation (45), the error estimate is:

$$\parallel b_{N_2}^{\delta }-b\parallel \le [C_5{({\tau }_2-\sqrt{2(M+1)})}^{-\frac{2}{p+2}}+C_6{({\tau }_2+\sqrt{2(M+1)})}^{\frac{p}{p+2}}]{\delta }^{\frac{p}{p+2}}{E}_2^{\frac{2}{p+2}},$$

(46)

In this, $C_5={(T\overline{C})}^{\frac{2}{p+2}}(C_4-1)$, $C_6={(\frac{1}{T\underline{C}})}^{\frac{p}{p+2}}$.

5. Numerical Implementation

In this section, different numeric instances are selected to illustrate the effectiveness of the truncated regularization technique for inversion of initial value functions concerning the time-fractional nonhomogeneous wave equation.

By resolving the direct question of Equation (1), we get the measurement information $g(x)$,

$$\left\{\begin{array}{cc}{D}_t^{\alpha }u=u_{xx}(x,t)+F(x,t),\hfill & 0<x<\pi ,0<t<T,\hfill \\ u(0,t)=u(\pi ,t)=0,\hfill & 0\le t\le T,\hfill \\ u(x,0)=a(x),\hfill & 0\le x\le \pi ,\hfill \\ u_t(x,0)=b(x),\hfill & 0\le x\le \pi ,\hfill \end{array}\right.$$

(47)

The finite difference method is used to discretize the above equation; the step for time and space are $\Delta t=\frac{T}{N}$ and $\Delta x=\frac{\pi }{M}$. Thus, we apply $t_n=n\Delta t(n=0,1,2,\cdots ,N)$ to represent the grid point in time domain $[0,T]$ and we apply $x_i=i\Delta x(n=0,1,2,\cdots ,M)$ to represent the grid point in space domain $[0,\pi ]$. The approximate value of each grid point u is denoted $u_i^n\approx u(x_i,t_n)$. Referring to [30], we can write the discrete form of the time-fractional derivative:

$$D_t^{\alpha }{u}_i^n=\frac{{(\Delta t)}^{-\alpha }}{\Gamma (3-\alpha )}(u_i^n-u_i^{n-1})-\frac{{(\Delta t)}^{-\alpha }}{\Gamma (3-\alpha )}\sum _{j=1}^{n-1}(B_{n-j-1}-B_{n-j})(u_i^j-u_i^{j-1})-\frac{{(\Delta t)}^{1-\alpha }}{\Gamma (3-\alpha )}{B}_{n-1}b(x_i),$$

(48)

here $B_k={(k+1)}^{2-\alpha }-k^{2-\alpha },k=0,1,2,\cdots .$ The final value $g(x)$ is calculated by combining the measured non-homogeneous term $F(x,t)$. The $Randn(·)$ function is used to add random perturbations to $g(x)$ and $F(x,t)$:

$$g^{\delta }=g+\epsilon (2randn(size(g))-1),F^{\delta }=F+\epsilon (2randn(size(F))-1).$$

(49)

where, the absolute error among the analytic solution and the regularized solution are:

$$e_{a_1}=\sqrt{\frac{1}{(M+1)}\sum _{i=1}^{M+1}(a_i-a_i^{N_1,\delta }{)}^2},e_{a_2}=\sqrt{\frac{1}{(M+1)}\sum _{i=1}^{M+1}(b_i-b_i^{N_2,\delta }{)}^2}.$$

(50)

Because the a priori bound is not well defined, only the numerical results under the rules of selecting a posteriori regularization parameters are given here. In the next instances, $M=100,N=50$, $F(x,t)=(1+x^2)exp(t)$ and $T=1$; for the sake of inverting the unknown function $a(x)$, we order $b(x)=sin(6\pi x)$; for the sake of inverting the unknown $b(x)$, we order $a(x)=xsin(x)$.

Figure 1 shows the absolute error between the exact solution $a(x)$ and the regularized solution under $\alpha =1.2$ and various $\epsilon$ for Instance 1. Figure 2 shows the exact solution $a(x)$ and the regularized solution for Instance 1 with difference $\alpha$. Figure 3 shows the exact solution $a(x)$ and the regularized solution for Instance 2 with difference $\alpha$ Figure 4 shows the exact solution $a(x)$ and the regularized solution for Instance 3 with difference $\alpha$ Figure 5 shows the exact solution $b(x)$ and the regularized solution for Instance 4 with difference $\alpha$ Figure 6 shows the exact solution $b(x)$ and the regularized solution for Instance 5 with difference $\alpha$ Figure 7 shows the exact solution $b(x)$ and the regularized solution for Instance 6 with difference $\alpha$.

Instance 1.

Let:

$$a(x)=xsin(x).$$

Table 1. Numerical results of Instance 1 for different $\alpha$ with $\epsilon =0.01$.

$\alpha$	1.1	1.3	1.5	1.7	1.9
$e_{a_1}$	0.0599	0.0935	0.1406	0.5209	1.0032

shows the relative error between the exact solution and the regularization solution as α increases. From Table 1, we can see the relative error between the exact solution and the regularization solution increases as α increases.

Instance 2.

Let:

$$a(x)=2x(0<x<\frac{\pi }{2})+2(\pi -x)(\frac{\pi }{2}\le x<\pi ).$$

Instance 3.

Let:

$$a(x)=0(0<x<\frac{\pi }{3})+1(\frac{\pi }{3}\le x<\frac{2\pi }{3})+0(\frac{2\pi }{3}\le x<\pi ).$$

Instance 4.

Let:

$$b(x)=2sin(x).$$

Instance 5.

Let:

$$b(x)=0(0<x\le \frac{\pi }{4})+4(x-\frac{\pi }{4})(\frac{\pi }{4}<x\le \frac{\pi }{2})-4(x-\frac{3\pi }{4})(\frac{\pi }{2}<x\le \frac{3\pi }{4})+0(\frac{3\pi }{4}<x\le \pi ).$$

Instance 6.

Let:

$$b(x)=0(0<x\le \frac{\pi }{2})+1(\frac{\pi }{2}<x\le \pi ).$$

As can be seen from the numerical examples above, the smaller the time order $\alpha$ and noise level $\epsilon$, the better the numerical results will be. Examples 1 and 4 outperform Examples 2, 3, 5 and 6, because their functions are better chosen.

6. Conclusions

In this paper, we study the question of the inversion of initial information for the time-fractional nonhomogeneous diffusion-wave formula, and choose truncated regularization means to resolve this question. The exact solution, conditional stability and a priori and a posteriori error estimates are given; some suitable numerical instances are selected to show the validity of this method. Compared with other scholars’ work, the convergence estimates in this paper are not saturated, and the non-homogeneous terms are both related to time and space. At the same time, combined with graphs and tables, the numerical results illustrate that the smaller the order $\alpha$ and noise level $\epsilon$, the better the numerical results will be.