Skip to main content

Sampling theories of boundary value problems with several internal points of discontinuity

Abstract

In this paper, our boundary value problem is a Dirac system with transmission conditions at several points of discontinuity. The main purpose of this paper is to derive the sampling theorems of this boundary value problem. To derive the sampling theorems including the construction of the Green’s matrix as well as the vector-valued eigenfunction expansion theorem, we briefly study the spectral analysis of the problem as in Levitan and Sargsjan (Introduction to Spectral Theory: Selfadjoint Ordinary Differential Operators, Translations of Mathematical Monographs, vol. 39, 1975; Sturm-Liouville and Dirac Operators, 1991) in a way similar to that of Fulton (Proc. R. Soc. Edinb., Sect. A 77:293-308, 1977). We derive sampling representations for transforms whose kernels are either solutions or the Green’s matrices of the problem. In the special case when our problem has one point of discontinuity, the obtained results coincide with the corresponding results in Tharwat et al. (Numer. Funct. Anal. Optim. 34:323-348, 2013).

1 Introduction

Let \(\mathfrak{H}(D)\) be a class of complex-valued functions defined on D, where D is a subset of \(\mathbb{C}\), which may coincide with \(\mathbb{C}\). We say that a sampling theorem holds for the class \(\mathfrak{H}(D)\) if there are two sequences \(\{\lambda_{k}\} _{k=1}^{\infty}\subset D\) and \(\{S_{k}(\lambda)\}_{k=1}^{\infty}\subset\mathfrak{H}(D)\), such that

$$ F(\lambda)=\sum_{k=1}^{\infty}F( \lambda_{k})S_{k}(\lambda),\quad \lambda\in D, F\in \mathfrak{H}(D), $$
(1.1)

where series (1.1) satisfies some sort of convergence. The points \(\{\lambda_{k}\}_{k=1}^{\infty}\), which are for practical reasons preferred to be real, will be called the sampling points and the functions \(\{S_{k}(\lambda)\}_{k=1}^{\infty}\) are called the sampling (reconstructing) functions. The name sampling comes about since the elements of \(\{S_{k}(\lambda)\} _{k=1}^{\infty}\), when a sampling theorem holds, can be recovered from their values at \(\{\lambda_{k}\}_{k=1}^{\infty}\). When a sampling theorem holds for a class \(\{S_{k}(\lambda)\}_{k=1}^{\infty}\), this may have many applications when the discrete (digital) state is preferable, like the role played by the Whittaker-Shannon-Kotel’nikov (WSK) sampling theorem [59] in signal processing.

The point now concerns the classes \(\mathfrak{H}(D)\), for which sampling theorems are guaranteed. In most of the known cases in sampling theory the classes of the sampled functions are integral transforms defined on a \(L^{2}\)-space. Now, in the following, we give an example.

Example 1.1

Let \(\mathfrak{H}(\mathbb{C})\) be the space of all bandlimited functions, i.e., functions of the form of finite Fourier transforms,

$$ F(\lambda)= \frac{1}{\sqrt{2\pi}} \int_{-\pi}^{\pi} f(x)e^{ix\lambda} \,dx,\quad f(\cdot) \in L^{2}(-\pi,\pi), \lambda\in\mathbb{C}. $$
(1.2)

Thus, \(\mathfrak{H}(\mathbb{C})\) is the Paley-Wiener space \(\mathit{PW}_{\pi }^{2}\) of entire functions of exponential type at most π, which belong to \(L^{2}(\mathbb{R})\) when restricted to \(\mathbb{R}\), cf. [5, 10, 11]. Elements of \(\mathfrak{H}(\mathbb{C})\), in this case, have sampling expansions of the form (1.1) with

$$ \lambda_{k}=k; \qquad S_{k}(\lambda)= \frac{\sin\pi(\lambda -k)}{\pi(\lambda-k)}, \quad k\in\mathbb{Z}, $$
(1.3)

where the sampling expansions converge absolutely on \(\mathbb{C}\) and uniformly on both \(\mathbb{R}\) and compact subsets of \(\mathbb{C}\), [12], p.262. This result is the well-known WSK sampling theorem [59] which has many applications in signal processing.

The above example is a special case of a generalized sampling theorem for generalized integral transforms. It is known as Kramer’s sampling theorem [5, 1316] and reads as follows.

Theorem 1.1

Let \(\{\lambda_{k}\}_{k=-\infty}^{\infty}\) be a sequence of real numbers. Let \(I\subseteq\mathbb{R}\) and \(\mathcal{K}(\cdot,\lambda ):I\times\mathbb{C}\to\mathbb{C}\) be a function such that

  1. 1.

    \(\mathcal{K}(\cdot,\lambda)\in L^{2}(I)\) for all \(\lambda\in \mathbb{C}\),

  2. 2.

    the sequence \(\{\mathcal{K}(\cdot,\lambda_{k})\}_{k=-\infty }^{\infty}\) forms a complete orthogonal set in \(L^{2}(I)\).

Let \(\mathfrak{H}(\mathbb{C})\) be the set of integral transforms

$$ F(\lambda)= \int_{I} f(x)\mathcal{K}(x,\lambda) \,dx, \quad f(\cdot)\in L^{2}(I). $$
(1.4)

Then elements of \(\mathfrak{H}(\mathbb{C})\) admit the sampling expansions

$$ F(\lambda)=\sum_{k=-\infty}^{\infty} F( \lambda_{k})\frac{\int_{I} \mathcal{K}(x,\lambda)\overline{\mathcal {K}(x,\lambda_{k})} \,dx}{\int_{I}|\mathcal{K}(x,\lambda_{k})|^{2}\,dx} ,\quad F\in\mathfrak{H}(\mathbb{C}). $$
(1.5)

Series (1.4) converges uniformly wherever \(\Vert \mathcal{K}(\cdot, \lambda)\Vert_{L^{2}(I)}\) is bounded.

We can see that Example 1.1 above is a special case of Theorem 1.1. There are many problems motivated by Kramer’s theorem. The first concerns the situation for which we can obtain the kernel and the sampling points such that the theorem holds. The second concerns the relationship between Kramer’s sampling expansions and Lagrange-type interpolation ones. Notice that the sampling series in Example 1.1 is in the form of a Lagrange-type interpolation expansion. A partial answer to the first question was given by Weiss [17], before the existence of Kramer’s theorem. In this short note [17], the kernel is extracted from the second order boundary value problem

$$\begin{aligned}& -y''(x)+q(x)y(x)=\lambda y(x),\quad - \infty< a\le x\le b< \infty, \end{aligned}$$
(1.6)
$$\begin{aligned}& \cos\alpha y(a)+\sin\alpha y'(a)=0, \end{aligned}$$
(1.7)
$$\begin{aligned}& \cos\beta y(b)+\sin\beta y'(b)=0, \end{aligned}$$
(1.8)

where \(q(\cdot)\) is continuous and real valued on \([a, b]\) and \(\alpha ,\beta\in[0,\pi)\). If we take \(u(\cdot,\lambda)\) to be the solution of (1.6) determined by the initial conditions

$$ u(a,\lambda)=\sin\alpha, \qquad u'(a,\lambda)=-\cos \alpha,\quad \lambda\in\mathbb{\mathbb{C}}, $$
(1.9)

then \(u(\cdot,\lambda)\) satisfies Kramer’s theorem with respect to the sequence of the eigenvalues, i.e., the sequence of the eigenvalues, \(\{\lambda_{k}\} _{k=0}^{\infty}\), will be taken as the sampling points. For the second question, it is shown in many articles that Kramer’s expansions are nothing more than Lagrange-type interpolation ones when the kernels of the sampled integral transforms are solutions of certain self adjoint eigenvalue problems or are to be expressed in terms of the Green’s functions of these problems; see [1820]. Using the connection between Kramer’s theorem and the eigenvalue problems most of the known sampling expansions can easily be obtained.

In [21, 22], the authors discussed sampling theorems of regular one-dimensional continuous Dirac systems; see also [23]. In the case of boundary value problems with one point of discontinuity, Tharwat et al. [4], discussed sampling theorems of Dirac systems; see also [24, 25]. In the present paper we derive sampling theorems associated with Dirac systems with transmission conditions at several points of discontinuity. That is, we define classes of transforms associated with these problems and derive sampling expansions for these classes. We use both solutions and the Green’s matrices to define the sampled transforms. In all cases the obtained sampling expansions are written in the form of a Lagrange interpolation series. The uniform convergence of the obtained sampling expansion as well as the analytic and growth properties of the sampled transforms are given. However, sampling theories associated with Dirac systems, which have several points of discontinuity, do not exist as far as we know. An illustrative example is given in detail in the last section.

2 Preliminaries

In this section we introduce some of the notations and relations that will be used in the sequel; then we prove some useful lemmas and theorems. We consider the Dirac system

$$\begin{aligned}& \left . \textstyle\begin{array}{l} r(x) u_{2}^{\prime}(x)-q_{1}(x)u_{1}(x)=\lambda u_{1}(x), \\ r(x)u_{1}^{\prime}(x)+q_{2}(x)u_{2}(x)=-\lambda u_{2}(x), \end{array}\displaystyle \right \} \quad x\in[a,c_{1})\cup(c_{1},c_{2}) \cup(c_{2},c_{3})\cup\cdots\cup(c_{m},b], \end{aligned}$$
(2.1)
$$\begin{aligned}& \mathcal{U}_{1}(u):=\sin\alpha u_{1}(a)+ \cos\alpha u_{2}(a)=0, \end{aligned}$$
(2.2)
$$\begin{aligned}& \mathcal{U}_{2}(u):=\sin\beta u_{1}(b)+\cos \beta u_{2}(b)=0, \end{aligned}$$
(2.3)

and the transmission conditions

$$ \mathcal{L}_{k}(u):= \left[ \textstyle\begin{array}{@{}c@{}} \gamma_{k} u_{1}(c_{k}^{-}) \\ \gamma_{k}^{\prime}u_{2}(c_{k}^{-}) \end{array}\displaystyle \right]- \left[ \textstyle\begin{array}{@{}c@{}} \delta_{k} u_{1}(c_{k}^{+}) \\ \delta_{k}^{\prime} u_{2}(c_{k}^{+}) \end{array}\displaystyle \right]=0, \quad k=1,2,\ldots,m, $$
(2.4)

where \(\lambda\in\mathbb{C}\);

$$r(x):=\left \{ \textstyle\begin{array}{l@{\quad}l} r_{1}, &x\in[a, c_{1}), \\ r_{k}, & x\in(c_{k-1}, c_{k})\ (k=2,3,\ldots,m), \\ r_{m+1}, & x\in(c_{m}, b], \end{array}\displaystyle \right . $$

\(r_{k}>0\) (\(k=1,2,\ldots,m+1\)) are real numbers; the real-valued functions \(q_{1}(\cdot)\) and \(q_{2}(\cdot)\) are continuous in \([a,c_{1})\), \((c_{k-1},c_{k})\), and \((c_{m},b]\) (\(k=2,3,\ldots,m\)), and have finite limits \(q_{1}(c_{k}^{\pm}):=\lim _{x\rightarrow c_{k}^{\pm}}q_{1}(x)\), \(q_{2} (c_{k}^{\pm}):=\lim_{x\rightarrow c_{k}^{\pm}}q_{2}(x)\) (\(k=1,2,\ldots,m\)); \(\gamma_{k}, \delta_{k}, \gamma_{k}^{\prime}, \delta _{k}^{\prime}\in \mathbb{R}\), \(\gamma_{k}, \delta_{k}, \gamma _{k}^{\prime}, \delta_{k}^{\prime}\neq0\) and \(\alpha,\beta\in [0,\pi)\); see [26, 27].

Let \(\mathcal{H}\) be the Hilbert space

$$ \mathcal{H}:= \left \{u(x)=\left ( \textstyle\begin{array}{@{}c@{}} u_{1}(x) \\ u_{2}(x) \end{array}\displaystyle \right ), u_{1}(x), u_{2}(x)\in\bigoplus _{k=1}^{m+1} L^{2}(c_{k-1},c_{k}), a=c_{0}, b=c_{m+1} \right \}. $$
(2.5)

The inner product of \(\mathcal{H}\) is defined by

$$\begin{aligned} \bigl\langle u(\cdot),v(\cdot)\bigr\rangle _{\mathcal{H}} :=& \frac{1}{r_{1}} \int^{c_{1}}_{a}u^{\top}(x)\overline{v}(x) \,dx+\sum_{k=2}^{m}\frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}u^{\top}(x)\overline{v}(x) \,dx \\ &{}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}u^{\top }(x)\overline{v}(x) \,dx, \end{aligned}$$
(2.6)

where \(D_{i}=\frac{\delta_{i}\delta_{i}^{\prime}}{\gamma _{i}\gamma_{i}^{\prime}}\), \(D_{i}>0\), \(i=1,2,\ldots,m\), and denotes the matrix transpose,

$$u(x)=\left ( \textstyle\begin{array}{@{}c@{}} u_{1}(x) \\ u_{2}(x) \end{array}\displaystyle \right ) , v(x)= \left ( \textstyle\begin{array}{@{}c@{}} v_{1}(x) \\ v_{2}(x) \end{array}\displaystyle \right ) \in\mathcal{H}, \qquad u_{i}(\cdot), v_{i}(\cdot)\in\bigoplus _{k=1}^{m+1} L^{2}(c_{k-1},c_{k}), \quad i=1,2. $$

For vector-valued functions \(u(x)\), which are defined on \([a,c_{1})\cup (c_{1},c_{2})\cup(c_{2},c_{3})\cup\cdots\cup(c_{m},b]\) and have a finite limit \(u(c_{k}^{\pm}):=\lim_{x\rightarrow c_{k}^{\pm}}u(x)\) (\(k=1,2,\ldots,m\)), by \(u_{(k)}(x)\) (\(k=1,2,\ldots,m+1\)) we denote the functions

