Это не официальный сайт wikipedia.org 01.01.2023

Задача Штурма — Лиувилля — Википедия

Задача Штурма — Лиувилля

(перенаправлено с «Задача Штурма-Лиувилля»)

Задача Шту́рма — Лиуви́лля, названная в честь Жака Шарля Франсуа Штурма и Жозефа Лиувилля, состоит в отыскании нетривиальных (то есть отличных от тождественного нуля) решений на промежутке ( a , b ) уравнения Штурма — Лиувилля

L [ y ] = λ ρ ( x ) y ( x ) ,

удовлетворяющих однородным краевым (граничным) условиям

α 1 y ( a ) + β 1 y ( a ) = 0 , α 1 2 + β 1 2 0 ; α 2 y ( b ) + β 2 y ( b ) = 0 , α 2 2 + β 2 2 0 ;

и значений параметра λ , при которых такие решения существуют.

Оператор L [ y ] здесь — это действующий на функцию y ( x ) линейный дифференциальный оператор второго порядка вида

L [ y ] d d x [ p ( x ) d y d x ] + q ( x ) y ( x )

(оператор Штурма — Лиувилля или оператор Шрёдингера), x  — вещественный аргумент.

Функции p ( x ) , p ( x ) , q ( x ) , ρ ( x ) предполагаются непрерывными на ( a , b ) , кроме того функции p ( x ) , ρ ( x ) положительны на ( a , b ) .

Искомые нетривиальные решения называются собственными функциями этой задачи, а значения λ , при которых такое решение существует — её собственными значениями (каждому собственному значению соответствует собственная функция).

Постановка задачиПравить

Вид уравненияПравить

Если функции ρ   и p   дважды непрерывно дифференцируемы и положительны на отрезке [ a , b ]   и функция q   непрерывна на [ a , b ]  , то уравнение Штурма — Лиувилля вида

( p ( x ) y ) q ( x ) y + λ ρ ( x ) y = 0  

при помощи преобразования Лиувилля приводится к виду[1][2]

y + q 1 ( x ) y = λ y . ( 1 )  

Поэтому часто рассматривают уравнение Штурма — Лиувилля в виде (1), функцию q 1 ( x )   называют потенциалом[3][4]. Изучаются задачи Штурма — Лиувилля с потенциалами из разных классов функций: непрерывными, L   (суммируемыми), L 2   и других.

Виды краевых условийПравить

  • Условия Дирихле y ( a ) = y ( b ) = 0.  
  • Условия Неймана y ( a ) = y ( b ) = 0  
  • Условия Робена y ( a ) h y ( a ) = 0 , y ( b ) + H y ( b ) = 0.  
  • Смешанные условия: условия разных видов в разных концах отрезка [ a , b ]  .
  • Распадающиеся краевые условия общего вида
α 1 y ( a ) + β 1 y ( a ) = 0 , α 1 2 + β 1 2 0 ; α 2 y ( b ) + β 2 y ( b ) = 0 , α 2 2 + β 2 2 0.  
  • Периодические условия y ( a ) = y ( b ) , y ( a ) = y ( b )  .
  • Антипериодические условия y ( a ) = y ( b ) , y ( a ) = y ( b )  .
  • Общие краевые условия
a i 1 y ( a ) + a i 2 y ( a ) + a i 3 y ( b ) + a i 4 y ( b ) = 0 , i = 1 , 2.  

В последнем случае обычно накладываются дополнительные условия регулярности на коэффициенты a i j  .[3][5]

Для удобства произвольный отрезок [ a , b ]   часто переводят в отрезок [ 0 , l ]   или [ 0 , π ]   с помощью замены переменной.

Оператор Штурма — ЛиувилляПравить

Оператор Штурма — Лиувилля

L y = 1 ρ ( x ) ( d d x [ p ( x ) d d x y ] q ( x ) y )  

представляет собой частный случай линейного дифференциального оператора[6]

p 0 ( x ) y ( n ) + p 1 ( x ) y ( n 1 ) + + p n ( x ) y .  

Область определения оператора L   состоит из дважды непрерывно дифференцируемых на отрезке [ a , b ]   функций y  , удовлетворяющих краевым условиям задачи Штурма — Лиувилля. Таким образом, задачу Штурма — Лиувилля можно рассматривать как задачу на собственные значения и собственные функции оператора L  : L y = λ y  . Если функции p ,   q ,   ρ   и коэффициенты краевых условий вещественные, оператор L   является самосопряжённым в гильбертовом пространстве L 2 ( [ a , b ] , ρ ( x ) d x )  . Следовательно, его собственные значения вещественны и собственные функции ортогональны с весом ρ ( x )  .

Решение задачиПравить

ПримерПравить

