We consider the Cauchy problem for a modified Helmholtz equation, where the Cauchy data is given at and the solution is sought in the interval . A spectral method together with choice of regularization parameter is presented and error estimate is obtained. Combining the method of lines, we formulate regularized solutions which are stably convergent to the exact solutions.
The Cauchy problem for Helmholtz equation arises from inverse scattering problems. Specific backgrounds can be seen in the existing literature; we can refer to [1–6] and so forth. A number of numerical methods for stabilizing this problem are developed. Several boundary element methods combined with iterative, conjugate gradient, Tikhonov regularization, and singular value decomposition methods are compared in . Cauchy problem for elliptic equations is well known to be severely ill-posed ; that is, the solution does not depend continuously on the boundary data, and small errors in the boundary data can amplify the numerical solution infinitely; hence it is impossible to solve Cauchy problem of Helmholtz equation by using classical numerical methods and it requires special techniques, for example, regularization methods. Although theoretical concepts and computational implementation related to the Cauchy problem of Helmholtz equation have been discussed by many authors [8–11], there are many open problems deserved to be solved. For example, many authors have considered the following Cauchy problem of the Helmholtz equation [8–11]:
However, the boundary condition is very strict. If the boundary condition is replaced by , where is a function of , then their methods cannot be applied easily. Therefore, in this paper, we consider the following Cauchy problem for the Helmholtz equation:
We want to seek the solution in the interval from the Cauchy data pairs located at . Of course, since and are assumed to be measured, there must be measurement errors, and we would actually have noisy data function , for which the measurement errors and are small. Here and in the following sections, denotes the norm. Thus (1.2) is a noncharacteristic Cauchy problem with appropriate Cauchy data given on the line .
However, the ill-posedness is caused by high frequency. By introducing a "cutoff'' frequency we can obtain a well-posed problem. This method has been used for solving inverse heat conduction problem , and sideways heat equation , sideways parabolic equation . An error estimate for the proposed method can be found in Section 2. The implementation of the numerical method is explained.
2. Regularization and Error Estimate
First, let be the Fourier transform of :
Taking Fourier transformation for (1.2), we have a family of problems parameterized by :
The solution can easily be verified to be
Following the idea of Fourier method, we consider (2.2) only for by cutting off high frequency and define a regularized solution:
where is the characteristic function of the interval . The solution can be found by using the inverse Fourier transform. Define the regularized solution with measured data by . The difference between the exact solution and the regularized solution can be divided into two parts:
We rewrite (2.2) as a system of ordinary differential equations:
Letting , we can rewrite (2.7) as
where the matrix is the one in (2.7). The reason for using instead of in the definition of is that with this choice, the matrix is normal and hence diagonalized by a unitary matrix. The eigenvalues of are . Thus we can factorize as
where is a unitary matrix, it follows that the solution of the system (2.8) can be written as
Since is unitary, and therefore
where denotes both the Euclidean norm in the complex vector space and the subordinate matrix norm, remembering in (2.11), we obtain
This inequality is valid for all and we can integrate over and use the Parseval theorem to obtain estimates for and in the -norm. First we will prove a bound on the difference between any two regularized solutions (2.4). We have errors in the measured and . These two cases are treated separately.
Assume that one has two regularized solutions and defined by (2.4), with the Cauchy data and . Then
The function satisfies the differential equation; thus solves (2.2), with initial data given by
and thus satisfies inequality(2.12). We have
If , then
holds. By inserting (2.15) into (2.16) we get
If we insert and integrate over , then we obtain
Thus (2.13) holds.
Assume that one has two regularized solutions and defined by (2.4), with the Cauchy data and . Then
The function satisfies the differential equation, with initial data given by
and thus inequality (2.12) holds. It follows that
and thus, for ,
By integrating over the interval , we get
Since is equal to zero outside the interval , we can extend the integrals. Inserting we get
This is precisely (2.19).
Next we prove that, for the regularized problem, we have stable dependence of the data. By using the two previous lemmas, we get the following.
Lemma 2.3 (stability).
Assume that is the regularized solution (2.4), with exact data and is the regularized solution with noisy data ; then
Let be a regularized solution defined by (2.4), with data . Then by Lemma 2.1,
Using Lemma 2.2, we have
By the triangle inequality,
This completes the proof.
By Lemma 2.3 the regularized solution depends continuously on the data. Next we derive a bound on the truncation error when we neglect frequencies in (2.3). So far we have not used any information about for we, assume that the Helmholtz equation is valid in a large interval . By imposing a priori bounds on the solution at and at , we obtain an estimate of the difference between the exact solution and a regularized solution with "cutoff'' level . This is a convergence result in the sense that as for the case of exact data . The following estimate holds.
Lemma 2.4 (convergence).
Suppose that is the solution of the problem (2.2), and that the Helmholtz equation is valid for . Then the difference between and a regularized solution can be estimated:
where the constant is defined by
The solution of (2.2) can be written in the floowing form:
where and can be determined from the boundary conditions:
Solving for and we find that the solution can be written:
We make the observation that
Using the expression (2.33) and the triangle inequality, we get
The first term on the right-hand side satisfies
Similarly, the second term can be estimated:
By combining these two expressions, we have
Thus the proof is complete.
The constant is well defined, but its value has to be estimated. From a numerical computation we conclude that .
When solving (2.2) numerically we need Cauchy data , along the line . The most natural way to obtain this is to use two thermocouples, located at , and , and compute by solving a well-posed problem for the Helmholtz equation in the interval . Hence it is natural to assume knowledge about at a second point.
Let us summarize what we have so far. The constant in Lemma 2.4 is unchanged. The propagated data error is estimated using Lemma 2.3:
and the truncation error is estimated using Lemma 2.4:
These two results can be combined into an error estimate for the spectral method. This is demonstrated in two examples.
Suppose that , and that we have an estimate of the noise level, . If and if we choose , with this choice, then by expression (2.42),
where we have assume that and used the bound . By expression (2.41),
Thus we obtain an error estimate:
Note that, under these assumptions, can be used as a rule for selecting the regularization parameter.
Suppose that the Helmholtz equation is valid in the interval and that we have a priori bounds and . Furthermore, we assume that the measured data satisfies and . Then we have the estimates:
By balancing these two components and , we can find a suitable value for the regularization parameter.
3. Numerical Implementation
In this section, we use the method of Lines. Let
and then we can obtain
The matrix approximates the second-order derivative by spectral cutoff method (see ). This ordinary differential equation system can easily be solved by various numerical methods.
It is easy to verify that the function is the exact solution of (1.2). Now we need to seek the solution , where from the Cauchy data pairs and .
In our numerical experiment, we give the comparison between the exact solution of problem (2.2) and its approximation by spectral cutoff method for different noise levels and different locations . Please see Figures 1, 2, 3, and 4. In the experiment, the regularization parameter is chosen according to the Remark of Remark 2.5, where we take . From the above results, it is easy to see that the numerical effect of the spectral method works well. Moreover, we can also see that the lower the noise level is, the better the approximate effect is; the closer to the boundary the location is, the better the approximate effect is. These accord with the theory in Section 2.
The author wants to express his thanks to the referee for many valuable comments. This work is supported by the Educational Commission of Hubei Province of China (Q20102804, T201009).
Harari, I, Barbone, PE, Slavutin, M, Shalom, R: Boundary infinite elements for the Helmholtz equation in exterior domains. International Journal for Numerical Methods in Engineering. 41, 1105–1131 (1998). Publisher Full Text
Marin, L, Elliott, L, Heggs, PJ, Ingham, DB, Lesnic, D, Wen, X: An alternating iterative algorithm for the Cauchy problem associated to the Helmholtz equation. Computer Methods in Applied Mechanics and Engineering. 192(5-6), 709–722 (2003). Publisher Full Text
Marin, L, Elliott, L, Heggs, PJ, Ingham, DB, Lesnic, D, Wen, X: Conjugate gradient-boundary element solution to the Cauchy problem for Helmholtz-type equations. Computational Mechanics. 31(3-4), 367–377 (2003)
Regińska, T, Regiński, K: Approximate solution of a Cauchy problem for the Helmholtz equation. Inverse Problems. 22(3), 975–989 (2006). Publisher Full Text
Marin, L, Elliott, L, Heggs, PJ, Ingham, DB, Lesnic, D, Wen, X: Comparison of regularization methods for solving the Cauchy problem associated with the Helmholtz equation. International Journal for Numerical Methods in Engineering. 60(11), 1933–1947 (2004). Publisher Full Text
Xiong, XT, Fu, CL: Two approximate methods of a Cauchy problem for the Helmholtz equation. Computational & Applied Mathematics. 26(2), 285–307 (2007). PubMed Abstract | Publisher Full Text | PubMed Central Full Text
Fu, CL, Feng, XL, Qian, Z: The Fourier regularization for solving the Cauchy problem for the Helmholtz equation. Applied Numerical Mathematics. 59(10), 2625–2640 (2009). Publisher Full Text
Qin, HH, Wei, T: Modified regularization method for the Cauchy problem of the Helmholtz equation. Applied Mathematical Modelling. 33(5), 2334–2348 (2009). Publisher Full Text
Qian, AL, Xiong, XT, Wu, YJ: On a quasi-reversibility regularization method for a Cauchy problem of the Helmholtz equation. Journal of Computational and Applied Mathematics. 233(8), 1969–1979 (2010). Publisher Full Text
Berntsson, F: A spectral method for solving the sideways heat equation. Inverse Problems. 15(4), 891–906 (1999). Publisher Full Text
Xiong, XT, Fu, CL, Cheng, J: Spectral regularization methods for solving a sideways parabolic equation within the framework of regularization theory. Mathematics and Computers in Simulation. 79(5), 1668–1678 (2009). Publisher Full Text