$$ \begin{aligned} &u_{(1)}(x)=\left \{ \textstyle\begin{array}{l@{\quad}l} u(x), & x\in[a,c_{1}) , \\ u(c_{1}^{-}), & x=c_{1} , \end{array}\displaystyle \right .\qquad u_{(i)}(x)=\left \{ \textstyle\begin{array}{l@{\quad}l} u(c_{i-1}^{+}), & x=c_{i-1}, \\ u(x), & x\in(c_{i-1},c_{i}), \\ u(c_{i}^{-}), &x=c_{i}, \end{array}\displaystyle \right .\quad i=2,3,\ldots,m, \\ &u_{(m+1)}(x)=\left \{ \textstyle\begin{array}{l@{\quad}l} u(x), & x\in(c_{m},b], \\ u(c_{m}^{+}), & x=c_{m}, \end{array}\displaystyle \right . \end{aligned} $$
(2.7)

which are defined on \(\mathcal{I}_{1}:=[a,c_{1}]\), \(\mathcal {I}_{k}:=[c_{k-1},c_{k}]\) (\(k=2,3,\ldots,m\)) and \(\mathcal {I}_{m+1}:=[c_{m},b]\), respectively.

In the following lemma, we will prove that the eigenvalues of the problem (2.1)-(2.4) are real.

Lemma 2.1

The eigenvalues of the problem (2.1)-(2.4) are real.

Proof

Assume to the contrary that \(\lambda_{0}\) is a non-real eigenvalue of problem (2.1)-(2.4). Let \(\bigl ({\scriptsize\begin{matrix}{} u_{1}(x) \cr u_{2}(x)\end{matrix}} \bigr )\) be a corresponding (non-trivial) vector-valued eigenfunction. By (2.1), we have

$$\begin{aligned}& r(x)\frac{d}{dx} \bigl\{ u_{1}(x)\overline{u}_{2}(x)- \overline{u}_{1}(x)u_{2}(x) \bigr\} \\& \quad =(\overline{ \lambda}_{0}-\lambda_{0}) \bigl\{ \bigl\vert u_{1}(x)\bigr\vert ^{2}+\bigl\vert u_{2}(x) \bigr\vert ^{2} \bigr\} ,\quad x\in[a,c_{1})\cup (c_{1},c_{2})\cup(c_{2},c_{3})\cup \cdots\cup(c_{m},b]. \end{aligned}$$

Integrating the above equation through \([a, c_{1})\), \((c_{k-1}, c_{k})\) (\(k=2,3,\ldots,m\)), and \((c_{m}, b]\), we obtain

$$\begin{aligned}& \frac{\overline{\lambda}_{0}-\lambda_{0}}{r_{1}} \biggl[ \int_{a}^{c_{1}} \bigl(\bigl\vert u_{1}(x) \bigr\vert ^{2}+\bigl\vert u_{2}(x)\bigr\vert ^{2} \bigr) \,dx \biggr] \\& \quad =u_{1}\bigl(c_{1}^{-} \bigr)\overline{u}_{2}\bigl(c_{1}^{-}\bigr)- \overline{u}_{1}\bigl(c_{1}^{-} \bigr)u_{2}\bigl(c_{1}^{-}\bigr)- \bigl[u_{1}(a)\overline{u}_{2}(a)- \overline{u}_{1}(a)u_{2}(a) \bigr], \end{aligned}$$
(2.8)
$$\begin{aligned}& \frac{\overline{\lambda}_{0}-\lambda_{0}}{r_{k}} \biggl[ \int_{c_{k-1}}^{c_{k}} \bigl(\bigl\vert u_{1}(x) \bigr\vert ^{2}+\bigl\vert u_{2}(x)\bigr\vert ^{2} \bigr) \,dx \biggr] \\& \quad =u_{1}\bigl(c_{k}^{-} \bigr)\overline{u}_{2}\bigl(c_{k}^{-}\bigr)- \overline{u}_{1}\bigl(c_{k}^{-} \bigr)u_{2}\bigl(c_{k}^{-}\bigr) - \bigl[u_{1}\bigl(c_{k-1}^{+}\bigr) \overline{u}_{2}\bigl(c_{k-1}^{+}\bigr)- \overline{u}_{1}\bigl(c_{k-1}^{+} \bigr)u_{2}\bigl(c_{k-1}^{+}\bigr) \bigr], \end{aligned}$$
(2.9)

for \(k=2,3,\ldots,m\), and

$$\begin{aligned}& \frac{\overline{\lambda}_{0}-\lambda_{0}}{r_{m+1}} \biggl[ \int_{c_{m}}^{b} \bigl(\bigl\vert u_{1}(x) \bigr\vert ^{2}+\bigl\vert u_{2}(x)\bigr\vert ^{2} \bigr) \,dx \biggr] \\& \quad =u_{1}(b)\overline{u}_{2}(b)- \overline{u}_{1}(b)u_{2}(b)- \bigl[u_{1} \bigl(c_{m}^{+}\bigr)\overline{u}_{2} \bigl(c_{m}^{+}\bigr)- \overline{u}_{1} \bigl(c_{m}^{+}\bigr)u_{2}\bigl(c_{m}^{+} \bigr) \bigr]. \end{aligned}$$
(2.10)

Then from (2.2), (2.3), and the transmission conditions, we have, respectively,

$$\begin{aligned}& u_{1}(a)\overline{u}_{2}(a)- \overline{u}_{1}(a)u_{2}(a)=0, \\& u_{1}(b)\overline{u}_{2}(b)- \overline{u}_{1}(b)u_{2}(b)=0, \end{aligned}$$

and

$$\begin{aligned}& u_{1}\bigl(c_{k}^{-}\bigr)\overline{u}_{2} \bigl(c_{k}^{-}\bigr)- \overline{u}_{1} \bigl(c_{k}^{-}\bigr)u_{2}\bigl(c_{k}^{-} \bigr) \\& \quad =\frac{\delta_{k}\delta _{k}^{\prime}}{\gamma_{k}\gamma_{k}^{\prime }}\bigl[u_{1}\bigl(c_{k}^{+} \bigr)\overline{u}_{2}\bigl(c_{k}^{+}\bigr)- \overline{u}_{1}\bigl(c_{k}^{+} \bigr)u_{2}\bigl(c_{k}^{+}\bigr)\bigr],\quad k=1,2, \ldots,m. \end{aligned}$$

Since \(\lambda_{0}\neq\overline{\lambda}_{0}\), it follows from the last three equations and (2.8), (2.9), (2.10) that

$$\begin{aligned}& \frac{1}{r_{1}} \int_{a}^{c_{1}} \bigl(\bigl\vert u_{1}(x) \bigr\vert ^{2}+\bigl\vert u_{2}(x)\bigr\vert ^{2} \bigr) \,dx +\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int _{c_{k-1}}^{c_{k}} \bigl(\bigl\vert u_{1}(x)\bigr\vert ^{2}+\bigl\vert u_{2}(x) \bigr\vert ^{2} \bigr) \,dx \\& \quad {} +\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int_{c_{m}}^{b} \bigl(\bigl\vert u_{1}(x) \bigr\vert ^{2}+\bigl\vert u_{2}(x)\bigr\vert ^{2} \bigr) \,dx=0. \end{aligned}$$
(2.11)

Then \(u_{i}(x)=0\), \(i=1,2\), and this is a contradiction. Consequently, \(\lambda_{0}\) must be real. □

Lemma 2.2

Let \(\lambda_{1}\) and \(\lambda_{2}\) be two different eigenvalues of the problem (2.1)-(2.4). Then the corresponding vector-valued eigenfunctions \(u(x,\lambda_{1})\) and \(v(x,\lambda_{2})\) of this problem satisfy the following equality:

$$\begin{aligned} \begin{aligned}[b] &\frac{1}{r_{1}} \int_{a}^{c_{1}}u^{\top}(x, \lambda_{1})v(x,\lambda _{2}) \,dx+\sum _{k=2}^{m}\frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int _{c_{k-1}}^{c_{k}}u^{\top}(x, \lambda_{1})v(x,\lambda_{2}) \,dx \\ &\quad {}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int_{c_{m}}^{b}u^{\top }(x, \lambda_{1})v(x,\lambda_{2}) \,dx=0. \end{aligned} \end{aligned}$$
(2.12)

Proof

By (2.1) we obtain

$$\begin{aligned}& r(x)\frac{d}{dx} \bigl\{ u_{1}(x,\lambda_{1})v_{2}(x, \lambda_{2})- u_{2}(x,\lambda_{2})v_{1}(x, \lambda_{1}) \bigr\} \\& \quad =(\lambda _{2}-\lambda_{1}) \bigl\{ u_{1}(x,\lambda_{1})v_{1}(x, \lambda_{2}) +u_{2}(x,\lambda_{1})v_{2}(x, \lambda_{2}) \bigr\} . \end{aligned}$$

Integrating the above equation through \([a, c_{1})\), \((c_{k-1}, c_{k})\) (\(k=2,3,\ldots,m\)), and \((c_{m}, b]\), and taking into account that \(u(x,\lambda_{1})\) and \(v(x,\lambda_{2})\) satisfy (2.2)-(2.4), we obtain (2.12), where \(\lambda_{1}\neq\lambda_{2}\). □

Now we will construct a special fundamental system of solutions of (2.1) for λ not being an eigenvalue. By virtue of Theorem 1.1 in [1], we will define the two solutions of (2.1)

$$\varphi(\cdot,\lambda)=\left ( \textstyle\begin{array}{@{}c@{}} \varphi_{1}(\cdot,\lambda) \\ \varphi_{2}(\cdot,\lambda) \end{array}\displaystyle \right ), \qquad \chi(\cdot,\lambda)=\left ( \textstyle\begin{array}{@{}c@{}} \chi_{1}(\cdot,\lambda) \\ \chi_{2}(\cdot,\lambda) \end{array}\displaystyle \right ), $$

where

$$\begin{aligned}& \begin{aligned} &\varphi_{1}(x,\lambda)=\left \{ \textstyle\begin{array}{l@{\quad}l} \varphi_{1,1}(x,\lambda), & x\in[c_{0},c_{1}) , \\ \varphi_{1,2}(x,\lambda), & x\in(c_{1},c_{2}) , \\ \ldots, \\ \varphi_{1,m+1}(x,\lambda), & x\in(c_{m},c_{m+1}] , \end{array}\displaystyle \right . \\ &\varphi_{2}(x,\lambda)=\left \{ \textstyle\begin{array}{l@{\quad}l} \varphi_{2,1}(x,\lambda), & x\in[c_{0},c_{1}) , \\ \varphi_{2,2}(x,\lambda), & x\in(c_{1},c_{2}) , \\ \ldots, \\ \varphi_{2,m+1}(x,\lambda), & x\in(c_{m},c_{m+1}] , \end{array}\displaystyle \right . \end{aligned} \end{aligned}$$
(2.13)
$$\begin{aligned}& \begin{aligned} &\chi_{1}(x,\lambda)=\left \{ \textstyle\begin{array}{l@{\quad}l} \chi_{1,1}(x,\lambda), & x\in[c_{0},c_{1}) , \\ \chi_{1,2}(x,\lambda), & x\in(c_{1},c_{2}) , \\ \ldots, \\ \chi_{1,m+1}(x,\lambda), & x\in(c_{m},c_{m+1}] , \end{array}\displaystyle \right . \\ &\chi_{2}(x,\lambda)=\left \{ \textstyle\begin{array}{l@{\quad}l} \chi_{2,1}(x,\lambda), & x\in[c_{0},c_{1}) , \\ \chi_{2,2}(x,\lambda), & x\in(c_{1},c_{2}) , \\ \ldots, \\ \chi_{2,m+1}(x,\lambda), & x\in(c_{m},c_{m+1}] , \end{array}\displaystyle \right . \end{aligned} \end{aligned}$$
(2.14)

as follows, where \(c_{0}=a\) and \(c_{m+1}=b\). By virtue of Theorem 1.1 in [1] the initial-value problem

$$\begin{aligned}& \begin{aligned} &r(x)u_{2}^{\prime}(x)-q_{1}(x)u_{1}(x)= \lambda u_{1}(x), \\ &r(x) u_{1}^{\prime}(x)+q_{2}(x)u_{2}(x)=- \lambda u_{2}(x), \quad x\in (a,c_{1}), \end{aligned} \end{aligned}$$
(2.15)
$$\begin{aligned}& u_{1}(a)=\cos\alpha, \qquad u_{2}(a)=-\sin \alpha, \end{aligned}$$
(2.16)

has a unique solution \(u=\bigl ({\scriptsize\begin{matrix}{} \varphi_{1,1}(x,\lambda)\cr \varphi_{2,1}(x,\lambda)\end{matrix}} \bigr )\), which is an entire function of \(\lambda\in\mathbb{C}\) for each fixed \(x\in[a,c_{1}]\). Similarly, employing the same method as in proof of Theorem 1.1 in [1], we see that the problem

$$\begin{aligned}& \begin{aligned} &r(x)u_{2}^{\prime}(x)-q_{1}(x)u_{1}(x)= \lambda u_{1}(x), \\ &r(x)u_{1}^{\prime}(x)+q_{2}(x)u_{2}(x)=- \lambda u_{2}(x), \quad x\in(c_{m},b), \end{aligned} \end{aligned}$$
(2.17)
$$\begin{aligned}& u_{1}(b)=\cos\beta, \qquad u_{2}(b)=-\sin \beta, \end{aligned}$$
(2.18)

has a unique solution \(u=\bigl ( {\scriptsize\begin{matrix}{} \chi_{1,m+1}(x,\lambda)\cr \chi_{2,m+1}(x,\lambda) \end{matrix}} \bigr )\), which is an entire function of the parameter λ for each fixed \(x\in[c_{m},b]\).

Now the functions \(\varphi_{i,k+1}(x,\lambda)\) and \(\chi _{i,k}(x,\lambda)\) are defined in terms of \(\varphi_{i,k}(x,\lambda)\) and \(\chi_{i,k+1}(x,\lambda)\), \(i=1,2\), \(k=1,2,\ldots,m\), respectively, as follows. The initial-value problem

$$\begin{aligned}& \begin{aligned} &r(x)u_{2}^{\prime}(x)-q_{1}(x)u_{1}(x)= \lambda u_{1}(x), \\ &r(x)u_{1}^{\prime}(x)+q_{2}(x)u_{2}(x)=- \lambda u_{2}(x), \quad x\in (c_{k},c_{k+1}), \end{aligned} \end{aligned}$$
(2.19)
$$\begin{aligned}& u_{1}(c_{k})=\frac{\gamma_{k}}{\delta_{k}}\varphi _{1,k}\bigl(c_{k}^{-},\lambda\bigr),\qquad u_{2}(c_{k})=\frac{\gamma_{k}^{\prime}}{\delta_{k}^{\prime }}\varphi_{2,k} \bigl(c_{k}^{-},\lambda\bigr), \quad k=1,2,\ldots,m, \end{aligned}$$
(2.20)

has a unique solution \(u=\bigl ( {\scriptsize\begin{matrix}{} \varphi_{1,k+1}(x,\lambda)\cr \varphi_{2,k+1}(x,\lambda) \end{matrix}} \bigr )\) for each \(\lambda\in \mathbb{C}\).

Similarly, the following problem also has a unique solution \(u=\bigl ( {\scriptsize\begin{matrix}{} \chi_{1,k}(x,\lambda)\cr \chi_{2,k}(x,\lambda)\end{matrix}} \bigr )\):

$$\begin{aligned}& \begin{aligned} &r(x)u_{2}^{\prime}(x)-q_{1}(x)u_{1}(x)= \lambda u_{1}(x), \\ &r(x)u_{1}^{\prime}(x)+q_{2}(x)u_{2}(x)=- \lambda u_{2}(x), \quad x\in (c_{k-1},c_{k}), \end{aligned} \end{aligned}$$
(2.21)
$$\begin{aligned}& u_{1}(c_{k})=\frac{\delta_{k}}{\gamma_{k}}\chi _{1,k+1}\bigl(c_{k}^{+},\lambda\bigr), \qquad u_{2}(c_{k})=\frac{\delta _{k}^{\prime}}{\gamma_{k}^{\prime}} \chi_{2,k+1} \bigl(c_{k}^{+},\lambda\bigr),\quad k=1,2,\ldots,m. \end{aligned}$$
(2.22)

By virtue of equations (2.20) and (2.22) these solutions satisfy both transmission conditions (2.4). These functions are entire in λ for all \(x\in[a,c_{1})\cup (c_{1},c_{2})\cup(c_{2},c_{3})\cup\cdots\cup(c_{m},b]\).

Let \(W(\varphi,\chi)({\cdot,\lambda})\) denote the Wronskian of \(\varphi(\cdot,\lambda)\) and \(\chi(\cdot,\lambda)\) defined in [2], p.194, i.e.,

$$W(\varphi,\chi) (\cdot,\lambda):=\left \vert \textstyle\begin{array}{@{}c@{\quad}c@{}} \varphi_{1}(\cdot,\lambda) & \varphi_{2}(\cdot,\lambda) \\ \chi_{1}(\cdot,\lambda) & \chi_{2}(\cdot,\lambda) \end{array}\displaystyle \right \vert . $$

It is obvious that the Wronskians

$$\begin{aligned} \omega_{k}(\lambda) :=&W(\varphi,\chi) (x,\lambda) \\ =& \varphi _{1,k}(x,\lambda) \chi_{2,k}(x,\lambda)- \varphi_{2,k}(x,\lambda) \chi_{1,k}(x,\lambda), \quad x\in \mathcal{I}_{k}, k=1,2,\ldots,m+1, \end{aligned}$$
(2.23)

are independent of \(x\in\mathcal{I}_{k}\) and are entire functions. Taking into account (2.20) and (2.22), a short calculation gives

$$\omega_{1}(\lambda)=D_{1}\omega_{2}(\lambda)= \cdots=\prod_{i=1}^{m}D_{i} \omega_{m+1}(\lambda), $$

for each \(\lambda\in\mathbb{C}\).

Corollary 2.3

The zeros of the functions \(\omega_{k}(\lambda)\) (\(k=1,2,\ldots,m+1\)) coincide.

Then we may introduce into the consideration the characteristic function \(\omega(\lambda)\) as

$$ \omega(\lambda):=\omega_{1}(\lambda)=D_{1} \omega_{2}(\lambda )=\cdots=\prod_{i=1}^{m}D_{i} \omega_{m+1}(\lambda). $$
(2.24)

Lemma 2.4

All eigenvalues of problem (2.1)-(2.4) are just zeros of the function \(\omega(\lambda)\).

Proof

Since the functions \(\varphi_{1}(x,\lambda)\) and \(\varphi _{2}(x,\lambda)\) satisfy the boundary condition (2.2) and the transmission conditions (2.4), to find the eigenvalues of the (2.1)-(2.4) we have to insert the functions \(\varphi _{1}(x,\lambda)\) and \(\varphi_{2}(x,\lambda)\) in the boundary condition (2.3) and find the roots of this equation. □

In the following lemma, we show that all eigenvalues of the problem (2.1)-(2.4) are simple.

Lemma 2.5

The eigenvalues of the boundary value problem (2.1)-(2.4) form an at most countable set without finite limit points. All eigenvalues of the boundary value problem (2.1)-(2.4) (of \(\omega(\lambda)\)) are simple.

Proof

The eigenvalues are the zeros of the entire function occurring on the left-hand side, see equation (2.24), in

$$\prod_{i=1}^{m}D_{i}\bigl[\sin \beta\varphi_{1,m+1}(b,\lambda)+ \cos\beta\varphi_{2,m+1}(b,\lambda) \bigr]=0. $$

We have shown (see Lemma 2.1) that this function does not vanish for non-real λ. In particular, it does not vanish identically. Therefore, its zeros form an at most countable set without finite limit points.

By (2.1) we obtain for \(\lambda,\mu\in\mathbb{C}\), \(\lambda\neq\mu\),

$$r(x)\frac{d}{dx} \bigl\{ \varphi_{1}(x,\lambda) \varphi_{2}(x,\mu)- \varphi_{1}(x,\mu)\varphi_{2}(x, \lambda) \bigr\} =(\mu-\lambda) \bigl\{ \varphi_{1}(x,\lambda) \varphi_{1}(x,\mu) +\varphi_{2}(x,\lambda) \varphi_{2}(x,\mu) \bigr\} . $$

Integrating the above equation through \([a, c_{1})\), \((c_{k-1}, c_{k})\) (\(k=2,3,\ldots,m\)) and \((c_{m}, b]\), we obtain

$$\begin{aligned}& \frac{\mu-\lambda}{r_{1}} \biggl[ \int_{a}^{c_{1}} \bigl(\varphi_{1,1}(x, \lambda)\varphi_{1,1}(x,\mu) +\varphi_{2,1}(x,\lambda) \varphi_{2,1}(x,\mu) \bigr) \,dx \biggr] \\& \quad = \varphi_{1,1}\bigl(c_{1}^{-},\lambda\bigr) \varphi_{2,1}\bigl(c_{1}^{-},\mu\bigr)- \varphi_{1,1}\bigl(c_{1}^{-},\lambda\bigr) \varphi_{2,1}\bigl(c_{1}^{-},\mu\bigr) \\& \qquad {}- \bigl(\varphi_{1,1}(a,\lambda)\varphi_{2,1}(a,\mu)- \varphi_{1,1}(a,\lambda)\varphi_{2,1}(a,\mu) \bigr), \end{aligned}$$
(2.25)
$$\begin{aligned}& \frac{\mu-\lambda}{r_{k}} \biggl[ \int_{c_{k-1}}^{c_{k}} \bigl(\varphi_{1,k}(x, \lambda)\varphi _{1,k}(x,\mu) +\varphi_{2,k}(x,\lambda) \varphi_{2,k}(x,\mu) \bigr) \,dx \biggr] \\& \quad =\varphi_{1,k}\bigl(c_{k}^{-},\lambda\bigr) \varphi _{2,k}\bigl(c_{k}^{-},\mu\bigr)- \varphi_{1,k}\bigl(c_{k}^{-},\lambda\bigr) \varphi_{2,k}\bigl(c_{k}^{-},\mu\bigr) \\& \qquad {}- \bigl[\varphi_{1,k}\bigl(c_{k-1}^{+}, \lambda\bigr)\varphi _{2,k}\bigl(c_{k-1}^{+},\mu \bigr)- \varphi_{1,k}\bigl(c_{k-1}^{+},\lambda\bigr) \varphi_{2,k}\bigl(c_{k-1}^{+},\mu \bigr) \bigr], \end{aligned}$$
(2.26)

for \(k=2,3,\ldots,m\), and

$$\begin{aligned}& \frac{\mu-\lambda}{r_{m+1}} \biggl[ \int_{c_{m}}^{b} \bigl(\varphi_{1,m+1}(x, \lambda)\varphi _{1,m+1}(x,\mu) +\varphi_{2,m+1}(x,\lambda) \varphi_{2,m+1}(x,\mu) \bigr) \,dx \biggr] \\& \quad =\varphi_{1,m+1}(b,\lambda)\varphi_{2,m+1}(b,\mu) - \varphi_{1,m+1}(b,\lambda)\varphi_{2,m+1}(b,\mu) \\& \qquad {}- \bigl[\varphi _{1,m+1}\bigl(c_{m}^{+}, \lambda\bigr)\varphi_{2,m+1}\bigl(c_{m}^{+},\mu\bigr)- \varphi_{1,m+1}\bigl(c_{m}^{+},\lambda\bigr) \varphi_{2,m+1}\bigl(c_{m}^{+},\mu \bigr) \bigr]. \end{aligned}$$
(2.27)

Then from (2.16) and the transmission conditions (2.20), we have

$$\begin{aligned}& \frac{\mu-\lambda}{r_{1}} \biggl[ \int_{a}^{c_{1}} \bigl(\varphi_{1,1}(x, \lambda)\varphi_{1,1}(x,\mu) +\varphi_{2,1}(x,\lambda) \varphi_{2,1}(x,\mu) \bigr) \,dx \biggr] \\& \qquad {} +(\mu-\lambda)\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \biggl[ \int_{c_{k-1}}^{c_{k}} \bigl(\varphi_{1,k}(x, \lambda)\varphi _{1,k}(x,\mu) +\varphi_{2,k}(x,\lambda) \varphi_{2,k}(x,\mu) \bigr) \,dx \biggr] \\& \qquad {} +(\mu-\lambda) \biggl[\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int_{c_{m}}^{b} \bigl(\varphi_{1,m+1}(x, \lambda)\varphi _{1,m+1}(x,\mu) +\varphi_{2,m+1}(x,\lambda) \varphi_{2,m+1}(x,\mu) \bigr) \,dx \biggr] \\& \quad =\prod_{i=1}^{m}D_{i} \bigl[\varphi_{1,m+1}(b,\lambda)\varphi _{2,m+1}(b,\mu)- \varphi_{1,m+1}(b,\lambda)\varphi_{2,m+1}(b,\mu)\bigr]. \end{aligned}$$
(2.28)

Dividing both sides of (2.28) by \((\lambda-\mu)\) and by letting \(\mu \rightarrow \lambda\), we arrive at the relation

$$\begin{aligned}& -\frac{1}{r_{1}} \biggl[ \int_{a}^{c_{1}} \bigl(\bigl\vert \varphi_{1,1}(x,\lambda)\bigr\vert ^{2}+\bigl\vert \varphi _{2,1}(x,\lambda)\bigr\vert ^{2} \bigr) \,dx \biggr] \\& \qquad {}-\sum _{k=2}^{m}\frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \biggl[ \int_{c_{k-1}}^{c_{k}} \bigl(\bigl\vert \varphi_{1,k}(x,\lambda)\bigr\vert ^{2} +\bigl\vert \varphi_{2,k}(x,\lambda)\bigr\vert ^{2} \bigr) \,dx \biggr] \\& \qquad {} -\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \biggl[ \int_{c_{m}}^{b} \bigl(\bigl\vert \varphi_{1,m+1}(x,\lambda)\bigr\vert ^{2} +\bigl\vert \varphi_{2,m+1}(x,\lambda)\bigr\vert ^{2} \bigr) \,dx \biggr] \\& \quad =\prod_{i=1}^{m}D_{i} \biggl(\varphi_{2,m+1}(b,\lambda)\frac {\partial \varphi_{1,m+1}(b,\lambda)}{\partial\lambda}- \varphi_{1,m+1}(b, \lambda)\frac{\partial \varphi_{2,m+1}(b,\lambda)}{\partial\lambda} \biggr). \end{aligned}$$
(2.29)

We show that the equation

$$ \omega(\lambda)=-\prod_{i=1}^{m}D_{i} \bigl(\sin\beta\varphi _{1,m+1}(b,\lambda)+ \cos\beta\varphi_{2,m+1}(b, \lambda) \bigr)=0 $$
(2.30)

has only simple roots. Assume the converse, i.e., equation (2.30) has a double root \(\lambda^{\ast}\), say. Then the following two equations hold:

$$\begin{aligned}& \sin\beta\varphi_{1,m+1}\bigl(b,\lambda^{\ast} \bigr)+ \cos\beta\varphi_{2,m+1}\bigl(b,\lambda^{\ast}\bigr)=0, \end{aligned}$$
(2.31)
$$\begin{aligned}& \sin\beta\frac{\partial\varphi_{1,m+1}(b,\lambda^{\ast })}{\partial\lambda} +\cos\beta\frac{\partial\varphi_{2,m+1}(b,\lambda^{\ast })}{\partial\lambda}=0. \end{aligned}$$
(2.32)

Equations (2.31) and (2.32) imply that

$$ \varphi_{2,m+1}\bigl(b,\lambda^{\ast}\bigr) \frac{\partial \varphi_{1,m+1}(b,\lambda^{\ast})}{\partial\lambda}- \varphi_{1,m+1}\bigl(b,\lambda^{\ast}\bigr) \frac{\partial \varphi_{2,m+1}(b,\lambda^{\ast})}{\partial\lambda}=0. $$
(2.33)

Combining (2.33) and (2.29), with \(\lambda=\lambda^{\ast}\), we obtain

$$\begin{aligned}& \frac{1}{r_{1}} \biggl[ \int_{a}^{c_{1}} \bigl(\bigl\vert \varphi_{1,1}(x,\lambda)\bigr\vert ^{2}+\bigl\vert \varphi _{2,1}(x,\lambda)\bigr\vert ^{2} \bigr) \,dx \biggr] \\& \quad {} +\sum_{k=2}^{m}\frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \biggl[ \int_{c_{k-1}}^{c_{k}} \bigl(\bigl\vert \varphi_{1,k}(x,\lambda)\bigr\vert ^{2} +\bigl\vert \varphi_{2,k}(x,\lambda)\bigr\vert ^{2} \bigr) \,dx \biggr] \\& \quad {} +\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \biggl[ \int_{c_{m}}^{b} \bigl(\bigl\vert \varphi_{1,m+1}(x,\lambda)\bigr\vert ^{2} +\bigl\vert \varphi_{2,m+1}(x,\lambda)\bigr\vert ^{2} \bigr) \,dx \biggr]=0. \end{aligned}$$
(2.34)

It follows that \(\varphi_{1}(x,\lambda^{\ast})=\varphi_{2}(x,\lambda^{\ast})=0\), which is impossible. This proves the lemma. □

Here \(\{\varphi(\cdot,\lambda_{n})\}_{n=-\infty}^{\infty}\) will be a sequence of vector-valued eigenfunctions of (2.1)-(2.4) corresponding to the eigenvalues \(\{\lambda_{n}\}_{n=-\infty}^{\infty}\). Since \(\chi(\cdot,\lambda)\) satisfies (2.3) and (2.4), then the eigenvalues are also determined via

$$ \sin\alpha\chi_{1,1}(a,\lambda)+\cos\alpha\chi_{2,1}(a,\lambda )=\omega(\lambda). $$
(2.35)

Therefore \(\{\chi(\cdot,\lambda_{n})\}_{n=-\infty}^{\infty}\) is another set of vector-valued eigenfunctions which is related by \(\{\varphi(\cdot,\lambda_{n})\}_{n=-\infty}^{\infty}\) with

$$ \chi(x,\lambda_{n})=\tau_{n}\varphi(x, \lambda_{n}),\quad x\in[a,c_{1})\cup(c_{1},c_{2}) \cup(c_{2},c_{3})\cup\cdots\cup (c_{m},b], n \in\mathbb{Z}, $$
(2.36)

where \(\tau_{n}\neq0\) are non-zero constants, since all eigenvalues are simple. Since the eigenvalues are all real, we can take the vector-valued eigenfunctions to be real valued.

3 Asymptotic formulas of eigenvalues and eigenvector functions

In this section, we derive the asymptotic formulas of the eigenvalues \(\{\lambda_{n}\}_{n=-\infty}^{\infty}\) and the vector-valued eigenfunctions \(\{\varphi(\cdot,\lambda_{n})\}_{n=-\infty}^{\infty}\). We shall transform equations (2.1), (2.13), (2.16), and (2.20) into the integral equations, see [2],

$$\begin{aligned}& \begin{aligned}[b] \varphi_{1,1}(x,\lambda)={}&\cos \biggl[ \frac{\lambda( x-a)}{r_{1}}-\alpha \biggr]-\frac{1}{r_{1}} \int_{a}^{x}\sin \biggl[\frac{\lambda(x-t)}{r_{1}} \biggr]q_{1}(t)\varphi_{1,1}(t,\lambda ) \,dt \\ &{}-\frac{1}{r_{1}} \int_{a}^{x}\cos \biggl[\frac{\lambda (x-t)}{r_{1}} \biggr]q_{2}(t) \varphi_{2,1}(t,\lambda) \,dt, \end{aligned} \end{aligned}$$
(3.1)
$$\begin{aligned}& \begin{aligned}[b] \varphi_{2,1}(x,\lambda)={}&\sin \biggl[ \frac{\lambda( x-a)}{r_{1}}-\alpha \biggr]+\frac{1}{r_{1}} \int_{a}^{x}\cos \biggl[\frac{\lambda(x-t)}{r_{1}} \biggr]q_{1}(t)\varphi_{1,1}(t,\lambda ) \,dt \\ &{}-\frac{1}{r_{1}} \int_{a}^{x}\sin \biggl[\frac{\lambda (x-t)}{r_{1}} \biggr]q_{2}(t) \varphi_{2,1}(t,\lambda) \,dt, \end{aligned} \end{aligned}$$
(3.2)
$$\begin{aligned}& \begin{aligned}[b] \varphi_{1,k+1}(x,\lambda)={}&\frac{\gamma_{k}}{\delta_{k}} \varphi _{1,k}\bigl(c_{k}^{-},\lambda \bigr)\cos \biggl[\frac{\lambda(x-c_{k})}{r_{k+1}} \biggr]-\frac{\gamma _{k}^{\prime}}{\delta_{k}^{\prime}}\varphi_{2,k} \bigl(c_{k}^{-},\lambda \bigr)\sin \biggl[\frac{\lambda(x-c_{k})}{r_{k+1}} \biggr] \\ &{}-\frac{1}{r_{k+1}} \int_{c_{k}}^{x}\sin \biggl[\frac{\lambda (x-t)}{r_{k+1}} \biggr]q_{1}(t) \varphi_{1,k+1}(t,\lambda) \,dt \\ &{}- \frac{1}{r_{k+1}} \int_{c_{k}}^{x}\cos \biggl[\frac{\lambda (x-t)}{r_{k+1}} \biggr]q_{2}(t)\varphi_{2,k+1}(t,\lambda) \,dt, \end{aligned} \end{aligned}$$
(3.3)
$$\begin{aligned}& \begin{aligned}[b] \varphi_{2,k+1}(x,\lambda)={}&\frac{\gamma_{k}}{\delta_{k}} \varphi _{1,k}\bigl(c_{k}^{-},\lambda \bigr)\sin \biggl[\frac{\lambda(x-c_{k})}{r_{k+1}} \biggr]+\frac{\gamma _{k}^{\prime}}{\delta_{k}^{\prime}}\varphi_{2,k} \bigl(c_{k}^{-},\lambda \bigr)\cos \biggl[\frac{\lambda(x-c_{k})}{r_{k+1}} \biggr] \\ &{}+\frac{1}{r_{k+1}} \int_{c_{k}}^{x}\cos \biggl[\frac{\lambda (x-t)}{r_{k+1}} \biggr]q_{1}(t) \varphi_{1,k+1}(t,\lambda) \,dt \\ &{}- \frac{1}{r_{k+1}} \int_{c_{k}}^{x}\sin \biggl[\frac{\lambda (x-t)}{r_{k+1}} \biggr]q_{2}(t)\varphi_{2,k+1}(t,\lambda) \,dt, \end{aligned} \end{aligned}$$
(3.4)

where \(k=1,2,\ldots,m\).

For \(|\lambda|\rightarrow \infty\) the following estimates hold uniformly with respect to x, \(x\in[a,c_{1})\cup(c_{1},c_{2})\cup (c_{2},c_{3})\cup\cdots\cup(c_{m},b]\), cf. [1], p.55:

$$\begin{aligned}& \varphi_{1,1}(x,\lambda)=\cos \biggl[\frac{\lambda( x-a)}{r_{1}}- \alpha \biggr]+\mathcal{O} \biggl(\frac{1}{\lambda } \biggr), \end{aligned}$$
(3.5)
$$\begin{aligned}& \varphi_{2,1}(x,\lambda)=\sin \biggl[\frac{\lambda( x-a)}{r_{1}}- \alpha \biggr]+\mathcal{O} \biggl(\frac{1}{\lambda } \biggr), \end{aligned}$$
(3.6)
$$\begin{aligned}& \begin{aligned}[b] \varphi_{1,k+1}(x,\lambda)={}&\frac{\gamma_{k}}{\delta_{k}}\varphi _{1,k}\bigl(c_{k}^{-},\lambda \bigr)\cos \biggl[ \frac{\lambda(x-c_{k})}{r_{k+1}} \biggr]-\frac{\gamma _{k}^{\prime}}{\delta_{k}^{\prime}}\varphi_{2,k} \bigl(c_{k}^{-},\lambda \bigr)\sin \biggl[\frac{\lambda(x-c_{k})}{r_{k+1}} \biggr] \\ &{}+\mathcal {O} \biggl(\frac{1}{\lambda} \biggr), \quad k=1,2,\ldots,m, \end{aligned} \end{aligned}$$
(3.7)
$$\begin{aligned}& \begin{aligned}[b] \varphi_{2,k+1}(x,\lambda)={}&\frac{\gamma_{k}}{\delta_{k}}\varphi _{1,k}\bigl(c_{k}^{-},\lambda \bigr)\sin \biggl[ \frac{\lambda(x-c_{k})}{r_{k+1}} \biggr]+\frac{\gamma _{k}^{\prime}}{\delta_{k}^{\prime}}\varphi_{2,k} \bigl(c_{k}^{-},\lambda \bigr)\cos \biggl[\frac{\lambda(x-c_{k})}{r_{k+1}} \biggr] \\ &{}+\mathcal {O} \biggl(\frac{1}{\lambda} \biggr), \quad k=1,2,\ldots,m. \end{aligned} \end{aligned}$$
(3.8)

Now we will find an asymptotic formula for the eigenvalues. Let \(\gamma _{k}\delta_{k}^{\prime}-\gamma_{k}^{\prime}\delta_{k}=0\), \(k=1,2,\ldots,m\). Since the eigenvalues of the boundary value problem (2.1)-(2.4) coincide with the roots of the equation

$$ \prod_{k=1}^{m}D_{k} \bigl[\sin\beta\varphi_{1,m+1}(b,\lambda)+ \cos\beta\varphi_{2,m+1}(b, \lambda)\bigr]=0, $$
(3.9)

then from the estimates (3.7) and (3.8) (with \(k=m\)) in (3.9) we get

$$\begin{aligned}& \sin\beta \biggl[\frac{\gamma_{m}}{\delta_{m}}\varphi _{1,m}\bigl(c_{m}^{-}, \lambda \bigr)\cos \biggl[\frac{\lambda(b-c_{m})}{r_{m+1}} \biggr]-\frac{\gamma _{m}^{\prime}}{\delta_{m}^{\prime}} \varphi_{2,m}\bigl(c_{m}^{-},\lambda \bigr)\sin \biggl[\frac{\lambda(b-c_{m})}{r_{m+1}} \biggr] \biggr] \\& \quad {} +\cos\beta \biggl[\frac{\gamma_{m}}{\delta_{m}}\varphi _{1,m} \bigl(c_{m}^{-},\lambda \bigr)\sin \biggl[\frac{\lambda(b-c_{m})}{r_{m+1}} \biggr]+\frac{\gamma _{m}^{\prime}}{\delta_{m}^{\prime}}\varphi_{2,m}\bigl(c_{m}^{-}, \lambda \bigr)\cos \biggl[\frac{\lambda(b-c_{m})}{r_{m+1}} \biggr] \biggr] \\& \quad {}+\mathcal{O} \biggl( \frac{1}{\lambda} \biggr)=0, \end{aligned}$$

which can be written as

$$\begin{aligned}& \frac{\gamma_{m}}{\delta_{m}}\varphi_{1,m}\bigl(c_{m}^{-}, \lambda \bigr)\sin \biggl[\frac{\lambda(b-c_{m})}{r_{m+1}}+\beta \biggr]+\frac {\gamma_{m}^{\prime}}{\delta_{m}^{\prime}} \varphi _{2,m}\bigl(c_{m}^{-},\lambda \bigr)\cos \biggl[\frac{\lambda(b-c_{m})}{r_{m+1}}+\beta \biggr] \\& \quad {}+\mathcal{O} \biggl(\frac{1}{\lambda} \biggr)=0. \end{aligned}$$
(3.10)

Then, from (3.5), (3.6), (3.7), (3.8), and \(\gamma_{k}\delta_{k}^{\prime}-\gamma_{k}^{\prime }\delta_{k}=0\), equation (3.10) has the form

$$\begin{aligned}& \sin \biggl[\frac{\lambda(c_{1}-a)}{r_{1}}+\frac{\lambda (c_{2}-c_{1})}{r_{2}}+ \frac{\lambda(c_{3}-c_{2})}{r_{3}}+\cdots +\frac{\lambda(c_{m}-c_{m-1})}{r_{m}}+\frac{\lambda (b-c_{m})}{r_{m+1}}-\alpha+\beta \biggr] \\& \quad {}+\mathcal{O} \biggl(\frac {1}{\lambda} \biggr)=0, \end{aligned}$$
(3.11)

which can be rewritten as

$$ \sin [\Theta\lambda-\alpha+\beta ]+\mathcal{O} \biggl( \frac{1}{\lambda} \biggr)=0,\quad \Theta:=\sum_{k=0}^{m} \frac {(c_{k+1}-c_{k})}{r_{k+1}}. $$
(3.12)

For large \(|\lambda|\) equation (3.12) obviously has solutions which, as is not hard to see, have the form, cf. [1], p.57,

$$ \Theta\lambda_{n}-\alpha+\beta=n\pi+ \delta_{n},\quad n=0,\pm1,\pm2,\ldots. $$
(3.13)

Inserting these values in (3.12), we find that \(\sin\delta_{n}=\mathcal{O} (\frac{1}{n} )\), i.e., \(\delta_{n}=\mathcal{O} (\frac{1}{n} )\). Thus we obtain the following asymptotic formula for the eigenvalues:

$$ \lambda_{n}=\frac{n\pi +\alpha-\beta}{\Theta}+\mathcal{O} \biggl( \frac{1}{n} \biggr),\quad n=0,\pm1,\pm2,\ldots. $$
(3.14)

Using equations (3.14), we obtain the following asymptotic formulas for the vector-valued eigenfunctions \(\boldsymbol{\varphi}(\cdot,\lambda_{n} )\), where \(\gamma _{k}\delta_{k}^{\prime}-\gamma_{k}^{\prime}\delta_{k}=0\):

$$ \varphi(x,\lambda_{n})=\left \{ \textstyle\begin{array}{l@{\quad}l} \Bigl( {\scriptsize\begin{matrix}{} { \varphi_{1,1}(x,\lambda_{n})=\cos [\frac {\lambda_{n}( x-a)}{r_{1}}-\alpha ]+\mathcal{O} (\frac{1}{n} )} \cr { \varphi_{2,1}(x,\lambda_{n})=\sin [\frac {\lambda_{n}( x-a)}{r_{1}}-\alpha ]+\mathcal{O} (\frac{1}{n} ) } \end{matrix}} \Bigr), & x\in[c_{0},c_{1}), \\ \Bigl( {\scriptsize\begin{matrix}{} { \varphi_{1,k+1}(x,\lambda_{n})=\prod_{i=1}^{k}\frac {\gamma_{k}}{\delta_{k}}\cos [\lambda_{n}\sum_{i=1}^{k} \frac{(c_{i}-c_{i-1})}{r_{i}}+\frac{\lambda _{n}(x-c_{k})}{r_{k+1}}-\alpha ]+\mathcal{O} (\frac {1}{n} ) } \cr { \varphi_{2,k+1}(x,\lambda_{n})=\prod_{i=1}^{k}\frac {\gamma_{k}}{\delta_{k}}\sin [\lambda_{n}\sum_{i=1}^{k} \frac{(c_{i}-c_{i-1})}{r_{i}}+\frac{\lambda _{n}(x-c_{k})}{r_{k+1}}-\alpha ]+\mathcal{O} (\frac {1}{n} ) } \end{matrix}} \Bigr), & x\in(c_{k},c_{k+1}), \end{array}\displaystyle \right . $$
(3.15)

where \(k=1,2,\ldots,m\).

4 Green’s matrix and expansion theorem

Let \(f(\cdot)=\bigl ( {\scriptsize\begin{matrix}{} f_{1}(\cdot)\cr f_{2}(\cdot) \end{matrix}} \bigr )\) be a continuous vector-valued function. Now, we derive the Green’s matrix of problem (2.1)-(2.4). Consider the inhomogeneous eigenvalue problem consisting of the differential system

$$\begin{aligned}& \left . \textstyle\begin{array}{l} r(x)u_{2}^{\prime}(x)- \{q_{1}(x)+\lambda \} u_{1}(x)=f_{1}(x), \\ r(x)u_{1}^{\prime}(x)+ \{q_{2}(x)+\lambda \} u_{2}(x)=-f_{2}(x), \end{array}\displaystyle \right \} \\& \quad x\in[a,c_{1})\cup(c_{1},c_{2}) \cup(c_{2},c_{3})\cup\cdots\cup (c_{m},b], \end{aligned}$$
(4.1)

and the boundary conditions (2.2)-(2.4) with λ is not an eigenvalue of problem (2.1)-(2.4).

Now, we can represent the general solution of (4.1) in the following form:

$$ u(x,\lambda)=\left \{ \textstyle\begin{array}{l@{\quad}l} \mathcal{A}_{1}(x,\lambda)\bigl( {\scriptsize\begin{matrix}{} \varphi_{1,1}(x,\lambda) \cr \varphi_{2,1}(x,\lambda) \end{matrix}} \bigr)+\mathcal{B}_{1}(x,\lambda)\bigl( {\scriptsize\begin{matrix}{} \chi_{1,1}(x,\lambda) \cr \chi_{2,1}(x,\lambda) \end{matrix}} \bigr), & x\in[a,c_{1}), \\ \mathcal{A}_{k}(x,\lambda)\bigl( {\scriptsize\begin{matrix}{} \varphi_{1,k}(x,\lambda) \cr \varphi_{2,k}(x,\lambda) \end{matrix}} \bigr)+\mathcal{B}_{k}(x,\lambda)\bigl( {\scriptsize\begin{matrix}{} \chi_{1,k}(x,\lambda) \cr \chi_{2,k}(x,\lambda) \end{matrix}} \bigr), & x\in(c_{k-1}, c_{k})\ (k=2,3, \ldots,m), \\ \mathcal{A}_{m+1}(x,\lambda)\bigl( {\scriptsize\begin{matrix}{} \varphi_{1,m+1}(x,\lambda) \cr \varphi_{2,m+1}(x,\lambda) \end{matrix}} \bigr) \\ \quad {}+\mathcal{B}_{m+1}(x,\lambda)\bigl( {\scriptsize\begin{matrix}{} \chi_{1,m+1}(x,\lambda) \cr \chi_{2,m+1}(x,\lambda) \end{matrix}} \bigr), & x\in(c_{m},b]. \end{array}\displaystyle \right . $$
(4.2)

We applied the standard method of variation of the constants to (4.2), thus, the functions \(\mathcal{A}_{k}(x,\lambda)\) and \(\mathcal{B}_{k}(x,\lambda)\) (\(k=1,2,\ldots,m+1\)) satisfy the linear system of equations

$$\begin{aligned}& \left . \textstyle\begin{array}{l} \mathcal{A}_{1}^{\prime}(x,\lambda)\varphi_{2,1}(x,\lambda) +\mathcal{B}_{1}^{\prime}(x,\lambda)\chi_{2,1}(x,\lambda)=\frac {f_{1}(x)}{r_{1}}, \\ \mathcal{A}_{1}^{\prime}(x,\lambda)\varphi_{1,1}(x,\lambda) +\mathcal{B}_{1}^{\prime}(x,\lambda)\chi_{1,1}(x,\lambda)=-\frac {f_{2}(x)}{r_{1}}, \end{array}\displaystyle \right \} \quad x\in[a,c_{1}), \end{aligned}$$
(4.3)
$$\begin{aligned}& \left . \textstyle\begin{array}{l} \mathcal{A}_{k}^{\prime}(x,\lambda)\varphi_{2,k}(x,\lambda) +\mathcal{B}_{k}^{\prime}(x,\lambda)\chi_{2,k}(x,\lambda)=\frac {f_{1}(x)}{r_{k}}, \\ \mathcal{A}_{k}^{\prime}(x,\lambda)\varphi_{1,k}(x,\lambda) +\mathcal{B}_{k}^{\prime}(x,\lambda)\chi_{1,k}(x,\lambda)=-\frac {f_{2}(x)}{r_{k}}, \end{array}\displaystyle \right \} \quad x\in(c_{k-1},c_{k})\ (k=2,3,\ldots,m), \end{aligned}$$
(4.4)

and

$$ \left . \textstyle\begin{array}{l} \mathcal{A}_{m+1}^{\prime}(x,\lambda)\varphi_{2,m+1}(x,\lambda) +\mathcal{B}_{m+1}^{\prime}(x,\lambda)\chi_{2,m+1}(x,\lambda )=\frac{f_{1}(x)}{r_{m+1}}, \\ \mathcal{A}_{m+1}^{\prime}(x,\lambda)\varphi_{1,m+1}(x,\lambda) +\mathcal{B}_{m+1}^{\prime}(x,\lambda)\chi_{1,m+1}(x,\lambda )=-\frac{f_{2}(x)}{r_{m+1}}, \end{array}\displaystyle \right \} \quad x\in (c_{m},b]. $$
(4.5)

Since λ is not an eigenvalue and \(\omega(\lambda)\neq0\), each of the linear system in (4.3), (4.4), and (4.5) has a unique solution which leads to

$$\begin{aligned}& \left . \textstyle\begin{array}{l} \mathcal{A}_{1}(x,\lambda)=\frac{1}{r_{1}\omega_{1}(\lambda )} \int_{x}^{c_{1}}\chi^{\top}(\xi,\lambda)f(\xi) \, d\xi+\mathcal{A}_{1}, \\ \mathcal{B}_{1}(x,\lambda)=\frac{1}{r_{1}\omega_{1}(\lambda )} \int_{a}^{x}\varphi^{\top}(\xi,\lambda)f(\xi) \, d\xi+\mathcal{B}_{1}, \end{array}\displaystyle \right \} \quad x\in[a, c_{1}), \end{aligned}$$
(4.6)
$$\begin{aligned}& \left . \textstyle\begin{array}{l} \mathcal{A}_{k}(x,\lambda)=\frac{1}{r_{k}\omega_{k}(\lambda )} \int_{x}^{c_{k}}\chi^{\top}(\xi,\lambda)f(\xi) \, d\xi+\mathcal{A}_{k}, \\ \mathcal{B}_{k}(x,\lambda)=\frac{1}{r_{k}\omega_{k}(\lambda )} \int_{c_{k-1}}^{x}\varphi^{\top}(\xi,\lambda)f(\xi )\, d\xi+\mathcal{B}_{k}, \end{array}\displaystyle \right \} \quad x\in(c_{k-1}, c_{k})\ (k=2,3,\ldots,m), \end{aligned}$$
(4.7)
$$\begin{aligned}& \left . \textstyle\begin{array}{l} \mathcal{A}_{m+1}(x,\lambda)=\frac{1}{r_{m+1}\omega_{m+1}(\lambda )} \int_{x}^{b}\chi^{\top}(\xi,\lambda)f(\xi) \, d\xi+\mathcal{A}_{m+1}, \\ \mathcal{B}_{m+1}(x,\lambda)=\frac{1}{r_{m+1}\omega_{m+1}(\lambda )} \int_{c_{m}}^{x}\varphi^{\top}(\xi,\lambda)f(\xi ) \, d\xi+\mathcal{B}_{m+1}, \end{array}\displaystyle \right \}\quad x\in(c_{m}, b], \end{aligned}$$
(4.8)

where \(\mathcal{A}_{k}\) and \(\mathcal{B}_{k}\), \(k=1,2,\ldots,m+1\), are arbitrary constants, and

$$\begin{aligned}& \varphi(\xi,\lambda)=\left \{ \textstyle\begin{array}{l@{\quad}l} \bigl( {\scriptsize\begin{matrix}{} \varphi_{1,1}(\xi,\lambda) \cr \varphi_{2,1}(\xi,\lambda) \end{matrix}} \bigr), & \xi\in [a,c_{1}), \\ \bigl( {\scriptsize\begin{matrix}{} \varphi_{1,k}(\xi,\lambda) \cr \varphi_{2,k}(\xi,\lambda) \end{matrix}} \bigr), & \xi\in(c_{k-1},c_{k}) \ (k=2,3,\ldots,m), \\ \bigl( {\scriptsize\begin{matrix}{} \varphi_{1,m+1}(\xi,\lambda) \cr \varphi_{2,m+1}(\xi,\lambda) \end{matrix}} \bigr), & \xi\in(c_{m},b], \end{array}\displaystyle \right . \\& \chi(\xi,\lambda)=\left \{ \textstyle\begin{array}{l@{\quad}l} \bigl( {\scriptsize\begin{matrix}{} \chi_{1,1}(\xi,\lambda) \cr \chi_{2,1}(\xi,\lambda) \end{matrix}} \bigr), & \xi\in [a,c_{1}), \\ \bigl( {\scriptsize\begin{matrix}{} \chi_{1,k}(\xi,\lambda) \cr \chi_{2,k}(\xi,\lambda) \end{matrix}} \bigr), & \xi\in(c_{k-1},c_{k}) \ (k=2,3,\ldots,m), \\ \bigl( {\scriptsize\begin{matrix}{} \chi_{1,m+1}(\xi,\lambda) \cr \chi_{2,m+1}(\xi,\lambda) \end{matrix}} \bigr), & \xi\in(c_{m},b]. \end{array}\displaystyle \right . \end{aligned}$$

Substituting equations (4.6), (4.7), and (4.8) into (4.2), we obtain the solution of (4.1),

$$ u(x,\lambda)=\left \{ \textstyle\begin{array}{l} \frac{\varphi(x,\lambda)}{r_{1}\omega_{1}(\lambda)} \int_{x}^{c_{1}}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi+\frac{\chi(x,\lambda)}{r_{1}\omega_{1}(\lambda)} \int_{a}^{x}\varphi^{\top}(\xi,\lambda)f(\xi) \,d\xi \\ \quad {}+ \mathcal{A}_{1}\varphi(x,\lambda)+\mathcal{B}_{1}\chi(x,\lambda),\quad x\in[a,c_{1}), \\ \frac{\varphi(x,\lambda)}{r_{k}\omega_{k}(\lambda)} \int_{x}^{c_{k}}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi+\frac{\chi(x,\lambda)}{r_{k}\omega_{k}(\lambda)} \int_{c_{k-1}}^{x}\varphi^{\top}(\xi,\lambda)f(\xi )\,d\xi \\ \quad {}+ \mathcal{A}_{k}\varphi(x,\lambda)+\mathcal{B}_{k}\chi(x,\lambda), \quad x\in(c_{k-1}, c_{k})\ (k=2,3,\ldots,m), \\ \frac{\varphi(x,\lambda)}{r_{m+1}\omega_{m+1}(\lambda)} \int_{x}^{b}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi+\frac{\chi(x,\lambda)}{r_{m+1}\omega_{m+1}(\lambda)} \int_{c_{m}}^{x}\varphi^{\top}(\xi,\lambda)f(\xi) \,d\xi \\ \quad {}+ \mathcal{A}_{m+1}\varphi(x,\lambda) +\mathcal{B}_{m+1}\chi(x,\lambda), \quad x\in(c_{m}, b]. \end{array}\displaystyle \right . $$
(4.9)

Then from (2.2), (2.3), and (2.4), we get

$$\begin{aligned}& \left \{ \textstyle\begin{array}{l} \mathcal{A}_{1}=\frac{1}{r_{2}\omega_{2}(\lambda)} \int_{c_{1}}^{c_{2}}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi+\frac{1}{r_{3}\omega_{3}(\lambda)} \int _{c_{2}}^{c_{3}}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi \\ \hphantom{\mathcal{A}_{1}={}}{}+\cdots+\frac{1}{r_{m+1}\omega_{m+1}(\lambda)} \int_{c_{m}}^{b}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi, \\ \mathcal{A}_{2}=\frac{1}{r_{3}\omega_{3}(\lambda)} \int_{c_{2}}^{c_{3}}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi+\frac{1}{r_{4}\omega_{4}(\lambda)} \int _{c_{3}}^{c_{4}}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi \\ \hphantom{\mathcal{A}_{2}={}}{}+\cdots+\frac{1}{r_{m+1}\omega_{m+1}(\lambda)} \int_{c_{m}}^{b}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi, \\ \ldots, \\ \mathcal{A}_{m}=\frac{1}{r_{m+1}\omega_{m+1}(\lambda )} \int_{c_{m}}^{b}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi, \\ \mathcal{A}_{m+1}=0, \end{array}\displaystyle \right . \end{aligned}$$
(4.10)
$$\begin{aligned}& \left \{ \textstyle\begin{array}{l} \mathcal{B}_{1}=0, \\ \mathcal{B}_{2}=\frac{1}{r_{1}\omega_{1}(\lambda)} \int_{a}^{c_{1}}\varphi^{\top}(\xi,\lambda)f(\xi) \,d\xi, \\ \mathcal{B}_{3}=\frac{1}{r_{1}\omega_{1}(\lambda)} \int_{a}^{c_{1}}\varphi^{\top}(\xi,\lambda)f(\xi) \,d\xi+\frac{1}{r_{2}\omega_{2}(\lambda)} \int _{c_{1}}^{c_{2}}\varphi^{\top}(\xi,\lambda)f(\xi) \,d\xi, \\ \ldots, \\ \mathcal{B}_{m+1}=\frac{1}{r_{1}\omega_{1}(\lambda)} \int_{a}^{c_{1}}\varphi^{\top}(\xi,\lambda)f(\xi) \,d\xi+\frac{1}{r_{2}\omega_{2}(\lambda)} \int _{c_{1}}^{c_{2}}\varphi^{\top}(\xi,\lambda)f(\xi) \,d\xi \\ \hphantom{\mathcal{B}_{m+1}={}}{}+\cdots+\frac{1}{r_{m}\omega_{m}(\lambda)} \int _{c_{m-1}}^{c_{m}}\varphi^{\top}(\xi,\lambda)f(\xi) \,d\xi. \end{array}\displaystyle \right . \end{aligned}$$
(4.11)

Equations (4.10) and (4.11) can be rewritten as follows:

$$ \mathcal{A}_{k}:=\left \{ \textstyle\begin{array}{l@{\quad}l} \sum_{i=k}^{m} \frac{1}{r_{i+1}\omega_{i+1}(\lambda )} \int_{c_{i}}^{c_{i+1}}\chi^{\top}(\xi,\lambda )f(\xi) \,d\xi, & k=1,2,\ldots,m , \\ 0, & k=m+1 \end{array}\displaystyle \right . $$
(4.12)

and

$$ \mathcal{B}_{k}:=\left \{ \textstyle\begin{array}{l@{\quad}l} 0, & k=1 , \\ \sum_{i=1}^{k-1}\frac{1}{r_{i}\omega_{i}(\lambda )} \int_{c_{i-1}}^{c_{i}}\varphi^{\top}(\xi,\lambda )f(\xi) \,d\xi, & k=2,3,\ldots,m+1 , \end{array}\displaystyle \right . $$
(4.13)

respectively. Then, from (4.12) and (4.13), (4.9) can be written as

$$ u(x,\lambda)=\left \{ \textstyle\begin{array}{l} \frac{\varphi(x,\lambda)}{r_{1}\omega_{1}(\lambda)} \int_{x}^{c_{1}}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi+\frac{\chi(x,\lambda)}{r_{1}\omega_{1}(\lambda)} \int_{a}^{x}\varphi^{\top}(\xi,\lambda)f(\xi) \,d\xi \\ \quad {}+ \sum_{i=1}^{m} \frac{\varphi(x,\lambda )}{r_{i+1}\omega_{i+1}(\lambda)} \int _{c_{i}}^{c_{i+1}}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi, \quad x\in[a,c_{1}), \\ \frac{\varphi(x,\lambda)}{r_{k}\omega_{k}(\lambda)} \int_{x}^{c_{k}}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi+\frac{\chi(x,\lambda)}{r_{k}\omega_{k}(\lambda)} \int_{c_{k-1}}^{x}\varphi^{\top}(\xi,\lambda)f(\xi )\,d\xi \\ \quad {}+ \sum_{i=k}^{m} \frac{\varphi(x,\lambda )}{r_{i+1}\omega_{i+1}(\lambda)} \int _{c_{i}}^{c_{i+1}}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi \\ \quad {}+ \sum_{i=1}^{k-1}\frac{\chi(x,\lambda)}{r_{i}\omega _{i}(\lambda)} \int_{c_{i-1}}^{c_{i}}\varphi^{\top }(\xi,\lambda)f(\xi) \,d\xi, \quad x\in(c_{k-1}, c_{k}) \ (k=2,3,\ldots,m), \\ \frac{\varphi(x,\lambda)}{r_{m+1}\omega_{m+1}(\lambda)} \int_{x}^{b}\chi^{\top}(\xi,\lambda)f(\xi) \,d\xi+\frac{\chi(x,\lambda)}{r_{m+1}\omega_{m+1}(\lambda)} \int_{c_{m}}^{x}\varphi^{\top}(\xi,\lambda)f(\xi) \,d\xi \\ \quad {}+ \sum_{i=1}^{m}\frac{\chi(x,\lambda)}{r_{i}\omega _{i}(\lambda)} \int_{c_{i-1}}^{c_{i}}\varphi^{\top }(\xi,\lambda)f(\xi) \,d\xi,\quad x\in(c_{m}, b], \end{array}\displaystyle \right . $$
(4.14)

which can be written as

$$\begin{aligned} u(x,\lambda) =& \frac{1}{r_{1}} \int^{c_{1}}_{a}G(x,\xi,\lambda)f(\xi) \,d\xi+\sum _{k=2}^{m}\frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}G(x,\xi,\lambda)f(\xi) \,d\xi \\ &{}+ \frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}G(x,\xi ,\lambda)f(\xi) \,d\xi, \end{aligned}$$
(4.15)

where

$$ G(x,\xi,\lambda)=\frac{1}{\omega(\lambda)} \left \{\textstyle\begin{array}{l@{\quad}l} \chi(x,\lambda)\varphi^{\top}(\xi,\lambda), & a\le\xi\le x\le b, x,\xi\neq c_{k}, k=1,2,\ldots,m, \\ \varphi(x,\lambda)\chi^{\top}(\xi,\lambda), & a\le x\le\xi\le b, x,\xi\neq c_{k}, k=1,2,\ldots,m. \end{array}\displaystyle \right . $$
(4.16)

Expanding (4.16) we obtain the concrete form

$$\begin{aligned}& G(x,\xi,\lambda) \\& \quad =\frac{1}{\omega(\lambda)} \left \{\textstyle\begin{array}{l@{\quad}l} \bigl( {\scriptsize\begin{matrix}{} \varphi_{1}(\xi,\lambda)\chi_{1}(x,\lambda) & \varphi_{2}(\xi ,\lambda) \chi_{1}(x,\lambda) \cr \varphi_{1}(\xi,\lambda)\chi_{2}(x,\lambda) & \varphi_{2}(\xi ,\lambda) \chi_{2}(x,\lambda) \end{matrix}} \bigr), & a\le\xi\le x\le b, x,\xi\neq c_{k}, k=1,2,\ldots,m, \\ \bigl( {\scriptsize\begin{matrix}{} \varphi_{1}(x,\lambda)\chi_{1}(\xi,\lambda) & \varphi _{1}(x,\lambda) \chi_{2}(\xi,\lambda) \cr \varphi_{2}(x,\lambda)\chi_{1}(\xi,\lambda) & \varphi _{2}(x,\lambda) \chi_{2}(\xi,\lambda) \end{matrix}} \bigr), & a\le x\le\xi\le b, x,\xi \neq c_{k}, k=1,2,\ldots,m. \end{array}\displaystyle \right . \end{aligned}$$
(4.17)

The matrix \(G(x,\xi,\lambda)\) is called the Green’s matrix of problem (2.1)-(2.4). Obviously \(G(x,\xi,\lambda)\) is a meromorphic function of λ, for every \((x,\xi)\in ([a,c_{1})\cup (c_{1},c_{2})\cup(c_{2},c_{3})\cup\cdots\cup(c_{m},b] )^{2}\), which has simple poles only at the eigenvalues. Although the Green’s function looks as simple as that of Dirac systems, cf. e.g. [1, 2], it is rather complicated because of the transmission conditions; see the example at the end of this paper.

The next theorem is a vector-valued eigenfunction expansion theorem. The proof is exactly similar to that of Levitan and Sargsjan derived in [1], pp.67-77; see also [2].

Theorem 4.1

  1. (i)

    For \(u(\cdot)\in\mathcal{H}\),

    $$ \bigl\| u(\cdot)\bigr\| _{\mathcal{H}}^{2}=\sum _{n=-\infty}^{\infty}\bigl\vert \bigl\langle u(\cdot), \psi_{n}(\cdot)\bigr\rangle _{\mathcal{H}}\bigr\vert ^{2}. $$
    (4.18)
  2. (ii)

    For \(u(\cdot)\in\mathcal{H}\),

    $$ u(x)=\sum_{n=-\infty}^{\infty}\bigl\langle u(\cdot),\psi_{n}(\cdot) \bigr\rangle _{\mathcal{H}} \psi_{n}(x), $$
    (4.19)

    where \(\psi_{n}(\cdot)=\frac{\varphi(\cdot,\lambda_{n})}{\|\varphi (\cdot,\lambda_{n})\|_{\mathcal{H}}}\) is the normalized vector-valued eigenfunction of the problem (2.1)-(2.4). The series is absolutely and uniformly convergent on \([a,c_{1})\cup(c_{1},c_{2})\cup(c_{2},c_{3})\cup\cdots\cup (c_{m},b]\).

5 The sampling theorems

In this section we introduce two sampling theorems associated with problem (2.1)-(2.4) above. In the first one we define an integral transform whose kernel is \(\varphi(\cdot,\lambda)\), while in the second one the sampled transform will be defined in terms of the Green’s matrix. The first sampling theorem of this section associated with the boundary value problem (2.1)-(2.4) is the following theorem.

Theorem 5.1

Let \(f(x)= \bigl ({\scriptsize\begin{matrix}{} f_{1}(x)\cr f_{2}(x)\end{matrix}}\bigr ) \in\mathcal{H}\). For \(\lambda\in\mathbb{C}\) let

$$\begin{aligned} \mathcal{F}(\lambda) =&\frac{1}{r_{1}} \int^{c_{1}}_{a}{f}^{\top }(x)\varphi(x,\lambda) \,dx+\sum_{k=2}^{m}\frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}{f}^{\top}(x)\varphi(x, \lambda) \,dx \\ &{}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}{f}^{\top }(x)\varphi(x,\lambda) \,dx, \end{aligned}$$
(5.1)

where \(\varphi(\cdot,\lambda)\) is the solution defined in (2.13). Then \(\mathcal{F}(\lambda)\) is an entire function of exponential type that can be reconstructed from its values at the points \(\{\lambda_{n}\}_{n=-\infty}^{\infty}\) via the sampling formula

$$ \mathcal{F}(\lambda)=\sum_{n=-\infty}^{\infty} \mathcal{F}(\lambda _{n})\frac{\omega(\lambda)}{ (\lambda-\lambda_{n})\omega^{\prime}(\lambda_{n})}. $$
(5.2)

The series (5.2) converges absolutely on \(\mathbb{C}\) and uniformly on any compact subset of \(\mathbb{C}\). Here \(\omega(\lambda)\) is the entire function defined in (2.24).

Proof

Equation (5.1) can be rewritten in the form

$$\begin{aligned} \mathcal{F}(\lambda) =&\bigl\langle f(\cdot),\varphi(\cdot,\lambda )\bigr\rangle _{\mathcal{H}} \\ =&\frac{1}{r_{1}} \int^{c_{1}}_{a}{f}^{\top}(x)\varphi(x,\lambda) \,dx+\sum_{k=2}^{m}\frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}{f}^{\top}(x)\varphi(x, \lambda) \,dx \\ &{}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}{f}^{\top }(x)\varphi(x,\lambda) \,dx, \quad \lambda\in\mathbb{C}. \end{aligned}$$
(5.3)

Since both \(f(\cdot)\) and \(\varphi(\cdot,\lambda)\) are in \(\mathcal{H}\), they have the Fourier expansions

$$ f(x)=\sum_{n=-\infty}^{\infty} \hat{{f}}(n)\frac{\varphi (x,\lambda_{n})}{ \|\varphi(\cdot,\lambda_{n})\|^{2}_{\mathcal{H}}},\qquad \varphi(x,\lambda)=\sum _{n=-\infty}^{\infty} \bigl\langle \varphi(\cdot,\lambda), \varphi(\cdot,\lambda_{n}) \bigr\rangle _{\mathcal{H}} \frac{\varphi(x,\lambda_{n})}{\|\varphi(\cdot,\lambda_{n})\| ^{2}_{\mathcal{H}}}, $$
(5.4)

where \(\lambda\in\mathbb{C}\) and \(\hat{{f}}(n)\) are the Fourier coefficients

$$\begin{aligned} \hat{f}(n) =&\bigl\langle f(\cdot),\varphi(\cdot, \lambda_{n})\bigr\rangle _{\mathcal{H}} \\ =&\frac{1}{r_{1}} \int^{c_{1}}_{a}{f}^{\top}(x)\varphi(x,\lambda _{n}) \,dx+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}{f}^{\top}(x)\varphi(x, \lambda_{n}) \,dx \\ &{}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}{f}^{\top }(x)\varphi(x, \lambda_{n}) \,dx. \end{aligned}$$
(5.5)

Applying Parseval’s identity to (5.3), we obtain

$$ \mathcal{F}(\lambda)=\sum_{n=-\infty}^{\infty} \mathcal{F}(\lambda_{n}) \frac{\langle \varphi(\cdot,\lambda), \varphi(\cdot,\lambda_{n})\rangle _{\mathcal{H}}}{\|\varphi (\cdot,\lambda_{n})\|^{2}_{\mathcal{H}}},\quad \lambda\in\mathbb{C}. $$
(5.6)

Now we calculate \(\langle \varphi(\cdot,\lambda),\varphi (\cdot,\lambda_{n})\rangle _{\mathcal{H}}\) and \(\|\varphi(\cdot,\lambda_{n})\|_{\mathcal{H}}\) of \(\lambda\in\mathbb{C}\), \(n\in\mathbb{Z}\). To prove expansion (5.2) we need to show that

$$ \frac{\langle \varphi(\cdot,\lambda),\varphi(\cdot,\lambda _{n})\rangle _{\mathcal{H}}}{ \|\varphi(\cdot,\lambda_{n})\|^{2}_{\mathcal{H}}}= \frac{\omega(\lambda)}{ (\lambda-\lambda_{n})\omega^{\prime}(\lambda)}, \quad n\in \mathbb{Z}, \lambda\in\mathbb{C}. $$
(5.7)

Indeed, let \(\lambda\in\mathbb{C}\) and \(n\in\mathbb{Z}\) be fixed. By the definition of the inner product of \(\mathcal{H}\), we have

$$\begin{aligned} \bigl\langle \varphi(\cdot,\lambda),\varphi(\cdot, \lambda_{n}) \bigr\rangle _{\mathcal{H}} =&\frac{1}{r_{1}} \int^{c_{1}}_{a}\varphi^{\top}(x,\lambda) \varphi (x,\lambda_{n}) \,dx+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}\varphi^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx \\ &{}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}\varphi ^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx. \end{aligned}$$
(5.8)

From the Green’s identity, see [1], p.51, we have

$$\begin{aligned}& (\lambda_{n}-\lambda) \Biggl(\frac{1}{r_{1}} \int ^{c_{1}}_{a}\varphi^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}\varphi^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx \\& \qquad {}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}\varphi ^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx \Biggr) \\& \quad = \varphi_{1,1}\bigl(c_{1}^{-},\lambda\bigr) \varphi_{2,1}\bigl(c_{1}^{-},\lambda_{n} \bigr)- \varphi_{1,1}\bigl(c_{1}^{-},\lambda\bigr) \varphi_{2,1}\bigl(c_{1}^{-},\lambda _{n} \bigr) \\& \qquad {}- \bigl[\varphi_{1,1}(a,\lambda)\varphi_{2,1}(a, \lambda_{n})- \varphi_{1,1}(a,\lambda)\varphi_{2,1}(a, \lambda_{n}) \bigr] \\& \qquad {} +\sum_{k=2}^{m}\prod _{i=1}^{k-1}D_{i} \bigl[\varphi _{1,k}\bigl(c_{k}^{-},\lambda\bigr) \varphi_{2,k}\bigl(c_{k}^{-},\lambda_{n} \bigr)- \varphi_{1,k}\bigl(c_{k}^{-},\lambda\bigr) \varphi_{2,k}\bigl(c_{k}^{-},\lambda _{n} \bigr) \bigr] \\& \qquad {} -\sum_{k=2}^{m}\prod _{i=1}^{k-1}D_{i} \bigl[\varphi _{1,k}\bigl(c_{k-1}^{+},\lambda\bigr) \varphi_{2,k}\bigl(c_{k-1}^{+},\lambda_{n} \bigr)- \varphi_{1,k}\bigl(c_{k-1}^{+},\lambda\bigr) \varphi_{2,k}\bigl(c_{k-1}^{+},\lambda _{n} \bigr) \bigr] \\& \qquad {} +\prod_{i=1}^{m}D_{i} \bigl[\varphi_{1,m+1}(b,\lambda)\varphi _{2,m+1}(b, \lambda_{n})- \varphi_{1,m+1}(b,\lambda)\varphi_{2,m+1}(b, \lambda_{n})\bigr] \\& \qquad {} -\prod_{i=1}^{m}D_{i} \bigl[\varphi_{1,m+1}\bigl(c_{m}^{+},\lambda \bigr) \varphi_{2,m+1}\bigl(c_{m}^{+},\lambda_{n} \bigr)- \varphi_{1,m+1}\bigl(c_{m}^{+},\lambda\bigr) \varphi_{2,m+1}\bigl(c_{m}^{+},\lambda _{n} \bigr) \bigr]. \end{aligned}$$
(5.9)

