Eigenvalues of Vectorial Sturm-Liouville Problems with Parameter Dependent Boundary Conditions
We generalize the regularized sampling method introduced in 2005 by the author to compute the eigenvalues of scalar Sturm-Liouville problems (SLPs) to the case of vectorial SLP with parameter dependent boundary conditions. A few problems are worked out to illustrate the effectiveness of the method and show by the same token that we have indeed a general method capable of handling with ease very broad classes of SLPs, whether scalar or vectorial.
In  we introduced the regularized sampling method, a method to compute the eigenvalues of scalar Sturm-Liouville problems (SLPs) with parameter dependent boundary conditions. We subsequently used this method to compute the eigenvalues of singular and non-self-adjoint Sturm-Liouville problems. The scope of the method was further extended to include the computation of the eigenvalues of discontinuous/impulsive, nonlocal ( and the references therein), and two-parameter SLPs . Continuing our effort we will tackle in this paper vectorial SLP with parameter dependent nonseparated boundary conditions. Vectorial Sturm-Liouville problems have been considered in [4–13] and the references therein while corresponding inverse problems appeared in [14–17].
2. The Characteristic Function
Consider the vectorial Sturm-Liouville problem, where is an matrix function, and are real matrix functions of the parameter such that the matrix has full rank.
Let , be the solutions of the Sturm-Liouville matrix equation subject to the initial conditions , and , , respectively, being the identity matrix and being the zero matrix.
The general solution of the Sturm-Liouville equation is given by with arbitrary constant vectors and . Replacing in the boundary conditions, we get To have a nontrivial solution a necessary and sufficient condition is that where the characteristic function is where
The eigenvalues of (1) are the square of the zeroes of . It is well known that the multiplicities of these eigenvalues are at most .
3. Main Results
Let be the Paley-Wiener space and recall the celebrated Whittaker-Shannon-Kotel’nikov theorem .
Theorem 1. Let ; then where the series converges uniformly on compact subset of and in .
It is known that, in the case of scalar Sturm-Liouville problems, is an entire function of for each fixed . is in a Paley-Wiener space as a function of for each only in the Dirichlet case. So, we had to subtract some terms from to make the difference fall in an appropriate space. We had even to subtract terms involving multiple integrals to get sharper results when it comes to computing of the eigenvalues. The regularized sampling method has been introduced recently  to overcome this problem; we do not have to subtract any term involving any (multiple) integration. In fact we multiplied and by an appropriate function of and got the eigenvalues with much greater precision at a reduced cost. Here and are known simple functions.
For the vectorial Sturm-Liouville problem at hand, we will use the regularized sampling method to recover the matrices , , , and from which we obtain , the characteristic function whose zeroes are the square roots of the sought eigenvalues of the problem.
Consider the compatible vector and matrix norms given by where , . In the following we will make use of the standard estimate.
Lemma 2. Consider where is some constant (we may take ).
To cover both cases (, and , ) we will consider the following initial value problem: where and are matrices or -vector. We have
Our first result is the following theorem.
Theorem 3. is an entire matrix function of for each fixed and satisfies the growth conditions for some positive constants , , , and .
Proof. From (11) and using standard arguments, we conclude that is an entire matrix function of for each in . Its derivative with respect to , is also an entire matrix function of for each in . Going back to (11) we get at once Multiplying by , using Gronwall’s lemma, and multiplying back by we get where . Now, using the above estimate in (10), we get where . Likewise we have where . As in the scalar case, is in a Paley-Wiener space only in the Dirichlet case; however, is. As for , it is not; nor is since they are not square integrable over the reals for fixed in . Also, where .
We get at once the following corollaries.
Corollary 4. , , , , , and are entire matrix functions of for each fixed and satisfy the growth conditions for some generic positive constants , , , and .
Corollary 5. The functions, belong to the Paley-Wiener space as functions of and thus can be recovered from their samples at , using the WSK series.
Theorem 6. Let be positive real number and a positive integer. Consider belong to the Paley-Wiener space where + as functions of for each fixed for and satisfy the growth condition , where is some positive constant ().
Proof. It is enough to note that is an entire function of and satisfies the estimate in the above Lemma and the fact that is the product of two entire functions thus entire.
Remark 7. To avoid the first singularity of we will take .
The use of the WSK theorem allows us to recover and as where , , , and .
Hence, or can be recovered as In practice, we take for some positive integer , large enough, so that can be reconstructed whose zeros are the square roots of the sought eigenvalues.
Since and are in , Jagerman’s result  is applicable and yields the following better estimate.
Lemma 8 (truncation error). Let denote the truncation of , . Then, for , where .
Lemma 9. Consider , where .
The approximation of and by and , respectively, induces an approximation of the characteristic function by , whose zeros are the square root of the eigenvalues of the problem.
Let denote an eigenvalue of the problem; then independent eigenfunctions associated can be obtained using basis vectors of the null space of the matrix as initial conditions to the differential equation , .
4. Numerical Examples
In this section we will illustrate the power of the regularized sampling method as applied to vectorial Sturm-Liouville problems with parameter dependent boundary conditions. We will take , and a precision of for the first three examples involving two dimensional SLPs. We will also work out two three-dimensional SLPs one of them involving parameter dependent boundary conditions. In these last two examples we take different values of , namely, , and , and take and a precision of . The reported multiplicities of the eigenvalues are just the dimensions of the null space of the corresponding matrices .
Example 1 (Chanane , 1D-version taken from fom Binding and Browne ). Consider where . The first three eigenvalues were obtained as , , and putting them at about from the exact eigenvalues. All these are double eigenvalues.
Example 2. Consider
The first four eigenvalues were obtained as , , , and . Their multiplicity is two. Figure 1 illustrates the graph of the characteristic function.
In the next example we change the boundary conditions in Example 2 to a parameter dependent one.
Example 3. Consider
The first ten eigenvalues were obtained as , , , , , , , , , and . All these are simple eigenvalues. Figure 2 illustrates the graph of the characteristic function.
Next we consider three-dimensional vectorial SLPs, with different boundary conditions.
Example 4. Consider
Here, we will take , , and , and a precision of . Figure 3 illustrates the characteristic function over the range , while Figures 4 and 5 zoom into the regions containing the eigenvalues. Note that, in the range of interest , the graphs of , , and , are on the top of each other. A precision on can be obtained with just . It appears clearly that in this example we have a simple eigenvalue and a double eigenvalue , followed by a simple eigenvalue and a double eigenvalue (Figures 7, 8, 9, and 10). To obtain the double eigenvalues we look for the roots of and then evaluate which happened to be in each case of the order of . Table 1 illustrates these as function of , the number of sampling points.
The first few coefficients in the cardinal series expansion of are given as follows: The above data have been reported with only a few digits.
Example 5. Consider
Here, we will take , , and , and a precision of . Figure 6 illustrates the characteristic function over the range . In this range, the graphs of , , and , are on the top of each other. In this example the first seven (07) eigenvalues , , are all simple. Tables 2(a) and 2(b) illustrate these as function of , the number of sampling points.
In this paper we have extended the domain of application of the regularized sampling method to the case of vectorial Sturm-Liouville problems with parameter dependent boundary conditions. We have presented the theoretical foundation of the method and worked out a few examples to illustrate the method and shown by the same token that we have indeed a general method capable of handling with ease very broad classes of SLPs, whether scalar or vectorial, and providing the results at a reduced cost.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
The author thanks KFUPM for its continuous support, the referees for their valuable comments, and the editor for handling the paper.
B. Chanane, “Computation of the eigenvalues of Sturm-Liouville problems with parameter dependent boundary conditions using the regularized sampling method,” Mathematics of Computation, vol. 74, no. 252, pp. 1793–1801, 2005.View at: Publisher Site | Google Scholar | MathSciNet
B. Chanane, “Computing the eigenvalues of a class of nonlocal Sturm-Liouville problems,” Mathematical and Computer Modelling, vol. 50, no. 1-2, pp. 225–232, 2009.View at: Publisher Site | Google Scholar | MathSciNet
B. Chanane and A. Boucherif, “Computation of the eigenpairs of two-parameter Sturm-Liouville problems using the Regularized Sampling Method,” Abstract and Applied Analysis, vol. 2014, Article ID 695303, 6 pages, 2014.View at: Publisher Site | Google Scholar | MathSciNet
J. D. Pryce, “Classical and vector sturm—liouville problems: recent advances in singular-point analysis and shooting-type algorithms,” Journal of Computational and Applied Mathematics, vol. 50, no. 1–3, pp. 455–470, 1994, Proceedings of the Fifth International Congress on Computational and Applied Mathematics (Leuven, 1992).View at: Google Scholar
M. Marletta, “Automatic solution of regular and singular vector Sturm-Liouville problems,” Numerical Algorithms, vol. 4, no. 1-2, pp. 65–99, 1993.View at: Publisher Site | Google Scholar | MathSciNet
R. Carlson, “Large eigenvalues and trace formulas for matrix Sturm-Liouville problems,” SIAM Journal on Mathematical Analysis, vol. 30, no. 5, pp. 949–962, 1999.View at: Publisher Site | Google Scholar | MathSciNet
J. M. Calvert and W. D. Davison, “Oscillation theory and computational procedures for matrix Sturm-Liouville eigenvalue problems, with an application to the hydrogen molecular ion,” Journal of Physics A: Mathematical and General, vol. 2, no. 3, pp. 278–292, 1969.View at: Publisher Site | Google Scholar
B. I. Bandyrskii, I. P. Gavrilyuk, I. I. Lazurchak, and V. L. Makarov, “Functional-discrete method (FDmethod) for matrix Sturm-Liouville problems,” Computational Methods in Applied Mathematics, vol. 5, no. 4, pp. 362–386, 2005.View at: Google Scholar
H. Volkmer, “Matrix Riccati equations and matrix Sturm-Liouville problems,” Journal of Differential Equations, vol. 197, no. 1, pp. 26–44, 2004.View at: Publisher Site | Google Scholar | MathSciNet
H. I. Dwyer and A. Zettl, “Eigenvalue computations for regular matrix Sturm-Liouville problems,” Electronic Journal of Differential Equations, no. 5, pp. 1–13, 1995.View at: Google Scholar
M. Jodeit Jr. and B. M. Levitan, “Isospectral vector-valued sturm-liouville problems,” Letters in Mathematical Physics, vol. 43, no. 2, pp. 117–122, 1998.View at: Publisher Site | Google Scholar
D. P. John, Numerical Solution of Sturm-Liouville Problems, Monographs on Numerical Analysis, Oxford Science Publications, The Clarendon Press, Oxford University Press, New York, NY, USA, 1993.
I. L. Makarov, “Exact and truncated difference schemes for the vector Sturm-Liouville problem,” Vychislitel'naya i Prikladnaya Matematika, no. 49, pp. 72–84, 1983 (Russian).View at: Google Scholar
C.-L. Shen, “Some inverse spectral problems for vectorial Sturm-Liouville equations,” Inverse Problems, vol. 17, no. 5, pp. 1253–1294, 2001.View at: Publisher Site | Google Scholar | MathSciNet
N. P. Bondarenko, “Necessary and sufficient conditions for the solvability of the inverse problem for the matrix Sturm-Liouville operator,” Functional Analysis and its Applications, vol. 46, no. 1, pp. 53–57, 2012.View at: Publisher Site | Google Scholar | MathSciNet
V. Yurko, “Inverse problems for the matrix Sturm-Liouville equation on a finite interval,” Inverse Problems, vol. 22, no. 4, pp. 1139–1149, 2006.View at: Publisher Site | Google Scholar | MathSciNet
Y. V. Mykytyuk and N. S. Trush, “Inverse spectral problems for Sturm-Liouville operators with matrix-valued potentials,” Inverse Problems, vol. 26, no. 1, Article ID 015009, 2010.View at: Publisher Site | Google Scholar | MathSciNet
A. I. Zayed, Advances in Shannon's Sampling Theory, CRC Press, Boca Raton, Fla, USA, 1993.View at: MathSciNet
P. A. Binding and P. J. Browne, “Oscillation theory for indefinite Sturm-Liouville problems with eigenparameter-dependent boundary conditions,” Proceedings of the Royal Society of Edinburgh: Section A Mathematics, vol. 127, no. 6, pp. 1123–1136, 1997.View at: Publisher Site | Google Scholar | MathSciNet