1. Introduction
The inverse Sturm–Liouville theory was originated in 1929 by Ambarzumian [
1] and further developed in [
2,
3,
4,
5,
6].
The
th order, nonsingular, self-adjoint eigenvalue problem (EVP) or Sturm–Liouville problem (SLP) is given by:
together with some boundary conditions at
a and
b, the functions
, and
being continuous on the finite closed interval
, and
having a continuous derivative. In the inverse SLP, the coefficient functions
,
need to be reconstructed, given suitable valid spectral data.
For a discussion of analytical methods and numerical methods for inverse SLP, see [
7,
8], respectively. Iterative methods [
9,
10], Rayleigh–Ritz method [
11], finite difference approximation [
12], Quasi-Newton method [
13], shooting method [
14], interval Newton’s method [
15], finite-difference method [
16], boundary value methods [
17,
18,
19,
20], Numerov’s method [
21,
22,
23], least-squares functional [
24], generalized Rundell–Sacks algorithm [
25,
26], spectral mappings [
27], Lie-group estimation method [
28], Broyden method [
29,
30], decent flow methods [
31], modified Numerov’s method [
32], Newton-type method [
33], Fourier–Legendre series [
34], and Chebyshev polynomials [
35] are of particular importance among the existing methods to solve inverse SLP.
Numerical algorithms to solve the inverse fourth-order Sturm–Liouville problem (FSLP) are proposed in [
36,
37,
38].
In [
39], authors implemented Barcilon [
9]’s iterative algorithm to solve inverse SLPs successfully. The present paper extends the results for the case
by considering the general-order inverse SLP algorithm as described in [
36].
2. Materials and Methods
Inverse SLP Algorithm for the General Order
Setting
in (
1) and taking the domain
, the EVP reads as:
According to [
40],
distinct spectra are required in order to determine
. In other words,
different EVPs with distinct boundary conditions (BC) are needed. In addition, the spectra need to be interlaced, i.e.,
Assume that m boundary conditions are common to all the given eigenvalue problems, say, conditions at the right end point .
Combine the eigenvalue problems into an equivalent linear differential equation:
- 1.1.
Define a vector function
:
where
and
are two solutions of Equation (
2).
- 1.2.
Differentiate each , times.
- 1.3.
Obtain linear equations for , where , .
- 1.4.
Express , , , in terms of and their first derivatives.
- 1.5.
Differentiate
,
again, to obtain
m linear coupled differential equations of order
, which is in matrix form:
Obtain the corresponding boundary conditions for the linear differential equation using the boundary conditions of the eigenvalue problems:
- 2.1.
Assume that and satisfy the m boundary conditions common to all the given eigenvalue problems.
- 2.2.
Then from 1.3, find the boundary conditions at , and the boundary conditions at for .
Solve Equation (
4) using above boundary conditions, to obtain the solutions
.
Solve the adjoint system of equations to the above system, and denote the solutions by .
Find the bi-orthogonal set of functions to , using the relation .
Some remarks on Algorithm 1 follow:
- Line (Data:)
Infinite eigenvalue (EV) sequence is replaced with the first
N EVs. The truncated EV sequences
,
are obtained, using the Magnus method [
39] to solve Equation (
2).
- Line (Result:)
Output : matrix of function values, n: the number of subdivisions in the x-axis and m: the number of potentials.
- Line (1)
We take initial guess as .
- Line (3)
while loop exit condition is replaced with a maximum number of iterations: and in each iteration the condition (: a fixed tolerance) is checked to exit the loop.
- Line (4)
denotes the integration w.r.t.
x from 0 to 1 and
.
and
are kept fixed at
and
, which are the solutions by setting
,
. Solutions
and
are calculated using Mathematica [
41]. Matlab [
42] built-in functions:
trapz() and
griddedInterpolant() with
pchip (piecewise cubic Hermite interpolating polynomial) option are used to approximate integrals and to have
p as a function, respectively. The latter is required since,
p as a function is input to the Magnus method. Again, the eigenvalue sequence
is calculated using
and the iteration repeats.
You may refer Section 2.6 of [
39] for the specific computational details of the Algorithm 1 for the case
and next section for the case
.
Algorithm 1: Solving inverse Sturm–Liouville problem (SLP). |
|
3. Inverse FSLP Algorithm
Setting
in Equation (
2) we obtain FSLP:
Three spectra
,
and
required to reconstruct the unknown coefficients:
p and
q will be obtained from the following EVPs:
Although, Equations (
9) and (
10) are self-adjoint, Equation (
11) is not. Its adjoint equation is:
The eigenvalues
,
and
(for
) are given by the solutions of the following equations:
Furthermore,
so that the spectra
,
and
are interlaced, i.e.,
Since the eigenvalues are interlaced, they can be coalesced into a single one by defining
Let
and
be two solutions of (
8) satisfying the common boundary conditions of Equations (
9)–(
12), namely, (
, etc.)
By differentiating
and
five times each, and eliminating
, and
, 9th order linear differential equation for
(or a coupled system of linear differential equations of order 5 for
and
) can be obtained. When
, this reduces to
or
in the form of the Equation (
4), and the boundary conditions read as:
The solution of the system of Equations (
21)–(
23) using Mathematica is:
with
. The adjoint of system of Equations (
21)–(
23) is
with
The bi-orthogonal set for
given by
is:
4. Results for Direct Singular and Inverse SLPs
Here, we present some numerical examples of direct and inverse SLPs of orders 2 and 4.
Algorithm 1 was implemented using MATLAB (2014) [
42]. The reference solutions of the EVPs are computed using Wolfram Mathematica 11 [
41]. Direct problems were solved using the Magnus expansion method [
39].
For the numerical calculations, n is the number of subdivisions in the interval , m is the number of subdivisions in the interval ; being the maximum eigenvalue searching, L is the number of multisection steps used to calculate each eigenvalue in the characteristic function and M is the number of inverse algorithm steps.
Error values are stated using max-norm which is defined as .
4.1. Singular SLPs
Example 1 (Singular SLP: infinite domain)
. This example extends the Magnus method, as described in [39] to the infinite domain. Consider the harmonic oscillator problemtogether with Dirichlet and Neumann BCs, having exact eigenvalues , [43].The domain is truncted to , and using the parameters , , the first three eigenvalues are computed (see Table 1). It can be deduced that the Magnus method has comparable performance to other methods. Moreover, fixed parameters can be changed to obtain better accuracy. Example 2 (Singular SLP with asymmetric potential)
. Consider the Morse potential problemtogether with Dirichlet and Neumann BCs. When this has a finite number of EV: , with , , where [43].The domain is truncated to , and using the parameters , , and , the first three eigenvalues are computed (see Table 2). Magnus method is more accurate when the truncated interval is small. Furthermore, it can be observed that initial eigenvalues are more accurate than the later ones. Example 3 (Singular SLP in half axis)
. Consider the hydrogen atom equationtogether with Dirichlet and Neumann BCs. This has a finite number of eigenvalues lying between , given by , with , [44].The domain is truncated to , and using the parameters , , the first few eigenvalues are computed (see Table 3). 4.2. Inverse SLPs
Example 4 (Inverse SLP with a different domain)
. This example extends Barcilon’s algorithm in [36] by changing the domain to :Here, is symmetric and normalized. In the inverse Algorithm 1 of [39], and need to be changed into , which are the basic solutions in the new domain. The truncated eigenvalue sequences at , for Dirichlet boundary conditions (DDBCs) and Dirichlet– Neumann boundary conditions (DNBCs) are calculated using the Magnus method (with , , ) and inverse algorithm steps are used. Figure 1 shows the reconstructed and exact potential, and the log absolute errors in the reconstructed potential, respectively. From Figure 1a, it is obvious that the potential is converging towards the exact one. The max-norm of difference between the reconstructed p and the actual p is ≈ and the max-norm of difference between the reconstructed eigenvalues and the actual ones is ≈. Figure 2 shows reconstructed p starting with perturbed eigenvalueswhere δ is the noise level and , , , , . It is obvious that reconstruction is possible even in the presence of a significant noise.
Example 5 (Inverse SLP with a non-smooth potential)
. This example extends Barcilon’s algorithm in [36] by reconstructing a non-smooth potential:with : symmetric and normalized, but non-smooth. The truncated eigenvalue sequences at , for DDBCs and DNBCs are calculated using the Magnus method (with , , ), and used in the reconstruction of the potential.Figure 3 shows the reconstructed and exact potential, and the log absolute errors in the reconstructed potential, respectively, using one inverse algorithm step. From Figure 3a, it is obvious that the potential is converging towards the exact one. The max-norm of difference between the reconstructed p and the actual p is ≈ and the max-norm of difference between the reconstructed eigenvalues and the actual ones is ≈. As anticipated, the error is maximum at the point of non-differentiability and at the boundaries. Example 6 (FSLP)
. Considerwith and . Using the three spectra corresponding to the sets of BCs: , and two potential functions , are reconstructed starting with zero initial guesses and , , (Figure 4). The max-norm of difference between the reconstructed and actual is ≈ and the max-norm of difference between the reconstructed eigenvalues and the actual ones is ≈. Example 7 (FSLP)
. Considerwith and . Using the three spectra corresponding to the sets of BCs: , and two potential functions , are reconstructed starting with zero initial guesses and , , (Figure 5). The max-norm of difference between the reconstructed and the actual is ≈ and the max-norm of difference between the reconstructed eigenvalues and the actual ones is ≈. 5. Discussion
In the first few examples, the Magnus method [
39] is further extended to solve some direct singular SLPs on infinite domain, half-axis and for an asymmetric potential. Later, Barcilon’s algorithm is extended to solve inverse SLPs in a non–unit domain and for a non–smooth potential. The first two examples, which extend the Magnus method in [
39] to infinite interval, illustrates it’s applicability to singular SLPs effectively. Furthermore, the present method does not require transforming the singular SLP to a regular one or re-scaling the eigenvalues. It effectively finds the required eigenvalues, just by truncating the interval to a finite one.
Next, the inverse SLP is solved with smooth and non-smooth potential, even in the presence of noise. It is observed that the method is successful even in the presence of significant noise, provided that the assumptions of the algorithm are satisfied.
The last few examples solve the inverse FSLP using Barcilon’s algorithm with the initial knowledge of three spectra. A simplified FSLP is solved keeping one of the unknown potential functions zero.
6. Conclusions
More than 40 years after Barcilon’s paper [
36], this paper gives a concrete implementation of the inverse SLP algorithm proposed therein, to our knowledge for the first time. Furthermore, computational feasibility and applicability of this algorithm for solving inverse SLPs of higher order is verified successfully in this paper (for
and
).
In future work, it is possible to extend the results to solve even higher order SLPs.
In conclusion, this work provides a method that can be adapted successfully for solving a direct (regular/singular) or inverse SLP of an arbitrary order with arbitrary BCs.
Author Contributions
All authors contributed equally. All authors have read and agreed to the published version of the manuscript.
Funding
This research was partially funded by UNIVERSITY OF KELANIYA, SRI LANKA.
Conflicts of Interest
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.
References
- Ambarzumian, V. Über eine frage der eigenwerttheorie. Z. Phys. A Hadron. Nucl. 1929, 53, 690–695. [Google Scholar]
- Tikhonov, A. On the uniqueness of the solution of the problem of electromagnetic sounding. Dokl. Akad. Nauk SSSR 1949, 69, 797–800. [Google Scholar]
- Levinson, N. The Inverse Sturm-Liouville Problem. Mat. Tidsskr. B 1949, 1949, 25–30. [Google Scholar]
- Marčenko, V. Sturm-Liouville Operators and Their Applications; Birkhauser Verlag: Basel, Switzerland, 1986. [Google Scholar]
- Levitan, B. Inverse Sturm-Liouville Problems; VNU Science Press: Utrecht, The Netherlands, 1987. [Google Scholar]
- Kravchenko, V.V. Direct and Inverse Sturm-Liouville Problems: A Method of Solution; Springer Nature: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
- McLaughlin, J.R. Analytical methods for recovering coefficients in differential equations from spectral data. SIAM Rev. 1986, 28, 53–72. [Google Scholar] [CrossRef]
- Freiling, G.; Yurko, V.A. Inverse Sturm-Liouville Problems and Their Applications; NOVA Science Publishers: New York, NY, USA, 2001. [Google Scholar]
- Barcilon, V. Iterative solution of the inverse Sturm-Liouville problem. J. Math. Phys. 1974, 15, 429–436. [Google Scholar] [CrossRef]
- Sacks, P.E. An iterative method for the inverse Dirichlet problem. Inverse Probl. 1988, 4, 1055–1069. [Google Scholar] [CrossRef]
- Hald, O.H. The inverse Sturm-Liouville problem and the Rayleigh-Ritz method. Math. Comput. 1978, 32, 687–705. [Google Scholar] [CrossRef]
- Paine, J. A numerical method for the inverse Sturm-Liouville problem. SIAM J. Sci. Stat. Comput. 1984, 5, 149–156. [Google Scholar] [CrossRef]
- Rundell, W.; Sacks, P.E. Reconstruction techniques for classical inverse Sturm-Liouville problems. Math. Comput. 1992, 58, 161–183. [Google Scholar] [CrossRef]
- Lowe, B.D.; Pilant, M.; Rundell, W. The recovery of potentials from finite spectral data. SIAM J. Math. Anal. 1992, 23, 482–504. [Google Scholar] [CrossRef]
- Neher, M. Enclosing solutions of an inverse Sturm-Liouville problem with finite data. Computing 1994, 53, 379–395. [Google Scholar] [CrossRef] [Green Version]
- Fabiano, R.H.; Knobel, R.; Lowe, B.D. A finite-difference algorithm for an inverse Sturm-Liouville problem. IMA J. Numer. Anal. 1995, 15, 75–88. [Google Scholar] [CrossRef]
- Ghelardoni, P. Approximations of Sturm-Liouville eigenvalues using boundary value methods. Appl. Numer. Math. 1997, 23, 311–325. [Google Scholar] [CrossRef]
- Aceto, L.; Ghelardoni, P.; Magherini, C. BVMs for Sturm-Liouville eigenvalue estimates with general boundary conditions. J. Numer. Anal. Ind. Appl. Math 2009, 4, 113–127. [Google Scholar]
- Ghelardoni, P.; Magherini, C. BVMs for computing Sturm-Liouville symmetric potentials. Appl. Math. Comput. 2010, 217, 3032–3045. [Google Scholar] [CrossRef]
- Aceto, L.; Ghelardoni, P.; Magherini, C. Boundary value methods for the reconstruction of Sturm-Liouville potentials. Appl. Math. Comput. 2012, 219, 2960–2974. [Google Scholar] [CrossRef]
- Andrew, A.L. Numerov’s method for inverse Sturm-Liouville problems. Inverse Probl. 2005, 21, 223–238. [Google Scholar] [CrossRef]
- Andrew, A.L. Numerical solution of inverse Sturm-Liouville problems. ANZIAM J. 2004, 45, 326–337. [Google Scholar] [CrossRef]
- Andrew, A.L. Computing Sturm-Liouville potentials from two spectra. Inverse Probl. 2006, 22, 2069–2081. [Google Scholar] [CrossRef]
- Röhrl, N. A least-squares functional for solving inverse Sturm–Liouville problems. Inverse Probl. 2005, 21, 2009–2017. [Google Scholar] [CrossRef] [Green Version]
- Rafler, M.; Böckmann, C. Reconstruction method for inverse Sturm–Liouville problems with discontinuous potentials. Inverse Probl. 2007, 23, 933–946. [Google Scholar] [CrossRef]
- Efremova, L.; Freiling, G. Numerical solution of inverse spectral problems for Sturm-Liouville operators with discontinuous potentials. Open Math. 2013, 11, 2044–2051. [Google Scholar] [CrossRef]
- Ignatiev, M.; Yurko, V. Numerical Methods for Solving Inverse Sturm-Liouville Problems. Results Math. 2008, 52, 63–74. [Google Scholar] [CrossRef]
- Liu, C.S. Solving an inverse Sturm-Liouville problem by a Lie-group method. Bound. Value Probl. 2008, 2008, 749865. [Google Scholar] [CrossRef] [Green Version]
- Böckmann, C.; Kammanee, A. Broyden method for inverse non-symmetric Sturm-Liouville problems. BIT Numer. Math. 2011, 51, 513–528. [Google Scholar] [CrossRef]
- Böckmann, C.; Rattana, A. An inverse fourth order sturm-liouville problem. AIP Conf. Proc. 2015, 1648, 850030. [Google Scholar] [CrossRef]
- Gao, Q. Decent flow methods for inverse Sturm–Liouville problem. Appl. Math. Model. 2012, 36, 4452–4465. [Google Scholar] [CrossRef]
- Gao, Q.; Cheng, X.; Huang, Z. Modified Numerov’s method for inverse Sturm–Liouville problems. J. Comput. Appl. Math. 2013, 253, 181–199. [Google Scholar] [CrossRef]
- Drignei, M.C. A Newton-type method for solving an inverse Sturm-Liouville problem. Inverse Probl. Sci. Eng. 2014, 23, 851–883. [Google Scholar] [CrossRef]
- Kravchenko, V.V. On a method for solving the inverse Sturm-Liouville problem. J. Inverse Ill Posed Probl. 2018. [Google Scholar] [CrossRef]
- Neamaty, A.; Akbarpoor, S.; Yilmaz, E. Solving Symmetric Inverse Sturm-Liouville Problem Using Chebyshev Polynomials. Mediterr. J. Math. 2019, 16, 74. [Google Scholar] [CrossRef]
- Barcilon, V. On the solution of inverse eigenvalue problems of high orders. Geophys. J. Int. 1974, 39, 143–154. [Google Scholar] [CrossRef] [Green Version]
- McLaughlin, J.R. An Inverse Eigenvalue Problem of Order Four. SIAM J. Math. Anal. 1976, 7, 646–661. [Google Scholar] [CrossRef]
- McLaughlin, J.R. An Inverse Eigenvalue Problem of Order Four—An Infinite Case. SIAM J. Math. Anal. 1978, 9, 395–413. [Google Scholar] [CrossRef]
- Perera, U.; Böckmann, C. Solutions of Direct and Inverse Even-Order Sturm-Liouville Problems Using Magnus Expansion. Mathematics 2019, 7, 544. [Google Scholar] [CrossRef] [Green Version]
- Barcilon, V. On the uniqueness of inverse eigenvalue problems. Geophys. J. R. Astron. Soc. 1974, 38, 287–298. [Google Scholar] [CrossRef] [Green Version]
- Mathematica; Wolfram Research, Inc.: Champaign, IL, USA, 2016.
- Matlab; The MathWorks, Inc.: Natick, MA, USA, 2014.
- Taşeli, H. Accurate numerical bounds for the spectral points of singular Sturm–Liouville problems over −∞ < x < ∞. J. Comput. Appl. Math. 2000, 115, 535–546. [Google Scholar] [CrossRef]
- Marletta, M.; Pryce, J.D. Automatic solution of Sturm-Liouville problems using the Pruess method. J. Comput. Appl. Math. 1992, 39, 57–78. [Google Scholar] [CrossRef] [Green Version]
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).