Then from the initial conditions (2.16) and (2.20), equation (5.9) becomes

$$\begin{aligned}& (\lambda_{n}-\lambda) \Biggl(\frac{1}{r_{1}} \int ^{c_{1}}_{a}\varphi^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}\varphi^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx \\& \qquad {}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}\varphi ^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx \Biggr) \\& \quad = \prod_{i=1}^{m}D_{i} \bigl[\varphi_{1,m+1}(b,\lambda)\varphi _{2,m+1}(b, \lambda_{n})- \varphi_{1,m+1}(b,\lambda)\varphi_{2,m+1}(b, \lambda_{n})\bigr]. \end{aligned}$$
(5.10)

From (2.36) and (2.18), we have

$$\begin{aligned}& \varphi_{1,m+1}(b,\lambda)\varphi_{2,m+1}(b, \lambda_{n})- \varphi_{2,m+1}(b,\lambda)\varphi_{1,m+1}(b, \lambda_{n}) \\& \quad =\tau_{n}^{-1}\bigl[\varphi_{1,m+1}(b,\lambda) \chi_{2,m+1}(b,\lambda_{n})- \varphi_{2,m+1}(b,\lambda) \chi_{1,m+1}(b,\lambda_{n})\bigr] \\& \quad =\tau_{n}^{-1}\bigl[-\sin\beta\varphi_{1,m+1}(b, \lambda)- \cos\beta\varphi_{2,m+1}(b,\lambda)\bigr] \\& \quad =\tau_{n}^{-1}\omega_{m+1}(\lambda). \end{aligned}$$
(5.11)

