INSTITUTO TECNOLÓGICO

SUPERIOR DE VILLA LA VENTA

CARRERA: ING INDUSTRIAS ALIMENTARIAS

ECUACIONES DIFERENCIALES

DOCENTE:MARIA ISABEL SOBERANO MAYO

ALUMNO: DANIEL HERNANDEZ GOMEZ

INVESTIGACIÓN
Transformada de Laplace

3.1 Teoría preliminar

3.1.1 Definición de la transformada de Laplace

Las transformadas de Laplace fueron formuladas para transformar una

ecuación diferencial que contiene las diferenciales de una función indefinida,
a partir de una ecuación t-espaciada hacia una ecuación s-espaciada que
puede ser resuelta con mucha facilidad. También pueden ser utilizadas para
resolver un sistema de ecuaciones diferenciales. Están dirigidas a una amplia
gama de problemas de valor inicial.

La transformada de Laplace se denomina a veces transformada operacional;

esto esporque transforma las operaciones de integración en simples
operaciones algebraicas que son mucho más convenientes de resolver.
Después de la cual la aplicación de la técnica de la transformada inversa
produce la solución exacta para la ecuación diferencial dada.

3.1.2 Condiciones suficientes de existencia para la transformada de Laplace

Antes de establecer las condiciones para la existencia de la transformada de

Laplace es esencial entender dos conceptos fundamentales que constituyen
la base de la transformada de Laplace. Estos son:

1. Función continua a trozos: Se dice que una función es a trozoso

seccionalmente continua en un intervalo finito a <= t <= b si el intervalo se
puede subdividir en un número finito de subintervalos, en cada uno de los
cuales f(t) es continua y tiene límites izquierdos, así como límites derechos.

Considera una función f(t) que es continua a trozos en [a, b], pero presenta
discontinuidades en algunos puntos.

3.2 Transformada directa

Uno de los algoritmos aritméticos más ampliamente utilizados es la

transformada rápida de Fourier, un medio eficaz de ejecutar un cálculo
matemático básico y de frecuente empleo. La transformada rápida de Fourier
es de importancia fundamental en el análisis matemático y ha sido objeto de
numerosos estudios. La aparición de un algoritmo eficaz para esta operación
fue una piedra angular en la historia de la informática.

Las aplicaciones de la transformada rápida de Fourier son múltiples. Es la

base de muchas operaciones fundamentales del procesamiento de señales,
donde tiene amplia utilización. Además, proporciona un medio oportuno para
mejorar el rendimiento de los algoritmos para un conjunto de problemas
aritméticos comunes.

3.3 Transformada inversa

Si la transformada de Laplace de una función f(t) es F(s), esto es,

Entonces f(t) se denomina la transformada inversa de Laplace de F(s) y se

escribe como,

Aquí L-1 es llamado el operador de la transformada inversa de Laplace.

Aunque existe una fórmula de inversión compleja la cual proporciona un
medio directo para determinar la transformada inversa de Laplace de una
función, esta implica un conocimiento amplio de la integración compleja. Sin
embargo, no es posible encontrar una transformada inversa de Laplace de
todas las funciones por este camino.

Un punto interesante a destacar aquí es que la transformada inversa de

Laplace de una función puede no ser única. Por ejemplo, sabemos que la
transformada de Laplace de f(t) = 1 es 1/ s. Sin embargo, existe una función
más para la que la transformada de Laplace resulta en el mismo término,
esta es,
f(t) = 1 si 0 < t < 3

−8 si t = 3

1 si t > 3

Por lo tanto, se puede concluir que la transformada inversa de Laplace de

1/s resulta para dos funciones como se describe anteriormente. Sin embargo,
entre las dos, sólo f(t) = 1 es continua, por lo tanto, f(t) = 1 es la única función
que tiene la transformada de Laplace como 1/ s.

También es posible escribir una transformada inversa de Laplace en la forma

de integración, la cual es llamada integral de Fourier Millen o integral
Bromwich o fórmula inversa de Millen. Se trata de una integral de línea y es
denotada como,

Propiedades

Propiedades de la Transformada Inversa de Laplace

Similar a la transformada de Laplace, la transformada inversa de Laplace
también tiene sus propias propiedades. Algunas de las propiedades más
importantes se discuten a continuación:

1. Propiedad de Linealidad: Si L{f(t)} = F(s) y L{g(t)} = G(s), entonces para

dos constantes cualesquiera c1 y c2 tenemos,

L-1{c1 f(t) + c2 g(t)} = c1 L-1{f(t)} + c2 L-1{g(t)}

= c1f(t) + c2 g(t)

Esto puede probarse como,

Por la propiedad de linealidad de la transformada de Laplace conocemos

que,

L{c1 f(t) + c2 g(t)} = c1 F(s) +c2 G(s)

Ahora tomando la transformada inversa de Laplace en ambos lados

obtenemos, c1f(t) + c2 g(t) = L-1{c1 F(s) +c2 G(s)} f(t) = L-1{F(s)} y, g(t) =

L-1{G(s)}

Transformada de Laplace de funciones definidas por tramos

Una función continua a trozos es aquella que es continua y está dividida en

pedazos. Antes de sumergirnos en el concepto primero debemos entender
lo que es una función continua y una función a trozos. Una función continua
es aquella que no se divide en su gráfico, es decir, se define continuamente
a lo largo de todo el dominio de la función. Un ejemplo de tal función sería =
y2, esta es la ecuación de una parábola.

Como podemos ver en el gráfico anterior, esta no se rompe.

Y, una función a trozos es aquella que no está definida continuamente, esto

es, la gráfica de la función está dividida. Unejemplo de
estoseríaunafunciónescalonada.
Como podemos observar en la imagen de arriba, la función no es continua.

Una función continua a trozoses aquella que combina los dos tipos de
funciones. Esta es una combinación muy interesante,porque¿Cómo puede
una función ser continua y no continua al mismo tiempo?.El siguiente gráfico
haría el concepto claro.
Como se puede observar en el ejemplo anterior, el gráficono es continuo,
pero el valor para la función está definido en cada punto, esto se debe a que
en varios puntos la gráfica tiene dos valores en lugar de uno, lo cual se divide
el gráfico aunque lo mantiene continuo.

Claramente, f(t) es continuo en los intervalos (a, A), (A, B), (C, y), (y, D) y (D,
b). También los límites derechose izquierdos de A son,

f(A + t) = f(A + 0) = f(A)

f(A - t) = f(A - 0) = f(A)

Aquí el valor de t siempre es positivo.

Por lo tanto, podemos concluir que una función continua a trozos se compone
de números finitos de secciones de continuidad para cada sub-intervalo finito
[0, t] en la gráfica de la función dada. Y el límite de la función dada es finitoa
medida que la función se aproxima a cada uno de los puntos continuos.

La transformada de Laplace de tal función está dada como,

3.4 Función escalón unitario

