О восстановлении интегродифференциальных операторов по функции Вейля Текст научной статьи по специальности «Математика»

УДК 517.984

О ВОССТАНОВЛЕНИИ ИНТЕГРОДИФФЕРЕНЦИАЛЬНЫХ ОПЕРАТОРОВ ПО ФУНКЦИИ ВЕЙЛЯ

М. Ю. Игнатьев1, С. Ю. Советникова2

1 Игнатьев Михаил Юрьевич, кандидат физико-математических наук, доцент кафедры математической физики и вычислительной математики. Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, 410012, Россия, Саратов, Астраханская, 83, IgnatievMU@info.sgu.ru

2Советникова Светлана Юрьевна, старший преподаватель кафедры математической физики и вычислительной математики. Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, 410012, Россия, Саратов, Астраханская, 83, sovetnikovasy@mail.ru

Исследуются обратные спектральные задачи для интегродифференциальных операторов второго порядка, которые являются возмущением оператора Штурма-Лиувилля интегральным вольтерровским оператором. Основное внимание уделяется нелинейной обратной задаче восстановления потенциала по заданной функции Вейля при условии, что ядро интегрального оператора известно априори. Получены свойства спектральных характеристик и функции Вейля, приведен алгоритм решения обратной задачи и установлена единственность решения. Для решения обратной задачи используется метод эталонных моделей.

Ключевые слова: интегродифференциальные операторы, обратные спектральные задачи, единственность, алгоритм.

001:10.18500/1816-9791 -2017-17-3-276-284

1. Исследуются обратные задачи спектрального анализа для следующего инте-гродифференциального оператора:

Вводится и изучается функция Вейля для этого класса интегродифференциальных операторов. Получен алгоритм построения потенциала д по заданной функции Вейля при условии, что ядро М известно априори. Обратные спектральные задачи состоят в восстановлении операторов по их спектральным характеристикам. Такие задачи часто возникают в математике, механике, физике, электронике, геофизике и других областях естествознания и техники. Обратная задача для дифференциальных операторов достаточно полно изучена (см. монографии [1-6] и библиографию в них). Для интегродифференциальных и других классов нелокальных операторов обратные задачи являются более трудными для исследования; для таких операторов классические методы (метод оператора преобразования и метод спектральных отображений [1-6]) либо совсем не работают, либо требуют существенной модификации. Поэтому для интегродифференциальных операторов теория решения обратных

х

1у := -у''(х) + д(х)у(х) + М(х, г)у(г) (г, х е [0, т]

0

на конечном интервале с краевыми условиями:

у'(0) - Ну(0) = у'(Т) + Ну(Т) =0.

задач еще не построена. В то же время нелокальные, в частности интегродифферен-циальные операторы, вызывают интерес исследователей, так как они имеют много приложений (см., например, [7]). Отметим, что некоторые аспекты теории решения обратных задач для интегродифференциальных операторов исследовались в [8-11] и др. работах.

2. Объектом исследования является краевая задача L = L(q, M, h, H) для инте-гродифференциального уравнения:

x

ly := -y''(x) + q(x)y(x) + y M(x,t)y(t) dt = Ay(x), x G [0,T], (1)

0

с краевыми условиями:

U(y) := y'(0) - hy(0) = 0, V(y) := y'(T) + Hy(T) = 0,

где q(x) и M(x, t) — интегрируемые комплекснозначные функции, а h и H — комплексные числа. Пусть функции C(x,A), S(x,A) и ^(x,A) являются решениями уравнения (1) при начальных условиях ц(0, A) = C(0, A) = S'(0,A) = 1, C'(0, A) = S(0, A) = 0, ц'(0, A) = h. При каждом фиксированном x G [0,T] функции C)(x,A), S)(x, A) и )(x, A), v = 0,1, являются целыми по A порядка 1/2. Ясно, что

<p(x, A) = C (x, A) + hS (x, A).

Обозначим A(A) := V(ц) = ц'(Т, A)+Нц(Т, A). Нули {An}n^0 целой функции A(A) совпадают с собственными значениями краевой задачи L. Функция A(A) называется характеристической функцией для L.

Пусть функция Ф^^) является решением уравнения (1) при условиях U(Ф) = 1, V(Ф) = 0. Обозначим N(A) := Ф(0, A). Тогда

Ф^, A) = S (x,A) + N (A)^x, A), N (A) = -Ai(A)/A(A), (2)

