Sarcina de a fugi este Liouville

Problema lui Sturm - Liouville. numit după Jacques Charles François Sturm și Joseph Liouville. constă în găsirea de soluții nontrivială (adică distincte de identitatea zero) pe interval (a, b; b) Ecuațiile Sturm-Liouville







\ alpha _1 y '(a) + \ beta _1 y (a) = 0, \ qquad \ alpha ^ 2_1 + \ beta ^ 2_1 \ ne 0; \\ alfa _2 y '(b) + \ beta _2 y (b) = 0, \ qquad \ alpha ^ 2_2 + \ beta ^ 2_2 \ ne 0; \\ \ end și valorile parametrilor \ lambda, la care există astfel de soluții.

operator L [y] aici este funcția funcției y (x) operatorul diferențial liniar secundar al formularului

(operatorul Sturm-Liouville sau operatorul Schrodinger); x Este un adevărat argument.

funcții p (x), \; p '(x), \; q (x), \\ rho (x) se presupune că sunt continuu (a, b; b), în plus funcții p (x), \\ rho (x) pozitiv pentru (a, b; b).

Soluțiile nontriviale necesare se numesc funcțiile proprii ale acestei probleme și valorile \ lambda, sub care există o astfel de soluție - prin propriile sale valori (fiecare valoare proprie corespunde funcției proprii).

Formularea problemei

Forma ecuației

Dacă funcțiile \ rho și p sunt de două ori continuu diferențiate și pozitive pe interval [a, b] și funcția q continuu pe [a, b], apoi ecuația Sturm-Liouville a formei

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

cu ajutorul transformării Liouville este redusă la forma

Prin urmare, considerăm adesea ecuația Sturm-Liouville în forma (1), funcția q (x) numit potențial Template: Sfn Template: Sfn. Problemele Sturm-Liouville cu potențiale din diferite clase de funcții sunt studiate: continuu. L (Summable) L_2 și altele.

Tipuri de condiții limită

  • Condițiile Dirichlet y (a) = y (b) = 0.
  • Condiții Neumann y '(a) = y' (b) = 0
  • Condițiile lui Robin y '(a) - h y (a) = 0, \ quad y' (b) + H y (b) = 0.
  • Condiții mixte: condiții de diferite specii la capete diferite ale segmentului [a, b].
  • Condițiile limită de degradare ale formei generale
\ începe

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

  • Condiții periodice y (a) = y (b), \ quad y '(a) = y' (b).
  • Condiții antiperiodice y (a) = -y (b), \ quad y '(a) = -y' (b).
  • Condiții generale de frontieră
a_y (a) + a_y '(a) + a_y (b) + a_y' (b) = 0, \ quad i = 1, 2.

În ultimul caz, de obicei, impunem condiții suplimentare de regularitate asupra coeficienților a_. Format: Sfn Format: Sfn

Pentru confort, un segment arbitrar [a, b] adesea traduse într-un segment [0, l] sau [0, \ pi] prin modificarea variabilei.

Operatorul Sturm-Liouville

Operatorul Sturm-Liouville

