Задача Штурма - Ліувілля

Задача Штурма - Ліувілля. названа на честь Жака Шарля Франсуа Штурма та Жозефа Ліувілля. полягає у знаходженні нетривіальних (тобто відмінних від тотожного нуля) рішень на проміжку (A, \; b) рівняння Штурма - Ліувілля

\ Alpha _1 y '(a) + \ beta _1 y (a) = 0, \ qquad \ alpha ^ 2_1 + \ beta ^ 2_1 \ ne 0; \\ \ alpha _2 y '(b) + \ beta _2 y (b) = 0, \ qquad \ alpha ^ 2_2 + \ beta ^ 2_2 \ ne 0; \\ \ end і значень параметра \ lambda, при яких такі рішення існують.

оператор L [y] тут - це діючий на функцію y (x) лінійний диференційний оператор другого порядку виду

(Оператор Штурма - Ліувілля або оператор Шредінгера), x - речовинний аргумент.

функції p (x), \; p '(x), \; q (x), \; \ rho (x) передбачаються безперервними на (A, \; b), крім того функції p (x), \; \ rho (x) позитивні на (A, \; b).

Шукані нетривіальні рішення називаються власними функціями цього завдання, а значення \ lambda, при яких таке рішення існує - її власними значеннями (кожному власному значенню відповідає власна функція).

Постановка задачі

вид рівняння

якщо функції \ rho і p двічі безперервно діфференцируєми і позитивні на відрізку [A, b] і функція q неперервна на [A, b], то рівняння Штурма - Ліувілля виду

