Acessibilidade / Reportar erro

Uniqueness Properties of The Solution of The Inverse Problem for The Sturm-Liouville Equation With Discontinuous Leading Coefficient

Abstract

The present paper studies uniqueness properties of the solution of the inverse problem for the Sturm-Liouville equation with discontinuous leading coefficient and the separated boundary conditions. It is proved that the considered boundary-value is uniquely reconstructed, i.e. the potential function of the equation and the constants in the boundary conditions are uniquely determined by given Weyl function or by the given spectral data.

Key words
Asymptotic formulas for eigenvalues; boundary value problems; inverse problems; spectral; analysis of ordinary differential operators; Sturm-Liouville theory; transformation operator

INTRODUCTION

This paper is concerned with the uniqueness theorems for the solution of some inverse spectral problems for the boundary value problem

- y ′′ + q ( x ) y = λ 2 ρ ( x ) y , 0 x π (1)
y ( 0 ) - h y ( 0 ) = 0 , y ( π ) + h 1 y ( π ) = 0 , (2)

where q(x) is real-valued function in L2(0,π), λ is a complex parameter, h ,h1 are real numbers,

ρ ( x ) = { 1 , 0 x a α 2 , a < x < π

with a(0,π), α1.

Inverse spectral problems consist in recovering differential operators from their spectral characteristics (see Marchenko 201120 MARCHENKO VA. 2011. Sturm-Liouville operators and applications, AMS Chelsea Publishing Volume 373, 393 p., Levitan 198714 LEVITAN BM. 1987. Inverse Surm-Liouville problems, VNU Sci. Press, Utrecht, 240 p.). Such problems arise in many areas of science and engineering (see Hald 19809 HALD O. 1980. Inverse eigenvalue problems for the mantle. Geophys J R Astron Soc 62: 41-48., Krueger 198213 KRUEGER R. 1982. Inverse problems for nonabsorbing media with discontinuous material properties. J Math Phys 23: 396-404. , Willis 198419 WILLIS C. 1984. Inverse problems for torsional modes. Geophys J R Astron Soc 78: 847-853.). The goal of this work is to prove the uniqueness theorems for the solution of the inverse problem which determines the potential function q(x) and the constants h, h1 by the Weyl function or by the spectral data of the boundary value problem (1)-(2).

In the classical case ρ(x)=1, the direct and inverse problems for the Sturm-Liouville operators have been completely studied (see Marchenko 2011, Levitan 1987, Levitan and Gasymov 196415 LEVITAN BM AND GASYMOV MG. 1964. Determination of a differential equation by two spectra. Uspehi Mat Nauk 19: 3-63., Freiling and Yurko 20016 FREILING G AND YURKO V. 2001. Inverse Sturm-Liouville problems and their applications, Nova Science Publishers, Inc., Huntington: NY, 305 p., Hryniv and Mykytyuk 200311 HRYNIV RO AND MYKYTYUK YAV. 2003. Inverse spectral problems for Sturm-Liouville operators with singular potentials. Inverse Probl 19: 665-684. , 200412 HRYNIV RO AND MYKYTYUK YAV. 2004. Half-inverse spectral problems for Sturm-Liouville operators with singular potentials. Inverse Probl 20: 1423-1444. , and the references therein). Direct and inverse problems for the discontinuous Sturm-Liouville boundary-value problems in different settings have been studied in Hald (198410 HALD O. 1984. Discontinuous inverse eigenvalue problems. Commun Pure Appl Math 37: 539-577. ), Andersson (19884 ANDERSSON L. 1988. Inverse eigenvalue problems for a Sturm-Liouville equation in impedance form. Inverse Probl 4: 929-971.), Guseinov and Pashaev (20028 GUSEINOV IM AND PASHAEV RT. 2002. On an inverse problem for a second-order differential equation. Russian Math Surveys 57: 597-598.), Carlson (19945 CARLSON R. 1994. An inverse spectral problem for Sturm-Liouville operators with discontinuous coefficients. Proc Amer Math Soc 120: 475-484.), Yurko (200021 YURKO VA. 2000. On boundary value problems with discontinuity conditions inside an interval. Differ Equ 36: 1266-1269.), Gasymov (19777 GASYMOV MG. 1977. The direct and inverse problem of spectral analysis for a class of equations with a discontinuous coefficient. Non-Classical Methods in Geophysics, 37-44, Nauka, Novosibirsk, Russia.), Amirov (20063 AMIROV RKH. 2006. On Sturm-Liouville operators with discontinuity conditions inside an interval. J Math Analysis Appl 317: 163-176.), Mamedov (200616 MAMEDOV KHR. 2006. Uniqueness of the solution of the inverse problem of scattering theory for Sturm-Liouville operator with discontinuous coefficient. Proc Inst Math Mech Natl Acad Sci Azerb 24: 163-172. , 201017 MAMEDOV KHR. 2010. On an inverse scattering problem for a discontinuous Sturm-Liouville equation with a spectral parameter in the boundary condition. Bound Value Probl ID 171967, 17 p.), Mamedov and Palamut (200918 MAMEDOV KHR AND PALAMUT N. 2009. On a direct problem of scattering theory for a class of Sturm-Liouville operator with discontinuous coefficient. Proc Jangjeon Math Soc 12: 243-251. )) and other works. Note that, the direct and inverse spectral problem for the equation(1) with simple boundary conditions on the interval (0,π) recently has been investigated in Akhmedova and Huseynov (201020 AKHMEDOVA EN AND HUSEYNOV HM. 2010. On inverse problem for Sturm-Liouville operator with discontinuous coefficients. Proc of Saratov University, New ser. Ser Math Mech and Inf 10: 3-9.) by using a new integral representations of the special solutions of Eq. (1).