Решение задачи Штурма — Лиувилля с нулевым потенциалом:

y = λ y , ( 2 )  
y ( 0 ) = y ( l ) = 0  

может быть найдено в явном виде[7]. Пусть λ = ρ 2  . Общее решение уравнения (2) при каждом фиксированном λ   имеет вид

y ( x ) = A sin ρ x ρ + B cos ρ x ( 3 )  

(в частности, при ρ = 0   (3) дает y ( x ) = A x + B  ). Из y ( 0 ) = 0   следует B = 0  . Подставляя (3) в краевое условие y ( l ) = 0  , получаем A sin ρ l ρ = 0  . Так как мы ищем нетривиальные решения, то A 0  , и мы приходим к уравнению на собственные значения

sin ρ l ρ = 0.  

Его корни ρ n = π n l  , следовательно, искомые собственные значения имеют вид

λ n = ( π n l ) 2 , n = 1 , 2 , 3 ,  

а соответствующие им собственные функции суть

y n ( x ) = sin π n l x , n = 1 , 2 , 3 ,  

(с точностью до постоянного множителя).

Общий случайПравить

В общем случае любое решение уравнения Штурма — Лиувилля

y + q ( x ) y = λ y ( 4 )  

представимо в виде линейной комбинации

y ( x ) = A S ( x , λ ) + B C ( x , λ ) ( 5 )  

его решений S ( x , λ )   и C ( x , λ )  , удовлетворяющих начальным условиям

S ( 0 , λ ) = C ( 0 , λ ) = 0 , S ( 0 , λ ) = C ( 0 , λ ) = 1  .

Решения S ( x , λ )   и C ( x , λ )   образуют фундаментальную систему решений уравнения (4) и являются целыми функциями по λ   при каждом фиксированном x  . (При q ( x ) 0   S ( x , λ ) = sin ρ x  , C ( x , λ ) = cos ρ x  , ρ = λ  ). Подставляя (5) в краевые условия y ( 0 ) = y ( π ) = 0  , получаем, что собственные значения совпадают с нулями характеристической функции

Δ ( λ ) = S ( π , λ ) ,  

аналитической во всей λ  -плоскости.[4]

В общем случае собственные значения и собственные функции не могут быть найдены в явном виде, однако для них получены асимптотические формулы:

λ n = n + c n + O ( 1 n 2 ) , c = 1 2 π 0 π q ( τ ) d τ ,  
y n ( x ) = sin n x + O ( 1 n 2 ) ,  

(в случае непрерывного на [ 0 , π ]   потенциала q ( x )  ).[8] При больших n   собственные значения и собственные функции близки к собственным значениям и собственным функциям задачи из примера с нулевым потенциалом.

Свойства собственных значений и собственных функцийПравить

  • Существует бесконечное счетное множество собственных значений: λ 1 < λ 2 < < λ n <  
  • Каждому собственному значению λ n   соответствует единственная с точностью до постоянного множителя собственная функция y n  .
  • Все собственные значения вещественны.
  • В случае граничных условий y ( a ) = y ( b ) = 0   и при выполнении условия q ( x ) 0   все собственные значения положительны λ n > 0  .
  • Собственные функции y n ( x )   образуют на [ a , b ]   ортогональную с весом ρ ( x )   систему { y n ( x ) }  :
a b y n ( x ) y m ( x ) ρ ( x ) d x = 0 , n m .  

Численные методы решенияПравить

Применение к решению уравнений в частных производныхПравить

Задачи Штурма — Лиувилля возникают при решении уравнений в частных производных методом разделения переменных.

В качестве примера рассмотрим краевую задачу для уравнения гиперболического типа:

ρ ( x ) u t t = ( k ( x ) u x ) x q ( x ) u , 0 < x < l , t > 0 , ( 6 )  
( h 1 u x h u ) | x = 0 = 0 , ( H 1 u x + H u ) | x = l = 0 , ( 7 )  
u | t = 0 = Φ ( x ) , u t | t = 0 = Ψ ( x ) . ( 8 )  

Здесь x   и t   — независимые переменные, u ( x , t )   — неизвестная функция, ρ  , k  , q  , Φ  , Ψ   — известные функции, h ,   h 1 ,   H ,   H 1   — вещественные числа.[15] Будем искать не равные тождественно нулю частные решения уравнения (6), удовлетворяющие краевым условиям (7) в виде

u ( x , t ) = Y ( x ) T ( t ) . ( 9 )  

Подстановка вида (9) в уравнение (6) дает

( k ( x ) Y ( x ) ) q ( x ) Y ( x ) ρ ( x ) Y ( x ) = T ( t ) T ( t ) .  

