1-Teoría - 2021 Matemática Avanzada
1-Teoría - 2021 Matemática Avanzada
3. Análisis de Fourier 75
3.1. Introducción y motivación . . . . . . . . . . . . . . . . . . . . . . 75
3.2. Sucesiones ortogonales . . . . . . . . . . . . . . . . . . . . . . . . 77
3.3. Series de Fourier generalizadas . . . . . . . . . . . . . . . . . . . 79
3.4. Series trigonométrica y exponencial de Fourier . . . . . . . . . . 83
3.5. Respuesta permanente a funciones periódicas . . . . . . . . . . . 90
3.6. Simetrı́a de ondas periódicas . . . . . . . . . . . . . . . . . . . . 91
3.7. Espectro discreto de frecuencias . . . . . . . . . . . . . . . . . . . 94
3.8. Propiedades del espectro discreto . . . . . . . . . . . . . . . . . . 99
3.9. Espectro de potencia . . . . . . . . . . . . . . . . . . . . . . . . . 102
3.10. Transformada de Fourier . . . . . . . . . . . . . . . . . . . . . . . 103
3.11. Transformada de algunas señales sencillas . . . . . . . . . . . . . 105
3.12. Transformada de Fourier de funciones absolutamente integrables 109
3.13. Propiedades de la transformada de Fourier . . . . . . . . . . . . . 111
3.14. Espectro de energı́a . . . . . . . . . . . . . . . . . . . . . . . . . . 117
3.15. Transmisión de señales a través de filtros lineales (SLIT) . . . . . 118
3.16. Transformada de Fourier de funciones no absolutamente integrables120
3.17. Transformadas de Fourier seno y coseno . . . . . . . . . . . . . . 124
1
ÍNDICE GENERAL 2
3
ÍNDICE GENERAL 4
Integrantes de la cátedra
Dra. Gabriela Messineo (responsable de cátedra)
Ing. Carlos Chiuro (JTP)
Dra. Vanesa Galli
Ing. Mg. Alberto López
Régimen de promoción
La evaluación de la asignatura a lo largo del cursado se realizará en las siguientes
instancias:
2 exámenes parciales teórico prácticos que serán calificados en la escala de
1 a 100 puntos. La ausencia a un parcial significa cero como nota.
Seminarios con Matlab u Octave (software libre), que deberán realizar-
se antes de cada parcial. Serán de asistencia obligatoria, siendo ésta la
condición indispensable para poder rendir el parcial correspondiente.
1 único recuperatorio que puede utilizarse para reemplazar la nota del
parcial con menor puntaje o para compensar una inasistencia justificada
a alguna de las instancias de evaluación regulares.
1 coloquio para definir la calificación final en los casos en que la cátedra
lo considere necesario.
El cierre de notas es posterior a la finalización de la cursada y hasta esa fecha
la cátedra podrá realizar evaluaciones.
- Las inasistencias se consideran justificadas si se presentan certificados médicos
o de trabajo.
- NO están justificadas la inasistencias por viajes.
Programa analı́tico
1. Variable Compleja
a) Transformación conforme Transformación lineal. Transformación
inversa. Transformación bilineal.
b) Integración. Polos y residuos. Integrales en el campo complejo.
Singularidades aisladas. Definición de residuos. Teorema de Cauchy-
Goursat. Consecuencias. Teorema de los residuos. Fórmula de la in-
tegral de Cauchy y de la Derivada de la integral de Cauchy.
c) Series de potencias. Series de Taylor. Series de Laurent. Méto-
dos prácticos. Desarrollos en serie. Convergencia. Fórmulas para el
cálculo de residuos.
2. Análisis de señales y sistemas en el dominio natural.
a) Clasificación y propiedades de Señales. Análisis de señales en
el dominio del tiempo. Señales periódicas y no periódicas. Señales de
energı́a y de potencia. Tipos de señales. Función impulso. Función
impulso como lı́mite de otra función. Propiedades.
b) Sistemas lineales e invariantes en el tiempo. Ejemplos de Mecáni-
ca y Electrotecnia. Función operacional del sistema. Solución de ecua-
ciones diferenciales en el dominio del tiempo. Solución transitoria y
permanente de un sistema estable. Aplicaciones a sistemas lineales.
Diagramas en bloques.
c) Convolución y respuesta al impulso. Cálculo y propiedades de la
integral de convolución. Cálculo de la respuesta al impulso. Solución
libre y forzada de un sistema lineal e invariante. Respuesta a funciones
exponenciales. Estabilidad de un sistema. Relación entre la respuesta
al escalón y la respuesta al impulso. Ejercicios de aplicación.
3. Análisis de señales y sistemas en el dominio transformado.
a) Series de Fourier. Sistemas ortogonales. Serie generalizada de Fou-
rier. Aproximación cuadrática. Coeficientes de Fourier. Identidad de
Parseval. Series trigonométrica y exponencial de Fourier. Simetrı́a de
la forma de onda. Integración y diferenciación de las series de Fou-
rier. Espectro de frecuencia discreta. Espectro de amplitud y de fase.
Propiedades. Espectro de potencia. Respuesta de un sistema lineal a
una función periódica.
ÍNDICE GENERAL 6
Bibliografı́a
E.KREYSZIG: Matemáticas avanzadas para Ingenierı́a. Ed. Limusa.
R.V.CHURCHILL; J.W.BROWN: Variable compleja y aplicaciones. Ed. Mc
Graw Hill.
M.R.SPIEGEL: Variable compleja. Serie Schaum. Ed. Mc Graw Hill.
R.A.GABEL, R.A. ROBERTS: Señales y Sistemas lineales. Ed. Limusa S.A.
C.D.MC GILLEN, H.R.COOPER: Continuos and Discrete Signal and System
Analysis. Ed.Holt-Rinehart and Winston.
R.V.CHURCHILL: Series de Fourier y problemas de contorno. Ed. Mc Graw
Hill.
H.R.HSU: Análisis de Fourier. Fondo educativo interamericano S.A.
M.R.SPIEGEL: Transformada de Laplace. Serie Schaum. Mc Graw Hill.
A.PAPOULIS: The Fourier Integral and its Aplications. Mc Graw Hill.
M.R.SPIEGEL: Análisis de Fourier. Serie Schaum. Ed. Mc Graw Hill.
OPPENHEIN –WILLSKY: Señales y Sistemas. Prentice Hall.
C.L.PHILLIPS - J.M.PARR: Signals, Systems, and Transforms. Prentice Hall.
http://www.jhu.edu/ signals
Capı́tulo 1
1.1. Introducción
Función analı́tica
La función ω = f (z), definida para los números complejos z = x+jy, es analı́tica
en un punto dado z0 ∈ D si la misma es derivable tanto en el propio punto z0
como en un cierto entorno del mismo.
Es condición necesaria (aunque no suficiente) para la derivabilidad de este tipo
de funciones que se verifiquen las Condiciones de Cauchy-Riemann (C-R), dos
ecuaciones diferenciales parciales, básicas en el análisis de funciones complejas
de variable compleja.
Recordamos las condiciones de C-R para f (z) = x(u, v) + jy(u, v):
∂u ∂v ∂v ∂u
= ; =− ,
∂x ∂y ∂x ∂y
Singularidades
El punto z0 es una singularidad de f (z) si f (z) es analı́tica en algún punto de
cierto entorno de z0 pero no lo es en el propio punto.
7
CAPÍTULO 1. TEMAS DE VARIABLE COMPLEJA 8
ϕ(z)
f (z) =
(z − z0 )m
lı́m f (z)(z − z0 )m = k; k 6= 0.
z→z0
Singularidad esencial
z = z0 es una singularidad esencial de f (z) si no existe lı́m f (z)
z→z0
Singularidad evitable
z = z0 es una singularidad evitable de f (z) si lı́m f (x) = L; L 6= 0.
z→z0
Ejercicios
Singularidades
Analizar las singularidades de cada función:
sen(z)
1. f (z) = z en z0 = 0; Rta.: Singularidad evitable
sen(z)
2. f (z) = z3 en z0 = 0; Rta.: Polo de segundo orden
CAPÍTULO 1. TEMAS DE VARIABLE COMPLEJA 9
1−e−z
3. f (z) = z en z0 = 0; Rta.: Singularidad evitable
1
4. f (z) = e z en z0 = 0; Rta.: Singularidad esencial
1−cos(z)
5. f (z) = z7 en z0 = 0; Rta.: Polo de quinto orden
Ceros
Determinar el orden de los ceros de las siguientes funciones:
6. f (z) = (z − a)3 ; Rta.: z = a es un cero de tercer orden de f (z).
7. f (z) = 1 − ez ; Rta.: z = 2πkj son ceros simples de f (z).
C1
C1*
a
C2
Plano Z
C2*
Plano W
Transformacion lineal
Por lo tanto:
u + jv = x + jy + br + jbi = x + br + j(y + bi ),
u = x + br
de donde salen la parte real y la parte imaginaria de w:
v = y + bi
Las ecuaciones de arriba son las denominadas coordenadas de transforma-
ción y permiten ver que esta transformación implica un corrimiento de br
unidades en el eje de abscisas y bi unidades en el eje de ordenadas.
B = 0; w = Az
Como tenemos el producto de dos número complejos (z y A), en este caso
conviene expresarlos en forma polar.
z = rejθ , A = Rejα .
Además, w = ρejφ .
Por lo tanto:
ρejφ = Rejα rejθ = Rrej(α+θ) ,
Transformación inversa
1 1
Su forma general es w = z oz= w.
En coordenadas polares
ρ = 1r
1
ρejφ = ⇒
rejθ φ = −θ
Puede obsevarse que existe una simetrı́a con respecto al eje real, ya que las fases
son opuestas. Entonces, si la fase de cierto punto z0 es, por ejemplo, π/6, la
CAPÍTULO 1. TEMAS DE VARIABLE COMPLEJA 11
En coordenadas cartesiana
1
w=
z
x
1 x − jy u = x2 +y2
u + jv = = 2 ⇒
x + jy x + y2
v = x2−y
+y 2
1
Análogamente, para z = w,
u −v
x= ;y = 2
u2 +v 2 u + v2
a(x2 + y 2 ) + bx + cy + d = 0
u2 + v 2 u v
a 2 2 2
+b 2 2
−c 2 +d=0
(u + v ) u +v u + v2
a + bu − cv + d(u2 + v 2 ) = 0
Puede encontrarse una transformación bilineal que transforme tres puntos cua-
lesquiera, z1 , z2 , z3 en tres puntos w1 , w2 , w3 , respectivamente, mediante la fórmu-
la:
(z − z1 )(z2 − z3 ) (w − w1 )(w2 − w3 )
=
(z − z3 )(z2 − z1 ) (w − w3 )(w2 − w1 )
Prueba:
La fórmula anterior puede escribirse como:
Ejercicios
Transformación lineal
1. Se desea encontrar la imagen de la región 0 ≤ x ≤ 1; 0 ≤ y ≤ 2 mediante
la transformación w = (1 + j)z + 2 − j. ¿Cómo afectan los coeficientes A
y B a la región original?
Respuesta: A = 1 + j, B = 2 − j
- Rotación dada por el argumento de A: π/4
√
- Magnificación dada por el módulo de A: 2
- Traslación según B: 2 unidades a la derecha y 1 para abajo.
Transformación inversa
2. Transformar x2 + y 2 − 4x + 2y = 0 mediante la transformación inversa.
Respuesta: v = −2u + 21 .
3. Hallar analı́tica y gráficamente mediante la transformación inversa, la ima-
gen de la región sombreada:
Respuesta:
u≤0
u2 + (v + 1)2 ≤ 1
u + (v + 12 )2 ≥ 21
2
Respuesta: w = 4.
Transformación bilineal
5. Encontrar la transformación que mapea: z1 = 1, z2 = 0, z3 = −1 en los
puntos w1 = j, w2 = 1, w3 = ∞. (Sugerencia: utilice la definición del punto
del infinito)
1+z(2j−1)
Respuesta: w = z+1
CAPÍTULO 1. TEMAS DE VARIABLE COMPLEJA 14
6. Mediante w = z−j
z+j , transformar x > 0, y > 0.
v<0
Respuesta:
2
u + v2 < 1
7. Hallar la transformación bilineal del cı́rculo |z| < 5 en el cı́rculo |w| < 1,
de tal modo que los puntos z1 = −5, z2 = 4 + 3j y z3 = 5 se transformen
en los puntos w1 = −1, w2 = j y w3 = 1, respectivamente.
2z−5
Respuesta: w = 10−z
Teorema de Green
Figura 1.3
Consideremos una función f (z) = u(x, y)+jv(x, y), analı́tica en todos los puntos
interiores y sobre un contorno cerrado C, y es tal que es continua allı́. Se quiere
evaluar la integral curvilı́nea: Z
f (z)dz.
C
Es fácil ver que las integrales de lı́nea complejas pueden expresarse en términos
de integrales de lı́nea reales, en función de sus componentes. Si se sustituye
f (z) = u(x, y) + jv(x, y) y dz = dx + jdy en la integral y se opera, resulta:
Z Z Z
(u(x, y)+jv(x, y))(dx+jdy) = u(x, y)dx−v(x, y)dy+j v(x, y)dx+u(x, y)dy.
C C C
Teorema de Cauchy-Goursat
Sea f (z) una función analı́tica en un dominio simplemente conexo D. Entonces,
dado un contorno cerrado simple C contenido en D, tenemos:
Z
f (z)dz = 0.
C
Demostración:
Recordemos que podemos escribir a la integral compleja como:
Z Z Z
f (z)dz = u(x, y)dx − v(x, y)dy + j v(x, y)dx + u(x, y)dy.
C C C
Además, por ser f (z) analı́tica, u, v y sus derivadas parciales son continuas en
D. Por lo tanto, es posible aplicar el Teorema de Green y la ecuación anterior
queda:
Z ZZ ZZ
∂u(x, y) ∂v(x, y) ∂u(x, y) ∂v(x, y)
f (z)dz = − ( + )dxdy+j ( − )dxdy.
C R ∂y ∂x R ∂x ∂y
Esta demostración, basada en el teorema de Green, exige que f 0 (z) sea continua
en R, ya que de lo contrario no podrı́amos aplicar dicho teorema.
Cauchy obtuvo en 1814 por primera vez este resultado valiéndose de una fórmula
equivalente, ya que Green no habı́a aún explicitado su teorema.
Otra demostración, menos restrictiva, fue formulada a fines del siglo XIX por
Goursat, que no requiere que f 0 (z) sea continua. Estas deducciones se conocen
como Teorema de Cauchy-Goursat, o a veces únicamente como teorema de la
integral de Cauchy.
Consecuencia 1
Supongamos que f (z) es analı́tica en una región comprendida entre dos curvas
C y C1, es decir, en un dominio doblemente conexo como el que se muestra
sombreado en la figura.
En ese caso puede probarse que:
Z Z
f (z)dz = f (z)dz
C C1
CAPÍTULO 1. TEMAS DE VARIABLE COMPLEJA 17
Demostración:
El dominio doblemente conexo puede transformarse en uno simplemente conexo
con un corte que una los puntos A y B, y entonces se puede aplicar el teorema
de Cauchy-Goursat:
Z Z A Z Z B
f (z)dz + f (z)dz + f (z)dz + f (z)dz = 0
C B −C1 A
Obsérvese que el contorno del dominio simplemente conexo está dado por la
¯ la
curva C recorrida en sentido positivo a partir del punto B, el segmento BA,
¯
curva C1 recorrida en sentido negativo y el segmento AB.
RA RB R R
Como B f (z)dz = − A f (z)dz y −C1 f (z)dz = − C1 f (z)dz, podemos escri-
bir Z Z B Z Z B
f (z)dz − f (z)dz − f (z)dz + f (z)dz = 0,
C A C1 A
En ese caso
Z Z Z Z
f (z)dz = f (z)dz + f (z)dz + ... + f (z)dz
C C1 C2 Cn
Consecuencia 2
La segunda consecuencia del teorema de Cauchy-Goursat establece el principio
de independencia de la trayectoria.
R R
Luego, como BA≡−C1
f (z)dz = − AB≡C1 f (z)dz, tenemos que
Z Z
f (z)dz = f (z)dz.
AB≡C2 AB≡C1
Residuos
Definición
Sea f (z) una función analı́tica en un contorno cerrado C simplemente conexo y
en todo punto del interior de C, salvo en z0 . Se denomina residuo de la función
f (z) en una singularidad aislada z = z0 al número definido por:
Z
1
Resz=z0 f (z) = f (z)dz
2πj C
Toda función tiene un residuo en cada uno de sus puntos singulares aislados.
Sin embargo, ese residuo puede ser cero.
- El residuo de una singularidad evitable es siempre cero.
- Cuando las singularidades son polos el cálculo de los residuos puede hacerse
utilizando las fórmulas que se dan abajo, y que demostraremos luego de estudiar
las series de potencias.
Polos simples
Si z0 es un polo simple tenemos dos fórmulas para calcular el residuo de la
función en ese punto.
Primero, recordemos que en este caso la función puede escribirse como f (z) =
φ(z)
z−z0 , con φ(z) analı́tica en z0 y φ(z0 ) 6= 0. Probaremos más adelante que
P (z)
La segunda fórmula puede usarse cuando f (z) = Q(z) , con P (z) y Q(z) analı́ticas
0
en z0 , P (z0 ) 6= 0, Q(z0 ) = 0 y Q (z0 ) 6= 0. Entonces
P (z0 )
Resz=z0 f (z) =
Q0 (z0 )
Polos múltiples
Si z0 es un polo de multiplicidad m, entonces la función puede escribirse como
φ(z)
f (z) = (z−z 0)
m , con φ(z) analı́tica en z0 . Entonces puede probarse que
φ(m−1) (z0 )
Resz=z0 f (z) =
(m − 1)!
Demostración:
Encerramos a cada uno de los puntos singulares zi en un cı́rculo Ci con radio
pequeño para que queden separados todos esos m cı́rculos y C. Entonces f (z)
es analı́tica en el dominio multiplemente conexo limitado por C, C1 , C2 , ..., Cm
y sobre la frontera.
Figura 1.6
1
R
Además, vimos que Resz=zi f (z) = 2πj Ci
f (z)dz. Por lo tanto,
Z
f (z)dz = 2πjResz=z1 f (z) + 2πjResz=z2 f (z) + ... + 2πjResz=zm f (z),
C
El teorema de los residuos junto con los diferentes métodos para el cálculo de
residuos nos permiten calcular integrales en el plano complejo cuando la función
integrando no cumple con las condiciones requeridas por el teorema de Cauchy-
Goursat.
CAPÍTULO 1. TEMAS DE VARIABLE COMPLEJA 21
Figura 1.7
φ(z)
La función z−z 0
es analı́tica en todos los puntos interiores y sobre C, excepto en
el punto z0 . Por lo tanto, es analı́tica en la región doblemente conexa entre ambos
contornos y, según la primera consecuencia del teorema de Cauchy-Goursat:
Z Z
ϕ(z) ϕ(z)
dz = dz,
C z − z0 C0 z − z0
ϕ(z) − ϕ(z0 )
Z Z Z
ϕ(z) 1
dz = dz + ϕ(z0 ) dz.
C z − z0 C0 z − z0 C0 z − z0
CAPÍTULO 1. TEMAS DE VARIABLE COMPLEJA 22
Puesto que puede ser tan pequeño como se desee, el valor absoluto de la
integral también puede hacerse arbitrariamente pequeño. Reducir equivale
simplemente a reducir el radio de C0 .
En cuanto a la segunda integral, vamos a escribir los números en coordenadas
polares, con lo que la integral queda en la variable θ.
ϕ(z) − ϕ(z0 )
Z Z Z
ϕ(z) 1
dz = dz + ϕ(z0 ) dz
C z − z0 C0 z − z0 C0 z − z0
Z
ϕ(z)
dz = 0 + 2πjϕ(z0 )
C z − z0
Z
1 ϕ(z)
ϕ(z0 ) = dz.
2πj C (z − z0 )
Ejercicios
ez dz = 0 para toda trayectoria cerrada en el plano complejo.
R
1. Probar que C
z 2 −2z+3
2. Determinar el residuo de f (z) = z−2 en sus singularidades.
Respuesta: Resz=2 f (z) = 3.
R cos(z)
Calcule C z−1 dz; Respuesta: cos(1)
R cos(z)
Calcule C z+1 dz; Respuesta: 0
8. Reconsidere la función del ejercicio 6 y calcule la integral a lo largo del
cı́rculo de radio 2 con centro en 0. Discuta el resultado en relación al
obtenido en el ejercicio 6 cuando a = j.
Respuesta: 0.
9. Considere ahora la misma curva que en el ejercicio 8 para integrar dos
nuevas funciones:
z
R
C z 2 −1
dz; Respuesta: 2πj
R 3z2 +1
C z(z−1)
dz; Respuesta: 6πj
CAPÍTULO 1. TEMAS DE VARIABLE COMPLEJA 24
2z+1
R
10. Calcule C:|z|=4 (z−1)(z+2)(z−3)
dz; Respuesta: 0
ch(z)
Respuesta: − πj
R
12. C:|z|=2 (z+1) 3 (z−1) dz; 2e
Se sabe que las series de potencias convergen para ciertos valores de la variable.
Las series de variables real convergen para un cierto intervalo. En cambio, en las
series de variable compleja lo que tenemos es una región de convergencia (RC)
en el plano complejo.
La región de convergencia queda determinada como un cı́rculo de centro a y
radio R, tal que en |z − a| < R la serie es convergente y en |z − a| > R es
divergente. El cı́rculo se llama cı́rculo de convergencia y su radio recibe el nombre
de radio de convergencia.
El radio se calcula a partir de los coeficientes de la serie, y según la fórmula de
Cauchy-Hadamard:
1
R= ,
L
donde
cn+1
L = lı́m .
n→∞ cn
Serie de Taylor
Sea f (z) analı́tica en todos los puntos interiores a una circunferencia C0 con
centro en z0 y radio R. En cada punto z interior a C0 la función queda
representada por una serie de Taylor, que tiene la forma:
f 00 (z0 ) f n (z0 )
f (z) = f (z0 ) + f 0 (z0 )(z − z0 ) + (z − z0 )2 + ... + (z − z0 )n + ... + Rn
2! n!
Serie de Laurent
Si una función f (z) no es analı́tica en z0 , no podemos aplicar el teorema de
Taylor en ese punto. No obstante, muchas veces es posible hallar una repre-
sentación de f (z) en forma de una serie que contiene tanto potencias positivas
como negativas de z − z0 . Es la denominada serie de Laurent.
y Z
1 f (z)
bn = dz; n = 1, 2, ...
2πj C (z − z0 )−n+1
ϕ(z)
Si z0 es un polo de orden m, entonces f (z) = (z−z0 )m , con ϕ(z) analı́tica en z0 .
Reescribiendo:
Z
1 ϕ(z) 1
bn = dz; n = 1, 2, ...,
2πj C (z − z0 )m (z − z0 )−n+1
vemos que cuando n > m los coeficientes se hacen cero porque el integrando es
una función analı́tica en z = z0 y en ese caso se verifica el teorema de Cauchy-
Goursat.
Singularidad esencial en z = z0
Decimos que el punto z0 es un punto singular esencial de la función, si la parte
principal contiene un número infinito de términos.
Singularidad evitable en z = z0
Si la Serie de Laurent carece de parte principal, y la función no es analı́tica en
z0 pero puede hacerse analı́tica utilizando una adecuada definición de la misma,
entonces el punto z0 es un punto singular evitable.
ϕ(z)
f (z) = ,
(z − z0 )m
lı́m f (z)(z − z0 )m = k, k 6= 0.
z→z0
CAPÍTULO 1. TEMAS DE VARIABLE COMPLEJA 28
Demostración
Si z0 es un polo de orden m, existe lı́mz→z0 f (z)(z − z0 )m = k, k 6= 0.
La serie de Laurent de f (z) es:
∞
X b1 b2 bm−1 bm
f (z) = an (z−z0 )n + + +...+ + .
n=0
(z − z0 ) (z − z0 )2 (z − z0 )m−1 (z − z0 )m
lı́m f (z)(z − z0 )m = bm , bm 6= 0.
z→z0
lı́m f (z)(z − z0 )m = k, k 6= 0,
z→z0
lı́m |f (z)| = ∞.
z→z0
Resz=z0 f (z) = b1 .
P (z0 )
Resz=z0 f (z) = .
Q0 (z0 )
P (z)
Resz=z0 f (z) = limz→z0 f (z)(z − z0 ) = limz→z0 (z − z0 )
Q(z)
P (z) P (z)(z − z0 )
limz→z0 (z−z0 ) = limz→z0 00 .
Q(z) Q(z0 ) + Q0 (z0 )(z − z0 ) + Q 2!(z0 ) (z − z0 )2 + ...
P (z) P (z0 )
Resz=z0 f (z) = limz→z0 Q00 (z0 )
= ,
Q0 (z0 ) + (z − z0 ) + ... Q0 (z0 )
2!
ϕ(m−1) (z0 )
Resz=z0 f (z) = b1 = .
(m − 1)!
Ejercicios
Determinar los radios y las regiones de convergencia de las siguientes series:
P∞ zn
1. n=1 n Respuesta: R = 1; |z| < 1
P∞ n(z−1)n
2. n=1 2n Respuesta: R = 2; |z − 1| < 2
Serie de Taylor
1
3. Determinar el radio de convergencia de la serie de Taylor para f (z) = z−1 ,
siendo z0 = j el centro de la serie.
√
Respuesta: R = 2
2
4. Hallar la serie de Taylor de f (z) = z−3 en potencias de z − 1. (Ver desa-
rrollo en el video)
Serie de Laurent (Ver desarrollos en el video)
2
5. Hallar la serie de Laurent de f (z) = z−3 en potencias de z − 1 en la zona
lejana
1
6. Hallar todas las series de Laurent de f (z) = − (z−1)(z−2) con centro en
z0 = 0
1
7. Encontrar todas las series de Laurent de f (z) = 1−z 2 con centro en z0 = 1
Residuos. Determinar el residuo de cada una de las siguientes funciones
en sus respectivas singularidades.
CAPÍTULO 1. TEMAS DE VARIABLE COMPLEJA 31
z−2
8. f (z) = z 2 −4 Respuesta: Resz=2 f (z) = 0; Resz=−2 f (z) = 1
2
9. f (z) = e1/z Respuesta: Resz=0 f (z) = 0
e−z
10. f (z) = (z−1)2 Respuesta: Resz=1 f (z) = −e−1
1
11. f (z) = z 4 +1
Respuesta:
e−j3π/4 e−jπ/4
Resz=ejπ/4 f (z) = ; Resz=ej3π/4 f (z) =
4 4
ejπ/4 ej3π/4
Resz=ej5π/4 f (z) = ; Resz=ej7π/4 f (z) =
4 4
4−3z
12. f (z) = z 2 −z Respuesta: Resz=0 f (z) = −4; Resz=1 f (z) = 1
13. f (z) = z 3 sen z12 Respuesta: Resz=0 f (z) = 0
Análisis en el dominio
natural
2.1. Señales
Una señal es la abstracción de una cantidad medible y puede ser representada
por una función de una o más variables independientes. Puede también inter-
pretarse a una señal como una perturbación que experimenta un medio (cambio
en su estado). Lo relevante de una señal es que ese cambio puede desplazarse,
es decir, la señal puede propagarse. Según la naturaleza del cambio es posible
distinguir diferentes dominios de propagación.
Llamamos dominio natural de la señal a aquel en el que el fenómeno fı́sico
medido está ocurriendo. Es decir, podemos tener transmisión de una señal que
dependa del tiempo, del espacio (en cualquier dimensión) o señales que varı́en
tanto en el tiempo como en el espacio.
Una señal tiene como una de sus propiedades relevantes la capacidad de co-
municar y/o transmitir información. La información a que nos referimos
podrı́a ser también la respuesta de un sistema a una solicitación. Esa informa-
ción toma la forma de perturbación y las perturbaciones pueden ser de diferente
ı́ndole. Por ejemplo, las señales de radio son perturbaciones electromagnéticas
que se propagan en el espacio. En este caso se trata de un estado de tensión
eléctrica del medio que se propaga. Las señales acústicas son perturbaciones
de presión que se propagan en un medio material. Las señales telegráficas son
perturbaciones eléctricas que se propagan a través de un conductor etc.
Para fı́sicos e ingenieros las señales tienen además un propósito no menos im-
portante que el anterior. En muchas oportunidades aplicamos señales para es-
tudiar un fenómeno de la naturaleza. En este caso lo relevante es que la señal
es modificada por el fenómeno que se estudia y aquı́ el énfasis no se hace en
evitar que la señal se distorsione sino en que la distorsión de ésta sea originada
por el fenómeno relevante que nos interesa estudiar. Entonces, resulta de gran
importancia tener información acerca del mecanismo de distorsión de la señal
generada.
32
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 33
Se podrı́a afirmar que el mundo moderno está repleto de señales. La mayor parte
de ellas no es posible percibirlas con los receptores de que está dotado el hombre
en forma natural y muchas, aún captadas con ayudas creadas por la inventiva
del ser humano, necesitan ser analizadas para obtener de ellas información útil.
De eso trata el análisis de señales. Esta disciplina es un conjunto de ideas
y recursos que permiten la interpretación de las señales naturales o artificiales
que inundan nuestro universo. Los recursos son de dos tipos, unos que modifican
las señales para su interpretación y otros que las toman y las ponen de forma
que resulten evidentes sus principales caracterı́sticas. En este curso daremos
una introducción a la primera parte de estas técnicas, es decir, aprenderemos a
analizar señales para obtener de ellas la información que nos interesa.
Desde tiempos muy antiguos los seres humanos han empleado señales de dife-
rentes tipos para comunicar acontecimientos importantes o dar voces de alarma.
Como ejemplos podemos citar:
- Haces luminosos que empleaban los griegos y romanos para fines militares.
- Señales de humo tradicionalmente empleadas por indios para enviar mensajes.
- Señales acústicas de tambores para comunicarse en sitios de difı́cil acceso
- Señales producidas por faros para guı́a de barcos y comunicación de noticias.
A fines del siglo dieciséis, Inglaterra empleó un sistema de faros para alertar
sobre la proximidad de la armada española. En esa época se acuñó el término
de “señal” para denotar un signo o noticia perceptible por el oı́do o la vista,
destinada a advertir, transmitir una información o comunicar alguna noticia. El
sistema de semáforos (del griego: “sema” = seña; “phoro” = portador) estaba
en Inglaterra tan perfeccionado en el año 1806, que era posible transmitir una
señal desde Plymouth a Londres obteniendo una confirmación en tan solo tres
minutos. En el año 1852, Morse inventó un código, el que junto a la invención del
telégrafo produjo un gran avance en las comunicaciones, tanto desde el punto
de vista de su rapidez como el de su confiabilidad.
Las señales de radar comienzan a ser aplicadas durante la segunda guerra mun-
dial para detectar aviones y alertar sobre la posibilidad de bombardeo. El sonar,
inventado por Langevin (1917), permitió aplicar señales acústicas a los sistemas
de detección de submarinos. Estas señales atraviesan una porción del espacio
cuyas condiciones cambian de forma azarosa, perturbándolas e impidiendo, a
veces, que estas cumplan con los objetivos para los que fueron generadas. Esta
circunstancia hizo que se desarrollara la teorı́a de señales desde el punto de vista
de su generación, mejorando la electrónica asociada para hacer las señales más
robustas. Desde el punto de vista de la detección, se desarrolló la herramien-
ta estadı́stica para buscar las caracterı́sticas relevantes de la información que
llegaba alterada a los sistemas de detección.
En el área de las comunicaciones, entonces, la misión de la teorı́a de señales
cumple múltiples propósitos: debe mejorar su generación para hacerlas inmunes
a las perturbaciones del medio, facilitar su recuperación e incluso hacer que la
transmisión de éstas sea más económica.
Una señal puede ser también un proyectil cuya trayectoria se altera por la pre-
sencia de un obstáculo, o una modificación en la amplitud y en el contenido
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 34
de frecuencia de una onda que se hace incidir sobre una zona afectada por un
fenómeno determinado, como una turbulencia, por ejemplo. Estas son las ideas
que subyaceb bajo el funcionamiento del radar, del sonar o de la ecografı́a ul-
trasónica.
Para fines de investigación las señales las generamos mediante dispositivos que
transforman el fenómeno de interés en tensiones eléctricas y viceversa. Un siste-
ma que transforma señales de un tipo en otras se llama transductor. Hay tantos
transductores como las interacciones que gobiernan su capacidad de transfor-
mación, ası́ tenemos transductores eléctrodinámicos (un buen ejemplo de estos
son los parlantes de las radios), magnéticos, electromagnéticos, piezoeléctricos,
termoeléctricos, etc.
Las señales, desde el punto de vista temporal se pueden dividir en continuas
y discretas, y ambas pueden ser periódicas o no periódicas. En este curso nos
ocuparemos de señales continuas.
En la próxima sección describimos las señales periódicas y sus principales carac-
terı́sticas, la mayorı́a de las cuales nos son familiares: perı́odo, amplitud, forma
de onda. También consideraremos señales no periódicas, más ligadas a even-
tos no repetitivos. Estas señales no son tan fáciles de describir. Sin embargo,
basándose en la forma en que éstas se despliegan en el tiempo, es posible dar
algunas de sus más importantes caracterı́sticas.
f (x) = f (x + mT ),
con m = 1, 2, 3...
Demostración:
Sea f (x + mT ) = f (x + T ) = f (x) con m = 1, entonces,
si m=2:
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 35
f (x + 2T ) = f ((x + T ) + T ) = f (x + T ) = f (x),
si m=n:
f (x + nT ) = f (x + (n − 1)T + T ) = f (x + (n − 1)T ) = ... = f (x).
2. Si f (x) es periódica y tiene perı́odo T , entonces f (ax) tiene periodo T /a.
Demostración:
Sea g(x) = f (ax).
Es evidente que g(x) es periódica. Entonces, supongamos que τ es el pe-
riodo de g(x).
Luego, g(x) = g(x + τ ) y f (ax) = f (a(x + τ )).
Como f (ax) = f (ax + aτ ), entonces aτ es el perı́odo de f (ax) y, por lo
tanto aτ = T ⇒ τ = Ta .
3. Si f (x) es periódica de perı́odo T1 y g(x) es periódica de perı́odo T2 ,
entonces si existe un valor T = aT1 = bT2 , con a/b un número racional, una
nueva función definida como la suma de las anteriores, y(x) = f (x) + g(x),
será periódica, de perı́odo T .
Esto es equivalente a decir que la relación entre los perı́odos de las señales
sumando debe ser un número racional.
Valores medios
En ciertas oportunidades resulta muy útil describir las señales, periódicas o no,
mediante un número limitado de parámetros que reflejan magnitudes más fáciles
de interpretar desde el punto de vista fı́sico. Estos valores están relacionados con
las propiedades que las funciones tienen en promedio.
Valor medio o valor PN promedio. Magnitud que es muy familiar en su
forma discreta: N1 i=1 fxi , donde fxi es el valor de cada muestra (punto
discreto) de la señal, y N el número de muestras. Cuando se trata de una
señal continua esta expresión se transforma en la bien conocida expresión
para el valor medio de una señal continua f (x) definida en a ≤ x ≤ b:
Z b
1
Vm = f (x)dx.
b−a a
Una señal tiene energı́a finita si E < ∞. En ese caso, caracterizaremos a f (t)
como una señal de energı́a.
Diremos que la señal f (t), definida en (−∞, ∞), es una señal de potencia si:
Z L
1
0 < lı́m |f (t)|2 dt < ∞.
L→∞ 2L −L
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 37
Puede probarse también que toda señal periódica es siempre una señal de po-
tencia, y que toda señal acotada en tiempo y amplitud de las llamadas “pulso”,
es siempre de energı́a.
Hay señales que no pertenecen a ninguna de las familias anteriores, que son las
de potencia infinita. Para estas señales:
Z L
1
lı́m |f (t)|2 dt = ∞.
L→∞ 2L −L
con
1
fp (t) = (f (t) + f (−t)])
2
y
1
fi (t) = (f (t) − f (−t)).
2
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 38
Figura 2.3
Sumando f (t) y f (−t) se obtiene una función par, y restándolas se obtiene una
función impar.
Luego, puede verse que sumando las dos funciones obtenidas se vuelve a tener
Figura 2.4
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 39
la f (t) dada.
Para tales funciones existen los lı́mites por izquierda y por derecha en el intervalo
[a,b], es decir lı́mt→a+ f (t) = f (a) y lı́mt→b− f (t) = f (b). Además, existen los
lı́mites por derecha y por izquierda en cada punto de discontinuidad con salto
finito.
Q[a, b] es el espacio de funciones seccionalmente continuas en [a, b], y Q0 [a, b]
simboliza el espacio de funciones seccionalmente continuas en [a, b], diferencia-
bles con continuidad por tramos en [a, b].
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 40
Escalón unitario
También llamada función de Heaviside.
0 si t<0
us (t) =
1 si t≥0
2.
0 si t < −2
us (t + 2) =
1 si t ≥ −2
3.
0 si t<2
−0, 7us (t − 2) =
−0,7 si t≥2
dr(t)
=0 t<0
dt
y
dr(t)
=1 t ≥ 0.
dt
Por lo tanto, la derivada de la rampa unitaria es el escalón unitario
dr(t)
= us (t).
dt
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 42
Ver Fig.2.12.
Función sinc
sen(t)
si t 6= 0
sinc(t) = t
1 si t=0
La función en t = 0 es el
sen(t)
lı́m = 1,
t→0 t
ya que t = 0 es una singularidad evitable de sinc(t).
La sinc(t) cumple con las siguientes propiedades:
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 44
sen(t) sen(t)
1. lı́mt→∞ t = 0 y lı́mt→−∞ t = 0 pues,
|sen(t)|
lı́mt→∞ | sen(t)
t | = lı́mt→∞ |t| ≤ lı́mt→∞ 1
|t| = 0,
sen(t)
⇒ lı́mt→∞ t =0
De la misma manera se analiza el caso para t → −∞.
2. Los ceros de la sinc(t) ocurren en t = kπ, k ∈ Z − {0}.
sen(t)
sinc(t) = 0 ⇔ = 0 ⇔ sen(t) = 0 ⇔ t = kπ; k ∈ Z.
t
El valor de k = 0 se excluye porque, como se mostró arriba, la sinc(t)
vale 1 en t = 0.
3. La función sinc(t) es par.
sen(−t) −sen(t) sen(t)
sinc(−t) = = = = sinc(t).
−t −t t
siempre que Z ∞
f (t)dt = 1.
−∞
y considerando que
lı́m f (t) = us (t).
a→0
Entonces,
dus (t) d lı́ma→0 f (t)
= .
dt dt
Finalmente, dado que puede intercambiarse el lı́mite con la derivada, y
que
df (t) 1
= Pa (t − a/2)
dt a
podemos escribir
df (t) 1
lı́m = lı́m Pa (t − a/2) = δ(t).
a→0 dt a→0 a
Rt
8. De la propiedad anterior surge que −∞ δ(t)dt = us (t).
Rt
9. t12 f (t)δ (k) (t − t0 )dt = (−1)k f (k) (t0 ), con t1 < t0 < t2 .
10. δ(t) = δ(−t).
11. Toda señal x(t) se puede expresar como un continuo de impulsos ponde-
rados.
Demostración. Según la propiedad 3, x(t)δ(t − t0 ) = x(t0 )δ(t − t0 ), que
si hacemos t0 = τ queda x(t)δ(t − τ ) = x(τ )δ(t − τ ). Integrando ambos
miembros: Z ∞ Z ∞
x(t)δ(t − τ )dτ = x(τ )δ(t − τ )dτ
−∞ −∞
Z ∞ Z ∞
x(t) δ(t − τ )dτ = x(τ )δ(t − τ )dτ,
−∞ −∞
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 48
R∞
y considerando que −∞
δ(t − τ )dτ = 1, queda:
Z ∞
x(t) = x(τ )δ(t − τ )dτ.
−∞
Luego, como se dijo, la señal x(t) quedó expresada como una suma conti-
nua de impulsos ponderados.
Ejercicios
1. Calcular la energı́a de la señal
0 si t<0
f (t) =
−bt
ae si t≥0
a2
Respuesta: E = 2b si b ∈ R+ .
2. Probar que la señal f (t) = e−at con t ∈ (−∞, ∞) no es una señal de
energı́a ni de potencia.
3. Escriba como combinación de escalones y grafique P4 (t), P4 (2t), P4 (2 − t)
y P4 (t/2).
4. Verificar analı́ticamente y gráficamente que la rampa finita se puede ex-
presar como
x(t) = r(t) − r(t − 1).
Las entradas xi (t), i = 1, 2, ...n y las salidas yj (t), j = 1, 2, ...m serán en general
señales de tiempo. Es decir, cualquier variable fı́sica que varı́e con el tiempo.
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 50
En este curso se analizarán sistemas de una entrada y una salida (SISO: Single
Input Single Output), sin perturbaciones externas.
Clasificación de sistemas
De acuerdo a determinadas caracterı́sticas, existen diferentes clasificaciones para
los sistemas.
Con memoria (dinámico) y sin memoria (estático).
Un sistema en el cual la salida en el instante t depende exclusivamente
de la entrada en ese instante es llamado sistema estático o sin memoria.
En contraposición, un sistema dinámico es uno en el cual la salida en el
instante t depende también de valores pasados y/o futuros de la entrada,
además del valor en t. Se lo llama también sistema con memoria.
Ejemplo. El sistema formado por un capacitor C con la entrada definida
como la corriente i(t) que circula por el mismo tiene memoria de duración
infinita, ya que la salida
1 t
Z
v(t) = i(t)dt
C −∞
indica que la salida depende de la entrada en todo el intervalo (−∞, t).
Causales y no causales.
Un sistema es causal si su respuesta a una entrada no depende de valores
futuros de esa entrada y/o valores futuros de salidas. Un sistema que no
verifica esta propiedad es llamado no causal. En particular, un sistema se
dice anticausal si su respuesta a una entrada depende exclusivamente de
valores futuros de esa entrada y/o valores futuros de salidas
De tiempo continuo y de tiempo discreto.
Un sistema es un sistema de tiempo continuo si opera sobre señales defini-
das en un intervalo continuo de tiempo y produce como respuesta señales
también continuas en el tiempo. Por otra parte, si el sistema procesa
señales que están definidas únicamente en instantes particulares de tiempo
(generalmente equiespaciados), es llamado de tiempo discreto. Para indi-
car el valor de estas señales en los instantes de tiempo tk se usa la notación
xk o x(k).
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 51
Lineal y no lineal.
Un sistema es lineal si verifica el principio de superposición, tanto para
entradas como para condiciones iniciales. Este principio se definirá más
adelante. Un sistema que no verifica el principio de superposición es no
lineal.
Invariante y variante.
Un sistema es invariante (estacionario) si su salida es siempre la misma
cada vez que se aplica la misma entrada para las mismas condiciones
iniciales, sin importar el instante en que se aplique la entrada.
Los sistemas invariantes en el tiempo se modelan mediante ecuaciones
diferenciales con coeficientes constantes.
Una caracterización sencilla de los sistemas invariantes en el tiempo se
tiene desplazando en el tiempo a la señal de entrada. El sistema es inva-
riante si x(t) → y(t) y x(t − T ) → y(t − T ). Es decir, que la misma entrada
desplazada genera la misma salida desplazada.
Análisis de sistemas
Una de las clasificaciones de sistemas los separa en lineales y no lineales. Y
aunque un sistema fı́sico nunca es completamente lineal, en ciertas áreas de
aplicación un modelo lineal se puede usar para representar apropiadamente al
sistema. El análisis de sistemas lineales es importante pues existe una gran
cantidad de teorı́as matemáticas que permiten analizarlos.
Por el contrario, los sistemas no lineales requieren análisis particulares, es decir,
cada sistema no lineal debe estudiarse como un caso único. No hay métodos
generales ni soluciones generales.
Como ingenieros, no sólo interesa el análisis, sino también la sı́ntesis o diseño de
sistemas, que constituye la parte creativa de la ingenierı́a. De aquı́ que, como en
otras actividades creativas, para abordar el diseño de sistemas primero se debe
aprender a analizarlo.
El análisis de sistemas puede dividirse en tres fases:
El desarrollo de un modelo matemático apropiado para el problema fı́sico
que se trate. En esta parte se deben obtener las ecuaciones de movimien-
to, condiciones de frontera o iniciales, valores de parámetros, etc. Es la
denominada etapa de modelado.
Después de obtener un modelo apropiado, se resuelven las ecuaciones re-
sultantes para encontrar la solución. Para ello, pueden aplicarse diversos
métodos.
Luego, la solución del modelo matemático se relaciona o interpreta en fun-
ción del problema fı́sico. Es evidente que lo descripto en el primer recuadro
debe ser tan exacto que permita hacer predicciones significativas concer-
nientes al sistema fı́sico. Esta última fase es la denominada validación
del modelo.
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 52
Principio de superposición
Un sistema es lineal si satisface el principio de superposición, que engloba las
propiedades de homogeneidad (escalado) y aditividad.
Homogeneidad
Si para la entrada x(t) se obtiene la salida y(t), y para la entrada αx(t)
se obtiene αy(t), cualquiera sea la constante α, entonces decimos que se
cumple la propiedad de homogeneidad. Es decir:
x(t) → y(t)
αx(t) → αy(t)
Aditividad
Si
x1 (t) → y1 (t)
y
x2 (t) → y2 (t),
entonces
x1 (t) + x2 (t) → y1 (t) + y2 (t).
y(t) = Lx(t),
L(αx1 (t) + βx2 (t)) = αLx1 (t) + βLx2 (t) = αy1 (t) + βy2 (t).
Modelado de SLIT
Cuando la función de entrada x(t) y la función de salida y(t) de un sistema
dinámico están relacionadas por una ecuación diferencial lineal ordinaria a co-
eficientes constantes (EDO), entonces puede probarse que ese sistema es lineal
e invariante en el tiempo (SLIT). Además, como la EDO involucra derivadas, el
sistema es dinámico o con memoria. Para sistemas de una sola entrada-salida,
la ecuación general puede ser:
an y (n) (t) + ... + a1 y 0 (t) + a0 y(t) = bm x(m) (t) + ... + b1 x0 (t) + b0 x(t),
dk y(t)
donde y (k) (t) = dtk .
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 53
donde
B(D)
H(D) =
A(D)
es la función operacional del sistema, un operador que actúa sobre la función
de entrada para producir la salida.
Finalmente, la ecuación diferencial escrita en forma operacional es:
y(t) = H(D)[x(t)].
Ecuación diferencial:
Z t
0 1
v(t) = Li (t) + Ri(t) + i(t)dt
C −∞
1
v(t) = Lq 00 (t) + Rq 0 (t) + q(t)
C
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 54
Ecuación diferencial:
Z t
0 1 1
i(t) = Cv (t) + v(t) + v(t)dt
R L −∞
Ecuación diferencial:
Vibraciones de torsión
Ecuación diferencial:
Ecuación diferencial:
1
P (t) = M V 00 (t) + Ra V 0 (t) + V (t)
Ca
Para cierto sistema dado, no hay un único diagrama en bloques que lo represente.
Veamos que si reescribimos la ecuación diferencial anerior como
1 dy(t)
y(t) = (bx(t) − ),
a dt
A(D)[y(t)] = x(t).
donde los coeficientes son los mismo del operador A(D). Esa ecuación tiene
n raı́ces, que son las que determinan la forma de las funciones yk (t).
1. Raı́ces reales distintas r1 , r2 , ..., rn , entonces
c1 eat cos (bt) + c2 eat sen(bt) + c3 teat cos (bt) + c4 teat sen(bt)
+... + c2m−1 tm−1 eat cos (bt) + c2m tm−1 eat sen(bt)
Solución.
- Solución homogénea. En primer lugar, buscamos las raı́ces de la ecua-
ción caracterı́stica, que en este caso es:
f (r) = r + 1,
y que tiene como única raiz (por ser sistema de primer orden) a r = −1.
La forma de la solución homogénea será:
yh (t) = c1 e−t .
yp (t) = c2
yp (t) = 3.
y(0) = c1 e0 + 3
1 = c1 + 3,
de donde c1 = −2.
Y ası́, vemos que
y(t) = −2e−t + 3.
aún cuando
y(t) = yl (t) + yf (t) = yh (t) + yp (t).
yh (t) = −2e−t ;
yp (t) = 3
y la solución total resultó:
y(t) = −2e−t + 3.
yl (t) = e−t .
Observemos que con este método se aplican las condiciones iniciales a la solución
libre, ya que la forzada tiene condiciones inciales nulas.
- Solución forzada. Esta respuesta se halla mediante la operación de con-
volución que estudiaremos más adelante. Entonces, es esta instancia vamos a
proponer una forma para la solución forzada en base a la entrada que se aplica
(que es una constante) y a las caracterı́sticas propias del sistema, con condiciones
iniciales nulas. Se propone:
yf (t) = k2 + k3 e−t .
Ası́,
yf (t) = 3 − 3e−t .
Ahora, sumando las respuestas libre y forzada se llega a la misma solución total
que sumando la homogénea y la particular:
y(t) = −2e−t + 3.
Estabilidad y causalidad
Sistema causal
Se dice que un sistema es causal o no anticipativo si la salida en el instante
t0 solo depende de la entrada para t ≤ t0 . Informalmente, podrı́amos decir
que la entrada es la causa de la salida. En otras palabras, la salida actual
no depende de valores futuros de la entrada. Intuitivamente, es claro que
todo sistema mecánico o eléctrico de existencia fı́sica en el mundo real es
inherentemente causal.
Sin embargo, actualmente hay una cantidad de situaciones ingenieriles que
requieren la formulación de sistemas no-causales. Por ejemplo, la “inteli-
gencia artificial” se basa justamente en sistemas anticipativos. En muchos
casos, se puede prevenir una salida, o anticipar la decisión. Esto sucede
en edificios inteligentes, seguimiento de satélites, calefacción de ambientes,
etc.
Para pensar: denominamos como ”derivador puro” a un sistema que se
define por la relación: y(t) = dx(t)
dt ¿Será éste un sistema causal?
Sitema estable
Una descripción intuitiva del concepto de estabilidad es la siguiente. Su-
ponga un sistema en reposo, sin entrada ni salida. En un instante, se pro-
porciona una entrada a ese sistema, que podrá ser una energı́a eléctrica,
mecánica, acústica o electromagnética. Si la señal de entrada desaparece
luego de un cierto tiempo o se mantiene en un valor constante, el sistema
será estable siempre que la señal de salida desaparezca o se mantenga en
un valor finito, respectivamente, luego de cierto tiempo. Si, por el contra-
rio, la señal de salida crece y crece cuando la señal de entrada se mantiene
acotada, entonces diremos que el sistema es inestable.
Ası́, un criterio de estabilidad dice que, si la entrada de un sistema es
acotada, la salida es acotada (BIBO: Bounded-Input Bounded-Output).
Entonces, el sistema es estable si:
para |x(t)| < M , ∀t resulta que |y(t)| < K, ∀t, con M y K constantes
reales.
Hemos visto que
y(t) = yh (t) + yp (t),
donde yp (t) tiene la misma forma que la entrada. Ası́, si la entrada está
acotada en amplitud, la solución particular también lo estará. Entonces,
para que la solución total esté acotada la solución homogénea debe hacerse
cero.
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 64
Recordemos que
|ck erk t | = |ck ||e(σk +jωk )t | = |ck ||eσk t ||ejωk )t | = |ck ||eσk t |.
Ejercicios
1. Sea un sistema con relación de entrada-salida dada por la ecuación lineal
y(t) = ax(t) + b. Verificar si el sistema es lineal o no.
2. Verificar si es lineal el sistema dado por la relación Lx(t) = 3x(t)(t + 4).
La integral de convolución
Consideremos un sistema LIT con condiciones iniciales nulas al que se le aplica
un impulso de peso k. La salida será kh(t) porque el sistema es lineal:
kδ(t) → kh(t).
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 66
kδ(t − t0 ) → kh(t − t0 ).
Por otro lado, y como ya vimos antes, cualquier función puede escribirse como
una suma continua de funciones impulso:
Z ∞
x(t) = x(τ )δ(t − τ )dτ.
−∞
La salida del sistema con condiciones iniciales nulas queda expresada por una
integral, dadas la entrada x(t) y la respuesta al impulso h(t), llamada integral
de convolución.
La notación abreviada para convolución es:
|x(t)| < M, ∀t
Interesa analizar la cota de la salida del sistema. Para ello, se aplica el concepto
de valor absoluto a ambos miembros de la igualdad dada por la integral de
convolución, y se obtiene que:
Z ∞ Z ∞ Z ∞
|y(t)| = | x(t − τ )h(τ )dτ | < |x(t − τ )||h(τ )|dτ < M |h(τ )|dτ.
−∞ −∞ −∞
R∞
Puede verse que y(t) está acotada si −∞
|h(τ )|dτ < ∞, es decir, si h(t) es
absolutamente integrable.
Para hallar la respuesta al impulso, h(t), se tendrá que x(t) = δ(t) y que las
condiciones iniciales del sistema son nulas:
donde las condiciones iniciales anteriores se emplean para calcular los coeficien-
tes de h(t).
Prueba. Consideremos el caso de la ecuación A(D)[y(t)] = x(t). Su solución,
cuando x(t) y h(t) son causales, puede expresarse como:
Z t
y(t) = x(τ )h(t − τ )dτ.
0
d t
Z Z t
d
y 0 (t) = f (t, τ )dτ = f (t, τ )|τ =t + f (t, τ )dτ,
dt 0 0 dt
con f (t, τ ) = x(τ )h(t − τ ).
Ası́, queda que:
d t
Z Z t
d
y 0 (t) = x(τ )h(t − τ )dτ = x(t)h(t − τ )|τ =t + (x(τ )h(t − τ ))dτ
dt 0 0 dt
Z t
0
y (t) = x(t)h(0) + x(τ )h0 (t − τ )dτ.
0
Como h(t) debe ser tal que las condiciones iniciales sobre la salida del sistema
sea cero, Z 0
0
y (0) = x(0)h(0) + x(τ )h0 (0 − τ )dτ = 0,
0
se ve que
h(0) = 0.
d t
Z Z t
(n) (n−1) (n−1)
y (t) = x(τ )h (t − τ )dτ = x(t)h (0) + x(τ )h(n) (t − τ )dτ.
dt 0 0
La ecuación diferencial que modela la relación entre ambas señales, para ciertos
valores de los parámetros, es:
Se vio que las condiciones inciales son ĥ(0) = 0; ĥ0 (0) = 1, por lo tanto:
ĥ0 (t) = (−c2 e−t sen(t) + c2 e−t cos(t))us (t) + c2 e−t sen(t)δ(t),
ĥ0 (0) = c2 → c2 = 1
y
ĥ(t) = e−t sen(t)us (t).
h(t) = (−e−t sen(t) + e−t cos(t))us (t) + e−t sen(t)δ(t) + e−t sen(t)us (t).
Por propiedad de la función impulso, el término que contiene δ(t) se hace cero,
y entonces:
h(t) = e−t cos(t)us (t).
Entonces, para este sistema con condiciones iniciales nulas y con una cierta
entrada x(t) causal, la salida y(t) puede expresarse como:
Z t
y(t) = x(τ )e−(t−τ ) cos(t − τ )dτ.
0
A(D)[h(t)] = B(D)[δ(t)],
D2 [h(t)] = (e−t cos (t) + 2e−t sen(t) − e−t cos (t))us (t)+
+(−e−t cos (t) − e−t sen(t))δ(t) + δ 0 (t)
D2 [h(t)] = 2e−t sen(t)us (t) − δ(t) + δ 0 (t).
d t
Z Z t Z t
dg(t) d
= h(τ )us (t−τ )dτ = h(τ ) us (t−τ )dτ = h(τ )δ(t−τ )dτ = h(t).
dt dt 0 0 dt 0
Ejemplo
Sea la respuesta al impulso de cierto sistema h(t) = e−3t us (t). La respuesta del
sistema cuando la entrada es un escalón será:
Z t Z t
−3τ 1 − e−3t
g(t) = e us (τ )dτ = e−3τ dτ = us (t).
0 0 3
Ejercicios
1. Determinar si es causal y/o estable el sistema cuya respuesta al impulso
está dada por h(t) = e−3t us (t) .
2. Sea x(t) = us (t) y h(t) = e−t us (t). Halle y(t) = x(t) ∗ h(t).
Respuesta: y(t) = 1 − e−t , t ≥ 0.
3. Hallar la convolución entre x(t) = us (t) − us (t − 4) y h(t) = e−2t us (t).
0
si t<0
1 1 −2t
Respuesta: y(t) =
2 − 2e si 0 ≤ t ≤ 4
1 −2t 8
2e (e − 1) si t>4
CAPÍTULO 2. ANÁLISIS EN EL DOMINIO NATURAL 74
Análisis de Fourier
Aquı́ se trata de representar una señal de forma que se pueda manejar con
formas potenciales. Una alternativa consiste en suponer que es posible ajustar
un polinomio para que sus valores representen los de la señal. Para esto se debe
obtener un sistema de ecuaciones con las que se determinan los coeficientes del
polinomio. Mientras mayor sea el número de puntos del polinomio que coincidan
con los de la curva que se pretende ajustar, mayor será el grado del polinomio
a considerar. Esto sugiere que el método es bastante potente y que se puede
conseguir tanta precisión como se quiera en el citado ajuste. Sin embargo, que el
sistema propuesto presenta ciertas debilidades para el caso de señales periódicas.
75
CAPÍTULO 3. ANÁLISIS DE FOURIER 76
El sentido del método es bastante claro y para ciertas curvas resulta suficiente,
pero la función resultante de este ajuste no es periódica, por lo tanto, es aplicable
solo a un ciclo de una señal periódica (o parte de él). Por esta razón sólo da
valores suficientemente aproximados para los puntos que se ha tomado para
hacer el cálculo (y su entorno) y la precisión de la representación se va perdiendo
rápidamente a medida que nos alejamos de los puntos de referencia.
b) Serie de potencias
∞
X
f (t) = an (t − c)n ,
n=−∞
con m, n = 1, 2, ...
Puede probarse que los productos escalares de las funciones del conjunto
tomadas de a pares dan cero, por lo tanto podemos asegurar que el sistema
es ortogonal.
Además, el correspondiente sistema ortonormal, B 0 , se puede deducir di-
vidiendo a cada función de B por su correspondiente norma:
( r r )
0 1 2 2
B = √ , cos (nω0 t), sen(nω0 t)
T T T
Sistema exponencial
Sea Q = [− T2 , T2 ], y B = {ejnω0 t }, con ω0 = 2π T . Puede probarse que B es
ortogonal, ya que:
R T /2
hejnω0 t , ejmω0 t i = −T /2 ejnω0 t e−jmω0 t dt = 0 si m 6= n.
Por otro lado, cuando m = n, se tiene la norma al cuadrado de las funcio-
R T /2
nes: hejnω0 t , ejnω0 t i = −T /2 ejnω0 t e−jnω0 t dt = T .
Ası́ que el sistema ortonormal que se obtiene a partir de B es
0 1 jnω0 t
B = √ e
T
.
Rb R b P∞
hf (t), φm (t)i = a
f (t)φm (t)dt = a n=1 kn φn (t)φm (t)dt
Intercambiando el orden de la sumatoria con la integral queda:
Rb Rb Rb
hf (t), φm (t)i = a k1 φ1 φm dt + a k2 φ2 φm dt + ... + a km φm φm dt + ...
Por brevedad, en la expresión anterior se omitió la dependencia de t de las
funciones.
Como la sucesión de funciones de la base es ortogonal, el único término de la
sumatoria que no se hace cero es el que corresponde a n = m:
Rb
hf (t), φm (t)i = km a φm φm dt = km kφm (t)k2 .
Entonces, la fórmula para calcular los coeficientes de una serie de Fourier que
represente a la función f (t) usando un conjunto de funciones ortogonales φn (t)
es:
hf (t), φn (t)i
kn = .
kφn (t)k2
kn = hf (t), φn (t)i.
entonces
lı́m εm = 0.
n→∞
Por lo tanto, si se aproxima una función f (t) por una serie finita de Fourier
Sm (t), dicha aproximación tiene la propiedad de tener el mı́nimo error cuadráti-
co medio.
CAPÍTULO 3. ANÁLISIS DE FOURIER 81
Si la condición: Z b
lı́m |f (t) − Sm (t)|2 dt = 0,
n→∞ a
es satisfecha por cualquier función f en nuestro espacio funcional, diremos que
el sistema ortonormal {φn (t)} es cerrado en el sentido de la convergencia en
media.
Desarrollando el integrando de la ecuación anterior y teniendo presente la defi-
nición de kn , se llega a que:
(Z m
)
b X
2 2 2
lı́m |f (t)| dt − |kn | kφn k = 0,
n→∞ a n=1
Aproximación cuadrática
Las series de Fourier, como se dijo, desempeñan un papel fundamental en la
aproximación de funciones por medio de funciones más simples. Esta área es
conocida como teorı́a de aproximaciones. A continuación se demostrará que los
coeficientes de la serie de Fourier hallados en el apartado anterior son aquellos
que proporcionan la aproximación con menor error cuadrático medio.
Sea f (t) una función que puede representarse por una serie de Fourier en el espa-
cio de funciones seccionalmente continuas en [a, b], y sean φ1 (t), φ2 (t), ..., φm (t),
m funciones de una sucesión ortonormal {φn (t)}, (n = 1, 2, ...) en dicho inter-
valo, y Km es una combinación lineal de las mismas:
Desde luego, debe definirse primero qué se entiende por el error E de esta
aproximación. Elegimos una definición que mida la bondad de la concordancia
entre f (t) y Km (t) en el intervalo, considerando el caso particular en que {φn (t)}
y f (t) son funciones reales solo a los efectos de simplificar la demostración.
Se elige:
Z b
E= [f (t) − Km (t)]2 dt,
a
que puede considerarse una medida de error, y se desea que sea lo más pequeña
posible. Esta es una aproximación a f (t) por mı́nimos cuadrados. Obsérvese que
E es el cuadrado de la distancia generalizada kf − Km k entre las funciones f (t)
y Km (t).
El error puede expresarse entonces como:
Z b
E= [f (t) − γ1 φ1 (t) − γ2 φ2 (t) − ... − γm φm (t)]2 dt,
a
Sean kn = hf (t), φn (t)i los coeficientes de Fourier de f (t), con {φn (t)} ortonor-
mal, entonces, reemplazando en el segundo término del integrando anterior a
f (t) por su representación en serie y considerendo todos los productos internos
que se anulan, la expresión queda simplemente:
Z b
E= {[f (t)]2 dt − 2k1 γ1 − 2k2 γ2 − ... − 2km γm φm + γ12 + γ22 + ... + γm
2
.
a
Entonces, podemos asegurar que los coeficientes de Fourier de una función f (t)
con respecto a las funciones φ1 (t), φ2 (t), ....., φm (t) de un sistema ortonormal
son aquellos para los cuales una combinación lineal de las m funciones resulta
ser la aproximación cuadrática óptima a f (t) en el intervalo.
El error cuadrático total de la aproximación, respecto de f (t) en el intervalo
fundamental es mı́nimo si y sólo si los coeficientes de dicha aproximación son
los coeficientes de Fourier de f (t):
hf (t), φn (t)i
kn = .
kφn (t)k2
CAPÍTULO 3. ANÁLISIS DE FOURIER 83
y
∞
X
f (t) = cn ejnω0 t .
n=−∞
o
fk (t) = c−k e−jkω0 t + ck ejkω0 t .
Una señal periódica que no cumple esta condición, y por lo tanto no puede
representarse con una serie de Fourier, es la que se muestra en la siguiente
figura:
1
f (t) = ; 0 < t ≤ 1; f (t) = f (t + 1)∀t.
t
T /2
hf (t), 1i
Z
a0 1
= = f (t)dt
2 k1k2 T −T /2
CAPÍTULO 3. ANÁLISIS DE FOURIER 86
T /2
hf (t), sen(nω0 t)i
Z
2
bn = = f (t)sen(nω0 t)dt, n∈N
ksen(nω0 t)k2 T −T /2
T /2
hf (t), ejnω0 t i
Z
1
cn = = f (t)e−jnω0 t dt, n ∈ Z.
kejnω0 t k2 T −T /2
Entonces:
1 R T /2 a0
- Para n = 0, c0 = f (t)dt =
T −T /2 2
a0
c0 =
2
- Para n ∈ N,
1 T /2 1 T /2 1 T /2
Z Z Z
cn = f (t)e−jnω0 t dt = f (t) cos (nω0 t)dt−j f (t)sen(nω0 t)dt
T −T /2 T −T /2 T −T /2
CAPÍTULO 3. ANÁLISIS DE FOURIER 88
an bn an − jbn
cn = −j = .
2 2 2
- De manera análoga puede probarse que
an bn an + jbn
c−n = +j = .
2 2 2
1 T /2
Z
a0 1
= δ(t)dt = ,
2 T −T /2 T
Z T /2
2 2
an = δ(t) cos (nω0 t)dt =
T −T /2 T
y
Z T /2
2
bn = δ(t)sen(nω0 t)dt = 0.
T −T /2
CAPÍTULO 3. ANÁLISIS DE FOURIER 89
Ejercicios
1. Hallar la serie de Fourier correspondiente a la función
0 si −π < x ≤ 0
f (t) =
t si 0 < x ≤ π
f (t) = f (t + T ); T = 2d.
a) Graficar la función.
b) Hallar ambas series de Fourier de f (t).
c) Escribir ambas expresiones de la identidad de Parseval.
3. Verificar todas las expresiones del apéndice C.
y(t) = H(jω)ejωt .
Para hallar la respuesta y(t) de un SLIT a una entrada periódica que puede
representarse por una serie de Fourier,
∞
X
x(t) = cn ejnω0 t ,
n=−∞
Ejemplo.
f (t) = et = cosh (t) + senh(t), siendo cosh (t) una función par y senh(t) una
función impar.
et + e−t et − e−t
f (t) = et = cosh (t) + senh(t) = + = et .
2 2
Observar que, para f (t) y g(t) funciones definidas en [−T /2, T /2] se cumple que:
Si ambas son pares, entonces su producto es par.
Si ambas son impares, entonces su producto es par.
Si una es par y la otra es impar, entonces su producto es impar.
Z T /2 Z T /2 Z T /2
1 −jnω0 t 1 1
cn = f (t)e dt = f (t) cos (nω0 t)dt−j f (t)sen(nω0 t)dt
T −T /2 T −T /2 T −T /2
y entonces
Z T /2
1
cn = f (t) cos (nω0 t)dt ∈ R ∀n ∈ Z
T −T /2
CAPÍTULO 3. ANÁLISIS DE FOURIER 92
y entonces
Z T /2
1
cn = −j f (t)sen(nω0 t)dt es imaginario puro ∀n ∈ Z
T −T /2
Z T /2 Z T /2
2 4
an = f (t) cos (nω0 t)dt = f (t) cos (nω0 t)dt
T −T /2 T 0
bn = 0
La serie trigonométrica es entonces:
∞
a0 X
f (t) = + an cos (nω0 t).
2 n=1
Aplicación
El uso más obvio de las simetrı́as de las ondas es la simplificación en los cálculos
de los coeficientes. Pero además, nos permite hallar de manera eficiente desa-
rrollos en serie de Fourier para señales no periódicas definidas entre [0, t0 ].
Una señal cualquiera, definida en un intervalo, puede expresarse en ese intervalo
mediante diferentes series, según se proponga lo que llamaremos una extensión
par de la misma o una extensión impar.
CAPÍTULO 3. ANÁLISIS DE FOURIER 93
Figura 3.7: Señal definida en [0, T /2] y sus extensiones par e impar.
La serie trigonométrica de Fourier que representa a las extensiones son las ha-
lladas antes para señales con alguna simetrı́a. Pero estas series son las corres-
pondientes a las funciones periódicas fˆ(t). Para representar a la señal original
(no periódica) debe hacerse cero todo lo que está fuera del intervalo [0, T /2]:
∞
!
a0 X
f (t) = + an cos (nω0 t) PT /2 (t − T /4) extensión par
2 n=1
o
∞
!
X
f (t) = bn sen(nω0 t) PT /2 (t − T /4) extensión impar
n=1
CAPÍTULO 3. ANÁLISIS DE FOURIER 94
Puede verse en las expresiones anteriores que cada función armónica está afec-
tado por el respectivo coeficientes y, justamente, esos coeficientes son los que
componen el espectro discreto de frecuencias de la señal que, para una función
f (t) los denominaremos como Fn (nω0 ):
Z T /2
1
F (nω0 ) = f (t)e−jnω0 t dt.
T −T /2
Entonces tendremos
espectro discreto de amplitud: |F (nω0 )| y
espectro discreto de fase: θ(nω0 ).
Ejemplo 1. Se desea encontrar el espectro discreto de frecuencias y graficar el
espectro de amplitud y el de fase de la señal definida como:
0 si −π < t ≤ 0
x(t) =
t si 0<t≤π
Operando, se obtiene:
cos (nπ) − 1 cos (nπ)
X(nω0 ) = +j n 6= 0,
2n2 π 2n
que es
1
X(nω0 ) = j n par
2n
−1 1
X(nω0 ) = −j n impar.
n2 π 2n
Además, para n = 0 tenemos el valor medio de la señal:
Z π
1 π
X(0) = tdt =
2π 0 4
Calculando la fase y la amplitud de X(nω0 ) se obtienen los espectros discretos
que se muestran debajo:
Puede verse que el espectro de una señal periódica rectangular está dada por
una sinc discreta. Esto significa que los valores del espectro, a diferencia del
ejemplo anterior, son números reales positivos o negativos. Estudiemos algunos
casos particulares de valores de d y T .
d = 1/10; T = 1/2; ω0 = 4π ⇒ d/T = 1/5
En este caso el espectro es:
A nπ
F (nω0 ) = sinc ∀n ∈ Z.
5 5
El valor medio de la señal es F (0) = A/5.
La otra caracterı́stica importante de este espectro son los valores de fre-
cuencias en los que se hace cero. Los ceros de esta función se dan para:
nπ
= kπ ∀k ∈ Z − {0},
5
n = 5k ∀k ∈ Z − {0},
Esto significa que entre cada cero del espectro hay 5ω0 . En este caso los
ceros están en ω = nω0 = 5kω0 = 20kπ.
Como el espectro discreto está formado por números reales, es posible
graficarlo en un único gráfico como el que se muestra en la Fig. 3.12.
Figura 3.12: Espectro de una señal rectangular con relación d/T = 1/5 y fre-
cuencia fundamental ω0 = 4π.
CAPÍTULO 3. ANÁLISIS DE FOURIER 98
También puede observarse en el gráfico que el espacio entre las barras del
espectro es ω0 = 4π.
Aunque el espectro esté formado por números reales, también pueden ha-
cerse los gráficos de los espectros de amplitud y de fase. El espectro de
amplitud es sencillamente (Fig.3.13).
A nπ
|F (nω0 )| = sinc ∀n ∈ Z.
5 5
Esto significa que entre cada cero del espectro hay ahora 10ω0 . Y como
ω0 = 2π, los ceros caen nuevamente en ω = nω0 = 10kω0 = 20kπ.
En el siguiente gráfico se observa cómo los ceros están en el mismo lugar,
pero las barras están más cercanas entre sı́, ya que la frecuencia de esta
señal es menor:
Figura 3.15: Espectro de una señal rectangular con relación d/T = 1/10 y
frecuencia fundamental ω0 = 2π.
F (−nω0 ) = F (nω0 )
θ(−nω0 ) = −θ(nω0 ).
(jnω0 )k F (nω0 ), n ∈ Z.
Prueba
Si
∞
X
f (t) = F (nω0 )ejnω0 t ,
n=−∞
entonces
∞
X
0
f (t) = jnω0 F (nω0 )ejnω0 t ,
n=−∞
CAPÍTULO 3. ANÁLISIS DE FOURIER 101
5. Desplazamiento en frecuencia.
Sea f (t) una función periódica de perı́odo T y a ∈ R, entonces la función
g(t) = ejaω0 t f (t)
es periódica de perı́odo T y su espectro de frecuencia discreta es:
G(nω0 ) = F (nω0 − aω0 ) = F ((n − a)ω0 ), n ∈ Z.
6. Convolución en frecuencia.
Si f (t) y g(t) son funciones periódicas de perı́odo T y h(t) = f (t)g(t) es
una función de perı́odo T > 0, entonces el espectro de frecuencia discreta
de h(t) es:
∞
X
H(nω0 ) = F (nω0 ) ∗ G(nω0 ) = F (kω0 )G((n − k)ω0 ), n ∈ Z.
k=−∞
7. Convolución en el tiempo.
Si f (t) y g(t) son funciones periódicas de perı́odo T , pertenecientes al
Q[−T /2, T /2], entonces la función
1 T /2
Z
h(t) = f (τ )g(t − τ )dτ
T −T /2
es continua en (−T /2, T /2), es periódica de perı́odo T y su espectro de
frecuencias discreta es:
H(nω0 ) = F (nω0 )G(nω0 ).
CAPÍTULO 3. ANÁLISIS DE FOURIER 102
∞
X ∞
X
|F− n|2 + |F0 |2 + |Fn |2 .
n=1 n=1
Luego, recordando que si f (t) es real, |F− n| = |Fn |, entonces la potencia media
de la señal puede expresarse en el dominio de la frecuencia como:
Z T /2 ∞
1 X
|f (t)|2 dt = |F0 |2 + 2 |Fn |2 .
T −T /2 n=1
Esta ecuación indica que la potencia media total de f (t) se puede calcular su-
mando todas potencias medias asociadas a cada componente de frecuencia nω0
de f (t).
El conjunto de las potencias de las componentes, |Fn |2 , como función de nω0 se
llama
espectro de potencia de f (t). O sea, el espectro de potencia de f (t)
es |Fn |2 n ∈ Z .
Al igual que con el espectro de frecuencia discreta de una señal de perı́odo T ,
aquı́ también tenemos el espectro de lı́nea de potencia media de f (t) y la curva
de ecuación η = |F (ω)|2 es la envolvente del espectro de potencia media.
Ejercicio
Para una señal rectangular como la del Ejemplo 2, con d = 1/40 y T = 1/8, se
pide:
1. Hallar analı́tica y gráficamente el espectro de potencia media, e indicar el
correspondiente espectro discreto y su envolvente.
CAPÍTULO 3. ANÁLISIS DE FOURIER 103
Figura 3.16
Esta última expresión puede interpretarse como una suma continua de expo-
nenciales complejas, ejωt . Y las funciones f (t) y F (ω) forman lo que se llama
un par transformado, que se representa mediante la notación
f (t) ←→ F (ω).
Demostración
Z ∞ Z ∞Z ∞
1 1
F −1 {F{f (t)}} = F (ω)ejωt dω = f (x)e−jωx dxejωt dω =
2π
−∞ 2π −∞ −∞
Z ∞Z ∞ Z ∞ Z ∞
1 1
= f (x)e−jω(x−t) dxdω = f (x) e−jω(x−t) dωdx =
2π −∞ −∞ 2π −∞ −∞
Z ∞
1 e−jΩ(x−t) − ejΩ(x−t)
= f (x) lı́m dx =
2π −∞ Ω→∞ −j(x − t)
Z ∞
1
f (x) lı́m 2Ωsinc (Ω(x − t)) dx.
2π −∞ Ω→∞
y entonces
Z ∞
1
F −1 {F{f (t)}} = f (x)2πδ(x − t)dx = f (t).
2π −∞
La función pulso rectangular que se muestra tiene una ecuación dada por:
1 si |t| ≤ a/2
f (t) = Pa (t) =
0 otros t
Solución
1) Su transformada, o espectro continuo de frecuencia F (ω), está dada por
la transformada de Fourier de Pa (t).
Z a/2
e−jωa/2 − ejωa/2 2 aω
F (ω) = e−jωt dt = = sen ,
−a/2 −jω ω 2
que puede escribirse como
aω
F (ω) = a sinc .
2
Obsérvese que por ser el espectro una función real de ω puede verse en
un único gráfico. Si fuese una función compleja tendrı́a que mostrarse el
CAPÍTULO 3. ANÁLISIS DE FOURIER 107
y el espectro de fase:
0 si F (ω) > 0
θ(ω) = π si F (ω) < 0, ω > 0
−π si F (ω) < 0, ω < 0
Puede probarse que la integral de esa sinc da 2π/a, con lo que se muestra
que
Pa (t = 0) = 1.
A
F (ω) = .
a + jω
CAPÍTULO 3. ANÁLISIS DE FOURIER 109
Ejercicio
1. Grafique la señal en el dominio temporal.
2. Demuestre que el espectro continuo de frecuencias es Absinc2 (ωb/2).
3. Grafique el espectro encontrado y analice algunas de sus caracterı́sti-
cas.
R∞
X(ω) = − −∞
f (t)sen(ωt)dt es una función impar.
Demostración
R∞ R∞
X(−ω) = − −∞ f (t)sen(−ωt)dt = −∞ f (t)sen(ωt)dt = −X(ω), por ser
el seno una función impar.
3. F (−ω) = F (ω)
Demostración
A partir de lo que se demostró en 2.:
1 ∞ 1 ∞
Z Z
f (t) = R(ω) cos (ωt)dω − X(ω)sen(ωt)dω.
π 0 π 0
y entonces Z ∞
F (ω) = 2 f (t) cos (ωt)dt, ∈ R.
0
Y, de acuerdo a lo visto en 5.:
Z ∞
1
f (t) = R(ω) cos (ωt)dω.
π 0
y entonces Z ∞
F (ω) = −2j f (t)sen(ωt)dt, ∈ Imag.
0
Y, de acuerdo a lo visto en 5.:
Z ∞
1
f (t) = − X(ω)sen(ωt)dω.
π 0
f1 (t) ←→ F1 (ω)
y
f2 (t) ←→ F2 (ω),
entonces
af1 (t) + bf2 (t) ←→ aF1 (ω) + bF1 (ω).
Prueba
En Z ∞
1
f (t) = F (ω)ejωt dω
2π −∞
3. Escala
f (t) ←→ F (ω)
1 ω
f (at) ←→ F
|a| a
Prueba:
Z ∞ Z ∞
−jωt dτ
F{f (at)} = f (at)e dt = f (τ )e−jωτ /a ,
−∞ −∞ a
Demostración
Z ∞ Z ∞ Z ∞
Y (ω) = y(t)e−jωt dt = x(τ )h(t − τ )dτ e−jωt dt
−∞ −∞ −∞
Z ∞ Z ∞ Z ∞ Z ∞
= x(τ ) h(t − τ )e−jωt dtdτ = x(τ ) h(λ)e−jω(λ+τ ) dλdτ
−∞ −∞ −∞ −∞
Z ∞ Z ∞ Z ∞
−jωλ −jωτ
= x(τ ) h(λ)e dλe dτ = x(τ )H(ω)e−jωτ dτ =
−∞ −∞ −∞
Z ∞
= H(ω) x(τ )e−jωτ dτ = H(ω)X(ω).
−∞
5. Convolución en frecuencia
Si x(t) ←→ X(ω) y g(t) ←→ G(ω), entonces
1
x(t)G(t) ←→ X(ω) ∗ G(ω).
2π
f (t) ←→ F (ω)
f (t − t0 ) ←→ F (ω)e−jt0 ω
Prueba
Z ∞
F{f (t − t0 )} = f (t − t0 )e−jωt dt,
−∞
f (t)ejω0 t ←→ F (ω − ω0 )
Prueba Z ∞
F{f (t)e jω0 t
}= f (t)ejω0 t e−jωt dt =
−∞
Z ∞
= f (t)e−j(ω−ω0 )t dt = F (ω − ω0 ).
−∞
8. Diferenciación en el tiempo
Si
f (t) ←→ F (ω)
dn f (t)
←→ (jω)n F (ω)
dtn
Demostración: partiendo de la fórmula de inversión se deriva con respecto a
t ambos miembros de la ecuación. El procedimiento se repite y la propiedad
se demuestra por inducción.
9. Diferenciación en frecuencia
Si
f (t) ←→ F (ω)
dn F (ω)
(−jt)n f (t) ←→
dω n
Demostración: es similar a la de la propiedad 8., pero se parte de la fórmula
de la transformada y se llega, también por inducción, derivando ambos
miembros con respecto a ω.
10. Integración
Si f (t) y g(t) son dos funciones que tienen transformada de Fourier y,
además, Z t
g(t) = f (t)dt,
−∞
F (ω)
G(ω) = + πF (0)δ(ω)
jω
Como f (t) es real, ya se vió que F (−ω) = F (ω). Luego F (ω)F (ω) = |F (ω)|2 y
se tiene la relación de Parseval para señales de energı́a:
Z ∞ Z ∞
1
E= f 2 (t)dt = |F (ω)|2 dω.
−∞ 2π −∞
|F (ω)|2
.
2π
Esta es una función par y real de ω, por lo cual puede escribirse la integral dos
veces entre cero e infinito, con lo que resulta:
Z ∞ Z ∞
1 2
E= |F (ω)| dω = S(ω)dω.
π 0 0
|F (ω)|2
S(ω) = ω ≥ 0,
π
que permite calcular la cantidad de energı́a comprendida en un cierto rango de
frecuencias.
Cuando se analizaron las funciones periódicas se asoció cierta cantidad de po-
tencia con cada armónica. En el caso de las señales de energı́a, se asocia la
CAPÍTULO 3. ANÁLISIS DE FOURIER 118
|F (ω)|2 1 16
S(ω) = = .
π π (4 + ω 2 )2
Supongamos que queremos saber cuál es la energı́a que aportan las frecuencias
comprendidas entre ω1 = 2 rad/s y ω2 = 3 rad/s. Entonces, debemos resolver
la integral: Z 3
1 16
Eb = 2 2
dω.
2 π (4 + ω )
2
En el gráfico siguiente se muestran la señal y su transformada, la función |F (ω)|
2π
y la densidad espectral de energı́a. También se señala la banda de frecuencias
cuya energı́a se desea calcular, energı́a que está dada por el área bajo la función
comprendida entre ω1 = 2 rad/s y ω2 = 3 rad/s.
Y1 (ω) = Y (ω)e−jt0 ω
y1 (t) = y(t − t0 ),
es decir, una entrada desplazada genera la misma salida que la señal sin
desplazar, pero desplazada.
4. Relación entre la densidad espectral de energı́a de la entrada y la salida
de un sistema. La densidad espectral de energı́a de la salida del sistema
es:
|Y (ω)|2
Sy (ω) = , ω ≥ 0.
π
Si la ponemos en función de la transformada de la entrada tendremos:
|H(ω)X(ω)|2
Sy (ω) = = |H(ω)|2 Sx (ω), ω ≥ 0.
π
f (t) = δ(t − t0 )
δ(t − t0 ) ←→ e−jωt0 ,
δ(t) ←→ 1.
* En el dominio de la frecuencia
F (ω) = δ(ω − ω0 ).
CAPÍTULO 3. ANÁLISIS DE FOURIER 121
La transformada inversa es
Z ∞
1 1 jω0 t
f (t) = δ(ω − ω0 )ejωt dt = e .
2π −∞ 2π
ejω0 t ←→ 2πδ(ω − ω0 ),
1 ←→ 2πδ(ω).
Como en este caso la función tiene valor medio cero, puede usarse la di-
ferenciación en t para calcular la transformada. La derivada de la función
signo es
df (t)
= 2δ(t) = g(t).
dt
Luego, por la propiedad de derivación sabemos que
df (t)
F = jωF (ω) = G(ω),
dt
que en este caso es
jωF (ω) = G(ω) = 2.
Por lo tanto, la transformada de la función signo es
2
F (ω) = .
jω
Transformada de Laplace
4.1. Introducción
La transformada de Laplace es una de las transformadas más importantes en
el análisis de sistemas lineales. Se la utiliza para descomponer una señal en la
suma continua de funciones exponenciales de la forma est , donde s = σ + jω es
una frecuencia compleja. La misma ofrece ventajas interesantes.
El análisis de la respuesta a entradas arbitrarias de sistemas fı́sicos que pue-
dan considerarse lineales e invariantes en tiempo continuo (SLIT) involucra la
solución de ecuaciones diferenciales ordinarias o, si se conoce la respuesta al
impulso del sistema, la evaluación de la integral de convolución entre la en-
trada y la respuesta al impulso. Estas técnicas de análisis suelen implicar la
realización de tediosas operaciones matemáticas aún para el caso de señales de
entrada relativamente simples. El uso de la Transformada de Laplace simplifica
considerablemente la solución del problema ya que convierte las ecuaciones dife-
renciales ordinarias en ecuaciones algebraicas en las transformadas, y la integral
de convolución en un producto de transformadas.
Otra ventaja del uso de la Transformada de Laplace para el análisis de SLIT
es que posibilita la obtención de la respuesta total (respuesta libre + respuesta
forzada), y que permite incluir automáticamente las condiciones iniciales en la
solución del problema.
Esta transformada, como la de Fourier, también ofrece transformada inversa:
f (t) ←→ F (s)
Z ∞
F (s) = L{f (t)} = f (t)e−st dt
0
Z C+j∞
1
f (t) = L−1 {F (s)} = F (s)est ds
2πj C−j∞
125
CAPÍTULO 4. TRANSFORMADA DE LAPLACE 126
lı́m |f (t)e−σt | = 0.
t→∞
RC : {s ∈ C/Re(s) > σc , σc ∈ R}
∞ ∞ L
e−(s−a)t
Z Z
F (s) = eat e−st dt = e−(s−a)t dt = lı́m .
0 0 L→∞ −(s − a) 0
Si se escribe la exponencial separando su esponente real de su exponente
imaginaro, se tiene:
L
e−(σ−a)t e−jωt e−(σ−a)L e−jωL 1
F (s) = lı́m = lı́m − ,
L→∞ −(s − a) 0
L→∞ −(s − a) −(s − a)
y se ve que la integral converge, es decir, la transformada de la Laplace
existe si σ > a (y entonces la exponencial e−(σ−a) L se hace cero cuando
L → ∞).
CAPÍTULO 4. TRANSFORMADA DE LAPLACE 128
Entonces,
1
L{e−at us (t)} = ,
s−a
converge para Re(s) > a.
Analice como se modificarı́a lo anterior si a fuese imaginario y si fuese
complejo.
2. Función escalón unitario
f (t) = us (t).
Z ∞ L L
e−st e−σt e−jωt e−σL e−jωL 1
F (s) = e−st dt = lı́m = lı́m = lı́m + .
0 L→∞ −s L→∞ −s L→∞ −s s
0 0
4.3. Propiedades
1. Linealidad
Si f1 (t) y f2 (t) son funciones que admiten transformada de Laplace, y a
y b son constantes reales o complejas, entonces resulta que:
Ejemplo
Hallar la transformada de Laplace de f (t) = cosh (at)us (t).
El coseno hiperbólico puede escribirse como
eat + e−at
cosh (at) =
2
y entonces, la transformada del coseno hiperbólico puede calcularse como
1 1 1 1 1 s
F (s) = L{eat } + L{e−at } = + = 2 ,
2 2 2 s−a s+a s − a2
2. Cambio de escala
1 s
L{f (at)} =
F ,
a a
RC: Re(s) > aσf , con a > 0 y σf la abscisa de convergencia de F (s).
Demostración Z ∞
L{f (at)} = f (at)e−st dt.
0
Haciendo la sustitución τ = at, queda
1 ∞
Z
s 1 s
L{f (at)} = f (τ )e− a τ dτ = F .
a 0 a a
3. Traslación en t
Sea t0 una constante real positiva,
Demostración
Z ∞
L{f (t − t0 )us (t − t0 )} = f (t − t0 )us (t − t0 )e−st dt,
0
5. Convolución en t
Sean f (t) y g(t) funciones causales. Se sabe que la convolución entre ellas
Rt
es f (t) ∗ g(t) = 0 f (τ )g(t − τ )dτ .
CAPÍTULO 4. TRANSFORMADA DE LAPLACE 131
6. Diferenciación en t
L{f (n) (t)} = sn F (s) − sn−1 f (0) − sn−2 f 0 (0) − ... − f (n−1) (0).
Prueba
Para la derivada de primer orden,
Z ∞
df (t) −st
L{f 0 (t)} = e dt,
0 dt
que se resuelve por partes con u = e−st y dv = dfdt(t) dt, con du = −se−st dt
y v = f (t). Entonces, queda:
Z ∞
∞
L{f 0 (t)} = f (t)e−st 0 + s f (t)e−st dt
0
Z t
F (s)
L f (x)dx =
0 s
A pesar de estar estudiando la transformada unilateral de Laplace, que se
aplica a funciones causales, también podemos mostrar que
Z t
F (s) kf
L f (x)dx = + ,
−∞ s s
Comprobación
a) Por definición
Z t Z ∞ Z t
L f (x)dx = f (x)dxe−st dt.
0 0 0
Rt
Se resuelve por parte, con u = 0 f (x)dx y dv = e−st dt, y entonces
−st
du = f (t)dt y v = e−s . Por lo tanto, queda:
t t ∞ ∞
e−st
Z Z Z
1
L f (x)dx =− f (x)dx + f (t)e−st dt =
0 s 0 0 s 0
!
Z 0 Z L
1 0 −Ls F (s)
= e f (x)dx − lı́m e f (x)dx + .
s 0 L→∞ 0 s
El primer término de la expresión anterior es cero porque, para
Re(s) > 0, la exponencial es cero en el lı́mite superior, ya que
Rt
0
f (x)dx puede probarse que es también una función de orden ex-
ponencial. También en el lı́mite inferior la integral es cero. Entonces
Z t
1 F (s) F (s)
L f (x)dx = (0 − 0) + =
0 s s s
Rt
b) Cuando se quiere encontrar la transformada de −∞
f (x)dx se tiene:
Z t Z 0 Z t
L f (x)dx =L f (x)dx + L f (x)dx .
−∞ −∞ 0
R0
Como por definición kf = −∞
f (x)dx, entonces
Z t
F (s)
L f (x)dx = L {kf } + .
−∞ s
dk F (s)
L (−t)k f (t) =
,
dsk
con RC igual a la de F (s)
Analiticidad de F (s)
Si e−σt f (t) es absolutamente integrable, F (s) = L{f (t)}, con t ≥ 0, es analı́tica
en la región del plano Re(s) > σ y lı́ms→∞ F (s) = 0 para Re(s) ≥ σ.
F (s) es analı́tica en su región de convergencia.
Se tendrá entonces:
∞ ∞
( )
X X
F (s) = L{f (t)} = L fn (t) = L {fn (t)}
n=1 n=1
∞
X
F (s) = Fn (s).
n=1
Además,
Fn (s) = L{fn (t)} = L{f1 (t − (n − 1)T )} = F1 (s)e−(n−1)T s .
Entonces:
∞
X ∞
X ∞
X
F (s) = F1 (s)e−(n−1)T s = F1 (s)e−nT s = F1 (s) e−nT s .
n=1 n=0 n=0
Demostración
La propiedad de diferenciación en el tiempo establece que
Z ∞
L{f 0 (t)} = f 0 (t)e−st dt = sF (s) − f (0).
0
Demostración
Si f (t) es discontinua en el origen, entonces f 0 (t) contiene un impulso de peso
f (0+ ) − f (0− ). Es decir, en el origen f 0 (t) = [f (0+ ) − f (0− )]δ(t). La pregunta
que surge es, entonces, cómo considerar el lı́mite inferior en la integral, si en
t = 0+ o en t = 0− .
Tomando t = 0− como lı́mite inferior de integración:
Z ∞
lı́m f 0 (t)e−st dt = lı́m (sF (s) − f (0− ))
s→∞ 0− s→∞
Z 0+ Z ∞
lı́m [f (0+ ) − f (0− )]δ(t)e−st dt + lı́m f 0 (t)e−st dt =
s→∞ 0− s→∞ 0+
de donde sale
lı́m sF (s) = f (0+ )
s→∞
Demostración
Se parte nuevamente de transformada de f 0 (t)
Z ∞
L{f 0 (t)} = f 0 (t)e−st dt = sF (s) − f (0)
0
I n
1 st
X
f (t) = F (s)e ds = Ress=sk F (s)esk t ,
2πj C0 k=1
st
con sk los polos de F (s)e .
Raramente usaremos esta expresión para calcular en la práctica una trans-
formada inversa. Generalmente, el procedimiento que se utilizará consistirá
en manipular la expresión de la transformada dada hasta obtener formas
de transformadas que aparecen en la tabla.
CAPÍTULO 4. TRANSFORMADA DE LAPLACE 138
b0 + b1 s + b2 s2 + ... + bm sm
F (s) = .
a0 + a1 s + a2 s2 + ... + an sn
b0 + b1 s + b2 s2 + ... + bm sm c1 c2 cn
F (s) = = + +...+ ,
an (s − p1 )(s − p2 )...(s − pn ) s − p1 s − p2 s − pn
ck = (s − pk ) F (s)|s=pk ,
con k = 1, 2, ..., n.
Polos reales múltiples
El desarrollo debe incluir, para cada polo múltiple, tantos términos
como la multiplicidad del polo. Por ejemplo, si el polo pi tiene mul-
tiplicidad r, se incluyen por ese polo los términos:
c1 c2 cr
+ + ... + .
s − pi (s − pi )2 (s − pi )r
Polos complejos
Las ecuaciones anteriores son válidas también para polos comple-
jos. Sin embargo, cuando se tienen polos complejos conjugados como
pc1,2 = α±jβ, por cada par complejo conjugado se tendrá un término
de la forma:
c1 s + c2
(s − α)2 + β 2 .
Las constantes del numerador se resuelven, generalmente, planteando
un sistema de ecuaciones de dos incógnitas.
CAPÍTULO 4. TRANSFORMADA DE LAPLACE 139
n
X P (pk ) −pk t
f (t) = e ,
Q0 (pk )
k=1
1 − e−5s
F (s) = .
s(1 − e−10s )
an (sn Y (s)−sn−1 y(0)−...−y (n−1) (0))+an−1 (sn−1 Y (s)−sn−2 y(0)−...−y (n−2) (0))+...
... + a1 (sY (s) − y(0)) + a0 Y (s) = X(s).
Agrupando por un lado los términos del primer miembro que involucran a la
función Y (s) y, por otro lado a los que contienen condiciones iniciales, se tendrá:
(an sn + an−1 sn−1 + ... + a1 s + a0 )Y (s) − (an sn−1 + an−1 sn−2 + ... + a1 )y(0)−
−(an sn−2 + an−1 sn−3 + ... + a2 )y 0 (0) − ... − an y (n−1) (0) = X(s).
Si se agrupan todos los términos que incluyen condiciones iniciales como C(s),
se tendrá:
X(s) C(s)
Y (s) = + .
an sn + an−1 sn−1 + ... + a1 s + a0 an sn + an−1 sn−1 + ... + a1 s + a0
(an sn + an−1 sn−1 + ...a1 s + a0 )Y (s) = (bm sm + am−1 sm−1 + ...b1 s + b0 )X(s).
y entonces
Bm (s)
Y (s) = X(s) = H(s)X(s).
An (s)
Bm (s) Y (s)
H(s) = =
An (s) X(s)
es la función de transferencia del sistema, una función racional que re-
presenta la relación entrada-salida en el dominio s, obtenida cuando el sistema
tiene condiciones iniciales nulas.
Si la entrada que se aplica al sistema es de la forma es0 t , la salida será: H(s0 )es0 t .
La demostración es totalmente análoga a la hecha para las funciones ejω0 t , caso
CAPÍTULO 4. TRANSFORMADA DE LAPLACE 141
que está incluido en el anterior, ya que son equivalentes con σ0 = 0. Todo esto
es cierto para todo s del plano complejo.
Además, la función de transferencia se puede encontrar haciendo la transforma-
da de Laplace de la respuesta al impulso del sistema, y entonces, para el caso
de sistemas causales será:
H(s) = L{h(t)}.
H(ω) = H(s)|s=jω
La relación anterior se cumple para toda señal, no solamente para h(t). Luego,
para una señal cualquiera x(t), podrá afirmarse que:
X(ω) = X(s)|s=jω
Bm (s) Bm (s)
H(s) = =
An (s) (s − p1 )(s − p2 )...(s − pn )
y
Bx (s) Bx (s)
X(s) = = .
Ax (s) (s − q1 )(s − q2 )...(s − qr )
Se sabe que la salida del sistema, con condiciones iniciales nulas, en el dominio
transformado es:
Bm (s)Bx (s)
Y (s) = .
(s − p1 )(s − p2 )...(s − pn )(s − q1 )(s − q2 )...(s − qr )
4.8. Estabilidad
La definición de sistema estable dice que si la entrada es acotada, la salida que
produce también debe estar acotada (BIBO).
Hay diferentes formas de asegurar la estabilidad en este sentido:
Con respecto a su región de convergencia
Diremos que un sistema es estable si y solo si el eje imaginario (σ = 0)
está contenido en la región de convergencia de H(s). Esto es equivalente
a decir que todos los polos de H(s) deben estar en el semiplano izquierdo.
Con respecto a su ecuación caracterı́stica.
Las raı́ces de esta ecuación deben tener parte real negativa.
Con respecto a su respuesta al impulso
Si h(t) es absolutamente integrable, entonces el sistema es estable.
Por otro lado, el requerimiento de que el grado del polinomio del denomi-
nador de H(s) es mayor o igual al grado del polinomio del numerador está
relacionado con la causalidad del sistema.
Observe la forma que se usó para descomponer la salida del sistema en dos
sumatorias en la sección anterior. En los sistemas estables, los términos de la
primera sumatoria tienden a cero cuando se incrementa t, mientras que los
de la segunda sumatoria pueden contener términos que son distintos de cero
indefinidamente. Esta parte, contiene las singularidades de la entrada. Por lo
tanto, si la entrada es acotada, esta parte también lo será. Por consiguiente, la
salida será no acotada si al menos uno de los términos de la respuesta natural
(ci epi t ) se hace no acotado. Esto ocurre solamente si la parte real de al menos
uno de los polos pi es positivo o cero.
Ejercicio.
Proponga un sistema cuya función de transferencia tenga polos en el semiplano
izquierdo del plano complejo. Encuentre la respuesta al impulso y analice la
estabilidad de acuerdo al tercero de los criterios de arriba.
Repita para el caso en que los polos están en el semiplano derecho.
Para ver ejemplos y explicaciones:
Transformada de Laplace 2: https://youtu.be/ZLBqf4v2oGo
Capı́tulo 5
Variables de estado
Concepto de estado
Se define estado del sistema en algún tiempo t0 a la información que, junto con
todas las entradas para todos los tiempos subsiguientes a t0 , determina el com-
portamiento del sistema para t ≥ t0 . Es decir, es la información suficiente acerca
del sistema en algún instante t0 , en el sentido de que permite el cálculo de las
salidas del sistema para todo tiempo posterior a t0 . Las variables que contienen
143
CAPÍTULO 5. VARIABLES DE ESTADO 144
x1 (t) = i(t);
x2 (t) = vc (t).
que al despejar las derivadas de las variables de estado, permite llegar al siguiente
sistema acoplado de ecuaciones diferenciales de primer orden:
R 1 1
ẋ1 (t) = − x1 (t) − x2 (t) + vi (t);
L L L
1
ẋ2 (t) = x1 (t).
C
CAPÍTULO 5. VARIABLES DE ESTADO 145
Para completar el modelo se escribe la ecuación que relacione la salida del siste-
ma con las variables de estado. En este caso, la salida es vc (t). Luego, la ecuación
que falta es:
y(t) = x2 (t).
Estas dos ecuaciones matriciales son las ecuaciones de estado para el circuito
RLC. Esta representación no es la única posible ya que las variables de estado
pueden definirse arbitrariamente. Entonces, dependiendo de las variables que se
elijan, se tendrán diferentes modelos en variables de estado.
Ecuaciones de estado
En general, la forma estándar para las ecuaciones de estados de un SLIT de
múltiples entradas y múltiples salidas es:
Ejemplo 2
Considere ahora el sistema de la figura formado por dos resortes y dos masas,
al que se le aplica una fuerza externa f (t). La fuerza provoca el desplazamiento
de las masas m1 y m2 , z1 (t) y z2 (t), respectivamente.
Puede verificarse que el sistema queda representado por dos ecuaciones diferen-
ciales de segundo orden:
d2 z1 (t)
m1 = −k1 z1 (t) + k2 (z2 (t) − z1 (t))
dt2
CAPÍTULO 5. VARIABLES DE ESTADO 146
d2 z2 (t)
m2 = −k2 (z2 (t) − z1 (t)) + f (t).
dt2
La representación en variables de estado de este sistema estará dada por un
modelo con cuatro variables de estado.
Una elección clásica de variables para este tipo de sistemas es:
x1 (t) = z1 (t)
k1 + k2 k2
ẋ2 (t) = − x1 (t) + x3 (t)
m1 m1
ẋ3 (t) = x4 (t)
k2 k2
ẋ4 (t) = x1 (t) − x3 (t) + f (t)
m2 m2
2 0 −2 0 1
Para completar el modelo es necesario definir la salida del sistema. Puede to-
marse, por ejemplo, y(t) = z1 (t). Entonces, la ecuación faltante del modelo es:
y(t) = 1 0 0 0 x(t).
CAPÍTULO 5. VARIABLES DE ESTADO 147
x1 (t) = y(t)
ẋn (t) = y n (t) = −an−1 y (n−1) (t) − ... − a1 y 0 (t) − a0 y(t) + b0 u(t)
y(t) = x1 (t).
con
0 1 0 ... 0 0
0 0 1 ... 0 0
...
A= ...
; B=
0 0 0 ... 1 0
−a0 −a1 −a2 ... an−1 b0
C= 1 0 0 ... 0 ; D=0
Se pide:
1. Verificar que la ecuación diferencial es
C = 2 − 15 4 − 20 = −13 −16 ; D=5
Observe la relación directa que hay entre los coeficientes de la ecuación diferen-
cial y los elementos de las matrices A, B, C y D.
Las variables de estado son las salidas de los integradores, y las ecuaciones de
estados quedan:
ẋ1 (t) = x2 (t)
ẋ2 (t) = x3 (t)
...
ẋn (t) = −a0 x1 (t) − a1 x2 (t) − ... − an−1 xn (t) + u(t)
Y la ecuación de la salida es:
y(t) = (b0 − bn a0 )x1 (t) + (b1 − bn a2 )x2 (t) + ... + (bn−1 − bn an−1 )xn (t) + bn u(t)
Solución
El modelo representa al mismo sistema si tiene la misma ecuación diferen-
cial. Entonces, se trata de hallar la ecuación diferencial que corresponde
a este diagrama en bloques.
Para eso, se nombran dos variables auxiliares a la salida de ambos ambos
integradores: v1 y v2 . Luego se escriben las ecuaciones de los dos sumadores
usando el operador D para representar las derivadas.
CAPÍTULO 5. VARIABLES DE ESTADO 151
Dv2 (t) = 2u(t) − 3v1 (t) ⇒ v2 (t) = D−1 [2u(t) − 3v1 (t)],
Además,
y(t) = v1 (t)
Por lo tanto, la ecuación diferencial que corresponde a este modelo es
que es la misma que corresponde al modelo del punto 4.1. Ese modelo es
el de la segunda forma canónica y tiene su correspondiente diagrama en
bloques. Se verá ahora que este diagrama en bloques corresponde a otro
modelo en variables de estado del mismo sistema.
Como es usual, consideremos que las variables de estado son las funcio-
nes que están a la salida de los integradores, v1 y v2 . Y en este caso las
ecuaciones de estado se obtienen de manera directa, ya que v̇1 y v̇2 son las
salidas de ambos sumadores. Entonces:
Con el ejemplo anterior se mostró que un mismo sistema puede tener más de
un modelo en variables de estado y que cada modelo queda representado por su
correspondiente diagrama en bloques.
CAPÍTULO 5. VARIABLES DE ESTADO 152
ẋ(t) = Ax(t).
Se propone una solución que depende del estado del sistema en el tiempo inicial,
t0 :
xl (t) = eA(t−t0 ) x(t0 ),
donde x(t0 ) es el valor inicial de los estados.
Si t0 = 0 la ecuación es:
xl (t) = eAt x(0),
con eAt = φ(t) la matriz de transición de los estados.
Por definición:
∞ k
X t
eAt = Ak ,
k!
k=0
k
donde A = AA...A, k veces.
Una caracterı́stica de esta matriz es que
deAt
= AeAt = eAt A
dt
Demostración
∞ ∞ ∞ ∞
deAt d X tk k X d tk k X k k−1 k X tk−1
= A = A = t A = AAk−1 =
dt dt k! dt k! k! (k − 1)!
k=0 k=0 k=0 k=1
∞ k
X t
=A Ak = AeAt .
k!
k=0
CAPÍTULO 5. VARIABLES DE ESTADO 153
Ahora puede verificarse que xl (t) = eAt x(0) es solución de ẋ(t) = Ax(t).
Para eso, simplemente se reemplazan xl (t) y su derivada en la ecuación diferen-
cial, dando la identidad:
x(t) = φ(t)x(0),
x(0) = φ(0)x(0),
Para probar que φ(t1 + t2 ) = φ(t1 )φ(t2 ) puede seguirse el mismo procedi-
miento que antes pero yendo primero de t = 0 a t = t2 , y luego de t = t2
a t = t1 + t2 .
3. φ−1 (t) = φ(−t)
Demostración
Otra vez se considera la solución
Estabilidad
Un SLIT es estable si los autovalores de la matriz de estados (A) tienen parte
real negativa. Para fundamentar esta propiedad debe analizarse la relación de los
autovalores con la ecuación caracterı́stica del sistema. Este análisis se completará
cuando se estudien los modelos en variables de estado en el dominio de Laplace.
Ejemplo 5
Analice la estabilidad del SLIT cuyo modelo está dado por las siguientes matri-
ces:
0 1 0
A= ; B=
−2 −3 1
CAPÍTULO 5. VARIABLES DE ESTADO 156
C= 4 5 ; D=0
det(λI − A) = 0
λ 0 0 1
det − =0
0 λ −2 −3
λ −1
det =0
2 λ+3
λ2 + 3λ + 2 = 0,
de donde sale que los autovalores son λ1 = −2 y λ2 = −1, que son valores reales
negativos y, por lo tanto, el sistema es estable.
La respuesta al impulso se calcula de acuerdo a la fórmula obtenida antes:
ẋ1 (t) = a11 x1 (t) + a12 x2 (t) + ... + a1n xn (t) + b11 u1 (t) + ... + b1r ur (t)
ẋ2 (t) = a21 x1 (t) + a22 x2 (t) + ... + a2n xn (t) + b21 u1 (t) + ... + b2r ur (t)
donde aij y bij son los elementos de las matrices A y B, respectivamente.
Ahora, aplicando la transformada de Laplace en ambos miembros de las dos
ecuaciones resulta:
sX1 (s) − x1 (0) = a11 X1 (s) + a12 X2 (s) + ... + a1n Xn (s) + b11 U1 (s) + ... + b1r Ur (s)
sX2 (s) − x2 (0) = a21 X1 (s) + a22 X2 (s) + ... + a2n Xn (s) + b21 U1 (s) + ... + b2r Ur (s)
Entonces, para hallar la matriz de transición de los estados, aunque hay otros
métodos (ver Apéndice D), el método más sencillo consiste en hacer:
Ejemplo 6
Verificar que la matriz de transición de los estados del modelo del Ejemplo 5 es
la dada en ese enunciado.
Solución −1
−1 s −1
φ(s) = (sI − A) =
2 s+3
T
1 s+3 −2 1 s+3 1
φ(s) = =
(s2 + 3s + 2) 1 s (s2 + 3s + 2) −2 s
Luego, !
s+3 1
(s+2)(s+1) (s+2)(s+1)
φ(s) = −2 s .
(s+2)(s+1) (s+2)(s+1)
Ejemplo 7
Halle la salida para el sistema anterior si u(t) = us (t) y las condiciones iniciales
de los estados son x1 (0) = 2 y x2 (0) = −1.
Solución
1. Solución libre de los estados
Como ya se tiene φ(t), esta solución puede calcularse directamente en el
dominio temporal:
xl (t) = φ(t)x(0)
2e − e−2t
−t
e−t − e−2t
2
xl (t) =
−2e−t + 2e−2t −e−t + 2e−2t −1
4e−t − 2e−2t − e−t + e−2t
−t
3e − 1e−2t
xl (t) = =
−4e−t + 4e−2t + e−t − 2e−2t −3e−t + 2e−2t
Puede observarse que la solución de los estados es un vector (2x1) en este
caso, ya que tenemos un sistema de segundo orden con una entrada.
CAPÍTULO 5. VARIABLES DE ESTADO 159
5s + 4 2 1 3
Yf (s) = = + −
(s + 2)(s + 1)s s s+1 s+2
yf (t) = 2 + e−t − 3e−2t ; t≥0
v(t) = Qx(t),
4. Hv (s) = H(s)
Como es de suponer, si la transformación de semejanza no cambia la re-
lación entrada-salida del sistema, entonces la función de transferencia no
debe verse modificada.
Para la siguiente demostración se tiene en cuenta la propiedad de las
matrices: (MN)−1 = N−1 M−1 .
Demostración
Ejemplo 8
Para el mismo sistema lineal de segundo orden considerado en el ejemplo 4.1,
suponga la transformación dada por:
y entonces
2 −1
P= .
−1 1
3 3 2 5s + 4
Hv (s) = − + =
s + 2 (s + 1)(s + 2) s + 1 (s + 1)(s + 2)
Si se hace el desarrollo en fracciones simples queda:
6 1
Hv (s) = − ,
s + 2 (s + 1)
por lo que
h(t) = (6e−2t − e−t )us (t).
(λ + 2)(λ + 1) = λ2 + 3λ + 2
CAPÍTULO 5. VARIABLES DE ESTADO 164
2 = 2 = (−1)(−2)
−3 = −3 = −2 − 1
Hv (s) = H(s)
Quedó verificado en el punto 3.
Transformación de autovalores
La transformación de autovalores es un caso particular de transformación de
semejanza en el que la matriz de estados transformada es una matriz diagonal,
que en este caso se denomina Λ.
Si A tiene n autovalores distintos, entonces puede diagonalizarse usando la
transformación
Λ = P−1 AP,
siendo P una matriz cuyas columnas son los autovectores de la matriz A, resul-
tando la matriz Λ una matriz diagonal cuyos elementos son los autovalores de
A.
Prueba
Para cada autovalor puede escribirse la ecuación
Aχ1 = λ1 χ1
Aχ2 = λ2 χ2
...
Aχn = λn χn
donde χi es el autovector correspondiente al autovalor λi , con i = 1, 2, ..., n.
Ahora, las n ecuaciones de arriba, pueden escribirse en forma matricial como:
AP = PΛ,
Λ = P−1 AP.
CAPÍTULO 5. VARIABLES DE ESTADO 165
Ejemplo 9.
Encuentre el modelo de variables de estado semejante al dado en el ejemplo 4.1,
realizando la transformación por los autovectores.
Solución
Para este sistema ya se halló un modelo semejante en el ejemplo 8. En este
caso vamos a encontrar otro modelo semajante aplicando la transformación de
autovalores.
Ya se conocen los autovalores de este sistema: λ1 = −1 y λ2 = −2. Lo siguiente
es calcular los autovectores, χi , que conformarán la mátriz de transformación y
que responden a la ecuación:
χ11 = −χ12
y su inversa,
−1 2 1
P = .
−1 −1
CAPÍTULO 5. VARIABLES DE ESTADO 166
Hv (s) = Cv φv (s)Bv + Dv
1
s+1 0
1 = −1 −6 1
Hv (s) = −1 −6 s+1 s+2
1 −1 −1
0 s+2
6 1 5s + 4
Hv (s) = − = ,
s+2 s+1 (s + 1)(s + 2)
que es la misma función de transferencia del resto de los modelos semejantes
que representan al mismo sistema (misma relación entrada-salida).
Operador anulador
con cierta entrada x(t) y condiciones iniciales y(0), y 0 (0), ..., y (n−1) (0).
El procedimiento requiere, como primer paso, encontrar un operador, LA , que
anule la función particular dada, x(t). Es decir,
LA [x(t)] = 0.
LA {A(D)[y(t)]} = LA [x(t)],
167
APÉNDICE A. OPERADOR ANULADOR 168
y obtenemos
LA {A(D)[y(t)]} = 0.
Ahora se tiene una ecuación diferencial homogénea pero de mayor orden que la
ecuación original. Por lo tanto, la nueva ecuación caracterı́stica asociada tendrá,
además de las n raı́ces originales, nuevas raı́ces que le imponga el operador
anulador. Entonces, la solución propuesta será:
c1 y1 (t) + c2 yn (t) + ... + cn yn (t) + cL1 yL1 (t) + cL2 yL2 (t) + ... + cLr yLr (t).
LA {A(D)[y(t)]} = LA {A(D)[cL1 yL1 (t) + cL2 yL2 (t) + ... + cLr yLr (t)]} = 0.
Los términos que quedan son los que provienen de la función particular de
excitación (entrada) y corresponden a la solución particular estudiada antes.
En este método no es necesario proponer una forma para la solución particular.
Sı́ es necesario, por el contrario, encontrar la forma del operador anulador.
Ejemplo. Considere la ecuación diferencial A(D)[y(t)] = (D2 + 1)[y(t)] = sen(t),
con condiciones iniciales y(0) = 1; y 0 (0) = 0.
Solución. En primer lugar se busca el operador que anule a la función x(t). En
este caso:
LA = D2 + 1.
(r2 + 1)(r2 + 1) = 0.
Claramente las raı́ces son j y −j, dobles. Entonces, la solución que se propone
es la que corresponde a raı́ces complejas conjugadas múltiples:
(D2 + 1){c1 cos (t) + c2 sen(t)} + (D2 + 1){c3 t cos (t) + c4 tsen(t)} = sen(t).
APÉNDICE A. OPERADOR ANULADOR 169
D2 {c3 t cos (t) + c4 tsen(t)} = −2c3 sen(t) − c3 t cos (t) + 2c4 cos (t) − c4 tsen(t).
Ası́,
(D2 + 1){c3 t cos (t) + c4 tsen(t)} = sen(t)
puede ponerse como:
y(0) = c1 = 1;
1 1
y 0 (0) = c2 − = 0 ⇒ c2 = .
2 2
Y ası́ se obtiene la solución completa de esta ecuación diferencial con la entrada
y las condiciones inciales dadas:
1 1
y(t) = cos (t) + sen(t) − t cos (t).
2 2
Apéndice B
Fórmula de Leibniz
170
Apéndice C
171
APÉNDICE C. ALGUNAS FORMAS DE ONDA 172
Apéndice D
Cálculo de φ(t)
Método de Cayley-Hamilton
Este método permite el cálculo directo de la matriz de transición, sin pasar por
el dominio de la transformada de Laplace. Sin embargo, por lo general tiene una
mayor dificultad operativa.
Para A diagonal, con elementos diferentes aii , i = 1, 2, ..n, puede demos-
trarse que φ(t) siempre será también una matriz nxn diagonal, cuyos
elementos son eaii t , i = 1, 2, ..n.
Para A cualquiera, debe calcularse mediante alguna de las siguientes ex-
presiones:
∞ k
X t k
eAt = A
k!
k=0
n−1
X
eAt = γi (t)Ai
i=0
173
APÉNDICE D. CÁLCULO DE φ(T ) 174
con j = 1, 2, ..., n.
El método se aplica directamente cuando todos los autovalores son diferentes.
Ejemplo D.1.
Hallar la matriz de transición de los estados del modelo cuya matriz de estados
es
−3 1 0
A = 1 −3 0
0 0 −3
λ3 + 9λ3 + 26λ + 24 = 0,
de donde salen los tres autovalores;
con respecto de λj :
teλj t = γ1 (t) + 2λj γ2 (t),
lo que nos da la tercera ecuación: