Ecuaciones Diferenciales Parciales

(EDPs)
Escuela Profesional: Ingeniería Civil

Coaquira Cárdenas Víctor

[email protected]

FACULTAD DE INGENIERÍA DE MINAS, GEOLOGÍA y CIVIL

Dapartamento Académico de Matemática y Física

Ayacucho, April 10, 2025

Introducción a las EDPs
Preliminares

Consideremos x = (x1 , x2 , . . . , xN ) un elemento del espacio eu-

clideo RN , se definen:
Definición (Bola abierta)
El conjunto formado por todos puntos x ∈ RN que se
encuentran a una distancia menor que r > 0 de un punto fijo
xo ∈ RN es llamado bola abierta centrado en xo de radio r,
denotado por B(xo , r), esto es
B(xo , r) = {x ∈ RN : d(x, xo ) < r} = {x ∈ RN : |x − xo | < r}.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Definición (Conjunto abierto)
Se dice que un subconjunto Ω de RN es abierto si, dado
cualquier xo ∈ Ω, existe B(xo , r) totalmente contenida en Ω.

Definición (Conjunto cerrado)

Un subconjuto F de RN , se dice que, es cerrado si, su
complemento Fc es abierto, donde
Fc = RN − F = {x ∈ RN : x ∈
/ F}

Definición (Clausura de un conjunto)

Dado un conjunto A ⊆ RN , el menor conjunto cerrado que lo
contiene a A, se le llama la clausura de A, denotado
generalmente por A.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Definición (Interior de un conjunto)
El interior de un conjunto A, denotado por Å, es el mayor
conjunto abierto contenido en A.

Definición (Frontera de un conjunto)

El borde o frontera de un conjunto A, denotado por ∂A o F rA, esta
definida por
∂A = {x ∈ A : x ∈
/ Å}

Consideremos el abierto Ω ⊆ RN y la función escalar u : Ω −→ R,

tales que u = u(x), x = (x1 , x2 , . . . , xN ), con los conceptos de
continuidad y diferenciabilidad, las derivadas parciales de primer
orden de u respecto a una de las variables xi puden denotarse por
∂u
= ∂xi u = Dxi u = uxi
∂xi
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 En forma análoga, las derivadas parciales de segundo orden podrían
denotarse por
∂  ∂u  ∂2u
= = ∂x2 u = Di2 u = uxi xi
∂xi ∂xi ∂xi 2 i

∂  ∂u  ∂2u
= = ∂xj ∂xi u = Dj Di u = uxi xj
∂xj ∂xi ∂xj ∂xi

Ejemplo: Si u = u(x1 , x2 ) = −x2 cos x1 + ϕ(x2 ), entonces

∂u
= ux1 = −x2 (−sen x1 ) + 0 = x2 sen x1
∂x1
∂u
= ux2 = − cos x1 + ϕ′ (x2 )
∂x2
∂  ∂u 
= ux1 x1 = x2 cos x1
∂x1 ∂x1
∂  ∂u  ∂  ∂u 
= ux1 x2 = sen x1 ; = ux2 x1 = −(−sen x1 )
∂x2 ∂x1 ∂x1 ∂x2
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
¿Qué son las ecuaciones diferenciales parciales?

Definición (EDP)
Una ecuación diferencial parcial o ecuación en derivadas
parciales (EDP), es una ecuación de la forma
∂u ∂u ∂ 2 u ∂2u
 
F x1 , · · · , xN , u, ,··· , , ,··· , ,··· =0 (1)
∂x1 ∂xN ∂x21 ∂x1 ∂xN

que permite relacionar las variables independientes x1 , · · · , xN ,

la variable dependiente u y sus derivadas parciales. Donde
u : Ω ⊆ RN −→ R tal que u = u(x) con x = (x1 , · · · , xN ) ∈ Ω.

La EDP estará definida y planteada en la región abierta Ω ⊆ RN

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Definición
El Orden de una EDP está dada por mayor orden de las derivadas
presentes en la ED.

Definición
Se dice que una EDP es lineal si la función incógnita u y todas sus
derivadas parciales presentes en la ED son de primer grado. Caso
contrario se dice que la EDP es no lineal
La forma más general de una EDP lineal de primero y segundo orden
respectivamente son:
N
X ∂u
αi (x) + β(x)u + γ(x) = 0 ; donde x = (x1 , . . . , xN ) (2)
∂xi
i=1

N N
X ∂2u X ∂u
δij (x) + αi (x) + β(x)u + γ(x) = 0 (3)
∂xi ∂xj ∂xi
i,j=1 i=1

.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Para el caso particular de las EDPs (2) y (3), cuando
u = u(x) = u(x, y) se tienen:
∂u ∂u
α1 (x, y) + α2 (x, y) + β(x, y)u + γ(x, y) = 0 (4)
∂x ∂y

∂2u ∂2u ∂2u ∂u ∂u

δ1 (x) 2
+ 2δ2 (x) + δ3 (x) 2 + α1 (x) + α2 (x) + β(x)u + γ(x) = 0
∂x ∂x∂y ∂y ∂x ∂y
(5)

Definición
Se dice que una EDP lineal es homogenea si el término que no
contiene a la variable dependiente es identicamente nulo.
Por ejemplo, las EDPs (2), (3), (4) y (5) serán homogeneas si y sola-
mente si la función γ es identicamente nula (es decir, si γ(x) = 0).
Tener en cuenta que, una función identicamente nula es siempre solu-
ción de una EDP lineal homogenea.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 La parte de una EDP que contiene las derivadas de mayor orden
se llaman parte principal de la EDP. Por ejemplo, las partes
principales de las EDPs (4) y (5), respectivamente son

∂u ∂u
α1 (x, y) + α2 (x, y) (6)
∂x ∂y

∂2u ∂2u ∂2u

δ1 (x, y) + 2δ2 (x, y) + δ 3 (x, y) (7)
∂x2 ∂x∂y ∂y 2

Definición (EDP semilineal)

Se dice que una EDP no lineal es semilineal si su parte
principal es lineal.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Por ejemplo, las EDPs de la forma

∂u ∂u
α1 (x, y) + α2 (x, y) = ψ(x, y, u) (8)
∂x ∂y

son EDPs semilineales.