где A1(A) := V(S) = S'(T, A) + HS(T, A). Нули {pn}n^0 целой функции A1(A) совпадают с собственными значениями краевой задачи Li для уравнения (1) с краевыми условиями y(0) = V (y) = 0. Функция N (A) называется функцией Вейля. Из (2) следует, что функция N(A) является мероморфной по A с полюсами {An}n^0 и нулями {pn}n^0. Пусть M(x,t),h и H известны априори. Обратная задача ставится следующим образом.

Обратная задача 1. Дана N(A), построить q(x).

Эта обратная задача является аналогом классической обратной задачи восстановления оператора Штурма - Лиувилля по заданной функции Вейля [3].

3. Пусть A = р2, т = Imр. Известным методом (см., например, [3]) получаем, что при |р| ^ го имеют место асимптотические формулы:

x

sinpx cospx f ,, 1 /exp(|t|x)\

s (xA)=-f-- j q(t) dt+0 v^^;,

0

С(х, Л) = ео—рх + ——— [ д(г) М + о е р(|т|х)

2Р .) \ Р

о

Б'(х, Л) = еоэ рх + —-—— [ д(г) М + о е Р(|т|х)

2Р .) \ Р

о

х

С (х,Л) = — р бш рх + ео—Р д(г) (г + о(ехр(|т |х))

и, следовательно,

Д(Л) = — рт + ш ео— рт + о(ехр(|т |Т)),

л т ш1—1прТ (ехр(|г|т)^ (3)

Д1 (Л) = ео— рТ +---Ь о 1 1

где

т т

ш := н + н +1У д(г) (г, ш1 := н +1У д(г) (г.

о о

Используя (3), вычисляем

/— ПП Ш /1\ - П / 1 \ Ш1 /1\

рп := V Лп = + - + о - Ь рп1 := Л/ДП = ^ п + о + - + о " Ь

Т пп \и) Т \ 2/ пп \и 1

Кроме того, задание чисел {Лп}п^0 и {дп}п^0 однозначно определяет характеристические функции Д(Л) и Д1(Л) по формулам [3]:

те

Л Л Л

Д(Л)= Т(Л0 — Л)П^ ' Д'(Л)^((П +^П/2)П/Т)2 . (4)

Учитывая (2) и (4), заключаем, что обратная задача 1 эквивалентна следующей обратной задаче.

Обратная задача 2. Даны два спектра {Лп}п^0 и {дп}п^0, построить д(х). Обозначим

-к

Г 2 := — У '(х)+ д(.ф (х)+/М (г)

х

Очевидно, что

! 1у(х) • г (х) = 0 (у г' — у'г) + ^ у (х) • Г г (х). (5)

00

Если у(х, Л) и г(х, д) — решения уравнений 1у = Лу и £*г = дг соответственно, то (5) дает

(Л — д) J у(х,Л)г(х,д) (х = (уг' — у'г). (6)

о

(8)

Пусть С*(ж,А),5*(ж,Л) и ф*(ж, А) — решения уравнения

I* г = Аг (7)

при начальных условиях:

С*(Т, А) = -5*(Т, А) = ф*(Т, А) = 1, С*(п, А) = Я(п, А) = 0, ф*(Т, А) = -Я. Тогда

V (ф*) = 0, ф* (ж, А) = С*(ж, А) + Я5*(ж, А). Из (6) при д = А вытекает, что

5 (п, А) = 5*(0, А), (п, А) = С* (0, А), С (п, А) = -5* (0, А), С' (п, А) = -С* (0, А).

Обозначим

Д*(А) := -и(ф*) = -ф*(0, А) + Дф*(0, А), А*(А) := ф*(0, А).

Используя (8), вычисляем

А1(А) = С*(0, А) + Я5* (0, А) = 5' (Т,А)+ Я5 (Т,А) = Ах (А), А*(А) = -С*(0, А) - Я5*(0, А) + ДС*(0, А) + ДЯ5*(0, А) = = С' (Т, А) + ЯС (Т, А) + (Т, А) + ДЯ5 (Т, А) = ^ (Т, А) + Я^(Т, А) = Ах (А)

и, следовательно,

А1(А) = Ах (А), А* (А) = А(А). (9)

Пусть Ф*(ж,А) — решение уравнения (7) при условиях и(Ф*) = 1, V(Ф*) = 0. Обозначим N * (А) := Ф*(0,А). Тогда

Ф (ж А) = ф-^ А) ¡и*(А) _ ф*(0,А) _ А1 (А

Вместе с (9) это дает N * (А) = N (А).

Известно, что существует фундаментальная система решений (ж,р),у2(ж,р)}, 1тр ^ 0, ж е [0,п], уравнения (1) такая, что при |р| ^ го, V = 0,1, имеем:

(ж,р) = (гр)^ ехр(грж) + 0(р-1), у2^(ж,р) = (—¿р)^ ехр(-¿рж)(1 + 0(р-1)).

Аналогично существует фундаментальная система решений (ж,р),г2(ж,р)}, 1тр ^ 0, ж е [0,п], уравнения (7) такая, что при |р| ^ го, V = 0,1, имеем:

\ж,р) = (ip)v exp(ipx)(1 + O(p 1)), p) = (—ip)v exp(-ipx) + O(p-1 exp(-ipn)).

