Ecuaciones Diferenciales de Derivadas P
Ecuaciones Diferenciales de Derivadas P
Ecuaciones Diferenciales de Derivadas P
d e O r d e n Sup e r i or
Pa r te V I
Ecua ci on e s D i fe r e n ci a l e s
D e D e r i va d a s Pa r ci a l e s
I n g . Ra m ón Ab a sca l
2006
Ecuaciones Diferenciales de Orden Superior. 6.3
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
•
dos o más variables independientes, y
una o más derivadas parciales de la función respecto de las variables independientes.
z ( x, y ) = a ∂ z ∂z
Ejemplos:
∂x ∂y
+ b
∂2 z ( x, y )
∂x
P
P
+ cy = 0
2
P P
z ( x, y ) = ∫ ϕ ( x, y ) d x + φ(y)
Como preparación a su estudio, veremos aquí cómo generar una ecuación diferencial ordinaria(1): T
P T
P
Ejemplo 1: Sea una circunferencia con centro en el origen de coordenadas, y radio r. Su ecuación
en el plano es, como sabemos:
x2 + y 2 = r 2
P P P P P r = cte.
P (6.2)
Si derivamos esta ecuación respecto de x, tenemos:
2x + 2y dy = 2 x + 2 y y' = 0
dx
De aquí surge la ecuación diferencial ordinaria:
x + y y' = 0 (6.3)
Pero también
2 x dx + 2y dy = 0
Esta es la ecuación diferencial de la familia. Al integrar, vemos que efectivamente, la (6.2) es la
solución de esta ecuación:
2 ∫ x dx + 2 ∫ y dy = x2 + y2 + C = 0
P P P P donde C = - r2
P P
La (6.2) es, obviamente, solución de dicha ecuación, cualquiera sea el valor de la constante r; es
decir que la ecuación (6.3) tiene infinitas soluciones que difieren entre sí precisamente en el valor
de la constante r: En este caso, las infinitas circunferencias de centro en el origen constituyen en
sentido estricto una familia muy particular, pues lo único que tienen en común es precisamente el
centro en el origen. El único parámetro común que comparten es por tanto el punto (0,0). El valor
de "r" identifica a cada uno de los miembros de la familia.
∴ y' = − x - a
y-b
(1) Es usual esta forma de introducir el estudio de las ecuaciones diferenciales, porque proporciona un panorama
más completo en el que confluyen el concepto de familia de ecuaciones, y su solución.
Ecuaciones Diferenciales de Orden Superior. 6.5
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
que es también una ecuación diferencial ordinaria, de primer orden, que caracteriza a una nueva
familia de circunferencias: las de centro en P. Los parámetros de la misma son obviamente a y b.
También en este caso "r" es el elemento que individualiza a los diferentes miembros de aquella.
P1 B B
S
P2 B B
z1 B B
z2
y
B B
(x1, y1)
B B B B
(x2, y2)
B B B B
x
Decimos "uno o más" valores de z, porque la superficie puede eventualmente tener una, o más de
una, "hojas". De otra forma, si x e y varían en forma continua, el punto P describirá una hoja de la
superficie S.
Si la función no tiene término independiente, la ecuación homogénea:
ϕ ( x, y, z ) = 0 (6.4)
representa a una superficie cónica que contiene al origen de coordenadas. Esto quiere decir que
las directrices de dicha superficie cónica pasan por el origen. En efecto, podemos verificar lo
dicho con ayuda de un ejemplo cualquiera.
Ecuaciones Diferenciales de Orden Superior. 6.6
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
Sea la ecuación:
F ( x, y, z ) = 2 x + 3 y - z = 0
∴ z = 2x + 3y
Entonces, si x = 0, estamos en presencia de la recta directriz z = 3 y,
situada en el plano z y y que pasa por el origen. Lo mismo, si y = 0, entonces z = 2x
es una recta del plano z x. Finalmente, si z = 0, y por tanto,
y = - 2 x,
3
tenemos una recta ubicada en el plano xy
La figura siguiente muestra las tres rectas directrices, correspondientes a cada plano cartesiano, y
obtenidas haciendo nula cada una de las tres variables x, y, z:
z
x=0
z = 3y
y=0
z = 2x
z=0
y = - 2x/3
B 0 1 3 zF B B
C 0 2 6 zC B B
D 1 0 2
E 1 1 5
F 1 2 8 zH
B B
G 2 0 4 z=6 zE
B B
H 2 1 7 zB B B
I 2 2 10
zG
B B
z = 12
x
Ecuaciones Diferenciales de Orden Superior. 6.8
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
Esta última figura, por razones de claridad, está limitada al primer cuadrante, únicamente. El
plano representado en ella por las proyecciones de las líneas de nivel para z = constante,
obviamente se extiende, lo mismo que las líneas de nivel, a los cuatro cuadrantes.
Las líneas AZG y AZC, en la primera figura, y sus prolongaciones en ambos sentidos,
B B B B
Una segunda forma de definir una superficie es la siguiente: Si las variables x, y, z, pueden ser
expresadas en función de un parámetro t,
x = x(t)
y = y(t)
z = z(t)
al variar t en forma tal que x, y, y z varíen en forma continua, las sucesivas posiciones de un
punto P ( x, y, z ) describen una superficie. Por lo que podemos decir que las tres ecuaciones del
sistema anterior representan también la superficie, pero en forma paramétrica.
Si derivamos z respecto de t, recordando que es función de x e y, obtendremos la ecuación:
dz = ∂z dx + ∂z dy
dt ∂x dt ∂y dt
O, multiplicando por d t:
∂ z d x + ∂ z d y = z' d x + z' d y
∂x ∂y
dz = x B
y B B B
Para simplificar su escritura, se acostumbra designar con "u" y "v" a las derivadas:
z'x = u
B B z'y = v
B B (6.5)
De tal forma, la ecuación anterior se escribirá así:
dz = u dx + v dy
Una familia de superficies está representada por una ecuación similar a la (6.4), pero incluyendo
un parámetro "a", característico de la familia, y que por consiguiente, la define:
ϕ ( x, y, z, a ) = 0
Ecuaciones Diferenciales de Orden Superior. 6.9
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
Despejando z en esta ecuación, volvemos a encontrar la relación (6.1), que liga z con las
variables independientes x e y:
z = f ( x, y )
Esta igualdad puede escribirse también, en forma implícita, como:
z - f ( x, y ) = 0 (6.6)
F ( x, y, z, α, β ) = 0
La función
representa asimismo una familia de superficies. En este caso es función de dos parámetros, α y β,
que la caracterizan. Esta ecuación, que para distinguirla de la anterior, en lugar de familia
llamaremos "haz de superficies", asociada a la (6.1), nos permite escribir el sistema siguiente:
F ( x, y, z, α, β ) = 0
(6.7)
z = z ( x, y )
Derivando F con respecto de x e y, se obtienen, respectivamente, las dos ecuaciones siguientes:
∂F + ∂F ∂z
∂x ∂z ∂x
= 0
∂F + ∂F ∂z
(6.8)
∂y ∂z ∂y
= 0
Justificación:
Descompongamos la función F en otras dos, F1 y F2, en la primera de las cuales no aparece la
B B B B
(6.10)
Derivando respecto de x:
∂ F1 = - ∂ F2 ∂z
∂x ∂z ∂x
B B B B
Ejemplo 1: Sea
F ( x, y, z, α, β ) = sen 4x - x -2 + α y 2 - 2 z + β = 0 P P P P
Entonces: z = 1 ( sen 4x - x -2 + α y 2 + β )
P P P P
∂F
Derivando
= 4 cos 4x + 2 x -3
∂x
P P
∂F
∂z
Lo mismo = - 2
∂z = 1 ( 4 cos 4x + 2 x -3 )
∂x
y P P
2
∂F + ∂F ∂z = 4 cos 4x + 2 x -3 - ( 4 cos 4x + 2 x -3 ) = 0
∂x ∂z ∂x
Por tanto: P P P P
∂F + ∂F ∂z = 2αy - 2αy = 0
∂y ∂z ∂y
Por tanto:
F ( x, y, z, α, β ) = e -x + e -y + α z - β x y = 0
Ejemplo 2:
P P P P
Aquí es:
z = 1 ( β x y - e -x - e -y )
α
P P P P
∂F = - e -x - β y
∂x
Entonces P P
∂F = α
∂z
∂z 1 ( β y + e -x )
∂x α
y = P P
∂F + ∂F ∂z
Por tanto:
= - e -x - β y + β y + e -x = 0
∂x ∂z ∂x
P P P P
∂z = 1 ( β x + e -y )
∂y α
y P P
∂F + ∂F ∂z = - e -y - β x + - e -y - β x = 0
∂y ∂z ∂y
Por tanto: P P P P
F ( x, y, z, α, β ) = 0
(6.7)
z = z ( x, y )
Las diferenciales totales de las funciones "F" y "z" son, respectivamente:
∂F dx + ∂F dy + ∂F dz
∂x ∂y ∂z
dF = (6.11)
Ecuaciones Diferenciales de Orden Superior. 6.12
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
∂z dx + ∂z dy
∂x ∂y
y dz = (6.12)
∂F + ∂F ∂z ∂F ∂F ∂z
∂x ∂z ∂x ∂y ∂z ∂y
dF = dx + + dy = 0
Esta última ecuación es igual a cero, pues los términos entre paréntesis son ambos nulos, según
acabamos de probar.
De todo lo anterior resulta, entonces, un sistema formado por tres ecuaciones, las dos que hemos
reunido en el sistema (6.8 ) más la igualdad:
F ( x, y, z, α, β ) = 0
Notemos que las ecuaciones (6.8) contienen las mismas variables que F más, respectivamente, las
derivadas parciales de z respecto de x e y. Por lo tanto, se pueden expresar como funciones
implícitas de las variables mencionadas, haciendo:
∂F + ∂F ∂z ∴ Φ1 x, y, z, α, β, ∂ z
∂x ∂z ∂x ∂x
= 0 B B = 0
∂F + ∂F ∂z ∴ Φ2 x, y, z, α, β, ∂ z
∂y ∂z ∂y ∂y
= 0 B B = 0
el sistema de tres ecuaciones mencionado más arriba puede ser formulado así:
Φ1 ( x, y, z, α, β, u ) = 0
Φ2 ( x, y, z, α, β, v ) = 0
B B
F ( x, y, z, α, β ) = 0
B B
F = α x + β y + α β − z = 0,
Si hacemos:
∂z = - y ∂z = - x
x e y:
∂x ∂y
y
w = υ u2 + η v2 + τ t 2 (6.15)
donde υ, η y τ son los parámetros variables que identifican a cada función de la familia.
P P P P P P
Ω = υ u2 + η v2 + τ t 2 - w = 0
implícita:
P P P P P P
Ecuaciones Diferenciales de Orden Superior. 6.14
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
∂w = w'v = 2 η v
∂v
B B
∂w = w't = 2 τ t
∂t
B B
∴
∂u ∂v ∂t
u
dϕ = ∂ϕ df + ∂ϕ dg = 0
∂f ∂g
(6.16)
De modo enteramente similar a como hicimos en el parágrafo anterior a partir de las ecuaciones
(6.11) y (6.12), obtenemos, en este caso:
∂f + ∂f ∂z ∂f + ∂f ∂z
∂x ∂z ∂x ∂y ∂z ∂y
df = dx + dy = 0
∂g ∂g ∂z ∂g + ∂g ∂z
y
∂x ∂z ∂x ∂y ∂z ∂y
dg = + dx + dy = 0
Ecuaciones Diferenciales de Orden Superior. 6.15
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
Reemplazando en la (6.16):
dϕ= ∂ϕ ∂ f + ∂f ∂z ∂f + ∂f ∂z
∂f ∂x ∂z ∂x ∂y ∂z ∂y
dx + dy +
+ ∂ϕ ∂g + ∂g ∂z ∂g + ∂g ∂z
∂g ∂x ∂z ∂x ∂y ∂z ∂y
dx + dy = 0
dϕ= ∂ϕ ∂f + ∂f ∂z + ∂ϕ ∂g + ∂g ∂z
∂f ∂x ∂z ∂x ∂g ∂x ∂z ∂x
dx +
∂ϕ ∂f + ∂f ∂z + ∂ϕ ∂g + ∂g ∂z
∂f ∂y ∂z ∂y ∂g ∂y ∂z ∂y
+ dy = 0
Esta igualdad debe cumplirse para cualquier valor de las variables x e y, pues no se ha hecho
ninguna salvedad a su respecto. Por lo tanto es condición necesaria que cada una de las sumas
encerradas entre llaves sea igual a cero.
∂ϕ ∂f + ∂f ∂z + ∂ϕ ∂g + ∂g ∂z
∂f ∂x ∂z ∂x ∂g ∂x ∂z ∂x
Es decir: = 0
∂ϕ ∂f + ∂f ∂z + ∂ϕ ∂g + ∂g ∂z
y también
∂f ∂y ∂z ∂y ∂g ∂y ∂z ∂y
= 0
De aquí, despejando los primeros miembros en ambas igualdades, obtenemos las dos ecuaciones
∂ϕ ∂f + ∂f ∂z = − ∂ϕ ∂g + ∂g ∂z
siguientes:
∂f ∂x ∂z ∂x ∂g ∂x ∂z ∂x
∂ϕ ∂f + ∂f ∂z = − ∂ϕ ∂g + ∂g ∂z
∂f ∂y ∂z ∂y ∂g ∂y ∂z ∂y
y
∂y ∂z ∂y ∂y ∂z ∂y
Recordemos que en una razón, el producto de los medios es igual al producto de los extremos.
Aplicando aquí esa propiedad y restando ambos productos entre sí, obtenemos la ecuación:
Ecuaciones Diferenciales de Orden Superior. 6.16
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
∂ f ∂g − ∂g ∂f + ∂z ∂f ∂g − ∂f ∂g +
∂x ∂y ∂x ∂y ∂x ∂z ∂y ∂y ∂z
+ ∂z ∂ f ∂g − ∂g ∂f + ∂z ∂z ∂f ∂g − ∂f ∂g =0
∂y ∂x ∂z ∂x ∂z ∂x ∂y ∂z ∂z ∂z ∂z
Simplificando el último paréntesis, que es nulo, como se ve, y despejando el primer sumando,
obtenemos:
∂ f ∂g − ∂g ∂f =
∂x ∂y ∂x ∂y
= ∂z ∂f ∂g − ∂f ∂g + ∂z ∂ f ∂g − ∂f ∂g
∂x ∂y ∂z ∂z ∂y ∂y ∂z ∂x ∂x ∂z
Si tenemos en cuenta que tanto el primer miembro como los términos entre paréntesis del
segundo son respectivamente los Jacobianos(1)
∂ ( f, g ) U = ∂ ( f, g ) ∂ ( f, g )
P
T P
T
∂ ( x, y ) ∂ ( y, z ) ∂ ( z, x )
R = V =
R = U ∂z + V ∂z
∂x ∂y
que es la ecuación diferencial buscada.
En lo que sigue analizaremos algunos ejemplos de familias de soluciones en las que el vínculo
consiste en que contienen determinadas funciones comunes a la familia.
(2) Sean las funciones u = u (x,y) y v = v (x,y). Se denomina Jacobiano, en homenaje al matemático alemán Carl
∂ (u,v) ∂u ∂v ∂u ∂v
Jacobi, a la expresión:
u'x u'y B B B B B B B B B B
Ejemplo 2: Encontrar una ecuación de derivadas parciales cuya solución sea una dada función de
z - x y de z - y:
Φ ( z - x, z - y ) = 0
Llamemos:
f = z-x y g = z-y
De aquí podemos derivar las siguientes ecuaciones diferenciales:
∂f = -1 ∂f = 1 ∂g = -1 ∂g = 1
∂x ∂z ∂y ∂z
Aplicando la ecuación (6.17), resulta la siguiente igualdad:
- 1 + z'x B B
= 0 + z'x B B
0 + z'y B B - 1 + z'y B B
∴ z'x + z'y = 1
B B B B
Esta es la ecuación buscada. Veremos a continuación algunos ejemplos dentro de esta familia de
ecuaciones, que certifican lo dicho:
Ejemplo 2.1:
Sea Φ ( z - x, z - y ) = ( z - x )2 - ( z - y )2 = 0 P P P P
Operando, obtenemos:
z2 - 2 x z + x 2 - z 2 + 2 z y - y 2 = 0
P P P P P P P P
x2 - y 2 - 2 z ( x - y ) = 0
P P P P
∴ z = 1 x2 - y 2P
P
P
P
= 1 ( x + y )
2 x - y 2
∂z + ∂z = 1 + 1 = 1
∂x ∂y
De aquí:
2 2
∂z + ∂z
Por lo tanto:
= 1 + 0 = 1
∂x ∂y
Ejemplo 2.3:
sen ( z - x ) + sen ( z - y ) = 0
∴ sen ( z - x ) = - sen ( z - y ) = sen ( y - z )
De aquí: z - x = y - z
∴ 2z = x + y
∴ z = 1 ( x + y )
2
Es decir que este caso se reduce al ejemplo 2.1 anterior.
Φ ( z - x2 , z + y2 ) = 0
P P P P
Hagamos:
f = z - x2 P P y g = z + y2 P P
= 0 + z'x B B
0 + z'y B B 2 y + z'y B B
∴ − 4 x y + 2 y z'x − 2 x z'y = 0
B B B B B B B B B B B B
B B B B
Φ = f ( x, y ) + g ( x, y ) = z - x2 + ( z + y2 ) = 0
Ejemplo 3.1:
P P P P
De aquí, resulta:
Ecuaciones Diferenciales de Orden Superior. 6.19
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
2 z = x2 - y 2 P P P P
∴ z = 1 ( x2 - y 2 ) P P P P
2
Derivando respecto de x e y,
∴ ∂z = x ∂z = - y
∂x ∂y
Y por tanto, reemplazando en la ecuación de la familia, que acabamos de obtener (6.18):
y z'x − x z'y − 2 x y = 0
B B B B
2
Entonces, e z-x = 1
P
P
Para que las dos funciones extremas sean iguales, deben serlo sus argumentos respectivos, es
decir:
z - x2 = - z - y 2
∴
P P P P
2 z = x2 - y 2 P P P P
Ecuaciones Diferenciales de Orden Superior. 6.20
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
Ejemplo N° 4:
Nos proponemos con este ejemplo encontrar la ecuación diferencial de la familia de todos los
planos que pasan por el origen de coordenadas. La ecuación correspondiente a esta familia es:
F(x,y,z) = ax + by + cz = 0
Despejando la variable dependiente, z, obtenemos:
z = − a x − b y (6.19)
c c
Y, derivando F respecto de x e y, resulta el sistema:
∂F + ∂F ∂z = a + c -a
∂x ∂z ∂x
= 0
c
∂F + ∂F ∂z = b + c -b
∂y ∂z ∂y
= 0
c
Por otra parte, como
∂z = - a
∂x c
∂z = - b
∂y
y
c
Ecuaciones Diferenciales de Orden Superior. 6.21
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
la ecuación (6.19) de la familia de planos puede por lo tanto expresarse como una ecuación
diferencial:
∂z + y ∂z
∂x ∂y
z = x
∂z + y ∂z
∂x ∂y
o bien: x - z = 0
Ejemplo N° 5:
Trataremos en este caso de hallar la ecuación diferencial de la familia de esferas de radio R cuyo
centro está ubicado sobre el eje de coordenas x. La ecuación algebraica correspondiente es, como
sabemos:
( x - a ) 2 + y 2 + z 2 = R2 P P P P P P P P
∂F + ∂F ∂z = 2 y + 2 z . z 'y = 0
∂y ∂z ∂y
B B (6.21)
-2ax a2
P P
Simplificando, hallamos:
F = y 2 + z 2 + z 2 . z'x2 - R2 = 0
P P P P B B B B P P B P
B P P P
que es la ecuación diferencial de la familia, y que puede reformularse más sencillamente sacando
factor común z2:
z 2 ( z'x2 + z'y2 + 1 ) − R2 = 0
P P
P P B B
P P B P
B P P P
z2 ∂z 2
∂z 2
− R2 = 0
∂x ∂y
P P
P
P
+ P P
P
P
+ 1 P P
∴ ∂z a - x
∂x
= P P
z
También, a partir de la (6.21), podemos ver que:
∂z - y
∂y
=
z
Al reemplazar estos valores en la ecuación de la familia, obtenemos:
z2 P P
(a - x)2 + ( - y)2 P
P
P P
P
P
P
P
+ 1 − R2 = 0 P P
z2 z2 P P P P
( x - a ) 2 + y 2 + z 2 - R2 = 0 P P P P P P P P
∂z = A x2 y
∂y
P P
Integrando z = A ∫ x2 y dy P P
Ecuaciones Diferenciales de Orden Superior. 6.23
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
∂ ( x 2 y2 ) = 2 x 2 y
∂y
Y como P P P P P P
podemos poner:
z = A ∫ ∂ ( x 2 y2 ) d y = A x2 y2 + A φ ( x )
∂y
P P P P P P P P
2 2
∂z = ∂ A x2 y2 + A φ ( x ) = 2 A x2 y + 0 = A x 2 y
∂y ∂y
P P P P P P P P
2 2
Ejemplo 2: En este ejemplo procederemos a analizar una ecuación diferencial de segundo orden,
como la siguiente:
∂2 z
∂x
P
P
+ Cy = 0
2
P P
Integrando dos veces, aparecerá cada vez una función que no depende de x, es decir:
z = - ∫∫ C . y dx dx + ∫ φ1 ( y ) dx + φ2 ( y ) =
∫ ∫ ∫
B B B B
= - Cy dx dx + φ1 ( y ) B B dx + φ2 ( y ) B B
Por lo tanto, calculando las integrales que contiene esta última expresión, obtendremos la
∫
solución buscada:
z = - C.y x dx + x φ1 ( y ) + φ2 ( y ) B B B B
y, por fin:
z = - C x2 y + x φ1 ( y ) + φ2 ( y )
P P B B B B
2
Comprobaremos seguidamente que ésta es la solución. Efectivamente, derivando la función z
hallada, dos veces con respecto a x, obtenemos, sucesivamente:
∂ − C x2 y + x φ1 ( y ) + φ2 ( y ) = - C x y + φ1 ( y )
∂x
P P B B B B B B
2
∂2 = ∂
P
∴ − C x2 y + x φ1 ( y ) + φ2 ( y ) - C x y + φ1 ( y )
∂x ∂x
P
P
P P B B B B B B = -Cy
2
2
con lo que queda demostrado. En todos los casos, las funciones φ se determinan a partir de las
P P P
condiciones de contorno.
Ecuaciones Diferenciales de Orden Superior. 6.24
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
R
y-C
C
x-a
x
a ai
B B
Cada miembro de la familia estará caracterizado por el valor de un parámetro a, en tanto que R y
C son constantes que identifican a la familia. La ecuación genérica de la misma es:
( x - a ) 2 + ( y - C ) 2 = R2
P P P P P P
y' = - ( x - a )
(y- C)
diferente de a, que la individualiza. Sea por ejemplo la circunferencia Γi, para la cual el
Como hemos dicho más ariba, a cada circunferencia de la familia le corresponde un valor
B B
y' ( y - C ) = - ( x - ai ) B B
Veamos qué ocurre cuando x = ai. En tal caso, como x - ai = 0, aparecen dos posibilidades para
B B B B
la variable y:
Ecuaciones Diferenciales de Orden Superior. 6.25
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
y = C
o bien y' = 0
En el primer caso, la ecuación corresponde a la recta horizontal C, lugar geométrico de los
centros de las circunferencias de la familia. Por su parte, cuando y' = 0, la ecuación (6.19) queda
simplificada así:
F = ( y - C ) 2 - R2 = 0
P P P P
Despejando, se obtiene:
( y - C ) 2 = R2
P P P ∴
P y - C = ± R
De donde surgen dos valores posibles para la variable y:
y1 =
B B R + C ó y2 = - R + C
B B (6.24)
Derivando y respecto de x, comprobamos nuevamente que:
y' = 0
Esta es la ecuación de una recta horizontal, y ocurre justamente cuando se cumplen las dos
condiciones siguientes:
Que la tangente a la circunferencia es horizontal, y a su vez,
Que x es igual a ai. B B
circunferencia dentro de la familia, coincide con el valor de la abscisa de su centro. Esto se puede
entender mejor observando la figura de más abajo, que representa algunas circunferencias de la
familia, y en la que hemos supuesto, para fijar ideas, que
R = 2 y C = 1
Ecuación de la Familia
( x - a ) 2 + ( y - 1) 2 = 4
P P P P
y
y1 = C + R
B B
y2 = C - R
B B
Resumiendo, las ecuaciones (6.24) de las envolventes son solución de la ecuación diferencial
paramétrica (6.23), de la familia.
Ecuaciones Diferenciales de Orden Superior. 6.26
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
Primera forma:
Como ya hicimos en la Sección 6.2, ecuación (6.5), llamaremos por simplicidad
z'x = u
B B
y z'y = v
B B
es una solución de la ecuación diferencial (6.25). En tal caso, se trata de una solución general
pues contiene un único parámetro, a.
Por el contrario, una solución de la forma:
F ( x, y, z, a, b ) = 0 B B (6.27)
en la que b es una función del parámetro a,
b = ϕ ( a)
F [ x, y, z, a, ϕ ( a ) ] = 0,
B B
o sea: B B (6.28)
es una solución completa, puesto que contiene dos funciones, a y b.
Ecuaciones Diferenciales de Orden Superior. 6.27
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
∂ F [ x, y, z, a, ϕ ( a ) ] = ∂ F + ∂ F db ∂ F + ∂ F ϕ' ( a ) = 0
∂a ∂a ∂b ∂a ∂b
B B = B B
da
∂F = - ∂F
Por lo tanto:
∴ ϕ' ( a )
∂a ∂b
B B (6.29)
El valor del parámetro "a" se obtiene por eliminación entre las ecuaciones (6.28) y (6.29).
Segunda forma:
Partamos de la ecuación (6.26):
F ( x, y, z, u, v ) = 0
B B (6.30)
Supóngase que por algún procedimiento particular hemos hallado una solución que contiene una
función arbitraria, a la que llamaremos "a". En esta situación, es posible escribir expresiones
como las siguientes:
Φ ( x, y, z, u, v ) = a
B B
u = u ( x, y, z, a ) B B
y v = v ( x, y, z, a )
En efecto, para obtener la primera de ellas, simplemente hemos expresado "a" en función de las
demás variables, de forma idéntica a lo que hicimos al despejar v, a partir de la (6.27) en el
ejemplo anterior. En cuanto a las otras dos, las obtenemos a partir de la ecuación
Φ ( x, y, z, u, v ) - a = 0
B B
Pero
v = v(y)
Por tanto
u = u [ x, y, z, v ( y ), a ] = u ( x, y, z, a ) B B B B
Es obvio que si elegimos expresar estas ecuaciones en esta forma, es únicamente por razón de
conveniencia.
Prosiguiendo con nuestro propósito, como z es función de x e y, su diferencial total responde a la
fórmula:
Ecuaciones Diferenciales de Orden Superior. 6.28
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
dz = ∂ z dx + ∂ z
∂x ∂y
dy = u . dx + v . dy (6.31)
Si integramos ahora, aparecerá una segunda función constante. Por lo tanto, la integral que resulta
es una solución completa de la ecuación diferencial.
También, a partir de
∂z + ∂z dy
∂x ∂y
= 0
dx
∂z
Despejando, obtenemos:
dy ∂
= - x ∂u
∂z ∂v
= -
∂y
dx
Si reemplazamos en la (6.33)
∂F ∂F dx
∂u ∂v
- = 0
dy
Ecuación que también podemos escribir así:
F'u - F'v
B B B B
dx = 0
dy
de donde:
Ecuaciones Diferenciales de Orden Superior. 6.29
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
dx = dy
F'u
B B F'v
B
Por simplicidad, utilizaremos un mismo nombre común, dt, para designar a estas dos razones,
idénticas entre sí. Es decir, llamaremos
dx = dy = dt (6.34)
F'u
B B F'v
B
dt B
(6.35)
dy = F'v B B
dt B
Detengámonos un momento para analizar lo realizado hasta ahora. De hecho, hemos expresado,
siempre de modo formal, las derivadas parciales
∂F ∂F
∂u ∂v
y
como derivadas totales respecto de un mismo parámetro auxiliar, que hemos denominado t.
Volvamos atrás a la ecuación (6.31), y dividamos sus dos miembros por dt. A continuación, en el
resultado obtenido, reemplazaremos las (6.35), alcanzando una nueva igualdad:
dz = u dx + v dy = u F'u + v F'v
B B B B
dt dt dt
Asimismo, partiendo de la ecuación (6.26), F ( x, y, z, z'x, z'y ) = F ( x, y, z, u, v ) = 0
B B B B
y calculando sus derivadas con respecto a las variables x e y, hallamos las igualdades:
∂F + ∂F ∂z + ∂F ∂u + ∂F ∂v
∂x ∂z ∂x ∂u ∂x ∂v ∂x
= 0
∂F + ∂F ∂z + ∂F ∂u + ∂F ∂v
∂y ∂z ∂y ∂u ∂y ∂v ∂y
y = 0 (6.36)
du = ∂ u dx + ∂ u dy
∂x ∂y
dv = ∂ v dx + ∂ v dy
∂x ∂y
y
Ecuaciones Diferenciales de Orden Superior. 6.30
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
dt
dv = ∂v F'u + ∂ v ∂v ∂F + ∂v ∂F
y (6.37)
∂x ∂y ∂x ∂u ∂y ∂v
B B F'v =
B B
dt
u = ∂z
∂x
Como
∴ ∂u = ∂2z
∂y ∂x ∂y
P
P
v = ∂z
∂y
Igualmente
∴ ∂v = ∂2z
∂x ∂x ∂y
P
P
∂F + ∂F ∂z
Vemos que
+ du
∂x ∂z ∂x
= 0
dt
Ecuaciones Diferenciales de Orden Superior. 6.31
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
du = - ( F'x + u F'z ) B B B B B ∴
B dt = - du
dt F'x + u F'z
B B B
dv = - ( F'y + v F'z ) B B B B B ∴
B dt = - dv
dt F'y + v F'z
B B B
También:
dt = dx = dy
F'u
B B F'v
B
dx = dy = dz = - du = - dv
F'u
B B F'v B u F'u + v F'v
B B B B F'x + u F'z
B B B B F'y + v F'z
B B B B
que posibilita llegar a determinar la solución completa de una ecuación diferencial en derivadas
parciales. En los ejemplos que siguen veremos cómo utilizarla para dicho fin:
u F'u + v F'v = 6 u2 + v y
B B B B P P
dx = dy = dz = du = dv
2
6u y B 6u + vy
B P P u 0
Esta ecuación sólo es posible si dv es igual a cero, lo que implica que v es una constante:
dv = 0 ∴ v = cte. = C1 B B
Por tanto
z = 3 u 2 + C1 y P P B B
o bien:
z - C1 y = 3 u2
∴
B B P P
d ( z - C1 y ) = 6 u du B B
dx = 6 u du = 6 du
u
Al integrar en ambos miembros hallamos:
∫ d x = 6 ∫ du
∴ x = 6 u + C2 B B
Por tanto
3 u 2 = z - C1 y =
P
( x - C2 ) 2 B B P
P B B
12
De donde
( x - C2 ) 2 + C y
z =
B B P P
1 B B P
12
Esta es la solución completa de la ecuación diferencial, porque contiene dos funciones arbitrarias,
C1 y C2.
B B B B
Para encontrar la solución general, basta expresar una de ellas como función explícita de la otra.
Pongamos, por ejemplo:
C1 = C y C2 = ϕ ( C )
z = ( x + ϕ ( C )) + C y = x2 + 2 x ϕ ( C ) + ϕ 2 ( C ) + C y
B B B B
2
P P P P P P P P
12 12
Verifiquemos este resultado en la ecuación diferencial:
z = 3 u2 + v y = 3 u 2 + C yP P P P
u = ∂z x + ϕ(C)
∂x
=
6
Ecuaciones Diferenciales de Orden Superior. 6.33
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
v = ∂z = C
∂y
y
P P
12
Ejemplo 2:
Sea la ecuación:
z = ( 2 - u ) x + v y = 2 x - x z'x + y z'y B B B B
F'u = - x
B B y F'v = y
B B
e integramos:
∫ dx = ∫ dy
-x y B B
ln y + ln x = ln C1 B B
o bien y x = C1 B B
A partir de la ecuación original, teniendo en cuenta los valores de F'u y F'v, podemos hacer: B B B B
u . F'u + v . F'v = - u x + v y = v y - u x = z - 2 x
B B B B (6.38)
Yendo de nuevo a la ecuación subsidiaria,
dz = dz = − dx = − dx
u . F'u + v . F'v
B B B B vy - ux F'u B x B B
vy - ux = z - 2x = − x dz
dx
Operando en los dos últimos términos de esta igualdad, hallamos:
2 x dx = z dx + x dz = d ( x z )
Integrando nuevamente podemos hacer:
∫ d(zx) = 2 ∫ x dx + C2 B B
De donde:
z x = x 2 + C2 P P B B ó C2 = z x - x2
B B P P
B B B B
z = x + C2 = x + 1 B B
ϕ ( C1 ) B B
x x
Donde ϕ ( C1 ) es una función arbitraria de x e y. Por ejemplo, si ϕ ( C1 ) = exy, entonces:
B B B B P P
z = x + exy P
P
B B
xy
u = z'x = 1 + x y e - e
xy
P
P
P
P
v = z'y = B B
x exy = exy
P
P
P P
x x2 P P x
Reemplazando en la ecuación implícita, verificamos que es solución de la misma. En efecto:
P P
P
P
x x
dx = du ∴ du = 5 dx
x 5 x B
Lo mismo, dy = - dv ∴ dv = −2 dy
y 2 y B
v = - 2 ( ln y + C2 ) B B
u = dz = 5 ( ln x + C1 ) + 5 - 5 = 5 ( ln x + C1 )
B B B B
dx
y v = dz = - 2 ( ln y + C2 ) - 2 + 2 = - 2 ( ln y + C2 )
B B B B
dy
La figura siguiente nos facilitará definir las condiciones del problema. Pero debe quedar bien en
claro que, aunque por necesidades del dibujo no lo parezca, el desplazamiento de la cuerda en
sentido transversal, es decir, su apartamiento de la horizontal, es de un orden varias veces inferior
a la longitud de la cuerda. En otras palabras, h es un infinitésimo frente a L.
y
F2
F1 dl
B B
B B
y h
x
dx L
0 x
F2 sen (α − ∆α)
que actúan sobre un
F2
∆l, de la cuerda.
elemento diferencial, F
B B
B B
α - ∆α
F2 cos (α − ∆α) = F cos α
∆l
α ∆y
F1 cos α = - F cos α
B B
∆x
B B
α
F1 = - F
B B
F1 sen αB B
de lo que sigue, téngase en cuenta que el ángulo α es muy pequeño, y por lo tanto, el ángulo ∆α
Nuevamente aquí las limitaciones de dibujo distorsionan la realidad. Para una mejor comprensión
Como queda dicho más arriba, las fuerzas que actúan en sentido longitudinal satisfacen la
F1 cos α + F2 cos α = 0
ecuación:
B B B B
cos ∆α ∼ 1 y sen ∆α ∼ ∆ α
Reemplazando arriba, se alcanza la igualdad:
∆ F ∼ F sen α - F2 sen α + F2 ∆α cos α
Si dividimos ambos miembros por ∆l:
B B B B
fuerza aplicada, F, que suponemos constante. Teniendo en cuenta esto, el cociente incremental de
la fuerza por unidad de longitud de la cuerda será:
dF ∼ 1 ( F sen α - F sen α + F ∂ α cos α )
dl ∂l
O, simplificando:
dF ∼ ∂α F cos α
∂l
(6.39)
dl
Por otra parte, sabemos que:
∂ sen α ∂α
= cos α
∂l ∂l
Reemplazando en la (6.9), nos queda:
∼ F ∂ sen α
dF
dl ∂l
En la figura puede verse también que:
Ecuaciones Diferenciales de Orden Superior. 6.38
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
sen α = ∂y
∂l
Reemplazando nuevamente:
∂ ∂y
∼ F = F ∂ y
2
dF
∂l ∂l ∂l2
P
P
dl P P
Como hemos dicho que la separación de la cuerda respecto de la horizontal es muy pequeña,
podemos afirmar que
∂l ∼ ∂x
dF ∼ F ∂2y
∂x2
y por tanto P
P
dl P P
dF ∼ F ∂ y dl
2
∂x2
P
P
(6.40)
P P
Al llegar aquí, recordemos la ley de la dinámica: Fuerza igual a masa por aceleración; en
símbolos
F = m.a
P P
dF = a.dm = µ ∂2y dl
∂t2
P
P
(6.41)
P P
Por tanto, y por carácter transitivo, las ecuaciones (6.40) y (6.41) conducen a a la deducción
F ∂ y = µ ∂ y
siguiente:
2 2
∂x ∂t2
P P
P P
2
P P P P
µ
Si definimos
= c2 P P
F
Ecuaciones Diferenciales de Orden Superior. 6.39
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
obtenemos la
∂2y = c2 ∂2y
∂x2 ∂t2
P P
P P
P P
P P P P
que, como puede verse, es una ecuación entre derivadas parciales, en la que la variable
dependiente, "y", es a su vez función de otras dos: "x" y "t".
Esta ecuación es idéntica, salvo por el valor de la constante c, a la
∂2ϕ = ν2 ∂ ϕ
2
∂t2 ∂x2
P P
P P
P P
P P P P
Empecemos por suponer que la posición inicial de la cuerda describe un triángulo isósceles, es
decir que sujetamos la cuerda exactamente por su punto medio, y la mantenemos así hasta el
momento de dejarla en libertad.
La figura siguiente ilustra la posición inicial descrita más arriba.
Posición de la
Cuerda en el
instante t = 0
L
U L
2
Si llamamos:
X" = - λ2 P P
X" = - λ2 X
de donde
P
λ2
P
y T" = - T
P P
2
P c
P P P
(1) Este método es conocido también como Método del Producto, e incluso Método de Fourier, por razones que se
comprenderán al estudiarlo. Con variantes, el procedimiento es valido también para otras posiciones iniciales
diferentes de la que consideramos aquí.
Ecuaciones Diferenciales de Orden Superior. 6.41
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
X" + λ2 X = 0 P
T" + λ T = 0
P
2
P P
c2 P P
λ
y = sen λ x . cos t
c
λ
y = sen λ x . sen t
c
λ
y = cos λ x . cos t
c
λ
y = cos λ x . sen t
c
A partir de las condiciones de contorno dadas, se deduce que únicamente la primera de estas
y = sen λ x . cos λ t
cuatro igualdades,
(6.42)
c
satisface el sistema. En efecto, solamente esta función cumple con la condición:
y ( 0, t ) = 0, para cualquier t,
así como la de velocidad inicial nula:
∂
∂t
y ( x, 0 ) = 0 P
La primera condición se cumple puesto que para x = 0 es sen x = 0, condición que descarta las
dos últimas soluciones de la lista. En cuanto a la segunda condición, derivando "y" respecto de
"t" obtenemos:
∂ λ t = − λ sen λ x . sen λ t
sen λ x . cos
∂t c c c
ecuación que es igual a cero para t = 0:
Pero, para que y (x, t) sea la solución buscada, la (6.42) debe satisfacer además la condición
inicial:
y ( L, 0 ) = 0
Para ello es necesario que λ L sea múltiplo de π:
λx = λL = nπ
x=L
En realidad, necesitamos que este valor sea igual a h, lo que a primera vista parecería una
incongruencia. Esto se resuelve multiplicando el valor de "y" por un factor, Cn. Como n es B B
y ( x, t ) = Σ Cn sen n π x . cos
B B
nπt (6.44)
n L cL
En el instante inicial, t = 0, la serie toma la forma simplificada:
y ( x, 0 ) = Σ Cn sen
B B
nπx = f(x)
n L
Esta función y ( x, 0 ) = f ( x ), que define la posición de la cuerda en el instante t = 0, coincide
con un pulso triangular, que podemos desarrollar por tanto como serie de Fourier. Dicha serie,
por tratarse de una función impar como muestra la figura siguiente, estará conformada por un
desarrollo en senos:
Ecuaciones Diferenciales de Orden Superior. 6.43
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
h
x
L
U L
2
Es decir que los coeficientes Cn coinciden por lo tanto con los coeficientes bn del desarrollo de
B B B B
Fourier:
∞
Cn =B B
2 ∫ f ( x ) sen n π x dx (6.45)
L 1 L
Es de notar que el período de la función a desarrollar es igual a 2 L, como se puede ver en la
misma figura.
Para calcular los coeficientes Cn, observemos que f ( x ) está constituida por dos tramos,
B B
2hx
L
2h
h
x
L L
2
Las ecuaciones que definen ambos tramos son, respectivamente (Ver figura):
Para
0 < x < L 2h
es : y = x
2 L
y para:
L < x < L es : y = 2h - 2h x = 2h (L - x)
2 L L
Reemplazamos ahora estos valores en la (6.45), y obtenemos:
Ecuaciones Diferenciales de Orden Superior. 6.44
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
∫ x sen n π x d x + ∫ 2 h ( L - x ) sen n π x d x
L/2 L
Cn = 2 2h
B B
L 0 L L L/2 L L
∫ x sen n π x d x + ∫ ( L - x ) sen n π x d x
L/2 L
4h
Cn =B B (6.46)
2
L P P 0 L L/2 L
Para calcular por partes la primera de estas integrales
∫ nπx dx
L/2
x sen
0 L
d v = sen n π x
haremos:
u = x y
L
cos n π x
De donde:
L
nπ
du = dx y v = -
L
∫ n π x d x = - L x cos n π x ∫ cos n π x d x
L/2 L/2 L/2
L
nπ nπ
x sen +
0 L L 0 0 L
Y, efectuando operaciones,
∫ n π x d x = - L cos nπ + sen n π x
L/2 L/2
2 2
L
2nπ n2 π2
x sen P P P P
0 L 2 P P L P P 0
∫ n π x d x = - L2 cos nπ + nπ
L/2
L2
2nπ n π
x sen sen (6.47)
P P P P
2 2
0 L 2 P P P P 2
Procediendo de modo idéntico con la segunda integral de la (6.46):
∫ ( L - x ) sen n π x d x
L
L/2 L
haremos
Ecuaciones Diferenciales de Orden Superior. 6.45
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
u = x- L y d v = sen n π x
L
cos n π x
L
nπ
Por tanto: du = dx y v = -
L
Integrando, hallamos:
∫ nπx dx = ∫ ( x - L ) sen n π x
L L
( L - x ) sen - dx =
L/2 L L/2 L
( x - L ) cos n π x nπx
L L
L L2
nπ n π
= - sen
P P
2 2
L L/2 P P L L/2
nπ - sen n π - sen n π
P P
L2 L2
2nπ n2 π2
= cos
P P P P
2 P P 2
cos n π + L2 sen n π
P P
L2
2nπ n2 π2
= (6.48)
P P P P
2 2
En esta última ecuación hemos simplificado sen n π = 0. Reemplazando ahora las igualdades
P P P P
B B
L2 P P P P P P 2 P P P P 2
Recordando, (6.44), que
∞
Y como sen
2
es alternadamente igual a + 1 y - 1, llegamos finalmente a la solución de la ecuación de la cuerda
vibrante, cuando la posición inicial es un triángulo isósceles:
πx cos π t - 1 sen 3 π x 3πt
y ( x, t ) = 8 h
π2
sen cos + ...
P P L cL 9 L cL
Ecuaciones Diferenciales de Orden Superior. 6.46
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
P P
2
P P P P
la que, poniendo
α = 1
c
podemos modificar así:
∂2y 1 ∂2y = α2 ∂ y
2
= P
P
P P
P
P
(6.49)
P P c2 P P P P P P
∂ r ∂r ∂s
De aquí:
= α = - α
∂x ∂t ∂t
= 1 y
P P P P P P
P P P P P
α −α
2 2
P P P
∂ 2 y = α2 ∂2y + ∂ y ∂2y
P P P
− 2
2
P P
P
P
P
P
P
P
(6.50)
2
P P P P P P
P P P P P
+ + 2
∂x2 ∂r2 ∂s2 ∂r∂s
P
P
P
P
P
P
P
P
(6.51)
P P P P P P
α2 ∂2y + ∂ y
2
− 2 ∂ y
2
= α2 ∂2y + ∂ y
2
+ 2 ∂ y
2
P P P P
P P P P P P P P
= 0
Nos encontramos por tanto con una nueva ecuación diferencial, cuya solución es de la forma:
y = Φ(r) + Ψ (s)
∂y ∂Φ(r)
En efecto,
∂r ∂r
=
∂2y ∂ ∂Φ(r)
y, por lo tanto:
∂r∂s ∂s ∂r
P
P
= = 0
Al reemplazar en ella los valores asignados a las variables r y s, la solución queda reformulada
y ( x, t ) = Φ ( x + α t ) + Ψ ( x - α t )
así:
(6.52)
Esta es la solución propuesta por D'alembert.
Φ' ( x ) = 1 A' ( x ) + 1 B ( x )
2 α
De manera similar, se obtiene:
Ψ' ( x ) = 1 1
α
A' ( x ) - B(x)
2
Si integramos Φ' (x) desde t = 0 hasta un instante t, recordando que Φ es función de r = x + α t:
x + αt
Φ ( x, t ) = 1 ∫ 1
α
A' ( x ) + B(x) dx
2 x
x + αt x + αt
Φ ( x, t ) = 1 ∫ 1 ∫
2α
A' ( x ) d x + B(x) dx
2 x x
De forma idéntica:
x - αt x - αt
Ψ ( x, t ) = 1 ∫ 1 ∫
2α
A' ( x ) d x - B(x) dx
2 x x
y ( x, t ) = Φ ( x , t ) + Ψ ( x, t )
Y como
x - αt x - αt
1 ∫ 1 ∫
2α
+ A' ( x ) d x - B(x) dx
2 x x
x + αt
∫ ∫
x
∴ y ( x, t ) = 1 A' ( x ) d x - 1 A' ( x ) d x +
2 x 2 x - αt
x + αt
∫ ∫
x
1 1
2α 2α
+ B(x) dx + B(x) dx
x x - αt
Finalmente
x + αt
1 [ A ( x + αt ) + A ( x - αt ) ] + 1 ∫
2α
y ( x, t ) = B(x) dx
2 x - αt
Dada una función de dos variables Φ ( x, y ), se conoce como Ecuación de Laplace a aquella cuya
expresión simbólica es:
Ecuación de
Laplace
∇2 Φ ( x, y ) = 0
P P
P P P P
y ρ = √ x2 + y 2 P P P P
∂Φ = ∂Φ ∂ρ + ∂Φ ∂θ
∂x ∂ρ ∂x ∂θ ∂x
∂2Φ = ∂2Φ ∂ρ 2
+ ∂Φ ∂ ρ + ∂ Φ
2 2
∂θ 2
+ ∂Φ ∂ θ
2
∂ x2 ∂ρ2 ∂x ∂ ρ ∂ x2 ∂θ2 ∂x ∂ θ ∂ x2
P P P P P P P
P P P P P P P
P P P P P P P P P P P P P P
(6.53)
∂2Φ = ∂2Φ ∂ρ + ∂Φ ∂ ρ + ∂ Φ ∂θ + ∂Φ ∂ θ
e igualmente:
2 2 2 2 2
∂ y2 ∂ρ2 ∂y ∂ ρ ∂ y2 ∂θ2 ∂y ∂ θ ∂ y2
P P P P P P P
P P P P P P P
P P P P P P P P P P P P P P
(6.54)
Por otra parte, derivando a partir de
ρ = √ x2 + y 2 = ( x2 + y 2 ) 2
1
P P P P P P P P
obtenemos:
∂ρ = 1 2x x
√ x + y
= (6.55)
∂x ρ
P
2 2
2 P P P P P P
∂ρ = 1 2y y
√ x + y
y = (6.56)
∂y 2 P P
2
P P
2
P ρ
P
=
P P P P P
∂2ρ ρ2 - x2 =
P P P P P
= y2
∂x ρ3 ρ
P
P
P
P
P
P
P
P
(6.57)
2 3
P P P P P P
= sen θ
Además, como
y
x cos θ
si derivamos también esta expresión respecto de x, obtenemos:
Ecuaciones Diferenciales de Orden Superior. 6.51
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
∂ y 1 cos2 θ ∂ θ + sen2 θ ∂ θ 1 ∂θ
∂x cos θ ∂x ∂x cos θ ∂x
= P P P P =
2 2
x P P P P
Comparando las dos últimas igualdades, por carácter transitivo surge que
∂θ y
∂x ρ2
= - (6.58)
P
= = (6.59)
2 2 2 2 2 2
P P x + y P P P P P P (x + y ) P P P P P P P P P P
∂ x 1 − sen2 θ ∂ θ − cos2 θ ∂ θ 1 ∂θ
∂y sen θ ∂y ∂y sen θ ∂y
= P P P P = -
2 2
y P P P P
∂2 θ ∂ x -2xy = − 2xy
∂y ∂y ρ4
P
P
= = (6.61)
2 2 2 2 2 2
P P x + y
P P P P P P (x + y ) P P P P P P P P P P
Ahora, reemplazando las igualdades (6.55) a (6.61), en las ecuaciones (6.53) y (6.54), resulta:
∂2Φ = ∂2Φ x2 + ∂Φ y2 + ∂ Φ
2
y2 + ∂Φ 2xy
∂ x2 ∂ρ2 ρ2 ∂ρ ρ3 ∂θ2 ρ ∂θ ρ4
P P P P P P
P P P P P P
4
P P P P P P P P P P P P P P
∂2Φ = ∂2Φ y2 + ∂Φ x2 + ∂ Φ
2
x2 − ∂Φ 2xy
∂ y2 ∂ρ2 ρ2 ∂ρ ρ3 ∂θ2 ρ4 ∂θ ρ4
y P
P
P
P
P
P
P
P
P
P
P
P
P P P P P P P P P P P P P P
P P P P P P P P P P P P P P
P P
P P P P P P P P P P
∇2Φ = 1 ∂ ρ ∂Φ + 1 ∂2Φ
ρ ∂ρ ∂ρ ρ2 ∂θ2
P P
P P
P P
P P P
∂ ρ ∂Φ ∂2Φ = 0
en coordenadas
+
1 1 polares
ρ ∂ρ ∂ρ ρ ∂θ2
P
P
P
P
(6.62)
2
P P P
Φ ( ρ, θ ) = R ( ρ) . Τ ( θ )
De esta última igualdad se deducen, por derivación, las siguientes:
Ecuaciones Diferenciales de Orden Superior. 6.53
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
∂Φ = T dR
∂ρ dρ
P
P
(6.63)
∂Φ = R dT
∂θ dθ
P
∂2Φ = R d T
2
∂θ2 dθ2
y P
P
P
P
P
P
Reemplazaremos a continuación estos valores en la ecuación de Laplace (6.62), con lo que ésta
P P P P
adopta la forma:
1 ∂ ρT dR + R d2T
ρ ∂ρ dρ ρ dθ
P
P
P
P
= 0
2 2
P P P
= 0
2
P P P P P
1 d2T = − ρ dR − ρ2 d2R
dθ R dρ dρ2
P P P
P P P
2
T P P R P P P
n2 = − 1 d2T
dθ
P P
P
P
(6.64)
2
T P P
n2 = ρ dR + ρ2 d2R
dρ dρ2
P P
P P
P P
R R P P
Al expresar estas dos últimas ecuaciones en forma implícita hallamos el sistema, homogéneo,
siguiente:
d2T + Tn2 = 0
dθ2
P
P
P P
d2R + 1 dR − n2
P P
dρ2 ρ dρ ρ2
P
P
P
P
R = 0 (6.65)
P P P P
T ( θ ) = An cos n θ + Bn sen n θ B B B B
Ecuaciones Diferenciales de Orden Superior. 6.54
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
dT = n ( - An sen n θ + Bn cos n θ )
dθ
En efecto, B B B B
P P B B B B (6.66)
2
P P
R ( ρ ) = C n ρ n + Dn ρ - n B B P P B B P P
Efectivamente:
d R = n . C ρ n -1 - n D ρ - n -1
dρ
n n B B P P B B P P
d 2 R = n ( n - 1 ) C ρ n-2 - n ( - n - 1 ) D ρ -n-2
d ρ2
y P
P
n n B B P P B B P P
P P
Reemplazando en la (6.65):
d2R + 1 dR − n2
dρ2 ρ dρ ρ
P
P
P
P
R = 0 (6.65)
2
P P P P
( n 2 - n ) Cn ρ n - 2 + ( n2 + n ) Dn ρ - n - 2 + n Cn ρ n - 2 - n Dn ρ - n - 2 -
obtenemos:
P P P P P P
- n 2 Cn ρ n - 2 - n 2 Dn ρ - n - 2 =
P B B P P B B P B B P B B P
P P P P
= Cn ρ n - 2 ( n2 - n + n - n2 ) + Dn ρ - ( n + 2 ) ( n2 + n - n - n2 ) = 0
P B B P P B B P
B B P P P P P P B B P P P P P P
Otra forma de verificarlo consiste en reemplazar directamente esta última igualdad en la ecuación
de Laplace
1 ∂ ρ ∂Φ + 1 ∂2Φ = 0
ρ ∂ρ ∂ρ ρ ∂θ2
P
P
P
P
(6.62)
2
P P P
Si en ésta multiplicamos ambos miembros por ρ 2, podemos expresarla bajo una nueva forma: P P
ρ ∂ ρ ∂Φ + ∂ Φ = 0
2
∂ρ ∂ρ ∂θ2
P
P
(6.67)
P
∂Φ dR
∂ρ dρ
= T
∂
resulta:
T ρ + Tρ d R
2
dR = T dR
∂ρ dρ dρ dρ2
P
P
P P B B B B P P
2
P P
∂2 Φ
de donde:
= R d2 T = - n2 RT
∂θ dθ
P P P
P P
P P P
2 2
P P P P
T ρ + Τ ρ2 d R - n2 RT = 0
2
dR
dρ dρ2
P
P
P P P P
n Cn ρ n - n Dn ρ - n + n ( n - 1 ) Cn ρ n - n ( - n - 1 ) Dn ρ - n - n 2 R =
= ( n + n 2 - n ) Cn ρ n + ( - n + n 2 + n ) Dn ρ - n - n 2 ( Cn ρ n + Dn ρ - n ) =
B B P P B B P P B B P P B B P P P P
= n 2 Cn ρ n + n 2 Dn ρ - n - n 2 ( Cn ρ n + Dn ρ - n ) = 0
P P B B P P P P B B P P P P B B P P B B P P
P P B B P P P P B B P P P P B B P P B B P P
y B B B B
Φ0 ( ρ, θ ) = R0 ( ρ) . Τ0 ( θ ) = ( C0 + D0 ) A0
Si n = 0:
B B B B B B B B B B B B
B B B B B B P P B B B B B B P P
Ecuaciones Diferenciales de Orden Superior. 6.56
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
n=1
αn = An Cn βn = Bn Cn γn = An Dn δn = Bn Dn
donde
B B B B B B B B B B B B B B B B B
y B B B B B B B
∂ 2 ϕ ( x, t ) + b ( x ) ∂ ϕ ( x, t ) + c ( x ) ϕ ( x, t ) +
∂x2 ∂x
a(x) P
P
P P
∂ 2 ϕ ( x, t ) + g ( x ) ∂ ϕ ( x, t )
∂t2 ∂t
+ d(x) P
P
= 0 (6.68)
P P
= P
P
P
P
P
P
P P 0 P P P P
Observemos que para el cálculo de esta transformada, hemos tenido en cuenta en primer lugar
que la operación de transformar involucra el dominio del tiempo, t, y por lo tanto es
independiente respecto de la variable x. En el último paso, además, hemos invertido el orden de
integración y derivación.
De modo semejante se obtiene la transformada del segundo término. O sea
∞
b ( x ) ∂ ϕ ( x, t ) ∫ b(x) ∂ ϕ e dt = b(x) ∂ L ϕ ( x, t )
-st
∂x ∂x ∂x
L = P
P
∞
d(x) ∂ ϕ ∫ d(x) ∂ ϕ e dt =
2 2 -st
∂t2 ∂t2
L P
P
= P
P
P
P
P P 0 P P
∴ d(x) ∂ ϕ s 2 L ϕ ( x, t ) - s ϕ ( x, 0 ) - ∂ ϕ ( x, 0 )
2
∂t2 ∂t
L P
P
= d(x) P P
P P
A su vez
∞
g(x) ∂ϕ ∫ g(x) ∂ϕ g dt=g(x) s L ϕ ( x, t ) - ϕ ( x, 0 )
-st
∂t ∂t
L = P
P
Para simplificar las ecuaciones que siguen, utilizaremos en adelante la forma usual de denominar
a la transformada como una función de s, es decir:
L { f ( x, t )} = F ( s )
En este caso:
∞
Φ ( x, s ) = L { ϕ ( x, t ) } = ∫ ϕ ( x, t ) e
-st
P
P
dt
0
Con esta simbología, y reemplazando en la ecuación diferencial (6.68) los valores obtenidos hasta
aquí, la misma queda representada por la igualdad:
a(x) ∂ Φ ( x, s ) + b ( x ) ∂ Φ ( x, s ) + c ( x ) Φ ( x, s ) +
∂x2 ∂x
+ d ( x ) s 2 Φ ( x, s ) - s ϕ ( x, 0 ) - ∂ ϕ ( x, 0 )
P P
∂t
P P +
+ e(x) s Φ ( x, s ) - ϕ ( x, 0 ) = 0
La única derivada respecto de t que aparece aquí es nula, porque la primitiva, ϕ ( x, 0), es función
de x solamente. Todas las demás derivadas lo son respecto de una única variable, x.
Por tanto, suprimida aquella, podemos considerar a las restantes como derivadas totales respecto
de x, con lo cual la ecuación diferencial se convierte, en el dominio transformado, en una
ecuación diferencial ordinaria, que se deberá resolver como tal:
a ( x ) d Φ ( x, s ) + b ( x ) d Φ ( x, s ) + c ( x ) Φ ( x, s ) +
dx2 P dx
s 2 Φ ( x, s ) - s ϕ ( x, 0 ) + e ( x ) s Φ ( x, s ) - ϕ ( x, 0 )
P
+ d(x) P P = 0
Ecuaciones Diferenciales de Orden Superior. 6.58
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
Ejemplo 1:
Como primer ejemplo analizaremos las ecuaciones de una línea de transmisión, que fueron
históricamente conocidas como Ecuaciones Telegráficas, porque definen el comportamiento de la
tensión y corriente eléctricas, en función de los parámetros físicos de un par de conductores de
cobre, como los utilizados en los sistemas alámbricos de comunicaciones.
Las relaciones entre la tensión y la intensidad de la corriente a lo largo de una tal línea de
transmisión de energía eléctrica satisfacen (Ver la figura a continuación), el sistema de
ecuaciones siguiente:
i ( x, t ) r .d x l .d x
Modelo matemático de
una línea de
transmisión de Energía v ( x, t ) g .d x c .d x
Eléctrica
dx
∂ v ( x, t ) = r . i ( x, t ) + l ∂ i ( x, t )
∂x ∂t
∂ i ( x, t ) = g . v ( x, t ) + c ∂ v ( x, t )
∂x ∂t
∂i g.v + c ∂v
∂x ∂t
= (6.69)
Ecuaciones Diferenciales de Orden Superior. 6.59
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
El problema consiste en hallar las funciones v ( x, t ) e i ( x, t ), que dan los valores de la tensión
y corriente en cada punto de la línea, y en cada instante.
Si bien estas ecuaciones pueden resolverse sin problemas en el campo de las variables x y t, aquí
recurriremos, a manera de ejemplo, al auxilio de la transformada de Laplace. La idea consiste en
redefinir las ecuaciones de modo que aparezca una sola variable independiente en cada una de
ellas. Empezamos por derivar la primera respecto de x y la segunda respecto de t:
∂2v = r ∂i + l ∂ i
2
∂x2 ∂x ∂t∂x
P P
P P
∂2i g ∂v + c ∂ v
P
∂x∂t ∂t ∂t2
P
P
= P
P
P P
∂x2 ∂x ∂t ∂t2
P P
P P
P P P P
Tomando para
∂ i
∂x
el valor definido por la (6.69), obtendremos
∂2v = rg.v + rc ∂v ∂v ∂2v
∂x2 ∂t ∂t ∂t2
P
P
+ lg + lc P
P
P P P P
o, agrupando términos:
∂2v = rg.v + (rc + l g ) ∂v + l c ∂2v
∂x2 ∂t ∂t2
P
P
P
P
(6.70)
P P P P
(6.71)
P P P P
L { v ( x, t ) } = V ( x, s ) = V
y L { i ( x, t ) } = I ( x, s ) = I
Y también:
Ecuaciones Diferenciales de Orden Superior. 6.60
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
vo = v ( 0)
B B v'o = d v ( 0 )
B B
dt
Con tales simplificaciones la ecuación resultante en el caso de la diferencia de potencial, V (x,t),
queda así:
d 2 V = r g . V + ( r c + l g ) ( s V - v ) + l c ( s 2 V - s v - v' )
P
P
o o o B B P P B B B B
dx2 P P
o o P P B B B B
dx2 P P
∴ d 2 V - ( r + l s ) ( g + c s ) V + l ( g + c s ) v + l c v' = 0
P
P
o o B B B B
dx2 P P
Esta ecuación, en el dominio de la variable de Laplace, es como vemos una ecuación diferencial
ordinaria. Para resolverla, bastará hallar la primitiva correspondiente.
o o B B B B
dx2 P P
Cuando la línea puede ser considerada sin pérdidas, es decir, cuando la resistencia específica
(Resistencia por unidad de longitud) del conductor es despreciable, y la aislación entre ambos
conductores es muy grande, es decir
r ≅ 0 y g ≅ 0
las ecuaciones anteriores se simplifican así:
d 2 V - l c ( s 2 V - s v - v' ) = 0
P
P
o o P P B B B B
2
d x P P
d2I P
P
- l c ( s 2 I - s i o - i 'o ) = 0
P P B B B B
2
d x P
Ejemplo 2:
Ecuación de las ondas (Ecuación de D'Alembert):
∂2ϕ ∂2ϕ
La ecuación:
= ν2
∂t2 ∂x2
P P
P P
P P
P P P P
Ecuaciones Diferenciales de Orden Superior. 6.61
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
debida a D'Alembert, es conocida como Ecuación de las ondas. Más arriba hemos resuelto en la
forma tradicional la ecuación de la cuerda vibrante, enteramente similar(1) a la que ahora T
P T
P
dt
0
∂t ∂x
P
P
e P
P
P P
P
P
e P
P
dt
2 2
0 0
P
P P P
∂t ∂x
P
P
e P
P
P P
P
P
P
P
dt
2 2
0 0
P P
P P
Obsérvese que este último paso está plenamente justificado porque, al no ser e- s t función de x, a P P
función.
Podemos ahora formular la siguiente igualdad:
(1)
P
T Similitud lógica si consideramos que el movimiento de la cuerda vibrante, al ser dejada en libertad, consiste en
P
T
s2 Φ ( x, s ) - s ϕ ( x, 0 ) - ∂ ϕ ( x, 0 ) = ν 2 ∂ 2 Φ ( x, s )
∂t ∂ x2
P
P
P P P P
P P P P
- P
P
2 2
P P P P
Φ ( x, s ) = A e v
s s
v x
P P + Be P P
Para obtener los valores de A y B, recurrimos a las condiciones de contorno que hemos
especificado para el problema, es decir:
ϕ ( 0, t ) = ν ( t )
cuya transformada es
Φ ( 0, s ) = L { ν ( t ) } = V ( s )
y ϕ ( ∞, t ) = 0
Φ ( ∞, s ) = 0
de donde:
Φ ( ∞, s ) = A e ∞ + B e - ∞ = A e ∞ = 0
P P
y P P P P P P
Φ ( x, s ) = V ( s ) e
s
v x
P P
ϕ ( x, t ) = ν ( t ) [ δ ( t − x ) . H ( t − x )]
ν ν
*
Finalmente, de acuerdo con los valores que toma la función de Heaviside, según que el
argumento sea positivo o negativo, llegamos al resultado siguiente:
= ν(t) δ(t − x ) t ≥ x
ν ν
*
si
ϕ ( x, t ) =
t < x
ν
= 0 si
6.14 - Problemas.
∂F = b + ∂z ∴ ∂z
∂y ∂y ∂y
= 0 = -b
∂z + ∂z = - b (x+y)
expresarla como ecuación diferencial.
∂x ∂y
R:
a
Ecuaciones Diferenciales de Orden Superior. 6.64
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
∂y ∂z ∂y ∂y ∂z ∂y
Aquí:
2 + z'x B B
= 0 + z'x B B
∴ - 4 y - 2 y z'x + 2 z'y = 0
B B B B
0 + z'y - 2 y + z'y
∂z - y ∂z
B B B B
∂y ∂x
R: = 2y
Solución:
f2 - g2 = z2 + 4xz + 4x2 - z2 +2zy2 - y4 = 0
P P P P P P P P P P P P P P
∴
4 2
z = y - 4x = 1
P
P
P
P
y4 - 4x2 =
P
P
P
P
y2 - 2x
P
P
4x + 2y2 2 y2 + 2x 2
∂z - y ∂z
P P P P
∴
∂x ∂y
= y - y(-1) = 2y
como la de la figura, conductora del calor, en uno de cuyos bordes el potencial calórico ϕ
6.14.4 - Ecuación de Laplace, o del Potencial. Sea una placa rectangular, delgada y uniforme,
aplicado responde a una función senoidal, mientras que los otros tres se mantienen a potencial
cero:
y
ϕ=0
l
ϕ=0 ϕ=0
ϕ = sen x π
x
0
∇2 ϕ ( x, y ) = ∂2ϕ + ∂ ϕ
2
∂x ∂y2
P P
P
P
P
P
= 0
2
P P P P
Y" = - X" = λ2 P P
Y X
Y a partir de aquí podemos establecer un sistema de ecuaciones diferenciales, cada una de ellas
de una única variable:
X" + λ2 X = 0
Y" - λ2 Y = 0
P P
P P
Por tanto:
ϕ ( x, y ) = [ A sen ( λ x ) + B cos ( λ x )] . ( C e λ y + D e - λ y )
ϕ ( x, y ) = ( M e λ y + N e - λ y ) sen ( λ x ) + ( P e λ y + Q e - λ y ) cos ( λ x )
P P P P
P P P P P P P P
Las condiciones de contorno indican que el potencial aplicado varía como la función seno
ϕ ( x, 0 ) = sen x
Por esto el término en coseno de x debe ser nulo. También es: λ = 1. La solución se reduce
entonces a:
ϕ ( x, y ) = ( M e y + N e - y ) sen x P P P P
Resta determinar los valores de las constantes M y N, para lo cual recurriremos de nuevo a las
condiciones de contorno:
De ϕ ( x, 0 ) = sen x resulta, al reemplazar y = 0 en la solución :
ϕ ( x, 0 ) = ( M + N ) sen x
Por lo que M + N = 1 ∴ Ν = 1− Μ
Por su parte, de ϕ ( x, l ) = 0, deducimos que
ϕ ( x, l ) = - ( M el + N e - l ) sen x = 0 P P P P
Y de aquí: M e l + ( 1 - M ) e -l = 0
P P P P
Despejamos M
M ( e l - e -l ) = - e -l
P P P P P P
∴ Μ =
l
e-P
P
2 sh l
y N = − el P
P
2 sh l
6.14.5 - Resolver el problema anterior pero con las condiciones de contorno siguientes:
a) ϕ ( x, 0 ) = sen λ x
ϕ ( x, l ) = 0
ϕ ( 0, y ) = 0
ϕ ( π/λ, y ) = 0
Ecuaciones Diferenciales de Orden Superior. 6.67
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
b) ϕ ( x, 0 ) = sen x - sen 3x
ϕ ( x, l ) = 0
ϕ ( 0, y ) = 0
ϕ ( π, y ) = 0
R: ϕ ( x, y ) = e (y - l ) - e - sen x −
(y - l )
P
P
P
P
P
P
e 3 (y - l ) - e - 3 (y - l )
P
P
P
P
P
P
sen 3 x
2 sh l 2 sh 3 l
a) ∂ z + ∂ z = 0
∂x ∂y
Solución:
z ( x, y ) = X ( x ) . Y ( y )
∴ ∂z ∂z
∂x ∂y
= YX' y = XY'
Nota: La simetría del sistema original implica que también es solución z = d e λ (y - x). Por tanto, P P
∇2 z ( x, y ) = ∂2z + ∂2z
∂x ∂y
P P
P
P
P
P
= 0
2 2
P P P P
R: z ( x, y ) = Σ ( M k e kλ y + Nk e - kλ y ) sen kλ x + ( Pk e
B B P P B B P P B B P
kλ y
P + Q k e - kλ y ) cos kλ x
B B P P
∂z ∂z
∂x ∂y
c) - y = 0
R: z ( x, y ) = k . y α . e α x P P P P
Aplica cion e s M a t la b.
Solución de ecuaciones con derivadas parciales respecto de una sola variable independiente.
» % Ejemplo : dz/ dx = 2 *y^ 2 *x
» syms x y
» z = 2 *int('y^ 2 *x',x)
z = y^ 2 *x^ 2
» % Para o btener la so lució n debemo s ag reg ar una funció n arbitraria de y (C o nstante
respecto de x):
» % z = y^ 2 *x^ 2 + 2 *f(y)
Time: 1 .2 5
Time: 1 .5
Time: 1 .7 5
Time: 2
Time: 2 .2 5
Time: 2 .5
Time: 2 .7 5
Time: 3
Time: 3 .2 5
Time: 3 .5
Time: 3 .7 5
Time: 4
Time: 4 .2 5
Time: 4 .5
Time: 4 .7 5
Time: 5
4 2 8 successful steps
6 2 failed attempts
9 8 2 functio n evaluatio ns
1 partial derivatives
1 4 2 LU deco mpo sitio ns
9 8 1 so lutio ns o f linear systems
»
» % Las filas de la matriz "u" so n la so lució n de la ecuació n. También, en este caso , la
derivada en el mo mento inicial co incide co n la fo rma de la o nda al pro mediar el tiempo .
Po r tanto , no s da una imag en de la misma.
» % La representamo s a co ntinuació n.
» ut0 = 3 *sin(pi*x).*exp(sin(pi/ 2 *y));
» pdesurf(p,t,ut0 )
» g rid
Ecuaciones Diferenciales de Orden Superior. 6.72
Parte VI - Ecuaciones Diferenciales de Derivadas Parciales.
Matlab puede ser una excelente herramienta para verificar la solución de ecuaciones en
general. Ejemplo, Verificar el resultado de la ecuación: dz/dx - y*dz/dy = 0. Problema
6.14.6c
»
» % Verificació n:
»
» diff(diff(f,x),x)
ans =
- (exp(1 / 2 *l/ sinh(l))*exp(y)+ (1 -exp(1 / 2 *l/ sinh(l)))*exp(-y))*sin(x)
» diff(diff(f,y),y)
ans =
(exp(1 / 2 *l/ sinh(l))*exp(y)+ (1 -exp(1 / 2 *l/ sinh(l)))*exp(-y))*sin(x)
» diff(diff(f,x),x) + diff(diff(f,y),y)
ans = 0