Substituting from (5.10), (5.11) into (5.8), we get

$$ \bigl\langle \varphi(\cdot,\lambda),\varphi(\cdot, \lambda_{n}) \bigr\rangle _{\mathcal{H}}=\tau_{n}^{-1} \frac{\prod_{i=1}^{m}D_{i} \omega_{m+1}(\lambda )}{\lambda_{n}-\lambda}=\tau_{n}^{-1} \frac{\omega(\lambda)}{\lambda_{n}-\lambda}. $$
(5.12)

Letting \(\lambda\rightarrow\lambda_{n}\) in (5.12) and since the zeros of \(\omega(\lambda)\) are simple, we get

$$ \bigl\langle \varphi(\cdot,\lambda_{n}),\varphi(\cdot, \lambda _{n})\bigr\rangle _{\mathcal{H}}= \bigl\Vert \varphi(\cdot, \lambda_{n})\bigr\Vert ^{2}_{\mathcal{H}} =- \tau_{n}^{-1}\omega^{\prime}(\lambda_{n}). $$
(5.13)

Since \(\lambda\in\mathbb{C}\) and \(n\in\mathbb{Z}\) are arbitrary, (5.12) and (5.13) hold for all \(\lambda\in \mathbb{C}\) and all \(n\in\mathbb{Z}\). Therefore from (5.12) and (5.13) we get (5.7). Hence (5.2) is proved with a pointwise convergence on \(\mathbb{C}\). Now we investigate the convergence of (5.2). First we prove that it is absolutely convergent on \(\mathbb{C}\). Using Cauchy-Schwarz’ inequality for \(\lambda\in\mathbb{C}\),