En particular:
∂u  x2 + y 2  ∂u
xy + + u2 = 0, es semilineal,
∂x 2 ∂y
∂u ∂2u
+ u 2 = 0, no es semilineal,
∂x ∂y
∂ 2 u ∂u
− = −sen u, es semilineal.
∂x2 ∂y
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
Ejemplos

Tabla para completar

Nombres EDPs EDPs Orden lineal homog Semi Cuasi
Laplace uxx + uyy = 0 2 Si Si No No
Calor ut = k2 uxx
Onda utt = k2 uxx
Burger ut + uux = kuxx
Poisson uxx + uyy = h(x, y)
Euler-Bernoulli utt + k2 uxxxx = 0
Schrodinger ut − ikuxx = 0
Helmholtz uxx + uyy + c2 u = 0
Tricomi uxx + xuyy = 0
Korleneg de Vries ut +k1 uux +k2 uxxx = 0
Sine - Gordon utt − uxx + sen u = 0
‡ xux − yuy = sen(xy)
† ux + uy + uz = 0
† xux + 2yuy = ku
† xuy + yux = u
† xux + yuy = 2
‡ uxx + uyy + utt = 0
‡ u2 2 2
x + uy + ut = f (x, y, z)

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
Solución de una EDP

Definición
Consideremos la EDP definida en (1) de orden k, se llama
solución de dicha EDP en cierta región Ω a una función escalar
cualquiera u : Ω ⊆ RN =⇒ R tal que al sustituir u y todas sus
derivadas parciales en dicha ecuación obtenemos una identidad.

Generalmente se exigirá que u = u(x1 , x2 , . . . , xN ) ∈ C k (Ω) sea continua o

diferenciable en ∂Ω, que tenga todas las derivadas parciales continuas en Ω
y que se cumpla la ecuación en el interior de Ω.
Siendo C k (Ω) el conjunto de funciones continuas en la región Ω junto con


todas las derivadas de hasta orden k .

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Definición
Una solución clásica de una EDP es una función lo
suficientemente regular que la satisface puntualmente.

Bajo esta teoría clásica se agrupan diversos métodos, que varían de acuerdo
con el tipo de la ecuación diferencial, y que guían los pasos a seguir para
encontrar una solución explícita.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Ejemplos #1
Una solución clásica a la ecuación del calor

∂u ∂2u
= (9)
∂t ∂x2

es la función u = u(t, x) definida sobre un dominio D ⊂ R2 tal que

todas las funciones u, ut , ux , utt , uxt , utx , uxx están bien definidas y
son continuas en cada punto (t, x) ∈ D, de modo que u ∈ C 2 (D), y,
además, (9) se cumple en cada (t, x) ∈ D.

Por ejemplo,
1 2
u= x +t (10)
2
es una solución de la ecuación del calor que está definida en todo el
dominio D = R2 , porque es C 2 , y además, ut = 1, ux = x y uxx = 1.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Ejemplo #2
Otra solución, más complicada pero extremadamente
importante, de la ecuación (9) es

2
e−x /4t
u= √ , (11)
2 πt

donde fácilmente se comprueba que u ∈ C 2 y resuelve la

ecuación del calor en el dominio D = {t > 0} ⊂ R2 .

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Ejemplo #3
La fución

√
u = e−t cos x + ie−t sin x donde i = −1 (12)

define una solución de valor complejo para la ecuación (9).

Esto se puede verificar directamente, ya que las reglas para

diferenciar exponenciales complejas son idénticas a las de sus
contrapartes reales.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
Solución de EDPs de primer y segundo orden
Consideremos las EDPs de primer y segundo orden respectiva-
mente  ∂u ∂u 
F x, y, u, , = 0.
∂x ∂y
 ∂u ∂u ∂ 2 u ∂ 2 u ∂ 2 u 
F x, y, u, , , , , = 0.
∂x ∂y ∂x2 ∂x∂y ∂y 2

La solución de estas EDP son las funciones escalares de la forma

u : Ω ⊆ R2 −→ R / u = u(x, y),
funciones que no contienen derivadas y satisfacen identicamente
a la ecuación diferencial.
Anotaremos algunos ejemplos atravez de las siguiente actividad
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Actividades

Example # 1
Verificar que las funciones u1 (x, t) = x2 + 2t y u2 (x, t) = e−t sen x son
soluciones de la EDP ut − uxx = 0.

Example # 2
Determine la relación entre a y b de tal modo que la función
u(x, y) = eay sen bx satisface la EDP uxx − uy = 0.

Example # 3
Compruebe que la función u(x, y) = ax3 + by 3 , donde a y b son
constantes, es la solución general de la EDP xux + yuy = 3u.
Además mensione dos soluciones particulares.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
EDPs que surgen de la eliminación de funciones
arbitrarias
Ya que las soluciones de las EDPs hacen intervenir funciones
arbitrarias, parece lógico que se obtengan EDPs por el proceso
inverso de eliminar tales funciones.
Example # 1
Encontrar una EDP de primer orden que tenga como solución general


la siguiente función u(x, y) = x − 6y − 3x3 f y siendo f diferenciable.

Example # 2
Encontrar una EDP de segundo orden que tenga como solución
general la siguiente función u(x, y) = xf (y) + yg(x) siendo f y g
funciones diferenciable.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Example # 3
Encontrar una EDP de primer orden que tenga como solución
y

general la siguiente función u(x, y) = xλ f 2 siendo f
x
diferenciable.

Example # 4
Encontrar una EDP de primer orden que tenga como solución
general la siguiente función u(x, y) = (x + y)f x2 − y 2 siendo f

diferenciable.

Example # 5
Compruebe que la función u(x, y, z) = h(x − y, y − z), con h
diferenciable es solución general de la EDP ux + uy + uz = 0.
Además mensione dos soluciones particulares.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Linealidad y Superposición .

Consideremos la EDP definida en (3) de orden k = 2. Observe

que, si k = 1, se tendría que δij = 0 cualesquiera que sean i, j ∈
{1, . . . , N } y existe i, 1 ≤ i ≤ N , tal que αi ̸= 0.
Reescribiendo la ecuación (3) de tal modo que

N N
X ∂2 X ∂
δij (x) u+ αi (x) u + β(x)u = −γ(x) (13)
i,j=1
∂xi ∂xj i=1
∂x i