La función escalón unitario o función escalón unitario de Heavisideu(t - a)

está definida como,

Aquí el valor de a es siempre mayor o igual que cero.

El gráfico de una función escalón unitario es parecido al siguiente,

Entonces, al observar el gráfico de la función, este se puede comparar con

un interruptor que se encuentra cerca de un tiempo en particular, que abre
por un tiempo y luego vuelve a cerrar.

Existen ciertas propiedades de una función escalón unitario relacionadas

con la transformada de Laplace. Estasse analizan a continuación:
3.5 Propiedades de la transformada de Laplace
(linealidad, teoremas de traslación)

Propiedades de la Transformada Inversa de Laplace

Similar a la transformada de Laplace, la transformada inversa de Laplace

también tiene sus propias propiedades. Algunas de las propiedades más
importantes se discuten a continuación:

1. Propiedad de Linealidad: Si L{f(t)} = F(s) y L{g(t)} = G(s), entonces para

dos constantes cualesquiera c1 y c2 tenemos,
L-1{c1 f(t) + c2 g(t)} = c1 L-1{f(t)} + c2 L-1{g(t)}

= c1f(t) + c2 g(t)

Esto puede probarse como,

Por la propiedad de linealidad de la transformada de Laplace conocemos

que,

L{c1 f(t) + c2 g(t)} = c1 F(s) +c2 G(s)

Ahora tomando la transformada inversa de Laplace en ambos lados

obtenemos, c1f(t) + c2 g(t) = L-1{c1 F(s) +c2 G(s)} f(t) = L-1{F(s)} y, g(t) =

L-1{G(s)}
3.6 Transformada de funciones multiplicadas por tn, y divididas entre t

3.7 Transformada de derivadas (teorema)

Al igual que en una función ordinaria, la transformada de Laplace también

puede aplicarse al diferencial de una función. En tal situación, colocamosen
la fórmula el diferencial de la función en el lugar de la función real para derivar
la transformada de Laplace, que es,

Sin embargo, mientras nos ocupamos de los diferenciales de la función,

necesitamos modificar el límite inferior de integración y colocar un valor
mayor que cero en el lugar de cero, como el límite inferior de integración.
Esto se hace principalmente porque el cero no manipula la solución obtenida
a partir de la integración, y de esta forma nos limitamos a la función clásica.

Existen, sin embargo, ciertas condiciones para que esto sea verdadero. La
función real debe ser definida para la variable tiempo ty el diferencial debe
existir para todos los valores mayores que cero. Asimismo, la función debe
ser definida de forma continua en el intervalo [0, ). De igual manera, el
diferencial de esta función debe ser una función continua a trozos para el
mismo intervalo, este es, [0, ).
Y, por último, tanto la función real, así como el diferencial de la función real
deben ser de orden exponencial cuando el valor de t tiende al infinito. Esto
significa que deben existir dos números reales positivos M y , y un número T
tal que,

Aquí el valor de t debe ser siempre mayor o igual que T.

3.8 Transformada de integrales (teorema)

Hasta ahora hemos estudiado la forma de determinar la transformada de

Laplace de una función dada. Pero, como sabemos, existen varias
operaciones que pueden realizarse en una determinada función. Una de las
principales operaciones entre ellas es la integración. Como sabemos, la
integración de una función nos da otra función. Por lo tanto, es esencial saber
si la técnica de la transformada de Laplacepuede aplicarse a la integral de
una función real.

La respuesta a la pregunta anterior es afirmativa, hasta cierto punto. La

cláusula de “hasta cierto punto”, se añade aquí porque esto no es cierto para
todas las integrales. Existen ciertos pre-requisitos que deben ser verdaderos
para obtener la transformada de Laplace de la integral de la función real.
La función real debe estar definida para la variable de tiempot, también la
función debe ser definida de forma continua en el intervalo [0, ). Asimismo,
el integral de esta función debe ser una función continua a trozos para el
mismo intervalo, esto es, [0, ).

Y, por último, ambas, la función real como el diferencial de la función real

deben ser de orden exponencial cuando el valor de t tiende al infinito. Esto
significa que deben existir dos números reales positivos M y , un número T
tal que,

Aquí el valor de t debe ser siempre mayor o igual que T.

Asumiendo que las condiciones anteriores son válidas para la función de

entrada y L{f(t)} = F(s), entonces la transformada de Laplace de la integral
de la función real puede darse como,

Teorema de la convolución

Si L-1{F(s)} = f(t) y L-1{G(s)} = g(t)entonces,

O esta puede ser reescrita como,

Aquí f * g se llama la convolución de F y G. esto es llamado el teorema de

convolución de la transformada de Laplace. Es una propiedad importante de
la transformada de Laplace. El teorema anterior indica que puede probarse
como,

A partir de la definición de la transformada de Laplace sabemos que,

L{ f(u) g(t – u) du} = e-st { f(u) g(t – u) du} dt

= e-stf(u) g(t – u) du dt

Aquí la región de integración es la parte sombreada de la siguiente figura,

Ahora, cambiando el orden de integración, tenemos la siguiente parte
sombreada como la región de integración.

Entonces obtenemos,

L{ f(u) g(t – u) du} = e-stf(u) g(t – u) du dt

= f(u) { e-st g(t – u) dt} du

= fijando t – u = v obtenemos, dt = dv

= f(u) { e-s(u + v) g(v) dv} du

= { e-suf(u) du}. e-sv g(v) dv}

O, L{ f(u) g(t – u) du} = F(s) G(s)

Invirtiendo ambos lados de la ecuación obtenemos,

L-1{F(s) G(s)} = f(u) g(t – u) du

L-1{F(s) G(s)} = f * g

Existe una cantidad amplia de problemas que pueden resolverse con la

ayuda del teorema de convolución. Uno de estos problemas se da aquí para
hacer el concepto más claro.

Usa el teorema de convolución para determinar L-1{1/ [s2 (s + 1)2]}

Sea F(s) = 1/ s2

Y G(s) = 1/ (s + 1)2entonces,

f(t) = L-1{F(s)} = L-1(1/ s2) = t

g(t) = L-1{G(s)} = L-1(1/ (s +

1)2)

= e-t L-1(1/ s2)

= t e-t

Por lo tanto, con la ayuda del teorema de convoluciónpodemos escribir,

L-1{1/ [s2 (s + 1)2]} = f(u) g(t – u) du

= u (t – u) e-(t – u) du

= e-t [t u eudu - u2eu du]

= e-t [t (u eu – eu) - (u2eu) + 2 u eu du]

= e-t [t {(t – 1) et + 1} – t2 et + 2(ueu – eu) ]

= e-t [t (t – 1) et + t + t2 et + 2(tet – et + 1)]

= e-t [t2 et - tet+ t - t2 et + 2tet - 2et + 2]

= e-t [tet+ t- 2et + 2]

= t + tet– 2 + 2et
El teorema de convolución tiene amplias aplicaciones en la práctica. También
se utiliza en la teoría de circuitos para calcular la respuesta al impulso de un
circuito concreto.

Aquí x(t) es la entrada del sistema, y(t) es la salida del sistema y h(t) es la
respuesta al impulso del sistema.

Por consiguiente, la salida del sistema calculado con la ayuda de la operación

de convolución está dada por,

y(t) = x(t) * h(t)

Aquí la función anterior está en el dominio de t, por lo tanto, los cálculos

pueden ser algo crípticos,por esto, con el propósito de conveniencia en los
cálculos, esta puede ser transformada en el dominio s. La operación de
convolución en el dominio s se convierte en la operación de multiplicación.

L{a(t) * b(t)} = A(s) B(s)

3.10Transformada de Laplace de una función periódica

Se dice que una función f(t) es una función periódica de período a> 0 si,
Esto significa que la gráfica de tal función a repetirá su forma para cada
intervalo (na, (n + 1)a). Un ejemplo de tal función es el seno ( ),el cual es una
función periódica del período 2 .

El valor de la función debe convertirse en cero en la porción negativa de la

recta numérica real.

Si f(t) es una función periódica con período a entonces,

Esto puede reorganizarse como,

En términos simples, podemos decir que para la función periódica f(t) con
período a, la transformada de Laplace es equivalente a la transformada de
Laplace en un período único de esafunción dividida por el término (1 - e-as).

También existe una prueba del teorema indicado arriba. Dado que f(t) es una

función periódica con período a, f(t + a) = f(t), f(t + 2a) = f(t) y así

sucesivamente.
Ahora, L{f(t)} = e-st f(t) dt

= e-stf(t) dt + e-st f(t) dt + e-st f(t) dt + …

Estableciendo t = (u + a) en la segunda integral, t = (u + 2a) en el tercer

integral etc. Entonces obtenemos,

= e-stf(t) dt + e-s(u + a) f(u + a) du + e-s(u + 2a) f(u + 2a) du + …

= e-suf(u) du + e-su f(u + a) du + e-2sa e-su f(u) du + …

= (1 + e-as + e-2as+ …) e-su f(u) du

= (1 - e-as)−1 e-stf(t) dt [Dado que 1 + x + x2 + … = (1 – x)−1]

L{f(t)} = [1/ (1 - e-as)] e-st f(t) dt

Por consiguiente, podemos ver que para llevar a cabo la operación de

transformación de Laplace a una función periódica necesitamos romper esa
función en los sub-intervalos, lo cual depende de los intervalos para los
cuales la función dada estádefinida. La sumatoria de todas las integrales
produce la transformada de Laplace para esa función.

Sin embargo, es posible aplicar directamente el

teorema discutido anteriormentecon propósitos de conveniencia para la
solución de problemas. Se provee un ejemplo ilustrando el uso de este teorema
para hacer más claro los conceptos.
Muestra que si f(t + a) = -f(t), entonces,

L{f(t)} = [1/ (1 + e-as)] e-st f(t) dt

Dado quef(t + 2a) = f[(t + a) + a)] = -f(t + a)

= -[-f(t)] = f(t)

La función f(t) dada es una función periódica con período 2a.

Ahora, utilizando el teorema para la transformada de Laplace de funciones

periódicas con el período areemplazado por 2a tenemos que,

L{f(t)} = [1/ (1 - e-2as)] e-st f(t) dt

= [1/ (1 - e-2as)] [ e-stf(t) dt + e-st f(t) dt]

Estableciendo t = (u + a) en la segunda integral, tenemos que

L{f(t)} = [1/ (1 - e-2as)] [ e-st f(t) dt + e-s(u + a) f(u + a) du]

= [1/ (1 - e-2as)] [ e-stf(t) dt – e-as e-su f(u) du]

= [(1 - e-as)/ (1 - e-2as)] e-stf(t) dt

L{f(t)} = [1/ (1 + e-as)] e-st f(t) dt

Función delta Dirac.

Las funciones deltade Dirac son las funciones que ejercen una enorme
cantidad de fuerza sobre un objeto, por una gran cantidad de tiempo. Aunque
a veces una función escalonada unitaria es comparada con una función
fuerza, la comparación no es muy adecuada dado que la cantidad de fuerza
ejercida por ellas es muy limitada. Una función delta de Dirac es una
diferencial de la función escalón unitario. Esta puede entenderse como
secuencias delta de funciones de fuerza generalizadas.

Esto implica que la función delta de Dirac no es una función real sino que es
una distribución que se extiende por un intervalo definido para la función
dada. También es llamada una función singular. Como tal, no existe una
definición formal de esta función. Pero puede ser definida mediante utilizar
la propiedad de la propia función, la cual es,

En términos simples, podemos decir que una función delta de Dirac es

aquella cuya salida se calcula a cero para cada valor del argumento de
entrada, excepto cuando el valor del argumento de la función en sí es igual
a cero. Aquí el argumento de la función es un parámetro valorado real. La
integral de la función en el rango de parámetros (- , ) es uno. A la luz de la
afirmación anterior, podemos concluir que esta es una función real desde el
punto de vista matemático, ya que para cualquier función real cuyo valor es
constante, excepto en un punto, el valor de la integral debe calcularse a cero,
el cual no es este caso.

La gráfica de la función se vería algo así como:

Debido a esta propiedad de la función, es ampliamente utilizada para

modelar el sistema que experimenta fuerzas extremasrepentinas. Una

propiedad muy importante de esta función es,

En este caso, sabemos que la función (t) toma el valor de cero para todos
los valores de t, excepto en t = 0. Esto implica que el valor de la función f(t)
también se vuelve insignificante, excepto cuando el argumento tde la función
se convierte en cero. En tal situación, tenemos el valor del integrando f(0) (t),
f(0)que puede tomarse fuera dado que se convierte en una constante,
haciéndolo de esta manera obtenemos el lado derecho de la ecuación.
Por lo tanto, podemos pensar en (t) dt como el operador funcional que saca
el valor de la función cuando el argumento de la función es igual a cero.

Otra forma popular de definir una función delta de Dirac es una medida, ya
que la función delta de Dirac tiene como su argumento un subconjunto de los
números realesS, es decir, S R. Aquí, el valor de S es cero cuando la salida
de la función es uno o infinita y en otro caso, el subconjunto S puede tomar
elementosinfinitos.

Veamos ahora un ejemplo de la función delta de Dirac.

Resuelve y’ + 2y’ – 15y = 6 (t – 9) dados y(0) = −5 ey’(0) = 7

Aplicando la transformada de Laplace para la función dada obtenemos,

s2Y(s) – sy(0) – y’(0) + 2(sY(s) – y(0)) – 15Y(s) = 6e-9s

= (s2 + 2s – 15)Y(s) + 5s + 3 = 6e-9s

Resolviendo la ecuación anterior obtenemos,

Y(s) = 6e-9s F(s) – G(s)

Ahora haciendo uso de las fracciones parciales para obtener la transformada

inversa de Laplace como,

F(s) = 1/ [(s + 5) (s – 3)]

= [(1/ 8)/ (s – 3)] – [(1/ 8)/ (s + 5)]

f(t) = (1/8) e3t - (1/8)e-5t

G(s) = [5s + 3/(s + 5) (s – 3)]

= [(9/4)/ (s – 3)] + [(11/4)/ (s +

5)] f(t) = (9/4) e3t - (11/4)e-5t

Por lo tanto, la solución es

Y(s) = 6e-9s F(s) – G(s)

3.11 Transformada de Laplace de la función delta Dirac

Una función delta de Dirac es una función especial cuyo valor es cero en todos
los puntos, excepto en un punto, este es cuando el argumento de la función es
igual a cero. Esto se denota como,
Cambiemos la función delta de Dirac por una constante, digamos c. Entonces
ahora la definición de esta función delta de Dirac desplazada es,

(t – c) = 0, t <> c

= , t=c

Esto es sólo una pseudodefinición de la función. Ahora bien, si derivamos el

área de la función para los límites de integración (- , ), y resulta ser uno, esto

es, (t – c) dt = 1

Se trata de una derivación importante y esta también nos da la noción de

pseudoinfinidad,como en la definición función delta de Diracdesplazada.
Aquí, la palabra pseudo se utiliza ya que pueden existir diferentes medidas
del infinito mediante tomar un producto del infinito con un número entero.
Para entenderlo, integremos el producto de la función delta de Dirac y un
entero, digamos dos para los mismos límites de integración, esto es, (- , ).

2 (t – c) dt

Uno podría suponer que la salida de la integración debería ser igual a dos,
ya que,
= 2 (t – c) dt

= (2) (1)

=2

Es decir, si la función delta de Dirac se multiplica por dos, el infinito sería dos
veces más grande que antes.

Ahora, multiplicando la función delta de Dirac desplazada por alguna otra

función,digamos f(t) y tomando la transformada de Laplace de esta, es decir,

L{ (t – c) f(t)}

En el caso de que uno desee determinar únicamente la transformada de

Laplace de la función delta de Dirac desplazada, asumimos que el valor de

f(t) es uno. Esto es, tenemos,

e-st f(t) (t – c) dt

Como sabemos, el proceso de integración nos da el área de la función que está

siendo integrada. Por lo tanto, primero dibujemos el área de las dos funciones
para averiguar qué área estamos determinando realmente. Mientras lo
hacemos,asume que f(t) es arbitrario. Por tanto, tenemos el gráfico de la función
como,
Aquí se dibuja una línea recta donde t = c ya que el valor de la función delta
de Dirac es siempre cero, excepto en t = c. Por lo tanto, el espacio común de
las dos curvas, cuyo valor será determinado por la operación de integración
viene a ser un solo punto, el cual es el punto de intersección de las dos
curvas, y el valor de la primera función en ese punto en seráe-sc f©. Este es
sólo un punto, el cual tiene un valor constante.

En consecuencia, tenemos un término constante dentro de la operación de

integración que se puede mover fuera y, por lo tanto, quedamos con, e-sc f©

(t – c) dt

Como sabemos, el valor de la integral (t – c) dtes uno, por esto, la

transformada de Laplace de la función delta de Dirac desplazada es e-sc f©.
Esto nos da la transformada de Laplace de la función delta de Dirac, donde
el valor de c = 0 y f(t) = 1, como,

L{ (t)} = e0 (1) = 1

L{ (t - c)} = e-cs (1) = e-cs

Solución de ecuaciones

La transformada de Laplace es especialmente útil para obtener la solución

de las ecuaciones diferenciales ordinarias no homogéneas con coeficientes
constantes, donde todas las condiciones de contorno se dan para la función
desconocida y sus diferencias en un solo punto. El procedimiento de trabajo
de la misma es la siguiente:

Sea el problema de valores iniciales dado como,

(i)

Donde y(0) = k0 e y’(0) = k1. Además a1, a2, k0, k1son todos constantes y
f(t) es función de t solamente.

1. Aplicando la transformada de Laplace en ambos lados de la ecuación

(i) tomando en cuenta que
L(d2y/ dt2) = s2Y(s) – sy(0) – y’(0) y,

L(dy/ dt) = sY(s) – y(0)

Donde, Y(s) = L{y(t)} e F(s) = L{f(t)}

Entonces, la ecuación (i) produce,

[s2Y(s) – sy(0) – y’(0)] + a1[sY(s) – y(0)] + a2Y = F(s)

Ahora, haciendo uso de las condiciones iniciales tenemos que,

(s2 + a1s + a2) Y(s) = F(s) + sk0 + k1 + a1k0 (ii)

2. Resuelve la ecuación (ii) y luego expresa el lado derecho como una

sumatoria de fracciones parciales.

Aplica la transformada inversa de Laplace a Y(s), obtenida en el paso

anterior. Esto dará la solución de la ecuación dada (i) con condiciones

iniciales, y(t) = L-1{Y(s)}

El procedimiento anterior también puede aplicarse a las ecuaciones

diferenciales ordinariasde orden superior.

Ahora demos un vistazo a un ejemplo ilustrativo en la categoríaanterior.

Resuelve y’’ + 2y’ + 5y = e-t sin (t) dados y(0) = 0 e y’(0) = 1.

Al tomar la transformada de Laplace de ambos lados conseguimos,

L{y’’} + L{2y’} + L{5y} = L{e-t sin (t)}

O, [s2Y(s) – sy(0) – y’(0)] + 2[sY(s) – y(0)] + 5Y = [1/ ((s + 1)2 + 1)]

Usando y(0) = 0 e y’(0) = 1 tenemos,

(s2 + 2s + 5) Y(s) – 1 = [1/ (s2 + 2s + 2)]

Y(s) = [1/ (s2 + 2s + 2) (s2 + 2s + 2)] + [1/ (s2 + 2s + 2)]

= (1/ 3) {[1/ (s2 + 2s + 2)] – [1/ (s2 + 2s + 2)]} + [1/ (s2 + 2s + 2)]

= (1/ 3) [1/ (s2 + 2s + 2)] + (2/ 3) [1/ (s2 + 2s + 2)]

= (1/ 3) [1/ ((s + 1)2 + 1)] + (2/ 3) [1/ ((s + 1)2 + 4)]

Invirtiendo ambos ladostenemos, y(t) = (1/ 3) L-1[1/ ((s

+ 1)2 + 1)] + (2/ 3) L-1[1/ ((s + 1)2 + 4)] = (1/ 3) e-t L-

1[1/ (s2 + 1)] + (2/ 3) e-t L-1[1/ (s2 + 4)] [Utilizando el

primer teorema de desplazamiento]

= (1/ 3) e-t sin (t) + (1/ 3) e-t sin (2t)

= (1/ 3) e-t (sin (t) + sin (2t))

La transformada de Laplace también puede utilizarse para resolver un
sistema de n ecuaciones diferenciales ordinarias simultáneas en n variables
dependientes, las cuales son funciones de la variable independiente t.

Considera un sistema de dos ecuaciones diferenciales ordinarias

simultáneas en dos variablesdependientes x e y, las cuales son funciones de
t.

(a1D2 + a2D + a3)x + (a4D2 + a5D + a6)y = f1(t)

(b1D2 + b2D + b3)x + (b4D2 + b5D + b6)y = f2(t)

Aquí D = (d/ dt) y las condiciones iniciales sonx(0) = c1, x’(0) = c2, y(0) = c3,
y’(0) = c4. También ai(i = 1, ), bi(i = 1, ), ci(i = 1, ) son constantes. El
procedimiento de trabajo de la misma es la siguiente:

1. Aplicando la transformada de Laplace a ambos lados de las dos

ecuaciones diferenciales ordinarias dadas, obtenemos

{a1[s2X(s) – sx(0) – x’(0)] + a2[sX – x(0)] + a3X} + {a4[s2Y(s) – sy(0) – y’(0)]

+ a5[sY – y(0)] + a6Y} = F1(s)

{b1[s2X(s) – sx(0) – x’(0)] + b2[sX – x(0)] + b3X} + {b4[s2Y(s) – sy(0) – y’(0)]

+ b5[sY – y(0)] + b6Y} = F2(s)

O, (a1s2 + a2s + a3)X + (a4s2 + a5s + a6)Y = F1(s) + s(a1c1 + a4c3) + (a1c2
+ a2c1 + a4c4 + a5c3)
(b1s2 + b2s + b3)X + (b4s2 + b5s + b6)Y = F2(s) + s(b1c1 + b4c3) + (b1c2 +
b2c1 + b4c4 + b5c3)

Donde, X(s) = L-1{x(t)}

Y(s) = L-1{y(t)}
Unidad IV: Sistemas de ecuaciones diferenciales lineales

4.1 Teoría preliminar

4.1.1 Sistemas de EDL

Los problemas de la vida real pueden representarse de mejor manera con

la ayuda de múltiples variables. Por ejemplo, piensaen el conteo de la
población representado con la ayuda de una sola variable. Pero, esta
depende del conteo de la población de depredadores así como también
de las condiciones climáticas y la disponibilidad de alimentos. Todas estas
condiciones en sí mismas forman una ecuación diferente definida en una
variable separada.

Por lo tanto, para estudiar las relaciones complejas, requerimos de varias

ecuaciones diferentes para definir diferentes variables. Tal sistema es el
sistema de ecuaciones diferenciales. Un sistema de ecuaciones
diferenciales lineales se puede denotar como,
Aquí xi (t) es una variable en términos de tiempo y el valor de i = 1, 2, 3,
…, n.
También A es una matriz que contiene todos los términos constantes,
como [ai,j].

Dado que los coeficientes de la matriz constante A no están definidos

explícitamente en términos de tiempo, por lo tanto, un sistema de
ecuaciones diferenciales lineales es llamadoa veces autónomo. La
notación convencional general para el sistema de ecuaciones diferenciales
lineales es,

dx/ dt = f(t, x, y)

dy/ dt = g(t, x,

y)

El sistema anterior de ecuaciones diferenciales tendrá numerosas

funciones para satisfacerla. Mediante la modificación de la variable tiempo
obtendremos un conjunto de puntos que se encuentran en el plano de dos
dimensiones x-y, los cuales se denominan trayectoria. La velocidad con
respecto a esta trayectoria, en algún tiempo t es,

= (dx/ dt, dy/ dt)

Un ejemplo de un sistema de ecuaciones diferenciales lineales es el

siguiente, dx1/ dt = −4×1 + 2×2 dx2/ dt = 0×1 + −2×2

Con el fin de determinar el conjunto completo de fórmulas para la variable
dependiente de tiempo xi(t) para todos los valores de i, es necesario
obtener primero los vectores propios y valores propios de la matriz
constante A. En el caso que la matriz constante A posea un conjunto de
valores propios repetidos para sus componentes, sería necesarioun vector
propio generalizado.

Este es t, toma en cuenta que los vectores propios y valores propios de la

matriz constante puede ser un subconjunto de los números reales o
también un subconjunto de los números complejos.

La representación de la matriz del problema anterior es la siguiente, dx/ dt

= A * x En este caso, A es la matriz constante que puedeserrepresentada
como,

A=

Y x(t)T es un vector de variables definidas en términos de tiempo, el cual

es representado como,

x(t)T =

dx/ dt =
En caso de que el vector propio de la matriz constante A sea un
subconjunto de los números reales para este ejemplo, podemos escribir,

A = S * D * S-1

Aquí D es la matriz diagonal de la matriz de vectores propios de la matriz

constante A y S es la matriz que contiene los vectores propios en forma
de columnas, en el mismo orden como los valores propios se escriben en
la matriz diagonal D.

En consecuencia, la forma de la matriz del ejemplo anterior se puede

escribir como,

dx/ dt = A * x

dx1/ dt

dx2/ dt =

−4 2 0 −4

* x1 x2

4.1.2 Sistemas de EDL homogéneos

Sabemos que una ecuación diferencial lineal es de la forma,

Si esta misma ecuación se transforma en la forma,

Obtenemos una ecuación diferencial lineal homogénea. Esta se da cuando

la función conocida no está presente en la ecuación diferencial lineal,
entonces se le llama una ecuación diferencial homogénea. Y si tenemos
una gran cantidad de tales ecuaciones juntas, de manera tal que dependen
unas de las otras, y definen colectivamente un problema común, entonces
se les llama un sistema de ecuaciones diferenciales lineal es homogéneo.

Tales sistemas pueden ser resueltos de manera eficiente con la ayuda de

las matrices, las cuales son denominadas matriz fundamental. Sean X1,
X2 … X3 las soluciones de la matriz fundamental del sistema de entrada
de ecuaciones diferenciales homogéneas, entonces puede representarse
de manera condensada como,

En la ecuación anterior, las soluciones del sistema de ecuaciones

diferenciales están definidas en algún intervalo, digamos I y la solución
general del sistema de ecuaciones diferenciales es este
En la ecuación anterior, los términos que se mantienen dentro de los
corchetes son los vectores fila, donde X1 = [xi1j], X2 = [xi2j] …Xn = [xinj].
Estas son las soluciones n fundamentales del sistema de entrada de
ecuaciones diferenciales lineales homogéneas para el intervalo dado I.
Entonces tenemos que la matriz fundamental para el sistema homogéneo
de ecuaciones diferenciales lineales para el intervalo dado como I es,

Los pasos para resolver un sistema homogéneo de ecuaciones

diferenciales lineales son los siguientes:

1.Construye la matriz de coeficientes para las ecuaciones del sistema

dado.

2.Determina los valores propios de esta matriz de coeficientes del sistema

dado de ecuaciones diferenciales lineales homogéneas.

3.Ahora, busca el vector propio inicial de este conjunto de valores propios

y nómbralo como EV1.

4.Determina la primera ecuación de este vector en términos de constantes

k1, k2.

5.Después de esto, determina el siguiente conjunto de valores propios, sus

vectores propios correspondientes y su ecuación.
6.Anota la solución general para las ecuaciones en términos de constantes
k1, k2.

7.Por último, deriva la solución general para el sistema de ecuaciones.

Aunque el procedimiento para solucionar el sistema de ecuaciones

diferenciales lineales homogéneo es bastante fácil, se da un ejemplo
ilustrativo que te ayudará a hacer los conceptos más claros.

dx/ dt = 2x + 3y

dy/ dt = 2x + y

Primero, escribamos la matriz constante para el sistema de ecuaciones

diferenciales homogéneas dado. Esto es,

La matriz columna de los valores propios construida a partir de esta matriz

de coeficientes es la siguiente,

Esto nos da 1 = −1 y 2 = 4. A partir de estos valores propios el vector propio

asociado se construye como,
Colocando el valor de 1 = −1 en lugar, el valor exacto deEV1se obtiene
como,

El determinante de este se obtiene

como, | EV1 | = 0.

La ecuación asociada de este vector propio es,

3k1 + 3k2 = 0

2k1 + 2k2 = 0

De manera similar, la otra ecuación para el segundo vector propio es,

-2k1 + 3k2 = 0

2k1 - 3k2 = 0

Cuando t = 0, c1 = 1 y c2 = 1.

X1 = e-t K1
Del mismo modo,

Esto nos da la solución general,

4.1.3 Solución general y solución particular de sistemas de EDL

Solución general de los sistemas de ecuaciones diferenciales lineales y

solución particular de los sistemas de ecuaciones diferenciales lineales
Para resolver un sistema de ecuaciones diferenciales lineales, es
absolutamente esencial conocer los conceptos de valor propio y vector
propio. Para una matriz M dada, es llamadolos valores propios, si la
condición es verdadera,

Aquí x se llama vector propio de la matriz M.

Es decir, los vectores propios son aquellos vectores que luego de ser
multiplicados por la matriz de entrada permanecen proporcionales a la
matriz de entrada o resultan cero.

Sea A la matriz que contiene los valores propios, como [ 1, 2, 3 … n]. A

continuación se indican los pasos para resolver un sistema de ecuaciones
diferenciales.

1. Calcula las ecuaciones del sistema de ecuaciones diferenciales y

construye la matriz que contiene los coeficientes de todas las ecuaciones
en el orden en que aparecen en el sistema de entrada de la ecuación.

2. Los valores propios se obtienen de esta matriz que contiene los

términos de los coeficientes.

3. Calcula todos los vectores propios de valores propios como los

obtenidos en el paso anterior. Nómbralos en la secuencia a medida que
son determinados como EV1, EV2, EV3 …EVn.

4. Calcula las ecuaciones correspondientes para cada uno de los

conjuntos de valores propios y los vectores propios asociados. Repite este
paso para cada par de valores y vectores propios.

5. Obtén la solución particular para un sistema de ecuaciones no

homogéneo como,
Aquí X(t) se define como,X(t) = [x1 x2 x3 … xn]

Y en caso de que el sistema de entrada sea una ecuación diferencial

homogénea, entonces la solución particular del sistema será dada de la
forma,

La ecuación anterior nos da la relación,

En la relación anterior, y son valores propios y vectores propios,

respectivamente.

6. Y la solución general para el sistema de ecuaciones diferenciales

lineales es dada de la forma,

En general podemos decir que la solución de un sistema de ecuación

diferencial es llamada solución general si los valores de las constantes no
se obtienen en la solución final. La misma solución puede convertirse en
una solución particular cuando tenemos el valor de las constantes
determinadas. Esto se hace en el caso que el sistema de entrada de la
ecuación diferencial sea un problema de valor inicial con las condiciones
iniciales establecidas para la determinación de los términos constantes.

El ejemplo siguiente aclarará el procedimiento para resolver un sistema de

ecuaciones diferenciales lineales.

Determina el conjunto de ecuaciones como xT(t) = [(x1(t), x2(t)] para el

sistema de ecuaciones dx/ dt = A * x con las condiciones iniciales
establecidas como x(0) = x0
= (x01, x02). El valor de la matriz A estádada como,

La ecuación característica de la matriz de coeficientes

arriba es, f( ) = 2 –trace(A) * + det(A) = 2 + 4 + 4

Y las raíces de esta ecuación nos dan repetidos valores propios de la

matriz como, 1 = 2 = −2.

Entonces, el vector propio de la matriz es dado de la forma,

La matriz tiene un solo vector propio ya que ambos valores propios son los
mismos. Por lo tanto, la solución general del problema se da como,

Por consiguiente, un vector propio generalizado puede ser calculado

como, v = vp
+s1* vh1 v1 v2 −1 0 +s1* 1 1

La solución particular de este problema sería vT(p) = (−1, 0) = v2, el cual

es el valor propio generalizado de esta matriz, junto con los valores propios
repetidos y vh es la solución homogénea dando el vector propiov1.

Y la solución general del problema es,

Estonos da, x1(0) = c1 e0 + c2 (0 –

e0) = c1 – c2 = x01 x1(0) = c1 e0 +

c2 (0) = c1 = x02
4.2 Métodos de solución para sistemas de EDL

Un sistema de diferenciales lineales puede resolver las ecuaciones. Al

igual que existen varias técnicas para resolver una ecuación diferencial
lineal, también las hay para un sistema de ecuaciones diferenciales
lineales. Como el método de eliminación de Gauss, método separable y
reducible etc. Sea un sistema de ecuaciones diferenciales lineales
representado como,

Entonces, la representación de la matriz equivalente de este sistema de

ecuaciones diferenciales lineales será,
Ahora, la técnica más conveniente para resolver este sistema de
ecuaciones diferenciales lineales es recurrir al método de multiplicación de
matrices. La fórmula para resolverlo está dada de la forma,

Aquí A es la matriz que contiene los términos coeficientes de toda la

ecuación del sistema, C, que es una matriz columna compuesta de
elementos no homogéneos y, finalmente, la matriz X es la que contiene los
elementos desconocidos. Cuando la matriz C es igual a cero, entonces el
sistema dado es un sistema homogéneo de ecuaciones diferenciales
lineales.

Todo el mundo está familiarizado con el hecho de que multiplicar unas

cuantas matrices es mucho mejor que solucionar las ecuaciones
algebraicas crípticas para cuatro variables desconocidas. Sin embargo, la
técnica anterior nos da la solución en todo momento. Por lo tanto, tenemos
que adoptar algunas otras técnicas para resolver las ecuaciones.

Dos sistemas de ecuaciones diferenciales lineales se llaman equivalentes

en la situación de que ambos produzcan las mismas soluciones para
variables desconocidas. Sin embargo, esto puede lograrse mediante
aplicar algunas operaciones elementales como la multiplicación con una
constante,la reordenación de las ecuaciones, etc.Sea un sistema de
ecuaciones diferenciales lineales dado como,

• x+y+z=0
• x – 2y + 2z = 4
• x + 2y – z = 2

4..3 Método de los operadores

Un operador es un objeto matemático que convierte una función en otra,

por ejemplo, el operador derivada convierte una función en una función
diferente llamada la función derivada. Podemos de nir el operador derivada
D que al actuar sobre una función diferenciable produce la derivada:
0f(x) = f(x) ; D1f(x) = f0(x) ; D2f(x) = f00(x) ; : : : ;Dnf(x) = f(n)(x) :

Es posible construir la siguiente combinaci on lineal con los operadores

diferenciales:
P

(D) = a0 + a1D + a2D2 + + anDn ; an 6= 0 : (1)

4.4Utilizando transformada de Laplace.

4.5 Aplicaciones

Los sistemas de ecuaciones diferenciales lineales encuentran sus

aplicaciones en varios problemas que surgen en el sistema del mundo real.
Algunos de estos problemas se discuten a continuación.

1. Problema mecánico del acoplamiento de los resortes: Dos cuerpos con

masa m1, m2, respectivamente, yacen sobre una mesa. La mesa está libre
de fricción. Los dos cuerpos están conectados entre sí con la ayuda de un
resorte. Esteresorte está en una posición no estirada. También cada uno
de estos cuerpos está conectado a una superficie estática con la ayuda de
los resortes. Una vez más, estos resortes no están estirados. La constante
elástica de cada uno de los resortes es k1, k2, k3, respectivamente. La
situación anterior puede ilustrarse como,
Aquí O1 es la posición inicial del primer cuerpo y O2 es la posición inicial
del segundo cuerpo. Los cuerpos pueden ser cambiados de su posición
de equilibrio mediante mover cualquiera delos cuerpos en cualquier
dirección y luego soltarlos. Unejemplo de estoes,

En la figura anterior, x1 es la cantidad de distancia recorrida por el primer

cuerpo cuando este se mueve desde la posición de equilibrio y x2 es la
cantidad de distancia recorrida por el segundo cuerpo cuando este se
mueve desde la posición de equilibrio. Esto implica que el primer resorte
se alarga desde la posición estática por una distancia de x1 y el segundo
resorte se alarga desde la posición estática por una distancia de x2 –
x1.Esto implica que dos fuerzas restauradoras están actuando sobre el
primer cuerpo, estas son:

• La fuerza del primer resorte la cual actúa en dirección izquierda. Esta

fuerza por la ley de Hookes igual ak1×1.
• La fuerza del segundo resorte que actúa en dirección derecha. Esta
fuerza es igual a k2(x2 – x1).
Unidad 5.- Introducción a las series de Fourier
5.1 Teoría preliminar
Una sinusoide es una señal de la forma

𝐴𝑠𝑒𝑛 (2𝜋𝑣𝑡 + 𝜑)

El número A > 0 es la amplitud, v > 0 es la frecuencia medida en ciclos por

segundo o Hercios (Hz), -π < φ 6 π (fase inicial), ω = 2πv es la frecuencia en
𝑟𝑎𝑑⁄𝑠, (que se llama a veces frecuencia angular). El periodo es el tiempo
que necesita la sinusoide para completar un ciclo completo, es decir, el
periodo es T = 1⁄𝑣 segundos.

𝐴𝑠𝑒𝑛 (2𝜋𝑣(𝑡 + 1⁄𝑣) + 𝜑) = 𝐴𝑠𝑒𝑛 (2𝜋𝑣𝑡 + 2𝜋 + 𝜑) = 𝐴𝑠𝑒𝑛 (2𝜋𝑣𝑡 + 𝜑) En

general, una función f: R C se dice que es periódica con periodo T si 𝑓(𝑡 +
𝑇) = 𝑓(𝑡) para todo
𝑡 ∈ 𝑅. En tal caso cualquier múltiplo entero de T es también un periodo f,
esto es,
𝑓(𝑡 + 𝑘𝑇) = 𝑓(𝑡) Para todo 𝑡 ∈ 𝑅, 𝑘 ∈ 𝑍.

Por convenio, una función constante se considera periódica con cualquier

periodo. Salvo este caso, cuando se dice que una función es periódica de
periodo T, se sobreentiende que T es el numero positivo, más pequeño que
verifica la igualdad 𝑓(𝑡 + 𝑇) = 𝑓(𝑡) para todo 𝑡 ∈ 𝑅. En la representación
gráfica de la señal 𝑓(𝑡) = 𝐴𝑠𝑒𝑛 (2𝜋𝑣𝑡 + 𝜑) se interpreta 𝑓(𝑡) como la
amplitud de señal en el instante t.

La amplitud A representa la máxima altura que alcanza dicha gráfica, esto

es, el máximo absoluto de la función f (el mínimo absoluto es –A). La
frecuencia es el número de veces (ciclos) que se repite la gráfica en un
segundo. El periodo es el tiempo necesario para que la gráfica complete un
solo ciclo.
5.2 Series de Fourier

La idea básica de las series de Fourier es que toda función periódica de

periodo T puede ser expresada como una suma trigonométrica de senos y
cosenos del mismo periodo T.

Algunas funciones periódicas de f (t) de periodo T pueden expresarse por la

siguiente serie, llamada “Serie trigonométrica de Fourier”

𝑓(𝑡) = 1 2 𝑎0 + ∑[ ∞ 𝑛=1 𝑎𝑛 cos(𝑛𝜔0𝑡) + 𝑏𝑛𝑠𝑒𝑛(𝑛𝜔0𝑡)]

Donde 𝜔0 = 2π/T y 𝑎0, 𝑎𝑛 y 𝑏𝑛 son los coeficientes de Fourier que toman

los valores:

𝑎0 = 2 𝑇 ∫ 𝑓(𝑡)𝑑𝑡 𝑇/2 −𝑇/2

𝑎𝑛 = 2 𝑇 ∫ 𝑓(𝑡) 𝑇/2 −𝑇/2 𝑐𝑜𝑠(𝑛𝜔0𝑡)𝑑𝑡 𝑛 = 1, 2, 3, …

𝑏𝑛 = 2 𝑇 ∫ 𝑓(𝑡) 𝑇/2 −𝑇/2 𝑠𝑒𝑛(𝑛𝜔0𝑡)𝑑𝑡 𝑛 = 1, 2, 3, …

Se dice que las funciones del conjunto {𝑓𝑘(t)} son ortogonales en el intervalo
a<t

∫ 𝑓𝑚(𝑡)𝑓𝑛(𝑡) 𝑏 𝑎 𝑑𝑡 = { 0 𝑝𝑎𝑟𝑎 𝑚 ≠ 𝑛 𝑟𝑛 𝑝𝑎𝑟𝑎 𝑚 = 𝑛 }

Forma compleja de serie de Fourier

Es posible obtener una forma alternativa usando las fórmulas de Euler:

𝑐𝑜𝑠(𝑛𝜔0𝑡) = 1 2 (𝑒 𝑖𝑛𝜔0 + 𝑒 −𝑖𝑛𝜔0)

𝑠𝑒𝑛(𝑛𝜔0𝑡) = 1 2𝑖 (𝑒 𝑖𝑛𝜔0 − 𝑒 −𝑖𝑛𝜔0)

Sustituyendo

𝑓(𝑡) = 1 2 𝑎0 + ∑[ ∞ 𝑛=1 𝑎𝑛 1 2 (𝑒 𝑖𝑛𝜔0 + 𝑒 −𝑖𝑛𝜔0) + 𝑏𝑛 1 2𝑖 (𝑒 𝑖𝑛𝜔0 −

𝑒 −𝑖𝑛𝜔0)]

Y usando el hecho de que 1/i = -i:

𝑓(𝑡) = 1 2 𝑎0 + ∑[ ∞ 𝑛=1 𝑎𝑛 1 2 (𝑎𝑛 − 𝑖𝑏𝑛 )𝑒 𝑖𝑛𝜔0𝑡 + 1 2 (𝑎𝑛 + 𝑖𝑏𝑛 )𝑒

−𝑖𝑛𝜔0𝑡 ]

Y definiendo,

𝐶0 ≡ 1 2 𝑎0

𝐶𝑛 ≡ 1 2 (𝑎𝑛 − 𝑖𝑏𝑛 )

𝐶−𝑛 ≡ 1 2 (𝑎𝑛 + 𝑖𝑏𝑛 )𝑒 −𝑖𝑛𝜔0

Queda:

𝑓(𝑡) = ∑ 𝐶𝑛𝑒 𝑖𝑛𝜔0𝑡 ∞ 𝑛=−∞

A esta se le llama forma compleja de la serie de Fourier y su coeficiente 𝐶𝑛

se puede obtener de lo siguiente:

𝐶𝑛 = 1 𝑇 ∫ 𝑓(𝑡) 𝑇 0 𝑒 −𝑖𝑛𝜔0𝑡𝑑

Convergencia de una serie de Fourier

Existe un teorema que a continuación se mencionara, que especifica las
condiciones suficientes de convergencia de una serie de Fourier en un punto.

Teorema condiciones de convergencia

Sean f y f´ continuas en tramos en el intervalo (-p, p) esto es, sean continuas

excepto con un numero finito de puntos en el intervalo y con discontinuidades
solo finitas en esos puntos. Entonces la seria de Fourier de f en el intervalo
converge hacia f(x) en un punto de continuidad. En un punto de
discontinuidad la serie de Fourier converge hacia el promedio f(x+)+f(x−) 2
en donde f(x+) y f(x) representan el límite de f en x, desde la derecha y la
izquierda, respectivamente.

Extensión periódica

Las funciones de conjunto básico tiene un periodo común de 2p, por eso se
dice que el lado derecho de la ecuación se le denomina periodo. Y con esto
podemos decir que una serie de Fourier no solo representa a la función en
el intervalo (-p, p), sino que también de la extensión periódica de f fuera en
el intervalo. Aplicando el teorema a la extensión periódica f, o desde el
principio suponer, que la función dada es periódica con un periodo de 2p. Y
cuando f es continúa por tramos y converge hacia el primero [f(p−)+f(p+)] 2
en esos extremos y hasta el valor extendido periódicamente a 3p, 5p, 7p,
etc.
5.3 Las series de Fourier en cosenos, senos y de medio intervalo

Series de Fourier en cosenos

Una serie de Fourier en cosenos no es más que extraer una función definida
en intervalos con una función par. La función lo es, si 𝑓(𝑥) = 𝑓(−𝑥). Ahora
que 𝑓(𝑥) es una función par en [−𝜋, 𝜋] y se extienda 𝑓(𝑥) a todo el intervalo
[−∞, ∞], originando que 𝑓(𝑥) de periodo 2π de tal manera que la serie de
Fourier de la función 𝑓(𝑥) sea:

𝑎0 = 1 𝜋 ∫ 𝑓(𝑥)𝑑𝑥 𝜋 −𝜋 = 0

𝑎𝑛 = 1 𝜋 ∫ [𝑓(𝑥) cos 𝑛 𝑥]𝑑𝑥 𝜋 −𝜋 = 0

La serie de Fourier, solo tiene términos en senos, es decir, será de la forma

𝑓(𝑥) = 𝑏1𝑠𝑒𝑛 𝑥 + 𝑏2𝑠𝑒𝑛 2𝑥 + 𝑏3𝑠𝑒𝑛 3𝑥 + ⋯

Series de Fourier en senos

Esto es extender el comportamiento de una función definida en medio de

intervalos como una función impar. Esta lo es si 𝑓(−𝑥) = −𝑓(𝑥)

𝑏𝑛 = 1 𝜋 ∫ [𝑓(𝑥)𝑠𝑒𝑛 𝑥] 𝑑𝑥 = 0 𝜋 −𝜋

La serie de Fourier correspondiente solo tiene el términos independiente y

términos en coseno, es decir, será de la forma

𝑓(𝑥) = 𝑎0 + 𝑎1 cos 𝑥 + 𝑎2 cos 2𝑥 + ⋯

Series de Fourier en medio intervalo

Definiendo la función f en –π < x < 0 como 𝑓(𝑥) = 𝑓(𝑥 + 𝜋) como se muestra
en la gráfica 3.

En la definición de serie de Fourier, de funciones pares o impares, solo se

utiliza la mitad del intervalo, es decir de 0 < x < π, por lo tanto en la práctica
no hay necesidad de reflejar la función haciéndola par o impar, se define la
función en la mitad del intervalo a partir del origen. Esto se conoce como
desarrollo en mitad del intervalo.

La función 𝑓(𝑥) es alternada: es decir 𝑓(𝑥 + 𝜋) = −𝑓(𝑥), las cuales son un

caso corriente en electrotecnia. En este caso la serie de Fourier
correspondiente solo tiene términos de senos y cosenos impares, ya que los
pares se anulan, en efecto

𝑎𝑛 = 1 𝜋 ∫ [𝑓(𝑥) cos 2𝑛 𝑥]𝑑𝑥 𝜋 −𝜋 = 1 𝜋 {∫ [𝑓(𝑥) cos 2𝑛 𝑥]𝑑𝑥 0 −𝜋 + ∫ [𝑓(𝑥)

cos 2𝑛 𝑥] 𝜋 0 𝑑𝑥}