4 1

Separación de variables. Repaso de EDOs

ut − ku xx = 0 X "" + λX = 0 utt − c2 u xx = 0 X "" + λX = 0 u xx + uyy = 0 X "" + λX = 0
Algunas EDOs de primer orden y" = f ( x, y) resolubles
u = XT → T " + λkT = 0 u = XT → T "" + λc2 T = 0 u = XY → Y "" − λY = 0
[Si f , f y son continuas en un entorno de ( xo , yo ) existe una única solución con y( xo ) = yo ].
u
urr + urr + urθθ2 = 0 , u = RΘ urr + 2ur r + urθθ2 + cos θ uθ
+ sen2φφθ r2 = 0 , u = RΘΦ → p( x ) y y
sen θ r2 Separables: y" = q(y) . y" = f ( x ) , z = x . y" = f ( ax + by) , z = ax + by .
Θ"" + λΘ = 0 r2 R"" + 2rR" − λR = 0 , Φ"" + µΦ = 0 , ' ' ' '
→ ! µ " ! µ "
Lineales: y" = a( x )y + f ( x ) → y = Ce a( x)dx + e a( x)dx e− a( x)dx f ( x )dx [ sgh + y p ].
r2 R"" + rR" − λR = 0 (sen θ Θ" )" + λ sen θ − sen 2 " "
θ Θ = 0 → ([1 − t ] Θ ) + λ − 1−t2 Θ = 0 $ %
t=cos θ Exactas: M ( x, y)+ N ( x, y) y" = 0 con M = Ux , N = Uy My ≡ Nx → U ( x, y) = C .
Ecuación de ondas / /
# [e] y"" + a( x )y" + b( x )y = f ( x ) , / y y /
utt − c2 u xx = F ( x, t) , x, t ∈ R EDOs lineales de segundo orden |W |( x ) = // 1" 2" //
La solución única de es: a , b , f continuas en I . y1 y2
u( x, 0) = f ( x ) , ut ( x, 0) = g( x )
$ % & x+ct & t& x+c[t−τ ] i] Si y1 , y2 son soluciones de la homogénea con wronskiano |W |(s) *= 0 para algún s ∈ I
1 1 1
u( x, t) = 2 f ( x + ct) + f ( x − ct) + g(s) ds + F (s, τ ) ds dτ e y p es solución particular de [e], la solución general de [e] es: y = c1 y1 + c2 y2 + y p .
2c x −ct 2c 0 x −c[t−τ ] ' y1 f ' y2 f
'
1 x 1 1
ii] Una solución particular de [e] es: y p = y2 |W |
dx − y1 |W |
dx .
Si G ( x ) ≡ 2c 0 g , 2 f ( x )− G ( x ) viaja hacia la derecha y 2 f ( x )+ G ( x ) hacia la izquierda.
  iii] Si xo ∈ I , existe solución única (definida en todo I ) con y( xo ) = yo , y" ( xo ) = y"o .
 utt − c2 u xx = F , x ≥ 0, t ∈ R  utt − c2 u xx = F , x ∈ [0,L], t ∈ R #
utt − c2 u xx = F ∗( x, t), x, t ∈ R µ1 *= µ2 reales → y = c1 eµ1 x + c2 eµ2 x
u( x, 0) = f , ut ( x, 0) = g , u( x, 0) = f , ut ( x, 0) = g → y"" + ay" + by = 0
  u( x, 0) = f ∗ ( x ), ut ( x, 0) = g∗ ( x ) Coeficientes constantes: µ doble (real) → y = (c1 + c2 x ) eµx
u(0, t) = 0 u(0, t) = u( L, t) = 0
µ2 + aµ + b = 0 µ = p ± iq → y = (c1 cos qx + c2 sen qx ) e px
f ∗ , g∗ , F ∗ extensiones impares de f ( x ), g( x ), F ( x, t) respecto a x = 0 o respecto a x = 0 y x = L .
$ % Método de coeficientes indeterminados para y"" + ay" + by = f ( x ) :
v( x, t) = 1 − Lx h0 (t)+ Lx h L (t) satisface v(0, t) = h0 (t) , v( L, t) = h L (t) .
+ $ % + Si f ( x ) = eµx pm ( x ) , pm polinomio de grado m , y µ no es autovalor hay solución
utt − c2 urr + 2r ur = 0 , r ≥ 0, t ∈ R v=ur vtt − c2 vrr = 0 , r ≥ 0, t ∈ R particular y p = eµx Pm ( x ) . Si µ es autovalor de multiplicidad r , hay y p = xr eµx Pm ( x ) .
u(r, 0) = f (r ) , ut (r, 0) = g(r )
−→ v(r, 0) = r f (r), v (r, 0) = rg(r), v(0, t) = 0 $ %
t Si f ( x ) = eµx p j ( x ) cos qx+qk ( x ) sen qx , p j , qk de grados j , k , y p±iq no es autovalor
 $ % $ %
