Usuari:Mcapdevila/temp

En matemàtiques, una equació de Sturm-Liouville, que pren el seu nom de Jacques Xerris François Sturm (1803-1855) i Joseph Liouville (1809-1882), és una equació diferencial lineal de segon ordre de la forma

(1) $-{\frac {d}{dx}}\left[p(x){\frac {dy}{dx}}\right]+q(x)y=\lambda w(x)y$

on les funcions $\ p(x),q(x),w(x)$ estan prescrites, i en el cas més simple són contínues en un interval finit tancat $[a,b]$ . El problema generalment ve formulat amb condicions de frontera, és a dir, valors específics de $y$ i/o ${\frac {dy}{dx}}$ en els extrems $a,b$ . La funció $w(x)$ és trucada funció de densitat o funció de pes.

El valor de λ; no s'especifica en l'equació; el trobar els valors λ; on hi hagi una solució no trivial de l'equació que satisfaci condicions de frontera es denomina el problema de Sturm-Liouville (S-L).

Tals valors de λ; són cridats valors propis o eigenvalores del problema de S-L que planteja (1) conjuntament amb les condicions de frontera. Les solucions corresponents són les funcions pròpies o les eigenfunciones o els eigenvectores del problema. Sota suposicions normals en els coeficients de les funcions $p(x),q(x),w(x)$ , aquestes indueixen operadors diferencials hermítics en algunes funcions definides per les condicions de frontera. La teoria resultant de l'existència i el comportament asimptòtic dels valors propis, la teoria qualitativa corresponent de les funcions pròpies i les seves funcions adequades completes es coneix com a teoria de Sturm-Liouville. Aquesta teoria és important en matemàtica aplicada, on els problemes S-L ocorren molt comunament, particularment en resoldre equacions diferencials parcials amb separació de variables.

Teoria de Sturm-Liouville

Quan les condicions de frontera són regulars de la forma

(2) $y(a)\cos \alpha -p(a)y^{\prime }(a)\sin \alpha =0$

(3) $y(b)\cos \beta -p(b)y^{\prime }(b)\sin \beta =0$

