УДК 517.927.2
DOI: 10.17277/vestnik.2018.04.pp.684-697
ИДЕНТИФИКАЦИЯ РАСПАДАЮЩИХСЯ КРАЕВЫХ УСЛОВИЙ ЗАДАЧИ ШТУРМА-ЛИУВИЛЛЯ
А. М. Ахтямов, Р. Ю. Галимов
Кафедра «Математическое моделирование», ФГБОУ ВО «Башкирский государственный университет», г. Уфа, Республика Башкортостан, Россия; [email protected]
Ключевые слова: задача Штурма-Лиувилля; идентификация краевых условий; обратная спектральная задача; собственные значения.
Аннотация: Рассмотрен вопрос идентификации распадающихся краевых условий задачи Штурма-Лиувилля по ее собственным значениям. Применены два метода идентификации краевых условий - неопределенной системы (метод миноров) и сравнения целых функций. В методе миноров используются только либо два, либо три собственных значения. Однако метод не позволяет ответить на фундаментальный вопрос: для каких классов обыкновенных дифференциальных уравнений и спектральных задач однозначно восстанавливаются виды и параметры краевых условий. Для ответа на данный вопрос использован метод сравнения целых функций. Показано, что в случае несимметрического потенциала задача идентификации краевых условий по всем собственным значениям имеет единственное решение, а в случае симметрического - два решения; при определенных условиях единственное решение может быть получено с использованием только трех собственных значений. Два решения могут быть получены с использованием только двух собственных значений.
Введение
На практике часто возникает задача диагностирования видов и параметров закрепления стержней и струн для того, чтобы определить, не нарушились ли заложенные в проект граничные условия. Кроме того, задача идентификации краевых условий важна при создании безопасных для здоровья человека технических систем. Это связано с тем, что некоторые частоты (инфразвуковые) находятся в опасном для здоровья человека диапазоне частот. Они совпадают с частотами важных органов человека - сердца, почек и т.п. Известно, что воздействие инфра-звуковых колебаний на определенной частоте может вызвать остановку сердца человека. Инфразвуковые частоты ухудшают самочувствие человека, могут вызывать недомогание и даже панические настроения. Поэтому при создании соответствующих технических систем со струнами и стержнями важно уходить от инфразвуковых частот, которые попадают в резонанс с частотами важных органов человека. Изложенные факты требуют создания таких закреплений элементов технических систем, которые давали бы нужный безопасный диапазон частот колебаний основных деталей. В математической постановке обе задачи (диагностирования и ухода от опасных частот) сводятся к задаче идентификации краевых условий по заданным собственным значениям.
Обозначим через Ь следующую задачу Штурма-Лиувилля:
1у = -у" + я(х)у = ^у =5 2 у; (1)
их(у ) = «пу(о)+«12 у'(0) = 0; (2)
и 2 (у) = «23 у(п) = «24 У (п) = 0 (3)
где д(х)е Ь1 (0, п) - вещественная функция; «у, / = 1,2, } = 1...4 - комплексные постоянные.
Близкие обратные задачи рассматривались в теории обратной задачи Штурма-Лиувилля, где по двум (нескольким спектрам) требуется восстановить потенциал ц(х) и краевые условия (2), (3). Обратная задача Штурма-Лиувилля для Ь с распадающимися краевыми условиями достаточно хорошо изучена в работах [1 - 21]. В данных работах восстанавливался потенциал ^(х) в уравнении (1) и краевые условия (2). При этом для восстановления использовались несколько спектров, спектральные данные (спектр и нормировочные числа), функция Вейля и т.п. То есть данные восстановления содержали в себе большую информацию, чем только сам спектр задачи Ь.
Восстановление только краевых условий (с известным дифференциальным уравнением) по собственным значениям началось, по-видимому, в 1990-х годах 20 века [22 - 24]. Оганесяном З. Б. исследовались несколько задач идентификации условий закрепления распределенных механических систем: на обоих концах стержня [23]; круговой [22] и прямоугольной [24] пластин. Однако им восстанавливались лишь коэффициенты канонических условий закрепления. Случай, когда не известен вид канонических условий (то есть, когда не известны все коэффициенты краевых условий) им рассмотрен не был.
В работе [25] (см. также библиографию к данной работе) изучалась идентификация краевых условий, в которых неизвестны все их коэффиценты (известен только ранг матрицы, составленной из коэффициентов краевых условий). Такая задача сводится к идентификации (с точностью до линейных преобразованиий строк) матрицы из коэффициентов краевых условий по ее минорам. В частности, в [26], [27] рассматривались частные случаи идентификации краевых условий задачи Штурма-Лиувилля. В работе [26] показано, что краевые условия Штурма (у'(0)-ку(00) = 0, у'(п) + Ну(п) = 0)задачи Штурма-Лиувилля (1) - (2) по собственным значениям восстанавливаются однозначно с точностью до перестановки к и Н местами. В [30], в случае ^(х) = 0 показано, что задача восстановления вида и параметров распадающихся краевых условий по двум собственным значениям имеет два решения.
Материалы и методы
Методы идентификации как вида, так и параметров краевых условий или условий сопряжения (дифференциальное уравнение считается известным) впервые разработаны А. М. Ахтямовым [25] и до настоящего времени являются оригинальными. В технической диагностике другими исследователями ранее по конечному числу собственных частот восстанавливалась лишь часть коэффициентов краевых условий определенного вида. А в теории обратных спектральных задач краевые условия определенного вида идентифицировались лишь попутно с коэффициентами самих дифференциальных уравнений. При этом для идентификации использовался не один, а два или несколько спектров или же спектр и дополнительно другие спектральные данные (весовые числа, функция Вейля, спектраль-
ная функция и т.п.). В отличие от этих направлений А. М. Ахтямовым предложено идентифицировать всю совокупность краевых условий спектральной задачи как линейную оболочку векторов, компонентами каждого из которых являются коэффициенты соответствующего краевого условия. Такой подход позволяет идентифицировать не только коэффициенты краевого условия определенного (канонического) вида, но и сам вид краевого условия (упругое закрепление, свободный конец, заделка и т.п.). Этот подход в идентификации краевых условий и условий сопряжения достаточно подробно изложен А. М. Ахтямовым и его учениками в многочисленных статьях и обобщающих их монографиях. В данных работах применяются авторские методы идентификации краевых условий: неопределенной системы (метод миноров); сравнения целых функций; введения дополнительных неизвестных; последовательного решения прямых задач и др.
Основными методами, применяемыми в статье, являются методы неопределенной системы (метод миноров) и сравнения целых функций.
Метод неопределенной системы (метод миноров) применяется для идентификации краевых условий по конечному числу собственных значений. На первом этапе данного метода строится характеристический определитель задачи, который представляет собой линейную однородную функцию от миноров Mkm матрицы A, составленной из коэффициентов краевых условий. Множителями при этих неизвестных минорах являются миноры Fkm некоторой известной матрицы F. Элементами данной известной матрицы F служат значения линейно независимых решений обыкновенного дифференциального уравнения. По пространственным переменным эти значения берутся в граничных точках, а спектральный параметр принимает значения, равные собственным. Подставив в разложение характеристического определителя собственные значения в количестве на единицу меньшей количества миноров (максимального порядка) матрицы A краевых условий, получим систему N - 1 (теперь уже линейных) уравнений от N неизвестных миноров Mkm. Пусть ранг данной системы равен N - 1. Тогда вектор из N неизвестных миноров определяется однозначно с точностью до ненулевого множителя. По данному вектору методами алгебраической геометрии и находится матрица краевых условий с точностью до линейных преобразований ее строк (то есть восстанавливаются краевые условия). В работе [28] данный метод уже был использован. Метод основан на восстановлении матрицы с точностью до линейного преобразования строк по ее минорам [28, 29].
Для метода миноров требуется проверка, что ранг некоторой матрицы F равен определенному числу. Проверка этого условия для конкретной краевой задачи не вызывает затруднений. Однако данный метод не позволяет ответить на фундаментальные вопросы: для каких классов обыкновенных дифференциальных уравнений и спектральных задач возможно однозначное восстановление видов и параметров краевых условий и решение задачи идентификации краевых условий двойственно и т.п. Для ответа на эти фундаментальные вопросы в статье используется метод сравнения целых функций, впервые предложенный для решения задач идентификации краевых условий [25], заключающийся в том, что, как правило, характеристический определитель рассматриваемых спектральных задач представляет собой целую функцию порядка одна вторая или четную целую функцию первого порядка. А такие функции восстанавливаются с точностью до ненулевого постоянного множителя по своим нулям (собственным значениям краевой задачи). Поскольку характеристический определитель представляет собой конечную сумму, в которой слагаемыми являются линейно независимые функции от спектрального параметра, то это позволяет доказать, что вектор, составленный из коэффициентов, которые стоят при этих линейно независимых функциях, восстанавливается по всему спектру однозначно с точностью до ненулевого множителя. По данному вектору можно однозначно восстановить краевые условия. Бесконеч-
ный набор собственных частот получить с помощью частотомеров невозможно. Однако данный метод применим к реальным динамическим системам. Он позволяет доказать, что необходимый для метода миноров конечный набор собственных значений найдется. В отличие от метода неопределенной системы, с помощью применения асимптотических формул для линейно независимых решений он позволяет описать классы спектральных задач, краевые условия которых можно однозначно идентифицировать. Полезен он и для предварительного анализа и выяснения того, какое минимальное число собственных значений необходимо для идентификации (это минимальное число частот связано с числом линейно независимых функций N в характеристическом определителе и числом соотношений Плюккера). Метод сравнения целых функций будет развит здесь и для случаев, когда нет однозначности восстановления краевых условий.
Результаты исследования
Обозначим матрицу, составленную из коэффициентов краевых условий (2) - (3), через А, а ее миноры, составленные из ¿-го и у'-го столбцов, через Мц
А =
а Ii 0
М23 =
а 12 0 а
а12
0 0
23 24 0
, м13 =
0
, м24 =
0
а11 0
0
а24
0
a23
, М14 =
аи 0
0
а24
(4)
, М12 = М34 = 0.
Будем считать, что ранг матрицы равен двум: гапкД = 2.
Условимся в дальнейшем задачу типа ¿, но с другими коэффицинтами в уравнении и параметрами в граничных формах, обозначать I. Всюду будем считать, что если некоторый символ обозначает объект из задачи ¿, то символ с волной ~ наверху обозначает аналогичный объект задачи Ь.
Определение 1. Краевые условия задач Ь и Ь назовем (к1 и шп)-смежными, если для миноров Мк1 и Мтп выполняются равенства Мк1 = СМтп, Мтп = СМк1, а для всех остальных миноров выполнены равенства Му = СМц.
Характеристический определитель задачи (1) - (3) (задачи Ь) имеет следующий вид [1]:
Д(Х) = М32 л(п, Х) + М42 у1(п, Х) + ЫХ3 у2 (п, Х) + Мму2 (п, X), (5)
где ух(х, Я) и у2(х,Я) —линейно независимые решения уравнения (1), удовлетворяющие условиям
У1(0Д) = 1, у£(0,Я) = 1, у2(0,Я) = 1, у^(0,Я) = 1 (6)
Справедливы следующие асимптотические формулы:
ух(х,Я) = cossx + О у2(хД) = jsinsx + О у((х,Я) = —s sinsx + 0(—),у2(х,Я) = cossx + О
(7)
для достаточно большого Я = б2 Е И [3, с. 62 - 65].
Обозначим задачу (1) - (3), но с другими коэффициентами краевых условий через I.
Д(Я) =М32У1 (п, Я) + М42 у^ (п, Я) + М13 у2 (п, Я) + М14 у'2 (п, Я). (8)
Из асимптотических представлений (7) для линейно независимых решений для функций ух(п, Я), у2(п,Я) получаем, что функции Д(Я) и Д(Я) являются целыми порядка 1/2. Поэтому из теоремы Адамара следует, что они связаны между собой тождеством
д(х) = сд(1).
(9)
Пусть q(x)Фq{п-x), ц^Ф <|(п - х). В этом случае у2(п,Я) ^ у^пД) [30, см. Лемма 4]. Тогда из (5) и линейной независимости ух(п, Я), у{(п, Я), у2 (п, Я), у2 (п, Я) получаем равенства:
М32 = СЛ?32, М42 = СМ42, М13 = СМ13, М14 = СЛ?14. (10)
Отсюда следует, что краевые условия задач Ъ и I совпадают (с точностью до линейных преобразований строк).
Пусть теперь q(x)=q(п — x), q(x)= <|(п — х). В этом случае у2(п,Я) = _1п,Я [30, см. Лемма 4]. Тогда из (5) и линейной независимости функций у1 (п, Я) = у2 (п, Я) , у{ (п, Я), у2 (п, Я) получаем равенства:
М32 + М14 = С(М32 + М14) , М42 = СМ42, М13 = СЛ?13. (11)
Для нахождения миноров воспользуемся тем, что произвольные числа не могут быть минорами матрицы. Для того, чтобы числа М13 ,М14 ,М23 ,М24 ,М34 были минорами матрицы необходимо и достаточно, чтобы выполнялись так называемые соотношения Плюккера [28]:
М12М34 - М13М24 + М14М23 = 0; (12)
М12М34 - М13М24 + М14М23 = 0. (13)
Тогда отсюда и из (11) получаем два набора равенств:
м13 = см13, м14 = см14, м32 = см32 , М42 = СМ42 (14)
и
м13 = см13, м14 = см32 , м32 = см14, м42 = см42. (15)
Отсюда следует, что краевые условия задач Ь и I либо совпадают (с точностью до линейных преобразований строк), либо являются (14 и 32) смежными.
Таким образом, верна следующая теорема.
Теорема 1. Пусть спектры задач Ь и Ь с распадающимися краевыми условиями совпадают с учетом их алгебраических кратностей, гапкД=2:
Тогда:
1. Если q(x)Ф ^(п — х), (¡^Ф с[(п — х), то матрицы коэффициентов краевых условий А = (сг^)2х4 и А = (сг^)2х4 совпадают с точностью до линейных преобразований строк. То есть, в случае несимметрического потенциала, задача идентификации распадающихся краевых условий по всем собственным значениям имеет единственное решение.
2. Если q(x)= ^(п — х), q(x)= <?(п — х), то либо матрицы коэффициентов краевых условий А = (а^)2х4 и А = (а^)2х4 совпадают с точностью до линейных преобразований строк, либо являются (14 и 32) смежными. То есть в случае симметрического потенциала задача идентификации распадающихся краевых условий по всем собственным значениям имеет два решения.
Рассмотрим теперь задачу идентификации краевых условий по конечному набору собственных значений.
Обозначим через Q матрицу следующего вида:
(16)
у1(п,Я1) у((п,Я1) у2(п,Я1) у2(п,Ях) Q = У1(п,Я2) у((п,Я2) у2(п,Я2) у2(п,Я2) У1(п,Я3) у((п,Я3) у2(п,Я3) у2(п,Я3) через Qj - минор матрицы Q, полученный вычеркиванием у-го столбца матрицы Q.
Теорема 2. Если д(х)Ф д(п — х), гапкД=2, три собственных значения
Я,-, (/ = 1,2,3) задачи Ь удовлетворяют условию:
mпkQ=3,
(17)
то задача идентификации краевых условий по этим трем собственным значениям имеет единственное решение, которое представляется формулами:
а) если M13 Ф 0, то A =
б) если M14 Ф 0, то A =
1 M 23 0
M13
0 M12 M13
1 M 24 M 34
M14 M14
0 M12 M13
M-
34
Mi
13
M,
14
то A = M13 1 0
M 23
- M12 0 M23
M14 1 M34
то A = M 24 M 24
- M12 0 M23
M
34
M 23
M 24
0
M24
Причем миноры в представлениях а) - г) для матрицы А даются следующими равенствами:
М12 = 0, М34 = 0, М32 = К}^ М42 = -М}2, М13 = £<2з, М14 = -1(}4, (18)
Доказательство. Так как у2(тс,Я) Ф ух(п,Я) [30, см. Лемма 4], и три собственных значения Я, (=1, 2, 3) задачи Ь являются корнями функции (5), то они удовлетворяют следующей системе уравнений относительно неизвестных
М32 , М42, М13, М14 :
Д(Х ] ) = М 32 уДп, X ]) + М 42 у^п, X ] ) + Мп 72 (п, X ] ) + М14У2 (п, X ] ) = 0. (19)
Матрица данной системы совпадает с матрицей Согласно условию теоремы, ранг матрицы 0 равен трем. Поэтому система уравнений (19) имеет единственное с точностью до ненулевого множителя t решение (18). По минорам (18) с помощью методов работ [25], [28] находится матрица А. В зависимости от того, какой из определителей отличен от нуля, она дается формулами а) - г). Теорема доказана.
Пример 1. Пусть собственные значения задачи (1), (2) с д(х)=2х+3, Я1 = 5,132512, Я2 = 8,59040, Я3 = 12,71359. Разложив линейно независимые решения у1 (п, Я) и у2 (п, Я) в ряд Тейлора по х и X и подставив частичную сумму ряда из первых 120 членов ряда в (16) и вычислив соответствующие миноры с точностью до семи значащих цифр, получим: Q1 = —18,28142, Q2 = 0, Q3 = -9.14071, = 0.
м12 = 0,М34 = 0, М23 = -М32 = -tQ1 = 18,28142 t, М24 = -М42 = -tQ2 = 0, М13 = tQ3 = -9,14071t, М14 = -tQ4 = 0.
(20)
Положив t = 18,28142 1, получим более простые представления для миноров:
м12 = М14 = М24 = М34 = 0, М23 = 1, м13 = --.
0
Так как М23 = 1^0 (случай в) из теоремы 2), то матрица А с точностью до линейных преобразований строк совпадает с матрицей
A =
М-
13
1
0
М-
34
М 23
- М12 0 М23
M23 М 24
1
0
0 0 10
R =
Таким образом, задача идентификации краевых условий в данном случае имеет единственное решение. Это решение представляет собой следующие краевые условия
у(0)- 2/(0) = 0, у(п) = 0.
Обозначим через Я матрицу следующего вида
у1(п,Я1) у2(л,Х1) (21)
у1(п,Я2) у[(п,Х2) у2(тс,Х2)
через Я,- обозначим минор матрицы Я, полученный вычеркиванием ¡-го столбца матрицы Я.
Теорема 3. Если д(х)=д(п — х), гапкА=2, два собственных значения Х,-(¡=1,2) задачи Ь удовлетворяют условию:
гапкй = 2, (22)
то задача идентификации краевых условий по этим двум собственным значениям имеет два решения, которые представляется формулами а) - г) из теоремы 2, в зависимости от того, какой из миноров отличен от нуля. Причем миноры в представлениях а) - г) для матрицы А даются следующими равенствами:
М12 = 0,М34 =Q,M32=2(fii±
М42 = tR2,M13 = tR3,M14 =2\Ri +
Äi + 4 R2R-
+ 4 R2R,
(23)
Доказательство. Два собственных значения Х^ , ¡=1, 2, задачи Ь являются корнями функции (5). Следовательно, они удовлетворяют следующей системе уравнений относительно неизвестных Х1 = М12 + М34, Х2 = М32, Х3 = М42, х4 = М13 х5 = М14 :
(М32 + М14)у1(п,А/) + М42 yi(n,A,-) + М13 у2(п,Я,) = 0.
(24)
Матрица этой системы совпадает с матрицей Я (21). Согласно условию теоремы, ранг матрицы Я равен двум. Поэтому система уравнений (24) имеет единственное с точностью до ненулевого множителя решение:
М32 + М14 = tR1, М42 = М13 = tR3. (25)
Из соотношений Плюккера (9), (4) и (25) получаем
М32М14 = -£2Д2Д3. (26)
Из равенств (25), (26) и обратной теоремы Виета следует, что М32 и М14 являются корнями квадратного уравнения
z2 - tR,z-t2RyR3 = 0.
Откуда
М32 =2\Ri±
Д? + 4 RyR,
,М14 =2lfii +
Д? + 4 RyR-
(27)
Из (27), (25) и вытекает утверждение теоремы.
2
Пример 2. Пусть собственные значения задачи (1), (2) с д(х)= 0, Я1 = 2,549144, Я2 = 6,560591. Линейно независимыми решениями задачи (1), (2)
1 (xVX)
VI
. Подставив
с q(x)= 0 являются функции ух(п,Я) = cos (xVX) и у2(пД) = :
их в (25) и вычислив соответствующие миноры Rj с точностью до семи значащих цифр, получим: R1 = -0,9187470, R2 = 0,R3 = -0,4593735. Отсюда, а также из (25) и (27) получаем
М12 = 0, М34 = 0;
М,
= 2(fii ± + 4R2R3 ) = 2 (—0,9187470 ± V(-0,9187470)2) =
(28)
= {-0,918747t или 0};
Мл
= 2 + = 2 (-0,9187470 + 7(-0,9187470)2) =
= {0 или - 0,918747t}; М24 = tR2 = 0, М13 = tR3 = -0,4593735t.
Положив t = —0,4593735_1, из (28) получим более простые представления для миноров:
М12 = 0, М34 =0, М32 = {2 или 0}, М14 = {0 или 2 },
М24 = 0, М13 = 1. (29)
Так как в первом случае М23 = —2^0 (случай в) из теоремы 2), а во втором М14 = 2^0 (случай б) из теоремы 2), то матрица А с точностью до линейных преобразований строк будет иметь два представления
M 13 1 0 M 34
= M 23 M 23
- M 12 0 M 23 M 24
1 M 24 M 34 0
A2 = M14 M14
0 M12 M13 M14
-1 1 2
0
0
0 - 2
10 0 0 0 0 12
Таким образом, задача идентификации краевых условий в данном случае имеет два решения. Эти решения представляет собой следующие краевые условия:
или
y(0)- 2 /(0) = 0, y(n) = 0 y(0) = 0, y(n) + 2y'(n) = 0.
Физически это означает, что закрепления на концах однородной струны находятся с точностью до перестановок их местами.
Заключение
Рассмотрен вопрос идентификации распадающихся краевых условий задачи Штурма-Лиувилля по ее собственным значениям. С помощью метода сравнения целых функций показано, что в случае несимметрического потенциала задача их идентификации распадающихся краевых условий по всем собственным значениям в случае несимметрического потенциала имеет единственное решение, а в случае симметрического потенциала - два решения.
Методом неопределенной системы доказано, что для задачи идентификации распадающихся краевых условий показано, что единственное решение может быть получено по трем собственным значениям, а два решения - по двум собственным значениям. Приведены соответствующие примеры восстановления краевых условий. Исследования могут быть продолжены и для восстановления краевых задач с более общими краевыми условиями
Работа выполнена при финансовой поддержке РФФИ и Правительства Республики Башкортостан (проекты 18-51-06002-Аз-а, 18-01-00250-а, 17-41-020230-р_а и Госзадание № АААА-А17-117040510250-6).
Список литературы
1. Марченко, В. А. Операторы Штурма-Лиувилля и их приложения / В. А. Марченко. - Киев : Наукова думка, 1977. - 331 с.
2. Левитан, Б. М. Обратные задачи Штурма-Лиувилля / Б. М. Левитан. -М. : Наука, 1984. - 240 с.
3. Наймарк, М. А. Линейные дифференциальные операторы / М. А. Наймарк. -2-е изд., перераб. и доп. - М. : Наука, 1969. - 528 с.
4. Коротяев, Е. Л. Обратная задача Штурма-Лиувилля со смешанными краевыми условиями / Е. Л. Коротяев, Д. С. Челкак // Алгебра и анализ. - 2009. - Т. 21, № 5. - С. 114 - 137.
5. Mamedov, Kh. R. A Uniqueness Theorem for a Sturm-Liouville Equation with Spectral Parameter in Boundary Conditions / Kh. R. Mamedov, F. A. Cetinkaya // Appl. Math. Inf. Sci. - 2015. - Vol. 9, No. 2. - P. 981 - 988.
6. Panakhov, E. S. Reconstruction Formula for the Potential Function of Sturm-Liouville Problem with Eigenparameter Boundary Condition / E. S. Panakhov, H. Koyunbakan, Ic. Unal // Inverse Problems in Science and Engineering. - 2010. -Vol. 18, No. 1. - P. 173 - 180.
7. Савчук, А. М. Обратные задачи для оператора Штурма-Лиувилля с потенциалами из пространств Соболева. Равномерная устойчивость / А. М. Савчук, А. А. Шкаликов // Функциональный анализ и его приложения. - 2010. - T. 44, № 4. - С. 34 - 53.
8. Mizrak, O. Characteristic Properties of Scattering Data of a Boundary Value Problem / O. Mizrak, Kh. R. Mamedov, A. M. Akhtyamov // Filomat. - 2017. - Vol. 31, No. 12. - P. 3945 - 3951.
9. Akhtyamov, A. M. Inverse Problem for the Diffusion Operator with Symmetric Functions and General Boundary Conditions / A. M. Akhtyamov, V. A. Sadovnichy, Ya. T. Sultanaev // Eurasian Mathematical Journal. - 2017. - Vol. 8, No. 1. - P. 10 - 22.
10. Del Rio, R. Inverse Problems for Jacobi Operators IV: Interior Mass-Spring Perturbations of Semi-Infinite Systems / R. Del Rio, L. O. Silva, M. Kudryavtsev // Inverse Problems. - 2017. - Vol. 33, No. 5. - P. 5 - 10. doi: 10.1088/1361-6420/aa6808
11. Monk, P. An Inverse Acoustic Waveguide Problem in the Time Domain / P. Monk, V. Selgas // Inverse Problems. - 2016. - Vol. 32, No. 5. - P. 5 - 10. doi: 10.1088/0266-5611/32/5/055001
12. Its, A. А Riemann-Hilbert Approach to the Inverse Problem for the Stark Operator on the Line / A. Its, V. A. Sukhanov // Inverse Problems. - 2016. - Vol. 32, No. 5. - P. 16 - 22. doi: 10.1088/0266-5611/32/5/055003
13. Kutsenko, A. A. Recovery of Defects from the Information at Detectors / A. A. Kutsenko // Inverse Problems. - 2016. - Vol. 32, No. 5. - P. 28 - 36. doi: 10.1088/0266-5611/32/5/055005
14. Horvath, M. Discrete Inverse Problems for the Schrödinger Operator on the Multi-Dimensional Square Lattice with Partial Cauchy Data / M. Horvath, Z. Marko // Inverse Problems. - 2016. - Vol. 32, No. 5. - P. 36 - 42. doi: 10.1088/0266-5611/32/5/ 055006
15. Rostamian, M. A Meshless Method for Solving 1D Time-Dependent Heat Source Problem / M. Rostamian, A. Shahrezaee // Inverse Problems in Science and Engineering. - 2018. - Vol. 26, No. 1. - P. 51 - 82.
16. Wei, Z. The Inverse Discrete Transmission Eigenvalue Problem for Absorbing Media / Z. Wei, G. Wei // Inverse Problems in Science and Engineering. - 2018. -Vol. 26, No. 1. - P. 83 - 99.
17. Mirsalimov, V. M. Inverse Problem of Fracture Mechanics for a Circular Disc under Mixed Boundary Conditions / V. M. Mirsalimov, S. H. Hasanov // Inverse Problems in Science and Engineering. - 2017. - Vol. 25, No. 11. - P. 1547 - 1559.
18. Fotopoulos, G. Inverse Scattering with Fixed Observation Angle Data in 2D / G. Fotopoulos, M. Harju // Inverse Problems in Science and Engineering. - 2017. -Vol. 25, No. 10. - P. 1492 - 1507.
19. Buterin, S. A. On the Half Inverse Spectral Problem for an Integra-Differential Operator / S. A. Buterin, M. Sat // Inverse Problems in Science and Engineering. -2017. - Vol. 25, No. 10. - P. 1508 - 1518.
20. Neamaty, A. Numerical Solution of Inverse Nodal Problem with an Eigenvalue in the Boundary Condition / A. Neamaty, S. Akbarpoor // Inverse Problems in Science and Engineering. - 2017. - Vol. 25, No. 7. - P. 978 - 994.
21. Qöl, A. Inverse Spectral Problem for Dirac Operator with Discontinuous Coefficient and Polynomials in Boundary Condition / A. Qöl // Inverse Problems in Science and Engineering. - 2016. - Vol. 24, No. 2. - P. 234 - 246.
22. Гнуни, В. Ц. Определение граничных условий круглой кольцевой пластинки по заданным частотам собственных колебаний / В. Ц. Гнуни, З. Б. Оганисян // Известия академии наук Армении. Механика. - 1991. - Т. 44, № 5. - С. 9 - 16.
23. Оганисян, З. Б. Об одной задаче восстановления граничных условий на концах стержня при заданном спектре частот собственных поперечных колебаний / З. Б. Оганисян // Сб. «Вопросы оптимального управления, устойчивости и прочности механических систем» (науч. тр. конф.). - Ереван, 1997. - С. 159 - 162.
24. Оганисян, З. Б. Об одной задаче восстановления граничных условий на краях пластинки при заданном спектре частот собственных поперечных колебаний / З. Б. Оганисян // Ученые записки ЕГУ. - 1991. - № 1. - С. 45 - 50.
25. Ахтямов, А. М. Теория идентификации краевых условий и ее приложения / А. М. Ахтямов. - М. : Физматлит, 2009. - 272 с.
26. Ахтямов, А. М. К единственности решения одной обратной спектральной задачи / А. М. Ахтямов // Дифференц. уравнения. - 2003. - Т. 39, № 8. - С. 1011 - 1015.
27. Ахтямов, А. М. Идентификация краевых условий на обоих концах струны по собственным частотам колебаний / А. М. Ахтямов, И. М. Утяшев // Акустический журнал. - 2015. - Т. 61, № 6. - С. 647 - 655. doi: 10.7868/S0320791915050019
28. Akhtyamov, A. M. On Reconstruction of a Matrix by its Minors / A. M. Akh-tyamov, M. Amram, A. Mouftakhov // International Journal of Mathematical Education in Science and Technology. - 2018. - Vol. 49, No. 2. - P. 268 - 281.
29. Cramer's Rule for Nonsingular m*n Matrices / A. M. Akhtyamov [et al.] // The Teaching of Mathematics. - 2017. - Vol. XX, No. 1. - P. 13 - 19.
30. Юрко, В. А. Обратная задача для дифференциальных операторов второго порядка с регулярными краевыми условиями / В. А. Юрко // Матем. заметки. -1975. - Т. 18, № 4. - С. 569 - 576.
Identification of Disintegrating Boundary Conditions of the Sturm-Liouville Problem
A. M. Akhtyamov, R. Yu. Galimov
Department of Mathematical Modeling, Bashkir State University, Ufa, Republic of Bashkortostan, Russia; [email protected]
Keywords: Sturm-Liouville problem; identification of boundary conditions; inverse spectral problem; own values.
Abstract: The problem of identifying decaying boundary conditions of the Sturm-Liouville problem by its eigenvalues is investigated. In this paper, we use two methods for identifying the boundary conditions - the method of an uncertain system (the minors method), and the method of comparing entire functions. In the minor method, only two or three eigenvalues are used. However, this method does give the answer to the fundamental question: what classes of ordinary differential equations and spectral problems are the uniquely restored types and parameters of boundary conditions. To answer this question, the article uses the method of comparing entire functions. It is shown that in the case of an asymmetric potential, the problem of identifying the boundary conditions for all eigenvalues has a unique solution, and in the case of a symmetric potential - two solutions. It is also shown that under certain conditions a single solution can be obtained using only three eigenvalues. Two solutions can be obtained using only two eigenvalues.
References
1. Marchenko V.A. Operatory Shturma-Liuvillya i ikh prilozheniya [Sturm-Liouville Operators and their Applications], Kiev: Naukova dumka, 1977, 331 p. (In Russ.)
2. Levitan B.M. Obratnyye zadachi Shturma-Liuvillya [Inverse Sturm-Liouville Problems], Moscow: Nauka, 1984, 240 p. (In Russ.)
3. Naymark M.A. Lineynyye differentsial'nyye operatory [Linear Differential Operators], Moscow: Nauka, 1969, 528 p. (In Russ.)
4. Korotyayev Ye.L., Chelkak D.S. [The Sturm-Liouville Inverse Problem with Mixed Boundary Conditions], Algebra i analiz [Algebra and Analysis], 2009, vol. 21, no. 5, pp. 114-137. (In Russ.)
5. Mamedov Kh.R., Cetinkaya F.A. A Uniqueness Theorem for a Sturm-Liouville Equation with Spectral Parameter in Boundary Conditions, Appl. Math. Inf. Sci., 2015, vol. 9, no. 2, pp. 981-988.
6. Panakhov E.S., Koyunbakan H., Unal Ic. Reconstruction Formula for the Potential Function of Sturm-Liouville Problem with Eigenparameter Boundary Condition, Inverse Problems in Science and Engineering, 2010, vol. 18, no. 1, pp. 173-180.
7. Savchuk A.M., Shkalikov A.A. [Inverse problems for the Sturm-Liouville Operator with Potentials from Sobolev Spaces. Uniform Stability], Funktsional'nyy analiz i yego prilozheniya [Functional Analysis and its Applications], 2010, vol. 44, no. 4, pp. 34-53. (In Russ.)
8. Mizrak O., Mamedov Kh.R., Akhtyamov A.M. Characteristic Properties of Scattering Data of a Boundary Value Problem, Filomat, 2017, vol. 31, no. 12, pp. 3945-3951.
9. Akhtyamov A.M., Sadovnichy V.A., Sultanaev Ya.T. Inverse Problem for the Diffusion Operator with Symmetric Functions and General Boundary Conditions, Eurasian Mathematical Journal, 2017, vol. 8, no. 1, pp. 10-22.
10. Del Rio R., Silva L.O., Kudryavtsev M. Inverse Problems for Jacobi Operators IV: Interior Mass-Spring Perturbations of Semi-Infinite Systems, Inverse Problems, 2017, vol. 33, no. 5, pp. 5 - 10, doi: 10.1088/1361-6420/aa6808
11. Monk P., Selgas V. An Inverse Acoustic Waveguide Problem in the Time Domain, Inverse Problems, 2016, vol. 32, no. 5, pp. 5-10, doi: 10.1088/02665611/32/5/055001
12. Its A., Sukhanov V.A. A Riemann-Hilbert Approach to the Inverse Problem for the Stark Operator on the Line, Inverse Problems, 2016, vol. 32, no. 5, pp. 16-22, doi: 10.1088/0266-5611/32/5/055003
13. Kutsenko A.A. Recovery of Defects from the Information at Detectors, Inverse Problems, 2016, vol. 32, no. 5, pp. 28-36, doi: 10.1088/0266-5611/32/5/055005
14. Horvath M., Marko Z. Discrete Inverse Problems for the Schrödinger Operator on the Multi-Dimensional Square Lattice with Partial Cauchy Data, Inverse Problems, 2016, vol. 32, no. 5, pp. 36-42, doi: 10.1088/0266-5611/32/5/055006
15. Rostamian M., Shahrezaee A. A Meshless Method for Solving 1D Time-Dependent Heat Source Problem, Inverse Problems in Science and Engineering, 2018, vol. 26, no. 1, pp. 51-82.
16. Wei Z., Wei G. The Inverse Discrete Transmission Eigenvalue Problem for Absorbing Media, Inverse Problems in Science and Engineering, 2018, vol. 26, no. 1, pp. 83-99.
17. Mirsalimov V.M., Hasanov S.H. Inverse Problem of Fracture Mechanics for a Circular Disc under Mixed Boundary Conditions, Inverse Problems in Science and Engineering, 2017, vol. 25, no. 11, pp. 1547-1559.
18. Fotopoulos G., Harju M. Inverse Scattering with Fixed Observation Angle Data in 2D, Inverse Problems in Science and Engineering, 2017,vol. 25, no. 10, pp. 1492-1507.
19. Buterin S.A., Sat M. On the Half Inverse Spectral Problem for an Integro-Differential Operator, Inverse Problems in Science and Engineering, 2017, vol. 25, no. 10, pp. 1508-1518.
20. Neamaty A., Akbarpoor S. Numerical Solution of Inverse Nodal Problem with an Eigenvalue in the Boundary Condition, Inverse Problems in Science and Engineering, 2017, vol. 25, no. 7, pp. 978-994.
21. Çôl A. Inverse Spectral Problem for Dirac Operator with Discontinuous Coefficient and Polynomials in Boundary Condition, Inverse Problems in Science and Engineering, 2016, vol. 24, no. 2, pp. 234-246.
22. Gnuni V.Ts., Oganisyan Z.B. [Determination of the Boundary Conditions of a Circular Annular Plate for Given Natural Frequencies], Izvestiya akademii nauk Armenii. Mekhanika [Proceedings of National Academy of Sciences of Armenia. Mechanics], 1991, vol. 44, no. 5, pp. 9-16. (In Russ.)
23. Oganisyan Z.B. Sbornik «Voprosy optimal'nogo upravleniya, ustoychivosti i prochnosti mekhanicheskikh sistem» (nauchnyye trudy konferentsii) [Collection "Questions of Optimal Control, Stability and Strength of Mechanical Systems" (Conference Proceedings)], Yerevan, 1997, pp. 159-162. (In Russ.)
24. Oganisyan Z.B. [On a Problem of Restoring Boundary Conditions at the Edges of a Plate for a Given Frequency Spectrum of Natural Transverse Oscillations], Uchenyye zapiski YEGU [Scholarly Notes YSU], 1991, no. 1, pp. 45-50. (In Russ.)
25. Akhtyamov A.M. Teoriya identifikatsii krayevykh usloviy i yeye prilozheniya [The Theory of Identification of Boundary Conditions and its Applications], Moscow: Fizmatlit, 2009, 272 p. (In Russ.)
26. Akhtyamov A.M. [To the Uniqueness of the Solution of One Inverse Spectral Problem], Differents. uravneniya [Differ. the Equations], 2003, vol. 39, no. 8, pp. 1011-1015. (In Russ.)
27. Akhtyamov A.M., Utyashev I.M. [Identification of Boundary Conditions at Both Ends of a String by the Eigenfrequencies of Vibrations], Akusticheskiy zhurnal [Akustichsky Journal], 2015, vol. 61, no. 6, pp. 647-655, doi: 10.7868/ S0320791915050019 (In Russ., abstract in Eng.)
28. Akhtyamov A.M., Amram M., Mouftakhov A. On Reconstruction of a Matrix by its Minors, International Journal of Mathematical Education in Science and Technology, 2018, vol. 49, no. 2, pp. 268-281.
29. Akhtyamov A.M., Amram M., Dagan M., Mouftakhov A. Cramer's Rule for Nonsingular m*n Matrices, The Teaching of Mathematics, 2017, vol. XX, no. 1, pp. 13-19.
30. Yurko V.A., [An Inverse Problem for Second-Order Differential Operators with Regular Boundary Conditions], Matematicheskiye zametki [Mathematical Notes], 1975, vol. 18, no. 4, pp. 569-576. (In Russ.)
Identifizierung der verteilenden Randbedingungen der Sturm-Liouville-Gleichung
Zusammenfassung: Es wird das Problem untersucht, verteilende Randbedingungen der Sturm-Liouville-Gleichung anhand ihrer eigenen Werte zu identifizieren. In der vorliegenden Arbeit verwenden wir zwei Methoden zur Identifizierung der Randbedingungen - die Methode eines unsicheren Systems (die Minor-Methode), die Methode des Vergleichs gesamter Funktionen. In der Minor-Methode werden nur zwei oder drei Eigenwerte verwendet. Dieses Verfahren erlaubt es jedoch nicht, die grundlegende Frage zu beantworten: für welche Klassen gewöhnlicher Differentialgleichungen und Spektralprobleme die Arten und Parameter der Randbedingungen eindeutig wiederhergestellt werden. Um diese Frage zu beantworten, wird in dem Artikel die Methode des Vergleichs ganzer Funktionen verwendet. Es wird gezeigt, dass das Problem der Identifizierung von Randbedingungen für alle Eigenwerte bei einem asymmetrischen Potential eine eindeutige Lösung hat, und bei einem symmetrischen Potential zwei Lösungen. Es wird auch gezeigt, dass unter bestimmten Bedingungen eine einzelne Lösung mit nur drei Eigenwerten erhalten werden kann. Zwei Lösungen können mit nur zwei Eigenwerten erhalten werden.
Identification des conditions de rupture de la théorie de Sturm-Liouville
Résumé: Est étudiée la tâche de l'identification des conditions de limite décroissantes de la théorie de Sturm-Liouvillepar ses propres valeurs. Dans cetarticle, sont utilisées deux méthodes d'identification des conditions de limite - la méthode d'un système non défini (méthode des mineurs) et la méthode de comparaison des fonctions complètes. Dans la méthode des mineurs sont utilisée seulement deux ou troistrois valeurs propres. Cependant, cette méthode ne permet pas de répondre à la question fondamentale: pour quelles classes d'équations différentielles ordinaires et des tâches spectrales les types et les paramètres des conditions de limite sont restaurés d'une
manière ou d'une autre. Pour répondre à cette question on utilise la méthode de comparaison des fonctions complètes. Il est démontré que, dans le cas d'un potentiel asymétrique, l'identification des conditions de limite pour toutes ses propres valeurs est la seule solution et, dans le cas d'un potentiel symétrique, il y a deux solutions. Il est également démontré que, dans certaines conditions, la seule solution peut être obtenue en utilisant seulement trois valeurs propres. Deux solutions peuvent être obtenues en utilisant seulement deux valeurs propres.
Авторы: Ахтямов Азамат Мухтарович - доктор физико-математических наук, профессор кафедры математического моделирования; Галимов Рустам Юмадилович - аспирант кафедры математического моделирования, ФГБОУ ВО «Башкирский государственный университет», г. Уфа, Республика Башкортостан, Россия.
Рецензент: Спивак Семён Израилевич - доктор физико-математических наук, профессор, заведующий кафедрой математического моделирования, ФГБОУ ВО «Башкирский государственный университет», г. Уфа, Республика Башкортостан, Россия.