$$\begin{aligned}& \sum_{k=-\infty}^{\infty}\biggl\vert \mathcal{F}(\lambda_{k}) \frac{\omega(\lambda)}{ (\lambda-\lambda_{k})\omega^{\prime}(\lambda_{k})}\biggr\vert \\ & \quad \le \Biggl(\sum_{k=-\infty}^{\infty} \frac{\vert \langle f(\cdot),\varphi(\cdot,\lambda_{k}) \rangle _{\mathcal{H}} \vert ^{2}}{ \|\varphi(\cdot,\lambda_{k})\|_{\mathcal{H}}^{2}} \Biggr)^{1/2} \Biggl(\sum _{k=-\infty}^{\infty}\frac{\vert \langle \varphi(\cdot ,\lambda), \varphi(\cdot,\lambda_{k})\rangle _{\mathcal{H}}\vert ^{2}}{ \|\varphi(\cdot,\lambda_{k})\|_{\mathcal{H}}^{2}} \Biggr)^{1/2}. \end{aligned}$$
(5.14)

Since \(f(\cdot), \varphi(\cdot,\lambda)\in \mathcal{H}\), the two series in the right-hand side of (5.14) converge. Thus series (5.2) converges absolutely on \(\mathbb{C}\). As for uniform convergence, let \(M\subset\mathbb{C}\) be compact. Let \(\lambda\in M\) and \(N>0\). Define \(\nu_{N}(\lambda)\) to be