Зафиксируем г e (0,п/2). Обозначим Q := {p : argp e — £]}. Используя эти фундаментальные системы решений, получаем

Ф(ж, А) = Ai(p)yi(x,p) + A2(p)y2(ж,р), Ф*(ж, А) = A*(p)zi(ж,р) + A*(p)z2(ж,р),

где Ак(р) и А*к (р) не зависят от х. Коэффициенты Ак (р) и А*к(р) вычисляются с помощью краевых условий и(Ф) = 1, V(Ф) = 0, и(Фк) = 1, V(Фк) = 0. Это дает следующую асимптотику при р е Q, |р| ^ ж, равномерно по х е [0,п — г]:

Ф(х, Л) = (гр)-1 ехр(грх)(1 + 0(р-1)) + 0(р-1),

Фк(х, Л) = (гр)-1 ехр(грх)(1 + 0(р-1)). (

4. Получим теперь алгоритм решения обратной задачи 1. Для этого наряду с Ь рассмотрим краевую задачу Ь := Ь(Ы,с[,к,Н) того же вида, но с другими коэффициентами. Условимся, что если некоторый символ а обозначает объект, относящийся к Ь, то а будет обозначать аналогичный объект, относящийся к Ь, и а := а — а.

Лемма 1. Пусть М(х, г) = М(х, г), Н = Н, Н = Н. Тогда

т

[ ((х)Ф(х, Л)Фк(х, Л) ¿х = й(Л), (11)

где ((х) = ((х) — ¡(х), ЩЛ) = N (Л) — N (Л). Доказательство. Имеем

X

—Ф''(х,л) + ((х)Ф(х,Л) + ! м (х,г)Ф(г,л) ¿г = лф(х,л),

0 т

— 4^ + ((х)фк(х,Л)+/ м ^ (г,л) Л = Лфк(х,Л).

X

Умножим первое соотношение на Фк(х, Л), затем вычтем второе соотношение, умноженное на Ф(х, Л), и проинтегрируем по х :

т т X

J ((х)Ф(х,л)Фк(х,л) йх + у Фк (х,л) ¿х J м (х,г)Ф(г,л) ¿г—

0 0 0 т т

Ф(х, л) м (г,х)Фк (г,л) ¿г = 0 (ф' (х,л)Фк(х,л) — Ф(х,л)Ф'к(х,л)).

0 X

Учитывая соотношение N к(Л) = N (Л), приходим к (11). □

Лемма 2. Пусть

хк / \ / £(х р) \ п(х р) г(х) = — (7 + р(х)), Н(х, р) = ехр(2грх) ( 1+--'—]+ехр(грх)-'—, х е [0,Т],

К\ V / V р / р

где р(х) е С[0,Т], р(0) =0, а функции £(х,р), п(х,р) являются непрерывными и ограниченными при х е [0,Т], р е Q, |р| > рк. Тогда при |р| ^ ж, р е Q, имеем:

т

1

г(х)Н(х, р) ¿х = (—2гр)к+1 (7 + о(1)).

0

Доказательство. Вычисляем

т

(-2гр)^+1 / г(ж)Я(ж, р) ¿ж = Д(р) + /2(р)+ 1з(р) + Ш:

где

т

/ж°

— ехр(2грж) ¿ж,

о

т

/ж°

—р(ж) ехр(2грж) ¿ж,

о т

Ш^ГУ/ гИехр^ (ж,р) ¿ж,

о

т

14(р) = (-2г)^+1р^У г(ж) ехр(грж)п(ж, р) ¿ж,

о

Так как

сю

/ж°

— ехр(2грж) ¿ж = 1, р е ф, к!

о

сю

/ж°

— ехр(2грж) ¿ж = о(1), р е ф, |р| ^ го, к!

т

то 11(р) - 7 ^ 0 при |р| ^ го, р е ф. Если р е ф, то существует го > 0 такое, что

11т р| > го|р| при р е ф. (12)

Возьмем г > 0 и выберем 6 = 6(г) так, чтобы |р(ж)| < 2£о+1 при ж е [0,6], где го определено в (12). Тогда, используя (12), вычисляем

£ т

г жо жо

(р)| < 2(2|р|го)°+Ч к, ехр(-2го|р|ж) ¿ж + (2|р|)°+М —|р(ж)| ехр(-2го|р|ж) ¿ж <

0 £ Г—£

< 2 + (2|р|)°+1 ехр(-2го|р|6^ |р(ж + 6)| ехр(-2го|р|ж) ¿ж.

о

В силу произвольности г получаем 12(р) ^ 0 при |р| ^ го, р е ф.

Так как |(7 + р(ж))£(ж, р) | < С, то при |р| ^ го, р е ф имеем

Гт

/жо С

— ехр(-2го|р|ж) ¿ж ,

к! Нго

о

о

следовательно, 13(р) ^ 0 при |р| ^ р G Q. Аналогично получаем 14(р) ^ 0 при |р| ^ го, р G Q. □

Для простоты предположим, что функция q(x) является аналитической на [0,Т]. Пусть при некотором фиксированном k ^ 0 коэффициенты Тейлора qj := q(j)(0), j = 0, k — 1 уже найдены. Выберем модельный потенциал q(x) так, чтобы первые k коэффициентов Тейлора функций q и q совпадали, т.е. qj = qj, j = 0, k — 1. Тогда, используя (10), (11) и лемму 2, можно вычислить следующий коэффициент Тейлора qk = q(k)(0). Точнее, справедливо следующее утверждение.

Лемма 3. Зафиксируем k. Пусть функции q(x) и q(x) являются аналитическими при x G [0, T] с (¡j : = qj — qj = 0 для j = 0, k — 1. Тогда

qk = (—2)k+1 lim (¿р)к+3 N(A). (13)

peQ

Таким образом, мы приходим к следующему алгоритму решения обратной задачи 1.

Алгоритм. Пусть задана функция Вейля N(A).

1. Вычисляем qk = q(k) (0), k ^ 0. Для этого последовательно выполняем следующие операции при k = 0,1,2,... : выбираем модельный потенциал q(x) так, чтобы qj = qj, j = 0, k — 1, и вычисляем qk = q(k) (0) by (13).

2. Строим функцию q(x) по формуле

те k

q(x) = qkk!, 0 < x < R, k=0 '

где

|qk К1/k\ -1

R = lim . ,,

Хк^те V k!

Если R < п, то при R < x <п функция q (x) строится по аналитическому продолжению.

Благодарности. Работа выполнена при финансовой поддержке РНФ (проект № 17-11-01193).

Библиографический список

1. Марченко В. А. Операторы Штурма - Лиувилля и их приложения. Киев : Наук. думка, 1977. 393 с.

2. Левитан Б. М. Обратные задачи Штурма - Лиувилля. М. : Наука, 1984. 246 с.

3. Freiling G., Yurko V. A. Inverse Sturm - Liouville Problems and their Applications. N. Y. : NOVA Sci. Publ., 2001. 305 p.

4. Beals R., Deift P., Tomei C. Direct and Inverse Scattering on the Line. (Math. Surveys and Monographs. Vol. 28). Providence, RI : Amer. Math. Soc., 1988. 209 p.

5. Yurko V. A. Method of Spectral Mappings in the Inverse Problem Theory. Inverse and Ill-posed Problems Series. Utrecht : VSP, 2002. 316 p.

6. Yurko V. A. Inverse Spectral Problems for Differential Operators and their Applications. Amsterdam : Gordon and Breach, 2000. 265 p.

7. Lakshmikantham V., Rama Mohana Rao M. Theory of Integro-differential Equations. (Stability and Control : Theory and Applications. Vol. 1). Singapure : Gordon and Breach, 1995. 308 p.

8. Юрко В. А. Обратная задача для интегро-дифференциальных операторов // Матем. заметки. 1991. Т. 50, вып. 5. С. 134-146.

9. Kuryshova Yu. An inverse spectral problem for differential operators with integral delay // Tamkang J. Math. 2011. Vol. 42, № 3. P. 295-303.

10. Buterin S. On an inverse spectral problem for a convolution integro-differential operator // Results in Math. 2007. Vol. 50. P. 173-181.

11. Бутерин С. А. О восстановлении сверточного возмущения оператора Штурма -Лиувилля по спектру // Дифференц. уравнения. 2010. Т. 46, № 1. С. 146-149.

Образец для цитирования:

Игнатьев М. Ю., Советникова С. Ю. О восстановлении интегродифференциальных операторов по функции Вейля // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2017. Т. 17, вып. 3. С. 276-284. ЭО!: 10.18500/1816-9791-2017-17-3-276-284.

On Recovering Integro-Differential Operators from the Weyl Function M. Yu. Ignatiev1, S. Yu. Sovetnikova2

1 Mikhail Yu. Ignatiev, ORCID: 0000-0002-4354-9197, Saratov State University, Astrakhanskaya Str., 83, Saratov, Russia, 410012, IgnatievMU@info.sgu.ru

2Svetlana Yu. Sovetnikova, ORCID: 0000-0002-5703-8195, Saratov State University, Astrakhanskaya Str., 83, Saratov, Russia, 410012, sovetnikovasy@mail.ru

We study inverse problems of spectral analysis for second order integro-differential operators, which are a perturbation of the Sturm - Liouville operator by the integral Volterra operator. We pay the main attention to the nonlinear inverse problem of recovering the potential from the given Weyl function provided that the kernel of the integral operator is known a priori. We obtain properties of the spectral characteristics and the Weyl function, provide an algorithm for constructing the solution of the inverse problem and establish the uniqueness of the solution. For solving the inverse problem we use the method of standard models.

Key words: integro-differential operators, inverse spectral problems, uniqueness result, algorithm.

Acknowledgements: This work was supported by the Russian Science Foundation (project no. 17-11-01193).

References

1. Marchenko V. A. Sturm - Liouville operators and their applications. Birkhauser, 1986. 393 p. (Russ. ed. : Kiev, Naukova Dumka, 1977. 393 p.)

2. Levitan B. M. Inverse Sturm - Liouville problems. Utrecht, VNU Sci. Press, 1987. 246 p. (Russ. ed. : Moscow, Nauka, 1984. 246 p.)

3. Freiling G., Yurko V. A. Inverse Sturm - Liouville Problems and their Applications. New York, NOVA Sci. Publ., 2001. 305 p.

4. Beals R., Deift P., Tomei C. Direct and Inverse Scattering on the Line. (Math. Surveys and Monographs, vol. 28), Providence, RI, Amer. Math. Soc., 1988. 209 p.

5. Yurko V. A. Method of Spectral Mappings in the Inverse Problem Theory. Inverse and Ill-posed Problems Series. Utrecht, VSP, 2002. 316 p.

6. Yurko V. A. Inverse Spectral Problems for Differential Operators and their Applications. Amsterdam, Gordon and Breach, 2000. 265 p.

7. Lakshmikantham V., Rama Mohana Rao M. Theory of integro-differential equations. (Stability and Control: Theory and Applications, vol.1), Singapure, Gordon and Breach, 1995. 308 p.

8. Yurko V. A. An inverse problem for integro-differential operators. Math. Notes, 1991, vol. 50, iss. 5-6, pp. 1188-1197.

9. Kuryshova Yu. An inverse spectral problem for differential operators with integral delay. TamkangJ. Math., 2011, vol. 42, no. 3, pp. 295-303.

10. Buterin S. A. On an inverse spectral problem for a convolution integro-differential operator. Results in Math., 2007, vol. 50, pp. 173-181.

11. Buterin S. A. On the reconstruction of a convolution perturbation of the Sturm - Liouville operator from the spectrum. Differential Equations, 2010, vol. 46, no. 1, pp. 150-154.

Cite this article as:

Ignatiev M. Yu., Sovetnikova S. Yu. On Recovering Integro-differential Operators From the Weyl Function. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2017, vol. 17, iss. 3, pp. 276-284 (in Russian). DOI: 10.18500/1816-9791-2017-17-3-276-284.