The spectral analysis of the boundary value problem (1)-(2) was examined in Adiloglu and Amirov (20131 ADILOGLU NA AND AMIROV RKH. 2013. On the boundary value problem for the Sturm-Liouville equation with the discontinuous coefficient. Math Methods Appl Sci 36: 1685-1700.) where useful integral representations for two linearly independent solutions of equation (1) were constructed (see also Akhmedova and Huseynov 2010), the asymptotic formulas for the eigenvalues and eigenfunctions were obtained, completeness and expansion theorems for the system of the eigenfunctions were proved. Using the results of Adiloglu and Amirov (2013), in the present paper, we investigate uniqueness properties of the solution of the inverse problems for the problem (1)-(2) and study some other types boundary value problems related with the equation (1). In section 2 we prove that the boundary-value problem (1)-(2) is uniquely reconstructed, i.e. the potential function q(x) and the constants h, h1 are uniquely determined by given Weyl function or by the given spectral data. We also show that in the special case the potential q(x) and the coefficient h can be determined by the one spectrum only. In section 3 we investigate the properties of the spectral characteristics of two boundary-value problems related to Eq. (1).

UNIQUENESS OF THE SOLUTION OF THE INVERSE PROBLEM

Let s(x,λ), c(x,λ) be solutions of Eq. (1) with initial conditions

s ( 0 , λ ) = c ( 0 , λ ) = 0 , s ( 0 , λ ) = c ( 0 , λ ) = 1

and φ(x,λ), ψ(x,λ) be solutions of (1) under conditions at π:

φ ( π , λ ) = ψ ( π , λ ) = - 1 , φ ( π , λ ) = ψ ( π , λ ) = 0 .

Then w1(x,λ)=hs(x,λ)+c(x,λ) and w2(x,λ)=ψ(x,λ)-h1φ(x,λ) are solutions of (1) with initial conditions

w 1 ( 0 , λ ) = 1 , w 1 ( 0 , λ ) = h

and

w 2 ( π , λ ) = - 1 , w 2 ( π , λ ) = h 1

respectively (see Adiloglu and Amirov 2013).

Denote by Φ(x,λ) the solution of equation (1) satisfying conditions

Φ ( 0 , λ ) - h Φ ( 0 , λ ) = 1 , Φ ( π , λ ) + h 1 Φ ( π , λ ) = 0 . (3)

We set M(λ)=Φ(0,λ) and consider the linearly independent solutions s(x,λ) and w1(x,λ) of equation (1). We have

Φ ( x , λ ) = s ( x , λ ) - ( 0 ) ( λ ) ( λ ) w 1 ( x , λ ) , (4)

where

( 0 ) ( λ ) = s ( π , λ ) + h 1 s ( π , λ ) ,
( λ ) = w 1 ( π , λ ) + h 1 w 1 ( π , λ ) = - w 2 ( 0 , λ ) + h w 2 ( 0 , λ ) .

Therefore,

M(λ)=Φ(0,λ)=-(0)(λ)(λ).(5)

Then Eq. (4) can be written as

Φ ( x , λ ) = s ( x , λ ) + M ( λ ) w 1 ( x , λ ) (6)

Additionally we see that the solution w2(x,λ) also satisfies the second one of the conditions (3). Consequently we obtain

Φ ( x , λ ) = - w 2 ( x , λ ) ( λ ) . (7)

We also have

w 1 ( x , λ ) , Φ ( x , λ ) := w 1 ( x , λ ) Φ ( x , λ ) - w 1 ( x , λ ) Φ ( x , λ ) 1 . (8)

The functions Φ(x,λ) and M(λ) are called the Weyl solution and the Weyl function of the boundary value problem (1)-(2), respectively. From equation (5) we have that the Weyl function M(λ) is a meromorphic function with simple poles at the points λ=λn, n0, where {λn2} are eigenvalues of the boundary value problem (1)-(2) (see Adiloglu and Amirov 2013).

We recall that the normalized numbers (see Adiloglu and Amirov 2013) αn of the boundary value problem (1)-(2) are defined as

α n = 0 π ρ ( x ) w 1 2 ( x , λ n ) d x , ( n = 0 , 1 , 2 , ) ,

the set {λn2,αn}n0 is the spectral data of the problem (1)-(2) . Note that there exists the sequence βn such that

β n α n = Δ ( λ n ) 2 λ n

(see Adiloglu and Amirov 2013).

The following theorem shows that the given spectral data uniquely determines the boundary- value problem (1)-(2).

Theorem 1.

The following formula holds M(λ)=n=01αn(λ-λn).

Proof.