Equivalentemente se tiene

N N
!
X ∂2 X ∂
δij (x) + αi (x) + β(x) u = f (x) (14)
i,j=1
∂xi ∂xj i=1 ∂xi
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 ∂ ∂2
Pero sabemos que = ∂xi y = ∂xj ∂xi = ∂x2i xj , entonces
∂xi ∂xj ∂xi
N N
X ∂2 X ∂
δij (x) + αi (x) + β(x) = L
∂xj ∂xi ∂xi
i,j=1 i=1

Reformulando la ecuación (14), tenemos

L(u) = f, (15)

siendo f (x) = −γ(x) y x = (x1 , . . . , xN ).

A cada función u (suficientemente diferenciable) le corresponde una única
función llamado Operador o Transformación.

L : C k (Ω) −→ C(Ω)
u 7−→ L(u)

donde el abierto Ω ⊂ RN , C(Ω) es el conjunto de funciones continuas y

C k (Ω) es el cunjunto de las funciones u : Ω −→ R k-veces continuamente
diferenciables.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Como en el caso de las EDOs, podemos asociarle a la EDP no homo-
genea (15) la EDP lineal homogenena asociada siguiente

Lu = 0. (16)

Por ejemplo:
Si L = ∂t2 − tx∂xt
2
, entonces la ecuación L(u) = sen x, expresado en su
∂2u ∂2u
forma estandar es utt − txuxt = sen x o 2
− tx = sen x .
∂t ∂t∂x
Si L = ∂t2 − t2 ∂xt
2
− y∂y , entonces la ecuación Lu = 0, expresado en su
∂2u ∂2u ∂u
forma estandar es utt − t2 uxt − yuy = 0 o 2
− t2 −y = 0.
∂t ∂t∂x ∂y
Ejercicio
N
X
Muestre que el operador L = αi (x)∂xi + β(x) es lineal.
i=1

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 ⋆
Consideremos un subconjunto abierto Ω de RN , las funciones reales de variable
vectorial u, v : Ω −→ R y el escalar λ ∈ R, donde u = u(x), v = v(x), β =
β(x), αi = αi (x), x = (x1 , · · · , xN ), entonces
N
X  N
X
L(u + λv) = αi ∂xi + β (u + λv) = αi ∂xi (u + λv) + β(u + λv)
i=1 i=1
∂ ∂ ∂
= α1 (u + λv) + α2 (u + λv) + · · · + αN (u + λv) + β(u + λv)
∂x1 ∂x 2 ∂x N
∂u ∂v ∂u ∂v
= α1 + λα1 + · · · + αN + λαN + βu + λβv
∂x1 ∂x 1 ∂x N ∂xN
∂u ∂u ∂v ∂v
= α1 + · · · + αN + βu + λα1 + · · · + λαN + λβv
∂x1 ∂xN ∂x 1 ∂xN
∂ ∂ ∂ ∂
   
= α1 + · · · + αN + β u + λ α1 + · · · + αN +β v
∂x 1 ∂xN ∂x 1 ∂x N
N
X  N
X 
= α i ∂x i + β u + λ αi ∂xi + β v = Lu + λLv
i=1 i=1

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Teorema 1.
Sea {ui } con 1 ≤ i ≤ r una familia de soluciones de la EDP lineal
r
X
homogenra (Lu = 0), entonces w = ci ui también es solución de la
i=1
misma, siendo los ci constantes y L un operador diferencial parcial lineal.

Demostración.
Como los ui son soluciones de Lu = 0 entonces Lui = 0, para todo
i = 1, . . . , r. Luego,
r
X
Lw = L ci ui = L(c1 u1 + c2 u2 + · · · + cr ur )
i=1
= L(c1 u1 ) + L(c2 u2 ) + · · · + L(cr ur )
= c1 L(u1 ) + c2 L(u2 ) + · · · + cr L(ur )
= c1 (0) + c2 (0) + · · · + cr (0)
= 0 + 0 + ··· + 0
= 0.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Teorema 2.
Sea up una solución particular de la EDP lineal no homogenea Lu = f ,
entonces la solución general da la EDP es u = up + ug , donde ug es la
solución general de la EDP homogenea Lu = 0.

Demostración.
1 Primero veamos que up + ug es solución de Lu = f .
En efecto: Tenemos que up es solución de Lu = f , esto es, Lup = f .
Además ug es solución de Lu = 0, es decir, Lug = 0.
Luego, L(up + ug ) = L(up ) + L(ug ) = f + 0 = f
∴ up + ug es solución de Lu = f .
2 Ahora, probemos que si u es solución de Lu = f entonces u es de la
forma up + ug , para alguna solución ug de Lu = 0.
En efecto: Sea u = up + v, de donde v = u − up entonces
Lv = L(u − up ) = Lu − Lup = f − f = 0. Esto prueba que, v es una
solución de Lu = 0. ∴ u siempre tendrá la forma up + ug .
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Proposición 1. (Principio de la Superposición)
Si L es un operador diferencial parcial lineal y {ui }, con 1 ≤ i ≤ r es un
conjunto de funciones que satisfacen las ecuaciones Lui = fi , con
r
X r
X
i = 1, . . . , r, entonces w = ci fi es solución de Lu = f , donde f = fi .
i=1 i=1

¿ Existe una forma infinita del Principio de Superposición?

Proposición 2. (Principio de la Superposición)

Sea L un operador diferencial lineal de orden k. Supongamos que un n≥1
,
k
 C (Ω) que satisfacen la EDP lineal
es una familia de funciones de clase
homogenea Lu = 0. Entonces si cn n≥1 es una sucesión de escalares tal
que
∞
X
u(x) = cn un (x)
n=1

es convergente y k-veces diferenciable término a término en Ω, u satisface