u − c 2 u +u +u = 0 hay y p = e px Pm ( x ) cos qx + Qm ( x ) sen qx con Pm y Qm de grado m = máx{ j, k } . Si

 tt xx yy zz u( x, y, z, t) = px
$ %
 $ t ' 2π' π % ' 2π' π % p ± iq es autovalor hay y p = xe Pm ( x ) cos qx + Qm ( x ) sen qx .
u( x, y, z, 0) = f ( x, y, z) t
 u ( x, y, z, 0 ) = g ( x, y, z )
→ ∂t∂ 4π 0 0 f (·, ·, ·) sen θ dθ dφ + 4π 0 0 g (·, ·, ·) sen θ dθ dφ ,
 t µ1 *= µ2 reales → y = c1 x µ1 + c2 x µ2

t ∈ R , ( x, y, z) ∈ R3 (·, ·, ·) = ( x + ct sen θ cos φ, y + ct sen θ sen φ, z + ct cos θ ) . x2 y"" + axy" + by = 0
Euler: µ doble (real) → y = [c1 + c2 ln x ] x µ
 $ % µ(µ − 1)+ aµ + b = 0
- & 2π& ct & 2π& ct . µ = p ± qi → y = [c1 cos(q ln x )+ c2 sen(q ln x )] x p
 utt − c2 u xx + uyy = 0 , t ∈ R u( x, y, t) =
1 ∂ r f (·, ·) drdθ
√ +
r g(·, ·) drdθ
√ ,
u( x, y, 0) = f ( x, y) , ( x, y) ∈ R 2 → 2πc ∂t 0 0 2 2
c t −r 2 0 0 2 2
c t −r 2 Haciendo x = es en x2 y"" + axy" + by = h( x ) se obtiene y"" (s)+( a − 1)y" (s)+ by(s) = h(es ) .

ut ( x, y, 0) = g( x, y) (·, ·) = ( x + r cos θ, y + r sen θ ) . Si b( x ) ≡ 0 , y" = v lleva [e] a lineal de primer orden en v .

Transformada de Fourier EDOs importantes resolubles utilizando series:

& ∞
La transformada de f ( x ) absolutamente integrable es: fˆ(k ) ≡ I f (k) ≡ √1 f ( x ) eikx dx . Legendre (1 − x2 )y"" − 2xy" + p( p + 1)y = 0 tiene soluciones polinómicas cuando p = n ∈ N :
2π −∞
& ∞ '1 '1
Si f ∈ C1 (R) y absolutamente integrable, f ( x ) = √1 fˆ( x ) e−ikx dk ∀ x ∈ R . P0 = 1 , P1 = x , P2 = 23 x2 − 12 , P3 = 25 x3 − 32 x , . . . ; −1 Pn Pm dx = 0 , si m *= n ; −1 Pn2 dx = 2n2+1 .
2π −∞
m/2 dm $ m2
%
Si f , f " , f "" ∈ C (R) y absolutamente integrables: I( f " ) = −ikI( f ) , I( f "" ) = −k2 I( f ) . Pnm ( x ) = (1 − x2 ) 2 ""
dx m Pn ( x ) satisfacen (1 − x ) y − 2xy + n ( n + 1)− 1− x2 y = 0 .
"

& ∞
Si ( f ∗ g)( x ) = √1 f ( x − s) g(s) ds , es f ∗ g = g ∗ f , y es I( f ∗ g) = I( f )I( g) . $ x %p ∞
soluciones acotadas en x = 0 , con
(−1)m [ x/2]2m
2π −∞ Bessel: x2 y"" +xy" +[ x2− p2 ]y = 0 . J p ( x ) ≡ ∑
2
m =0
infinitos ceros [los de J0 : 2.4.., 5.5.., ...].
m! Γ( p+m+1)
# I(e − ax2
) = √12a e −k2 /4a
.
! " 1 , x ∈ [ a, b] ikb ika Las soluciones linealmente independientes de ellas son no acotadas en 0 .
I−1 f (k ) eiak = f ( x − a) . Si h( x ) = , I(h) = e√ −e . 2 2 $ %
0 en el resto 2π ik 1
I (e
− − ak ) = √12a e−x /4a . Si p = 12 , 32 , . . . , solución expresable con funciones elementales p = 12 → y = c1 sen √ x + c2 cos
√x .
x x
#
ut − u xx = 0 , x ∈ R, t > 0 1
& ∞
2 2p $ %"
La solución única de es: u( x, t) = √ f (s) e−( x−s) /4t ds . Recurrencia: J p+1 = x J p − J p−1 . Derivadas: x p J p ( x ) = x p J p−1 ( x ) , J0" ( x ) = − J1 .
u( x, 0) = f ( x ) , u acotada 2 πt −∞
 2 3

