УДК 517.984
ОБРАТНАЯ ЗАДАЧА ДЛЯ ОПЕРАТОРОВ ШТУРМА-ЛИУВИЛЛЯ В КОМПЛЕКСНОЙ ПЛОСКОСТИ
А. А. Голубков
Голубков Андрей Александрович, доктор физико-математических наук, профессор кафедры физики, Специализированный учебно-научный центр (факультет) — школа-интернат имени А. Н. Колмогорова Московского государственного университета имени М. В. Ломоносова (СУНЦ МГУ), Россия, 121357, Москва, Кременчугская, 11, апс1ге]2501 @уапс1ех.ш
Впервые изучена обратная задача для стандартного уравнения Штурма-Лиувилля со спектральным параметром р и потенциалом, кусочно-целым на спрямляемой кривой 7 с С, у которой задана только начальная точка. Ограниченная на кривой 7 функция Q является кусочно-целой на ней, если 7 можно разбить конечным числом точек на участки, на которых Q совпадает с целыми функциями, различными на соседних участках. Точки разбиения, начальная и конечная точки кривой называются критическими точками. Ставится задача нахождения всех критических точек 7 и потенциала на ней по столбцу или строке передаточной матрицы Р вдоль 7. На основе полученной асимптотики Р при |р| ^ го доказано, что если хотя бы один её элемент ограничен при любых р € С, то 7 после удаления всех «невидимых петель» вырождается в точку («невидимая петля» -- такая петля кривой 7 с заданной кусочно-целой функцией, узел которой совпадает с двумя последовательными критическими точками). В статье доказана единственность решения поставленной обратной задачи для кривых без «невидимых петель». На примере обратной задачи для уравнения а! (яж) ах) + (я(х) - г(х)А2) у(х) = 0 с кусочно-целым потенциалом д(х) и кусочно-постоянной функцией г(х) = 0 на отрезке действительной оси показана полезность полученных результатов при исследовании обратных задач для обобщенных уравнений Штурма-Лиувилля, приводимых к изученному в статье типу.
Ключевые слова: уравнение Штурма-Лиувилля на кривой, кусочно-целый потенциал, передаточная матрица, асимптотика, обратная спектральная задача.
001:10.18500/1816-9791-2018-18-2-144-156
1. ПОСТАНОВКА ЗАДАЧИ И ОСНОВНЫЕ РЕЗУЛЬТАТЫ
Обратные задачи для уравнения Штурма - Лиувилля стандартного вида
и"(*) + (^) - А2)«(2) = 0 (1)
в случае вещественной переменной хорошо изучены для различных, в том числе комплекснозначных потенциалов [1-3]. Однако обратные спектральные задачи для обобщенных уравнений Штурма - Лиувилля
5Х (/(х)Ю + (<?(х) - Г(х)А2) и(х) = 0 (х € [0,1])
с кусочно-аналитическими комплекснозначными коэффициентами на отрезке в общем случае почти не исследованы [4]. Между тем такого типа обратные задачи имеют большое практическое значение, так как возникают, например, в спектроскопии одномерно неоднородных сред [5,6]. Одним из эффективных методов их изучения могло бы стать преобразование с помощью соответствующих подстановок [7]
обобщенных уравнений на отрезке в стандартные уравнения Штурма - Лиувилля на кривых в комплексной плоскости. Например, в случае f (х) = 1,д(х) = 0 полезность такого перехода продемонстрирована в статье [8]. К сожалению, несмотря на большое число работ, посвященных исследованию решений обыкновенных дифференциальных уравнений в комплексной плоскости (см. книги [9-11] и библиографию в них), асимптотическое поведение решений уравнений Штурма - Лиувилля (1) вдоль кривых при больших значениях модуля спектрального параметра хорошо изучено только при достаточно жестких ограничениях на форму кривой и (или) на потенциал [8, 11, 12]. Вероятно, именно поэтому среди обратных задач на комплексной плоскости детально исследована только задача о безмонодромных уравнениях Штурма - Лиувилля с мероморфным потенциалом на кусочно-гладкой кривой, являющейся границей некоторой выпуклой ограниченной области [13,14]. В данной работе впервые рассмотрена обратная спектральная задача для уравнения (1) с потенциалом кусочно-целым на спрямляемой кривой 7 с С произвольной формы.
Определение 1. Функцию Q будем называть кусочно-целой на кривой 7 с С, соединяющей точки ,zf и допускающей параметрическое задание непрерывной функцией z = z(£) (£ е [t0,tf(£0) = z0,z(tf) = Zf), если Q ограничена на этой кривой и совпадает на ней с конечным числом целых функций, т. е. существуют такие целое число N ^ 0 и набор чисел Т = ^^: ^ < ^ < ... < = tf, что
Q(z(t))= Qi(z), если t е (^,^+1) (г = 0^), (2)
где все Qi — целые функции и при N ^ 1 для всех п е {1,..., N} функции Qn и Qn-1 различны.
Замечание 1. Кусочно-целая на кривой функция может быть многозначной: принимать разные значения при прохождении кривой дважды через одну точку.
В силу (2) функция Q и все её производные вдоль кривой 7 имеют конечные односторонние пределы в точках
:= х) (о = 0ГЯГТ).
Если N ^ 1, то для всех п е {1,..., N} односторонние пределы функции Q или хотя бы одной её производной конечного порядка вдоль 7 различны при подходе к zn «слева» и «справа», так как по определению 1 функции Qn и Qn-1 разные. Будем называть такие точки точками обобщенного скачка функции Q на кривой 7. Точки обобщенного скачка, а также начальную и конечную точки кривой будем называть её критическими точками.
Определение 2. Назовём Ш := {^ ^{Qi} набором ключевых данных кривой 7 и заданной на ней с помощью (2) кусочно-целой функции Q.
Определение 3. Пусть и1 (7^),и2(7^) — непрерывно дифференцируемые решения (1) вдоль кривой 7, проходящей через точку zb, удовлетворяющие условиям:
и(7, Zb) = 1, и1 (7, Zb) = 0, и(7, Zb) = 0, и2(7, Zb) = 1. (3)
Назовём передаточной матрицей уравнения (1) матрицу
Пт,^) = ) ^0.
Лемма 1. Если функция Q является кусочно-целой на спрямляемой кривой 7, соединяющей точки и то передаточная матрица Р(ч, zf ) уравнения (1) существует и однозначно определяется заданием р = А2 и набором ключевых данных Ш. Матрица Р(7, zf, 20) является целой функцией параметра р; её определитель равен 1, и для всех 2 € 7 справедливы следующие соотношения:
Р(7,2,20 )= Р(7,^с )Р(7,2С ,20 ), (4)
Р(7,*0)= Р-1 (7,2,20)=( и2^ . (5)
Доказательство. Рассмотрим вспомогательные уравнения Штурма - Лиувилля:
Ж2 и«
Ж
2
+ ^ - А2)и^ =0, 2 € С, i = 0, N. (6)
В силу (2) и определения 3 решения и1 (7,2),и2(7,2) могут быть построены «сшиванием» непрерывно дифференцируемых решений уравнений (6) вдоль соседних участков кривой 7. Поэтому лемма 1 следует из известных свойств решений
линейных дифференциальных уравнений с голоморфными коэффициентами [10,11]. □
Следствие 1. Для расчета передаточной матрицы вместо исходной кривой 7 при решении (1) можно использовать любую спрямляемую кривую 7 (например, ломаную), последовательно соединяющую критические точки кривой 7, если для каждого i € {0,..., N} на участке 7, соединяющем точки и , полагать
Q(z) = Qi (2).
Определение 4. Петлёй кривой 7 с узлом в точке начинающийся и кончающийся в точке её самопересечения 2. Все точки петли, кроме её узла, будем называть внутренними.
Определение 5. Пусть на кривой 7 задана кусочно-целая функция Q. Петля кривой 7 называется «невидимой петлей» уровня один или просто «невидимой петлей», если её узел совпадает с двумя последовательными критическими точками. Петля кривой 7 называется «невидимой петлей» уровня М ^ 2, если её узел совпадает с двумя критическими точками и функция Q совпадает с одной целой функцией на всей петле, кроме конечного числа участков, как минимум один из которых является «невидимой петлёй» уровня М - 1, а остальные — «невидимыми петлями» меньшего уровня.
Лемма 2. Добавление к кривой или удаление из неё «невидимой петли» любого конечного уровня не меняет начальную и конечную точку кривой, а также передаточную матрицу вдоль неё.
Доказательство. В случае «невидимой петли» уровня один лемма следует из определения 5, формулы (4) и того, что передаточная матрица уравнения (1) с целым потенциалом вдоль любой петли равна единичной матрице [10,11]. Применяя метод индукции по уровню «невидимой петли» и пользуясь определением 5, получим, что лемма справедлива для «невидимой петли» любого конечного уровня. □
Определение 6. Кривую с заданной на ней кусочно-целой функцией будем называть простой, если на кривой отсутствуют «невидимые петли» любого уровня.
Определение 7. Набор ключевых данных Ш будем называть простым, если для него выполнены следующие условия:
Azi := Zi+i - Zi = 0 (i = 0,N). (7)
Упорядоченное множество точек {zj+1, входящих в простой набор ключевых данных, будем называть множеством базовых точек соответствующего потенциала Q.
Лемма 3. Кривая y с заданной на ней кусочно-целой функцией Q является простой, если и только если ей соответствует простой набор ключевых данных.
Доказательство. Лемма непосредственно следует из определений 5, 6 и 7, так как в силу определения 5 на кривой могут быть «невидимые петли» уровня M ^ 2, только если на ней есть «невидимые петли» предыдущего уровня, а выполнение всех условий (7) равносильно отсутствию «невидимых петель» первого уровня. □
В силу следствия 1 и леммы 2 при заданной точке z0 передаточная матрица P уравнения (1) может соответствовать бесконечно большому числу кривых с кусочно-целыми потенциалами, отличающихся в том числе количеством, расположением и формой «невидимых петель» и, значит, наборами ключевых данных. После удаления «невидимых петель» всех уровней исходная кривая либо переходит в простую кривую, либо вырождается в точку. В последнем случае очевидно, что P = I, где I — единичная матрица.
Теорема 1. Если хотя бы один элемент передаточной матрицы P уравнения (1) с кусочно-целым потенциалом на некоторой (не заданной) спрямляемой кривой y ограничен на всей комплексной плоскости параметра р, то y после удаления «невидимых петель» всех уровней вырождается в точку, и P = I.
Доказательство теоремы 1 приведено в парагр. З.Из теоремы 1 следует следующая теорема единственности, доказательство которой приведено в парагр. 4.
Теорема 2. Пусть z0^ = z02) = z0. Тогда простым спрямляемым кривым y(1) и y(2) с разными простыми наборами ключевых данных W(1) = |Njzj^ j , f rnЛм(1) 1 í f (2)(2)+i f í2)in(2) 1
|Q¿ | \ и W(2) = < N(2), |zj j , |Q« j г соответствуют передаточные матрицы P(1) и P(2), у которых соответствующие строки и столбцы различны хотя бы при одном значении параметра р.
Следствие 2. Пусть функция Q — кусочно-целая на спрямляемой кривой y, у которой задана только начальная точка z0. Тогда задание при всех р е C столбца или строки отличной от единичной передаточной матрицы уравнения (1) между точкой z0 и конечной точкой кривой y однозначно определяет набор ключевых данных кривой, получающейся из y удалением всех «невидимых петель».
Замечание 2. В силу леммы 1 все элементы передаточной матрицы являются целыми функциями р. Поэтому в следствии 2 и всех последующих утверждениях достаточно требовать их задания (совпадения) на множестве точек комплексной плоскости параметра р, имеющем хотя бы одну конечную предельную точку.
Следствие 3. Если два уравнения (1) имеют вдоль данной спрямляемой кривой 7 при всех р € С передаточные матрицы с одинаковым столбцом (строкой), то их потенциалы совпадают всюду на 7, за исключением, возможно, тех её однозначно определённых петель, среди внутренних точек которых нет базовых точек потенциала.
Доказательство. В силу определения 5 и лемм 1, 2 замена потенциала Q на петле, среди внутренних точек которой нет его базовых точек, на любую целую функцию превращает эту петлю в «невидимую», не меняя передаточную матрицу. Поэтому на таких петлях потенциалы в двух рассматриваемых уравнениях могут различаться (а могут и совпадать). Если после удаления всех таких петель кривая 7 вырождается в точку, то следствие доказано. Если же она превращается в некоторую кривую 7тт, то последняя будет заведомо простой, так как по построению не может иметь «невидимых петель» первого уровня и, значит, любого другого. Поэтому по следствию 2 потенциал Q однозначно определен на 7ть. □
Теорема 2 может быть полезной при изучении обратных задач для различных типов уравнений Штурма - Лиувилля, приводимых к уравнению (1) с кусочно-целым потенциалом. В качестве примера рассмотрим обратную спектральную задачу для уравнения
с кусочно-целым потенциалом д(х) и кусочно-постоянной отличной от нуля функцией г(х) на отрезке [х0,xf] действительной оси. Все коэффициенты в уравнении (8) комплекснозначные. Причем нас будут интересовать такие непрерывные решения (8), для которых функция у'(х)/г(х) также непрерывна. В качестве спектральных данных обратной задачи для уравнения (8) используем его обобщенную передаточную матрицу между концами отрезка [х0, xf]:
где у1 (х), у2(х) — решения (8) на отрезке [х0, xf] в указанном выше смысле, удовлетворяющие следующим условиям:
У1(х0) = 1, у1 (х0 )/г(х0) = 0, у2 (х0) = 0, у2(х0 )/г(х0) = 1.
Нетрудно убедиться, что если уравнение (8) имеет на отрезке [х0, xf] обобщенную передаточную матрицу Рг, то замена и(2) := у(х),
приводит его к уравнению (1) с кусочно-целым потенциалом, однозначно определенным на задаваемой (9) ломаной Ь, причём передаточная матрица этого уравнения (1) вдоль Ь равна Рг. Пользуясь этим, из теоремы 2 можно получить два следствия.
Следствие 4. Если два уравнения (8) с заданной кусочно-постоянной ком-плекснозначной функцией г(х) имеют на отрезке [х0, xf] обобщенные передаточные матрицы с одинаковым столбцом или одинаковой строкой, то их кусочно-целые потенциалы совпадают всюду, за исключением, возможно, однозначно
(9)
определенных строкой или столбцом передаточной матрицы «невидимых участков» отрезка [х0, х/], переходящих при отображении (9) в «невидимые петли» ломаной Ь.
Следствие 5. Пусть на некотором отрезке действительной оси хотя бы одно уравнение типа (8) имеет отличную от единичной обобщенную передаточную матрицу с заданным столбцом или заданной строкой. Тогда этим столбцом или строкой однозначно определена такая величина Ь(0) > 0, что среди бесконечного множества уравнений типа (8), имеющих на отрезке [0,1] комплексно-значный коэффициент г(х) с постоянным модулем и передаточную матрицу с заданным столбцом или заданной строкой, 1) существует единственное уравнение, у которого |г(х)| = Ь(0); 2) не существует уравнения, у которого |г(х)| = Ь(г) < Ь(0); 3) существует бесконечно много уравнений, у которых |г(х)| = Ь(г) > Ь(0).
Заметим, что следствие 4 существенно ослабляет полученные в статье [4] достаточные условия однозначности нахождения потенциала q(х) уравнения (8) по известному столбцу или известной строке обобщенной передаточной матрицы и известной функции г(х). В соответствии с результатами статьи [4] для этого достаточно, чтобы все значения функции г(х) лежали в открытой половине комплексной плоскости с границей, проходящей через ноль. По следствию 4 для этого достаточно, чтобы задаваемая (9) ломаная Ь не имела самопересечений. Правда, в [4] рассматривался более общий тип обобщенного уравнения Штурма - Лиувилля с тремя независимыми кусочно-аналитическими комплекснозначными коэффициентами, любые два из которых считались известными.
2. АСИМПТОТИКА ПЕРЕДАТОЧНОЙ МАТРИЦЫ
Лемма 4. Для любой ограниченной области С с С, любых К е N и г е {0,..., N} и любого набора чисел {д,,,}К существуют числа А,,К > 0, С^К > 0
и два непрерывно дифференцируемых решения Р+К(я, А) соответствующего уравнения (6), которые для всех А = 0 и я е С могут быть представлены в виде
(я, А) = (я, А) ехр{±А(я - я,)}, —^ = ±А£&(я, А) ехр{±А(я - я,)}, (10)
(я, А) := 1 + ±(±1)к^ + (±1)К+1 ^Акг^
Ак у Ак+1
к=1
(я, А) := 1 + ±(±1)кС*^(я) + (±1)к+1 ^АКР
к=1
С(г) := -(Сл-Л*) + Я<(я)С,,-1 (г))А С,,о(я) := 1, С,,(2,)= , (£ = Т7К).
(11) (12) (13)
Причём если |А| ^ А,то (я, А) — линейно независимые решения (6) в С и для всех я е Ь, (Ь, — отрезок, соединяющий точки я, и я,+1):
(я, А)| < С(,0К, (я, А)| < С(,0К, (14)
(я, А) - 1| < 1/4, (я, А) - 1| < 1/4. (15)
Доказательство. Формулы (10)-(13) проверяются подстановкой в (6), а оценки (14), (15) на отрезке Ь следуют из известных результатов по асимптотическому разложению при |А| ^ то решений обобщенных уравнений Штурма - Лиувилля на отрезке действительной оси [10]. □
Символами 0(1) и 0(1) будем обозначать соответственно функции и матрицы функций параметра А, конкретный вид которых для нас не важен, ограниченные при |А| > Асг, где Асг - конечная величина (разная для разных функций и матриц).
Лемма 5. Пусть функция Q является кусочно-целой на спрямляемой кривой 7, соединяющей точки 20 и zf. Тогда существует такое целое число К0 ^ 2, что для любого целого К ^ К0 существует такое конечное Ак > 0, что при |А| ^ Ак матрица , 20) может быть записана в виде
Р(7, zf, z0) = С^)Т)Т-1)... Т(2)Т(1)Т(0)4(0) . (16)
В (16) Т(0) = I (единичная матрица),
4(0) = ( АЕ-0К(«,, А) С-К («,, А)
TT(n) = ff ( ^УК^-и u j (n = T7Ñ При N > 1), (19)
0 V 0 exp(-AAzn-i) ) V У h
D? V AE+K(zo,A) —C+K(zo,A)
((f) = / C+K(zn+1, A)exp(AAzN) C-K(z^+i, A) exp(—AAzn) AE+K(zn+1, A) exp(AAzN) -AE-K)(zn+1, A) exp(—AAz^)
(n) ( exp(AAzn-i) 0
0
где Azi (i = 0, N) определены в (7). Л (17) и (19) соответственно
Df = -A(C+0¿(zo, A)E-0K(zo, A) + C-K (zo, A)E+oK(zo, A)) = 0,
f (n) = / 1 + C-K-1' (zn, A)W+K (A) C-K-1)(zn, A)W-nK (A) \ 0(1) o l -С+П-1) (zn, A)W+n)(A) 1 - С+П-1) (z„, A) W-n)(A) j + AK+1'
причём при |A| ^ AK функции W±k(A) можно представить в виде
W±n) (A) = ^(t^A) "+¿n(l + °AT^)=0 (n = T,N при N ^ 1) , (20)
где целые числа mn (mn G [0, Ko — 2]) и комплексные числа 5n = 0 не зависят от A.
Доказательство. По лемме 4 для всех K G N существуют такие числа Ai)K > 0 (i = 0,N), что при |A| ^ AijK решения u1 (7, z),u2(7, z) уравнения (1) на участке кривой y, соединяющем точки zi и zi+1, можно представить как линейную комбинацию решений F±K(z, A) c
{Co,fc(zo)}K = 0 (при i = 0), Ci,fc(zi) = Ci_i,fc(zi) (при N ^ 1, i = T,N, k = T7K).
Поэтому при |А| ^ Ак := шах{А^,к, I = 0, N}, пользуясь граничными условиями (3) и непрерывной дифференцируемостью и1 (7,2), и2 (7,2) в точках обобщенного
скачка потенциала, получим, что для Р(7, я/,я0) справедливо представление (16), в котором матрицы А(0), С(/) и Т(" (п = 1,N при N ^ 1) удовлетворяют соотношениям (17)-(19). При этом С^К и (п = 0, N) в силу (15) не равны нулю в точках и
^п+Ъ а
(А) = -¿У £(±1)*^) "Г"-1*) (п = при N > 1), (22)
где £>(" = -А(С+"К(я", А)Е-"К^, А) + С-"К(я", А)Е+"К(^, А)). Из (11)-(14) имеем:
= -2А ( 1 + =0 (г = 0^). (23)
Остаётся доказать справедливость (20). При N ^ 1 для всех п е {1,...,N} по определению 1 функции и ф" различны. Значит, существуют целое число
Ш" ^ 0 и комплексное число 5" = 0 такие, что при я =
dzm dzm dmQn dm Qn-i
- = ¿n = 0 при m = m„; (24)
= 0 при m = 0, mn — 1, если mn ^ 1. (25)
dzm dzm
Кроме того, для всех n G {1,..., N} в силу (13), (21), (24) и (25) /dm(Cn,i(z) — Cn-1,1 (z))\ io, 0 < m < mn,
(26)
V dzm ) z=zn y— ¿n/2, m = mn + 1.
При mn = 0 (20) непосредственно следует из (22), (23) и (26). Пользуясь (13), (21), (25), (26) и методом индукции по k, получим, что для mn ^ 1 и k = 1,mn + 1
idm (Cn,fc(z) — Cn-i,k(z)) \ = i0, 0 < m < mn + 1 — k,
V dzm !(—-)k ¿n, m = mn + 2 — k.
Полагая в последнем равенстве m = 1 и учитывая (22), (23), получим (20) при mn ^ 1 и K ^ mn + 2. Полагая K0 := 2 + max{mn,n = 1,N}, получаем, что (20) выполнено для любого целого K ^ K0 и |А| ^ Ак. □
Будем обозначать (N + 1)-мерные вектора буквой со стрелкой сверху, а их скалярное произведение — круглыми скобками. Например, z := (Az0, Az1,..., AzN); а := (a0s),ais), • • •, «N), где a(s) G {±1} (i = 0,N), а s = 1 + £]=0(1 + a(s))2i-i, т.е. s G {1,..., 2n+1} (как в двоичной системе счисления); (z, а.) = ^¿=0 a(s)Az^.
Следствие 6. В условиях леммы 5 при K ^ K0 и |А| ^ AK элементы матрицы p(y, z/,z0) могут быть записаны в виде
2n+1
Paß = £ dsß(А) exp {Ah.} (a, в G {1,2}), (27)
s=1
где h. = (z, a.) не зависят от А и все функции d0ß(А) можно представить в виде
(s)
daß (А) = (А) ¿Oß (1 + Op) = 0 (a, в G {1,2}, s = и**1),
где целые числа m0ß ^ 0 и комплексные числа ¿0ß = 0 не зависят от А.
Часть показателей экспонент в (27) может совпадать, и возникает вопрос: не сокращаются ли в (27) коэффициенты при экспонентах, наиболее быстро растущих с ростом |А|? Для простых кривых отрицательный ответ даёт следующая лемма.
Лемма 6. Пусть ^тах := шах{|^|, 5 = 1,2^+1}, где ^ := (2,а). Тогда существуют хотя бы два различных в0 е {1,..., 2^+1} таких, что 1 = ^тах, причём если 2 = 0, то ^тах > 0. Кроме того, при выполнении (7) = ^ для любого такого, что |^01 = ^тах и для всех в е {1,..., 2^+1 }\{в0}.
Доказательство. Существование хотя бы одного в0 е {1,..., 2^+1} : |^01 = ^тах следует из конечности множества {^, в = 1,2^+1}. Так как множество , в = 1, 2^+1} можно разбить на пары векторов, сумма которых равна нулю, то число таких в0, что |^01 = ^тах, всегда чётное. Допустим, что ^тах = 0. Тогда все ^ = 0, и компоненты вектора 2 удовлетворяют, в частности, системе линейных однородных уравнений вида: в уравнении под номером п (п = коэффициент при равен минус 1, а остальные коэффициенты равны 1. В последнем N + 1 уравнении все коэффициенты равны 1. Такая система имеет только тривиальное решение, поэтому ^тах = 0, только если 2 = 0.
Докажем вторую часть леммы. Пусть существует в1 е {1,..., 2^+1 }\{в0}: (г, (ав1 — а50)) = 0. В силу (7) = 0 для всех г = 0, N, поэтому а50 и ав1 имеют минимум две не совпадающие (противоположные по знаку) компоненты с номерами 0 ^ г1 < ... < г/ ^ N. Положим е0 := (ав1 + а50)/2. Тогда ^тах = |(2,е0)| > 0 (в силу (7) и первой части леммы) и у е0 не равны нулю (равны ±1) только компоненты с номерами, отличными от г1,г2,...,г/. Значит, 2 ^ J ^ N, что возможно, только если N ^ 2.
Рассмотрим все вектора вида вР = (в0,..., в^), где в« е {±1} при г е {г1,..., }, а остальные компоненты равны нулю (р = 1, ). Пусть для всех р е {1,..., } |(г,вР)| = 0. Тогда Дг«1,Дг«2,...,удовлетворяют, в частности, J линейным однородным уравнениям следующего вида: в уравнении под номером ] = 1, J — 1) коэффициент при Дг^. равен минус 1, а остальные коэффициенты равны 1. В последнем уравнении все коэффициенты равны 1. Очевидно, такая система имеет только тривиальное решение, что противоречит (7). Значит, среди векторов вР (р = 1, ) есть хотя бы один вР0: вР0)| > 0. Рассмотрим два вектора: т/± = е0±вР0. По построению они оба входят в множество , в = 1, 2^+1}. Причём (2, п±) = (2, е0) ± вР0) и |(^,вро)1 > 0. Поэтому шах{|(^, )|, П-)|} > 1(2, ¿0)| = ^тах (при геометрическом представлении комплексных чисел использованное неравенство эквивалентно тому, что медиана всегда меньше большей из двух сторон треугольника, между которыми она лежит). Полученное противоречие доказывает вторую часть леммы 6. □
Лемма 7. Пусть у передаточных матриц Р и Р на всей комплексной плоскости параметра р совпадают первые или вторые строки. Тогда хотя бы один из
элементов матрицы Р(0) := Р(Р)-1 равен нулю на всей комплексной плоскости р.
Доказательство. Перемножая матрицы, получим:
Р (0) = (Р22р11 — Р21р12 Р11 р12 — Р12ри\ \Р22Р21 — Р21Р22 Р11Р22 — Р12Р21) '
Откуда следует, что если у матриц Р и Р совпадают элементы первой строки, то р102 = 0, а если второй, то р21) =0. □
3. ДОКАЗАТЕЛЬСТВО ТЕОРЕМЫ 1
Пусть в условиях теоремы 1 после исключения всех «невидимых петель» кривая 7 превращается в простую кривую 7тщ. Тогда Р можно рассчитать вдоль 7т1п. Из леммы 6 следует, что в этом случае все показатели экспонент в правых частях равенств (27) с максимальными (при данном А = 0) модулями различны и не равны нулю, и существует минимум два противоположно направленных луча, исходящих из нуля комплексной плоскости параметра А (и, значит, минимум один луч, исходящий из нуля комплексной плоскости параметра р), такие, что среди 2]+1 слагаемых в (27) для каждого элемента Р существует ровно одно слагаемое, имеющее наибольший экспоненциальный рост при А ^ то (р ^ то) вдоль каждого из этих лучей, и значит, каждый из элементов матрицы Р при этом неограниченно возрастает. Последнее противоречит условию теоремы 1. Значит, после исключения всех «невидимых петель» кривая 7 вырождается в точку, и, следовательно, Р = I.
4. ДОКАЗАТЕЛЬСТВО ТЕОРЕМЫ 2
Докажем теорему 2 методом от противного, рассмотрев два случая.
Случай 1. Пусть у передаточных матриц Р(1) и Р(2) при всех р £ C совпадает одна из строк. Из (4), (5) следует, что передаточная матрица уравнения (1) вдоль кривой 7, получающейся при последовательном обходе сначала кривой 7(1) (с заданным на ней кусочно-целым потенциалом) от точки к точке 20, а затем кривой 7(2) от
точки 20 к точке +1, будет равна Р = Р(2) (Р(1))-1. Тогда в силу леммы 7 и теоремы 1 кривая 7 после удаления «невидимых петель» всех уровней вырождается в точку. По условию ¿01( = ¿02( = 20 и простые наборы ключевых данных Ш(1) и Ш(2) различны. Это означает, что существует целое число г0 ^ 1 такое, что для всех г = 0, г0 — 1 2(2) = 2(1) = 2 и при г0 ^ 2 для всех г = 0, г0 — 2 0(2) = 0(1), но при этом 2(02) = 2(1) или (и) = ^(01-1. Пусть для опредёленности N(1) ^ N(2).
Допустим, что г0 = N(1) + 1. Тогда после удаления всех невидимых петель у кривой 7 останется только участок простой кривой 7(2), соединяющий точки +1
и +1. Поскольку 7 должна вырождаться в точку, то N(1) = N(2), т.е. наборы ключевых данных Ш(1) и Ш(2) совпадают, что противоречит условию теоремы.
Пусть г0 < N(1) + 1, а функции О^1— и различны. Тогда после удаления у
кривой 7 «невидимых петель», образованных участками кривых 7(1) и 7(2) от точки 20 до точки -1, оставшаяся кривая 7т1п сначала будет совпадать с участком кривой 7(1) от точки г:]1!)+1 до точки -1, а затем с участком кривой 7(2) от точки -1 (2)
до точки (2)+1. Причём точка -1 будет точкой обобщенного скачка потенциала 0 на кривой 7т1п и, значит (так как по условию наборы ключевых данных Ш(1) и Ш(2) простые), кривой 7т1п будет соответствовать простой набор ключевых данных, что противоречит тому, что она должна вырождаться в точку. Следовательно,
С-! = 02-1 • (28)
Допустим, что г0 < N (1) + 1, выполнено (28), но 2(2) = 2(01). Тогда опять же кривой 7 будет соответствовать простой набор ключевых данных. Следовательно,
4Ч = • (29)
Формулы (28), (29) противоречат определению числа i0, что доказывает теорему в случае 1.
Случай 2. Обозначим := z^L)+1 (а £ {1,2}). Пусть у передаточных матриц Pи P(2) при всех р £ C совпадают элементы одного из столбцов. Тогда в силу (5) у передаточных матрицы (y(1),z0, zf1) и p^Y(2), z0,zf2)) совпадают элементы
одной из строк. Пусть т := z(1) — z(2). Сделав в уравнении (1) с потенциалом Q(2) (z) замену переменной: x = z + т, после обратного переобозначения x на z получим, что матрица P^i является передаточной для уравнения (1) с потенциалом Q(2) (z — т) между точками z(1) и z0 + т простой кривой, получающейся из 7(2) параллельным переносом. Пользуясь теоремой 2, доказанной выше для случая 1, получаем:
N <2> = N (1> = N, Qf(z — т) = Q« (z), zf + т = zf (i = ÖTN), z<2»1 + т = z<J+1. ( '
По условию теоремы z02) = z01), поэтому из (30) при i = 0 следует, что т = 0. Но при т = 0 (30) противоречит тому, что W(1) = W(2). Это противоречие доказывает теорему.
Благодарности. Автор благодарен В. А. Юрко и М. Ю. Игнатьеву за полезные обсуждения полученных результатов.
Библиографический список
1. Марченко В. А. Операторы Штурма - Лиувилля и их приложения. Киев : Наук. думка, 1977. 331 с.
2. Левитан Б. М. Обратные задачи Штурма - Лиувилля. М. : Наука, 1984. 240 с.
3. Yurko V. A. Method of Spectral Mappings in the Inverse Problem Theory. Inverse and Ill-posed Problems Series. Utrecht : VSP, 2002. 316 p.
4. Голубков А. А., Макаров В. А. Обратная спектральная задача для обобщенного уравнения Штурма - Лиувилля с комплекснозначными коэффициентами // Дифференц. уравнения. 2011. Т. 47, № 10. С. 1498-1502.
5. Голубков А. А., Макаров В. А. Восстановление координатной зависимости тензора диэлектрической проницаемости диагонального вида одномерно-неоднородной среды // Вестн. Моск. ун-та. Сер. 3. Физика. Астрономия. 2010. № 3. С. 32-36.
6. Ангелуц А. А., Голубков А. А., Макаров В. А., Шкуринов А. П. Восстановление спектра диэлектрической проницаемости однородной среды в терагерцовом диапазоне частот по измерениям угловых зависимостей коэффициентов пропускания плоскопараллельной пластины // Письма в ЖЭТФ. 2011. Т. 93, № 4. С. 209-212.
7. Левитан Б. М., Саргсян И. С. Операторы Штурма - Лиувилля и Дирака. М. : Наука, 1988. 432 с.
8. Ишкин Х. К. О необходимых условиях локализации спектра задачи Штурма—Лиувилля на кривой // Матем. заметки. 2005. Т. 78, № 1. С. 72-84. DOI : 10.4213/mzm2563
9. Федорюк М. В. Асимптотические методы для линейных обыкновенных дифференциальных уравнений. М. : Наука, 1983. 352 с.
10. Коддингтон Э. А., Левинсон Н. Теория обыкновенных дифференциальных уравнений. М. : Изд-во иностр. лит., 1958. 475 с.
11. Вазов В. Асимптотические разложения решений обыкновенных дифференциальных уравнений. М. : Мир, 1968. 465 с.
12. Ишкин Х. К. Критерий локализации спектра оператора Штурма — Лиувилля на кривой // Алгебра и анализ. 2016. Т. 28, № 1. С. 52-88.
13.
14.
Ишкин Х. К. О критерии однозначности решений уравнения Штурма - Лиувилля // Матем. заметки. 2008. Т. 84, № 4. С. 552-566. 001 : 10.4213/шгш4173 Ишкин Х. К. О критерии безмонодромности уравнения Штурма — Лиувилля // Матем. заметки. 2013. Т. 94, № 4. С. 552-568. 001 : 10.4213/тгт10319
Образец для цитирования:
Голубков А. А. Обратная задача для операторов Штурма - Лиувилля в комплексной плоскости // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2018. Т. 18, вып. 2. С. 144-156. 001: 10.18500/1816-9791-2018-18-2-144-156
Inverse Problem for Sturm - Liouville Operators in the Complex Plane
A. A. Golubkov
Andrey A. Golubkov, https://orcid.org/0000-0002-5265-1310, The Advanced Educational Scientific Center (faculty) — Kolmogorov's boarding school of Moscow State University, 11, Kremenchugskaya Str., Moscow, 121357, Russia, andrej2501 @yandex.ru
The inverse problem for the standard Sturm-Liouville equation with a spectral parameter p and a potential function, piecewise-entire on a rectifiable curve 7 c C, on which only the starting point is given, is studied for the first time. A function Q that is bounded on a curve 7 is piecewise-entire on it if 7 can be splitted by a finite number of points into parts on which Q coincides with entire functions, different in neighboring parts. The split points, the initial and final points of the curve are called critical points. The problem is to find all the critical points of the curve 7 and the potential on it by the column or row of the transfer matrix P along 7. On the basis of the obtained asymptotics of matrix P for |p| ^ ro, it is proved that if at least one of its elements is bounded for Vp e C, then the curve 7 degenerates to a point after removing all „invisible loops". An „invisible loop" is a loop of the curve 7 (with a given piecewise-entire function) whose knot coincides with two successive critical points. The uniqueness of the solution of the inverse problem for curves without „invisible loops" is proved. On the example of the inverse problem for the equation dx (rtxj dl) + W^) - r(x) A2) y(x) = 0 with a piecewise-entire function q(x) and a piecewise constant function r (x) = 0 on the segment of the real axis, the usefulness of the results obtained in the article is shown for the study of inverse problems for generalized Sturm - Liouville equations, which can be reduced to the type studied in the article.
Key words: Sturm - Liouville equation on a curve, piecewise-entire potential function, transfer matrix, asymptotics, inverse spectral problem.
Acknowledgements: The author is grateful to Vyacheslav A. Yurko and Mikhail Yu. Ignatiev for useful discussions of the results obtained.
References
1. Marchenko V. A. Sturm - Liouville operators and their applications. Birkhauser, 1986. 393 p. (Russ.ed. : Kiev, Naukova Dumka, 1977. 331 p.)
2. Levitan B. M. Inverse Sturm - Liouville problems. Utrecht, VNU Sci. Press, 1987. 246 p. (Russ.ed. : Moscow, Nauka, 1984. 240 p.)
3. Yurko V. A. Method of Spectral Mappings in the Inverse Problem Theory. Inverse and Ill-posed Problems Series. Utrecht, VSP, 2002. 316 p.
4. Golubkov A. A., Makarov V. A. Inverse spectral problem for a generalized Sturm--Liouville equation with complex-valued coefficients. Differential Equations, 2011, vol. 47, no. 10, pp. 1514-1519. DOI : 10.1134/S0012266111100156
5. Golubkov A. A., Makarov V. A. Reconstruction of the coordinate dependence of the diagonal form of the dielectric permittivity tensor of a one-dimensionally inhomogeneous medium. Moscow University Physics Bulletin, 2010, vol. 65, no. 3, pp. 189-194. DOI : 10.3103/S0027134910030070
6. Angeluts A. A., Golubkov A. A., Makarov V. A., Shkurinov A. P. Reconstruction of the spectrum of the relative permittivity of the plane-parallel plate from the angular dependences of its transmission coefficients. JETP Letters, 2011, vol. 93, no. 4, pp. 191— 194. DOI : 10.1134/S0021364011040047
7. Levitan B. M., Sargsyan I. S. Sturm - Liouville and Dirac Operators. Mathematics and Its Applications (Soviet Series), vol. 59. Dordrecht, Springer, 1990. 350 p. (Russ. ed.: Moscow, Nauka, 1988. 432 p.). DOI 10.1007/978-94-011-3748-5
8. Ishkin Kh. K. Necessary Conditions for the Localization of the Spectrum of the Sturm -Liouville Problem on a Curve. Math. Notes, 2005, vol. 78, no. 1, pp. 64-75. DOI: 10.1007/s11006-005-0100-5.
9. Fedoryuk M. V. Asymptotic analysis: linear ordinary differential equations. Berlin, Springer-Verlag, 1993. 363 p. (Russ. ed.: Moscow, Nauka, 1983. 352 p.)
10. Coddington E. A., Levinson N. Theory of ordinary differential equations. New York, McGraw-Hill, 1955. 429 p. (Russ. ed.: Moscow, Izd-vo inostr. lit., 1958. 475 p.)
11. Wasow W. Asymptotic expansions for ordinary differential equations. New York, Dover Publications, 1988. 384 p. (Russ. ed.: Moscow, Mir, 1968. 465 p.)
12. Ishkin Kh. K. Localization criterion for the spectrum of the Sturm - Liouville operator on a curve. St. Petersburg Mathematical Journal, 2017, vol. 28, no. 1, pp. 37-63. DOI: 10.1090/spmj/1438
13. Ishkin Kh. K. On the uniqueness criterion for solutions of the Sturm - Liouville equation. Math. Notes, 2008, vol. 84, no. 4, pp. 515—528. DOI : 10.1134/S000143460809023X
14. Ishkin Kh. K. On a trivial monodromy criterion for the Sturm - Liouville equation. Math. Notes, 2013, vol. 94, no. 4, pp. 508-523. DOI : 10.1134/S0001434613090216
Cite this article as:
Golubkov A. A. Inverse Problem for Sturm - Liouville Operators in the Complex Plane. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2018, vol. 18, iss. 2, pp. 144-156 (in Russian). DOI: 10.18500/1816-9791-2018-18-2-144-156