la ecuación Lu = 0, donde x = (x1 , . . . , xN ) ∈ Ω.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Demostración. (Prop. #1)
La demostración que como ejercicio para el estudiante.

Demostración. (Prop. #2)

Para la demostración consideremos k = 2 (orden del operador (??)-(15)).
En este caso, por hipótesis, cualesquiera que sean x ∈ Ω; 1 ≤ i, j ≤ N , se
tienen:
∞
X
u(x) = cn un (x) = c1 u1 (x) + c2 u2 (x) + · · · + cr ur (x) + · · ·
n=1
∞
X
∂xi u(x) = c1 ∂xi u1 (x) + c2 ∂xi u2 (x) + · · · = cn ∂xi un (x)
n=1
∞
X
∂xj ∂xi u(x) = c1 ∂xj ∂xi u1 (x) + c2 ∂xj ∂xi u2 (x) + · · · = cn ∂xj ∂xi un (x)
n=1
Observamos que estas series convergen.Por lo tanto, usando (??)-(15)
N N
X X
Lu(x) = δij (x)∂xj ∂xi u(x) + αi (x)∂xi u(x) + β(x)u(x)
i,j=1 i=1
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 P.P.S:
N ∞ N ∞
X X X X
Lu(x) = δij (x) cn ∂xj ∂xi un (x) + αi (x) cn ∂xi un (x)u(x)+
i,j=1 n=1 i=1 n=1
∞
X
β(x) cn un (x)
n=1
∞ N
X N 
X X
= cn δij (x)∂xj ∂xi un (x) + αi (x)∂xi un (x)u(x) + β(x)un (x)
n=1 i,j=1 i=1
∞
X

= cn Lun (x)
n=1
∞
X

= cn 0
n=1
=0.

Lo que demuestra la proposición 2. para el caso k = 2.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Operador Laplaciano
Definición (Gradiente o Vector Gradiente).
El gradiente de u (o vector gradiente u : Ω ⊆ RN −→ R) denotado ∇u
es un campo vectoria ∇u : Ω ⊆ RN −→ RN . El vector gradiente de u
evaluado en un punto genérico x ∈ Ω, definido por
 ∂u(x) ∂u(x)

∇u(x) = ,..., , indica la dirección en la cual el campo u
∂x1 ∂xN
varía más rápidamente y su módulo representa el ritmo de variación de u en
la dirección de dicho vector.

Definición (Laplaciano).
El Operador Laplaciano (o Laplaciano), es un operador diferencial
elíptico de segundo orden, denotado como ∆ o ∇2 . En coordenadas
cartesianas, si u : Ω ⊆ RN −→ R, entonces el Laplaciano de éste es,
N N
∂2u ∂2u ∂2u X ∂2u X 2
∆u = ∇ · ∇u = ∇2 u = 2
+ 2
+ ··· + 2
= = ∂xi u
∂x1 ∂x2 ∂xN ∂x2i
i=1 i=1

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Éste operador tiene ese nombre en reconocimiento a Pierre-Simon Laplace
que estudió soluciones de ecuaciones diferenciales en derivadas parciales en
las que aparecía dicho operador.
Por ejemplo:
la EDP uxx + uyy + uzz = 0 es equivalente a expresarlo como ∆u = 0,
denominado ecuación de Laplace en tres variables independientes,
∂2u ∂2u ∂2u
donde ∆u = 2
+ 2
+ .
∂x ∂y ∂z 2
∂u
la EDP: = c2 ∆u, es la ecuación del calor en dimensiones mayores,
∂t
donde u = u(x), x ∈ Ω ⊆ RN , t > 0, ∆ es el Laplaciano en RN .

Ejercicio
Demuestre que el operador laplaciano ∆ en las variables espaciales
x, y, z, es un operador diferencial parcial lineal.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Condiciones de frontera
Considerando la EDP (1) sobre un dominio Ω ⊆ RN , normalmente se nos
pide no sólo que encontremos una función u = u(x) que satisfaga dicha
EDP en todo punto de Ω, sino también que tal función cunpla una serie de
condiciones

∂u ∂u ∂ 2 u ∂2u
 
µi x, u, ,··· , , 2,··· , ,··· = 0, x ∈ ∂i , i = 1, . . . , r
∂x1 ∂xN ∂x1 ∂x1 ∂xN

donde los símbolos ∂i denotan partes de ∂Ω, y los µi son funciones dependi-
entes de las variables xi .
Condiciones de este tipo se denominan condiciones de contorno (o de
frontera).
Un problema consistente en resolver una EDP sobre un dominio Ω y un
conjunto de condiciones de contorno se denomina problema de contorno
(o problema de frontera)

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 ⋆
Ejemplo:
Se tiene la ecuación de Laplace ∆u = 0 en una caja rectangular
Ω = ⟨0, a⟩ × ⟨0, b⟩ × ⟨0, c⟩
Donde las condiciones de contorno son

u(x, y, c) = g(x, y) en Ω
u(x, y, 0) = 0, en ∂Ω
u(0, y, z) = 0, en ∂Ω
u(a, y, z) = 0, en ∂Ω
u(x, 0, z) = 0, en ∂Ω
u(x, b, z) = 0, en ∂Ω

La solución u = u(x, y, z) a este problema modela el potencial electrostático

en una caja con todas sus caras a potencial cero, mientras que la primera
condición tiene el potencial g(x, y).
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Actividad
Ejercicio # 1
Resolver el problema siguiente
∆u = 0 ; en Ω = ⟨0, 4⟩ × ⟨0, 4⟩
Donde las condiciones de contorno son

u(x, 4) = 4x − x2 ; 0 < x < 4

u(x, 0) = 0 ; 0<x<4
u(4, y) = 0 ; 0<y<4
u(0, y) = 0 ; 0<y<4

El problema tambien puede expresarse en la forma

⟨0, 4⟩ × ⟨0, 4⟩

 uxx + uyy = 0,
u(x, 4) = 4x − x2 ,



 0<x<4
u(x, 0) = 0, 0<x<4



 u(4, y) = 0, 0<y<4

u(0, y) = 0, 0<y<4
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Condiciones de Robin, Dirichlet & Neumann
Las condiciones de frontera aparecen de manera natural en la descripción de fenomenos
físicos estacionarios (sus características no varian con el tiempo) muchas veces en-
contramos condiciones del tipo

∂u
αu + β = f (x) sobre ∂Ω (17)
∂n

para algunas constantes no nulas α , β y una función dada f definida sobre ∂Ω.
∂u
Aquí, u es la solución desconocida definida sobre el conjunto abierto Ω, := n·∇u
∂n
es la derivada normal en la frontera ∂Ω y n es el vector normal unitario.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 † La condición (17) es llamado condición de frontera de Robin (o de
tercer tipo) denominado así, en honor a Victor Gustave Robin (1855-1897).
Por ejemplo, en una dimensión, si, Ω = [0, 1], la condición de Robin es:
αu(0) − βu′ (0) = f (0)


αu(1) + βu′ (1) = f (1)

donde se puede observar el cambio de signo en el frente que involucra la
derivada, esto es porque la normal a [0, 1] en 0 apunta en la dirección negativa,
mientras que en 1 apunta en dirección positiva.

† Si, en la condición (17) se tiene que solo β = 0, dicha condición se conoce

con el nombre de condición de Dirichlet (o de primer tipo), denominado
así en honor a Johann Peter Gustav Lejeune Dirichlet (1805-1859).
Por ejemplo, para una EDP ∆u + u = 0 sobre un dominio Ω ⊂ RN , las
condiciones de frontera de Dirichlet toman la forma: u(x) = f (x), ∀x ∈ ∂Ω.

† Si solo α = 0, en la condición (17) se tendrá la condición de frontera

de Neumann (o de segundo tipo), llamada así en alusión a Carl Gottfried
Neumann.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
Condición de frontera mixta
Una condición de frontera mixta para una EDP indica que se utilizan difer-
entes condiciones de frontera sobre partes diferentes de la frontera del dominio
de la ecuación.

Por ejemplo, si u es una solución a una EDP sobre el conjunto abierto Ω

con frontera ∂Ω suave a tramos, y ∂Ω está dividida en dos partes, Γ1 y Γ2 ,
una puede usar la condición de frontera de Dirichlet sobre Γ1 y una condición
de frontera de Neumann sobre Γ2 :

u = uo
Γ1




 ∂u

=g

∂n Γ2

donde u0 y g son funciones dadas definida sobre aquellas porciones de la

frontera.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Como la condición de frontera de Robin es una combinación lineal de las condiciones

de frontera de Dirichlet y Neumann, es otro tipo de condición de frontera híbrida.

Verde: condición de frontera de Neumann;

Púrpura: condición de frontera de Dirichlet.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Condiciones iniciales
Otro tipo de condiciones que suelen exigise a las soluciones de una EDP son
las denominadas condiciones iniciales respecto de una de las variables.
Como en las EDPs tenemos más de una variable dependiente (por ejemplo x y
t), es natural fijar una de las variables (por ejemplo t = 0) e imponer el valor
de la solución y de sus derivadas parciales respecto a la variable fijada como
función de las otras variables (por ejemplo u(x, 0) = f (x) y ut (x, 0) = g(x),
donde f y g serán dadas).
Observese que en el caso n = 2, con variables x y t, fijar t = 0 significa
inponer el valor de la solución y sus derivadas normales a lo largo de la curva
inicial t = 0.
De forma similar, en el caso n = 3, con variables x, y y t, fijar t = 0 significa
observar la solución (y sus derivadas normales si fuera el caso) a lo largo de
la superficie inicial t = 0. En estos casos el problema corrrespondiente es
denominado problema de Cauchy o problema de valor inicial.
En general, las condiciones iniciales no son condiciones de contorno ya que
también se consideran situaciones en las que el conjunto determinado por la
relación t = to puede estar en el interior de Ω.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Ejemplo:
Consideremos la ecuación del calor ut = uxx , en 2 dimensiones,
sobre el dominio Ω = {(x, t) ∈ R2 / t > 0 ; 1 < x < 2}.
Donde las condiciones iniciales son

u(x, 0) = (x − 1)(x − 2) en Ω
u(1, t) = 0, en ∂Ω
u(2, t) = 0, en ∂Ω

En este caso la condición inicial es también condición de contorno.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
Observación.
En problemas mixtos, la condición de contorno y la condición inicial no son
completamente independientes y por ende, es preciso que se satisfaga una
condición de compatibilidad para que exista solución.
Por ejemplo, El problema para la ecuación de onda en un intervalo acotado
utt = k2 uxx ,

 en Ω = ⟨0, ℓ⟩ × ⟨0, +∞⟩




 u(0, t) = 0, t≥0
u(ℓ, t) = 0, t≥0



 u(x, 0) = f (x), 0≤x≤ℓ


ut (x, 0) = g(x), 0≤x≤ℓ
también puede ser considerado como un problema mixto donde u(0, t) = 0 =
u(ℓ, t), t ≥ 0, son las condiciones de frontera y u(x, 0) = f (x) y u(x, 0) =
g(x), x ∈ [0, ℓ] son las condiciones iniciales. O como un problema de frontera
en la región no acotada [0, ℓ] × [0 + ∞⟩.
Para que exista solución en este problema es necesario que f satisfaga la
condición de compatibilidad f (0) = 0 = f (ℓ).
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 La generalización del problema anterior con dos variables espaciales estará
dada por

utt = k2 ∆u,


 en Ω × ⟨0, +∞⟩

 u(x, y, t) = 0, (x, y) ∈ ∂Ω, t ≥ 0
,

 u(x, y, 0) = f (x, y), (x, y) ∈ Ω

(x, y) ∈ Ω

ut (x, y, 0) = g(x, y),

donde Ω es un subconjunto abierto de R2 y ∆ es el laplaciano en R2 y la

condición de compatibilidad en este caso es f (x, y) = 0, ∀(x, y) ∈ ∂Ω.

Este problema describe el movimiento de una membrana vibrando (como en

un tambor) y también tiene una única solución si f satisface la condición de
compatibilidad.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Observación Importante.
Dado un problema que consiste en una EDP y sus condiciones de frontera
y/o iniciales, existen tres preguntas fundamentales.
1 existencia de soluciones,
2 unicidad de soluciones,
3 dependencia de la solución en los datos iniciales y/o de frontera.

Para estudiar la existencia de soluciones, es necesario especificar no

solamente la clase de funciones donde buscamos solución sino también
en que sentido las condiciones de frontera y/o iniciales son satisfechas.
Una vez obtenida la existencia, el significada de la unicidad es claro:
deseamos saber si la solución es única dentro de la clase especificada.
El estudio de la dependencia de la solución en los datos iniciales y/o
frontera es muy importante. Es decir, la solución depente
continuamente de los datos iniciales y/o frontera.
Un problema para el cual se dan la existencia, unicidad y dependencia con-
tinua de los datos iniciales y/o frontera se llama un problema bien puesto
(en el sentido de Hadamard). De lo contrario se dice que es un problema
mal puesto.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Actividad

Ejercicio # 1
Supongamos que se quiere obtener la temperatura u(x, y) correspondiente
al estado estacionario de una placa rectangular.
Ω = {(x, y) ∈ R2 / 0 < x < a ∧ 0 < y < b}.
Cuando no se pierde  calor a través de las caras laterales, el problema es.

 ∆u = 0, en Ω

u(x, b) = φ(x), 0<x<a




u(x, 0) = 0, 0<x<a

ux (0, y) = 0, 0<y<b





ux (a, y) = 0, 0<y<b


Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Ejercicio # 2
La siguiente EDP modela la variación de temperatura con el tiempo en una
varilla de longitud 1. Si u(t, x) es la temperatura de la varilla en el instante
t en la posición x de la varilla, entonces, podemos plantear el problema
como sigue 

 uxx = ut , 0<x<1 ; t>0


 u(0, x) = f (x), 0≤x≤1


 u(t, 0) = 0, t≥0

t≥0

u(t, 1) = 0,
La condición u(0, x)=f (x) es la temperatura en el instante inicial t = 0 y
representa una condición inicial del problema. Las condiciones u(t, 0) = 0 y
u(t, 1) = 0 son restricciones que se le imponen a los extremos de la varilla
en cualquier instante de tiempo y se denominan condiciones de contorno.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
Resolución.
Nos piden la temperatura u(t, x) de la varilla en el instante t en la posición
x de la varilla, donde
el problema está modelado por:

 uxx = ut , 0<x<1 ; t>0


0≤x≤1

 u(0, x) = f (x),

 u(t, 0) = 0,

 t≥0


u(t, 1) = 0, t≥0


La condición inicial u(0, x)=f (x) indica la distribución de temperatura en

la varilla en el instante inicial. Las condiciones de contorno u(t, 0) = 0
y u(t, 1) = 0 indican que la temperatura en los extremos de la varilla es
constante e igual a 0.
Supongamos que u(t, x) = g(t)h(x), obviamente la función nula (u = 0) es
solución de la EDP, sin embargo sólo será solución del problema si f (x) = 0
que es el caso trivial.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 Busquemos soluciones no triviales, para esto supondremos que h(x) ̸= 0 y
g(t) ̸= 0.
Tomando las condiciones de frontera, se tienen:
u(t, 0) = g(t)h(0) = 0 =⇒ h(0) = 0 y u(t, 1) = g(t)h(1) = 0 =⇒ h(1) = 0.
Además, si u(t, x) = g(t)h(x), entonces uxx = g(t)h′′ (x) y ut = g ′ (t)h(x).
h′′ (x) g ′ (t)
PPS: g(t)h′′ (x) = g ′ (t)h(x) ⇐⇒ =
h(x) g(t)
Observamos que, un lado de la igualdad depende sólo de x y el otro depende
sólo de t, por tanto, para que se dé la igualdad ambos miembros deben ser
constantes, es decir
h′′ (x) g ′ (t)
= = k, con k ∈ R
h(x) g(t)
′′ ′
De aquí, h (x) − kh(x) = 0 · · · (α) y g (t) − kg(t) = 0 · · · (β)
⋆ Si k = 0, en (α) y (β) respectivamente se tienen h′′ (x) = 0 y g ′ (t) = 0 cuya
solución general se obtiene integrando h(x) = c1 x + c2 y g(t) = c3
Luego, u(t, x) = g(t)h(x) = c3 (c1 x + c2 ) = Ax + B
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
⋆

Como u(t, 0) = 0 = A(0) + B, entonces B = 0, quedando u(t, x) = Ax,

además u(t, 1) = 0 = A(1), entonces A = 0, así u(t, x) = 0.

⋆ Si k > 0, podemos suponer que k = λ2 > 0; λ ̸= 0.

En (α) y (β) respectivamente se tienen h′′ (x)−λ2 h(x) = 0 y g ′ (t)−λ2 g(t) = 0.

2
cuyas soluciones son: h(x) = c4 cosh λx + c5 senh λx y g(t) = c6 eλ t

Como h(0) = 0 = h(1), es decir c4 cosh λ0 + c5 senh λ0 = 0, entonces c4 = 0,

quedando h(x) = c5 senh λx.

Además h(1) = 0, esto es c5 senh λ = 0, de aquí c5 ̸= 0, entonces senh λ = 0

pero, senh λ = 0 ⇐⇒ λ = 0 lo cual contradice a lo supuesto (que λ ̸= 0).

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
⋆

⋆ Si k < 0, podemos suponer que k = −λ2 < 0; λ ̸= 0.

En (α) y (β) respectivamente se tienen h′′ (x)+λ2 h(x) = 0 y g ′ (t)+λ2 g(t) = 0.

2
cuyas soluciones son: h(x) = c7 cos λx + c8 senλx y g(t) = c9 e−λ t

Como h(0) = 0, es decir c7 cos λ0 + c8 senλ0 = 0, entonces c7 = 0, quedando

h(x) = c8 senλx.

Además, h(1) = 0, es decir c8 senλ = 0, ya que c8 ̸= 0, se tiene senλ = 0.

senλ = 0 si solo si λ = nπ para todo n ∈ Z ∖ {0}.

2
Expresemos hn (x) = cn sen(nπx) y gn (t) = kn e−(nπ) t , para n = 1, 2, . . .
2 2
Luego, un (t, x) = kn e−(nπ) t cn sen(nπx) = (kn cn ) e−(nπ) t sen(nπx)

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 ⋆
Cualquier combinación lineal de estas funciones, es también solución (Por el
P∞
Principio de superposición u = n=1
an un ).

∞ ∞
X 2 X 2
u(t, x) = an (kn cn ) e−(nπ) t sen(nπx) = bn e−(nπ) t sen(nπx)
n=1 n=1

. Usado la condición inicial, u(0, x) = f (x), es decir

∞
X 2
bn e−(nπ) 0 sen(nπx) = f (x)
n=1

. ∞
X Z 1
f (x) = bn sen(nπx) ⇝⇝ bn = 2 f (x)sen(nπx)dx
n=1 0

.
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
Condición de frontera de Cauchy
Las condiciones de frontera de Cauchy en EDOs o en EDPs imponen valores
específicos a la solución de una ecuación diferencial que se toma de la frontera
del dominio y de la derivada normal a la frontera.

Esto es igual a imponer dos tipos de condiciones: la condición de frontera de

Dirichlet y la condición de frontera de Neumann. Su nombre hace honor al
prolífero matemático francés del siglo XIX Augustin Louis Cauchy.

Las condiciones de Cauchy son también llamadas condiciones de valor inicial

o simplemente valores de Cauchy.

Las condiciones de frontera de Cauchy son una generalización de estos tipos

de condiciones.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Definamos una EDP de segundo orden:

∆u = u(x, y)

Tenemos un dominio de dos dimensiones cuya frontera es una línea, que

puede ser descrita por las siguientes ecuaciones paramétricas:
(
x = ξ(s)
y = η(s)
así, de manera similar que en las ecuaciones diferenciales ordinarias de se-
gundo orden, necesitamos ahora conocer el valor de la función en la frontera
y su derivada normal a esta para hallar la solución de la ecuación diferencial
parcial, es decir, que:
 u(s)
 ∂u = n · ∇u(s)
∂n

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Observe que si bien las condiciones de frontera de Cauchy implican tener
tanto las condiciones de frontera de Dirichlet como las de Neumann, que no
es lo mismo del todo que tener una condición de frontera de Robin.

Una mezcla entre las condiciones de frontera de Dirichlet y Neumann está

dada por
∂u
a(s)u(s) + b(s) (s) = f (s) sobre ∂Ω
∂n
donde se entiende que a(s), b(s) y f (s) deben darse en la frontera (esto
contrasta con el término condiciones de frontera mixtas, que es generalmente
usado para referirse a condiciones de frontera de tipos diferentes en subgrupos
distintos de frontera).

En este caso la función y sus derivadas deben cumplir una condición dentro
de la misma ecuación de la condición de frontera.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
.

El problema

 uy = 0, en R2
 u(x, p(x)) = f (x), x ∈ R,

es un problema de Cauchy, donde p, f ∈ C 1 (R) son funciones dadas. Luego,

la función u se busca a lo largo de la curva inicial y = p(x) en el plano.
Como uy = 0, tenemos que u(x, y) = g(x), ∀(x, y) ∈ R2 donde g ∈ C 1 (R) es
arbitrario y u ∈ C 1 (R2 ).
Puesto que u(x, p(x)) = f (x), entonces g(x) = f (x), es decir, la solución del
problema es u(x, y) = f (x), ∀(x, y) ∈ R2 .

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Se ha probado entonces que, si u es una solución clasica del problema dado,
dicha solución está dada por u(x, y) = f (x), ∀(x, y) ∈ R2 . Recíprocamente,
si u está dada por u(x, y) = f (x), ∀(x, y) ∈ R2 , entonces u es solución del
problema dado. Por lo tanto, el problema dado tiene solución clásica única.

Observaciones:

Tenemos dependencia continua en los datos del problema dado, puesto

que la solución u(x, y) es igual a la condición inicial f (x). Por lo
tanto, el problema dado es bien puesto.

Si, en el problema dado, se cambia la condición u(x, p(x)) = f (x) por

u(0, y) = f (y), y ∈ R, el problema sigue siendo de Cauchy, la curva
inicial es el eje de las ordenadas. Así, el problema no tiene solución (si
f no es constante) o tiene infinitas soluciones (si f es constante).

Además, así propuesto el problema no es bien puesto, puesto que

ux (0, y) = f ′ (y), y por tanto f ′ (y) = 0, ∀y ∈ R.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Luego, para que exista solución, es preciso que f sea constante. Si f = k es
constante, como la solución general de la EDP es u(x, y) = g(x), ∀(x, y) ∈ R2 ,
donde g ∈ C 1 (R) es cualquier función que satisface g(0) = k.
Los problemas que involucren EDPs no homogeneas se resuelven de la
misma manera.
Por ejemplo, dado el problema

 uy = φ(x, y), (x, y) ∈ R2
 u(x, p(x)) = f (x), x ∈ R,

donde φ ∈ C(R2 ), p, f ∈ C 1 (R) serán conocidos, las curvas

características planas son las mismas, a lo largo de la recta x = xo .
la ecuación queda
d
u(xo , y) = φ(xo , y)
dy

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Para hallar la solución en el punto (xo , yo ) basta integral a lo largo del seg-
mento de recta determinado pot (xo , p(xo )) y (xo , yo ),

Por lo tanto:
Z yo Z yo
u(xo , yo ) = u(xo , p(xo )) + uy (xo , s)ds = f (xo ) + φ(xo , s)ds
p(xo ) p(xo )

osea, Z y
u(x, y) = f (x) + φ(x, s)ds, (18)
p(x)

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
Ejercicios:
Verificar si los siguientes problemas tienen solución si fueran así presentare-
mos dicha solución.
(
uy = xey , ∀(x, y) ∈ R2
1 2
u(y 2 , y) = ey + y 4 , y ∈ R.

uy = xey , ∀(x, y) ∈ R2
2
u(y 2 , y) = y 2 ey , y ∈ R.

ux = 2xy, ∀(x, y) ∈ R2
3
u(ey , y) = y 2 + 1, y ∈ R.

ux = 2xy, ∀(x, y) ∈ R2
4
u(x, x2 ) = 1, x ∈ R.

ux = 2xy, ∀(x, y) ∈ R2
5
u(x, x2 ) = x, x ∈ R.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
Método de las curvas características planas
Primero recordemos que, las soluciones generales dependen de una función
arbitraria.
1 si ux = 0 y uy = 0 entonces respectivamente se tienen las soluciones
u = ψ(y) y u = ϕ(x) con ψ, ϕ funciones arbitrarias.
2 Para la EDP ux = uy , el cambio de variables α = x + y y β = x − y
conduce a
α+β α−β
u(x, y) = u( , ) = w(α, β)
2 2
La EDP se transforma entonces en wβ = 0, cuya solución es,
w = f (α) = f (x + y), siendo f una función arbitraria.
3 La EDP que describe el movimiento ondulatorio unidireccional es

ut + cux = 0, (19)

donde la velocidad de propagación c es una constante.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Consideremos además la condición inicial u(x, 0) = f (x), siendo f una fun-
ción conocida.
La ecuación (17) representa la derivada de u en la dirección s dada por
 dt
 ds = 1,


(20)
 dx = c.


ds
dx
Además, = c, de donde se obtiene la familia de curvas (en este caso
dt
rectas) en el plano de las fases (x, t), descriptas por la ecuación x = ct + k.

El sistema (22) determina una familia de curvas en el plano xt. De modo

que t = t(s) y x = x(s).


A lo largo de estas curvas la función u = u(x, t) es u = u x(s), t(s) y verifica

du dx dt
= ux + ut = cux + ut = 0.
ds ds ds
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
 du
Como = 0, entonces se tiene que
ds
u=C (21)

con C constante (relativa) sobre cada recta x = ct + k, es decir u = C(k).

Una solución general entonces es

u = C(x − ct) (22)

y el problema quedará resuelto determinando C mediante la condición inicial:

u(x, 0) = f (x), es decir C(x − c0) = f (x), quedando f (x) = C(x), de donde
se obtiene
u(x, t) = f (x − ct) (23)

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Consideremos la EDP lineal no homogenea de primer orden

α(x, y)ux + β(x, y)uy + ψ(x, y)u = γ(x, y) (24)

La parte principal de la ecuación, representa la derivada de u en la

β(x, y)
dirección , en los puntos en los cuales no son simultáneamente
α(x, y)
nulos los coeficientes α(x, y) y β(x, y).

Si se consideran en el plano XY las curvas definidas por x = x(s) y

β(x, y) dy
y = y(s), cuyas tangentes tienen la dirección = de modo que
α(x, y) dx
dx dy
= α(x, y) y = β(x, y).
ds ds
Además
dx dy
= (25)
α(x, y) β(x, y)

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 du dx dy
Si u(x, y) = u(x(s), y(s)), entonces = ux + uy
ds ds ds
du dx dy
⋆ = ux + uy = ux α(x, y) + uy β(x, y) = γ(x, y) − ψ(x, y)u
ds ds ds
du dy β(x, y) ux α(x, y) + uy β(x, y)
⋆ = ux + uy = ux + uy =
dx dx α(x, y) α(x, y)
du γ(x, y) − ψ(x, y)u
=
dx α(x, y)
Quedando
du dx
= (26)
γ(x, y) − ψ(x, y)u α(x, y)
Finalmente, de (25) y (26) se tiene que

dx dy du
= = (27)
α(x, y) β(x, y) γ(x, y) − ψ(x, y)u

La relación (27) es llamado sistema auxiliar de Lagrange.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 Si en la ecuación (24) la función ψ es identicamente nula, queda

α(x, y)ux + β(x, y)uy = γ(x, y) (28)

el sistema auxiliar de Lagrange será

dx dy du
= = (29)
α(x, y) β(x, y) γ(x, y)

Si en la ecuación (24) la función γ es identicamente nula, queda una

EDP lineal homogénea

α(x, y)ux + β(x, y)uy + ψ(x, y)u = 0 (30)

el sistema auxiliar de Lagrange será

dx dy du
= = (31)
α(x, y) β(x, y) −ψ(x, y)u
Coaquira Cárdenas Víctor Alcides UNSCH
Universidad Nacional de San Cristóbal de Huamanga
Actividades
Ejercicio # 1
hallaremos una integral general de la EDP

∂u ∂u
xy − xy = (x − y)u .
∂x ∂y

Ejercicio # 2
Resolver el problema de valor inicial
(
yux + xuy − 2u = 0,
u(0, 1) = 1.

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
Resolución ejercicio #1

Consideremos

xyux − xyuy + (y − x)u = 0.

Por el sistema auxiliar de Lagrange se tiene

dx dy du
= =
xy −xy −(y − x)u

dx dy
i) = ⇝ dx + dy = 0 ⇝ x + y = c ⇝ x = c − y
xy −xy

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 ⋆

dy du dy du
ii) = ⇝ =
−xy −(y − x)u (c − y)y (2y − c)u
(2y − c)dy du du dy dy
= ⇝ + + =0
(c − y)y u u (y − c) y
Integrando

ln u + ln(y − c) + ln y = ln k1 ⇝ u(y − c)y = k1 ⇝ u(−x)y = k1

k
Despejando u(x, y) =
xy

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
Resolución ejercicio #2

Tenemos la EDP yux + xuy − 2u = 0.

Por el sistema auxiliar de Lagrange se tiene

dx dy du dx dy du
= = ⇝ = =
y x −(−2)u y x 2u

dx dy y 2 x2
i) = ⇝ y dy − x dx = 0 ⇝ − = c1 ⇝ y 2 − x2 = c
y x √ 2 2
consideremos y = c + x2 > 0 por la condición u(0, 1) = 1.
du dx du dx du 2dx
Z Z
ii) = ⇝ =√ ⇝ = √
2u y 2u c+x 2 u c + x2

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 ⋆
√ √ √
Si x = c tan θ ⇝ x2 + c = c sec2 θ ⇝ x2 + c = c sec θ
√
Además, dx √ = c sec2 θdθ , entonces
du 2 c sec2 θdθ
Z Z Z
= √ = 2 sec θdθ
u c sec θ
Integrando
ln |u| = 2 ln | sec θ +tan θ|+ln |k1 | ⇝ ln |u| = ln |(sec θ +tan θ)2 k1 |

u = ±(sec θ + 2
√tan θ) k1
 c + x2 x 2 k
u(x, y) = k √ +√ = (y + x)2
c c c
La condición inicial u(0, 1) = 1 determinará la constante
k
(1 + 0)2 = 1 ⇝ k = c .
c
Finalmente, u(x, y) = (x + y)2

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga
 ⋆

Ejercicio para Exposición

Resolver el problema de valor inicial

 (xy 3 − 2x4 )ux + (2y 4 − x3 y)uy + 9(y 3 − x3 )u = 0,



1
 u(x, ) = x3 .
x

Coaquira Cárdenas Víctor Alcides UNSCH

Universidad Nacional de San Cristóbal de Huamanga