EDPs de primer orden Problemas de contorno para EDOs

+
A( x, y) uy + B( x, y) u x = C ( x, y) u + D ( x, y)
dy
=
A( x, y)
→
características [ py" ]" − qy + λry = 0 p ∈ C1 [ a, b] , q, r ∈ C [ a, b] , p, r > 0 en [ a, b] ,
dx B( x, y) ξ ( x, y) = K (Ps ) ,
αy( a)− α" y" ( a) = βy(b)+ β" y" (b) = 0 |α|+|α" | , | β|+| β" | *= 0 .
# #
ξ = ξ ( x, y) ξ = ξ ( x, y)
→ A uη = C u + D ; → B uη = C u + D Los λn → ∞ . Las {yn } son espacio vectorial de dimensión 1 y cada yn tiene exactamente
η=y η=x 'b
n − 1 ceros en ( a, b) . +yn , ym , ≡ a r yn ym dx = 0 , si yn , ym están asociadas a λn *= λm .
Hay una única solución local u( x, y) con valores dados sobre una curva G = ( g(s), h(s))
Si α · α" ≥ 0 , β · β" ≥ 0 y q( x ) ≥ 0 en [ a, b] todos los λn ≥ 0 .
del plano xy si G no es tangente a las características ⇔ ∆(s) ≡ g" A( g, h)− h" B( g, h) *= 0 .
∞
+f,y ,
Si f es C1 a trozos en [ a, b] la serie de Fourier ∑ n
y ( x ) converge hacia f ( x ) en los
Clasificación y formas canónicas de EDPs de 2 orden o
n =1
+yn , yn , n
1
$ +
%
x ∈ ( a, b) en que f es continua y hacia 2 f ( x )+ f ( x ) en los x en que es discontinua.
−
Auyy + Bu xy + Cu xx + Duy + Eu x + Fu = G ∞ $ %
a
Serie en senos y cosenos en [− L, L]: f ( x ) = 2o + ∑ an cos nπx nπx
L + bn sen L ,
# uyy = uξξ ξ y2 + 2 uξη ξ y ηy + uηη ηy2 + uξ ξ yy + uη ηyy n =1
ξ = ξ ( x, y) 1
& L
1
& L
→ u xy = uξξ ξ y ξ x + uξη [ξ x ηy + ξ y ηx ] + uηη ηy ηx + uξ ξ xy + uη ηxy con: an = f ( x ) cos nπx
L dx , n = 0, 1, 2, ... y bn = f ( x ) sen nπx
L dx , n = 1, 2, ...
η = η ( x, y) L −L L −L
u xx = uξξ ξ 2x + 2 uξη ξ x ηx + uηη ηx2 + uξ ξ xx + uη ηxx
Problemas conocidos para y"" + λy = 0 :

> 0 , la EDP es hiperbólica c.c. y (0) = y ( L ) = 0 y " (0) = y " ( L ) = 0 y (0) = y " ( L ) = 0 y " (0) = y ( L ) = 0
Si en una región del plano es B2 ( x, y)− 4A( x, y)C ( x, y) =
< parabólica . n2 π 2 n2 π 2 [2n−1]2 π 2 [2n−1]2 π 2
elíptica λn L2
, n = 1, 2, ... L2
, n = 0, 1, ... 2 2 , n = 1, 2, ... 2 2 , n = 1, 2, ...
0 1 L 0 1 0 2 L [2n−1]πx 1 L 0 2 L [2n−1]πx 1 L
√ # yn , +yn , yn , sen nπx
L ,2 cos nπx
L , L ó L
2 sen 2L ,2 cos 2L ,2
dx B± B2 − 4AC ξ ( x, y) = K ξ = ξ ( x, y) ( n =0)
H: = → ; → uξη + · · · = G ∗ Problemas singulares conocidos:
dy 2A η ( x, y) = K η = η ( x, y)
# # "" # #$ %"
dx B ξ = ξ ( x, y) xy + 2y" + λxy = 0 xy"" + y" + λxy = 0 (1 − x2 )y" + λy = 0
P: = → ξ ( x, y) = K ; → uηη + · · · = G ∗ y acotada en x = 0, y(1) = 0 y acotada en x = 0, y(1) = 0 y acotada en x = ±1
dy 2A η = y (ó η = x ) 0 1 √ 0 1 0 1
√ # u = xy → senxnπx s = λ x → J0 (wn x ) , J0 (wn ) = 0 Pn ( x ) de Legendre
dx B ± i 4AC − B2 ξ = α( x, y)
E: = → α( x, y)+ iβ( x, y) = K ; → uξξ + uηη +· · · = G ∗ #$ %"
dy 2A η = β( x, y) p( x )y" + g( x )y = f ( x )
(Pf ) tiene 1 solución si y sólo si el homogéneo
+ √ + αy( a)− α" y" ( a) = βy(b)+ β" y" (b) = 0
B+ B2 −4AC # 2Ax − By
Si A , B , C ξ = x− y ξ = x− B
y ξ=√ (Ph ) tiene sólo la solución y ≡ 0 . Si (Ph ) tiene soluciones no triviales {yh } , (Pf ) tiene
H: √ 2A P: 2A E: 4AC − B2 'b
constantes: η = x− B− B2 −4AC
y η=y η=y infinitas soluciones o ninguna según sea = 0 ó *= 0 la integral: a f ( x ) yh ( x ) dx .
2A

[Cambios del tipo u = e py eqx w suprimen derivadas de menor orden]. Si (Ph ) tiene sólo la solución y ≡ 0 e y1 , y2 son soluciones no triviales de [ py" ]" + gy = 0
cumpliendo, respectivamente, αy( a)−α" y" ( a) = 0 , βy(b)+ β" y" (b) = 0 , la solución de (Pf ) es:

Tienen solución única:  y1 ( s ) y2 ( x ) , a ≤ s ≤ x
# # 'b p|W |(y ,y )
ut − ku xx = F ( x, t), x ∈ (0, L), t > 0 ∆u = F en D acotado [Si F ≡ 0 máximo y y( x ) = a G ( x, s) f (s)ds , con G ( x, s) = y ( x)y1 (s2) , función de Green.
,  1 2 , x≤s≤b
u( x, 0) = f ( x ), [Link] físicas u = f en ∂D mínimo se dan en ∂D ]. p|W |(y ,y ) 1 2
# && 2
∆u = F en D
Tiene unicidad salvo constante si se cumple F dxdy = f ds . Funciones de Green
un = f en ∂D D ∂D
# #
∆u = F en D ∆G = δ(ξ − x, η − y) en D
Algunas propiedades de funciones trigonométricas e hiperbólicas: (PD ) . Si G ( x, y; ξ, η ) satisface , para ( x, y) ∈ D
u = f en ∂D G = 0 en ∂D
sen a sen b = 12 [cos( a − b)− cos( a + b)] y vista como función de (ξ, η ) , se le llama función de Green y la solución de (PD ) es:
&& 2
sen ( a ± b) = sen a cos b ± cos a sen b
cos a cos b = 21 [cos( a + b)+ cos( a − b)] u( x, y) = G ( x, y; ξ, η ) F (ξ, η ) dξdη + Gn ( x, y; ξ, η ) f (ξ, η ) ds .
cos ( a ± b) = cos a cos b ∓ sen a sen b D ∂D
sen a cos b = 21 [sen( a + b)+ sen( a − b)] #
∆u = F (r, θ ) , r < R 1
$ % 1 $ 2 r2 σ2 %
sen2 a = 21 [1 − cos 2a] , cos2 a = 12 [1 + cos 2a] , sen3 a = 34 sen a − 14 sen 3a , cos3 a = 34 cos a + 14 cos 3a . → G (r, θ; σ, φ) = 4π ln σ2 +r2 −2rσ cos(θ − φ) − 4π ln R + R2 −2rσ cos(θ − φ)
u( R, θ ) = f (θ )
x −e− x x +e− x
' R' 2π 1
' 2π R2 −r 2
sh x = e 2 , ch x = e 2
sh x
, th x = ch x ; [sh x ]" = ch x , [ch x ]" = sh x , [th x ]" = 1 − th2 x = 1
. → u(r, θ ) = 0 0 σ G (r, θ; σ, φ) F (σ, φ) dφ dσ + 2π 0 R2 −2Rr cos(θ−φ)+r2 f ( φ ) dφ
ch2 x