$$ \nu_{N}(\lambda):=\Biggl\vert \mathcal{F}(\lambda)- \sum_{k=-N}^{N}\mathcal{F}( \lambda_{k}) \frac{\omega(\lambda)}{ (\lambda-\lambda_{k})\omega^{\prime}(\lambda_{k})}\Biggr\vert . $$
(5.15)

Using the same method as developed above

$$ \nu_{N}(\lambda)\le \Biggl(\sum _{k=-N}^{N}\frac{\vert \langle f(\cdot),\varphi (\cdot,\lambda_{k})\rangle _{\mathcal{H}}\vert ^{2}}{ \|\varphi(\cdot,\lambda_{k})\|_{\mathcal{H}}^{2}} \Biggr)^{1/2} \Biggl(\sum_{k=-N}^{N}\frac{\vert \langle \varphi(\cdot,\lambda ),\varphi (\cdot,\lambda_{k})\rangle _{\mathcal{H}}\vert ^{2}}{ \|\varphi(\cdot,\lambda_{k})\|_{\mathcal{H}}^{2}} \Biggr)^{1/2}. $$
(5.16)

Therefore

$$ \nu_{N}(\lambda)\le \bigl\| \varphi(\cdot, \lambda_{k})\bigr\| _{\mathcal{H}} \Biggl(\sum_{k=-N}^{N} \frac{\vert \langle f(\cdot),\varphi(\cdot,\lambda_{k})\rangle _{\mathcal{H}}\vert ^{2}}{ \|\varphi(\cdot,\lambda_{k})\|_{\mathcal{H}}^{2}} \Biggr)^{1/2}. $$
(5.17)

Since \([a,b]\times M\) is compact, cf. e.g. [28], p.225, we can find a positive constant \(C_{M}\) such that

$$ \bigl\| \varphi(\cdot,\lambda)\bigr\| _{\mathcal{H}}\le C_{M}, \quad \text{for all } \lambda\in M. $$
(5.18)

Then

$$ \nu_{N}(\lambda)\le C_{M} \Biggl(\sum _{k=-N}^{N} \frac{\vert \langle f(\cdot),\varphi(\cdot,\lambda_{k}) \rangle _{\mathcal{H}}\vert ^{2}}{ \|\varphi(\cdot,\lambda_{k})\|_{\mathcal{H}}^{2}} \Biggr)^{1/2} $$
(5.19)

uniformly on M. In view of Parseval’s equality,

$$\Biggl(\sum_{k=-N}^{N} \frac{\vert \langle f(\cdot),\varphi(\cdot,\lambda_{k}) \rangle _{\mathcal{H}}\vert ^{2}}{ \|\varphi(\cdot,\lambda_{k})\|_{\mathcal{H}}^{2}} \Biggr)^{1/2}\rightarrow 0 \quad \text{as } N\rightarrow \infty. $$

Thus \(\nu_{N}(\lambda)\rightarrow0\) uniformly on M. Hence (5.2) converges uniformly on M. Thus \(\mathcal{F}(\lambda)\) is an entire function. From the relation

$$\begin{aligned} \bigl\vert \mathcal{F}(\lambda)\bigr\vert \le& \frac{1}{r_{1}} \int_{a}^{c_{1}}\bigl\vert f_{1}(x)\bigr\vert \bigl\vert \varphi_{1,1}(x,\lambda )\bigr\vert \,dx + \frac{1}{r_{1}} \int_{a}^{c_{1}}\bigl\vert f_{2}(x)\bigr\vert \bigl\vert \varphi_{2,1}(x,\lambda )\bigr\vert \,dx \\ &{}+\sum_{k=2}^{m}\frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int _{c_{k-1}}^{c_{k}}\bigl\vert f_{1}(x)\bigr\vert \bigl\vert \varphi_{1,k}(x,\lambda)\bigr\vert \,dx \\ &{}+\sum _{k=2}^{m}\frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int _{c_{k-1}}^{c_{k}}\bigl\vert f_{2}(x)\bigr\vert \bigl\vert \varphi_{2,k}(x,\lambda)\bigr\vert \,dx +\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int _{c_{m}}^{b}\bigl\vert f_{1}(x)\bigr\vert \bigl\vert \varphi_{1,m+1}(x,\lambda)\bigr\vert \,dx \\ &{}+ \frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int _{c_{m}}^{b}\bigl\vert f_{2}(x)\bigr\vert \bigl\vert \varphi_{2,m+1}(x,\lambda)\bigr\vert \,dx,\quad \lambda\in\mathbb{C}, \end{aligned}$$

and the fact that \(\varphi_{1,k}(\cdot,\lambda)\) and \(\varphi _{2,k}(\cdot,\lambda)\), \(k=1,2,\ldots,m+1\), are entire functions of exponential type, see Lemma 7.2.1 in [2], we conclude that \(\mathcal{F}(\lambda)\) is of exponential type. □

Remark 5.2

If we replace \(\omega(\lambda)\) by the canonical product

$$ \widetilde{\omega}(\lambda)=(\lambda-\lambda_{0})\prod _{n=1}^{\infty} \biggl(1-\frac{\lambda}{\lambda_{n}} \biggr) \biggl(1-\frac{\lambda }{\lambda_{-n}} \biggr), $$
(5.20)

then expansion (5.2) is a Lagrange-type interpolation. Indeed, from Hadamard’s factorization theorem, see [10], there is an entire function with no zeros, \(h(\lambda)\), such that

$$\omega(\lambda)=h(\lambda)\widetilde{\omega}(\lambda). $$

Thus,

$$\frac{\omega(\lambda)}{\omega^{\prime}(\lambda_{n})} =\frac{h(\lambda)\widetilde{\omega}(\lambda)}{h(\lambda _{n})\widetilde{\omega}^{\prime}(\lambda_{n})} $$

and (5.1), (5.2) remain valid for the function \(\mathcal{F}(\lambda)/h(\lambda)\). Hence

$$ \mathcal{F}(\lambda)=\sum_{n=-\infty}^{\infty} \mathcal{F}(\lambda _{n})\frac{h(\lambda)\widetilde{\omega}(\lambda)}{ h(\lambda_{n})\widetilde{\omega}^{\prime}(\lambda_{n})(\lambda -\lambda_{n})}. $$
(5.21)

We may redefine (5.1) by taking the kernel \(\frac{\varphi(\cdot,\lambda)}{h(\lambda)}=\widetilde{\varphi }(\cdot,\lambda)\) to get

$$ \widetilde{\mathcal{F}}(\lambda)=\frac{\mathcal{F}(\lambda )}{h(\lambda)}=\sum _{n=-\infty}^{\infty}\widetilde{\mathcal{F}} ( \lambda_{n})\frac{\widetilde{\omega}(\lambda)}{ (\lambda-\lambda_{n})\widetilde{\omega}^{\prime}(\lambda_{n})}. $$
(5.22)

For more details as regards the argument of Remark 5.2, see [9], p.110.

The next theorem is devoted to vector-type interpolation sampling expansions associated with problem (2.1)-(2.4) for vector-valued integral transforms whose kernels are defined in terms of the Green’s matrix. As we see in (4.16), the Green’s matrix \(G(x,\xi,\lambda)\) of problem (2.1)-(2.4) has simple poles at \(\{\lambda_{k}\}_{k=-\infty}^{\infty}\). Define the function \(G(x,\lambda)\) to be \(G(x,\lambda):=\omega(\lambda )G(x,\xi_{0},\lambda)\), where \(\xi_{0}\in[a,c_{1})\cup(c_{1},c_{2})\cup(c_{2},c_{3})\cup\cdots \cup(c_{m},b]\) is a fixed point and \(\omega(\lambda)\) is the function defined in (2.24) or it is the canonical product (5.20).

Theorem 5.3

Let \(f(x)=\bigl ( {\scriptsize\begin{matrix}{} f_{1}(x) \cr f_{2}(x)\end{matrix}} \bigr ) \in\mathcal{H}\). Let \(\frak{F}(\lambda)=\bigl ( {\scriptsize\begin{matrix}{} \mathfrak{F}_{1}(\lambda)\cr \mathfrak{F}_{2}(\lambda)\end{matrix}} \bigr )\) be the vector-valued transform

$$\begin{aligned} \begin{aligned}[b] \frak{F}(\lambda)={}& \frac{1}{r_{1}} \int^{c_{1}}_{a}G(x,\lambda) f(x) \,dx+\sum _{k=2}^{m}\frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}G(x,\lambda) f(x) \,dx \\ &{}+ \frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int ^{b}_{c_{m}}G(x,\lambda) f(x) \,dx. \end{aligned} \end{aligned}$$
(5.23)

Then \(\frak{F}(\lambda)\) is a vector-valued entire function of exponential type that admits the vector-valued sampling expansion

$$ \frak{F}(\lambda)=\sum_{n= -\infty}^{\infty} \frak{F}(\lambda_{n})\frac{\omega(\lambda)}{ (\lambda-\lambda_{n})\omega^{\prime}(\lambda_{n})}. $$
(5.24)

The vector-valued series (5.24) converges absolutely on \(\mathbb{C}\) and uniformly on compact subsets of \(\mathbb{C}\). Here (5.24) means

$$ \begin{aligned} &\mathfrak{F}_{1}(\lambda)=\sum_{n=-\infty}^{\infty} \mathfrak {F}_{1}(\lambda_{n})\frac{\omega(\lambda)}{ (\lambda-\lambda_{n})\omega^{\prime}(\lambda_{n})}, \\ &\mathfrak{F}_{2}(\lambda)=\sum_{n=-\infty}^{\infty} \mathfrak {F}_{2}(\lambda_{n})\frac{\omega(\lambda)}{ (\lambda-\lambda_{n})\omega^{\prime}(\lambda_{n})}, \end{aligned} $$
(5.25)

where both series converge absolutely on \(\mathbb{C}\) and uniformly on compact sets of \(\mathbb{C}\).

Proof

Let \(\lambda\in\mathbb{C}\) such that \(\lambda\neq\lambda_{n}\) for \(n\in\mathbb{Z}\). Since \(\frak{F}(\lambda)\) is nothing but the unique solution (multiplied by \(\omega(\lambda)\)) of the inhomogeneous Dirac system (4.1), (2.2)-(2.4) when \(\xi=\xi_{0}\), see [1], p.77, \(\frak{F}(\lambda)\) has the expansion

$$\begin{aligned}& \frak{F}(\lambda)=\sum_{n=-\infty}^{\infty} \frac{a_{n}\omega(\lambda)}{\lambda_{n}-\lambda}\psi_{n}(\xi _{0}), \\& a_{n}= \frac{1}{r_{1}} \int^{c_{1}}_{a}{f}^{\top}(x) \psi_{n}(x) \,dx+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}{f}^{\top}(x) \psi_{n}(x) \,dx \\& \hphantom{a_{n}={}}{}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}{f}^{\top }(x) \psi_{n}(x) \,dx \end{aligned}$$
(5.26)

and \(\{\psi_{n}(\cdot) \}_{n=-\infty}^{\infty}\) is a complete orthonormal set of vector-valued eigenfunctions. Thus

$$ \frak{F}(\lambda)=\sum_{n=-\infty}^{\infty} \frac{a_{n}\omega(\lambda)}{\lambda_{n}-\lambda} \frac{\varphi(\xi_{0}, \lambda_{n})}{\|\varphi(\cdot,\lambda_{n})\|^{2}_{\mathcal{H}}}. $$
(5.27)

Now, we calculate \(\|\varphi(\cdot,\lambda_{n})\|^{2}_{\mathcal{H}}\) using the Green’s identity [1], p. 51. Let \(\lambda,\mu\in \mathbb{C}\) be different and \(\mu\neq\lambda_{n}\neq\lambda\); then, from (2.20) and (2.22),

$$\begin{aligned}& (\lambda-\mu) \Biggl(\frac{1}{r_{1}} \int^{c_{1}}_{a}\varphi^{\top}(x,\lambda ) \overline{\chi}(x,\mu) \,dx+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}\varphi^{\top}(x,\lambda) \overline{\chi}(x,\mu ) \,dx \\& \qquad {}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}\varphi ^{\top}(x,\lambda) \overline{\chi}(x,\mu) \,dx \Biggr) \\& \quad =\prod_{i=1}^{m}D_{i} \bigl[\varphi_{1,m+1}(b,\lambda)\overline{\chi }_{2,m+1}(b,\mu)- \varphi_{2,m+1}(b,\lambda)\overline{\chi }_{1,m+1}(b,\mu)\bigr] \\& \qquad {}- \bigl(\varphi_{1,1}(a,\lambda)\overline{\chi}_{2,1}(a,\mu)- \varphi _{2,1}(a,\lambda)\overline{\chi}_{1,1}(a,\mu)\bigr), \end{aligned}$$
(5.28)

where \(\varphi(x,\lambda)\) and \(\chi(x,\mu)\) are the solutions defined above. By the initial conditions (2.16) and (2.18), we have

$$\begin{aligned}& (\lambda-\mu) \Biggl(\frac{1}{r_{1}} \int^{c_{1}}_{a}\varphi^{\top}(x,\lambda ) \overline{\chi}(x,\mu) \,dx+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}\varphi^{\top}(x,\lambda) \overline{\chi}(x,\mu ) \,dx \\& \qquad {}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}\varphi ^{\top}(x,\lambda) \overline{\chi}(x,\mu) \,dx \Biggr) \\& \quad =-\prod_{i=1}^{m}D_{i} \bigl(\sin\beta\varphi_{1,m+1}(b,\lambda)+\cos \beta\varphi_{2,m+1}(b, \lambda)\bigr) \\& \qquad {}-\bigl(\cos\alpha\overline{\chi}_{2,1}(a,\mu)+\sin\alpha \overline {\chi}_{1,1}(a,\mu)\bigr). \end{aligned}$$
(5.29)

Letting \(\mu=\lambda_{n}\), for some n, and noting that \(\chi(x,\lambda_{n})\) is a real-valued eigenfunction, from (2.16) and (2.36), equation (5.29) becomes

$$\begin{aligned}& \tau_{n}(\lambda-\lambda_{n}) \Biggl( \frac{1}{r_{1}} \int^{c_{1}}_{a}\varphi^{\top}(x,\lambda ) \varphi(x,\lambda_{n}) \,dx+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}\varphi^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx \\& \qquad {}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}\varphi ^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx \Biggr) \\& \quad =-\prod_{i=1}^{m}D_{i} \bigl(\sin\beta\varphi_{1,m+1}(b,\lambda)+\cos \beta\varphi_{2,m+1}(b, \lambda)\bigr). \end{aligned}$$
(5.30)

Then from (2.24), we obtain

$$\begin{aligned}& \frac{1}{r_{1}} \int^{c_{1}}_{a}\varphi^{\top}(x,\lambda) \varphi (x,\lambda_{n}) \,dx+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}\varphi^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx \\& \quad {}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}\varphi ^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx=\frac{\tau_{n}^{-1}\omega(\lambda)}{\lambda-\lambda_{n}}. \end{aligned}$$
(5.31)

Taking the limit in (5.31) when \(\lambda\rightarrow\lambda_{n}\),

$$\begin{aligned} \bigl\Vert \varphi(\cdot,\lambda_{n})\bigr\Vert ^{2}_{\mathcal{H}} =&\frac {1}{r_{1}} \int^{c_{1}}_{a}\varphi^{\top}(x,\lambda) \varphi (x,\lambda_{n}) \,dx+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}\varphi^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx \\ &{}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int^{b}_{c_{m}}\varphi ^{\top}(x,\lambda) \varphi(x,\lambda_{n}) \,dx=\tau_{n}^{-1} \omega'(\lambda_{n}). \end{aligned}$$
(5.32)

Substituting in (5.27), we get the following series, which holds if λ is not an eigenvalue:

