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)

on les funcions estan prescrites, i en el cas més simple són contínues en un interval finit tancat . El problema generalment ve formulat amb condicions de frontera, és a dir, valors específics de i/o en els extrems . La funció é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 , 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-LiouvilleModifica

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

(2) 

(3) 

on   és diferenciable, les funcions   són contínues i les funcions   són positives sobre l'interval  , i els valors   estan en l'interval   la teoria ens indica que

  • Els valors propis   del problema de S-L, són valors reals i ben ordenats en el sentit de que  .
  • A cada valor propi   li correspon una única funció pròpia   i   té exactament   zeros a la frontera  .
  • Les funcions pròpies són mútuament ortogonals i satisfan la relació d'ortogonalitat

(4) 

on   és la funció de pes.

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

(5) 

on   és la Delta de kronecker.

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-LiouvilleModifica

L'equació diferencial

 

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.

ExemplesModifica

 

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

 

 

pot ser transformada fàcilment en una forma de Sturm-Lioville, ja que  ; així l'equació de Legendre equivalent és:

 

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

 

Si dividim per   tenim:

 

Multiplicand per un factor integrant:

 

ens dóna

 

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

 

que és equivalent a dir:

 

  • En general, donada una equació diferencial

 

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

 

Operadors diferencials Sturm-LiouvilleModifica

L'operador lineal:

(6) 

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

(7) 

Aquest és precisament un problema de valors propis; on es tracta de trobar valors propis λ; i vectors propis   de l'operador  . Tanmateix, també s'ha d'incloure les condicions de frontera. Com a exemple es dirà que avaluarem el problema en l'interval   i es pondra les condicions de frontera  .

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

 

en l'interval  

 

en 0 i 1.

Aquí   és la funció en l'espai  . Si una solució   existeix i és única, la hi pot escriure de la forma:

 

perquè la transformació de   a   ha de ser lineal. Ara s'observa que el trobar els vectors propis i els valors propis de   és essencialment el mateix que trobar els vectors i valors propis de  . Efectivament, si   és un vector propi de   amb valors propis   ha d'existir un   que també és vector propi de   amb valors propis  .


ExempleModifica

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

(8) 

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

 

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

 

és una solució amb eigenvalor λ; = −k 2. 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

 

amb   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

 

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 2, 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

 

Podríem haver trobat la solució amb antidiferenciació. Aquesta tècnica dóna u=(x32x)/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ó normalsModifica

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; xL 1, 0 &# 8804; iL 2. El desplaçament vertical W (x' y' t) de la membrana és governada per l'equació d'ona:

 

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/ c2) 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 "= X on   és constant. Les condicions de frontera ("sostinguda en un marc rectangular") són W =0 when x =0, L 1 o i =0, L 2 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

 

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

 

Les funcions Wmn 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  . Aquesta representació pot requerir una suma infinita convergent.

Operadors de Sturm-Liouville com operadors HermíticsModifica

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

 

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èricModifica

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ó

 

é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 λ; 0* (sovint λ; 0*=0; no ha de ser un eigenvalor) i amb qualsevol solució i 0 de (1) amb λ; =λ 0* 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 0 i λ; 0* apropiats.) Dues successions de funcions X(n) (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 02) i w'y 02 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,   per a n parell,

(10)  per a n senar,   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>

  i   .

Llavors per a qualsevol λ; (real o complex), o 0 i o 1 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 0.)

Ara s'escullen coeficients c 0, c 1 de manera que la combinació i =c 0o 0 + c 1o 1 satisfaci la primera condició de frontera (2). Això és fàcil de fer perquè X(n) (a) =0 i X~(n) (a) =0, per a n >0. Els valors de X(n) (b) i X~(n) (b) donen els valors de o 0(b) i o 1(b) i de les derivades o 0'(b) i o 1'(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 λ; = λ 0, 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·laModifica

El mateix mètode Spps pot usar-se per a trobar una solució inicial i 0. 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=0. Mentre no hi ha grantía que u0 o u1 no s'anul·li, la funció complexa y0= u0+ i' ' o 1 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 0 de (1) per al valor λ; 0=0. A la pràctica, si (1) té coeficients reals, les solucions basades en i 0 tindran parts imaginàries molt petites que s'hauran de descartar.

Vegeu tambéModifica

ReferènciesModifica

  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.
  • 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.