4
Separacin de variables para Laplace
Repaso de EDOs
Algunas EDOs de primer orden
dy d
= F en D acotado tiene solucin nica. [Si F 0 mximo y mnimo se dan en D ]. = en D = F en D Tiene unicidad salvo constante si se cumple D F ddy = D ds . n = en D + yy = 0 , = XT X + X = 0 , Y Y = 0 .
= (, y ) resolubles
+ 1 = 0 , = R + = 0 , r 2 R + rR R = 0 . rr + 1 r r r 2 En un crculo, n = cos n, sen n , Rn = {r n } (interior), Rn = {r n } (exterior), n = 0, 1, . . . 2 2 1 1 Si G(r,; ,) = 4 ln 2 + r 2 2r cos( ) 4 ln R2 + r R 2 2 r cos ( ) , la solucin de 2 r 2 2 R 2 = F (r, ) , r < R () d . es (r, ) = 0 0 G(r, ; , ) F ( , ) d d + R 2 0 R2 2Rr cos()+r 2
(R, ) = ()
[Si , y son continuas en un entorno de (o , yo ) existe solucin nica con y (o ) = yo ]. y p() y Separables: y = q(y) . y = , z = . y = ( + by ) , z = + by . Lineales: y = ()y + () y = Ce ()d + e ()d e ()d ()d [ sgh + yp ]. Exactas: M(, y )+ N(, y ) y = 0 con M = U , N = Uy My N U(, y ) = C . EDOs lineales de segundo orden [n] y + ()y + b()y = ()
, b , continuas en .
|W |() = y y . 1 2
y 1
y2
En el espacio:
r rr + 2r + cos + sen2 =0 , + r2 sen r 2 r2
En esfera con simetra, n = Pn (cos ) , Rn = {r n } (interior), Rn = {r n1 } (exterior), n = 0, 1, . . .
= R
r 2 R + 2rR R = 0 , + = 0 , s=cos sen + sen sen = 0 (con = 0 si no d [1 s2 ] d + 2 = 0 depende de ) ds ds
1s
. i] Si o , tiene una sola solucin (denida en todo ) con y (o ) = yo , y (o ) = yo
ii] Si y1 , y2 son soluciones de la homognea ( 0 ) con wronskiano |W |(s) = 0 para algn s , la solucin general de la homognea es: y = c1 y1 + c2 y2 . iv] Una solucin particular de [n] es: yp = y2 y + y + by = 0 2 + + b = 0 iii] Si yp es una solucin de [n], la solucin general de [n] es: y = c1 y1 + c2 y2 + yp . y y 2 1 d y1 |W d [fvc]. |W | |
Ecuacin de ondas
La solucin nica de (, t ) = Si F = 0 ,
Coecientes constantes: tt = F ( , t ) , , t R (, 0) = () , t (, 0) = g()
+ct ct
c2
es:
1 2c
1 1 ( + ct ) + ( ct ) + 2 2 c
g(s) ds +
+c[t ]
tt c2 = F, 0, t R tt c2 = F, [0,L], t R tt c2 = F (, t ) , , t R (, 0) = , t (, 0) = g , (, 0) = , t (, 0) = g (,0) = (), t (,0) = g() (0, t ) = 0 (0, t ) = (L, t ) = 0
1 1 () 2 g 2 c 0
viaja hacia la derecha y
1 1 () + 2 g 2 c 0
0 c[t ]
F (s, ) ds d
hacia la izquierda.
Mtodo de coecientes indeterminados para y + y + by = () : Si () = e pm () , pm polinomio de grado m , y no es autovalor hay solucin particular yp = e Pm () . Si es autovalor de multiplicidad r , hay yp = r e Pm () . Si () = e pj ( ) cos q+qk () sen q , pj ,qk de grados j , k , y pq no es autovalor hay yp = ep Pm () cos q + Qm () sen q con Pm y Qm de grado m = mx{j, k } . Si p q es autovalor hay yp = ep Pm () cos q + Qm () sen q . Euler: 1 = 2 reales y = c1 1 + c2 2 doble (real) y = [c1 + c2 ln ] ( 1)+ + b = 0 = p q y = c1 cos(q ln )+ c2 sen(q ln ) p Una yp de 2 y + y + by = h() se obtiene de la [fvc] con () = h()/ 2 , o utilizando que con = es se convierte en y (s)+( 1)y (s)+ by (s) = h es . 2 y + y + by = 0
1 = 2 reales y = c1 e1 + c2 e2 doble (real) y = (c1 + c2 ) e = p q y = (c1 cos q + c2 sen q) ep
, g , F extensiones impares de (), g(), F (, t ) respecto a = 0 a = 0 y = L .
Condiciones = 0 llevan a extensiones pares de los datos iniciales. tt c2 rr + 2 = 0 , r 0, t R =r tt c2 rr = 0 , r 0, t R r r (r, 0) = r (r ), t (r, 0) = rg(r ), (0, t ) = 0 (r, 0) = (r ) , t (r, 0) = g(r )
Si b() 0 , y = lleva [e] a lineal de primer orden en . Con el cambio = y se convierte y + 2y + y = 0 en + = 0 . EDOs importantes resolubles utilizando series:
Transformada de Fourier
La transformada de () absolutamente integrable es (k ) F ( )(k ) Si C1 (R) y absolutamente integrable, () =
1 2
1 2
( ) ek d .
Legendre (1 2 )y 2y + p(p + 1)y = 0 tiene soluciones polinmicas si p = n N :
Si , , C(R) y absolutamente integrables: F ( ) = k F ( ) , F ( ) = k 2 F ( ) . 2 2 F e = 1 ek / 4 . 1 , [, b] ekb ek 2 F 1 (k ) ek = ( ) . Si h() = , F (h) = . k 2 2 1 1 0 en el resto 2 k F e = e / 4 . Si ( g)() =
1 2
(k ) ek dk R .
P0 = 1 , P1 = , P2 = 3 2 1 , P3 = 5 3 3 , P4 = 35 4 15 2 + 3 , ... 2 2 2 2 8 4 8 1 1 2 Se cumple que: 1 Pn Pm d = 0 , si m = n ; 1 P2 d = . 2n+1 n m/ 2 dm m2 m 2 2 Pn () = (1 ) P () satisfacen (1 )y 2y + n(n + 1) 12 y = 0 . d m n Jp ()
p
( s) g(s) ds , es g = g , y es F ( g) = F ( )F (g) . (, t ) =
1 2 t
Bessel: 2 y + y +[2 p2 ]y = 0
m=0
( 1)m [ / 2]2m m! (p+m+1)
La solucin nica de
t = 0 , R, t > 0 es: (, 0) = () , acotada
( s ) e(s)
2 / 4t
ds .
soluciones acotadas en = 0 , con innitos ceros [los de J0 : 2.4.., 5.5.., ... ]. Las soluciones linealmente independientes de ellas son no acotadas en 0 . cos + c2 Si p = 1 , 3 , . . . , las soluciones son funciones elementales p = 1 y = c1 sen . 2 2 2 2p Recurrencia: Jp+1 = Jp Jp1 . Derivadas: p Jp () = p Jp1 () , J ( ) = J . 1 0
EDPs de primer orden
A(, y ) y + B(, y ) = C(, y ) + D(, y )
Problemas de contorno para EDOs
caractersticas = (, y ) = K d B(, y )
dy A(, y )
(Ps )
Hay una nica solucin local (, y ) con valores dados sobre una curva G = g(s), h(s) del plano y si G no es tangente a las caractersticas (s) g A(g, h) h B(g, h) = 0 .
= (, y ) A = C + D ; =y
= (, y ) B = C + D . =
p C1 [, b] , q, r C[, b] , p, r > 0 en [, b] , [py ] qy + ry = 0 , | | + | | , || + | | = 0 . y () y () = y (b)+ y (b) = 0
Si 0 , 0 y q() 0 en [, b] todos los n 0 .
, y n
Los n . Las {yn } son espacio vectorial de dimensin 1 y cada yn tiene exactamente b n 1 ceros en (, b) . yn , ym r yn ym d = 0 , si yn , ym estn asociadas a n = m . Si es C1 a trozos en [, b] la serie de Fourier yn () converge hacia () en los n=1 yn , yn 1 + (, b) en que es continua y hacia 2 ( )+ ( ) en los en que es discontinua.
Clasicacin y formas cannicas de EDPs de 2o orden
Ayy + By + C + Dy + E + F = G
yy = 2 + 2 y y + 2 + yy + yy y y = (, y ) y = y + [ y + y ] + y + y + y = (, y ) = 2 + 2 + 2 + + B2 (, y ) 4A(, y )C(, y )
Problemas conocidos para [Link] n yn , yn , yn y (0) = y (L) = 0
y + y = 0
:
[2n1]2 2 , n = 1, 2, ... 22 L2 [2n1] L sen , 2 2L
Si en una regin del plano es
hiperblica > = 0 , la EDP es parablica . < elptica
n2 2 , n = 1, 2, ... L2 L sen n , 2 L
n2 2 , n = 0, 1, ... L2 L cos n , L 2 L (n=0)
y (0) = y (L) = 0
y (0) = y (L) = 0
[2n1]2 2 , n = 1, 2, ... 22 L2 [2n1] L cos , 2 2L
y (0) = y (L) = 0
(, y ) = K = (, y ) + = G ; (, y ) = K = (, y ) d B = (, y ) P: = (, y ) = K ; + = G = y ( = ) dy 2A d B 4AC B2 = (, y ) E: = (, y ) (, y ) = K ; + + = G = (, y ) dy 2A H:
dy = 2A
B2 4AC
o Serie en senos y cosenos en [L, L]: () = 2 + n cos n + bn sen n , L L n=1
con: n = 1 L
Converge hacia () en los en que su extensin 2L-peridica es continua (y hacia 1 ( )+ (+ ) en los puntos de discontinuidad). Adems converge uniformemente 2 en todo intervalo cerrado que no contenga discontinuidades de la extendida. La serie en senos (cosenos) converge hacia la extensin impar (par) y 2L-peridica de . Problemas singulares conocidos: y + 2y + y = 0 y + y + y = 0 y acotada en = 0, y (1) = 0 y acotada en = 0, y (1) = 0 sen n = y s = J0 (n ) , J0 (n ) = 0
(1 2 )y + y = 0 y acotada en = 1 Pn () de Legendre
() cos n d , n = 0, 1, 2, ... y bn = 1 L L
() sen n d , n = 1, 2, ... L
Si A , B , C constantes:
H:
= =
B B+
B2 4AC 2A B2 4AC 2A
y y
P:
= =y
B 2A
E:
2ABy = 2
=y
4ACB
[Cambios del tipo = epy eq suprimen derivadas de menor orden]. Formas cannicas resolubles:
+ E + F = G ,
+ D = G .
(Pf )
Tienen solucin nica (0, t ) = h0 (t ) ; (L, t ) = hL (t )
t k = F (, t ), (0, L), t > 0 . (, 0) = (), [Link] fsicas (0, t ) (0, t ) = h0 (t ) , , b > 0 (L, t )+ b (L, t ) = hL (t )
[temperaturas o ujo dados o radiacin libre]
(Ph ) tiene slo la solucin y 0 . Si (Ph ) tiene soluciones no triviales {yh } , (Pf ) tiene b innitas soluciones o ninguna segn sea = 0 = 0 la integral: () yh () d .
p()y + g()y = () tiene 1 solucin si y slo si el homogneo y () y () = y (b)+ y (b) = 0
(0, t ) = h0 (t ) ; (L, t ) = hL (t )
Algunas propiedades de funciones trigonomtricas e hiperblicas: 1 sen sen b = 1 cos( b) cos( + b) cos = sen + , 1 + tn2 = cos 2 2 2 1 cos cos b = 2 cos( + b)+ cos( b) sen ( b) = sen cos b cos sen b cos ( b) = cos cos b sen sen b sen cos b = 1 sen( + b)+ sen( b) 2 sen2 =
1 2 cos 2 2
e
Si (Ph ) tiene slo la solucin y 0 e y1 , y2 son soluciones no triviales de [py ]+gy = 0 ( ) = 0 , y ( b )+ y ( b ) = 0 , la solucin de (P ) es: cumpliendo y1 () y1 2 f 2 y1 (s)y2 () , s b p|W |(y1 ,y2 ) , funcin de Green. y () = G(, s) (s)ds , con G(, s) = y ()y (s) 1 2 , sb p|W |(y ,y )
1 2
Separacin de variables para calor y ondas
. X + X = 0 , T + kT = 0 t k = 0 , = XT X + X = 0 , tt c2 = 0 , = XT T + c2 T = 0 rR + nR + rR = 0 , T + T = 0
tt rr + n = 0 , = RT r r
cos2
1 2 + cos 2 2
sen3
3 sen 3 sen 4 4
sh x
cos3
3 cos 3 + cos 4 4
sh = e
, ch = e
+e
sh , th = ch .
[sh ] = ch , [th ] = 1 th2 = 1 . ch2 [ch ] = sh
ch x
th x
Para problemas no homogneos se prueban desarrollos en autofunciones del homogneo. Si las condiciones de contorno no son homognas hay que hacer un cambio de variable. h0 (t )+ h (t ) satisface (0, t ) = h0 (t ) , (L, t ) = hL (t ) . En particular (, t ) = 1 L L L
sh( y ) = sh ch y ch sh y , ch( y ) = ch ch y sh sh y .