International Journal of Partial Differential Equations

Volume 2013 (2013), Article ID 590737, 6 pages

http://dx.doi.org/10.1155/2013/590737

*A Posteriori* Regularization Parameter Choice Rule for Truncation Method for Identifying the Unknown Source of the Poisson Equation

School of Science, Lanzhou University of Technology, Lanzhou, Gansu 730050, China

Received 30 April 2013; Accepted 19 August 2013

Academic Editor: Athanasios N. Yannacopoulos

Copyright © 2013 Xiao-Xiao Li and Dun-Gang Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We consider the problem of determining an unknown source which depends only on one variable in two-dimensional Poisson equation. We prove a conditional stability for this problem. Moreover, we propose a truncation regularization method combined with an *a posteriori* regularization parameter choice rule to deal with this problem and give the corresponding convergence estimate. Numerical results are presented to illustrate the accuracy and efficiency of this method.

#### 1. Introduction

Inverse source problems arise in many branches of science and engineering, for example, heat conduction, crack identification, electromagnetic theory, geophysical prospecting, and pollutant detection. In this paper, we consider the following inverse problem: to find a pair of functions satisfying where is the unknown source depending only on one spatial variable and is the supplementary condition. In applications, input data can only be measured, and there will be measured data function which is merely in and satisfies where the constant represents a bound on the measurement error.

For the heat equation, there has been a large number of research results for the different forms of heat source [1–6]. In [7], the authors identified the unknown source of the Poisson equation using the modified regularization method. In [8], the authors identified the unknown source of the Poisson equation using the truncation method. In [7, 8], the regularization parameters which depend on the noise level and the *a priori* bound are selected by the *a priori* rule. Generally speaking, there is a defect for any *a priori* method; that is, the *a priori* choice of the regularization parameter depends seriously on the *a priori* bound of the unknown solution. However, in general, the *a priori* bound cannot be known exactly in practice, and working with a wrong constant may lead to the bad regularized solution.

In the present paper, *a posteriori* choice of the regularization parameter will be given. To the authors’ knowledge, there are few papers for choosing the regularization parameter by the *a posteriori* rule for this problem.

The truncation regularization methods have been studied for solving various types of inverse problems. Eldén et al. [9] used the truncation method to analyze and compute one-dimensional IHCP, Xiong et al. [10] used it to consider the surface heat flux for the sideways heat equation, Fu et al. [11] used it to solve the BHCP, Qian et al. [12] used it to consider the numerical differentiation, and Regińska and Regiński [13] applied the idea of truncation to a Cauchy problem for the Helmholtz equation. In [9–13], the ill-posedness of the problem was caused by the high frequency components, and they all used the truncation method to eliminate all high frequencie and called this truncation method Fourier method. In [14], the authors ever identified the unknown source on Poisson equation on half unbounded domain using Fourier Transform. In this paper, we identified the unknown source on Poisson equation on half band domain using separation of variables and gave the numerical example to illustrate our methods. In [15], the authors ever used this method to consider a Cauchy problem for the Laplace equation in the bounded domain and obtained the error estimate which had a good convergence rate. The truncation method of this paper seems to be more direct and simple for considering the inverse problem in the bounded or half band domain.

This paper is organized as follows. The ill-posedness of this problem (1) is analyzed in Section 2. In Section 3, we give a conditional stable result. In Section 4, we propose the truncation method and give the convergence estimates under an *a priori* assumption for the exact solution and the *a posteriori* regularization parameter choice rule. In Section 5, a numerical example is proposed to show the effectiveness for this method. Section 6 puts an end to this paper with a brief conclusion.

#### 2. Ill-Posedness of Problem (1)

By separation of variables, we obtain that the solution of problem (1) has the form where is an orthogonal basis in , and Using the supplementary condition, we have Defining the operator , we obtain It is easy to see that is a linear compact operator, and the singular values of satisfy That is, Therefore,

Note that as ; thus, the exact data function must satisfy the property that decays rapidly as . As for the measured data function , we cannot expect that it has the same decay rate in . Thus, the problem (1) is ill-posed. It is impossible to gain the unknown source using classical methods. In the following section, we will use the truncation regularization method to deal with the ill-posed problem. Before doing that, we impose an a priori bound on the heat source; that is,
where is a constant, denotes the norm in Sobolev space, is defined by [16] as following:
For any ill-posed problem, the a-priori bound of the exact solution is very necessary. Otherwise, the regularization solution is not convergent or the rate of convergence is very slow. The detail of the a *priori* bound can be seen in [17].