Так как x   и t   — независимые переменные, то равенство возможно только если обе дроби равны константе. Обозначим эту константу через λ  . Получаем

T ( t ) + λ T ( t ) = 0 , ( 10 )  
( k ( x ) Y ( x ) ) + q ( x ) Y ( x ) = λ ρ ( x ) Y ( x ) , 0 < x < l . ( 11 )  

Подстановка вида (9) в краевые условия (7) дает

h 1 Y ( 0 ) h Y ( 0 ) = 0 , H 1 Y ( l ) + H Y ( l ) = 0. ( 12 )  

Нетривиальные решения (6) — (7) вида (9) существуют только при значениях λ  , являющихся собственными значениями задачи Штурма — Лиувилля (11) — (12) λ n  . Эти решения имеют вид T n ( t ) Y n ( x )  , где Y n ( x )   — собственные функции задачи (11) — (12), T n ( t )   — решения уравнения (10) при λ = λ n  . Решение задачи (6) — (8) находится в виде суммы частных решений (ряда Фурье по собственным функциям задачи Штурма — Лиувилля Y n ( x )  ):

u ( x , t ) = n = 1 T n ( t ) Y n ( x ) .  

Обратные задачи Штурма — ЛиувилляПравить

Обратные задачи Штурма — Лиувилля состоят в восстановлении потенциала q ( x )   оператора Штурма — Лиувилля y + q ( x ) y   и коэффициентов краевых условий по спектральным характеристикам.[8][3][4] Обратные задачи Штурма — Лиувилля и их обобщения имеют приложения в механике, физике, электронике, геофизике, метеорологии и других областях естествознания и техники. Существует важный метод интегрирования нелинейных эволюционных уравнений (например, уравнения КдФ), связанный с использованием обратной задачи Штурма — Лиувилля на оси ( < x <  ).

Одного спектра (множества собственных значений) как правило недостаточно для того, чтобы однозначно восстановить оператор. Поэтому в качестве исходных данных обратной задачи обычно используют следующие спектральные характеристики:

  1. Два спектра, соответствующие разным краевым условиям (задача Борга).
  2. Спектральные данные, включающие в себя собственные значения и весовые числа, равные квадратам норм собственных функций в пространстве L 2  .
  3. Функцию Вейля — мероморфную функцию, равную отношению двух характеристических функций разных краевых задач.

Каждый из наборов данных 1 — 3 однозначно определяет потенциал q ( x )  . Кроме того, задание функции Вейля равносильно заданию двух спектров или спектральных данных, поэтому обратные задачи по данным 1 — 3 эквивалентны. Существуют конструктивные методы решения обратных задач Штурма — Лиувилля, основанные на сведении нелинейных обратных задач к линейным уравнениям в некоторых банаховых пространствах.[4]

См. такжеПравить

ПримечанияПравить

  1. Левитан, Саргсян, 1988, с. 10.
  2. Юрко, 2010, с. 45.
  3. 1 2 3 Марченко, 1977.
  4. 1 2 3 4 Юрко, 2007.
  5. Наймарк, 1969, с. 72.
  6. Наймарк, 1969.
  7. Юрко, 2010, с. 25.
  8. 1 2 Левитан, Саргсян, 1988.
  9. Калиткин, 1978, с. 281.
  10. Калиткин, 1978, с. 284.
  11. Бахвалов Н. С., Жидков Н. П., Кобельков Г. М. Численные методы. — Бином, 2008. — ISBN 978-5-94774-815-4.
  12. Калиткин, 1978, с. 286.
  13. Калиткин, 1978, с. 287.
  14. Гельфанд И. М., Фомин С. В. Вариационное исчисление. — 1961.
  15. Юрко, 2010, с. 30.

ЛитератураПравить

  • Левитан Б. М., Саргсян И. С. Операторы Штурма — Лиувилля и Дирака. — М.: Наука. Гл. ред. физ.-мат. лит., 1988. — 432 с. — ISBN 5-02-013751-0.
  • Марченко В. А. Операторы Штурма — Лиувилля и их приложения. — Киев: Наукова думка, 1977.
  • Ахтямов А. М., Садовничий В. А., Султанаев Я. Т. Обратные задачи Штурма — Лиувилля с нераспадающимися краевыми условиями. — М.: Изд-во Московского университета, 2009. — 184 с. — ISBN 978-5-211-05557-5.
  • Юрко В. А. Уравнения математической физики. — Саратов, 2010.
  • Юрко В. А. Введение в теорию обратных спектральных задач. — М.: Физматлит, 2007. — 384 с. — ISBN 978-5-9221-07.
  • Наймарк М. А. Линейные дифференциальные операторы. — М.: Наука, 1969.
  • Калиткин Н. Н. Численные методы. — М.: Наука, 1978.