$$ \frak{F}(\lambda)=\sum_{n=-\infty}^{\infty} \frac{a_{n}\tau_{n}\omega(\lambda)}{\lambda_{n}-\lambda} \frac{\varphi(\xi_{0}, \lambda_{n})}{\omega'(\lambda_{n})}. $$
(5.33)

To prove (5.24), it suffices to prove that \(\frak{F}(\lambda_{n})=a_{n}\tau_{n}\varphi(\xi_{0}, \lambda_{n})\). Indeed, from (4.16) and (5.23), we have

$$ \frak{F}(\lambda_{n})=\left ( \textstyle\begin{array}{@{}c@{}} \mathfrak{F}_{1}(\lambda_{n}) \\ \mathfrak{F}_{2}(\lambda_{n}) \end{array}\displaystyle \right ), $$
(5.34)

where

$$ \mathfrak{F}_{1}(\lambda_{n})=\left \{ \textstyle\begin{array}{l} \frac{\varphi_{1}(\xi_{0},\lambda_{n})}{r_{1} } \int_{\xi_{0}}^{c_{1}}\chi^{\top}(x,\lambda _{n})f(x) \,dx+\frac{\chi_{1}(\xi_{0},\lambda_{n})}{r_{1} } \int_{a}^{\xi_{0}}\varphi^{\top}(x,\lambda_{n})f(x) \,dx \\ \quad {}+ \sum_{i=1}^{m} \frac{\varphi_{1}(\xi_{0},\lambda _{n})\prod_{j=1}^{i}D_{j}}{r_{i+1} } \int _{c_{i}}^{c_{i+1}}\chi^{\top}(x,\lambda_{n})f(x) \,dx, \quad \xi_{0}\in[a,c_{1}), \\ \frac{\varphi_{1}(\xi_{0},\lambda_{n})\prod_{i=1}^{k-1}D_{i}}{r_{k} } \int_{\xi_{0}}^{c_{k}}\chi^{\top}(x,\lambda _{n})f(x) \,dx \\ \quad {}+\frac{\chi_{1}(\xi_{0},\lambda_{n})\prod_{i=1}^{k-1}D_{i}}{r_{k} } \int_{c_{k-1}}^{\xi_{0}}\varphi^{\top}(x,\lambda _{n})f(x) \,dx \\ \quad {}+ \sum_{i=k}^{m} \frac{\varphi_{1}(\xi_{0},\lambda _{n})\prod_{j=1}^{i}D_{j}}{r_{i+1} } \int _{c_{i}}^{c_{i+1}}\chi^{\top}(x,\lambda_{n})f(x) \,dx \\ \quad {}+ \sum_{i=1}^{k-1}\frac{\chi_{1}(\xi_{0},\lambda _{n})\prod_{j=1}^{i}D_{j-1}}{r_{i} } \int _{c_{i-1}}^{c_{i}}\varphi^{\top}(x,\lambda_{n})f(x) \,dx, \\ \quad \xi_{0}\in(c_{k-1}, c_{k})\ (k=2,3,\ldots,m), D_{0}=1, \\ \frac{\varphi_{1}(\xi_{0},\lambda_{n})\prod_{i=1}^{m}D_{i}}{r_{m+1} } \int_{\xi_{0}}^{b}\chi^{\top}(x,\lambda_{n})f(x) \,dx \\ \quad {}+\frac{\chi_{1}(\xi_{0},\lambda_{n})\prod_{i=1}^{m}D_{i}}{r_{m+1} } \int_{c_{m}}^{\xi_{0}}\varphi^{\top}(x,\lambda _{n})f(x) \,dx \\ \quad {}+ \sum_{i=1}^{m}\frac{\chi_{1}(\xi_{0},\lambda _{n})\prod_{j=1}^{i}D_{j-1}}{r_{i} } \int _{c_{i-1}}^{c_{i}}\varphi^{\top}(x,\lambda_{n})f(x) \,dx,\quad \xi_{0}\in(c_{m}, b] \end{array}\displaystyle \right . $$
(5.35)

and

$$ \mathfrak{F}_{2}(\lambda_{n})=\left \{ \textstyle\begin{array}{l} \frac{\varphi_{2}(\xi_{0},\lambda_{n})}{r_{1} } \int_{\xi_{0}}^{c_{1}}\chi^{\top}(x,\lambda _{n})f(x) \,dx+\frac{\chi_{2}(\xi_{0},\lambda_{n})}{r_{1} } \int_{a}^{\xi_{0}}\varphi^{\top}(x,\lambda_{n})f(x) \,dx \\ \quad {}+ \sum_{i=1}^{m} \frac{\varphi_{2}(\xi_{0},\lambda _{n})\prod_{j=1}^{i}D_{j}}{r_{i+1} } \int _{c_{i}}^{c_{i+1}}\chi^{\top}(x,\lambda_{n})f(x) \,dx,\quad \xi_{0}\in[a,c_{1}), \\ \frac{\varphi_{2}(\xi_{0},\lambda_{n})\prod_{i=1}^{k-1}D_{i}}{r_{k} } \int_{\xi_{0}}^{c_{k}}\chi^{\top}(x,\lambda _{n})f(x) \,dx \\ \quad {}+\frac{\chi_{2}(\xi_{0},\lambda_{n})\prod_{i=1}^{k-1}D_{i}}{r_{k} } \int_{c_{k-1}}^{\xi_{0}}\varphi^{\top}(x,\lambda _{n})f(x) \,dx \\ \quad {}+ \sum_{i=k}^{m} \frac{\varphi_{2}(\xi_{0},\lambda _{n})\prod_{j=1}^{i}D_{j}}{r_{i+1} } \int _{c_{i}}^{c_{i+1}}\chi^{\top}(x,\lambda_{n})f(x) \,dx \\ \quad {}+ \sum_{i=1}^{k-1}\frac{\chi_{2}(\xi_{0},\lambda _{n})\prod_{j=1}^{i}D_{j-1}}{r_{i} } \int _{c_{i-1}}^{c_{i}}\varphi^{\top}(x,\lambda_{n})f(x) \,dx, \\ \quad \xi_{0}\in(c_{k-1}, c_{k})\ (k=2,3,\ldots,m), D_{0}=1, \\ \frac{\varphi_{2}(\xi_{0},\lambda_{n})\prod_{i=1}^{m}D_{i}}{r_{m+1} } \int_{\xi_{0}}^{b}\chi^{\top}(x,\lambda_{n})f(x) \,dx \\ \quad {}+\frac{\chi_{2}(\xi_{0},\lambda_{n})\prod_{i=1}^{m}D_{i}}{r_{m+1} } \int_{c_{m}}^{\xi_{0}}\varphi^{\top}(x,\lambda _{n})f(x) \,dx \\ \quad {}+ \sum_{i=1}^{m}\frac{\chi_{2}(\xi_{0},\lambda _{n})\prod_{j=1}^{i}D_{j-1}}{r_{i} } \int _{c_{i-1}}^{c_{i}}\varphi^{\top}(x,\lambda_{n})f(x) \,dx,\quad \xi_{0}\in(c_{m}, b]. \end{array}\displaystyle \right . $$
(5.36)

Now we prove that \(\mathfrak{F}_{1}(\lambda_{n})=a_{n}\tau_{n}\varphi_{1}(\xi_{0}, \lambda_{n})\) and \(\mathfrak{F}_{2}(\lambda_{n})=a_{n}\tau_{n}\varphi_{2}(\xi_{0}, \lambda_{n})\). From (2.36), (5.35), and (5.36), we get

$$\begin{aligned} \mathfrak{F}_{1}(\lambda_{n}) =&\tau_{n} \varphi_{1}(\xi_{0},\lambda_{n}) \Biggl( \frac{1}{r_{1}} \int^{c_{1}}_{a}\varphi^{\top }(x, \lambda_{n})f(x) \,dx \\ &{}+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}\varphi^{\top}(x, \lambda_{n})f(x) \,dx \\ &{}+\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int ^{b}_{c_{m}}\varphi^{\top}(x, \lambda_{n})f(x) \,dx\Biggr) \\ =&\tau_{n}\varphi_{1}(\xi_{0}, \lambda_{n}) \Biggl(\frac{1}{r_{1}} \int^{c_{1}}_{a}f^{\top }(x)\varphi(x, \lambda_{n}) \,dx \\ &{}+\sum_{k=2}^{m} \frac{\prod_{i=1}^{k-1}D_{i}}{r_{k}} \int ^{c_{k}}_{c_{k-1}}f^{\top}(x)\varphi(x, \lambda_{n}) \,dx \\ &{} +\frac{\prod_{i=1}^{m}D_{i}}{r_{m+1}} \int ^{b}_{c_{m}}f^{\top}(x)\varphi(x, \lambda_{n}) \,dx\Biggr) \\ =&a_{n}\tau_{n}\varphi_{1}(\xi_{0}, \lambda_{n}), \end{aligned}$$

for \(x,\xi_{0}\in[a,c_{1})\cup(c_{1},c_{2})\cup(c_{2},c_{3})\cup\cdots \cup(c_{m},b]\). The same for \(\mathfrak{F}_{2}(\lambda_{n})\) and the proof of (5.24) is complete. The convergence properties as well as the analytic and growth properties can be established as in Theorem 5.1 above. □

6 Examples

Example 6.1

The boundary value problem

$$\begin{aligned}& u_{2}^{\prime}-q(x)u_{1}=\lambda u_{1},\qquad u_{1}^{\prime}+q(x)u_{2}=- \lambda u_{2},\quad x\in[-1,0)\cup\biggl(0,\frac{1}{2} \biggr)\cup\biggl(\frac{1}{2},1 \biggr], \end{aligned}$$
(6.1)
$$\begin{aligned}& u_{1}(-1)=u_{1}(1)=0, \end{aligned}$$
(6.2)
$$\begin{aligned}& u_{1}\bigl(0^{-}\bigr)-2 u_{1} \bigl(0^{+}\bigr)=0,\qquad u_{2}\bigl(0^{-}\bigr)-2 u_{2}\bigl(0^{+}\bigr)=0, \end{aligned}$$
(6.3)
$$\begin{aligned}& 2u_{1}\biggl(\frac{1}{2}^{-}\biggr)- u_{1}\biggl(\frac{1}{2}^{+}\biggr)=0,\qquad 2u_{2}\biggl(\frac{1}{2}^{-}\biggr)- u_{2} \biggl(\frac{1}{2}^{+}\biggr)=0, \end{aligned}$$
(6.4)

is a special case of problem (2.1)-(2.4) when \(m=2\), \(r_{1}=r_{2}=r_{3}=1\), \(\gamma_{1}=\gamma_{1}^{\prime}=\delta_{2}=\delta_{2}^{\prime }=1\), \(\gamma_{2}=\gamma_{2}^{\prime}=\delta_{1}=\delta _{1}^{\prime}=2\), \(\alpha=\beta=\frac{\pi}{2}\) and \(q_{1}(x)=q_{2}(x)=q(x)\),

$$q(x)=\left \{ \textstyle\begin{array}{l@{\quad}l} -1, & {-}1\le x< 0 , \\ 0, & 0< x< \frac{1}{2}, \\ 1, & \frac{1}{2}< x\le1 . \end{array}\displaystyle \right . $$

Then \(\gamma_{k}\delta_{k}^{\prime}=\gamma_{k}^{\prime}\delta _{k}\), \(k=1,2\). In the notations of Section 2, the solutions \(\varphi(\cdot,\lambda)\) and \(\chi(\cdot,\lambda)\) are

$$\begin{aligned}& \varphi(x,\lambda)=\left \{ \textstyle\begin{array}{l@{\quad}l} \bigl[ {\scriptsize\begin{matrix}{} \varphi_{1,1}(x,\lambda) \cr \varphi_{2,1}(x,\lambda) \end{matrix}} \bigr]= \bigl[ {\scriptsize\begin{matrix}{} \sin[(\lambda-1)(x+1)] \cr -\cos[(\lambda-1)(x+1)] \end{matrix}} \bigr], & x\in [-1,0), \\ \bigl[ {\scriptsize\begin{matrix}{} \varphi_{1,2}(x,\lambda) \cr \varphi_{2,2}(x,\lambda) \end{matrix}} \bigr]= \bigl[ {\scriptsize\begin{matrix}{} -\frac{1}{2}\sin[1-\lambda(x+1)] \cr -\frac{1}{2}\cos[1-\lambda(x+1)] \end{matrix}} \bigr], & x\in(0,\frac{1}{2}), \\ \bigl[ {\scriptsize\begin{matrix}{} \varphi_{1,3}(x,\lambda) \cr \varphi_{2,3}(x,\lambda) \end{matrix}} \bigr]= \bigl[ {\scriptsize\begin{matrix}{} -\sin [\frac{3}{2}-\lambda-(\lambda+1)x ] \cr -\cos [\frac{3}{2}-\lambda-(\lambda+1)x ] \end{matrix}} \bigr], & x\in(\frac{1}{2},1 ], \end{array}\displaystyle \right . \end{aligned}$$
(6.5)
$$\begin{aligned}& \chi(x,\lambda)=\left \{ \textstyle\begin{array}{l@{\quad}l} \bigl[ {\scriptsize\begin{matrix}{} \chi_{1,1}(x,\lambda) \cr \chi_{2,1}(x,\lambda) \end{matrix}} \bigr]= \bigl[ {\scriptsize\begin{matrix}{} -\sin[\lambda+\frac{1}{2}+(1-\lambda)x] \cr -\cos[\lambda+\frac{1}{2}+(1-\lambda)x] \end{matrix}} \bigr], & x\in[-1,0), \\ \bigl[ {\scriptsize\begin{matrix}{} \chi_{1,2}(x,\lambda) \cr \chi_{2,2}(x,\lambda) \end{matrix}} \bigr]=\bigl[ {\scriptsize\begin{matrix}{} -\frac{1}{2}\sin[\lambda+\frac{1}{2}-\lambda x] \cr -\frac{1}{2}\cos[\lambda+\frac{1}{2}-\lambda x] \end{matrix}} \bigr], & x\in(0,\frac{1}{2} ), \\ \bigl[ {\scriptsize\begin{matrix}{} \chi_{1,3}(x,\lambda) \cr \chi_{2,3}(x,\lambda) \end{matrix}} \bigr]=\bigl[ {\scriptsize\begin{matrix}{} \sin[(\lambda+1)(x-1)] \cr -\cos[(\lambda+1)(x-1)] \end{matrix}} \bigr], & x\in(\frac{1}{2},1]. \end{array}\displaystyle \right . \end{aligned}$$
(6.6)

The eigenvalues are the solutions of the equation

$$ \omega(\lambda)=\sin\biggl[\frac{1}{2}-2\lambda\biggr]=0. $$
(6.7)

Hence the eigenvalues are

$$ \lambda_{n}=\frac{1-2n\pi}{4},\quad n\in\mathbb{Z}, $$
(6.8)

which is illustrated in Figure 1. The Green’s matrix of problem (6.1)-(6.4) is given by

$$ G(x,\xi,\lambda)=\frac{1}{\sin[\frac{1}{2}-2\lambda]} \left \{\textstyle\begin{array}{l@{\quad}l} G_{1}(x,\xi,\lambda),& {-}1\le\xi\le x< 0, \\ G_{2}(x,\xi,\lambda),& {-}1\le x\le\xi< 0, \\ G_{3}(x,\xi,\lambda),& {-}1\le\xi< 0, 0< x< \frac{1}{2}, \\ G_{4}(x,\xi,\lambda),& {-}1\le x< 0, 0< \xi< \frac{1}{2}, \\ G_{5}(x,\xi,\lambda),& {-}1\le\xi< 0, \frac{1}{2}< x\le1, \\ G_{6}(x,\xi,\lambda),& {-}1\le x< 0, \frac{1}{2}< \xi\le1, \\ G_{7}(x,\xi,\lambda) , &0< \xi\le x< \frac {1}{2}, \\ G_{8}(x,\xi,\lambda), & 0< x\le\xi< \frac {1}{2}, \\ G_{9}(x,\xi,\lambda),& 0< \xi< \frac{1}{2}, \frac{1}{2}< x\le1, \\ G_{10}(x,\xi,\lambda),& 0< x< \frac{1}{2}, \frac{1}{2}< \xi\le1, \\ G_{11}(x,\xi,\lambda) ,& \frac{1}{2}< \xi\le x\le 1, \\ G_{12}(x,\xi,\lambda),& \frac{1}{2}< x\le\xi \le1, \end{array}\displaystyle \right . $$
(6.9)

where

$$\begin{aligned}& G_{1}(x,\xi,\lambda) \\& \quad = \left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} -\sin[(\lambda-1)(\xi+1)]\sin[\lambda+\frac{1}{2}+(1-\lambda)x] & \cos[(\lambda-1)(\xi+1)] \sin[\lambda+\frac{1}{2}+(1-\lambda)x] \\ -\sin[(\lambda-1)(\xi+1)]\cos[\lambda+\frac{1}{2}+(1-\lambda)x] & \cos[(\lambda-1)(\xi+1)] \cos[\lambda+\frac{1}{2}+(1-\lambda)x] \end{array}\displaystyle \right ), \\& G_{2}(x,\xi,\lambda) \\& \quad = \left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} -\sin[(\lambda-1)(x+1)]\sin[\lambda+\frac{1}{2}+(1-\lambda)\xi] & -\sin[(\lambda-1)(x+1)] \cos[\lambda+\frac{1}{2}+(1-\lambda)\xi] \\ \cos[(\lambda-1)(x+1)]\sin[\lambda+\frac{1}{2}+(1-\lambda)\xi] & \cos[(\lambda-1)(x+1)] \cos[\lambda+\frac{1}{2}+(1-\lambda)\xi] \end{array}\displaystyle \right ), \\& G_{3}(x,\xi,\lambda) \\& \quad =\frac{1}{2}\left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} -\sin[(\lambda-1)(\xi+1)]\sin[\lambda+\frac{1}{2}-\lambda x] & \cos[(\lambda-1)(\xi+1)] \sin[\lambda+\frac{1}{2}-\lambda x] \\ -\sin[(\lambda-1)(\xi+1)]\cos[\lambda+\frac{1}{2}-\lambda x] & \cos[(\lambda-1)(\xi+1)] \cos[\lambda+\frac{1}{2}-\lambda x] \end{array}\displaystyle \right ), \\& G_{4}(x,\xi,\lambda) \\& \quad =\frac{1}{2}\left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} -\sin[(\lambda-1)(x+1)]\sin[\lambda+\frac{1}{2}-\lambda\xi] & -\sin[(\lambda-1)(x+1)] \cos[\lambda+\frac{1}{2}-\lambda\xi] \\ \cos[(\lambda-1)(x+1)]\sin[\lambda+\frac{1}{2}-\lambda\xi] & \cos[(\lambda-1)(x+1)] \cos[\lambda+\frac{1}{2}-\lambda\xi] \end{array}\displaystyle \right ), \\& G_{5}(x,\xi,\lambda) \\& \quad =\left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} \sin[(\lambda-1)(\xi+1)]\sin[(\lambda+1)(x-1)] & -\cos[(\lambda -1)(\xi+1)] \sin[(\lambda+1)(x-1)] \\ -\sin[(\lambda-1)(\xi+1)]\cos[(\lambda+1)(x-1)] & \cos[(\lambda -1)(\xi+1)] \cos[(\lambda+1)(x-1)] \end{array}\displaystyle \right ), \\& G_{6}(x,\xi,\lambda) \\& \quad =\left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} \sin[(\lambda-1)(x+1)]\sin[(\lambda+1)(\xi-1)] & -\sin[(\lambda-1)(x+1)] \cos[(\lambda+1)(\xi-1)] \\ -\cos[(\lambda-1)(x+1)]\sin[(\lambda+1)(\xi-1)] & \cos[(\lambda-1)(x+1)] \cos[(\lambda+1)(\xi-1)] \end{array}\displaystyle \right ), \\& G_{7}(x,\xi,\lambda) \\& \quad =\frac{1}{4}\left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} \sin[1-\lambda(\xi+1)]\sin[\lambda+\frac{1}{2}-\lambda x] & \cos [1-\lambda(\xi+1)] \sin[\lambda+\frac{1}{2}-\lambda x] \\ \sin[1-\lambda(\xi+1)]\cos[\lambda+\frac{1}{2}-\lambda x] & \cos [1-\lambda(\xi+1)] \cos[\lambda+\frac{1}{2}-\lambda x] \end{array}\displaystyle \right ), \\& G_{8}(x,\xi,\lambda) \\& \quad =\frac{1}{4}\left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} \sin[1-\lambda(x+1)]\sin[\lambda+\frac{1}{2}-\lambda\xi] & \sin [1-\lambda(x+1)] \cos[\lambda+\frac{1}{2}-\lambda\xi] \\ \cos[1-\lambda(x+1)]\sin[\lambda+\frac{1}{2}-\lambda\xi] & \cos [1-\lambda(x+1)] \cos[\lambda+\frac{1}{2}-\lambda\xi] \end{array}\displaystyle \right ), \\& G_{9}(x,\xi,\lambda) \\& \quad =\frac{1}{2}\left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} -\sin[1-\lambda(\xi+1)]\sin[(\lambda+1)(x-1)] & -\cos[1-\lambda (\xi+1)] \sin[(\lambda+1)(x-1)] \\ \sin[1-\lambda(\xi+1)]\cos[(\lambda+1)(x-1)] & \cos[1-\lambda (\xi+1)] \cos[(\lambda+1)(x-1)] \end{array}\displaystyle \right ), \\& G_{10}(x,\xi,\lambda) \\& \quad =\frac{1}{2}\left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} -\sin[1-\lambda(x+1)]\sin[(\lambda+1)(\xi-1)] & \sin[1-\lambda(x+1)] \cos[(\lambda+1)(\xi-1)] \\ -\cos[1-\lambda(x+1)]\sin[(\lambda+1)(\xi-1)] & \cos[1-\lambda(x+1)] \cos[(\lambda+1)(\xi-1)] \end{array}\displaystyle \right ), \\& G_{11}(x,\xi,\lambda) \\& \quad = \left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} -\sin[\frac{3}{2}-\lambda-(\lambda+1)\xi]\sin[(\lambda+1)(x-1)] & -\cos[\frac{3}{2}-\lambda-(\lambda+1)\xi] \sin[(\lambda+1)(x-1)] \\ \sin[\frac{3}{2}-\lambda-(\lambda+1)\xi]\cos[(\lambda+1)(x-1)] & \cos[\frac{3}{2}-\lambda-(\lambda+1)\xi] \cos[(\lambda+1)(x-1)] \end{array}\displaystyle \right ), \\& G_{12}(x,\xi,\lambda) \\& \quad = \left ( \textstyle\begin{array}{@{}c@{\quad}c@{}} -\sin[\frac{3}{2}-\lambda-(\lambda+1)x]\sin[(\lambda+1)(\xi-1)] & \sin[\frac{3}{2}-\lambda-(\lambda+1)x] \cos[(\lambda+1)(\xi-1)] \\ -\cos[\frac{3}{2}-\lambda-(\lambda+1)x]\sin[(\lambda+1)(\xi-1)] & \cos[\frac{3}{2}-\lambda-(\lambda+1)x] \cos[(\lambda+1)(\xi-1)] \end{array}\displaystyle \right ). \end{aligned}$$