on $p(x)$ és diferenciable, les funcions $p(x),q(x),w(x)$ són contínues i les funcions $\ p(x),w(x)$ són positives sobre l'interval $[a,b]$ , i els valors $\alpha ,\beta$ estan en l'interval $[0,\pi [$ la teoria ens indica que

Els valors propis $\lambda _{n}$ del problema de S-L, són valors reals i ben ordenats en el sentit de que $\lambda _{1}<\lambda _{2}<\lambda _{3}<\cdots <\lambda _{n}<\cdots \to \infty ;$ .
A cada valor propi $\lambda _{n}$ li correspon una única funció pròpia $\ y_{n}(x)$ i $\ y_{n}(x)$ té exactament $\ n-1$ zeros a la frontera $\ (a,b)$ .
Les funcions pròpies són mútuament ortogonals i satisfan la relació d'ortogonalitat

(4) $\int _{a}^{b}y_{n}(x)y_{m}(x)w(x)\,dx=0,m\neq n,$

on $\ w(x)$ és la funció de pes.

Un conjunt ortonormal pot ser format si el conjunt de funcions pròpies satisfà la relació d'ortogonalitat

(5) $\int _{a}^{b}y_{n}(x)y_{m}(x)w(x)\,dx=\delta _{mn},$

on $\ \delta _{mn}$ és la Delta de kronecker.

Els valors propis del problema de S-L poden ser caracteritzats pel Quocient de rayleigh

No s'ha pogut entendre (error de sintaxi): {\displaystyle \lambda_n = \frac{-p(x) y_{n}(x) i ' _{n}(x)|_a^b + \int_a^b [pàg. y'_{n}(x)^2 + q y_{n}(x)^2]\,dx}{\int_a^b y_{n}(x)^2 w(x)\, dx}}

.

Forma de Sturm-Liouville

L'equació diferencial

$-{d \over dx}\left[p(x){d \over dx}y(x)\right]+q(x)y(x)=\lambda w(x)y(x)$

es diu que és de la forma de S-L o de la forma autoadjunta. Tota equació diferencial ordinària lineal de segon ordre pot portada a aquesta forma en multiplicar-li per un factor integrant apropiat.

Exemples

L'Equació de bessel:

$\ x^{2}y''+xy'+(\lambda ^{2}x^{2}-\nu ^{2})y=0$

pot ser escrita en la forma de S-L així:

$(xy')'+(\lambda ^{2}x-\nu ^{2}/x)y=0.\,$

L'Equació de legendre:

$(1-x^{2})y''-2xy'+\nu (\nu +1)y=0\;\!$

pot ser transformada fàcilment en una forma de Sturm-Lioville, ja que

\scriptstyle (1-x^{2})'=-2x

; així l'equació de Legendre equivalent és:

$[(1-x^{2})y']'+\nu (\nu +1)y=0\;\!$

Un altre exemple simple és una equació diferencial de la forma:

$\ x^{3}y'''-xy'+2y=0$

Si dividim per

\ x^{3}

tenim:

$y'''-{x \over x^{3}}y'+{2 \over x^{3}}y=0$

Multiplicand per un factor integrant:

$e^{\int -{x/x^{3}}\,dx}=e^{\int -{1/x^{2}}\,dx}=e^{1/x},$

ens dóna

$e^{1/x}y'''-{e^{1/x} \over x^{2}}y'+{2e^{1/x} \over x^{3}}i=0$

que pot posar-se fàcilment en la forma de S-L així:

$De^{1/x}=-{e^{1/x} \over x^{2}}$

que és equivalent a dir:

$(e^{1/x}y')'+{2e^{1/x} \over x^{3}}i=0.$

En general, donada una equació diferencial

$P(x)y''+Q(x)y'+R(x)y=0\,$

dividida per a

\ P(x)

, multiplicada per un factor integrant tenim la forma de S-L:

$e^{\int {Q(x)/P(x)}\,dx},$

Operadors diferencials Sturm-Liouville

L'operador lineal:

(6) $Lo=-{d \over dx}\left[p(x){du \over dx}\right]+q(x)u$

pot ser vista com la transformació d'una funció $\ u$ en una altra funció $\ Lu$ . Es pot estudiar aquest operador lineal en el context de l'anàlisi funcional. Si posem $\ w=1$ en l'equació (1), podem escriure-la:

(7) $Lo=\lambda o\,.$

Aquest és precisament un problema de valors propis; on es tracta de trobar valors propis λ; i vectors propis $\ u$ de l'operador $\ L$ . Tanmateix, també s'ha d'incloure les condicions de frontera. Com a exemple es dirà que avaluarem el problema en l'interval $\ [0,1]$ i es pondra les condicions de frontera $\ u(0)=u(1)=0$ .

La importància de problemes de valors propis aquesta en el fet que ens ajuda a resoldre problemes associats inhomogenis:

$Lo=f\,$

en l'interval $\ [0,1]$

$u=0\,$

en 0 i 1.

Aquí $\ f$ és la funció en l'espai $\ L^{2}$ . Si una solució $\ u$ existeix i és única, la hi pot escriure de la forma:

$u=Af\,$

perquè la transformació de $\ f$ a $\ u$ ha de ser lineal. Ara s'observa que el trobar els vectors propis i els valors propis de $\ A$ és essencialment el mateix que trobar els vectors i valors propis de $\ L$ . Efectivament, si $\ u$ és un vector propi de $\ L$ amb valors propis $\lambda$ ha d'existir un $\ u$ que també és vector propi de $\ A$ amb valors propis ${\frac {1}{u}}$ .

Exemple

Es busca una funció o (x) que resolgui el següent problema S-L:

(8) $Lo={\frac {d^{2}u}{dx^{2}}}=\lambda u$

en on les incògnites són λ; and o (x) . Prenguem per exemple les condicions de frontera

u(0)=u(\pi )=0\,

i observem que si k és qualsevol enter, llavors la funció

u(x)=\operatorname {sen} kx\,

és una solució amb eigenvalor λ; = −k ². Sabem que les solucions d'un problema S-L formen una base ortogonal, i de la teoria de les sèries de Fourier sabem que aquest conjunt de funcions sinusoïdals és una base ortoganal. Ja que les bases otorgonales per definició són màximes, concloem que aquest problema S-L no té més eigenvectores.

Amb base en això resolguem el problema inhomogeni

Lu=x\,

amb $x\in (0,\pi )$ i amb les mateixes condicions de frontera. En aquest cas hem de posar f (x) = x en formes de sèrie de Fourier. El lector pot verificar, sigui integrand exp(ikx) x dx o consultant una taula de transformades de Fourier, que

Lo=\sum _{k=1}^{\infty }-2{\frac {(-1)^{k}}{k}}\sin kx.

Aquesta Sèrie de fourier és problemàtica per les seves males propietats de convergència: no està clar a priori si convergeix puntualment. Pel fet (de la teoria d'anàlisi de Fourier) que els coefientes són quadrat-sumables, la sèrie convergeix en L ², amb el qual n'hi ha prou per a la present discussió. Esmentem per al lector interessat que podem aplicar un resultat que diu que la sèrie de Fourier convergeix en cada punt de diferenciabilidad, i en els punts de salt (per exemple (la funció x, considerada com a funció periòdica, té un salt en n ) convergeix a la mitjana dels límits esquerre i dret.

Per tant, per la fórmula (8) obtenim que la solució és

u=\sum _{k=1}^{\infty }2{\frac {(-1)^{k}}{k^{3}}}\sin kx.

Podríem haver trobat la solució amb antidiferenciació. Aquesta tècnica dóna u=(x³-π²x)/6, la sèrie del qual de Fourier concorda amb la solució que trobem. La tècnica d'antidiferención generalment no és estri per a les equacions diferencials de diverses variables.

Aplicació a les maneres de vibració normals

Certes equacions en derivades parcials poden resoldre's amb l'ajuda de la teoria de S-L. Suposem que ens interessen les maneres de vibració d'una membrana prima, sostinguda en un marc rectangular, 0 &# 8804; x ≤L ₁, 0 &# 8804; i ≤L ₂. El desplaçament vertical W (x' y' t) de la membrana és governada per l'equació d'ona:

{\frac {\partial ^{2}W}{\partial x^{2}}}+{\frac {\partial ^{2}W}{\partial y^{2}}}={\frac {1}{c^{2}}}{\frac {\partial ^{2}W}{\partial t^{2}}}.

El mètode de separació de variables suggereix buscar primer solucions de la forma senzilla W = X(((( x)) ×; Y(i) ×; T(t) . Per a una tal funció Wl'equació en derivades parcials pot escriure's com X" / X+ Y"/ Y= (1/ c²) T" )) /T . Ja que els tres termes d'aquesta equació són funcions de x,i,t separadament, han de ser constants. Per exemple, el primer terme dóna X "= $\lambda$ X on $\lambda$ és constant. Les condicions de frontera ("sostinguda en un marc rectangular") són W =0 when x =0, L ₁ o i =0, L ₂ el qual defineix els problemes de S-L eigenvalue el més senzills possibles com en l'exemple i donen les "solucions de maneres normals" per a W amb dependència harmònica del temps

W_{mn}(x,y,t)=A_{mn}\sin \left({\frac {m\pi x}{L_{1}}}\right)\sin \left({\frac {n\pi y}{L_{2}}}\right)\cos \left(\omega _{mn}t\right)

on m,n són enters no nuls, A_mn són constants arbitràries i

\omega _{mn}^{2}=c^{2}\left({\frac {m^{2}\pi ^{2}}{L_{1}^{2}}}+{\frac {n^{2}\pi ^{2}}{L_{2}^{2}}}\right).

Les funcions W_mn formen una base de l'Espai d'hilbert de solucions (generalitzades) de l'equació d'ona; això és, una solució arbitrària W pot descompondre's com a suma d'aquestes maneres, que vibren a freqüències individuals $\omega _{mn}$ . Aquesta representació pot requerir una suma infinita convergent.

Operadors de Sturm-Liouville com operadors Hermítics

Moltes de les propietats dels operadors de S-L vénen del fet que aquests són operadors hermÍTIQUES respecte al producte intern:

$\langle o,v\rangle _{w}=\int _{a}^{b}w(x)u(x)v(x)dx.$

I així els valors propis dels operadors de S-L són reals i que les funcions pròpies corresponen a diferents valors propis són ortogonals.

Representació de Solucions i Càlcul Numèric

L'equació de S-L (1) amb condicions de frontera pot resoldre's a la pràctica per una varietat de mètodes numèrics. En casos difícils pot ser necessari dur a terme els càlculs intermedis amb diversos cents de dígits decimals per aconseguir els eigenvalores correctament a unes quantes xifres.

1. Mètodes de tret. ^[1]. ^[2]. Aquests mètodes funcionen endevinant un valor de λ;, resolent un problema en valors inicials definit per les condicions de frontera en un extrem, diguem a, de l'interval [a,b ], comparant el valor que aquesta solució pren en l'altre extrem b amb l'altra condicion de frontera desitjada i finalment augmentant o reduint λ; segun sigui necessari per corregir el valor original. Aquesta estratègia no és aplicable per trobar eigenvalores complexos.

2. Mètode de diferències finites.

3. El mètode de Sèries de Potències del Paràmetre Espectral < ref name="Kp " >. V. V. Kravchenko, R. M. Porter, "Spectral parameter power sèries for Sturm-Liouville problems, " Mathematical Methods in the Applied Sciences (Mmas) 33, 2010, 459-468 </ref>. (Spps per les seves sigles en anglès) aplica una generalització del seguiente fet sobre les equacions en derivades ordinàries de segon ordre: si i és una solució que no s'anul·la en cap punt de [a,b ], llavors la funció

y(x)\int _{a}^{x}{\frac {dt}{p(t)y(t)^{2}}}

és una solució de la mateixa equació i és linealment independent de i . A més, totes les solucions són combinacions lineals d'aquestes dues solucions. En l'algoritme Spps, un ha de començar amb un valor arbitrari λ; ₀^* (sovint λ; ₀^*=0; no ha de ser un eigenvalor) i amb qualsevol solució i ₀ de (1) amb λ; =λ ₀^* la qual no s'anul·li en [a,b ]. (Discussió #CONSTRUCCIÓN D'UNA SOLUCIÓ QUE no s'ANUL·LA|a baix sobre maneres de trobar i ₀ i λ; ₀^* apropiats.) Dues successions de funcions X⁽ⁿ⁾ (t), X ^~(n)(t) en [a,b ], que es cridaran integrals iterades, es defineixen de manera recursiva com segueix. Primer, quan n =0, es prenen iguals a 1 idènticament en [a,b]. Per obtenir les següents funcions es multipliquen alternativament per 1/ ( p' i ₀²) i w'y ₀² i després s'integren, específicament

(9)No s'ha pogut entendre (error de sintaxi): {\displaystyle X^{(n)}(t) = - \int_a^x X^{(n-1)}(t) p(t)^{-1} y_0(t)^{-y_0(t)^{-2}\,dt } per a n senar, $\int _{a}^{x}X^{(n-1)}(t)y_{0}(t)^{2}w(t)\,dt$ per a n parell,

(10) ${\tilde {X}}^{(n)}(t)=\int _{a}^{x}{\tilde {X}}^{(n-1)}(t)y_{0}(t)^{2}w(t)\,dt$ per a n senar, $-\int _{a}^{x}{\tilde {X}}^{(n-1)}(t)p(t)^{-1}y_{0}(t)^{-2}\,dt$ per a n parell,

quan n >0. Les integrals iterades així obtingudes s'apliquen ara com coeficients en les dues següents sèries de potències en λ; :

< dl>

u_{0}=y_{0}\sum _{k=0}^{\infty }(\lambda -\lambda _{0}^{*})^{k}{\tilde {X}}^{(2k)}

i

u_{1}=y_{0}\sum _{k=0}^{\infty }(\lambda -\lambda _{0}^{*})^{k}X^{(2k+1)}

.

Llavors per a qualsevol λ; (real o complex), o ₀ i o ₁ són solucions linealment dependents de l'equació (1) corresponent. (Les funcions pàg. (x) i q (x) participen en aquesta construcció de forma implícita per la seva influència en l'elecció de i ₀.)

Ara s'escullen coeficients c ₀, c ₁ de manera que la combinació i =c ₀o ₀ + c ₁o ₁ satisfaci la primera condició de frontera (2). Això és fàcil de fer perquè X⁽ⁿ⁾ (a) =0 i X^~(n) (a) =0, per a n >0. Els valors de X⁽ⁿ⁾ (b) i X^~(n) (b) donen els valors de o ₀(b) i o ₁(b) i de les derivades o ₀'(b) i o ₁'(b), llavors la segona condició de frontera (3) es converteix en una sèrie de potències en λ; . Per al treball numèric u pot truncar aquesta sèrie a un nombre finit de termes i obtenir un polinomi en λ; calculable, les arrels del qual són aproximacions dels eigenvalores que es busquen.

Quan λ; = λ ₀, això es redueix a la construcció original descrita a dalt per a una solució linealment independent a una donada. Les representacions (9),(10) també tenen aplicacions teòriques en la teoria de S-L. < ref name="Kp"/ >.

Construcció d'una solució que no s'anul·la

El mateix mètode Spps pot usar-se per a trobar una solució inicial i ₀. Es consideri l'equació (pi) ' = μ;'qy; i.e. q) ) w((, i λes reemplacen en (1) amb 0, - qi μrespectivament. Llavors la funció constant 1 és una solució que no s'anul·la, respecte a l'eigenvalor μ₀=0. Mentre no hi ha grantía que u₀ o u₁ no s'anul·li, la funció complexa y₀= u₀+ i' ' o ₁ mai s'anul·larà perquè dues solucions linealmentes d'una equació S-L regular no poden anul·lar-se simultàniament (com a conseqüència de la unicitat de solucions per a problemes de valors inicials). Aquest truc dóna una solució i ₀ de (1) per al valor λ; ₀=0. A la pràctica, si (1) té coeficients reals, les solucions basades en i ₀ tindran parts imaginàries molt petites que s'hauran de descartar.

Vegeu també

Referències

↑ . J. Sr. Pryce Numerical Solution of Sturm-Liouville Problems, Clarendon Press, Oxford, 1993.
↑ . V. Ledoux, M. Van Daele, G. Vanden Berghe, "Efficient computation of high index Sturm-Liouville eigenvalues for problems in physics, " Comput. Phys. Comm. 180, 2009, 532 ??554.

Pàg. Hartman Ordinary Differential Equations, Siam, Philadelphia, 2002 (2 nd edition). Isbn 978-0-89871-510-1
A. SR. Polyanin and V. F. Zaitsev, Handbook of Exact Solutions for Ordinary Differential Equations, Chapman & Hall/Crc Press, Boca Raton, 2003 (2 nd edition). Isbn 1-58488-297-2
G. Teschl, Ordinary Differential Equations and Dynamical Systems, http://www.mat.univie.ac.at/~gerald/ftp/book-ode/ (Chapter 5)
G. Teschl, Mathematical Methods in Quàntum Mechanics and Applications to Schrödinger Operators, http://www.mat.univie.ac.at/~gerald/ftp/book-schroe/ (see Chapter 9 for singular S-L operators and connections with quàntum mechanics)
A. Zettl, Sturm-Liouville Theory, American Mathematical Society, 2005. Isbn 0-8218-3905-5.

[1] . J. Sr. Pryce Numerical Solution of Sturm-Liouville Problems, Clarendon Press, Oxford, 1993.

[2] . V. Ledoux, M. Van Daele, G. Vanden Berghe, "Efficient computation of high index Sturm-Liouville eigenvalues for problems in physics, " Comput. Phys. Comm. 180, 2009, 532 ??554.

[1]

[2]