Since (0)(λ)=s(π,λ)+h1s(π,λ)=w2(0,λ) it follows from (Adiloglu and Amirov (2013)) (see formula (50), Adiloglu and Amirov 2013) that |(0)(λ)|Ce|Imλ|μ+(π), where μ+(π)=απ-αa+a. Since |(λ)|Cδ~|λ|e|Imλ|μ+(π), λGδ, |λ|>ρ0>0, where, Gδ={λ:|λ-λn0|δ} and Cδ~>0 for some δ>0, we have

| M ( λ ) | C δ | λ | , λ G δ , | λ | ρ 0 (9)

Further using the Lemma 1 (Adiloglu and Amirov (2013)) (see also the formula (44) there) we find

R e s λ = λ n M ( λ ) = - ( 0 ) ( λ n ) ( λ n ) = - β n ( λ n ) = 1 α n . (10)

Let

I N ( λ ) = 1 2 π i Γ N M ( p ) λ - p 𝑑 p ,
λ D N := { λ : | λ | < ( N + 1 2 ) 2 } , Γ N = D N .

Then by virtue of (9) we have limNIN(λ)=0. On the other hand by the residue theorem

I N ( λ ) = - M ( λ ) + n = 0 N 1 α n ( λ - λ n )

which gives the desired results as N+. Theorem is proved. ∎

Now let the Weyl function M(λ) of the boundary value problem (1)-(2) is given. In the following theorem we prove that the boundary problem (1)-(2) is uniquely reconstructed, i.e. the potential function q(x) and the constants h, h1 are uniquely determined by given Weyl function.

Let us denote the boundary value problem (1)-(2) by L=L(q(x),h,h1) and the similar boundary value problem with the potential q~(x) and boundary constants h~,h~1 by L~=L(q~(x),h~,h~1). Then the following theorem is satisfied.

Theorem 2.

Let the Weyl functions of the boundary value problem L and L~ are M(λ) and M~(λ) respectively. If M(λ)=M~(λ) then L=L~.

Proof.

Define the matrix P(x,λ)=[Pik(x,λ)]i,k=1,2 by the formula

P(x,λ)[w~1(x,λ)Φ~(x,λ)w~1(x,λ)Φ~(x,λ)]=[w1(x,λ)Φ(x,λ)w1(x,λ)Φ(x,λ)](11)

where w~1(x,λ) and Φ~(x,λ) are solutions of the equation of the boundary value problem L~ which is identical to the solutions w1(x,λ) and Φ(x,λ). Since w~1(x,λ),Φ~(x,λ)1 we have