#### 3. A Conditional Stability Estimate

We firstly establish a conditional stability estimate for the inverse source problem (1).

Lemma 1. *If , one has the following inequality:
*

Theorem 2. *Assume that is the solution of (1) and condition (11) is satisfied; then, the following estimate holds:
*

*Proof. *Using the Hölder inequality and (11), we obtain
So

*Remark 3. *If are the solutions of the inverse source problem with exact data , respectively, then there holds
It is obvious that if , then .

#### 4. An *A Posteriori* Regularization Parameter Choice Rule for the Truncation Method and Convergence Estimate

Noting (10), small errors in the components of large can blow up and completely destroy the solution. A nature way to stabilize the problem is to eliminate all the components of large from the solution and instead consider (10) only for . Then, we get a regularized solution Noting (18), if the parameter is large, is close to the exact solution . On the other hand, if the parameter is fixed, is bounded. So the positive integer plays the role of regularization parameter. We consider an a posteriori regularization parameter choice by the discrepancy principle. Define Due to the discrepancy principle, we will take as the solution of where is a constant. Before giving the main conclusion of this section, we first give an important lemma.

Lemma 4. *Let condition (2) hold. If is taken as the solution of (20), then one has the following inequality:
*

*Proof. *Due to (2) and (19), we obtain
So
On the other hand,
So
Combining (23) with (25), we obtain

The main conclusion of this paper is as the following.

Theorem 5. *Let be the regularized solution given by (18), and let be the exact solution given by (10). Let be the measured data at satisfying (2), and priori condition (11) holds for . If the regularization parameter is chosen as the solution of (20), then one gets the following estimate:
**
where .*

* Proof. *Due to the triangle inequality, we obtain
where . It is easy to know
Due to triangle inequality and (20), we have
Using Theorem 2, we obtain
Using (18), we obtain
Due to (13), (21), and (32), we obtain
Combining (28), (31) and (33), we obtain
where .

#### 5. Numerical Implementation

From (7), we know that

We use the trapezoidal rule to approach the integral and do an approximate truncation for the series by choosing the sum of the front terms. After considering an equidistant grid ), we get where The bisection method is used to solve (20) and .

*Example 6. *It is easy to see that the function and the function are the exact solution of the problem (1). Consequently, the data function is .

Adding a random distributed perturbation to each data function, we obtain vector ; that is,
where indicates the noise level of the measurement data and the function “” generates arrays of random numbers whose elements are normally distributed with mean , variance , and standard deviation . “rand (size())” returns an array of random entries that is the same size as . The bound on the measurement error can be measured in the sense of root mean square error (RMSE) according to

Using as data function, we obtain the computed approximation by (19). The relative error is given by
where is defined by (39).

From Figures 1, 2, 3, and 4, we find that the numerical results are quite satisfactory. Even with the noise level up to , the numerical solutions are still in good agreement with the exact solution. In addition, the numerical results are better with the increase of at first, but the numerical results are not so good after . This means that the numerical results are not so good for stronger “smoothness” assumptions on the exact solution which is consistent with the Tikhonov regularization method in [16]. These results are consistent with [8]. Comparing with [8], we can also easily find that the *a posteriori* parameter choice rule works better than the *a priori* parameter choice rule.

#### 6. Conclusions

In this paper, a truncation regularization method is used to identify the unknown source term depending only on the spatial variable for the Poisson equation. The *a posteriori* rule for choosing regularization parameter with strict theory analysis is presented. In practice, the *a priori* bound is unknown exactly, and the regularization parameters can not be obtained exactly. This is a defect of the *a priori* rule for choosing regularization parameter. Using Morozov’s discrepancy principle, we give *a posteriori* parameter choice rule which depends only on the measured data. For the *a posteriori* parameter choice rule, we obtain the Hölder type error estimate.

#### Acknowledgment

The project is supported by the Distinguished Young Scholars Fund of Lanzhou University of Technology (Q201015).

#### References