L x = - \ frac \ Bigl (\ frac \ left [p (x) \ fracție \ dreapta] - q (x)

este un caz particular al unui operator diferențial liniar

Domeniul operatorului L constă din două rapoarte continuu diferențiate pe interval [a, b] funcție y, satisfăcând condițiile limită ale problemei Sturm-Liouville. Astfel, problema Sturm-Liouville poate fi privită ca o problemă asupra valorilor proprii și a funcțiilor proprii ale operatorului L: L y = \ lambda y. Dacă funcțiile p, q, \ rho iar coeficienții condițiilor limită sunt reali. operator L se auto-adăpostește într-un spațiu Hilbert L_2 ([a, b], \ rho (x) \, dx). În consecință, valorile proprii sunt reale și funcțiile proprii sunt ortogonale cu greutate \ rho (x) .

Rezolvarea problemei

Soluția problemei Sturm-Liouville cu potențial zero:

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

y (0) = y (1) = 0 poate fi găsit într-o formă explicită. lăsa \ lambda = \ rho ^ 2. Soluția generală a ecuației (2) pentru fiecare măsură fixă \ lambda are forma

y (x) = A \ frac + B \ cos \ rho x \ qquad (3) (în special cu \ rho = 0 (3) y (x) = ax + B). din y (0) = 0 ar trebui să fie B = 0. Înlocuindu-se (3) în condiția limită y (l) = 0, avem A \ frac = 0. Din moment ce căutăm soluții non-triviale, atunci A \ ne 0, și ajungem la ecuația eigenvalue

\ frac = 0. Rădăcinile sale \ rho_n = \ frac, în consecință, valorile proprii dorite au forma







\ lambda_n = \ stânga (\ frac \ right) ^ 2, \ quad n = 1, 2, 3, \ dots Funcțiile personale corespunzătoare sunt

y_n (x) = \ sin \ fracx, \ quad n = 1, 2, 3, \ puncte (într-un multiplicator constant).

Caz general

În cazul general, orice soluție a ecuației Sturm-Liouville

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

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

deciziile sale S (x, \ lambda) și C (x, \ lambda), satisfacerea condițiilor inițiale

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

soluţii S (x, \ lambda) și C (x, \ lambda) formează un sistem fundamental de soluții de ecuație (4) și sunt funcții întregi în \ lambda pentru fiecare fix x. (Când q (x) \ equiv 0 S (x, \ lambda) = \ sin \ rho x, C (x, \ lambda) = \ cos \ rho x, \ rho = \ sqrt \ lambda). Înlocuind (5) condițiile de graniță y (0) = y (\ pi) = 0, găsim că valorile proprii coincid cu zerourile funcției caracteristice

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

În cazul general, valorile proprii și funcțiile proprii nu pot fi găsite în mod explicit, ci se obțin formule asimptotice pentru ele:

\ sqrac \ lambda_n = n + \ frac + O \ stânga (\ frac \ right), \ quad c = \ frac \ int_0 ^ q (\ tau) y_n (x) = \ sin n x + O \ stânga (\ frac \ dreapta),

(în cazul unui proces continuu [0, \ pi] potențial q (x)) Model: Sfn Pentru mari n valorile proprii și funcțiile proprii sunt apropiate de valorile proprii și de funcțiile proprii ale problemei din exemplul cu potențial zero.

Proprietățile valorilor proprii și ale funcțiilor proprii

  • Există un set infinit de numere proprii: \ lambda_1 <\lambda_2 <\dots <\lambda_n <\dots.
  • La fiecare valoare proprie \ lambda_n corespunde unei funcții proprii unice, până la un factor constant y_n.
  • Toate valorile proprii sunt reale.
  • În cazul condițiilor limită y (a) = y (b) = 0 și dacă este condiție q (x) \ geqslant 0 toate valorile proprii sunt pozitive \ lambda_n> 0.
  • Funcții personalizate y_n (x) forma pe [a, \; b] ortogonal cu greutate \ rho (x) sistemul \:
\ int \ limits_a ^ b y_n (x) y_m (x) \ rho (x) \, dx = 0, \ quad n \ neq m.
  • Teorema lui Steklov deține.

Metode numerice de rezolvare

  • Metoda de fotografiere. Pentru a rezolva problema Sturm-Liouville cu condițiile de graniță Dirichlet y (a) = y (b) = 0, putem lua pentru ecuația inițială problema Cauchy cu condițiile inițiale u (a) = 0, u '(a) = 1 și conduceți ajustarea parametrilor \ lambda înainte de executarea condiției de limită dreapta. Model: Sfn
  • Metoda diferențelor finite Format: Sfn [1]. Se construiește o aproximație cu diferență finită care ne permite să înlocuim problema Sturm-Liouville prin găsirea valorilor proprii ale matricei.
  • Metoda vectorială completă. Diferența de funcționare y = \ este completat de un y_ = \ lambda. În ceea ce privește vectorul augmentat, obținem un sistem neliniar care poate fi rezolvat prin metoda Newton. Modelul: Sfn
  • Metoda lui Galerkin Model: Sfn
  • Metode variate. [2]

Aplicarea la soluția ecuațiilor diferențiale parțiale

De exemplu, considerăm problema de valoare limită pentru o ecuație de tip hiperbolic.

\ r (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)

aici x și T Sunt variabile independente. u (x, t) Este o funcție necunoscută, \ rho, k, q, \ Phi, \ Psi - funcții cunoscute, h, h_1, H, H_1 (6), care îndeplinesc condițiile de graniță (7) sub formă

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

Substituția formei (9) în ecuația (6) dă

deoarece x și T - variabile independente, atunci egalitatea este posibilă numai dacă ambele fracții sunt egale cu o constantă. Denumiți această constantă prin -\ lambda. Avem

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

Înlocuirea formei (9) în condițiile de graniță (7) dă

(1) + H Y (l) = 0. \ qquad (12)

Soluțiile nontriviale (6) - (7) ale formularului (9) există numai pentru valori \ lambda, care sunt valorile proprii ale problemei Sturm-Liouville (11) - (12) \ lambda_n. Aceste soluții au forma T_n (t) Y_n (x), unde Y_n (x) Sunt funcțiile proprii ale problemelor (11) - (12), T_n (t) - soluții ale ecuației (10) cu \ lambda = \ lambda_n. Soluția problemei (6) - (8) este sub forma unei sume de soluții particulare (seria Fourier în funcțiile proprii ale problemei Sturm-Liouville Y_n (x)):

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

Probleme inverse ale Sturm-Liouville

Problemele inverse de la Sturm-Liouville constau în restabilirea potențialului q (x) operatorul Sturm-Liouville -y + q (x) y și coeficienții condițiilor de graniță prin caracteristicile spectrale Modelul Sfn Template Sfn Problemele inverse ale lui Sturm-Liouville și generalizările lor au aplicații în mecanică. fizica. electronica. geofizica. Meteorologie și alte domenii ale științei și tehnologiei. Există o metodă importantă de integrare a ecuațiilor evoluției neliniare (de exemplu, ecuația KdV) asociată cu utilizarea problemei inverse Sturm-Liouville pe axă -\ infty ).

Un singur spectru (un set de valori proprii) nu este, de obicei, suficient pentru recuperarea unică a operatorului. Prin urmare, ca date inițiale ale problemei inverse, se folosesc de obicei următoarele caracteristici spectrale:

  1. Două spectre care corespund unor condiții limită diferite (problema Borg).
  2. Datele spectrale incluzând valori proprii și numere de greutate egale cu pătratele normelor funcțiilor proprii în spațiu L_2.
  3. Funcția Weyl este o funcție meromorfă. egală cu raportul dintre cele două funcții caracteristice ale diferitelor probleme de valoare a limitelor.

Fiecare set de date 1-3 identifică în mod unic potențialul q (x). În plus, specificarea funcției Weil este echivalentă cu specificarea a două date spectrale sau spectrale, astfel încât problemele inverse ale datelor 1-3 sunt echivalente. Există metode constructive pentru rezolvarea problemelor inverse Sturm-Liouville bazate pe reducerea problemelor inverse neliniare la ecuațiile liniare în anumite spații Banach.

notițe

literatură







Trimiteți-le prietenilor: