Instituto Politécnico Nacional.

Escuela Superior de Fı́sica y Matemáticas del IPN

Notas de Ecuaciones diferenciales parciales

José Santos Ordoñez

6 de febrero de 2026

1. Series de Fourier
1.1. Relaciones de ortogonalidad de las funciones trigonométricas.
Las siguientes relaciones integrales son de gran utilidad en la solución de los problemas que involucran
soluciones de EDP’s.

Teorema 1.1

Z +L
nπ mπ
sin( x) sin( x)dx = Lδnm (1.1)
−L L L
Z +L
nπ mπ
cos( x) cos( x)dx = Lδnm (1.2)
−L L L
Z +L
nπ mπ
sin( x) cos( x)dx = 0 ∀n, m ∈ Z (1.3)
−L L L

donde δnm es la delta de Kronecker.

Notando que los productos de senos y de cosenos generan funciones par entonces se pueden reescribir como
sigue.

Teorema 1.2

Z +L
nπ mπ L
sin(
x) sin( x)dx = δnm (1.4)
0 L L 2
Z +L
nπ mπ L
cos( x) cos( x)dx = δnm (1.5)
0 L L 2

1.2. Series de Fourier para periodos arbitrarios.

Aunque se estudia primero el caso de funciones 2π-periódicas, este es un caso particular de las funciones
2p-periódicas.

1
 Definición 1.1: Serie de Fourier
Se conoce como serie de Fourier una serie trigonométrica de la forma:
∞
X nπ nπ
a0 + an cos( x) + bn sin( x) (1.6)
p p
n=1

El interés está en poder expresar a un función f (x) en forma de serie de Fourier, lo que ocurre y se establece
en el siguiente resultado.

Teorema 1.3
Sea f (x) una función suave por trozos. La serie de Fourier de f (x) está dada por
∞
X nπ nπ
a0 + an cos( x) + bn sin( x) (1.7)
p p
n=1

donde los coeficientes están dados por:

Z p
1
a0 = f (x)dx (1.8)
2p −p
1 p
Z
nπ
an = f (x) cos( x)dx (1.9)
p −p p
Z p
1 nπ
bn = f (x) cos( x)dx (1.10)
p −p p

con n = 1, 2, 3... La serie trigonométrica converge a f (x) si f es continua en x. Para cualquier otro caso
la serie converge al promedio

f (x+) + f (x−)
(1.11)
2

Algunos datos aspectos importantes sobre las series de Fourier son:

1. La serie de Fourier de una función, es en sı́ misma una función, por lo que la obedece las operaciones de
funciones conocidas.
2. Se pueden obtener series de Fourier a partir de una serie conocida por medio de cambios de variables o
con las operaciones entre funciones.
3. Al manipular series de Fourier se debe cuidar como puede verse alterado el intervalo de validez. Éste puede
cambiar al manipular las series.

1.3. Series de funciones Par e Impar

La paridad de una función es de gran utilidad porque simplifica considerablemente el cálculo de los coefi-
cientes de una serie de Fourier.

Definición 1.2: Función Par e Impar

1. Función Par. Una función f (x) se dice par si f (x) = f (−x)∀x ∈ D(f ).

2. Función Impar. Una función se dice impar si f (x) es tal que −f (x) = f (−x)∀x ∈ D(f )

donde D(f ) es el dominio de la función.

2
 Las funciones par e impar cumplen propiedades útiles de las cuales se pueden obtener ventajas en el cálculo
de las integrales.

Teorema 1.4
si f es par, entonces:
Z L Z L
f (x)dx = 2 f (x)dx (1.12)
−L 0

si f es impar, entonces:
Z L
f (x)dx = 0 (1.13)
−L

Usando las propiedades integrales de función par e impar y el teorema sobre cómo obtener la serie de Fourier
de una función se obtienen los siguientes resultados.

Teorema 1.5
Si f es una función 2p-periódica y además tiene serie de Fourier.

1. Si f es par entonces

∞
X nπ
f (x) = a0 + an cos( x) (1.14)
p
n=1
Z p
1
a0 = f (x)dx (1.15)
p 0
Z p
2 nπ
an = f (x) cos( x)dx (1.16)
p 0 p

2 Si f es impar entonces

∞
X nπ
f (x) = bn sin( x) (1.17)
p
n=1
Z p
2 nπ
bn = f (x) sin( x) (1.18)
p 0 p

1.4. Extensiones Par e Impar de una función y series de medio rango

Se refiere a que una función cualquiera definida en un intervalo finito 0 < x < p puede extenderse de manera
par o impar, hacerse periódica y en consecuencia obtener su serie de Fourier.
Una misma función puede tener tanto una serie de senos (para extensión impar) o una serie de cosenos
(extensión par).

3
 Definición 1.3: Extensión para e impar de una función

Sea f suave por partes definida en un intervalo finito 0 < x < p.

1. Su extensión par se define como

fp (x) = f (x) ∀x : 0 < x < p (1.19)

fp (−x) = f (x) ∀ − p < x < 0 (1.20)
fp (x) = f (x + 2p) (1.21)

2 Su extensión impar se define como:

fi (x) = f (x) ∀x : 0 < x < p (1.22)

fi (x) = −f (−x) ∀x : −p < x < 0 (1.23)
fi (x) = fi (x + 2p) (1.24)

Según sea la extensión que se elija la función tendrá una u otra representación en serie como las obtenidas
previamente.

Teorema 1.6
Sea f una función suave por trozos y definida en un intervalo finito 0 < x < p. Entonces la función tiene
dos representaciones en serie de Fourier.

1. Si f se extiende de manera par, entonces

∞
X nπ
f (x) = a0 + an cos( x) (1.25)
p
n=1
Z p
1
a0 = f (x)dx (1.26)
p 0
2 p
Z
nπ
an = f (x) cos( x)dx (1.27)
p 0 p

2 Si f se extiende de manera impar entonces

∞
X nπ
f (x) = bn sin( x) (1.28)
p
n=1
Z p
2 nπ
bn = f (x) sin( x) (1.29)
p 0 p

Las series obtenidas de esta forma se conocen como expansiones de medio rango de f.

1.5. Serie compleja de Fourier

2. Modelos que obedecen la ecuación de onda

La ecuación se obtiene modelando una cuerda de longitud L sujeta y fija por sus extremos en x = 0 y x = L.
Se supone que la cuerda experimenta vibración como consecuencia de haber movido la cuerda de su punto de
equilibrio y liberándola. Las suposiciones que se hacen sobre la cuerda son las siguientes:

4
 Figura 2.1:

1. la cuerda tiene densidad de masa lineal constante

2. la cuerda no ofrece resistencia a doblarse, es decir, es perfectamente elástica.

3. el movimiento de la cuerda ocurre en el plano u − x de manera transversal.

∂u
4. la tensión horizontal es cada punto de la cuerda es constante y la pendiente de la cuerda ∂x en cada punto
es muy pequeña.

5. el peso de la cuerda se considera despreciable.

A continuación se muestran diferentes modelos que implican la ecuación de onda o una variación de [Link]
modelos son:
1. Vibración libre de una cuerda.

2. Vibración forzada de una cuerda.

3. Vibración amortiguada de una cuerda.

4. Cadena colgante

5. Vibración libre de una membrana rectangular

2.1. Vibración libre de una cuerda

Se parte del modelo fı́sico mostrado en figura 2.1 con las siguientes suposiciones adicionales:

1. No hay fuerzas externas actuando sobre la cuerda

2. La tensión de la cuerda es constante en cada punto de la cuerda.

3. La tensión de la cuerda es tangente en cada punto.

Se usa la suma de fuerza sobre el eje u, la segunda ley de Newton y como se asumen ángulos pequeños por lo
que tan θ ≈ sin θ y además tan θ = ∂u(x,t)
∂x . Bajo estas consideraciones el modelo matemático que se obtiene es el
siguiente:

∂2u 2
2∂ u
= c (2.1)
∂t2 ∂x2
Observación:Cabe recalcar que el modelo obtenido a partir de la vibración libre de una cuerda de longitud L y
fija en los extremos 0 y L es el modelo de referencia y que comúnmente identificamos como ecuación de onda

5
2.2. Vibración forzada
El modelo fı́sico del cuál se deriva la vibración forzada de una cuerda aparece en la figura 2.2. La diferencia
con el modelo de la cuerda libre es que se considera una fuerza externa vertical F (x, t)∆x donde F (x, t) es la
fuerza por unidad de longitud. Derivar matemáticamente el modelo es igual que en el caso previo pero agregando
en la suma de fuerza la fuerza externa como se muestra en la figura 2.2. El modelo matemático obtenido es el

Figura 2.2:

siguiente:
∂2u ∂ 2 u F (x, t)
2
= c2 2 + (2.2)
∂t ∂x ρ
dos casos especiales que se derivan de la ecuación es cuando la fuerza externa se debe al campo gravitacional y
cuando la fuerza externa se debe a la resistencia la cuál se considera proporcional a la velocidad instantanea. Estos
modelos son los siguientes:
∂2u 2
2∂ u
= c −g (2.3)
∂t2 ∂x2
∂2u 2
2∂ u ∂u
2
= c 2
−k (2.4)
∂t ∂x ∂t

2.3. Cuerda colgante

El modelo fı́sico aparece en la figura 2.3. Al deducir el modelo de la cuerda colgante se asume:
1. La cuerda está fija en L
2. La tensión no es constante, por lo que T (x + ∆x) ̸= T (x)
3. Se mantiene la suposición de ángulos pequeños.
Es importante señalar que en el modelo fı́sico las tensiones aparecen apuntando hacia la izquierda, sin embargo,
para obtener el modelo correctamente las tensiones deben apuntar en direcciones opuestas, recomendando que la
tensión superior apunte hacia la derecha. El proceso es como en los casos previos: suma de fuerzas sobre u, 2da
ley de Newton y simplificación considerando ángulos pequeños ası́ como la definición de derivada.
El modelo matemático de la cadena colgante es:
∂2u ∂2u
 
∂u
=g +x 2 (2.5)
∂t2 ∂x ∂x

6
 con

Figura 2.3:

2.4. Vibración libre en una membrana rectangular

3. Separación de Variables y condiciones de frontera

El método de separación de variables se utiliza cuando se cumplen las siguientes condiciones a la vez:

1. La EDP es lineal, homogénea y no tiene derivadas cruzadas. O se puede homogeneizar por algún cambio
de variables

2. Sus condiciones de frontera son homogéneas o se pueden volver homogéneas por algún cambio de variable.

3. La región sobre la cuál se busca obtener la solución coincide con los ejes del sistema coordenado en el cuál
se está trabajando para las fronteras se puedan separar.

Las condiciones de frontera se clasifican en:

1. Condiciones de Dirichlet. Especifican el valor de la solución en la frontera del dominio.

2. Condiciones de Neumann. Especifican el valor de la derivada normal en la frontera del dominio.

3. Condiciones de Robin. Cuando se especifica el valor de la combinación lineal de la función y su derivada

normal en la frontera del dominio.

Adicionalmente, cuando se aplica la separación de variables se busca que en algunas de las variables se
obtenga un problema de Sturm Liouville, que permitirá expresar la solución en término de las eigenfunciones y
obtener sus eigenvalores.

4. Solución de la ecuación de onda en 1 dimensión y en coordenadas

rectangulares
La solución de la ecuación de onda

utt = c2 uxx (4.1)

dependerá de las condiciones de frontera, iniciales ası́ como de la región en donde se busca solucionar la ecuación.
En el caso de las condiciones iniciales, aparece la posición inicial u(x, 0) y la velocidad inicial de la cuerda
ut (x, 0).

7
 Teorema 4.1
La solución de la ecuación

utt = c2 uxx

sujeta a:

1. BC: u(0, t) = u(L, t) = 0

2. IC: u(x, 0) = f (x) y ut (x, 0) = g(x)

Por separación de variables, está dada por:

∞
X nπc nπc nπ
u(x, t) = (An cos( t) + Bn sin( t)) sin( x) (4.2)
L L L
n=1

donde los coeficientes están dados por:

2 L
Z
nπ
An = f (x) sin( x)dx (4.3)
L 0 L
Z L
2 nπ
Bn = g(x) sin( x)dx (4.4)
nπc 0 L

Donde al deducir esta expresión se usa constantemente algunas de las siguientes relaciones de ortogonalidad.
Z +L
nπx mπx
sin( ) sin( )dx = Lδnm (4.5)
−L L L
Z L
nπx mπx L
sin( ) sin( )dx = δnm (4.6)
0 L L 2

Esta solución corresponde al caso de una cuerda de longitud L sujeta por ambos extremos. Si estas condiciones
de frontera cambian es necesario hacer ajustes al problema, es decir homogeneizar las condiciones de frontera,
como se hace en el siguiente problema cuando las condiciones de frontera no son homogéneas, sino son iguales a
constantes.

8
 Teorema 4.2
La solución de la ecuación de onda

utt = c2 uxx (4.7)

sujeta a las siguientes condiciones frontera no homogéneas

1. B.C: u(0, t) = A y u(L, t) = B ambas constantes diferente de cero.

2. I.C: u(x, 0) = f (x) y ut (x, 0) = g(x)

se soluciona con el cambio:

w(x, t) = u(x, t) − u1 (x) (4.8)

donde u1 (x) = (B − A)/Lx + A. Con lo que se obtiene el problema equivalente:

wtt = c2 wxx (4.9)

w(0, t) = w(L, t) = 0 (4.10)
w(x, 0) = f (x) − u1 (x) = F (x), wt (x, 0) = g(x) (4.11)

cuya solución corresponde al problema de la cuerda con condiciones homogéneas. Entonces la solución al
problema con condiciones de forntera no homogéneas está dada por:

u(x, t) = w(x, t) + u1 (x) (4.12)

∞
X nπc nπc nπ B−A
u(x, t) = (An cos( t) + Bn sin( t)) sin( x) + x+A (4.13)
L L L L
n=1

y cuyos coeficientes están dados por:

Z L
2 nπ
An = (f (x) − u1 (x)) sin( x)dx (4.14)
L 0 L
Z L
2 nπ
Bn = g(x) sin( x)dx (4.15)
nπc 0 L

4.1. Nodos y modos normales

1. Nodos: Puntos cuyo desplazamiento es nulo durante la vibración de la cuerda. Es decir, los nodos son los
puntos tales que u(x, t) = 0, Por ejemplo, para el caso de la cuerda con condiciones homogéneas son los
puntos tales que sin( nπ
L x) = 0.

2. Modos normales: De las soluciones de la cuerda se llaman modos normales a los dados por un (x, t). De
modo que la vibración de la cuerda o el sistema es la superposición de los modos normales de vibración.

3. Modo Fundamental: Es el primer modo normal de vibración, es decir, n = 1. El resto de modos normales
se conocen como sobretonos.
p
4. Intensidad: La intensidad del sonido producido por un modo normal está dado por A2n + Bn2 para la
cuerda con BC homogéneas.
nπc
5. Frecuencia natural o ciruclar está dada por ωn = L .

6. Relación de la velocidad con la tensión. En el caso de la cuerda la velocidad

p de disipación se relaciona
con la tensión y con la densidad de la cuerda por la siguiente expresión c = τ /ρ.

9
4.2. Método de D’Alambert

5. Ecuación de calor en 1 dimensión y en coordenadas rectangulares

Se conoce como ecuación de calor o ecuación de difusión a la EDP de la forma:
∂u ∂2u
= c2 2 (5.1)
∂t ∂x
La solución de la ecuación de calor depende, como en el caso de la ecuación de onda, de las condiciones de
frontera, inicial y de la región donde se soluciona la EDP. Algunas soluciones clásicas son las que se muestran a
continuación.

Teorema 5.1: Ecuación de calor con BCH

La solución del siguiente problema

ut = c2 uxx (5.2)
u(0, t) = u(L, t) = 0 (B.C.H) (5.3)
u(x, 0) = f (x) (I.C) (5.4)

está dada por:

∞
X nπ nπc
u(x, t) = An sin( x) exp(− t) (5.5)
L L
n=1
Z L
2 nπ
An = f (x) sin( x) (5.6)
L 0 L

La ecuación de calor tiene una solución de estado estable la cuál es independiente del tiempo, es decir,
es la solución de uxx = 0 por tanto dicha solución es:

Teorema 5.2: Solución de estado estable

T2 − T1
u(x) = x + T1 (5.7)
L

Pero a partir de la solución de estado estable se puede formular una para condiciones de frontera no
homogéneas que se establece en el siguiente resultado.

10
 Teorema 5.3: Ecuación de calor con B.C.N.H
La solución del siguiente problema

ut = c2 uxx (5.8)
u(0, t) = T1 y u(L, t) = T2 (5.9)
u(x, 0) = f (x) (5.10)

Donde T1 y T2 son constantes. Entonces, el problema se transforma a un problema con BCH como en el
caso original por medio del siguiente cambio de variable

w(x, t) = u(x, t) − u1 (x) (5.11)

T2 − T1
u1 (x) = x + T1 (5.12)
L
de modo que:

w(0, t) = w(L, t) = 0 B.C.H (5.13)

w(x, 0) = F (x) = f (x) − u1 (x) I.C (5.14)

cuya solución es de la forma 5.5 con la diferencia de que el coeficiente está dado por:
Z L
2 nπ
An = (f (x) − u1 (x)) sin( x)dx (5.15)
L 0 L

Por tanto la solución del problema de la ecuación de calor sujeta a condiciones no homogéneas está dada
por:

u(x, t) = w( x, t) + u1 (x) (5.16)

donde w(x, t) es la solución del problema con condiciones de frontera homogéneas obtenido por el cambio
de variables.

1. Solución de estado estable: Es una función que depende solamente de x que se obtiene cuando t → ∞.

6. Ecuación de onda en 2 dimensiones y coordenadas rectangulares.

En dos dimensiones la ecuación de onda se define como:

Definición 6.1

∂2u ∂2u ∂2u

 
= c2 + 2 (6.1)
∂t2 ∂x2 ∂y

El problema que se soluciona es el de la membrana rectangular sometida a una vibración.

11
 Teorema 6.1
La solución del problema:

utt = c2 (uxx + uyy ) (6.2)

BC (6.3)
u(0, y, t) = u(a, y, t) = 0 (6.4)
u(x, 0, t) = u(x, b, t) = 0 (6.5)
IC (6.6)
u(x, y, 0) = f (x, y) (6.7)
ut (x, y, 0) = g(x, y) (6.8)

Tiene como solución:

∞ X
∞
X mπ nπ
u(x, y, t) = (Amn cos(λmn t) + Bmn sin(λmn t)) sin( x) sin( y)dxdy (6.9)
a b
n=1 m=1

con coeficientes
Z bZ a
4 mπ nπ
Amn = f (x, y) sin( x) sin( y)dxdy (6.10)
ab 0 0 a b
Z bZ a
4 mπ nπ
Bmn = g(x, y) sin( x) sin( y)dxdy (6.11)
abλmn 0 0 a b

donde
r 
m 2  n 2
λmn = cπ + (6.12)
a b

6.1. Lineas nodales

1. Linea nodal: Son la lineas formadas por el conjunto de puntos que permanecen fijos durante la vibración
de la membrana.

2. Las lineas nodales son los puntos (x, y) : un (x, y, t) = 0 ∀t > 0, 0 ≤ x ≤ a; 0 ≤ y ≤ b, es decir, son los puntos
que satisfacen la ecuación

mπ nπ
sin( x) sin( y) = 0 (6.13)
a b

12
7. Ecuación de calor en 2 dimensiones y en coordenadas rectangulares.

8. Ecuación de Laplace: caso particular y caso general.

9. Eigenfunciones de la ecuación de Laplace, Ecuación de Helmholtz y solución

por eigenfucniones para [Link] homogéneas.

10. Método por transformada de Fourier y principio de Duhamel para ecuación

de onda y ecuación de calor.

11. Método por transformada de Laplace.

12. Ecuación de Laplace en sistemas de coordenadas curvilineos

13. Ecuación de onda en sistemas curvilineos ortogonales

14. Ecuación de calor en sistemas curvilineos ortogonales

13