{w1(x,λ)=P11(x,λ)w~1(x,λ)+P12(x,λ)w~1(x,λ)Φ(x,λ)=P11(x,λ)Φ~(x,λ)+P12(x,λ)Φ~(x,λ)(12)

or

{ P 11 ( x , λ ) = w 1 ( x , λ ) Φ ~ ( x , λ ) - w ~ 1 ( x , λ ) Φ ( x , λ ) P 12 ( x , λ ) = w ~ 1 ( x , λ ) Φ ( x , λ ) - w 1 ( x , λ ) Φ ~ ( x , λ ) . (13)

Using the formula (7) in equation (13) we find

P 11 ( x , λ ) = 1 + 1 ( λ ) { w 2 ( x , λ ) [ w ~ 1 ( x , λ ) - w 1 ( x , λ ) ] -
= - w 1 ( x , λ ) [ w ~ 2 ( x , λ ) - w 2 ( x , λ ) ] } , (14)
P 12 ( x , λ ) = 1 ( λ ) { w 1 ( x , λ ) w ~ 2 ( x , λ ) - w ~ 1 ( x , λ ) w 2 ( x , λ ) } (15)

Now using the estimations (49),(50) and (84) in (Adiloglu and Amirov (2013)) we have

| P 11 ( x , λ ) - 1 | c δ | λ | , | P 12 ( x , λ ) | c δ | λ | , λ G δ , | λ | ρ 0 (16)

On the other hand from the equation (6) we obtain

P 11 ( x , λ )
= w 1 ( x , λ ) s ~ ( x , λ ) - w ~ 1 ( x , λ ) s ( x , λ ) +
+ [ M ~ ( λ ) - M ( λ ) ] w 1 ( x , λ ) w ~ 1 ( x , λ )
P 12 ( x , λ )
= w ~ 1 ( x , λ ) s ( x , λ ) - w 1 ( x , λ ) s ~ ( x , λ ) +
+ [ M ( λ ) - M ~ ( λ ) ] w 1 ( x , λ ) w ~ 1 ( x , λ )

Thus if M(λ)=M~(λ) then the functions P11(x,λ) and P12(x,λ) are entire in λ. Together with (16) this gives P11(x,λ)1, P12(x,λ)0. Then from (12) we obtain w1(x,λ)=w~1(x,λ), Φ(x,λ)=Φ~(x,λ) for all x and λ. Since w1(x,λ) and w~1(x,λ) satisfy the equations

- w 1 ′′ + q ( x ) w 1 = λ 2 ρ ( x ) w 1 and w ~ 1 ′′ + q ~ ( x ) w ~ 1 = λ 2 ρ ( x ) w ~ 1

correspondingly, substructing these equations we have

( q ( x ) - q ~ ( x ) ) w 1 ( x , λ ) = 0 .

Therefore q(x)=q~(x) a.e. on (0,π). Further, because of Φ(x,λ)=Φ~(x,λ) for all x and λ from the conditions (3) we also have that h=h~, h1=h1..~Consequently L=L~. Theorem is proved. ∎

Theorem 3.

Let S={λn2,αn}n0 and S~={λ~n2,α~n}n0 are the spectral data of the problems L and L~ respectively. If λn=λ~n, αn=α~n, n0 then L=L~.

Proof.

It is known (Adiloglu and Amirov (2013)) that the solutions c(x,λ) and s(x,λ) of equation (1) satisfying the conditions

c ( 0 , λ ) = s ( 0 , λ ) = 1 , c ( 0 , λ ) = s ( 0 , λ ) = 1

are expressed as

c ( x , λ ) = c 0 ( x , λ ) + 0 μ + ( x ) N + ( x , t ) cos λ t d t (17)
s ( x , λ ) = s 0 ( x , λ ) + 0 μ + ( x ) N - ( x , t ) sin λ t λ 𝑑 t (18)

respectively, where

c 0 ( x , λ )
= r + ( x ) cos λ μ + ( x ) + r - ( x ) cos λ μ - ( x )
s 0 ( x , λ )
= r + ( x ) sin λ μ + ( x ) λ + r - ( x ) sin λ μ - ( x ) λ
r ± ( x ) = 1 2 ( 1 ± 1 ρ ( x ) ) , μ ± ( x ) = ± x ρ ( x ) + a ( 1 ± ρ ( x ) )

and N±(x,)L1(0,μ+(x)) for each x[0,π]. Using (17) and (18) we have

w 1 ( x , λ ) = w 1 ( 0 ) ( x , λ ) + 0 μ + ( x ) W 1 ( x , t ) cos λ t d t , (19)

where

w 1 ( 0 ) ( x , λ ) = c 0 ( x , λ ) + h s 0 ( x , λ ) (20)
W 1 ( x , t ) = N + ( x , t ) + h t μ + ( x ) N - ( x , ξ ) 𝑑 ξ . (21)

Clearly W1(x,)L1(0;μ+(x)) for all x[0,π]. Note that the function W1(x,t) also satisfies the conditions (see Adiloglu and Amirov 2013)

d d x W 1 ( x , μ + ( x ) ) = 1 4 ρ ( x ) ( 1 + 1 ρ ( x ) ) q ( x ) , (22)
d d x { W 1 ( x , μ - ( x ) + 0 ) - W 1 ( x , μ - ( x ) - 0 ) } = 1 4 ρ ( x ) ( 1 - 1 ρ ( x ) ) q ( x ) . (23)

Since

s 0 ( x , λ ) = 0 μ + ( x ) r ( x , t ) cos λ t d t ,

where r(x,t)={1, 0<t<μ-(x)α+, μ-(x)<t<μ+(x)the equation (19) can be written as

w 1 ( x , λ ) = c 0 ( x , λ ) + 0 μ + ( x ) W ( x , t ) cos λ t d t , (24)

where

W ( x , t ) = W 1 ( x , t ) + h r ( x , t ) . (25)

It is clear that if 0xa then

w 1 ( x , λ ) = cos λ x + 0 x W ( x , t ) cos λ t d t , (26)

where W(x,t)=W1(x,t)+h is continuous kernel. Then we can see the equation (26) as a Volterra integral equation with respect to cosλx. From the theory of the Volterra integral equations, we know that the equation (26) is then uniquely solvable and the solution is

cos λ x = w 1 ( x , λ ) + 0 x W ^ ( x , t ) w 1 ( t , λ ) 𝑑 t (27)

where W1^(x,t) is a continuous kernel.

Let now x>a. In this case the equation (24) is written as

w 1 ( x , λ ) = α + cos λ μ + ( x ) + α - cos λ ( 2 a - μ + ( x ) ) + + 0 μ + ( x ) W ( x , t ) cos λ t d t (28)

where α±=12(1±1α). Here the kernel W(x,t) has a jump discontinuity at t=μ-(x). Clearly, μ+(x)>a and 0<μ-(x)<a when x>a and therefore (28) takes the form

w 1 ( x , λ ) = α + cos λ μ + ( x ) + a μ + ( x ) W ( x , t ) cos λ t d t + h ( x , λ ) , (29)

where

h ( x , λ ) = α - cos λ μ - ( x ) + 0 a W ( x , t ) cos λ t d t (30)

Now using (27) we obtain that

h ( x , λ ) = α - w 1 ( μ - ( x ) , λ ) + α - 0 μ - ( x ) W ^ ( μ - ( x ) , t ) w 1 ( t , λ ) 𝑑 t
- α - 0 μ - ( x ) ( a x U ( x , μ + ( s ) ) W ^ ( μ - ( s ) , t ) 𝑑 s ) w 1 ( t , λ ) 𝑑 t
+0a(W(x,t)+taW(x,s)W^(s,t)𝑑s)w1(t,λ)𝑑t,(31)

where U(x,μ+(t)) is a continuous kernel. Therefore the equation (29) is written as

w1(x,λ)-h(x,λ)=α+cosλμ+(x)+α+axW(x,μ+(t))cosλμ+(x)𝑑t,(32)

where the kernel W(x,μ+(t)) is continuous. Hence, we obtain the Volterra integral equation

w 1 ( x , λ ) - h ( x , λ ) = ( I + W ) α + cos λ μ + ( x ) , (33)

where

( I + W ) f ( x ) = f ( x ) + a x W ( x , μ + ( t ) ) f ( t ) 𝑑 t (34)

Solving the Volterra integral equation (33) with respect to α+cosλμ+(x), we find

α + cos λ μ + ( x ) = w 1 ( x , λ ) - h ( x , λ ) + a x U ( x , μ + ( t ) ) [ w 1 ( t , λ ) - h ( t , λ ) ] 𝑑 t . (35)

Taking into account, the expression for h(x,λ), we have

cos λ μ + ( x ) = 1 α + w 1 ( x , λ ) - α - α + w 1 ( μ - ( x ) , λ ) + 0 x W ^ ( x , t ) w 1 ( t , λ ) 𝑑 t (36)

where W^(x,t) is a kernel with jump at t=μ-(x). Consequently, we have

w 1 ( x , λ ) = w 1 ~ ( x , λ ) + 0 x H ( x , t ) w 1 ~ ( t , λ ) 𝑑 t , (37)

where H(x,t) is a kernel with jump at t=μ-(x).

Let f(x)L2(0,π). (37) implies that

0 π ρ ( x ) f ( x ) w 1 ( x , λ ) 𝑑 x = 0 π ρ ( x ) g ( x ) w 1 ~ ( x , λ ) 𝑑 x ,

where g(x)=f(x)+xπH(t,x)f(t)𝑑t. Hence for all n0

a n = b n ~ ,

where an:=0πρ(x)f(x)w1(x,λn)𝑑x and bn~:=0πρ(x)g(x)w1~(x,λn)𝑑x. Using the Parseval’s equality (Adiloglu and Amirov (2013)), we calculate

f L 2 ( 0 , π , ρ ) = g L 2 ( 0 , π , ρ ) . (38)

Now if we consider the operator

A f ( x ) = f ( x ) + x π H ( t , x ) f ( t ) 𝑑 t , (39)

we have Af=g and AfL2(0,π;ρ)=gL2(0,π;ρ). Then by (38) we obtain

A f L 2 = f L 2

for any f(x)L2(0,π). Then A=A-1 is possible if and only if H(t,x)0. Thus w1(x,λ)=w1~(x,λ), i.e. q(x)=q~(x) a.e. on (0,π) and h=h~, H=H~. Theorem is proved. ∎

Let q(x)=α2q(απ-αx) for x>a and h1=αh. We now show that in this case the potential q(x) and coefficient h can be determined by the spectrum {λn2}n0 only.

Theorem 4.

If q(x)=α2q(απ-αx), x>a, h1=αh, q~(x)=α2q~(απ-αx), x>a, h1~=αh~and λn=λn~, n0 then q(x)=q~(x) a.e. on (0,π) and h=h~.

Proof.

If q(x)=α2q(απ-αx), x>aand y(x) is a solution of the Eq. (1) for 0xa, then y1(x):=y(απ-αx) is a solution of (1) for x>a.

Indeed, it is easy to check that

- y 1 ′′ ( x ) = ρ ( x ) λ 2 y 1 ( x ) , i.e.

y1(x) is a solution for x>a. In particular, if we take the solution

y ( x ) = - w 1 ( x , λ ) , 0 x a

then y1(x)=-w1(απ-αx,λ), (x>a) is the solution of (1) satisfying the initial conditions y1(π)=-w1(0,λ)=-1, y1(π)=αw1(0,λ)=αh=h1.

Consequently, w2(x,λ)-w1(απ-αx,λ). Now using w2(x,λn)=βnw1(x,λn), βn0 we have

w 2 ( x , λ n ) = - 1 β n w 1 ( x , λ n )

which implies βn2=1. On the other hand βnw1(π,λn)=-1 and

w 1 ( π , λ n ) = c 0 ( π , λ n 0 ) + ξ n ( π ) n , | ξ n ( π ) | C , (see Adiloglu and Amirov 2013).

Consequently we have βn=(-1)n+1, hence from the formula βnαn=-Δ(λn)2λn we obtain

α n = ( - 1 ) n Δ ( λ n ) 2 λ n .

Since λn=λn~ we have αn=αn~. Then by the previous theorem q(x)=q~(x) a.e. on (0,π) and h=h~. Theorem is proved. ∎

BOUNDARY VALUE PROBLEMS L0 AND L1

(i) Consider the boundary value problem L1=L1(q(x),h) for equation (1) with the boundary conditions

y ( 0 ) - h y ( 0 ) = 0 , y ( π ) = 0 .

The characteristic function of the problem L1 is d(λ)=w1(π,λ) and the eigenvalues of L1 are the squares of zeros of the equation w1(π,λ)=0.

Note that as in the case of the problem L we can prove that the eigenvalues {μn2} of the problem L1 are real and simple.

Since

w 1 ( π , λ ) = c 0 ( π , λ ) + O ( e | Im λ μ + ( π ) | | λ | ) , | λ |

we have for the eigenvalues {μn2} the following asymptotic formula:

μ n = μ n 0 + ε n , (40)

where

μ n 0 = π μ + ( π ) ( n + 1 2 ) + θ n , ( sup n | θ n | < + )

are the roots of the equation c0(π,λ)=0 and εn=o(1), n.

Further it is easy to show that (see Adiloglu and Amirov 2013 also) the solution w1(x,λ) satisfies the asymptotic relation

w 1 ( x , λ ) = c 0 ( x , λ ) + h s 0 ( x , λ ) + 0 x g ( x , t , λ ) c 0 ( t , λ ) q ( t ) 𝑑 t + O ( e | Im λ μ + ( π ) | λ 2 ) , λ

where

g ( x , t , λ ) = s 0 ( x , λ ) c 0 ( t , λ ) - c 0 ( x , λ ) s 0 ( t , λ ) , 0 x π , 0 t π .

Therefore we obtain

d ( λ ) = w 1 ( π , λ ) = α + cos λ μ + ( π ) + α - cos λ μ - ( π )
+ b 1 α + sin λ μ + ( π ) λ + b 2 α - sin λ μ - ( π ) λ + k ( λ ) λ (41)

where

b 1 = h + 1 2 0 π q ( t ) ρ ( t ) 𝑑 t , (42)
b 2 = h + 1 2 0 π s g n ( a - t ) ρ ( t ) q ( t ) 𝑑 t , (43)
k ( λ ) = - μ + ( π ) μ + ( π ) Q ( t ) sin λ t d t + O ( e | Im λ μ + ( π ) | λ ) , λ (44)

with Q(t)L2(-μ+(π),μ+(π)).

Since d(μn)=0, equations (40), and (41) imply that

( μ n 0 + ε n ) ( α cos ( μ n 0 + ε n ) μ + ( π ) + α - cos ( μ n 0 + ε n ) μ - ( π ) )
+ b 1 α + sin ( μ n 0 + ε n ) μ + ( π ) + b 2 α - sin ( μ n 0 + ε n ) μ - ( π ) + k n = 0

and consequently

μ n 0 [ α + μ + ( π ) sin μ n 0 μ + ( π ) + α - μ - ( π ) sin μ n 0 μ - ( π ) ] ε n
= μ + ( π ) b 1 α + cos μ n 0 μ + ( π ) + μ - ( π ) b 2 α - cos μ n 0 μ - ( π )
+ b 1 α + sin μ n 0 μ + ( π ) + b 2 α - sin μ n 0 μ - ( π ) + k n , (45)

which implies

ε n = d n μ n 0 + k n n , (46)

where

d n = b 1 α + sin μ n 0 μ + ( π ) + b 2 α - sin μ n 0 μ - ( π ) α + μ + ( π ) sin μ n 0 μ + ( π ) + α - μ - ( π ) sin μ n 0 μ - ( π ) , (47)

knl2. Hence

μ n = μ n 0 + d n μ n 0 + k n n . (48)

Moreover if we define the normalized numbers (αn1) for the problem L1 as

α n 1 = 0 π ρ ( x ) w 1 2 ( x , μ n ) 𝑑 x (49)

then we have

α n 1 = α n 1 0 + t n 1 n , (50)

where

α n 1 0 = 0 π ρ ( x ) c 0 2 ( x , μ n ) 𝑑 x and t n 1 l 2 .

Since d(λ) is an entire function of order one by the H’Adamard’s theorem the function d(λ) is uniquely determined up to a multiplicative constant C by its zeros:

d ( λ ) = C n = 0 ( 1 - λ 2 μ n 2 ) . (51)

We also have

c 0 ( π , λ ) = α + cos λ μ + ( π ) + α - cos λ μ - ( π ) = n = 0 ( 1 - λ 2 ( μ n 0 ) 2 ) . (52)

Since

lim | λ | d ( λ ) c 0 ( π , λ ) = 1 , lim | λ | n = 0 ( 1 + μ n 2 - ( μ n 0 ) 2 ( μ n 0 ) 2 - λ 2 ) = 1 ,

we have

C = n = 0 μ n 2 ( μ n 0 ) 2 ,

therefore

d(λ)=n=0μn2(μn0)2(1-λ2μn2)=n=0(μn2-λ2(μn0)2).(53)

We have proved the following two theorems:

Theorem 5.

The boundary value problem L1= L1(q(x),h) has a countable set of eigenvalues {μn2} n1and for sufficiently large values of n,the asymptotic formula

μ n = μ n 0 + d n μ n 0 + k n n

are satisfied,where

μ n 0 = π μ + ( π ) ( n + 1 2 ) + θ n , ( sup n | θ n | < + ) ,
d n = b 1 α + sin μ n 0 μ + ( π ) + b 2 α - sin μ n 0 μ - ( π ) α + μ + ( π ) sin μ n 0 μ + ( π ) + α - μ - ( π ) sin μ n 0 μ - ( π ) ,

is bounded and kn2. Moreover, we have

α n 1 = α n 1 0 + t n 1 n ,

for the normalized numbers

α n 1 = 0 π ρ ( x ) w 1 2 ( x , μ n ) 𝑑 x

o f the problem L1, where

α n 1 0 = 0 π ρ ( x ) c 0 2 ( x , μ n 0 ) 𝑑 x and t n 1 l 2 .

Theorem 6.

The specification of the spectrum {μn2}n0uniquely determines the characteristic function d(λ)=w1(π,λ)by the formula

d ( λ ) = n = 0 ( μ n 2 - λ 2 ( μ n 0 ) 2 ) .

(ii) Consider the boundary value problem L0=L0(q(x),h1) for equation (1) with the boundary condition

y ( 0 ) = 0 , y ( π ) + h 1 y ( π ) = 0 .

The eigenvalues {(λn0)2}n0 of the problem L0 are simple and coincide with the squares of zeros of the characteristic function

Δ 0 ( λ ) := w 2 ( 0 , λ ) = s ( π , λ ) + h 1 s ( π , λ )

where s(x,λ) is the solution of Eq.(1) with the initial conditions s(0,λ)=0, s(0,λ)=1.

From the results of Adiloglu and Amirov (2013), we have

s ( x , λ ) = s 0 ( x , λ ) + 0 x g ( x , t , λ ) q ( t ) s ( t , λ ) 𝑑 t . (54)

Since

s ( x , λ ) = s 0 ( x , λ ) + O ( e | Im λ μ + ( x ) | λ 2 ) , | λ |
s ( x , λ )
= s 0 ( x , λ ) + α 0 a [ α + cos λ ( μ + ( x ) - t ) - α - cos λ ( μ - ( x ) - t ) ] sin λ t λ q ( t ) 𝑑 t
+ a x cos λ ( μ + ( x ) - μ + ( t ) ) q ( t ) 𝑑 t s 0 ( x , λ ) + O ( e | Im λ μ + ( x ) | λ 2 ) , | λ | ,

We have

Δ 0 ( λ ) = α α + 2 cos λ μ + ( π ) - α α - 2 cos λ μ - ( π )
+ α + w + 0 sin λ μ + ( π ) λ + α - w - 0 sin λ μ - ( π ) λ + k 0 ( λ ) λ (55)

where

w ± 0 = h 1 ± α 2 0 a q ( t ) 𝑑 t + 1 2 a π q ( t ) 𝑑 t (56)
k 0 ( λ ) = - α α + 2 0 π q ( t ) ρ ( t ) sin λ ( μ + ( π ) - 2 μ + ( t ) ) 𝑑 t + α α - 2 0 a q ( t ) sin λ ( μ - ( π ) - 2 t ) 𝑑 t
+ α α - 2 0 π q ( t ) sin λ ( 2 μ - ( t ) - μ - ( π ) ) 𝑑 t + O ( e | Im λ μ + ( π ) | | λ | ) . (57)

Therefore

Δ 0 ( λ ) = α 2 ( α + cos λ μ + ( π ) - α - cos λ μ - ( π ) )
+ α + w + 0 sin λ μ + ( π ) λ + α - w - 0 sin λ μ - ( π ) λ + k 0 ( λ ) λ (58)

From (58) we obtain the following expression for the roots {λn0} of the function Δ0(λ):

λ n 0 = γ n 0 + η n , (59)

where

γ n 0 = π μ + ( π ) ( n + 1 2 ) + θ n 0
sup n | θ n 0 | < + and η n = o ( 1 ) , n .

Since Δ0(λn0)=0 some simple transformations lead to

η n = d n 0 γ n 0 + k n 0 n , k n 0 l 2 , (60)
d n 0 = α + w + 0 sin μ + ( π ) γ n 0 + α - w - 0 sin μ - ( π ) γ n 0 α α + μ + ( π ) sin μ + ( π ) γ n 0 + α α - μ - ( π ) sin μ - ( π ) γ n 0 (61)

where (kn)l2. Therefore

λ n 0 = γ n 0 + d n 0 γ n 0 + k n 0 n , k n 0 l 2 (62)

Moreover we can obtain that

Δ 0 ( λ ) = n = 0 ( λ n 0 ) 2 - λ 2 ( γ n 0 ) 2 (63)

We can formulate the following theorem for the boundary- value problem L0:

Theorem 7.

The boundary value problem L0= L0(q(x),h1) has a countable set of eigenvalues {(λn0)2} n1and for sufficiently large values of nthe asymptotic formula (62)is satisfied. Moreover, the characteristic function

Δ 0 ( λ ) := w 2 ( 0 , λ ) = s ( π , λ ) + h 1 s ( π , λ )

is uniquely determined by the spectrum {(λn0)2}n0 via the formula (63).

Lemma 1.

The following relation holds

λ n < μ n < λ n + 1 , n 0 , (64)

i.e. the eigenvalues of two boundary problems L and L1 are alternating.

Proof.

Consider the characteristic functions Δ(λ) and d(λ) of the boundary value problems L and L1 respectively. Let λ2=z. Then

Δ ( z ) = w 1 ( π , z ) + h 1 w ( π , z ) , d ( λ ) = w 1 ( π , z ) .

Since

( z - μ ) 0 π ρ ( x ) w 1 ( x , z ) w 1 ( x , μ ) 𝑑 x
= ( Δ ( μ ) - Δ ( z ) ) d ( z ) - ( d ( μ ) - d ( z ) ) Δ ( z ) .

for μz we obtain

0 π ρ ( x ) w 1 2 ( x , z ) 𝑑 x = Δ ( z ) d ( z ) - Δ ( z ) d ( z ) , (65)

where

Δ ( z ) = d d z Δ ( z )

and

d ( z ) = d d z d ( z )

From (65) we have for z=λn2

α n = - Δ ( λ n ) d ( λ n ) 2 λ n (66)

and

1 d 2 ( z ) 0 π ρ ( x ) w 1 2 ( x , z ) 𝑑 x = - d d z ( Δ ( z ) d ( z ) ) ,

where λ=(z) is real and d(z)0.

Thus the function Δ(z)d(z) is monotonically decreasing when zμn2, n0 and

lim z ( μ n ± 0 ) 2 Δ ( z ) d - 1 ( z ) = ±

Then from asymptotic formulas for λn and μn, we arrive at (64).

ACKNOWLEDGMENTS

We thank the referees for their encouraging remarks and insightful comments. Also, this work is supported by the Scientific and Technological Research Council of Turkey (TÜBİTAK) with Project 113F366.

REFERENCES

  • 1
    ADILOGLU NA AND AMIROV RKH. 2013. On the boundary value problem for the Sturm-Liouville equation with the discontinuous coefficient. Math Methods Appl Sci 36: 1685-1700.
  • 20
    AKHMEDOVA EN AND HUSEYNOV HM. 2010. On inverse problem for Sturm-Liouville operator with discontinuous coefficients. Proc of Saratov University, New ser. Ser Math Mech and Inf 10: 3-9.
  • 3
    AMIROV RKH. 2006. On Sturm-Liouville operators with discontinuity conditions inside an interval. J Math Analysis Appl 317: 163-176.
  • 4
    ANDERSSON L. 1988. Inverse eigenvalue problems for a Sturm-Liouville equation in impedance form. Inverse Probl 4: 929-971.
  • 5
    CARLSON R. 1994. An inverse spectral problem for Sturm-Liouville operators with discontinuous coefficients. Proc Amer Math Soc 120: 475-484.
  • 6
    FREILING G AND YURKO V. 2001. Inverse Sturm-Liouville problems and their applications, Nova Science Publishers, Inc., Huntington: NY, 305 p.
  • 7
    GASYMOV MG. 1977. The direct and inverse problem of spectral analysis for a class of equations with a discontinuous coefficient. Non-Classical Methods in Geophysics, 37-44, Nauka, Novosibirsk, Russia.
  • 8
    GUSEINOV IM AND PASHAEV RT. 2002. On an inverse problem for a second-order differential equation. Russian Math Surveys 57: 597-598.
  • 9
    HALD O. 1980. Inverse eigenvalue problems for the mantle. Geophys J R Astron Soc 62: 41-48.
  • 10
    HALD O. 1984. Discontinuous inverse eigenvalue problems. Commun Pure Appl Math 37: 539-577.
  • 11
    HRYNIV RO AND MYKYTYUK YAV. 2003. Inverse spectral problems for Sturm-Liouville operators with singular potentials. Inverse Probl 19: 665-684.
  • 12
    HRYNIV RO AND MYKYTYUK YAV. 2004. Half-inverse spectral problems for Sturm-Liouville operators with singular potentials. Inverse Probl 20: 1423-1444.
  • 13
    KRUEGER R. 1982. Inverse problems for nonabsorbing media with discontinuous material properties. J Math Phys 23: 396-404.
  • 14
    LEVITAN BM. 1987. Inverse Surm-Liouville problems, VNU Sci. Press, Utrecht, 240 p.
  • 15
    LEVITAN BM AND GASYMOV MG. 1964. Determination of a differential equation by two spectra. Uspehi Mat Nauk 19: 3-63.
  • 16
    MAMEDOV KHR. 2006. Uniqueness of the solution of the inverse problem of scattering theory for Sturm-Liouville operator with discontinuous coefficient. Proc Inst Math Mech Natl Acad Sci Azerb 24: 163-172.
  • 17
    MAMEDOV KHR. 2010. On an inverse scattering problem for a discontinuous Sturm-Liouville equation with a spectral parameter in the boundary condition. Bound Value Probl ID 171967, 17 p.
  • 18
    MAMEDOV KHR AND PALAMUT N. 2009. On a direct problem of scattering theory for a class of Sturm-Liouville operator with discontinuous coefficient. Proc Jangjeon Math Soc 12: 243-251.
  • 19
    WILLIS C. 1984. Inverse problems for torsional modes. Geophys J R Astron Soc 78: 847-853.
  • 20
    MARCHENKO VA. 2011. Sturm-Liouville operators and applications, AMS Chelsea Publishing Volume 373, 393 p.
  • 21
    YURKO VA. 2000. On boundary value problems with discontinuity conditions inside an interval. Differ Equ 36: 1266-1269.

Publication Dates

  • Publication in this collection
    Dec 2017

History

  • Received
    05 Feb 2016
  • Accepted
    27 Jan 2017
Academia Brasileira de Ciências Rua Anfilófio de Carvalho, 29, 3º andar, 20030-060 Rio de Janeiro RJ Brasil, Tel: +55 21 3907-8100 - Rio de Janeiro - RJ - Brazil
E-mail: aabc@abc.org.br