(P (x) y ')' - q (x) y + \ lambda \ rho (x) y = 0

за допомогою перетворення Ліувілля приводиться до вигляду

Тому часто розглядають рівняння Штурма - Ліувілля в вигляді (1), функцію q (x) називають потенціалом Шаблон: Sfn Шаблон: Sfn. Вивчаються завдання Штурма - Ліувілля з потенціалами з різних класів функцій: безперервними. L (Сумовною), L_2 та інших.

Види крайових умов

  • умови Діріхле y (a) = y (b) = 0.
  • умови Неймана y '(a) = y' (b) = 0
  • умови Робена y '(a) - h y (a) = 0, \ quad y' (b) + H y (b) = 0.
  • Змішані умови: умови різних видів в різних кінцях відрізка [A, b].
  • Розпадаються крайові умови загального вигляду
\ begin

\ Alpha _1 y '(a) + \ beta _1 y (a) = 0, \ qquad \ alpha ^ 2_1 + \ beta ^ 2_1 \ ne 0; \\ \ alpha _2 y '(b) + \ beta _2 y (b) = 0, \ qquad \ alpha ^ 2_2 + \ beta ^ 2_2 \ ne 0. \\ \ end

  • періодичні умови y (a) = y (b), \ quad y '(a) = y' (b).
  • Антіперіодіческіе умови y (a) = -y (b), \ quad y '(a) = -y' (b).
  • Загальні крайові умови
a_ y (a) + a_ y '(a) + a_ y (b) + a_ y' (b) = 0, \ quad i = 1, 2.

В останньому випадку зазвичай накладаються додаткові умови регулярності на коефіцієнти a_. Шаблон: Sfn Шаблон: Sfn

Для зручності довільний відрізок [A, b] часто переводять в відрізок [0, l] або [0, \ pi] за допомогою заміни змінної.

Оператор Штурма - Ліувілля

Оператор Штурма - Ліувілля

L y = - \ frac \ Bigl (\ frac \ left [p (x) \ fracy \ right] - q (x) y \ Bigr)

являє собою окремий випадок лінійного диференціального оператора

Область визначення оператора L складається з двічі безперервно диференційовних на відрізку [A, b] функцій y, задовольняють крайовим умовам задачі Штурма - Ліувілля. Таким чином, завдання Штурма - Ліувілля можна розглядати як задачу на власні значення і власні функції оператора L: L y = \ lambda y. якщо функції p, q, \ rho і коефіцієнти крайових умов речові. оператор L є самосопряженним в гільбертовому просторі L_2 ([a, b], \ rho (x) \, dx). Отже, його власні значення речовинні і власні функції ортогональні з вагою \ Rho (x) .

Рішення завдання

Рішення задачі Штурма - Ліувілля з нульовим потенціалом:

-y = \ lambda y, \ qquad (2)

y (0) = y (l) = 0 може бути знайдено в явному відеШаблон: Sfn. нехай \ Lambda = \ rho ^ 2. Загальне рішення рівняння (2) при кожному фіксованому \ lambda має вигляд

y (x) = A \ frac + B \ cos \ rho x \ qquad (3) (Зокрема, при \ Rho = 0 (3) дає y (x) = Ax + B). з y (0) = 0 слід B = 0. Підставляючи (3) в крайове умова y (l) = 0, отримуємо A \ frac = 0. Так як ми шукаємо нетривіальні рішення, то A \ ne 0, і ми приходимо до рівняння на власні значення

\ Frac = 0. його коріння \ Rho_n = \ frac, отже, шукані власні значення мають вид

\ Lambda_n = \ left (\ frac \ right) ^ 2, \ quad n = 1, 2, 3, \ dots а відповідні їм власні функції суть

y_n (x) = \ sin \ fracx, \ quad n = 1, 2, 3, \ dots (З точністю до постійного множника).

Загальний випадок

У загальному випадку будь-яке рішення рівняння Штурма - Ліувілля

-y + q (x) y = \ lambda y \ qquad (4)

y (x) = A S (x, \ lambda) + B C (x, \ lambda) \ qquad (5)

його рішень S (x, \ lambda) і C (x, \ lambda), задовольняють початковим умовам

S (0, \ lambda) = C '(0, \ lambda) = 0, \ quad S' (0, \ lambda) = C (0, \ lambda) = 1.

рішення S (x, \ lambda) і C (x, \ lambda) утворюють фундаментальну систему рішень рівняння (4) і є цілими функціями по \ lambda при кожному фіксованому x. (При q (x) \ equiv 0 S (x, \ lambda) = \ sin \ rho x, C (x, \ lambda) = \ cos \ rho x, \ Rho = \ sqrt \ lambda). Підставляючи (5) в крайові умови y (0) = y (\ pi) = 0, отримуємо, що власні значення збігаються з нулями характеристичної функції

\ Delta (\ lambda) = S (\ pi, \ lambda),

У загальному випадку власні значення і власні функції не можуть бути знайдені в явному вигляді, проте для них отримані асимптотичні формули:

\ Sqrt \ lambda_n = n + \ frac + O \ left (\ frac \ right), \ quad c = \ frac \ int_0 ^ q (\ tau) \, d \ tau, y_n (x) = \ sin n x + O \ left (\ frac \ right),

(В разі безперервного на [0, \ pi] потенціалу q (x)) .Шаблон: Sfn При великих n власні значення і власні функції близькі до власних значень і власних функцій задачі з прикладу з нульовим потенціалом.

Властивості власних значень і власних функцій

  • Існує нескінченна рахункове безліч власних значень: \ lambda_1 <\lambda_2 <\dots <\lambda_n <\dots.
  • Кожному власному значенню \ lambda_n відповідає єдина з точністю до постійного множника власна функція y_n.
  • Всі власні значення речовинні.
  • У разі граничних умов y (a) = y (b) = 0 і при виконанні умови q (x) \ geqslant 0 всі власні значення позитивні \ Lambda_n> 0.
  • власні функції y_n (x) утворюють на [A, \; b] ортогональную з вагою \ Rho (x) систему \:
\ Int \ limits_a ^ b y_n (x) y_m (x) \ rho (x) \, dx = 0, \ quad n \ neq m.
  • Має місце теорема Стеклова.

Чисельні методи рішення

  • Метод стрільби. Щоб вирішити задачу Штурма - Ліувілля з крайовими умовами Діріхле y (a) = y (b) = 0, можна взяти для вихідного рівняння задачу Коші з початковими умовами u (a) = 0, u '(a) = 1 і вести пристрілювання параметра \ lambda до виконання правого крайового условія.Шаблон: Sfn
  • Метод кінцевих різниць Шаблон: Sfn [1]. Будується кінцево-різницева апроксимація, яка дозволяє замінити задачу Штурма - Ліувілля завданням знаходження власних значень матриці.
  • Метод доповненого вектора. Разностная власна функція y = \ доповнюється компонентом y_ = \ lambda. Щодо доповненого вектора виходить нелінійна система, яка може бути вирішена методом Ньютона .Шаблон: Sfn
  • Метод Гальоркіна .Шаблон: Sfn
  • Варіаційні методи. [2]

Застосування до вирішення рівнянь в приватних похідних

Як приклад розглянемо крайову задачу для рівняння гіперболічного типу.

\ Rho (x) u_ = (k (x) u_x) _x - q (x) u, \ quad 0 0, \ qquad (6) (H_1 u_x - h u) _ = 0, \ quad (H_1 u_x + H u) _ = 0, \ qquad (7) u_ = \ Phi (x), \ quad u_ = \ Psi (x). \ Qquad (8)

тут x і t - незалежні змінні. u (x, t) - невідома функція, \ rho, k, q, \ Phi, \ Psi - відомі функції, h, h_1, H, H_1 - речові числа .Шаблон: Sfn Будемо шукати нерівні тотожне нулю приватні рішення рівняння (6), що задовольняють крайовим умовам (7) у вигляді

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

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

Так як x і t - незалежні змінні, то рівність можливо тільки якщо обидві дробу рівні константі. Позначимо цю константу через -\ lambda. отримуємо

T (t) + \ lambda T (t) = 0, \ qquad (10) -(K (x) Y '(x))' + q (x) Y (x) = \ lambda \ rho (x) Y (x), \ quad 0

Підстановка виду (9) в крайові умови (7) дає

h_1 Y '(0) - h Y (0) = 0, \ quad H_1 Y' (l) + H Y (l) = 0. \ qquad (12)

Нетривіальні рішення (6) - (7) виду (9) існують тільки при значеннях \ lambda, є власними значеннями завдання Штурма - Ліувілля (11) - (12) \ lambda_n. Ці рішення мають вигляд T_n (t) Y_n (x), де Y_n (x) - власні функції задачі (11) - (12), T_n (t) - рішення рівняння (10) при \ Lambda = \ lambda_n. Рішення завдання (6) - (8) знаходиться у вигляді суми приватних рішень (ряду Фур'є за власними функціями задачі Штурма - Ліувілля Y_n (x)):

u (x, t) = \ sum_ ^ T_n (t) Y_n (x).

Зворотні задачі Штурма - Ліувілля

Зворотні задачі Штурма - Ліувілля складаються у відновленні потенціалу q (x) оператора Штурма - Ліувілля -y + q (x) y і коефіцієнтів крайових умов по спектральним характерістікам.Шаблон: Sfn Шаблон: Sfn Шаблон: Sfn Зворотні задачі Штурма - Ліувілля і їх узагальнення мають додатки в механіці. фізики. електроніці. геофізики. метеорології та інших областях природознавства і техніки. Існує важливий метод інтегрування нелінійних еволюційних рівнянь (наприклад, рівняння КдФ), пов'язаний з використанням оберненої задачі Штурма - Ліувілля на осі ( -\ infty ).

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

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

Кожен з наборів даних 1-3 однозначно визначає потенціал q (x). Крім того, завдання функції Вейля рівносильно завданням двох спектрів або спектральних даних, тому зворотні завдання за даними 1-3 еквівалентні. Існують конструктивні методи розв'язання обернених задач Штурма - Ліувілля, засновані на зведенні нелінійних обернених задач до лінійних рівнянь в деяких банахових просторах .Шаблон: Sfn

Примітки

література