By Theorem 5.1, the transform

$$\begin{aligned} \mathcal{F}(\lambda) =& \int_{-1}^{0} \bigl[f_{1}(x)\sin\bigl[( \lambda-1) (x+1)\bigr] -f_{2}(x)\cos\bigl[(\lambda-1) (x+1)\bigr] \bigr] \,dx \\ &{}-2 \int_{0}^{\frac{1}{2}} \bigl[f_{1}(x)\sin\bigl[1- \lambda(x+1)\bigr] +f_{2}(x)\cos\bigl[1-\lambda(x+1)\bigr] \bigr] \,dx \\ &{}- \int_{\frac{1}{2}}^{1} \biggl[f_{1}(x)\sin\biggl[ \frac{3}{2}-\lambda -(\lambda+1)x\biggr] +f_{2}(x)\cos\biggl[ \frac{3}{2}-\lambda-(\lambda+1)x\biggr] \biggr] \,dx \end{aligned}$$
(6.10)

has the following expansion:

$$ \mathcal{F}(\lambda)=\sum_{n=-\infty}^{\infty} \mathcal{F} \biggl(\frac{1-2n\pi}{4} \biggr) \frac{\sin [\frac{1}{2}+n\pi-2\lambda ]}{ \frac{1}{2}-n\pi-2\lambda}, $$
(6.11)

where \(\omega' (\frac{1-2n\pi}{4} )=2(-1)^{n+1}\).

Figure 1
figure 1

The eigenvalues of \(\pmb{\omega(\lambda)}\) .

In view of Theorem 5.3, the vector-valued transform

$$ \mathfrak{F}(\lambda)= \left \{ \textstyle\begin{array}{l@{\quad}l} \bigl[ {\scriptsize\begin{matrix}{} -\sin[\lambda+\frac{1}{2}+(1-\lambda)\xi_{0}] \cr -\cos[\lambda+\frac{1}{2}+(1-\lambda)\xi_{0}] \end{matrix}} \bigr]\mathcal{T}_{1}+\bigl[ {\scriptsize\begin{matrix}{} \sin[(\lambda-1)(\xi_{0}+1)] \cr -\cos[(\lambda-1)(\xi_{0}+1)] \end{matrix}} \bigr](\mathcal{T}_{2}+\mathcal{T}_{3}+ \mathcal{T}_{4}),& {-}1\le \xi_{0}< 0, \\ \frac{1}{2}\bigl[ {\scriptsize\begin{matrix}{} -\sin[\lambda+\frac{1}{2}-\lambda\xi_{0}] \cr -\cos[\lambda+\frac{1}{2}-\lambda\xi_{0}] \end{matrix}} \bigr]( \mathcal{T}_{5}+\mathcal{T}_{6})+\frac{1}{2} \bigl[ {\scriptsize\begin{matrix}{} -\sin[1-\lambda(\xi_{0}+1)] \cr -\cos[1-\lambda(\xi_{0}+1)] \end{matrix}} \bigr](\mathcal{T}_{7}+ \mathcal{T}_{4}),& 0< \xi_{0}< \frac{1}{2}, \\ \bigl[ {\scriptsize\begin{matrix}{} \sin[(\lambda+1)(x-1)] \cr -\cos[(\lambda+1)(x-1)] \end{matrix}} \bigr](\mathcal{T}_{5}+ \mathcal{T}_{8}+\mathcal{T}_{9})+\bigl[ {\scriptsize\begin{matrix}{} -\sin[\frac{3}{2}-\lambda-(\lambda+1)x] \cr -\cos[\frac{3}{2}-\lambda-(\lambda+1)x] \end{matrix}} \bigr]\mathcal{T}_{10},& 0< \xi_{0}\le1, \end{array}\displaystyle \right . $$
(6.12)

where

$$\begin{aligned}& \mathcal{T}_{1}= \int_{-1}^{\xi_{0}} \bigl[\sin\bigl[(\lambda-1) (x+1) \bigr]f_{1}(x) -\cos\bigl[(\lambda-1) (x+1)\bigr]f_{2}(x) \bigr] \,dx, \\& \mathcal{T}_{2}=- \int_{\xi_{0}}^{0} \biggl[\sin\biggl[\lambda+ \frac {1}{2}+(1-\lambda)x\biggr]f_{1}(x) +\cos\biggl[\lambda+ \frac{1}{2}+(1-\lambda)x\biggr]f_{2}(x) \biggr] \,dx, \\& \mathcal{T}_{3}=-\frac{1}{2} \int_{0}^{\frac{1}{2}} \biggl[\sin \biggl[\lambda+ \frac{1}{2}-\lambda x\biggr]f_{1}(x) +\cos\biggl[\lambda+ \frac{1}{2}-\lambda x\biggr]f_{2}(x) \biggr] \,dx, \\& \mathcal{T}_{4}= \int_{\frac{1}{2}}^{1} \bigl[\sin\bigl[(\lambda +1) (x-1) \bigr]f_{1}(x) -\cos\bigl[(\lambda+1) (x-1)\bigr]f_{2}(x) \bigr] \,dx, \\& \mathcal{T}_{5}= \int_{-1}^{0} \bigl[\sin\bigl[(\lambda-1) (x+1) \bigr]f_{1}(x) -\cos\bigl[(\lambda-1) (x+1)\bigr]f_{2}(x) \bigr] \,dx, \\& \mathcal{T}_{6}=-\frac{1}{2} \int_{0}^{\xi_{0}} \bigl[\sin \bigl[1-\lambda(x+1) \bigr]f_{1}(x) +\cos\bigl[1-\lambda(x+1)\bigr]f_{2}(x) \bigr] \,dx, \\& \mathcal{T}_{7}=- \int_{\xi_{0}}^{\frac{1}{2}} \biggl[\sin\biggl[\lambda + \frac{1}{2}+(1-\lambda)x\biggr]f_{1}(x) +\cos\biggl[\lambda+ \frac{1}{2}+(1-\lambda)x\biggr]f_{2}(x) \biggr] \,dx, \\& \mathcal{T}_{8}=-\frac{1}{2} \int_{0}^{\frac{1}{2}} \bigl[\sin \bigl[1-\lambda(x+1) \bigr]f_{1}(x) +\cos\bigl[1-\lambda(x+1)\bigr]f_{2}(x) \bigr] \,dx, \\& \mathcal{T}_{9}=- \int_{\frac{1}{2}}^{\xi_{0}} \biggl[\sin\biggl[\frac {3}{2}- \lambda-(\lambda+1)x\biggr]f_{1}(x) +\cos\biggl[\frac{3}{2}- \lambda-(\lambda+1)x\biggr]f_{2}(x) \biggr] \,dx, \\& \mathcal{T}_{10}= \int_{\xi_{0}}^{1} \bigl[\sin\bigl[(\lambda+1) (x-1) \bigr]f_{1}(x) -\cos\bigl[(\lambda+1) (x-1)\bigr]f_{2}(x) \bigr] \,dx. \end{aligned}$$

The vector-valued transform (6.12) has the following vector-valued expansion:

$$ \mathfrak{F}(\lambda)=\left [ \textstyle\begin{array}{@{}c@{}} \mathfrak{F}_{1}(\lambda)=\sum_{n=-\infty}^{\infty }\mathfrak{F}_{1} (\frac{1-2n\pi}{4} ) \frac{\sin [\frac{1}{2}+n\pi-2\lambda ]}{\frac{1}{2}-n\pi-2\lambda} \\ \mathfrak{F}_{2}(\lambda)=\sum_{n=-\infty}^{\infty }\mathfrak{F}_{2} (\frac{1-2n\pi}{4} ) \frac{\sin [\frac{1}{2}+n\pi-2\lambda ]}{\frac{1}{2}-n\pi-2\lambda} \end{array}\displaystyle \right ]. $$
(6.13)

References

  1. Levitan, BM, Sargsjan, IS: Introduction to Spectral Theory: Selfadjoint Ordinary Differential Operators. Translations of Mathematical Monographs, vol. 39. Am. Math. Soc., Providence (1975)

    MATH  Google Scholar 

  2. Levitan, BM, Sargsjan, IS: Sturm-Liouville and Dirac Operators. Kluwer Academic, Dordrecht (1991)

    Book  Google Scholar 

  3. Fulton, CT: Two point boundary value problems with eigenvalue parameter contained in the boundary conditions. Proc. R. Soc. Edinb., Sect. A 77, 293-308 (1977)

    Article  MATH  MathSciNet  Google Scholar 

  4. Tharwat, MM, Yildirim, A, Bhrawy, AH: Sampling of discontinuous Dirac systems. Numer. Funct. Anal. Optim. 34(3), 323-348 (2013)

    Article  MATH  MathSciNet  Google Scholar 

  5. Higgins, JR: Sampling Theory in Fourier and Signal Analysis: Foundations. Oxford University Press, Oxford (1996)

    MATH  Google Scholar 

  6. Kotel’nikov, V: On the carrying capacity of the ‘ether’ and wire in telecommunications. In: Material for the First All Union Conference on Questions of Communications, vol. 55, pp. 55-64. Izd. Red. Upr. Svyazi RKKA, Moscow (1933)

    Google Scholar 

  7. Shannon, C: Communication in the presence of noise. Proc. IRE 37, 10-21 (1949)

    Article  MathSciNet  Google Scholar 

  8. Whittaker, E: On the functions which are represented by the expansion of the interpolation theory. Proc. R. Soc. Edinb., Sect. A 35, 181-194 (1915)

    Article  Google Scholar 

  9. Zayed, AI: Advances in Shannon’s Sampling Theory. CRC Press, Boca Raton (1993)

    MATH  Google Scholar 

  10. Boas, R: Entire Functions. Academic Press, New York (1954)

    MATH  Google Scholar 

  11. Paley, R, Wiener, N: Fourier Transforms in the Complex Domain. Am. Math. Soc. Colloquium Publ. Ser., vol. 19. Am. Math. Soc., Providence (1934)

    MATH  Google Scholar 

  12. Lasser, R: Introduction to Fourier Series. Dekker, New York (1996)

    MATH  Google Scholar 

  13. Butzer, PL: A survey of the Whittaker-Shannon sampling theorem and some of its extensions. J. Math. Res. Exposition 3, 185-212 (1983)

    MathSciNet  Google Scholar 

  14. Butzer, P, Higgins, JR, Stens, RL: Sampling theory in signal analysis. In: Pier, JP (ed.) Development of Mathematics 1950-2000, pp. 193-234. Birkhäuser, Basel (2000)

    Chapter  Google Scholar 

  15. Butzer, P, Nasri-Roudsari, G: Kramer’s sampling theorem in signal analysis and its role in mathematics. In: Blackledge, JM (ed.) Image Processing: Mathematical Methods and Applications. Proc. IMA Conference, Cranfield University, UK, pp. 49-95. Clarendon, Oxford (1997)

    Google Scholar 

  16. Kramer, HP: A generalized sampling theorem. J. Math. Phys. 38, 68-72 (1959)

    Article  MATH  Google Scholar 

  17. Weiss, P: Sampling theorems associated with Sturm-Liouville systems. Bull. Am. Math. Soc. 163, 242 (1957)

    Google Scholar 

  18. Everitt, WN, Garcia, AG, Hernández-Medina, MA: On Lagrange-type interpolation series and analytic Kramer kernels. Results Math. 51, 215-228 (2008)

    Article  MATH  MathSciNet  Google Scholar 

  19. Garcia, AG, Littlejohn, LL: On analytic sampling theory. J. Comput. Appl. Math. 171, 235-246 (2004)

    Article  MATH  MathSciNet  Google Scholar 

  20. Higgins, JR: A sampling principle associated with Saitoh’s fundamental theory of linear transformations. In: Saitoh, S, Hayashi, N, Yamamoto, M (eds.) Analytic Extension Formulas and Their Applications. Kluwer Academic, Dordrecht (2001)

    Google Scholar 

  21. Abd-alla, MZ, Annaby, MH: Sampling of vector-valued transforms associated with Green’s matrix of Dirac systems. J. Math. Anal. Appl. 284, 104-117 (2003)

    Article  MATH  MathSciNet  Google Scholar 

  22. Zayed, A, Garcia, AG: Sampling theorems associated with a Dirac operator and the Hartley transform. J. Math. Anal. Appl. 214, 587-598 (1997)

    Article  MATH  MathSciNet  Google Scholar 

  23. Annaby, MH, Tharwat, MM: On sampling and Dirac systems with eigenparameter in the boundary conditions. J. Appl. Math. Comput. 36, 291-317 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  24. Tharwat, MM: Discontinuous Sturm-Liouville problems and associated sampling theories. Abstr. Appl. Anal. (2011). doi:10.1155/2011/610232

    MathSciNet  Google Scholar 

  25. Tharwat, MM: On sampling theories and discontinuous Dirac systems with eigenparameter in the boundary conditions. Bound. Value Probl. (2013). doi:10.1186/1687-2770-2013-65

    Google Scholar 

  26. Hira, F, Altinisik, N: Dirac system with discontinuities at two points. Abstr. Appl. Anal. 2014, Article ID 927693 (2014)

    Article  MathSciNet  Google Scholar 

  27. Hira, F, Altinisik, N: Expansion theorems for discontinuous Sturm-Liouville problems. J. Math. Comput. Sci. 4, 335-349 (2014)

    Google Scholar 

  28. Coddington, EA, Levinson, N: Theory of Ordinary Differential Equations. McGraw-Hill, New York (1955)

    MATH  Google Scholar 

Download references

Acknowledgements

The author is very grateful to the referees for carefully reading the paper and for their comments and suggestions, which have improved the paper.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mohammed M Tharwat.

Additional information

Competing interests

The author declares that he has no competing interests.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tharwat, M.M. Sampling theories of boundary value problems with several internal points of discontinuity. Bound Value Probl 2016, 4 (2016). https://doi.org/10.1186/s13661-015-0515-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-015-0515-1

MSC

Keywords