- J. R. Cannon and P. Duchateau, “Structural identification of an unknown source term in a heat equation,”
*Inverse Problems*, vol. 14, no. 3, pp. 535–551, 1998. View at Publisher · View at Google Scholar · View at Scopus - G. Li, “Data compatibility and conditional stability for an inverse source problem in the heat equation,”
*Applied Mathematics and Computation*, vol. 173, no. 1, pp. 566–581, 2006. View at Publisher · View at Google Scholar · View at Scopus - Z. Yi and D. A. Murio, “Source term identification in 1-D IHCP,”
*Computers and Mathematics with Applications*, vol. 47, no. 12, pp. 1921–1933, 2004. View at Publisher · View at Google Scholar · View at Scopus - L. Yan, F.-L. Yang, and C.-L. Fu, “A meshless method for solving an inverse spacewise-dependent heat source problem,”
*Journal of Computational Physics*, vol. 228, no. 1, pp. 123–136, 2009. View at Publisher · View at Google Scholar · View at Scopus - A. Farcas and D. Lesnic, “The boundary-element method for the determination of a heat source dependent on one variable,”
*Journal of Engineering Mathematics*, vol. 54, no. 4, pp. 375–388, 2006. View at Publisher · View at Google Scholar · View at Scopus - L. Yan, C.-L. Fu, and F.-L. Yang, “The method of fundamental solutions for the inverse heat source problem,”
*Engineering Analysis with Boundary Elements*, vol. 32, no. 3, pp. 216–222, 2008. View at Publisher · View at Google Scholar · View at Scopus - F. Yang and C.-L. Fu, “The modified regularization method for identifying the unknown source on Poisson equation,”
*Applied Mathematical Modelling*, vol. 36, no. 2, pp. 756–763, 2012. View at Publisher · View at Google Scholar · View at Scopus - F. Yang, “The truncation method for identifying an unknown source in the Poisson equation,”
*Applied Mathematics and Computation*, vol. 217, no. 22, pp. 9334–9339, 2011. View at Publisher · View at Google Scholar · View at Scopus - L. Eldén, F. Berntsson, and T. Regińska, “Wavelet and Fourier methods for solving the sideways heat equation,”
*SIAM Journal on Scientific Computing*, vol. 21, no. 6, pp. 2187–2205, 2000. View at Publisher · View at Google Scholar · View at Scopus - X.-T. Xiong, C.-L. Fu, and H.-F. Li, “Fourier regularization method of a sideways heat equation for determining surface heat flux,”
*Journal of Mathematical Analysis and Applications*, vol. 317, no. 1, pp. 331–348, 2006. View at Publisher · View at Google Scholar · View at Scopus - C.-L. Fu, X.-T. Xiong, and Z. Qian, “Fourier regularization for a backward heat equation,”
*Journal of Mathematical Analysis and Applications*, vol. 331, no. 1, pp. 472–480, 2007. View at Publisher · View at Google Scholar · View at Scopus - Z. Qian, C.-L. Fu, X.-T. Xiong, and T. Wei, “Fourier truncation method for high order numerical derivatives,”
*Applied Mathematics and Computation*, vol. 181, no. 2, pp. 940–948, 2006. View at Publisher · View at Google Scholar · View at Scopus - T. Regińska and K. Regiński, “Approximate solution of a Cauchy problem for the Helmholtz equation,”
*Inverse Problems*, vol. 22, pp. 975–989, 2006. View at Google Scholar - F. Yang and D. Li, “The fourier regularization method for identifying the unknown source on poisson equation,”
*International Journal of Mathematical Analysis*, vol. 4, no. 25–28, pp. 1357–1361, 2010. View at Google Scholar · View at Scopus - Z. Qian, C.-L. Fu, and Z.-P. Li, “Two regularization methods for a Cauchy problem for the Laplace equation,”
*Journal of Mathematical Analysis and Applications*, vol. 338, no. 1, pp. 479–489, 2008. View at Publisher · View at Google Scholar · View at Scopus - A. Kirsch,
*An Introduction to the Mathematical Theory of Inverse Problems*, Springer, New York, NY, USA, 1996. - H. W. Engl, M. Hanke, and A. Neubauer,
*Regularization of Inverse Problems*, Kluwer Academic, Boston, Mass, USA, 1996.