Numeros Complejos y Funciones Complejas Maquina Alfa Oficial SEGUNDA SEMANA PDF
Numeros Complejos y Funciones Complejas Maquina Alfa Oficial SEGUNDA SEMANA PDF
Numeros Complejos y Funciones Complejas Maquina Alfa Oficial SEGUNDA SEMANA PDF
y Cálculo Operacional
Teorı́a y Práctica usando MATLAB
William La Cruz
Universidad Central de Venezuela
Facultad de Ingenierı́a
Escuela de Ingenierı́a Eléctrica
Departamento de Electrónica, Computación y Control
c 2015 W. La Cruz
Contenido
Prólogo vi
I Variable Compleja 2
1 Números Complejos 3
1.1 Definición . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2 Operaciones Algebraicas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.3 Representación Geométrica . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.4 Valor Absoluto y Conjugado . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.5 Coordenadas Polares . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.5.1 Fórmula de Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.6 Potencias y Raı́ces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.7 Regiones en el Plano Complejo . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.8 Aritmética compleja con MATLAB . . . . . . . . . . . . . . . . . . . . . . . 16
1.9 Problemas Resueltos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.10 Problemas Propuestos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
ii
iii
8 Transformada z 237
8.1 Definición . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 237
8.1.1 Región de Convergencia . . . . . . . . . . . . . . . . . . . . . . . . . 238
8.1.2 Transformada z Inversa . . . . . . . . . . . . . . . . . . . . . . . . . 239
8.2 Propiedades de la Transformada z . . . . . . . . . . . . . . . . . . . . . . . 239
8.3 Algunos Pares de Transformadas . . . . . . . . . . . . . . . . . . . . . . . . 245
8.4 Cálculo de la Transformada z Inversa . . . . . . . . . . . . . . . . . . . . . . 250
8.4.1 Integración Compleja . . . . . . . . . . . . . . . . . . . . . . . . . . 250
8.4.2 Expansión en Serie de Potencias . . . . . . . . . . . . . . . . . . . . 251
8.4.3 Inversión por Tablas . . . . . . . . . . . . . . . . . . . . . . . . . . . 252
Contenido v
Bibliografı́a 269
Prólogo
“La persona que nunca se equivoca
nunca prueba algo nuevo”
Albert Einstein
Hoy en dı́a la creencia de los estudiantes, y en general de las personas, es que los proble-
mas asociados a matemáticas, ingenierı́a o a cualquier ciencia básica o aplicada, se pueden
resolver utilizando una computadora sin otro conocimiento previo o, simplemente, las so-
luciones de los mismos se pueden encontrar en Iternet. A pesar de que esta visión no está
muy alejada de la realidad, debido a los grandes avances en la ciencia de la computación,
no es cierta. Sólo basta con pensar y hacerse las siguientes preguntas: ¿cómo se diseñan
las computadoras? o ¿cómo funciona Internet? Las respuestas a estas preguntas encierra
un gran mundo de conceptos y métodos, todos ellos relacionados con la matemática y la
ingenierı́a. Por ejemplo, para diseñar un procesador de una computadora o diseñar una
red de comunicaciones, es necesario tener conocimiento de la suma, la multiplicación, la
derivación, la integración, entre otros conceptos y operaciones matemáticas.
El objetivo principal de este escrito es presentar de manera formal (usando solo concep-
tos matemáticos correspondientes al nivel del cuarto semestre de la carrera de Ingenierı́a)
los conceptos de Variable Compleja y del Cálculo Operacional. Especı́ficamente, en este
escrito se presentan los conceptos básicos del cálculo diferencial e integral en el conjunto de
los números complejos, comenzando con la definición de número complejo, pasando luego
con la descripción formal de funciones de variable compleja, enfatizando en los conceptos
de continuidad, diferenciación e integración.
También se presenta una breve introducción al cálculo operacional, en la cual se des-
criben las funciones elementales de dominios continuo y discreto, y se definen las trans-
formadas de Fourier y Laplace (para funciones de dominio continuo) y la transformada z
(para funciones de dominio discreto). Además, se describen y demuestran las propiedades
de la transformada de Fourier, la transformada de Laplace y la transformada z.
Para la comprensión de los conceptos expuestos en el escrito, es necesario que el lector
posea un nivel aceptable de conocimientos del Cálculo Elemental (particularmente, del
Cálculo Infinitesimal Real en una y dos variables); también debe manejar con destreza las
operaciones algebraicas en los conjuntos reales, racionales y enteros.
El escrito está estructurado en dos partes: I. Variable Compleja y II. Cálculo Operacio-
nal, ya que el mismo, en principio, está dirigido a los estudiantes del curso Variable Com-
pleja y Cálculo Operacional, que forma parte del Plan de Estudios de Ingenierı́a Eléctrica.
Es de notar que estas notas han evolucionado a medida que transcurren los semestres en
los cuales se ha dictado la asignatura Variable Compleja y Cálculo Operacional, agregando
un problema resuelto o una imagen que permita visualizar mejor un concepto. En todo
caso, siempre en mejora del aprendizaje.
Ahora creemos necesario agregar un enfoque computacional del análisis complejo, es
decir, en el escrito también se presenta una visión computacional de los conceptos, usando
vi
PRÓLOGO 1
para ello el software MATLAB que es empleado en numerosas universidades y posee mu-
chas facilidades gráficas de calidad, flexibilidad e interactividad. Para una introducción a
MATLAB, puede consultar el material Cálculo Cientı́fico con MATLAB, disponible en el
link: http://neutron.ing.ucv.ve/cruzw/_private/CalculoCientificoconMatlab.pdf o el libro de Higham
y Higham [4]. Al final de cada capı́tulo se presentan ejemplos prácticos codificados en
MATLAB. Tal iniciativa tiene el propósito de introducir al estudiante en el uso de una
herramienta computacional muy poderosa como MATLAB, pero dejándoles ver que tal
herramienta por sı́ sola no resuelve los problemas, sino que es necesario que el usuario
tenga la destreza de entender y definir correctamente el problema para luego resolverlo con
la computadora.
Parte I
Variable Compleja
2
1
Números Complejos
En este capı́tulo se introducen algunas propiedades del conjunto de los números com-
plejos. Tal conjunto de números es ampliamente utilizado en el desarrollo de las ideas
teóricas de Ingenierı́a, en especial, de Ingenierı́a Eléctrica.
1.1 Definición
Se dice que z es un número complejo si se expresa como
z = x + iy
o, de manera equivalente,
z = x + y i,
donde x e y son números reales (recuerde que R denota el conjunto de los números reales).
El sı́mbolo i se conoce como unidad imaginaria. El conjunto de los números complejos se
denota como C, y z ∈ C indica que z pertenece al conjunto de los números complejos.
Por otra parte, un número complejo z = x + i y también se puede definir como el par
ordenado z = (x, y). Esta definición equivalente de número complejo permite definir a i
como el número complejo dado por
i = (0, 1).
Más adelante se demostrará que i2 = −1.
Se denota con x = Re z la parte real del número complejo z, y con y = Im z la parte
imaginaria de z. Los números complejos de la forma x + i 0, se denominan reales puros
o, simplemente, reales. Además, los números complejos de la forma 0 + i y, se denominan
imaginarios puros.
Los números reales 0 y 1, también se pueden definir como números complejos. El cero
de los números complejos, denotado 0, se define como 0 = 0 + i 0. El 1 de los números
complejos, denotado 1, se define por 1 = 1 + i 0.
Definición 1.1. Se dice que dos números complejos z y w son iguales si, y sólo si sus partes
reales son iguales y sus partes imaginarias son iguales. En otras palabras, si z = x + i y y
w = u + i v, entonces z = w, si y sólo si Recuerde que
los números
x=u e y = v. reales se pueden
representar
como un punto
En particular, en una recta,
z = x + iy = 0 ⇔ x = y = 0. que se recorre
de izquierda a
derecha, con lo
Observación 1.1. No existe relación de orden en los números complejos. En los números cual se impone
un orden.
reales, por ejemplo, se tiene que 5 > 3, pero no tiene sentido afirmar que 1 + i < 2 + i 3.
3
4 1.2. OPERACIONES ALGEBRAICAS
z1 · z2 = z1 z2 = (x1 x2 − y1 y2 ) + i (x1 y2 + x2 y1 ).
1. Conmutativa
• z+w =w+z
• zw = wz
2. Asociativa
• z + (w + s) = (z + w) + s
• z(ws) = (zw)s
3. Elemento Neutro
• z+0=z
• 1·z =z
CAPÍTULO 1. NÚMEROS COMPLEJOS 5
4. Elemento Inverso
5. Distributiva
• z(w + s) = zw + zs.
a) z1 + z2 + z3
b) z1 · z2 · z3
c) z1 (z2 + z3 )
d) z2 /z3
e) z1 /(z2 · z1 )
f) z3 /(z1 /z2 )
z
= z w−1 .
w
z
Demostración. Demostremos que las partes real e imaginaria de los números complejos
w
y z w−1 6 0. Se tiene
son iguales. Tomemos z = x + i y, w = u + i v y asumamos que w =
que
z xu + yv yu − xv
= 2 2
+i 2 . (1.1)
w u +v u + v2
Ahora, veamos la forma que tiene w−1 . Asumamos que w−1 = a+i b. Por la Proposición 1.1
se tiene que w w−1 = 1. Por tanto, podemos escribir:
z = (x, y)
y b
x Eje real
Eje imaginario
z = x + iy
y
x Eje real
inicial es el origen (0, 0) y su extremo final es el punto (x, y). De esta forma, a un número ¿Por qué se
toma a z como
complejo lo podemos representar como un par ordenado o como un vector en el plano xy. el vector
orientado con
origen (0, 0) y
1.4 Valor Absoluto y Conjugado extremo final
(x, y), y no otro
vector?
En la sección anterior notamos que un número complejo z se puede representar como
un vector orientado. Este hecho permite que z herede algunas propiedades de los vectores,
entre ellas la magnitud o valor absoluto.
Definición 1.6 (Valor absoluto). El valor absoluto del número complejo z = x + i y,
denotado por |z|, se define como la longitud del vector z, la cual se calcula a través de la
fórmula p
|z| = x2 + y 2 .
Geométricamente, el valor absoluto |z| representa la distancia del punto (x, y) al origen
(0, 0).
Definición 1.7 (Conjugado). El conjugado del número complejo z = x + i y, denotado
por z, se define como
z = x + i (−y).
Geométricamente, el conjugado z es la reflexión del punto (x, y) con respecto al eje real.
En la Figura 1.3 se muestra un ejemplo gráfico del conjugado de un número complejo.
z = x + iy
y b
−y b
z = x − iy
2. z1 ± z2 = z1 ± z2 .
3. z1 z2 = z1 z2 .
z1 z1
4. = , siempre y cuando z2 6= 0.
z2 z2
8 1.4. VALOR ABSOLUTO Y CONJUGADO
5. |z1 | = |z1 |.
6. z1 z1 = |z1 |2 .
z1
7. z1−1 = , siempre y cuando z1 6= 0.
|z1 |2
8. |z1 | ≥ |Re z1 | ≥ Re z1 .
9. |z1 | ≥ |Im z1 | ≥ Im z1 .
pero
z1 z2 + z1 z2 = 2 Re (z1 z2 ) ≤ 2 |z1 z2 | = 2 |z1 | |z2 | = 2 |z1 | |z2 |,
luego
|z1 + z2 |2 ≤ |z1 |2 + 2 |z1 | |z2 | + |z2 |2 = (|z1 | + |z2 |)2 ,
de donde se deduce que
|z1 + z2 | ≤ |z1 | + |z2 |,
con lo cual queda establecida la desigualdad triangular.
z = x + iy
|z
=
r
Los valores de r y θ = arg z definen de manera única a z, es decir, para cada par (r, θ)
existe un único número complejo z que tiene como coordenadas polares a (r, θ). En cambio,
el número complejo z caracteriza de manera única a r, pero no a θ = arg z, esto es, dado
un número complejo z, entonces existe un único r > 0 tal que r = |z|, pero existen infinitos
valores de θ = arg z. Observe el siguiente ejemplo.
10 1.5. COORDENADAS POLARES
Ejemplo 1.1. Para el número complejo z = 1, se tiene que los valores de θ están dados
por:
θ = arg z ∈ {0, ±2π, ±4π, . . .}.
La existencia de infinitos valores de θ es una consecuencia de la periodicidad de las
funciones trigonométricas. En efecto, si z1 = r1 (cos θ1 + i sen θ1 ), entonces se tiene que
−π < Arg z ≤ π.
Note que la fórmula anterior considera la posición z en el plano complejo para calcular
el valor de Arg z, esto es, dependiendo del cuadrante donde se encuentre z, su argumento
principal se calcula de una forma muy particular, por supuesto, usando para ello el valor
de arctan(y/x). Observe el siguiente ejemplo.
Ejemplo 1.2. Utilizando el argumento principal, represente en forma polar los siguientes
números complejos: z1 = 1 + i , z2 = −1 + i , z3 = −1 − i, y z4 = 1 − i.
√ √
Solución. Se tiene que el valor absoluto de z1 es r = |z1 | = 11 + 12 = 2. Como
z1 = 1 + i está en el primer cuadrante, su argumento principal es
π
Arg z1 = arctan(1) = .
4
CAPÍTULO 1. NÚMEROS COMPLEJOS 11
Por un razonamiento similar, se puede verificar que la forma polar de los números complejos
z2 = −1 + i , z3 = −1 − i y z4 = 1 − i, es respectivamente:
√
z2 = 2(cos(3π/4) + i sen (3π/4)),
√
z3 = 2(cos(3π/4) − i sen (3π/4)),
√
z4 = 2(cos(π/4) − i sen (π/4)).
Por otra parte, pasemos a describir algunas propiedades del argumento, las cuales se
muestran en la siguiente proposición.
que es la forma polar de 1/z2 y arg (1/z2 ) = −θ2 ; luego, se cumple la identidad 2.
Como
z1 1
= z1 · ,
z2 z2
entonces, por las identidades 1 y 2, se prueba que la identidad 3 se cumple.
12 1.5. COORDENADAS POLARES
Ejemplo 1.3. La Proposición 1.6 se puede utilizar para calcular un valor del argumento
de z1 z2 o z1 /z2 . Por ejemplo, tomemos z1 = i y z2 = √ −1 + i. La forma polar de estos
números complejos es z1 = cos(π/2) + i sen (π/2) y z2 = 2(cos(3π/4) + i sen (3π/4)), por
supuesto, usando el argumento principal. Ası́, un valor del argumento de z1 z2 es
π 3π 5π
arg (z1 z1 ) = + = ,
2 4 4
Un ejercicio
interesante es lo cual efectivamente es un valor del argumento de z1 z2 = −1−i (verifique esta afirmación).
encontrar z1 y Pero, este valor no es el argumento principal de z1 z2 = −1 − i, aunque π/2 y 3π/4 sı́ son
z2 tales que
Arg (z1 z2 ) = los argumentos principales de z1 y z2 , respectivamente.
Arg z1 + Arg z2
El Ejemplo 1.3 nos permite asegurar que, en general,
z = reiθ
Observe que en la forma polar de los números dados se utilizó el argumento principal, pero
pudiera también emplearse cualquier valor del argumento.
CAPÍTULO 1. NÚMEROS COMPLEJOS 13
Definición 1.11 (Raı́ces n-ésimas). Sean z = r eiθ y n un entero positivo. Las raı́ces
n-ésimas de z, denotadas por z 1/n , se definen como
1/n √
n i θ+2kπ √
n
θ + 2kπ θ + 2kπ
z = re n = r cos + i sen ,
n n
√
para k = 0, 1, . . . , n − 1, donde n
r denota la raı́z n-ésima positiva del número real r.
√
2 3π √ i π +2kπ
4
(1 + i)3 = 23 ei (1 + i)1/5 =
10
4 y 2e 5 , k = 0, 1, . . . , 4.
b
z
b
z0
¿Todo conjunto
Definición 1.13 (Conjunto abierto). Se dice que un conjunto de números complejos S es
del plano que abierto, si para cada z ∈ S existe una vecindad de z, B(z, ε), tal que B(z, ε) ⊂ S.
no contiene sus
bordes es Ejemplo 1.6. Los siguientes conjuntos del plano complejo son ejemplos de conjuntos
abierto?
abiertos.
a) S1 = {z ∈ C : |z| < 1}
b) S2 es la región del plano complejo formada por los puntos interiores del cuadrado
cuyos vértices son los puntos z1 = 0, z2 = 1, z3 = 1 + i, y z4 = i.
c) S3 es la región del plano complejo formada por los puntos z = x + i y que satisfacen
el siguiente sistema de inecuaciones:
x − y > −1
x+y < 1
y > 0
1 1 S2 1
S1
S3
b
−1 1 1 −1 1
−1
a) b) c)
S z1
b
z2
b
z3
b
Ejemplo 1.8. Considere el conjunto S de la Figura 1.7. Aquı́ observamos que z1 y z2 son
puntos de acumulación de S, pero z3 no es un punto de acumulación de S, ya que existe
ε > 0 tal que B(z3 , ε) ∩ S = ∅.
¿Todo conjunto
del plano que
Definición 1.16 (Conjunto cerrado). Se dice que un conjunto de números complejos S es contiene sus
cerrado, si todos sus puntos de acumulación pertenecen a él. De forma equivalente, S es bordes es
cerrado?
cerrado si contiene a su frontera.
Ejemplo 1.9. Considere todos los conjuntos mostrados en las Figuras 1.5, 1.6 y 1.7.
Suponga que con estos conjuntos se forman nuevos conjuntos que contienen los puntos de
la frontera. Estos últimos conjuntos son todos cerrados.
S1 S2
Ejemplo 1.11. Todos los conjuntos mostrados en las Figuras 1.5, 1.6 y 1.7 son dominios.
Asimismo, el conjunto S1 de la Figura 1.8 también es un dominio.
S2
S1
>> z1 = 1+3 i , z2 = -1 -2 i
z1 =
1.0000 + 3.0000 i
z2 =
-1.0000 - 2.0000 i
>> z1 + z2
ans =
0 + 1.0000 i
>> z1 - z2
ans =
2.0000 + 5.0000 i
>> z1 * z2
ans =
5.0000 - 5.0000 i
>> z1 / z2
ans =
-1.4000 - 0.2000 i
1
Para una introducción a MATLAB, puede consultar el material Cálculo Cientı́fico con MATLAB dispo-
nible en el link: http://neutron.ing.ucv.ve/cruzw/wlacruz_cursos_archivos/CalculoCientificoconMatlab.pdf
o el libro de Higham y Higham [4].
CAPÍTULO 1. NÚMEROS COMPLEJOS 17
Los comandos: real(z), imag(z), conj(z), abs(z) y angle(z), se utilizan para cal- ¿Qué resultado
arrojarı́a el
cular la parte real, la parte imaginaria, el conjugado, el módulo y el argumento de z, comando
respectivamente. Observe el siguiente ejemplo. angle(i)?
>> z = 1+ i /2
z =
1.0000 + 0.5000 i
>> real ( z )
ans =
1
>> imag ( z )
ans =
0.5000
>> conj ( z )
ans =
1.0000 - 0.5000 i
>> abs ( z )
ans =
1.1180
>> angle ( z )
ans =
0.4636
Note que el valor angle(z) está expresado en radianes, además, es el argumento principal
de z. Observe el siguiente ejemplo.
>> z = 1+ i ; angle ( z )
ans =
0.7854
>> z = -1+ i ; angle ( z )
ans =
2.3562
>> z = -1 - i ; angle ( z )
ans =
-2.3562
>> z = 1 -i ; angle ( z )
Para más
ans =
detalles de
-0.7854 plot, teclee
help plot en el
escritorio de
MATLAB
El comando plot puede usarse para ver la representación gráfica de un número complejo
como un par ordenado. Por ejemplo, con las instrucciones
z = 1/3 + sqrt (2) * i ;
plot (z , ’ob ’ , ’ M a r k e r E d g eCo lor ’ , ’b ’ , ’ M a r k e r F a c eCol or ’ , ’b ’ ,...
’ MarkerSize ’ ,5)
text (0.4 ,1.6 , ’ $z = \ frac {1}{3} + \ sqrt {2}\ , i$ ’ ,...
’ FontSize ’ ,14 , ’ B a c k g r o u nd Col or ’ , ’w ’ , ’ Interpreter ’ , ’ latex ’)
xlabel ( ’x ’)
ylabel ( ’y ’)
grid on
axis square
18 1.8. ARITMÉTICA COMPLEJA CON MATLAB
2.5
1
√
z= 3 + 2i
1.5
y
1
0.5
0
−1 −0.5 0 0.5 1 1.5
x
1
√
Figura 1.10. Representación como un par ordenado del número z = 3 + 2i
√
se genera la gráfica de la Figura 1.10, donde se representa al número complejo z = 31 + 2i
como un par ordenado.
El comando compass(z) dibuja a z = x + i y como un vector con extremo inicia en el
origen y extremo final en el par ordenado (x, y). Las instrucciones
Para más zz = [1+i , -1 -i , -(1/2)+i , 1 - i ];
detalles de
compass, teclee
c o m p a s s( zz )
help compass
en el escritorio se utilizaron para representar como vectores a los números complejos z1 = 1+i, z2 = −1−i,
de MATLAB
z3 = − 12 + i, z4 = 1 − i, cuya gráfica se aprecia en la Figura 1.11.
90
1.5
120 60
150 30
0.5
180 0
210 330
240 300
270
Conocidas las coordenadas polares, (r, θ), de un número complejo z, se puede utilizar
a la función exp para hallar la forma cartesiana de z. Por ejemplo, la forma cartesiana,
z = x + iy, del número complejo z = 2 eiπ/3 , que está representado en forma polar usando
la fórmula de Euler, se halla de la siguiente manera:
El comando z^n permite calcular la potencia enésima del número complejo z. Observe
el siguiente ejemplo donde se calcula (1 + i)5 .
>> (1+ i ) ^5
ans =
-4.0000 - 4.0000 i
Ahora, el Listado 1.1 muestra la función raizn que calcula las n raı́ces enésimas de z.
f u n c t i o n szn = raizn (z , n )
% raizn Halla las n raı́ces de z .
% szn es un arrego que c o n t i e n e las n raı́ces de z
theta = angle ( z ) ;
r = abs ( z ) ;
k = 0:1: n -1;
szn = ( r ^(1/ n ) ) * exp (( theta + 2* k * pi ) * i / n ) ;
end
Ejercicio 1.3. Emplee convenientemente las funciones √ compass y raizn para obtener la
gráfica de las raı́ces séptimas de la unidad, esto es, 7 1.
√ √
Problema 1.2. Demuestre que |(2 z + 5)( 2 − i)| = 3|2z + 5|.
|z1 | |z1 |
≤ ,
|z2 + z3 | ||z2 | − |z3 ||
lo que implica
z1 |z1 |
z2 + z3 ≤ ||z2 | − |z3 || ,
Problema 1.5. Si z1 z2 6= 0, aplicar la forma polar para demostrar que Re (z1 z2 ) = |z1 | |z2 |
si, y sólo si θ1 + θ2 = 2nπ, (n = 0, ±1, ±2, . . .), donde θ1 = arg z1 y θ2 = arg z2 .
1 − z n+1
1 + z + z2 + · · · + zn = ,
1−z
se satisface para todo entero n ≥ 1 y todo z 6= 1.
a) (2 + 3i) + (4 + i) e) (8 + 6i)3
2
2 + 3i 3
b) f) 1 +
4+i 1+i
1 3 1 + 2i 2 − i
c) + g) +
i 1+i 3 − 4i 5i
d) (2 + 3i)(4 + i) h) (1 − i)4
a) z 2 = 3 − 4i b) (z + 1)2 = 3 + 4i
a) z 5 − 2 = 0 c) z 6 + 8 = 0
b) z 4 + i = 0 d) z 3 − 4 = 0
1.8. Encontrar todos los valores del argumento de z cuando el número complejo z está
dado por:
√ √
i 3 + 3 + (3 − 3)i
a) z = c) z =
i−1 3 + 3i
√ √
(2 3 + 1) + ( 3 − 2) i (3 + i)
b) z =
5 − 5i
2.1 Definición
Una función f de variable compleja es una regla de asignación que le hace corresponder
a un número complejo z = x + i y uno o varios números complejos w = u + i v. El
número w se llama valor o imagen de f en z y se designa por f (z), es decir, w = f (z) o,
equivalentemente, u + i v = f (x + i y).
Definición 2.1 (Funciones monovaluadas y multivaluadas). Sea f (z) una función de va-
riable compleja. Si a cada número complejo z la función f le hace corresponder una y sólo
una imagen w, se dice que f es monovaluada. Ahora, si la función f le hace corresponder ¿Todas las
a z más de una imagen, digamos w1 , w2 , . . ., se dice que f es multivaluada. funciones vistas
en los cursos de
En la Figura 2.1 se muestra una representación gráfica, a través de diagramas, de Cálculo son
monovaluadas?
funciones monovaluadas y multivaluadas.
w = f (z) w = f (z) b w1
b b b
b w2
z w z ..
.
b
wn
p(z) = a0 + a1 z + a2 z 2 + · · · + an z n , (an 6= 0)
se denomina polinomio complejo o, simplemente, polinomio de grado n.
23
24 2.2. LÍMITE Y CONTINUIDAD
Recuerde el Ejercicio 2.1. Demuestra que tanto el dominio como el rango de un polinomio de grado
teorema
fundamental n, es todo el plano complejo.
del Álgebra:
todo polinomio Definición 2.4 (Función racional). Sean p(z) y q(z) polinomios. La función r(z) dada
de una variable
no constante
por
con coeficientes p(z)
complejos tiene r(z) =
una raı́z q(z)
compleja
se denomina función racional y está definida en todo número complejo z, excepto donde
q(z) = 0, esto es, el conjunto {z ∈ C : q(z) 6= 0} que no contiene ninguna de las raı́ces del
polinomio q(z).
Solución. Se tiene que el dominio de r(z) es el conjunto de número complejos z tales que
r(z) está bien definida, es decir, el conjunto de números complejos que no produzcan una
división por 0, que en este caso es {z ∈ C : z 6= i}. Ahora, para determinar el rango de r(z)
se utiliza el siguiente procedimiento. Se toma w = r(z) y luego se despeja a z en función
de w. A la función obtenida (la inversa de r(z)), que depende de w, se le determina el
dominio. Este último conjunto de números complejos constituye el rango de r(z). Ası́, se
tiene
z+1
w= ,
z−i
despejando z en función de w, se obtiene
1 − iw
z= ,
w−1
de lo cual se deduce que el rango de r(z) es {w ∈ C : w 6= 1}. ♦
Ejemplo 2.2. Las funciones componentes, u(x, y) y v(x, y), de f (z) = z 2 + a, donde
a = a1 + i a2 , están dadas por
Note que la
2.2.2 Lı́mite definición de
lı́mite de una
Sea f (z) una función de variable compleja definida en un dominio D ⊂ C. Sea z0 un función de
variable
punto de D. compleja es una
extensión del
Definición 2.5 (Lı́mite). Sea f (z) una función definida en todos los puntos de cierta lı́mite de una
vecindad de z0 excepto, posiblemente en el mismo z0 . Se dice que w0 es un lı́mite de función de
variable real.
w = f (z), si para cada número positivo ε, existe un número positivo δ tal que
Denotamos con
lı́m f (z) = w0
z→z0
Algunas de las fórmulas de lı́mites que se estudiaron en el cálculo elemental tienen sus
homólogas en el caso de funciones de variable compleja. Las fórmulas equivalentes, que se
presentan aquı́, se demuestran usando la definición de lı́mite.
Teorema 2.1. Sean f (z) y g(z) funciones de variable complejas tales que lı́m f (z) = w0
z→z0
y lı́m g(z) = W0 . Entonces,
z→z0
w0
3. lı́m [f (z)/g(z)] = .
z→z0 W0
Demostración. Sea ε > 0. Como w0 y W0 son respectivamente los lı́mites de f (z) y g(z)
cuando z tiende a z0 , podemos escribir:
ε
|f (z) − w0 | < , siempre que 0 < |z − z0 | < δ1 ,
2
y
ε
|g(z) − W0 | < , siempre que 0 < |z − z0 | < δ2 ,
2
26 2.2. LÍMITE Y CONTINUIDAD
donde δ1 > 0 y δ2 > 0 dependen de ε. Ası́, por la desigualdad triangular tenemos que
|[f (z) + g(z)] − [w0 + W0 ]| ≤ |f (z) − w0 | + |g(z) − W0 | < ε,
siempre que 0 < |z − z0 | < δ, donde δ ≤ mı́n(δ1 , δ2 ). Por lo tanto, w0 + W0 es el lı́mite de
f (z) + g(z) cuando z tiende a z0 .
Supongamos que
r
ε
|f (z) − w0 | < , siempre que 0 < |z − z0 | < δ3 ,
2
y r
ε
|g(z) − W0 | < , siempre que 0 < |z − z0 | < δ4 ,
2
donde δ3 > 0 y δ4 > 0 dependen de ε. Como W0 (f (z) − w0 ) = W0 f (z) − W0 w0 y w0 (g(z) −
W0 ) = w0 g(z) − w0 W0 , entonces es claro que lı́m W0 f (z) = W0 w0 y lı́m w0 g(z) = w0 W0 ;
z→z0 z→z0
luego, por 1., lı́m (W0 f (z) − w0 g(z)) = 0. Ahora, tenemos que
z→z0
si, y sólo si
lı́m u(x, y) = u0 y lı́m v(x, y) = v0 .
(x,y)→(x0 ,y0 ) (x,y)→(x0 ,y0 )
Ahora, si
lı́m u(x, y) = lı́m f1 (z) = u0 y lı́m v(x, y) = lı́m f2 (z) = v0 ,
(x,y)→(x0 ,y0 ) z→z0 (x,y)→(x0 ,y0 ) z→z0
2.2.3 Continuidad
La continuidad es uno de los conceptos más importantes del análisis. A continuación
veremos que la continuidad de una función de variable compleja se puede ver como una
extensión del concepto de continuidad de una función de variable real.
Definición 2.6. Se dice que una función f (z) es continua en z0 , si satisface las dos
condiciones siguientes:
(z − 1)(z + 1)
f (z) = = z + 1 = (x + 1) + i y,
(z − 1)
luego, las funciones componentes de f (z) son u(x, y) = x + 1 y v(x, y) = y. Como el valor
de f (z0 ) está bien definido para todo z0 = x0 + i y0 6= 1, entonces por el Teorema 2.2
podemos escribir:
El siguiente teorema dice que las funciones definidas como la suma, multiplicación o di-
visión de funciones continuas, son también funciones continuas en su dominio de definición.
La demostración de este teorema se obtiene utilizando el Teorema 2.1.
Teorema 2.3. Sean f (z) y g(z) dos funciones de variable compleja continuas en un do-
minio D. Entonces:
Observación 2.1. El Teorema 2.3 nos permite asegurar que un polinomio es una función
continua en todo el plano complejo; también nos permite aseverar que una función racional
es continua en todo el plano complejo excepto donde su denominador se anule.
En el cálculo elemental vimos que la composición de funciones de variable real continuas
produce una nueva función continua. En el caso de las funciones de variable compleja
también se satisface este hecho, es decir, la composición de funciones de variable compleja
continuas es continua. Esta afirmación se plantea en el siguiente teorema, cuya demostra-
ción se obtiene directamente de la definición de continuidad.
Teorema 2.4. Sean f (z) y g(z) dos funciones de variable compleja definidas respectiva-
mente en los dominios D y E, tales que f (D) ⊆ E. Si f es continua en D y g es continua
en f (D), entonces la función h(z) = g(f (z)) es continua en D.
Por otra parte, dado que el lı́mite de una función de variable compleja se puede calcular
a través del lı́mite de sus funciones componentes (Teorema 2.2), es natural inferir que la
continuidad de una función de variable compleja se corresponde con la continuidad de
sus funciones componentes. Seguidamente se muestra un teorema que trata este aspecto,
cuya demostración es inmediata utilizando la definición de continuidad de funciones reales
conjuntamente con el Teorema 2.2.
Teorema 2.5. Sea f (z) = u(x, y) + i v(x, y) una función de variable compleja. Entonces,
f es continua en un punto z0 = x0 + i y0 si, y sólo si sus funciones componentes, u(x, y)
y v(x, y), son continuas en el punto (x0 , y0 ).
A continuación damos un ejemplo que muestra la aplicación de los Teoremas 2.3 y 2.5,
en el estudio de la continuidad de una función de variable compleja.
Ejemplo 2.5. Estudiemos la continuidad en z = 1 de la función
x − iy
f (z) = 2x + iy + .
x2 + y 2
Verifiquemos, a través de dos procedimientos, la continuidad de f (z). Primero, utilizando
el Teorema 2.3; después, empleando el Teorema 2.5. Escribamos a f (z) en función de z,
para ello consideremos las siguientes expresiones de x e y,
z+z z−z
x= , y= .
2 2i
Sustituyendo convenientemente estas expresiones en la fórmula de definición de f (z), ob-
tenemos:
z+z
z+z z−z + i z−z
f (z) = 2 +i + 2 2i
2 2i zz
z−z z 3 1 1
= z+z+ + = z+ z+ .
2 zz 2 2 z
Definamos tres funciones f1 (z), f2 (z) y f3 (z) como
3 1 1
f1 (z) = z, f2 (z) = z, f3 (z) = .
2 2 z
Las funciones f1 (z) y f3 (z) son, respectivamente, un polinomio de grado uno y una función
racional. Ası́, f1 (z) es continua en todo el plano, particularmente en z = 1, y f3 (z) es
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 29
2.3 Diferenciación
Inmediatamente damos la definición de derivada de una función de variable compleja. ¿Por qué en la
definición de
Definición 2.7 (Derivada). Sea f (z) una función de variable compleja definida en un derivada es
necesario que z0
dominio D. Sea z0 un punto de acumulación de D. La derivada de f (z) en z0 , denotada sea un punto de
d acumulación?
por f (z0 ) o f ′ (z0 ), es una función de variable compleja que se define mediante la ecuación
dz
d f (z0 + h) − f (z0 )
f (z0 ) = f ′ (z0 ) = lı́m , (2.1)
dz h→0 h
Demostración. Como la derivada de f (z) existe en z0 , entonces |f ′ (z0 )| < ∞ y f (z0 ) está
bien definido. Por tanto,
Observación 2.2. El Teorema 2.6 nos garantiza que toda función derivable es continua. Lo
que no es cierto es que toda función continua es derivable; por ejemplo, f (z) = |z|2 es
continua en todo el plano complejo, pero es derivable solamente en z = 0 (se deja al lector
verificar este hecho).
30 2.3. DIFERENCIACIÓN
4. [f1 (z) + f2 (z) + · · · + fn (z)]′ = f1′ (z) + f2′ (z) + · · · + fn′ (z).
8. Si f (z) = a0 + a1 z + a2 z 2 + · · · + am z m , entonces
′
f1 (z) f1′ (z)f2 (z) − f2′ (z)f1 (z)
9. = , siempre y cuando f2 (z) 6= 0.
f2 (z) (f2 (z))2
10. Regla de la Cadena. Sean f (z) derivable en z0 y g(w) derivable en f (z0 ). Entonces
la función h(z) = g(f (z)) es derivable en z0 , y
Note que todas las reglas de diferenciación se cumplen, si las funciones involucradas
son derivables. De esta forma, para aplicar alguna de estas reglas, es necesario primero
verificar que las funciones consideradas sean derivables. En el Cálculo Elemental, por lo
general, era relativamente sencillo verificar si la derivada de una función existı́a o no. Por
ejemplo, observando la gráfica de una función f : R → R se podı́a identificar los puntos
donde existı́a la derivada. Este procedimiento, que por demás es muy sencillo, no se puede
usar para verificar que una función de variable compleja es derivable. En la siguiente
sección damos un procedimiento para estudiar la derivabilidad de las funciones de variable
compleja.
Ejemplo 2.6. Comprobemos que la función f (z) = ex cos y + i ex sen y es derivable en todo
punto z del plano complejo. Las funciones componentes de f (z) están dadas por
∂u ∂u
= ex cos y, = −ex sen y,
∂x ∂y
∂v ∂v
= ex sen y, = ex cos y,
∂x ∂y
las cuales son continuas en todo R2 y, además, satisfacen las ecuaciones de Cauchy-Riemann
en todo z = x + i y. Entonces, por el Teorema 2.8 la función f (z) = ex cos y + i ex sen y es
derivable en todo el plano complejo.
Definición 2.8 (Función analı́tica). Se dice que una función f (z) es analı́tica en z0 ∈ C,
si la derivada de f (z) existe en z0 y en todo punto z de alguna vecindad de z0 . Si f (z) es
analı́tica en todos los puntos de un dominio D ⊂ C, se dice que f (z) es analı́tica en D.
es analı́tica en todo el plano complejo, excepto en las raı́ces del polinomio del denominador,
en otras palabras, salvo en los puntos z tales que b0 + b1 z + · · · + bm z m = 0. Estas raı́ces
son los puntos singulares de r(z).
∂2h ∂2h
∇2 h(x, y) = (x, y) + (x, y) = 0,
∂x2 ∂y 2
El siguiente teorema establece una relación entre las funciones armónicas y las funciones
analı́ticas.
∂2u ∂2v
= . (2.5)
∂x2 ∂x∂y
∂u ∂v
Ahora, derivando con respecto a y la otra ecuación de Cauchy-Riemann ∂y = − ∂x , obte-
nemos:
∂2u ∂2v
= − . (2.6)
∂y 2 ∂y∂x
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 35
∂2u ∂2u
∇2 u(x, y) = (x, y) + (x, y) = 0,
∂x2 ∂y 2
con lo cual se demuestra que u es armónica en D. Utilizando un razonamiento similar se
demuestra que v también es armónica. En consecuencia, u y v son funciones armónicas
que satisfacen las ecuaciones de Cauchy-Riemann en D, es decir, v es armónica conjugada
de u.
(⇐) Supongamos que v es armónica conjugada de u en D y demostremos que la función
f (z) = u(x, y) + i v(x, y) es analı́tica en D. Como v es armónica conjugada de u en D,
entonces u y v son armónicas en D y satisfacen las ecuaciones de Cauchy-Riemann; además,
sus primeras derivadas parciales son continuas en D. Ası́, por el Teorema 2.8, la función
f (z) = u(x, y) + i v(x, y) es analı́tica en D.
Observación 2.4. El Teorema 2.10 nos garantiza que una función f (z) = u(x, y) + i v(x, y)
es analı́tica si v es armónica conjugada de u. Pero, no es cierto que si u y v son armónicas,
entonces f (z) = u(x, y) + i v(x, y) es analı́tica; por ejemplo, si tomamos u(x, y) = x + y y
v(x, y) = x, observamos que u y v son armónicas en todo el plano, pero f (z) = (x + y) + i x
no es analı́tica en ningún punto del plano (se deja al lector demostrar esta afirmación).
Es bueno aclarar que si v es armónica conjugada de u en cierto dominio D, no es en ge-
neral cierto que u sea armónica conjugada de v en ese dominio. Por ejemplo, consideremos
las funciones u(x, y) = x2 − y 2 y v(x, y) = 2xy. Vemos que u y v son las partes real e ima-
Note que
ginaria, respectivamente, de f (z) = z 2 , que es una función entera. Por el Teorema 2.10, v g(z) =
es armónica conjugada de u en todo el plano, pero u no puede ser una armónica conjugada 2xy + i(x2 − y 2 )
de v, ya que la función g(z) = 2xy + i (x2 − y 2 ) no es analı́tica en ningún punto del plano, solo es derivable
en z = 0.
porque sus funciones componentes no satisfacen las ecuaciones de Cauchy-Riemann.
Por otra parte, es cierto que una armónica conjugada, cuando existe, es única, excepto
por una constante aditiva. Expliquemos con un ejemplo un método para obtener una
armónica conjugada de una función armónica dada.
Ejemplo 2.11. Determine v(x, y), la armónica conjugada de u(x, y), si u(x, y) = x + y.
Solución. Para que v sea armónica conjugada de u, la función v debe ser armónica
y, además, se deben satisfacer las ecuaciones de Cauchy-Riemann (2.2). El método para
generar v, una armónica conjugada de u, consiste en forzar que se cumplan las ecuacio-
nes de Cauchy-Riemann y luego resolver ciertas ecuaciones diferenciales. Apliquemos tal
metodologı́a al problema planteado.
Como u(x, y) = x + y, se tiene que u es armónica y
∂u ∂u
= 1, y = 1.
∂x ∂y
Ahora, utilizando una de las ecuaciones de Cauchy-Riemann, digamos
∂u ∂v
= ,
∂x ∂y
36 2.5. FUNCIONES ARMÓNICAS
podemos escribir:
∂v
= 1.
∂y
Integrando con respecto a y esta última ecuación obtenemos:
φ(x) = −x + c,
v(x, y) = y − x + c,
que es armónica (se deja al lector verificar esta afirmación). Es decir, hemos construido
una familia de funciones armónicas conjugadas de u que se diferencian entre sı́ por una
constante. Para obtener una de ellas es necesario un valor de v en algún punto del plano.
Por ejemplo, si v(0, 0) = 1, entonces la constante c correspondiente es c = 1 y la función v
es v(x, y) = y − x + 1. ♦
Ejemplo 2.12. Sea v(x, y) = x. Determine una función analı́tica f (z) = u(x, y) + i v(x, y)
tal que f (0) = −1.
∂u
= 0.
∂x
Integrando con respecto a x esta última ecuación obtenemos:
u(x, y) = φ(y),
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 37
∂u ∂v
donde φ : R → R. Ahora, utilizando la otra ecuación de Cauchy-Riemann, ∂y = − ∂x ,
podemos escribir:
∂u ∂v
φ′ (y) = =− = −1,
∂y ∂x
es decir,
φ′ (y) = −1.
Integrando con respecto a y la ecuación anterior obtenemos:
φ(y) = −y + c,
donde c es una constante real. Por lo tanto, considerando la expresión de φ(y), la función
u está dada por
u(x, y) = −y + c.
Por construcción, v es armónica conjugada de u, luego, por el Teorema 2.10 las funciones
f (z) están dadas por
f (z) = c − y + i x,
las cuales constituyen una familia de funciones enteras que se diferencian entre sı́ por una
constante real c. Ahora, la función deseada debe satisfacer f (0) = −1. Forzando esta
última condición y despejando c, se tiene que c = −1 y la función pedida es
f (z) = −1 − y + i x = iz − 1.
las cuales satisfacen las ecuaciones de Cauchy-Riemann y sus derivadas parciales son con-
tinuas en todo el plano complejo. Luego, la función exponencial es analı́tica en todo el
plano complejo (ver Ejemplo 2.6). Por lo tanto, ez es una función entera, cuya derivada
está dada por:
d z ∂u ∂v
e = +i = ex cos y + i ex sen y = ez ,
dz ∂x ∂x
que coincide con la propiedad deseada de la función exponencial real.
38 2.6. FUNCIONES ELEMENTALES
eiz − e−iz
sen z = (2.10)
2i
y la función coseno complejo como
eiz + e−iz
cos z = . (2.11)
2
Las funciones sen z y cos z son combinaciones lineales de las funciones enteras eiz y e−iz ,
por lo tanto, sen z y cos z son funciones enteras; además, sus derivadas son:
d i(eiz + e−iz )
sen z = = cos z
dz 2i
y
d i(eiz − e−iz )
cos z = = −sen z.
dz 2i
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 39
• Las funciones seno y coseno son periódicas con periodo 2π. (Se deja como ejercicio
para el lector la demostración).
• Los ceros de sen z son todos los números complejos z = nπ, para n = 0, ±1, ±2, . . ..
En efecto, para que z ∈ C sea un cero de sen z, se debe cumplir que sen z = 0, luego,
utilizando la ecuación (2.12) podemos escribir:
sen x cosh y = 0, cos x senh y = 0.
Ahora, utilizando la ecuación sen x cosh y = 0, se deduce que y = 0. Sustituyendo
este valor de y en la ecuación sen x cosh y = 0, obtenemos sen x = 0, de lo cual se
tiene que x = nπ, para n = 0, ±1, ±2, . . ... En consecuencia, los ceros de sen z son
todos los números complejos z = nπ, para n = 0, ±1, ±2, . . ..
• Los ceros de cos z son todos los números complejos z = (n + 21 )π, para n = 0, ±1, . . ..
(Se deja al lector verificar este hecho).
40 2.6. FUNCIONES ELEMENTALES
Tangente
sen z
tan z =
cos z
Cotangente
cos z
cot z =
sen z
Secante
1
sec z =
cos z
Cosecante
1
csc z =
sen z
Las funciones tan z, cot z, sec z y csc z, son analı́ticas en todos los números complejos z
donde no se anule el denominador de la expresión que las define, además, la derivada de
cada una de ellas es, respectivamente,
d
tan z = sec2 z,
dz
d
cot z = − csc2 z,
dz
d
sec z = tan z sec z,
dz
d
csc z = − cot z csc z.
dz
ez − e−z
senh z = (2.14)
2
y la función coseno hiperbólico complejo como
ez + e−z
cosh z = . (2.15)
2
Las funciones senh z y cosh z son combinaciones lineales de las funciones enteras ez y e−z ,
por lo tanto, senh z y cosh z son funciones enteras; además, sus derivadas son:
d
senh z = cosh z
dz
y
d
cosh z = senh z.
dz
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 41
• Las identidades hiperbólicas que satisfacen el seno y coseno hiperbólicos reales también
son válidas en el caso complejo; por ejemplo:
cosh2 z − senh2 z = 1,
senh(z1 + z2 ) = senh z1 cosh z2 + cosh z1 senh z2 ,
cosh(z1 + z2 ) = cosh z1 cosh z2 + senh z1 senh z2 ,
entre otras. (Se deja como ejercicio para el lector la comprobación de cada una de
las identidades).
y
cosh z = cosh x cos y + i senh x sen y. (2.17)
• Las funciones seno y coseno hiperbólicos son periódicas con periodo 2πi.
• Los ceros de senh z son todos los números complejos z = nπi, n = 0, ±1, ±2, . . .
• Los ceros de cosh z son todos los números complejos z = (n + 21 )πi, n = 0, ±1, ±2, . . .
Tangente hiperbólica
senh z
tanh z =
cosh z
Cotangente hiperbólica
cosh z
coth z =
senh z
Secante hiperbólica
1
sech z =
cosh z
Cosecante hiperbólica
1
csch z =
senh z
42 2.6. FUNCIONES ELEMENTALES
Las funciones tanh z, coth z, sech z y csch z, son analı́ticas en todos los números complejos
z donde no se anule el denominador de la expresión que las define, además, la derivada de
cada una de ellas es, respectivamente,
d
tanh z = sech 2 z,
dz
d
coth z = −csch 2 z,
dz
d
sech z = − tanh z sech z,
dz
d
csch z = − coth z csch z.
dz
donde ln · denota el logaritmo natural. Veamos que log z es una función multivaluada. Se
tiene que, para todo z 6= 0, el argumento de z se puede escribir como:
Observe que para cualquier z 6= 0, los valores de log z tienen la misma parte real y sus
partes imaginarias difieren en múltiplos enteros de 2π. Por lo tanto, log z es una función
multivaluada. El siguiente ejemplo ilustra este hecho.
π
Ejemplo 2.13. Calculemos todos los valores de log(1 + i). Se tiene que Arg (1 + i) = y
√ 4
|1 + i| = 2. Ası́,
√ π
log(1 + i) = ln 2 + i + 2nπ (n = 0, ±1, ±2, . . .).
4
Definición 2.13 (Valor principal del logaritmo). El valor principal de log z es el valor que
se obtiene de la fórmula (2.18) cuando se utiliza el argumento principal de z. Ese valor se
denota como Log z y se da por medio de la ecuación
• La función Log z es continua en el dominio (|z| > 0, −π < arg z < π). (Se deja al
lector la prueba).
• La función Log z es analı́tica en el dominio (|z| > 0, −π < arg z < π), y su derivada
está dada por
d 1
Log z = .
dz z
(Se deja al lector la verificación).
Definición 2.14 (Rama de una función multivaluada). Una rama de una función multi-
valuada f (z), es una función monovaluada F (z) que es analı́tica en cierto dominio D ⊂ C
y que coincide con f (z) en D, es decir, F (z) = f (z) para todo z ∈ D.
Definición 2.15 (Corte y punto ramal). Un corte ramal es una lı́nea o curva de puntos
singulares que se introducen al definir una rama de una función multivaluada. El punto
común a todos los cortes ramales de una función multivaluada se denomina punto ramal.
Observe que Log z es una rama de la función multivaluada log z, además, se tiene que
Log z = log z, para todo z 6= 0 tal que |z| > 0 y −π < arg z < π. La función Log z se
denomina rama principal del logaritmo. El corte ramal de Log z es el eje real negativo con
el origen (ver Figura 2.2).
donde α es un número real fijo. El corte ramal de esta rama es el rayo arg z = α (ver
Figura 2.3).
al
punto ramal r am
e
c ort
α
b
Ejemplo 2.14. Determinemos el valor de log(1 + i), si log z está definido por
log z = ln |z| + i arg z (|z| > 0, π/2 < arg z < 5π/2).
√ iπ/4
Se tiene que 1+i = 2e , pero π/4 6∈ (π/2, 5π/2). Luego, el argumento principal de 1+i
no es el valor indicado para calcular log(1 + i). Para encontrar arg (1 + i) ∈ (π/2, 5π/2),
se utiliza el argumento principal de 1 + i. Se tiene que el valor del argumento de 1 + i
buscado está dado por:
π 9π
+ 2π = ∈ (π/2, 5π/2);
4 4
por lo tanto, arg (1 + i) = 9π/4 es el valor indicado para calcular log(1 + i). Ası́,
√ 9π
log(1 + i) = ln 2+i .
4
z c = ec log z , (2.20)
donde |z| > 0 y α < arg z < α + 2π, con α un número real fijo. En otras palabras, las
ramas de la función exponente complejo se definen según la rama del logaritmo que se esté
utilizando. De esta forma, cuando α = −π, estaremos empleando la rama principal del
logaritmo para definir la rama principal de la función exponente complejo. Para obtener
la derivada de la función exponente complejo, se emplea la regla de la cadena. Ası́, la
derivada de cada una de las ramas de la función exponente complejo está dada por
d c
z = c z c−1 .
dz
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 45
para todo z ∈ C. De igual manera que la función exponente complejo, la función expo-
nencial de base c es generalmente multivaluada, cuyas ramas se definen según la rama
del logaritmo que se esté empleando; por ejemplo, si utilizamos el valor principal de log c,
entonces la rama que se obtiene es la rama principal de la función exponencial de base c.
Para cada rama de la función exponencial de base c, su derivada está dada por:
d z
c = cz log c.
dz
eiw − e−iw
z=
2i
o, equivalentemente,
e2iw − 2iz eiw − 1 = 0.
eiw = iz + (1 − z 2 )1/2 ,
donde (1−z 2 )1/2 es una función bivaluada de z. Si tomamos logaritmo en ambos miembros
y recordamos que w = sen −1 z, se llega a la siguiente fórmula cerrada de la función sen −1 z
sen −1 z = −i log iz + (1 − z 2 )1/2 . (2.22)
Las funciones inversas del coseno y la tangente se pueden obtener fácilmente realizando
un procedimiento similar al antes descrito. Las expresiones de cos−1 z y tan−1 z son,
respectivamente:
cos−1 z = −i log z + i(1 − z 2 )1/2 , (2.23)
i i+z
tan−1 z = log . (2.24)
2 i−z
En general, las funciones trigonométricas inversas son multivaluadas. Por ejemplo, para
la función sen −1 z podemos elegir dos valores igualmente válidos para la raı́z cuadrada que
aparece en su definición. Una vez que hemos elegido una valor, existe un número infinito
de valores posibles para el logaritmo de iz + (1 − z 2 )1/2 . En resumen, vemos que, debido a
la raı́z cuadrada y al logaritmo, hay dos conjuntos distintos de valores de sen −1 z, cada uno
de los cuales posee un número infinito de elementos. Este hecho sucede de forma semejante
para las restantes funciones trigonométricas inversas.
Las derivadas de las funciones trigonométricas inversas se obtienen directamente de sus
expresiones cerradas. Las derivadas de las funciones inversa del seno e inversa del coseno
dependen de los valores escogidos de la raı́ces cuadradas, a saber:
d 1
sen −1 z = ,
dz (1 − z 2 )1/2
d −1
cos−1 z = .
dz (1 − z 2 )1/2
En cambio, la derivada de la inversa de la tangente,
d 1
tan−1 z = ,
dz 1 + z2
no depende de la manera en que la función se haga monovaluada.
Ejercicio 2.3. Determinar las expresiones cerradas de las derivadas de sec−1 z, csc−1 z y
cot−1 z.
El siguiente ejemplo ilustra la aplicación de las expresiones de las funciones inversas en
la resolución de ecuaciones que involucran funciones trigonométricas.
Ejemplo 2.17. Resolver la ecuación
sen 2 z + 2i sen z − 1 = 0.
Solución. La ecuación dada se puede escribir equivalentemente como
(sen z + i)2 = 0,
la cual es cierta si, y sólo si la siguiente ecuación se cumple
sen z = −i.
Al tomar la inversa del seno en ambos miembros de la ecuación anterior, podemos escribir:
z = sen −1 (−i) = −i log i(−i) + (1 − (−i)2 )1/2
√
= −i log 1 + 21/2 = −i log 1 ± 2 .
Según el valor de la raı́z cuadrada empleado, obtenemos dos conjuntos de valores de z que
resuelven la ecuación dada. Ası́, uno de estos conjuntos de números complejos está dado
por: √ √
zk = −i log 1 + 2 = 2kπ − i ln 1 + 2 , para k = 0, ±1, ±2, . . .
El otro conjunto es:
√
√
zn = −i log 1 − 2 = (2n + 1)π − i ln 1 − 2 , para n = 0, ±1, ±2, . . .
♦
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 47
Las derivadas de las funciones senh−1 z y cosh−1 z, dependen de los valores escogidos
de las raı́ces cuadradas, a saber:
d 1
senh−1 z = ,
dz (z 2+ 1)1/2
d −1
cosh−1 z = .
dz (z − 1)1/2
2
2.7 Mapeos
En esta sección se estudia la transformación de regiones del plano complejo a través
de funciones de variable compleja, tales como polinomios de grado 1 y funciones racionales
obtenidas como el cociente de dos polinomios de grado 1. Recordemos que en el cálculo
elemental se podı́a obtener la gráfica de una función real y = f (x) y con ello estudiar su
comportamiento; en cambio para una función w = f (z) de una variable compleja, no es
tan fácil elaborar su gráfica. Para ello se requieren dos números x e y para definir un
valor z cualquiera y otros dos números para los valores de u y v correspondientes. Por lo
tanto, se requiere un espacio de cuatro dimensiones para representar w = f (z) en forma
gráfica. Evidentemente una gráfica de cuatro dimensiones no es un medio conveniente para
estudiar el efecto gráfico de una función de variable compleja. Es preciso recurrir a otros
medios para visualizar w = f (z).
Aquı́, visualizamos la relación w = f (z) como el efecto geométrico que tiene la función
f (z) sobre un conjunto de puntos complejos. A este proceso lo denominamos mapear o
transformar y a la función f (z) de variable compleja la denominamos transformación o
mapeo. El objetivo principal de esta sección es determinar analı́tica y gráficamente el efecto
que tiene un mapeo sobre un conjunto de puntos del plano complejo; en particular, estudia-
mos el efecto de los mapeos lineales, inversión y bilineales, sobre rectas y circunferencias.
Para esta parte necesitamos la siguiente notación:
48 2.7. MAPEOS
z Mapeo o w
b b
Transformación
x u
Plano z Plano w
2.7.1 Mapeo w = z + c
El mapeo del plano z en el plano w definido por la ecuación
w = z + c, (2.25)
b
Traslación
w0 = z0 + c
x u
b
c
El mapeo (2.25) es inyectivo; por tanto, posee mapeo inverso definido como
f −1 (w) = w − c. (2.26)
x + i y = f −1 (u + i v),
α(x2 + y 2 ) + βx + δy + γ = 0, (2.27)
x + i y = f −1 (u + i v) = u + i v − c1 + i c2 = (u − c1 ) + i (v − c2 ),
3
S
2
x
1 2 3
Ejemplo 2.18. Sea S el conjunto de números complejos que se muestra en la Figura 2.6.
Determinar el transformado de S bajo el mapeo w = f (z) = z − i.
2 f (S)
u
1 2 3
♦
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 51
2.7.2 Mapeo w = bz
El mapeo del plano z en el plano w definido por la ecuación
w = bz, (2.28)
z = reiθ
b
θ
b
b = leiβ θ+β
β
x u
El mapeo (2.28) es inyectivo; por tanto, posee mapeo inverso definido como
w
f −1 (w) = . (2.29)
b
Además, el mapeo (2.28) transforma rectas a rectas y circunferencias a circunferencias.
Ejercicio 2.5. Demuestre que el mapeo lineal (2.28) transforma rectas a rectas y cir-
cunferencias a circunferencias. Ayuda: utilice la ecuación general de una recta o una
circunferencia (2.27) conjuntamente con la expresión del mapeo inverso (2.29).
Ejemplo 2.19. Sea S el conjunto de números complejos que se muestra en la Figura 2.9.
Determinar el transformado de S bajo el mapeo w = f (z) = (1 − i)z.
y
2 S
1 b
x
1 2
v
2 f (S)
1 √
2
b
u
1 2 3
−1
2.7.3 Mapeo w = bz + c
El mapeo del plano z en el plano w definido por la ecuación
w = bz + c, (2.30)
donde b y c son números complejos, es una rotación en el ángulo Arg b y una expansión ó
contracción según sea el valor de |b|, seguida de una traslación en la dirección del vector c.
En efecto, el mapeo (2.30) se puede expresar a través de los siguientes mapeos sucesivos:
z / Z = bz / w =Z+c
Rotación y Traslación
Expansión ó
Contracción
w−c
f −1 (w) = , b 6= 0. (2.31)
b
Además, como (2.30) es la composición de una rotación y una traslación, entonces el mapeo
w = bz + c transforma rectas a rectas y circunferencias a circunferencias.
Ejemplo 2.20. Sea S el conjunto de números complejos que se muestra en la Figura 2.11.
Determinar el transformado de S bajo el mapeo w = f (z) = iz + 3 − i.
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 53
3
S
2
1 b
x
1 2 3
Solución. Se tiene que el mapeo inverso es f −1 (w) = −iw + 1 + 3i, además, el conjunto
S se puede escribir como S = D1 ∪ D2 , donde D1 y D2 son los conjuntos de números
complejos definidos respectivamente por:
D1 = {z = x + i y : 2x − y ≥ 1, 2x + y ≤ 7, y ≥ 1}
y
D2 = {z = x + i y : (x − 2)2 + (y − 1)2 ≤ 1, y ≤ 1}.
De esta forma, como el mapeo es inyectivo, se tiene que el transformado de S está dado
por f (S) = f (D1 ) ∪ f (D2 ). Ahora, utilizando la ecuación x + i y = f −1 (u + i y), tenemos:
x = v + 1, y = 3 − u. Sustituyendo los valores de x e y en las inecuaciones que describen
a los conjuntos D1 y D2 , obtenemos que los conjuntos transformados de D1 y D2 , bajo el
mapeo w = iz + 3 − i, están dados por
f (D1 ) = {w = u + i v : u + 2v ≥ 2, −u + 2v ≤ 2, u ≤ 2}
y
f (D2 ) = {w = u + i v : (u − 2)2 + (v − 1)2 ≤ 1, u ≥ 2}.
Ası́, la unión de f (D1 ) y f (D2 ) conforma el transformado de S bajo el mapeo w = iz +3−i,
cuya representación gráfica se aprecia en la Figura 2.12.
f (S)
2
1 b
u
1 2 3
♦
54 2.7. MAPEOS
1
w= , (2.32)
z
para todo z 6= 0, se denomina mapeo inversión, y establece una correspondencia uno a uno
entre los puntos de los planos z y w distintos de cero. Por lo tanto, el mapeo inverso de
f (z) = 1/z es f −1 (w) = 1/w, es decir, el mismo mapeo inversión.
Por otra parte, el mapeo (2.32) también puede escribirse como
z
w= , para todo z 6= 0.
|z|2
Por ello, el mapeo inversión se puede expresar a través de los siguientes mapeos sucesivos:
z /Z= 1 z /w=Z
2 |z|
γ(u2 + v 2 ) + βu − δv + α = 0. (2.34)
y 1 v
w=
z
x u
y 1 v
w=
z
x u
y 1 v
w=
z
x u
y 1 v
w=
z
x u
Ejercicio 2.6. Compruebe cada uno de los conjuntos transformados de los Ejemplos 2.21-
2.24.
Ejemplo 2.25. Sea S el conjunto de números complejos dado en el Ejemplo 2.18 y que se
muestra en la Figura 2.6. Comprobemos que el transformado de S bajo el mapeo inversión
es el conjunto de números complejos que se muestra en la Figura 2.13. Se tiene que
S = {z = x + i y : 2x − y ≥ 1, 2x + y ≤ 7, y ≥ 1},
además, la función inversa de f (z) = 1/z es f −1 (w) = 1/w. Ahora, utilizando la ecuación
x + i y = f −1 (u + i v) obtenemos: x = u/(u2 + v 2 ), e y = −v/(u2 + v 2 ). Sustituyendo
convenientemente estas últimas expresiones en el conjunto de inecuaciones que describe
a S, se tiene que los puntos w = u + i v del transformado de S bajo el mapeo inversión
satisfacen las siguientes inecuaciones
1 2 5 1 2 1 2 5 1 2 1
(u − 1)2 + v − ≤ , u− + v+ ≥ , u2 + v + ≤ .
2 4 7 14 196 2 4
1
b
u
b
1 2
−1
f (S)
az + b
w= , (2.35)
cz + d
Contracción, y Contracción, y
Traslación Traslación
En otras palabras, el mapeo (2.35) es la composición de mapeos lineales, por lo cual él
transforma circunferencias o rectas en el plano z a circunferencias o rectas en el plano w.
Por otra parte, cuando el determinante del mapeo es cero, ad − bc = 0, el mapeo (2.35)
adquiere la forma
a
w= ,
c
es decir, el mapeo (2.35) transforma todo el plano complejo en el punto w = a/c.
Ejemplo 2.26. Sea S el conjunto de números complejos dado en el Ejemplo 2.18 y que se
muestra en la Figura 2.6. Determinar el transformado de S bajo el mapeo bilineal
iz + 1 + i
w = f (z) = .
iz + i
Solución. Se tiene que
S = {z = x + i y : 2x − y ≥ 1, 2x + y ≤ 7, y ≥ 1}.
b
b
b u
1 2
−1
f (S)
>> syms x y
Los objetos simbólicos se pueden manipular siguiendo las reglas de las operaciones ma-
temáticas, por ejemplo:
>> x + y - 5* x + y ^2
ans =
y ^2 + y - 4* x
Una función de variable compleja f (z) también se puede crear como una función
simbólica. Primero se crea el objeto simbólico z, luego se define la función simbólica f(z)
según su forma algebraica. En el siguiente ejemplo se crea la función f (z) = z 2 + iz − 1.
>> syms z
>> f ( z ) = z ^2 + i * z - 1
f(z) =
z ^2 + z * i - 1
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 59
>> f ( i )
ans =
-3
>> syms z
>> f1 ( z ) = z ^2 - i , f2 ( z ) = (z -1) *(z -2)
f1 ( z ) =
z ^2 - i
f2 ( z ) =
( z - 1) *( z - 2)
>> g1 ( z ) = f1 ( z ) + f2 ( z ) , g2 ( z ) = f1 ( z ) * f2 ( z ) , g3 ( z ) = f1 ( z ) / f2 ( z )
g1 ( z ) =
( z - 1) *( z - 2) + z ^2 - i
g2 ( z ) =
( z ^2 - i ) *( z - 1) *( z - 2)
g3 ( z ) =
Para más
( z ^2 - i ) /(( z - 1) *( z - 2) ) detalles del
comando limit,
teclee help
Una operación importante en matemáticas es el cálculo de lı́mites, la cual también limit en el
se puede realizar con la herramienta simbólica de MATLAB. En particular, el comando escritorio de
MATLAB
limit(f,z,a) calcula el lı́mite lı́mz→a f (z). Observe el siguiente ejemplo donde se calcula
el lı́mite
z + z + 2 Re z (Im z + 1)i
lı́m = 1 + i.
z→0 z+z
>> syms z
>> f ( z ) = ( z + conj ( z ) +2* real ( z ) *( imag ( z ) +1) * i ) /( z + conj ( z ) )
f(z) =
( z + conj ( z ) + real ( z ) *( imag ( z ) + 1) *2* i ) /( z + conj ( z ) )
>> limit (f ,z ,0)
ans ( z ) =
1 + i
>> syms z
>> f ( z ) =(2* i - z ^4) /( z ^3+ i * z ^2 - z +1)
f(z) =
-( z ^4 - 2* i ) /( z ^3 + z ^2* i - z + 1)
>> diff (f , z )
ans ( z ) =
-(4* z ^3) /( z ^3+ z ^2*i - z +1) +(( z ^4 -2* i ) *(3* z ^2+ z *(2* i ) -1) ) /
( -z + z ^2* i + z ^3+1) ^2
>> diff (f ,z ,2)
ans ( z ) =
-(12* z ^2) /( z ^3+ z ^2*i - z +1) -(2*( z ^4 -2* i ) *(3* z ^2+ z *(2* i ) -1) ^2) /
( -z + z ^2* i + z ^3+1) ^3+(8* z ^3*(3* z ^2+ z *(2* i ) -1) ) /
( -z + z ^2* i + z ^3+1) ^ 2 + ( ( 6 *z + 2* i ) *( z ^4 -2* i ) ) /
( -z + z ^2* i + z ^3+1) ^2
>> diff (f ,z ,3)
ans ( z ) =
(6*( z ^4 -2* i ) ) /( - z + z ^2* i + z ^3+1) ^2 -(24*z ) /( z ^3+ z ^2*i - z + ) +
(6*( z ^4 -2* i ) *(3* z ^2+ z *(2* i ) -1) ^3) /( - z + z ^2* i + z ^3+1) ^4+
(12* z ^3*(6* z +2* i ) ) /( - z + z ^2* i + z ^3+1) ^ 2 + ( 3 6 *z ^2*(3* z ^2+
z *(2* i -1) ) /( - z + z ^2* i + z ^3+1) ^2 -(24*z ^3*(3* z ^2+ z *(2* i ) -1) ^2) /
( -z + z ^2* i + z ^3+1) ^3 -(6*(6* z +2* i ) *( z ^4 -2* i ) *
(3* z ^2+ z *(2* i ) -1) ) /( - z + z ^2* i + z ^3+1) ^3
En este caso, las expresiones obtenidas de la primera, segunda y tercera derivada de f (z)
no están simplificadas. Para buscar la forma más simple de una expresión simbólica se
utiliza el comando simple. Utilicemos tal comando para hallar las formas simplificadas de
la primera, segunda y tercera derivada de f (z) en el ejemplo anterior.
>> cos (1 - i )
ans =
0.8337 + 0.9889 i
>> log (1+ i )
ans =
0.3466 + 0.7854 i
>> exp ( pi * i )
ans =
-1.0000 + 0.0000 i
Con el comando cos(1-i) √ se obtiene el valor de cos(1 − i), con log(1+i) se obtiene el
valor de Log (1 + i) = ln 2 + i π4 , y con exp(pi*i) se obtiene eπi = −1.
Con las funciones elementales se pueden crear cualquier tipo de funciones simbólicas.
Observe el siguiente ejemplo donde se crea la función
sen (1 − Log z)
f (z) =
ez 2 +1 − 2
>> syms z
>> f ( z ) = sin (1 - log ( z ) ) /( exp ( z ^2+1) -2)
f(z) =
- sin ( log ( z ) - 1) /( exp ( z ^2 + 1) - 2)
>> abs ( f ( i ) )
ans =
( cosh ( pi /2) ^2* sin (1) ^2 + sinh ( pi /2) ^2* cos (1) ^2) ^(1/2)
62 2.8. FUNCIONES, LÍMITE, DERIVADA Y MAPEOS CON MATLAB
MATLAB también puede emplearse para hallar los conjuntos transformados de circun-
ferencias y polı́gonos, bajo el mapeo lineal w = bz + c. El Listado 2.1 muestra la función
MapeoPoligono que halla el transformado de un polı́gono con vértices z1 , z2 , . . . , zn ∈ C,
bajo el mapeo lineal w = bz + c.
f u n c t i o n w = M a p e o P o l i g o n o (b ,c , z )
% M a p e o P o l i g o n o Halla halla el t r a n s f o r m a d o de un p o l ı́ g o n o
% cuyos v é r t i c e s son los e l e m e n t o s del a r r e g l o z ,
% bajo el mapeo w = bz + c .
% w es un arrego que c o n t i e n e los v é r t i c e s del
% p o l ı́ g o n o t r a n s f o r m a d o
w = b*z + c;
s u b p l o t (1 ,2 ,1)
plot ([ z , z (1) ] , ’r ’ , ’ LineWidth ’ ,2)
xlabel ( ’x ’)
ylabel ( ’y ’)
grid on
xmin = min ([ min ( real ( z ) ) min ( real ( w ) ) ]) -.5;
xmax = max ([ max ( real ( z ) ) max ( real ( w ) ) ]) +.5;
ymin = min ([ min ( imag ( z ) ) min ( imag ( w ) ) ]) -.5;
ymax = max ([ max ( imag ( z ) ) max ( imag ( w ) ) ]) +.5;
axis ([ xmin xmax ymin ymax ])
axis equal
title ( ’ Polı́gono ’)
%
s u b p l o t (1 ,2 ,2)
plot ([ w , w (1) ] , ’b ’ , ’ LineWidth ’ ,2)
xlabel ( ’u ’)
ylabel ( ’v ’)
grid on
axis ([ xmin xmax ymin ymax ])
axis equal
texto = s p r i n t f( ’w =( % s ) z +% s ’ , n u m 2 s t r( b ) , n u m 2 s t r( c ) ) ;
texto = [ ’ P o l ı́ g o n o t r a n s f o r m a d o bajo el mapeo ’, texto ];
title ( texto )
end
5 5
4 4
3 3
2 2
y
v
1 1
0 0
−1 −1
−2 −2
−2 0 2 −2 −1 0 1 2 3
x u
4 4
3 3
2 2
1 1
0 0
y
−1 −1
−2 −2
−3 −3
−4 −4
−2 0 2 4 −2 0 2 4
x u
64 2.9. PROBLEMAS RESUELTOS
f u n c t i o n [ w0 , l ] = M a p e o C i r c u n f e r e n c i a (b ,c , z0 , r )
% M a p e o C i r c u n f e r e n c i a Halla el t r a n s f o r m a d o de una circunfe -
% rencia con centro z0 y radio r , bajo el
% mapeo w = bz + c .
% w0 es el centro y l es el radio de la c i r c u n f e r e n c i a
% transformada
w0 = b * z0 + c ;
l = abs ( b ) * r ;
t = l i n s p a c e (0 ,2* pi ) ;
s u b p l o t (1 ,2 ,1)
xt = r * cos ( t ) + real ( z0 ) ;
yt = r * sin ( t ) + imag ( z0 ) ;
plot ( xt , yt , ’r ’ , ’ LineWidth ’ ,2)
xlabel ( ’x ’)
ylabel ( ’y ’)
grid on
rmax = max ([ r l ]) ;
xmin = min ([ real ( z0 ) real ( w0 ) ]) - rmax -.2;
xmax = max ([ real ( z0 ) real ( w0 ) ]) + rmax +.2;
ymin = min ([ imag ( z0 ) imag ( w0 ) ]) - rmax -.2;
ymax = max ([ imag ( z0 ) imag ( w0 ) ]) + rmax +.2;
axis ([ xmin xmax ymin ymax ])
axis equal
title ( ’ C r i c u n f er encia ’)
%
s u b p l o t (1 ,2 ,2)
ut = l * cos ( t ) + real ( w0 ) ;
vt = l * sin ( t ) + imag ( w0 ) ;
plot ( ut , vt , ’b ’ , ’ LineWidth ’ ,2)
xlabel ( ’u ’)
ylabel ( ’v ’)
grid on
axis ([ xmin xmax ymin ymax ])
axis equal
texto = s p r i n t f( ’w =( % s ) z +% s ’ , n u m 2 s t r( b ) , n u m 2 s t r( c ) ) ;
texto = [ ’ C i r c u n f e r e n c i a t r a n s f o r m a d a bajo el mapeo ’, texto ];
title ( texto )
end
Solución. Demostremos que los lı́mites dados alcanzan valores distintos cuando z tiende
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 65
a z0 por los siguientes caminos: i) una recta horizontal que pasa por z0 , esto es, tomando
z = x + iy0 ; y ii) una recta vertical que pasa por z0 , esto es, tomando z = x0 + iy.
a) El lı́mite de z−z
z−z0 cuando z tiende a z0 por el camino i) está dado por:
0
z − z0 x + iy0 − x0 + iy0 x − x0
lı́m = lı́m = lı́m = 1.
z→z0 z − z0 z→z 0 x + iy0 − x0 − iy0 z→z 0 x − x0
z−z0
Ahora, el lı́mite de z−z0 cuando z tiende a z0 por el camino ii) está dado por:
z − z0 x0 + iy − x0 + iy0 i(−y + y0 )
lı́m = lı́m = lı́m = −1.
z→z0 z − z0 z→z0 x0 + iy − x0 − iy0 z→z0 i(y − y0 )
z−z0
Por lo tanto, el lı́mite lı́m no existe.
z→z0 z−z0
Re z−Re z0
b) El lı́mite de z−z0 cuando z tiende a z0 por el camino i) está dado por:
Re z−Re z0
El lı́mite de z−z0 cuando z tiende a z0 por el camino ii) está dado por:
Re z−Re z0
Por lo tanto, el lı́mite lı́m z−z0 no existe. ♦
z→z0
Problema 2.2. Sea f (z) una función de variable compleja definida como
2
z ,
Im z < Re z,
2
f (z) = 2(Re z) i, Im z = Re z,
(Im z)z, Im z > Re z.
Determine el conjunto de todos los números complejos z donde la función f (z) es continua.
Solución. Es claro que f (z) está bien definida en todo z ∈ C. Como f (z) = z 2 para
Im z < Re z y z 2 es un polinomio, entonces f (z) es continua en z ∈ C tal que Im z < Re z.
Ahora, para z = x + iy tal que Im z > Re z, se tiene que f (z) = (Im z)z = yx + iy 2 ; de lo
cual se deduce que las funciones componentes de f (z) son u(x, y) = yx y v(x, y) = y 2 , para
y > x, que son continuas en todo (x, y) ∈ R2 , en particular, en todo (x, y) tal que y > x.
De esta forma, f (z) es continua en todo z tal que Im z > Re z o Im z < Re z. Veamos ahora
que f (z) no es continua en z tal que Im z = Re z 6= 0. Sea z0 tal que Im z0 = Re z0 6= 0.
Calculemos el lı́mite lı́m f (z) en la dirección de la recta Im z = Re z0 − Re z + Im y0 , con
z→z0
Im z > Re z. Ası́,
lı́m f (z) = lı́m (Im z)z = (Im z0 )z0 = (Re z0 )2 + i(Re z0 )2 6= 2(Re z0 )2 i = f (z0 ),
z→z0 z→z0
Por lo cual f (z) no es continua en z tal que Im z = Re z 6= 0. Es claro que f (z) es continua
z0 = 0. De todo lo anterior, se tiene que el conjunto de todos los números complejos z
donde f (z) es continua está dado por:
C − {z ∈ C : Im z = Re z 6= 0}
♦
66 2.9. PROBLEMAS RESUELTOS
Problema 2.3. Si u y v se expresan en términos de las coordenadas polares (r, θ), muestre
que las ecuaciones de Cauchy-Riemann pueden escribirse de la forma, para r > 0,
∂u 1 ∂v
=
∂r r ∂θ
1 ∂u = − ∂v
r ∂θ ∂r
denominada forma polar de las ecuaciones de Cauchy-Riemann.
∂u ∂v
= −
∂y ∂x
Ahora, haciendo el cambio de variable x = r cos θ e y = r sen θ y aplicando la regla de la
cadena, obtenemos:
∂u ∂u ∂x ∂u ∂y ∂u ∂u
= + = cos θ + sen θ,
∂r ∂x ∂r ∂y ∂r ∂x ∂y
∂u ∂u ∂x ∂u ∂y ∂u ∂u
= + = −r sen θ + r cos θ,
∂θ ∂x ∂θ ∂y ∂θ ∂x ∂y
∂v ∂v ∂x ∂v ∂y ∂v ∂v
= + = cos θ + sen θ,
∂r ∂x ∂r ∂y ∂r ∂x ∂y
∂v ∂v ∂x ∂v ∂y ∂v ∂v
= + = −r sen θ + r cos θ.
∂θ ∂x ∂θ ∂y ∂θ ∂x ∂y
De esta forma, usando las ecuaciones anteriores conjuntamente con las ecuaciones de
Cauchy-Riemann en su forma cartesiana, podemos escribir:
∂u 1 ∂v ∂u ∂u 1 ∂v ∂v
− = cos θ + sen θ − −r sen θ + r cos θ
∂r r ∂θ ∂x ∂y r ∂x ∂y
∂u ∂v ∂u ∂v ∂u 1 ∂v
= − cos θ + + sen θ = 0 ⇒ = ,
∂x ∂y ∂y ∂x ∂r r ∂θ
1 ∂u ∂v 1 ∂u ∂u ∂v ∂v
+ = −r sen θ + r cos θ + cos θ + sen θ
r ∂θ ∂r r ∂x ∂y ∂x ∂y
∂u ∂v ∂u ∂v 1 ∂u ∂v
= − + cos θ + + sen θ = 0 ⇒ =− .
∂x ∂y ∂y ∂x r ∂θ ∂r
♦
Problema 2.4. Sea f (z) una función de variable compleja definida como
1
f (z) = log ,
z−1
donde log z = ln |z| + iarg z, |z| > 0 y −π/4 < arg z < 7π/4. Determine el conjunto de
todos los números complejos z donde f (z) es analı́tica.
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 67
Solución. Por las condiciones del enunciado del problema, se tiene que la función log z
es analı́tica en todo el plano complejo excepto en el conjunto {z ∈ C : Re z ≥ 0, Im z =
−Re z}. Ahora, como
1 Re z − 1 1 −Im z
Re = 2 2
e Im = ,
z−1 (Re z − 1) + (Im z) z−1 (Re z − 1)2 + (Im z)2
1
entonces la función f (z) = log z−1 es analı́tica en todo el plano complejo excepto en el
conjunto
Re z − 1 −Im z Re z − 1
z∈C: ≥ 0, =−
(Re z − 1)2 + (Im z)2 (Re z − 1)2 + (Im z)2 (Re z − 1)2 + (Im z)2
C − {z ∈ C : Re z ≥ 1, Im z = Re z − 1}
Determine el conjunto de todos los pares ordenados (x, y) donde u(x, y) es armónica,
además, construya una función f (z) = u(x, y) + iv(x, y) y el mayor dominio D ⊂ C tales
que f (z) es analı́tica en D y f (−1 − i) = 1 − i.
Solución. Se tiene que la función z 2 + 1 es entera; por tanto, sus funciones componentes
Re (z 2 + 1) e Im (z 2 + 1) son funciones armónicas en todo R2 . Como u(x, y) = Re (z 2 + 1) =
x2 − y 2 + 1, para cada z = x + iy tal que x < 0, entonces u(x, y) es armónica en el conjunto
{(x, y) ∈ R2 : x ≤ 0}. Ahora, para (x, y) ∈ R2 tal que x > 0, se tiene que u(x, y) = x−3 y 3 ;
ası́
∂2u ∂2u
∇2 u(x, y) = (x, y) + (x, y)
∂x2 ∂y 2
= 12x−5 y 3 + 6x−3 y
= 6x−3 y x−2 y 2 + 1 = 0 ⇔ y = 0,
Es decir, la función
es analı́tica en D y f (−1 − i) = 1 − i. ♦
Problema 2.6. Sean α un número real y z0 un número complejo tal que Im z0 > 0.
Demostrar que el mapeo bilineal
iα z − z0
w=e ,
z − z0
|w|2 ≤ 1,
1
2
-1 − 21 1
1 x
2
z−1
Ahora, despejemos z = x + iy en función w = u + iv, a partir de la ecuación w = .
i−z
Se tiene que
iw + 1
iw − wz = z − 1 ⇔ z(w + 1) = iw + 1 ⇔ z= , (2.38)
w+1
de donde se deduce
u−v+1 u2 + v 2 + u − v
x= , y= . (2.39)
(u + 1)2 + v 2 (u + 1)2 + v 2
|(1 − v) + iu|2 ≤ |(u + 1) + iv|2 ,
(1 − v)2 + u2 ≤ (u + 1)2 + v 2 ,
|(1 − v) + iu|2 ≥ 1
4 |(u + 1) + iv|2 , ⇔ 4 (1 − v)2 + u2 ≥ (u + 1)2 + v 2 ,
(u + 1 )2 + (v − 1 )2 ≥ 1 ,
u2 + v 2 + u − v ≥ 0,
2 2 2
v 2 − 2v + 1 + u2 ≤ u2 + 2u + 1 + v 2 ,
v ≥ −u,
4v 2 − 8v + 4 + 4u2 ≥ u2 + 2u + 1 + v 2 , ⇔ u2 + v 2 − 32 u − 83 v + 1 ≥ 0,
(u + 1 )2 + (v − 1 )2 ≥ 1 . (u + 21 )2 + (v − 12 )2 ≥ 12 .
2 2 2
..
4 f (A)
.
3
..
.
2
-4 -3 -2 -1 1 2 3 4 u
-1
..
.
♦
2.2. Encuentre las funciones componentes, u(x, y) y v(x, y), de las siguientes funciones
f (z), donde z = x + iy.
1 c) f (z) = z 2 − 3z + 2.
a) f (z) = .
z2
z+1 z+2
b) f (z) = . d) f (z) = .
2z − 5 2z − 1 + i
z+2
2.3. Dada la función f (z) = , establecer en cada caso el conjunto de puntos
2z − 1
z ∈ C tales que:
b) lı́m (2z 3 − 3z 2 + z + i) z 2 − 2i z
z→2i e) lı́m
z→2i z2 + 4
z2 z2 + 1
c) lı́m f) lı́m
πi
z→e 4
z4 +z+1 z→i z6 + 1
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 71
z
g) lı́m i) lı́m Arg (1/z)
z→(3−4i) z+z z→i
2.6. Determine el conjunto de números complejos donde las siguientes funciones son
continuas:
2 Re z
3z − 5iz + 2
, z 6= 2i;
, Re z 6= Im z, Im z 6= 0;
a) f (z) = z 2 + iz + 6
Im z
7 c) h(z) = 1, Re z = Im z, Im z 6= 0;
5, z = 2i.
2
z +1
, z 6= −i; Re z, Im z = 0.
b) g(z) = z + i
−2i, z = −i.
2.8. Pruebe que cada una de las siguientes funciones satisfacen las ecuaciones de
Cauchy-Riemann:
2.11. Determine el conjunto D de todos los números complejos donde cada una de las
siguientes funciones f (z) es analı́tica y calcule su derivada:
a) f (z) = (z + 1)2 z4
e) f (z) =
1 1 + z4
b) f (z) = z +
z Log (z + 4)
3z − 1 f) f (z) =
c) f (z) = z2 + i
3−z g) f (z) = ee
z
2z + 1
d) f (z) = h) f (z) = sen (ez )
z(z 2 + 1)
2.12. Comprobar que cada una de las siguientes funciones son enteras:
y 2 − x2
f (z) = + i|1 − xy|.
2
Determine el conjunto de todos los números complejos z = x+ iy donde la función
f (z) es analı́tica.
2.14. Determine en qué conjunto D ⊂ R2 cada una de las siguientes funciones son
armónicas, y encuentre de ser posible su armónica conjugada.
a) e(2±3πi) = −e2
b) ez+πi = −ez
c) exp(z) = exp(iz), para z ∈ C
d) senh(z + πi) = − senh z, para z ∈ C
e) cosh(z + πi) = − cosh z, para z ∈ C
f) senh 2z = 2 senh z cosh z, para z ∈ C
g) Log (ei) = 1 + (π/2)i
h) Log (1 + i) = 12 (ln 2 + (π/2)i)
i) log 1 = 2nπi, para n = 0, ±1, ±2, . . .
CAPÍTULO 2. FUNCIONES DE VARIABLE COMPLEJA 73
a) ii b) (1 − i)4i c) (−i)3πi
1
2.19. Sea f (z) = Log . Determine el conjunto D de todos los números comple-
z+i
jos z donde f (z) es analı́tica.
2.20. Sea f (z) = Log (z 2 + 1). Determine el conjunto D de todos los números complejos
z donde f (z) es analı́tica.
2.21. Determine la rama principal y su derivada, para cada una de las siguientes fun-
ciones multivaluadas:
√ 2
a) f (z) = ez + 1 d) f (z) = z log z ≡ e(log z)
b) f (z) = cos(log z) e) f (z) = icos z ≡ ecos z log i
c) f (z) = log(ez + 1) f) f (z) = z sen z ≡ esen z log z
donde se utiliza un rama tal que se satisface f ′ (iπ) = 3π/2. Determine f ′ (z) y
f ′ (πi/2).
2.23. Demostrar que la transformación w = iz + i mapea el semiplano x > 0 sobre el
semiplano v > 1.
2.24. Demostrar que si c1 < 0, la imagen del semiplano x < c1 bajo la transformación
w = 1/z es el interior de un cı́rculo. ¿Cuál es la imagen cuando c1 = 0?
2.25. Demostrar que la imagen del semiplano y > c2 bajo la transformación w = 1/z es
el interior de un cı́rculo, siempre y cuando c2 > 0. Encontrar la imagen cuando
c2 < 0; también cuando c2 = 0.
74 2.10. PROBLEMAS PROPUESTOS
2.26. Sean α un número real y z0 un número complejo tal que |z0 | ≤ 1. Demostrar que
la transformación bilineal
iα z − z0
w=e ,
zz 0 − 1
mapea el disco |z| ≤ 1 sobre disco |w| ≤ 1.
y
1
-1 1 x
-1 D
z+1+i
Determine el conjunto transformado de D bajo el mapeo w = .
z+i
2.28. Sea S el conjunto de números complejos que se muestra en la siguiente figura.
y
1 S
x
−1 1
lı́m zn = z.
n→∞
zn = xn + i yn , para n = 0, 1, 2, . . .
1
zn = .
(1 + i)n
75
76 3.1. SERIE DE NÚMEROS COMPLEJOS
si, y sólo si
∞
X ∞
X
xn = X y yn = Y.
n=0 n=0
Ası́,
N
X N
X N
X
SN = zn = xn + i yn = XN + i YN .
n=0 n=0 n=0
Como
S = X + i Y = lı́m SN = lı́m XN + i lı́m YN
N →∞ N →∞ N →∞
si, y sólo si
lı́m XN = X y lı́m YN = Y,
N →∞ N →∞
P
entonces,
P∞podemosPconcluir que la serie ∞n=0 zn converge a S = X + i Y si, y sólo si las
series n=0 xn y ∞ y
n=0 n convergen respectivamente a X e Y.
El siguiente teorema da una idea completa del campo de convergencia de las series de
potencias. La demostración de este teorema se puede apreciar en [8].
De este modo, cuando 0 < α < ∞, existe un cı́rculo con centro en el punto z = z0 ,
en el interior del cual la serie (3.1) es absolutamente convergente y en el exterior del cual
la serie es divergente. Éste se llama cı́rculo de convergencia de la serie de potencias y su
radio R = 1/α, radio de convergencia de la misma. Los casos α = ∞ y α = 0 se pueden
78 3.1. SERIE DE NÚMEROS COMPLEJOS
1
R= .
α
Esta última se llama fórmula de Cauchy-Hadamard. Para las aplicaciones de la fórmula de
Cauchy-Hadamard, en muchos casos suele ser útil la relación siguiente:
r
n n! 1
lı́m n
= . (3.2)
n→∞ n e
La demostración de esta relación se deja como ejercicio para el lector.
Ejemplo 3.3. Hallar el cı́rculo y el radio de convergencia de la serie
∞
X nn
zn.
n=1
n!
para |z − z0 | < r0 .
Observación 3.2.
P∞ f (n) (z0 )
• El Teorema de Taylor garantiza que la serie n=0 n! (z − z0 )n converge a f (z),
para todo z tal que |z − z0 | < r0 .
• La serie de potencias (3.3) se denomina desarrollo en serie de Taylor o, simplemente,
desarrollo de Taylor de f (z) centrado en el punto z0 .
• Si z0 = 0, entonces el desarrollo de Taylor adquiere la forma
∞
X f (n) (0)
f (z) = zn
n=0
n!
El siguiente teorema nos garantiza que el desarrollo de Taylor de una función f (z)
alrededor de z0 , es la única serie de potencias que converge a f (z) en un disco centrado en
z0 .
Teorema 3.5. El desarrollo en serie de Taylor de una función f (z) alrededor de z0 es la
única serie en potencias de (z − z0 ) que converge a f (z) en todo punto de un disco centrado
en z0 .
Demostración. Por el absurdo. Supongamos que f (z) posee dos desarrollos distintos en
serie de Taylor alrededor de z0 válidos en un mismo dominio D ⊂ C, a saber:
∞
X ∞
X
n
f (z) = an (z − z0 ) y f (z) = bn (z − z0 )n
n=0 n=0
Como estos desarrollos son distintos, entonces se debe satisfacer que an 6= bn , para algún
n ≥ 0. Ahora, también tenemos que
∞
X
0 = f (z) − f (z) = (an − bn )(z − z0 )n ⇐⇒ a n = bn , para z ∈ D y todo n ≥ 0.
n=0
Observación 3.3. Sean f (z) y g(z) funciones tales que su desarrollo de Taylor centrado en
z0 está dado respectivamente por
∞
X
f (z) = an (z − z0 )n , para z tal que |z − z0 | < r1 ,
n=0
y
∞
X
g(z) = bn (z − z0 )n , para z tal que |z − z0 | < r2 .
n=0
Entonces del desarrollo de Taylor de f (z) + g(z) centrado en z0 es
∞
X
f (z) + g(z) = (an + bn ) (z − z0 )n , para z tal que |z − z0 | < r,
n=0
donde r = mı́n(r1 , r2 ).
El siguiente teorema nos permite calcular el radio de convergencia del mayor cı́rculo
para el cual existe la serie de Taylor de una función f (z).
Teorema 3.6. Consideremos el desarrollo en serie de Taylor (3.3) de una función f (z)
alrededor de z0 . Entonces, |z − z0 | < r es el mayor cı́rculo dentro del cual la serie converge
a f (z) para cada z en éste cı́rculo, donde r es la distancia entre z0 y el punto singular de
f (z) más cercano.
Obsérvese que este teorema no afirma que la serie de Taylor no converja fuera de
|z − z0 | = r. Sólo asevera que éste es el mayor cı́rculo en todo punto del cual la serie
converge a f (z). El cı́rculo en todo punto del cual la serie de Taylor (3.3) converge a f (z)
y el cı́rculo en todo punto del cual la serie converge no son necesariamente iguales. El
segundo de ellos puede tener un radio mayor. No obstante, se puede mostrar que cuando el
punto singular más cercano a z0 es tal que |f (z)| se hace infinito, los dos cı́rculos coinciden.
Éste es el caso en la mayor parte de las funciones que consideraremos.
CAPÍTULO 3. SERIES DE POTENCIAS Y SINGULARIDADES AISLADAS 81
X1h ∞ i
1
i) = (−1)n + 2−(n+1) (z − 1)n , |z − 1| < 1.
z(3 − z) n=0 3
Solución. a) Se tiene que f (z) = ez es una función entera y f (n) (z) = ez . Ası́, por los
Teoremas 3.4 y 3.6, el desarrollo de Maclaurin de ez es
∞
X 1 n
ez = z , |z| < ∞.
n!
n=0
eiz − e−iz
b) Se tiene que f (z) = sen z es una función entera, además, sen z = . De esta
2i
forma, utilizando el desarrollo de Macalurin de ez podemos escribir:
∞ ∞
X 1 n
X 1
(iz) − (−iz)n ∞
eiz − e−iz n=0
n! n=0
n! 1 X 1 n
sen z = = = (i − (−i)n )z n .
2i 2i 2i n!
n=0
♦
CAPÍTULO 3. SERIES DE POTENCIAS Y SINGULARIDADES AISLADAS 83
donde la serie converge a f (z) en cierto dominio o región y cn ∈ C. Ası́ pues, un desarrollo
de Laurent, a diferencia del desarrollo de Taylor, puede contener uno o más términos con
(z − z0 ) elevado a una potencia negativa. También puede contener potencias positivas de
(z − z0 ).
Normalmente los desarrollos de Laurent se obtienen a partir de los desarrollos de Taylor.
Por ejemplo, para encontrar el desarrollo de Laurent de f (z) = e1/z centrado en z0 = 0, se
considera el desarrollo de Maclaurin de es ,
∞
s
X 1 n
e = s ,
n!
n=0
• Si f (z) es analı́tica en todos los puntos de la región |z−z0 | < r2 , entonces el desarrollo
de Laurent de f (z) centrado en z0 , se convierte en el desarrollo de Taylor de f (z)
centrado en z0 .
84 3.1. SERIE DE NÚMEROS COMPLEJOS
• Si f (z) es analı́tica en todos los puntos de la región |z−z0 | < r2 excepto en en el punto
z0 , entonces el desarrollo de Laurent es válido en toda la región 0 < |z − z0 | < r2 .
• Los coeficientes de la serie de Laurent se pueden definir como una integral que invo-
lucra a la función f (z). Este último hecho permite calcular integrales utilizando los
resultados de series de potencias, lo cual se estudiará en el Capı́tulo 4.
1
f (z) = .
1+z
Solución. Observamos que el único punto singular de f (z) es z = −1, entonces existen dos
regiones donde f (z) posee desarrollos de Laurent centrados en z0 = 0, a saber: a) |z| < 1,
y b) |z| > 1.
a) Como f (z) es analı́tica en todo punto del disco |z| < 1, el desarrollo de Laurent de f (z)
centrado en z0 = 0 válido en el dominio |z| < 1, coincide con su desarrollo de Maclaurin.
En otras palabras, el desarrollo de Laurent de f (z) centrado en z0 = 0 válido en el dominio
|z| < 1, es:
∞
X
f (z) = (−1)n z n , |z| < 1.
n=0
X∞ X∞
• f (z) = − (1 − i)n (z − i)−(n+1) − (2 − i)−(n+1) (z − i)n , para z ∈ C tal que
√ n=0 √ n=0
2 < |z − i| < 5,
∞
X √
• f (z) = [(2 − i)n − (1 − i)n ] (z − i)−(n+1) , para z ∈ C tal que |z − i| > 5.
n=0
CAPÍTULO 3. SERIES DE POTENCIAS Y SINGULARIDADES AISLADAS 85
ahora, como
d 1 1
= − 2,
dz z z
entonces utilizando el desarrollo de Taylor de 1/z centrado en z0 = 1, podemos escribir:
∞ ∞
1 d 1 X d n n
X
= − = − [(−1) (z − 1) ] = (−1)n+1 n (z − 1)n−1 , |z| < 1.
z2 dz z dz
| n=0 {z } n=0
|z|<1
Como
d 1 1
= ,
dz (1 − z) (1 − z)2
entonces utilizando el desarrollo de Maclaurin de 1/(1 − z), podemos escribir:
∞ ∞
z d 1 X d n X n
= z = z z = nz , |z| < 1.
(1 − z)2 dz (1 − z) dz
n=0 n=0
| {z }
|z|<1
♦
86 3.2. SINGULARIDADES AISLADAS
1
Ejemplo 3.11. Determine la parte principal de f (z) = en cada uno de sus puntos
z(z − 1)
singulares.
Los puntos singulares aislados se pueden caracterizar según la forma que adquiere la
parte principal de la función en tales puntos. Se distinguen tres tipos de singularidades
aisladas: polo de orden m, punto singular esencial y punto singular removible. A continua-
ción se describen cada una de ellas.
Como lı́mz→z0 |(z − z0 )m | = 0, entonces por (3.5) se tiene que lı́mz→z0 |f (z)| = ∞.
88 3.2. SINGULARIDADES AISLADAS
El teorema anterior no solo nos permite identificar si un punto es un polo, sino también
el orden del mismo. Basándose en este teorema y, particularmente, en la ecuación (3.5),
las siguientes reglas nos permiten identificar si un punto singular aislado z0 es un polo y,
además, calcular el orden del mismo.
1
Ejemplo 3.12. Sea f (z) = . Probar que 0 y 1 son polos simples de f (z).
z(z − 1)
Solución. Utilicemos el Teorema 3.8 y las Reglas I y II, para determinar que los puntos
0 y 1 son polos simples de f (z). Se tiene que
Es común encontrarse con problemas en los que se quiere determinar el orden de los
polos de una función de la forma f (z) = p(z)/q(z), por ello les dedicaremos una atención
especial. El siguiente teorema nos permite identificar si un punto es un polo, y además su
orden, para una función f (z) = p(z)/q(z).
Teorema 3.9. Sea z0 ∈ C. Sea f (z) una función tal que se puede escribir como
p(z)
f (z) = ,
q(z)
por tanto,
q(z) q (m) (z0 )
lı́m = 6= 0, ∞,
z→z0 (z − z0 )m m!
lo cual implica que
(z − z0 )m
lı́m (z − z0 )m f (z) = lı́m p(z) lı́m 6= 0, ∞.
z→z0 z→z0 z→z0 q(z)
Según la Regla I, lo anterior indica que z0 es un polo de orden m de f (z).
Ejemplo 3.13. Verificar que todos los puntos singulares aislados de la función
ez
f (z) =
sen z
son polos simples.
90 3.2. SINGULARIDADES AISLADAS
Solución. Se tiene que f (z) = p(z)/q(z), donde p(z) = ez y q(z) = sen z. Ahora, los
puntos singulares aislados de f (z) son: zn = nπ, para n = 0, ±1, ±2, . . .. Como p(zn ) 6= 0,
q(zn ) = 0 y q ′ (zn ) = cos(nπ) 6= 0, para todo n, entonces por el Teorema 3.9 todos los
puntos singulares de f (z) son polos simples. ♦
Ejercicio 3.3. Verifique la existencia de los polos y el orden de los mismos para cada una
de las funciones dadas.
1
a) El punto z0 = 0 es un polo de orden m = 2 de la función f (z) = z
.
z(e − 1)
(z + 1)
b) Los puntos z0 = 3i y z1 = −3i son polos simples de la función f (z) = .
(z 2 + 9)
z 2 − 2z + 3
c) El punto z0 = 2 es un polo simple de la función f (z) = .
(z − 2)
senh z
d) El punto z0 = 0 es un polo de orden m = 3 de la función f (z) = .
z4
Por lo tanto, para determinar si un punto singular aislado z0 es un punto singular removible
de f (z), basta con verificar que el lı́mz→z0 f (z) existe y es finito.
>> syms z n
>> symsum ( z ^n ,n ,0 , Inf )
ans =
p i e c e w i s e ([1 <= z , Inf ] , [ abs ( z ) < 1 , -1/( z - 1) ])
En esta caso, la respuesta que da MATLAB, piecewise([1 <= z, Inf], [abs(z) < 1,
-1/(z - 1)]), indica que la serie no converge para |z| ≥ 1, pero si converge a 1/(1 − z)
cuando |z| < 1. Veamos otros ejemplos. Seguidamente se calcula:
∞
X zn
• = −Log (1 − z), |z| < 1
n=1
n
92 3.3. SERIES DE POTENCIAS Y SINGULARIDADES CON MATLAB
∞
X n n z
• z = 2 , |z| < 2
n=1
2n 2 z
−1
2
∞
X n(z + 2)n 2z + 4
• = , |z + 2| < 2
2n z2
n=0
>> syms z n
>> symsum ( z ^ n /n ,n ,1 , Inf )
ans =
p i e c e w i s e ([1 <= z , Inf ] , [ abs ( z ) <= 1 and z ~= 1 , - log (1 - z ) ])
>> symsum (( n * z ^ n ) /2^n ,n ,1 , Inf )
ans =
p i e c e w i s e ([ abs ( z ) < 2 , z /(2*( z /2 - 1) ^2) ])
>> symsum (( n *( z +2) ^ n ) /2^ n ,n ,1 , Inf )
ans =
p i e c e w i s e ([ abs ( z /2 + 1) < 1 , (2*( z + 2) ) / z ^2])
Hemos visto que el comando symsum permite calcular a qué y en dónde converge una
serie de potencias. Ahora utilicemos el comando taylor para determinar el desarrollo
de Taylor de una función f (z). Con taylor(f,z,a) se calculan, por defecto, los cinco
primeros términos (f (a) + f ′ (a)(z − a) + f ′′ (a)(z − a)2 + f ′′′ (a)(z − a)3 + f (iv) (a)(z − a)4 )
del desarrollo de Taylor de f (z) al rededor de z0 = a. El desarrollo de Maclaurin de f (z)
se calcula usando taylor(f). Observe el siguiente ejemplo donde se calculan los cinco
primeros términos de los desarrollos de Maclaurin de ez , cos z y 1/(1 − z).
>> syms z
>> taylor ( exp ( z ) )
ans =
z ^5/120 + z ^4/24 + z ^3/6 + z ^2/2 + z + 1
>> taylor ( cos ( z ) )
ans =
z ^4/24 - z ^2/2 + 1
>> taylor (1/(1 -z ) )
ans =
z ^5 + z ^4 + z ^3 + z ^2 + z + 1
>> syms z
>> taylor (( z +1) /( z -5) ^6 , z , -1 , ’ Order ’ ,8)
ans =
z /46656 + ( z + 1) ^ 2 / 4 6 6 5 6 + (7*( z + 1) ^3) / 5 5 9 8 7 2 +
(7*( z + 1) ^4) / 1 2 5 9 7 1 2 + (7*( z + 1) ^5) / 3 3 5 9 2 3 2 +
(7*( z + 1) ^6) / 1 0 0 7 7 6 9 6 + (77*( z + 1) ^7) / 3 6 2 7 9 7 0 5 6 + 1 / 4 6 6 5 6
b1 z m + b2 z m−1 + b3 z m−2 + · · · + bm
f (z) = ,
a1 z n + a2 z n−1 + a3 z n−2 + · · · + an
con n > m, hemos diseñado la función LaurentRacional (ver Listado 3.1), que permite
hallar el desarrollo de Laurent de f (z) alrededor de uno de sus puntos singulares z0 ∈ C.
f u n c t i o n [ dsl ]= L a u r e n t R a c i o n a l (b ,a , z0 , m )
% L a u r e n t R a c i o n a l Halla el d e s a r r o l l o de L a u r e n t de la f u n c i ó n
% r a c i o n a l propia
% b (1) z ^ m + b (2) z ^(m -1) +...+ b ( m )
% f(z) = ------------------------------------
% a (1) z ^ n + a (2) z ^(n -1) +...+ a ( n )
% al r e d e d o r de uno de sus puntos s i n g u l a r e s z0
% m es el número de t é r m i n o s de los d e s a r r o l l o s
% de p o t e n c i a s o b t e n i d o s
% dls guarda el d e s a r r o l l o de L a u r e n t
syms z ;
dsl = sym (0) ; dsl1 = sym (0) ;
[r , p , k ] = r e s i d u e(b , a ) ;
np = length ( r ) ;
fi = sym ( zeros ( size ( r ) ) ) ;
fi (1) = sym ( r (1) ) /( z - sym ( p (1) ) ) ;
mr = 1;
if np >1
for t =2: np
if p ( t ) == p (t -1)
fi ( t ) = sym ( r ( t ) ) /( z - sym ( p ( t ) ) ) ^( sym ( mr +1) ) ;
mr = mr +1;
else
fi ( t ) = sym ( r ( t ) ) /( z - sym ( p ( t ) ) ) ;
mr = 1;
end
end
end
for t =1: np
if sym ( p ( t ) ) == z0
dsl = dsl + fi ( t ) ;
else
dsl1 = dsl1 + fi ( t ) ;
end
end
dsl = dsl + taylor ( dsl1 ,z , z0 , ’ Order ’ ,m ) ;
if ~ i s e m p t y( k )
dsl = dsl + sym ( k ) ;
end
end
94 3.3. SERIES DE POTENCIAS Y SINGULARIDADES CON MATLAB
z+1 z+1
f (z) = = 3
(z − 1)2 (z + 3) z + z 2 − 5z + 3
Por otra parte, el comando limit puede usarse para determinar la naturaleza de una
singularidad aislada, es decir, podemos identificar si es un polo, un punto singular remo-
vible o un punto singular esencial. A continuación verificamos cada una de las siguientes
afirmaciones:
a) Se tiene
>> syms z
>> f ( z ) = 1/( z *( exp ( z ) -1) )
f(z) =
1/( z *( exp ( z ) - 1) )
>> limit ( abs ( f ) ,z ,0)
ans ( z ) =
Inf
b) Se tiene
>> syms z
>> f ( z ) = exp (1/( z - i ) )
f(z) =
exp (1/( z - i ) )
>> limit (f ,z , i )
ans ( z ) =
NaN
En este caso, el sı́mbolo NaN indica que el lı́mite no existe. Por lo tanto, z0 = i es un punto
singular esencial de f (z) = e1/(z−i) .
c) Se tiene
>> syms z
>> f ( z ) = cos ( z ) /(2* z - pi )
f(z) =
- cos ( z ) /( pi - 2* z )
>> limit (f ,z , pi /2)
ans ( z ) =
-1/2
Solución. Es claro que la serie de potencias dada se puede escribir como la resta de dos
series,
∞ h
X i ∞
X ∞
X
(−1)n − 2−(n+1) (z − 1)n = (−1)n (z − 1)n − 2−(n+1) (z − 1)n .
n=0 n=0 n=0
an = (−1)n y z0 = 1.
P
Los coeficientes y el centro de la serie ∞ n=0 2
−(n+1) (z − 1)n son: b = 2−(n+1) y z = 1.
n 0
Ası́, aplicando el criterio del cociente obtenemos:
|bn+1 | 2−(n+2) 1
α = lı́m = lı́m −(n+1) = .
n→∞ |bn | n→∞ 2 2
P
Luego, |z − 1| < 2 es el disco de convergencia de la serie ∞ n=0 2
−(n+1) (z − 1)n .
P∞ n
−(n+1) (z − 1)n es
Por lo tanto, la región de convergencia de la serie n=0 (−1) − 2
|z − 1| < 1. ♦
2z + 1
f (z) = ,
z(2z − 1)
centrado en z0 = − 12 .
2 1
f (z) = 1 − .
z− 2
z
! ∞
1 −1 X 1 n 1 1
= 2 = (−2) 2n z + , z + < .
z 1 − 2(z + 21 ) n=0
2 2 2
∞ ∞
X 1 n X
n 1 n
f (z) = (−2) z+ − (−2) 2 z+
2 2
n=0 n=0
| {z } | {z }
|z+ 21 |<1 |z+ 21 |< 21
∞
X
n+1 1 n 1 1
= (2 − 2) z + , z + < .
n=0
2 2 2
♦
CAPÍTULO 3. SERIES DE POTENCIAS Y SINGULARIDADES AISLADAS 97
i
f (z) = ,
z2
centrado en z0 = −i.
Como f (z) = −g′ (z), entonces el desarrollo de Taylor de f (z) centrado en z0 = −i, está
dado por:
" ∞ # ∞
d X X
f (z) = − − i−n (z + i)n = i−n n (z + i)n−1 , |z + i| < 1.
dz n=0 n=0
♦
98 3.4. PROBLEMAS RESUELTOS
1
f (z) = ,
z(z + i)
centrados en z0 = 1.
i i
f (z) = − .
z+i z
i i
f (z) = −
z+i z
i 1 1
= −i
1 + i 1 + (1 + i)−1 (z − 1) 1 + (z − 1)
∞ ∞
i X X
= (−1)n (1 + i)−n (z − 1)n − i (−1)n (z − 1)n
1+i
|n=0 {z } |n=0 {z }
√
|z−1|< 2 |z−1|<1
∞
X h i
= i(−1)n (1 + i)−(n+1) − 1 (z − 1)n , |z − 1| < 1.
n=0
√
b) El desarrollo de Laurent de f (z) centrado en z0 = 1, válido en 1 < |z − 1| < 2, está
dado por:
i i
f (z) = −
z+i z
i 1 i 1
= −
1 + i 1 + (1 + i)−1 (z − 1) z − 1 1 + (z − 1)−1
∞
X ∞
X
= i(−1)n (1 + i)−(n+1) (z − 1)n − i(−1)n (z − 1)−(n+1)
|n=0 {z } |n=0 {z }
√
|z−1|< 2 |z−1|>1
∞ ∞
X X √
= i(−1)n (1 + i)−(n+1) (z − 1)n + i(−1)n+1 (z − 1)−(n+1) , 1 < |z − 1| < 2.
n=0 n=0
CAPÍTULO 3. SERIES DE POTENCIAS Y SINGULARIDADES AISLADAS 99
√
c) El desarrollo de Laurent de f (z) centrado en z0 = 1, válido en |z − 1| > 2, está dado
por:
i i
f (z) = −
z+i z
i 1 i 1
= −
z − 1 1 + (1 + i)(z − 1)−1 z − 1 1 + (z − 1)−1
∞
X ∞
X
n n −(n+1)
= i(−1) (1 + i) (z − 1) − i(−1)n (z − 1)−(n+1)
n=0
| {z } |n=0 {z }
√
|z−1|> 2 |z−1|>1
∞
X √
= i(−1)n [(1 + i)n − 1] (z − 1)−(n+1) , |z − 1| > 2.
n=0
1
f (z) = ,
z (z − 1) (z − 2)
centrados en z0 = 0.
1/2 1 1/2
f (z) = − + .
z−2 z−1 z
La función f (z) dada tiene tres desarrollos de Laurent centrados en z0 = 0, cuyos dominios
de validez son los siguientes anillos: a) 0 < |z| < 1, b) 1 < |z| < 2, y c) |z| > 2.
a) El desarrollo de Laurent de f (z) centrado en z0 = 0, válido en 0 < |z| < 1, está dado
por:
1/2 1 1/2
f (z) = − +
z − 2 z−1 z
1 1 1 1
= − −1
+ + z −1
4 1−2 z 1−z 2
∞ ∞
1 X −n n X n 1 −1
= − 2 z + z + z
4 n=0 |2 {z }
| {z } |n=0 {z } |z|>0
|z|<2 |z|<1
∞ h i
X 1 −1
= 1 − 2−(n+2) z n + z , 0 < |z| < 1.
n=0
2
100 3.4. PROBLEMAS RESUELTOS
b) El desarrollo de Laurent de f (z) centrado en z0 = 0, válido en |z| > 2, está dado por:
1/2 1 1/2
f (z) = − +
z − 2 z−1 z
1 1 −1 1 1
= − −1
−z −1
+ z −1
4 1−2 z 1−z 2
∞ ∞
1 X −n n X −(n+1) 1 −1
= − 2 z − z + z
4 |2 {z }
n=0
| {z } |n=0 {z } |z|>0
|z|<2 |z|>1
∞ ∞
X X 1 −1
= (−1)2−(n+2) z n + (−1)z −(n+1) + z , 1 < |z| < 2.
2
n=0 n=0
c) El desarrollo de Laurent de f (z) centrado en z0 = 0, válido en 1 < |z| < 2, está dado
por:
1/2 1 1/2
f (z) = − +
z−2 z−1 z
−1
z 1 −1 1 1
= −1
−z −1
+ z −1
2 1 − 2z 1−z 2
∞ ∞
X X 1
= 2n−1 z −(n+1) − z −(n+1) + z −1
n=0 |2 {z }
| {z } |n=0 {z } |z|>0
|z|>2 |z|>1
∞
X 1
= 2n−1 − 1 z −(n+1) + z −1 , |z| > 2.
2
n=0
1
f (z) = .
(z + 2)3 e1/(z+2)3
Por lo que la parte principal de f (z) en z0 = −2 posee infinitos términos distintos de cero;
por lo tanto, z0 = −2 es un punto singular esencial de f (z). ♦
CAPÍTULO 3. SERIES DE POTENCIAS Y SINGULARIDADES AISLADAS 101
p(z)
f (z) = ,
q(z)
donde ( sen z
, z 6= 0,
p(z) = z y q(z) = z 4 .
1, z = 0,
Además, se tiene que: p(z) y q(z) son analı́ticas en z0 = 0, p(0) = 1 6= 0, q(0) = q ′ (0) =
q ′′ (0) = q ′′′ (0) = 0 y q (iv) (0) = 24 6= 0. Entonces, z0 = 0 es un polo de orden 4 de f (z). ♦
p(z)
f (z) = ,
q(z)
2
donde p(z) = e−z y q(z) = z 2 . Además, se tiene que: p(z) y q(z) son analı́ticas en z0 = 0,
p(0) = 1 6= 0, q(0) = q ′ (0) = 0 y q ′′ (0) = 2 6= 0. Por lo tanto, z0 = 0 es un polo de orden 2
de f (z). ♦
∞ ∞
X n(z + 2)n X (−1)n−1
k) l) zn
n=0
2n n=1
n
3.3. Determine y clasifique todas las singularidades aisladas de las siguientes funciones
complejas.
1 cos z
a) f (z) = . e) f (z) = .
2 2
z (z − 4z + 5) z2
z+i f) f (z) = z 3 e1/z .
b) f (z) = 3 .
z + z 2 − 8z − 12
z−1 z
c) f (z) = 4 . g) f (z) = .
z − z 2 (1 + i) + i sen z
cos(z − 1) ez
d) f (z) = 4 . h) f (z) =
z − z 2 (1 + i) + i z(1 − e−z )
4
Integración Compleja
Es de notar que en los cursos de cálculo elemental, particularmente, en el tema de
integración de funciones definidas de R en R, elR concepto de integral tenı́a la siguiente in-
b
terpretación geométrica: el valor de la integral a h(x) dx corresponde al área comprendida
entre el eje real y la gráfica de la función h(x) en el intervalo [a, b]. Desafortunadamente,
esta interpretación geométrica no tiene sentido para funciones de variable compleja. En
sı́, la integración de funciones de variable compleja es una herramienta esencial para el
desarrollo teórico de las ideas del cálculo operacional, en particular, el estudio de las trans-
formadas de Laplace, Fourier y Z, que son temas indispensables para la compresión de
ciertos conceptos estudiados en Ingenierı́a Eléctrica. En este capı́tulo se describen los con-
ceptos básicos de la integración compleja, comenzando con la integral definida, pasando
luego por la integral de lı́nea, el Teorema de Cauchy-Goursat y la primitiva de una función,
para finalizar con el Teorema de los Residuos.
Z b Z b Z b
F (t) dt = U (t) dt + i V (t) dt
a a a
103
104 4.2. INTEGRACIÓN DE LÍNEA
y
C
β = z(b)
α = z(a) b b
b
z(t)
b b b
a t b x
−1 1 x −1 1 x
−1
(
z1 (t) = eit , 0 ≤ t ≤ 2π −t + i(1 + t), −1 ≤ t ≤ 0,
z2 (t) =
−t + i(1 − t), 0 < t ≤ 1.
Definición 4.3 (Contorno). Un contorno o curva suave a tramos, es una curva que consta
de un número finito de curvas suaves unidas por sus extremos.
Ejemplo 4.3. Seguidamente se aprecia una representación gráfica de un contorno confor-
mado por seis curvas suaves C1 , C2 , . . . , C6 .
C2
C3
C1
C4
C6 C5
Definición 4.4 (Contorno cerrado simple). Se dice que C es un contorno cerrado simple,
si C es un contorno y z(a) = z(b) y z(t1 ) 6= z(t2 ), para todo t1 6= t2 ∈ (a, b).
Ejemplo 4.4. En la siguiente gráfica se muestran dos curvas C y S. La curva C es un
contorno cerrado simple, en cambio S no lo es.
C z(a) = z(b)
S
z(t1 ) = z(t2 )
z(a) = z(b)
106 4.2. INTEGRACIÓN DE LÍNEA
se le asocia el contorno −C, cuya parametrización está definida por la ecuación z = z(−t),
donde −b ≤ t ≤ −a. Gráficamente, el contorno −C es el mismo contorno C pero recorrido
en sentido contrario.
con extremo inicial α = z(a) y extremo final β = z(b). Supongamos que f (z) = u(x, y) +
i v(x, y) es continua a trozos en C, es decir, las partes real e imaginaria, u(x(t), y(t)) y
v(x(t), y(t)), de f (z(t)) son funciones de t continuas a trozos. Bajo estas condiciones, se
define la integral de lı́nea de f (z) a lo largo de C como:
Z Z b
f (z) dz = f (z(t)) z ′ (t) dt,
C a
Z Z b
v) f (z) dz ≤
|f (z(t)) z ′ (t)| dt.
C a
CAPÍTULO 4. INTEGRACIÓN COMPLEJA 107
Ası́,
Z Z 2π Z 2π
1 1 it
dz = ie dt = i dt = 2πi.
|z|=1 z 0 eit 0
deben anularse para cualquier contorno cerrado simple C. En todos estos casos, el in-
tegrando es una función entera. Ahora, desde el punto de vista teórico, el Teorema de
Cauchy-Goursat permite introducir el concepto de primitiva de una función (lo cual se
estudiará más adelante), entre otros conceptos. R
Obsérvese que la dirección de integración en la integral C f (z) dz = 0 no afecta el
resultado, pues Z Z
f (z) dz = − f (z) dz.
−C C
En el siguiente ejemplo se verifica la validez del Teorema de Cauchy-Goursat.
Ejemplo 4.7. Verifique que Z
z n dz = 0,
C
donde n es un entero positivo y C es la circunferencia |z| = r, con r > 0.
Solución. Como n > 0, la función f (z) = z n es entera, luego por el Teorema de Cauchy-
Goursat Z
z n dz = 0.
C
Veamos que esto es efectivamente cierto. Una parametrización de |z| = r es:
Luego,
Z Z 2π Z 2π
z n dz = (r eit )n (ir eit ) dt = ir n+1 e(n+1)it dt
C 0 0
Z 2π
= ir n+1 [cos((n + 1)t) + i sen ((n + 1)t)] dt
0
n+1 sen ((n + 1)t) − i cos((n + 1)t) 2π
= ir = 0.
(n + 1) 0
D1 D2
C C1
C2
Solución. Sea R la región cerrada que consta de todos los puntos dentro y sobre |z| = 2
excepto los puntos interiores a |z + 1| = 1/2, |z| = 1/2 y |z − 1| = 1/2. El integrando es
analı́tico excepto en los puntos z = 0 y z = ±1, y estos tres puntos no pertenecen a R.
Por lo tanto, aplicando el Teorema 4.3 concluimos que
Z
dz
2 2
= 0.
B z (z − 1)
♦
110 4.4. INTEGRAL DEFINIDA
pero Z Z Z Z Z
f (z) dz = f (z) dz + f (z) dz = f (z) dz − f (z) dz,
C C1 −C2 C1 C2
por lo tanto, Z Z
f (z) dz = f (z) dz,
C1 C2
lo cual indica que la integral desde z1 hasta z2 es ası́ independiente del contorno seguido,
en tanto ese contorno se encuentre dentro de D.
Definición 4.7 (Integral indefinida o primitiva). Sea f (z) una función analı́tica en un
dominio simplemente conexo D ⊂ C. Sea z0 ∈ D. La función F (z) definida en D por
Z z
F (z) = f (s) ds + c, (4.1)
z0
para indicar todas las posibles primitivas de f (z). Note que la primitiva de cada una de las
funciones elementales complejas, son las mismas primitivas de las funciones elementales de
variable real; por ejemplo, la primitiva de ez es, por supuesto, ella misma, y la primitiva
de sen z es − cos z, etc.
Por otra parte, el valor de la constante c correspondiente a una primitiva especı́fica
Z z
f (s) ds
z0
Para determinar el valor de c observemos que el lado izquierdo de esta ecuación es cero
cuando z = π. Por lo tanto,
0 = −π cos π + sen π + c,
De la ecuación (4.1) se infiere que una integral definida se puede evaluar de igual forma
que en el cálculo elemental:
Z β
β
f (z) dz = F (β) − F (α) = F (z) ,
α α
es decir, usando la Regla de Barrow, que es una aplicación del Teorema Fundamental del
Cálculo.
112 4.5. FÓRMULA INTEGRAL DE CAUCHY
1 1
= .
z2 + 4 (z − 2i)(z + 2i)
Observamos que el factor (z − 2i) se anula dentro del contorno de integración y que (z + 2i)
no se anula ni sobre el contorno ni en su interior. Escribiendo la integral considerada en
la forma
1
Z
(z + 2i)
dz
|z−i|=2 (z − 2i)
1
Como f (2i) = , entonces el valor de la integral considerada es
4i
Z
1 π
2
dz = 2πif (2i) = .
|z−i|=2 z + 4 2
♦
CAPÍTULO 4. INTEGRACIÓN COMPLEJA 113
Pasemos ahora a ver que si una función es analı́tica en un punto, entonces sus derivadas
de todos los órdenes existen en ese punto y son también analı́ticas. Tal resultado se conoce
como extensión de la fórmula integral de Cauchy.
Tomando f (z) = 1/(z + 2i)2 y z0 = 2i, por el Teorema 4.6 podemos escribir:
Z Z
′ 1 f (z) 1 1
f (2i) = dz = dz.
2πi |z−i|=2 (z − 2i)2 2πi |z−i|=2 (z 2 + 4)2
−2
Como f ′ (z) = , entonces el valor de la integral considerada es
(z + 2i)3
Z
1 π
2 2
dz = 2πif ′ (2i) = .
|z−i|=2 (z + 4) 16
4.6 Residuo
Definición 4.8. Sean C un contorno cerrado simple y z0 ∈ C un punto interior a C. Sea
f (z) una función analı́tica sobre C y en todo punto de su interior, salvo en z0 . El residuo
de f (z) en z0 , denotado por Res [f (z)], se define como
z=z0
Z
1
Res [f (z)] = f (z) dz.
z=z0 2πi C
El siguiente teorema describe la relación que existe entre el residuo Res [f (z)] y una y
z=z0
sólo una serie de Laurent de f (z) centrada en z0 .
114 4.6. RESIDUO
Teorema 4.7. Sea z0 un punto singular aislado de una función f (z). Entonces, el residuo
de f (z) en z0 es igual al coeficiente de (z − z0 )−1 en la serie de Laurent que representa a
f (z) en el anillo
0 < |z − z0 | < r,
Demostración. Como z0 es un punto singular aislado de f (z), entonces existe r > 0 tal
que el desarrollo de Laurent de f (z) centrado en z0 , válido en el anillo 0 < |z − z0 | < r,
está dado por:
∞
X ∞
X
f (z) = an (z − z0 )n + bn (z − z0 )−n .
n=0 n=1
Sea C un contorno cerrado simple contenido en el anillo 0 < |z − z0 | < r tal que z0 es un
punto interior de C. Ası́, integrando en ambas partes de la ecuación anterior obtenemos:
Z Z "X
∞ ∞
X
#
n −n
f (z) dz = an (z − z0 ) + bn (z − z0 ) dz
C C n=0 n=1
∞
X Z X∞ Z
= an (z − z0 )n dz + bn (z − z0 )−n dz,
n=0 C n=1 C
además,
Z
(z − z0 )n dz = 0, para n ≥ 0,
C
Z
(z − z0 )−1 dz = 2πi,
ZC
(z − z0 )−n dz = 0, para n ≥ 2,
C
por lo tanto, Z
f (z) dz = b1 2πi,
C
c) Res z 4 sen (1/z)
z=0
CAPÍTULO 4. INTEGRACIÓN COMPLEJA 115
Teorema 4.8. Sea z0 un punto singular aislado de una función f (z). Si para cierto entero
positivo m, la función
φ(z) = (z − z0 )m f (z)
se puede definir en z0 de modo que sea analı́tica ahı́ y φ(z0 ) 6= 0, entonces f (z) tiene un
polo de orden m en z0 y, además,
φ(z0 ) = lı́m (z − z0 )m f (z), m = 1;
z→z0
Res [f (z)] = (4.3)
z=z0
φ(m−1) (z0 )
, m > 1.
(m − 1)!
Demostración. Como φ(z) = (z − z0 )m f (z) es analı́tica en z0 y φ(z0 ) 6= 0, entonces la
función f (z) se puede escribir como
p(z)
f (z) = ,
q(z)
donde p(z) = φ(z) y q(z) = (z − z0 )m , además,
q(z0 ) = q ′ (z0 ) = · · · = q (m−1) (z0 ) = 0, y q (m) (z0 ) = m! 6= 0.
Luego, por el Teorema 3.9 z0 es un polo de orden m de f (z).
Pasemos ahora a demostrar la fórmula (4.3). Como φ(z) es analı́tica en z0 ella posee
desarrollo de Taylor centrado en z0 , válido en cierto disco |z − z0 | < r, dado por
φ(m) (z0 )
φ(z) = (z − z0 )m f (z) = φ(z0 ) + φ′ (z0 )(z − z0 ) + · · · + (z − z0 )m + · · · .
m!
Ası́, en cada punto z del disco |z − z0 | < r excepto en z0 , se tiene que
φ(z0 ) φ′ (z0 ) φ(m−1) (z0 ) φ(m) (z0 )
f (z) = + + · · · + + + ··· ,
(z − z0 )m (z − z0 )m−1 (m − 1)!(z − z0 ) m!
que es el desarrollo de Laurent de f (z) centrado en z0 válido en el anillo 0 < |z − z0 | < r,
con
φ(m−1) (z0 )
b1 = .
(m − 1)!
Entonces, por el Teorema 4.7 la ecuación (4.3) es cierta.
116 4.6. RESIDUO
Log z
c) Res 4 .
z=1 z (z − 1)2
sen z − z
d) Res , para n = 0, ±1, ±2, . . .
z=nπi z senh z
z−2
Solución. Sea f (z) = . Los puntos singulares de f (z) son z0 = 0 y z1 = 1, que
(z − 1)z
son puntos interiores a C, además, z0 y z1 son polos simples de f (z). Luego,
b0 + b1 z + · · · + bM z M
f (z) = ,
a0 + a1 z + · · · + z N
Teorema 4.10. Sea f (z) una función racional propia dada por
b0 + b1 z + · · · + bM z M
f (z) = ,
a0 + a1 z + · · · + z N
(i) Si todos los polos de f (z) son simples, la expansión en fracciones parciales de f (z)
es:
A1 A2 AN
f (z) = + + ··· + ,
(z − p1 ) (z − p2 ) (z − pN )
donde los números complejos Ak , denominados coeficientes, se calculan como
Ak = Res [f (z)] ,
z=pk
para k = 1, 2, . . . , N .
(ii) Si todos los polos de f (z) son simples, excepto el polo pℓ que es de orden rℓ , la
expansión en fracciones parciales de f (z) es:
A1 A2 Aℓ−1
f (z) = + + ··· +
(z − p1 ) (z − p2 ) (z − pℓ−1 )
Aℓ,1 Aℓ,2 Aℓ,rℓ
+ + 2
+ ··· +
(z − pℓ ) (z − pℓ ) (z − pℓ )rℓ
Aℓ+1 Aℓ+2 AT
+ + + ··· + ,
(z − pℓ+1 ) (z − pℓ+2 ) (z − pT )
donde
Ak = Res [f (z)] ,
z=pk
para k = 1, 2, . . . , T , y k 6= ℓ; y
" #
1 d(rℓ −) rℓ
Aℓ, = · [(z − p ℓ ) f (z)] ,
(rℓ − )! dz (rℓ −)
z=pℓ
para = 1, 2, . . . , rℓ .
El teorema anterior se puede extender a funciones racionales propias que tienen dos o
más polos cuyo orden es mayor que 1. En el siguiente ejemplo se muestra tal extensión.
Ejemplo 4.15. Halle la expansión en fracciones parciales de
144 z 2 + 144 z + 144
f (z) = .
(z − 3)2 (z − 2)2 (z + 1)
Solución. Los puntos p1 = −1, p2 = 2 y p3 = 3 son los polos de f (z). Se observa que p1
es de orden 1 y p2 y p3 son de orden 2. De esta forma, la expansión en fracciones parciales
de f (z) es de la forma
A1 A2,1 A2,2 A3,1 A3,2
f (z) = + + 2
+ + .
(z + 1) (z − 2) (z − 2) (z − 3) (z − 3)2
Se tiene que
A1 = Res [f (z)] = 1,
z=−1
A2,2 = (z − 2)2 f (z) z=2 = 336,
d 2
A2,1 = [(z − 2) f (z)] = 800,
dz z=2
A3,2 = (z − 3)2 f (z) z=3 = 468,
d 2
A3,1 = [(z − 3) f (z)] = −801,
dz z=3
CAPÍTULO 4. INTEGRACIÓN COMPLEJA 119
z+1 z5
f (z) = sen (z + 1), g(z) = , h(z) = .
(z − 2)2 ez−1
>> syms z
>> f ( z ) = sin ( z +1) ; int (f , z )
ans ( z ) =
- cos ( z + 1)
>> g ( z ) = ( z +1) /( z -2) ^2; int (g , z )
ans ( z ) =
log ( z - 2) - 3/( z - 2)
>> h ( z ) = z ^5/ exp (z -1) ; int (h , z )
ans ( z ) =
- exp (1 - z ) *( z ^5 + 5* z ^4 + 20* z ^3 + 60* z ^2 + 120* z + 120)
>> syms z
>> f ( z ) = ( exp ( z ) + z ^3 - i ) / exp (z -1) ; int (f ,z , -1 -i , -1)
ans ( z ) =
exp (1) * i + exp (2) *( - 2 + i ) + cos (1) * exp (2) *(2 - 3* i )
+ exp (2) * sin (1) *(3 + 2* i )
120 4.7. INTEGRACIÓN Y RESIDUOS CON MATLAB
En MATLAB no existe una comando que permita calcular, en forma directa, una
integral de lı́nea. Nosotros hemos diseñado la función IntegralLinea (ver Listado 4.1),
que permite calcular la integral de lı́nea
Z Z b
f (z) dz = f (z(t)) z ′ (t) dt,
C a
donde C es una curva con parametrización z(t) = x(t) + iy(t), a ≤ t ≤ b y f (z) es una
función de variable compleja.
f u n c t i o n [ vi ]= I n t e g r a l L i n e a (f , zt ,a , b )
% I n t e g r a l L i n e a C a l c u l a el valor de la i n t e g r a l de f ( z ) a lo
% largo de la curva C con p a r a m e t r i z a c i ó n
% z ( t ) , a <= t <= b
% vi : valor de la i n t e g r a l de lı́nea
%
syms z t ;
gt = f ( zt ) ;
dzt = diff ( zt , t ) ;
vi = int ( gt * dzt ,t ,a , b ) ;
end
>> syms z t
>> f ( z ) = conj ( z ) ^3/ z ^2; zt =2*( cos ( t ) + i * sin ( t ) ) ;
>> vi = I n t e g r a l L i n e a (f , zt , pi /2 , pi )
vi =
0
>> f ( z ) = log ( z ) ; zt = t +(1 -2*t ) * i ;
>> vi = I n t e g r a l L i n e a (f , zt ,0 ,1)
vi =
pi *(1/4 - i /4) + log (2) *(1/2 - i /2) - 1 + 2* i
Por otra parte, MATLAB no cuenta con un comando que calcule directamente el residuo
de una función f (z) en z0 . Hemos diseñado la función ResiduoPolo (ver Listado 4.2), que
calcula el residuo Res [f (z)], cuando z0 es un polo de orden m. En el siguiente ejemplo, se
z=z0
utiliza la función ResiduoPolo para calcular los residuos:
ez
a) Res 2 2 = 1.
z=0 z (z + 1)
Log z
b) Res 4 = 0.
z=1 z (z − 1)2
CAPÍTULO 4. INTEGRACIÓN COMPLEJA 121
sen z − z π−senh(π)
c) Res = π .
z=πi z senh z
f u n c t i o n [ res ]= R e s i d u o P o l o(f , z0 , m )
% R e s i d u o P o l o C a l c u l a el r e s i u d o de f ( z ) en el polo
% z0 de orden m
% res : valor del r e s i d u o
%
syms z ;
phi = (z - z0 ) ^ m * f ;
if m ==1
res = limit ( phi ,z , z0 ) ;
else
dphi = diff ( phi ,z ,m -1) ;
res = limit ( dphi ,z , z0 ) / f a c t o r i a l(m -1) ;
end
end
>> syms z
>> f ( z ) = exp ( z ) /( z ^2*( z ^2+1) ) ; z0 = 0; m = 2;
>> res = R e s i d u o P o l o (f , z0 , m )
res ( z ) =
1
>> f ( z ) = log ( z ) /( z ^4*( z -1) ) ; z0 = 1; m = 2;
>> res = R e s i d u o P o l o (f , z0 , m )
res ( z ) =
0
>> f ( z ) = ( sin ( z ) -z ) /( z * sinh ( z ) ) ; z0 = pi * i ; m = 1;
>> res = R e s i d u o P o l o (f , z0 , m )
res ( z ) =
( pi - sinh ( pi ) ) / pi
Ahora bien, para calcular el residuo de una función f (z) en un punto singular esencial
z0 , es necesario hallar el desarrollo de Laurent para obtener tal residuo. Para ello podemos
usar convenientemente el comando taylor para determinar solo aquellos coeficientes del
desarrollo de Laurent necesarios para hallar Res [f (z)]. Tenga presente que el procedi-
z=z0
miento que explicaremos a continuación, se aplica a funciones f (z) que pueden expresarse
como f (z) = f1 (z)f2 (z) (o un producto finito de funciones), donde z0 es un punto singular
esencial de f2 (z) y f1 (z) es analı́tica en z0 . Suponga que se desea hallar el residuo de
f (z) = (1 + z 2 )e1/(z−2i) en z0 = 2i, que, por supuesto, es un punto singular esencial de
f (z). Es claro que
f (z) = f1 (z)f2 (z),
donde
f1 (z) = 1 + z 2 y f2 (z) = e1/(z−2i) .
Primero, hallamos el desarrollo de Taylor de f1 (z) centrado en z0 = 2i (en este caso se
hallan todos los términos, que son tres),
122 4.7. INTEGRACIÓN Y RESIDUOS CON MATLAB
>> syms z
>> f1 ( z ) = 1+ z ^2; dtf1 ( z ) = taylor ( f1 ,z ,2 i , ’ Order ’ ,3)
dtf1 ( z ) =
z *4* i + ( z - 2* i ) ^2 + 5
o, equivalentemente,
f1 (z) = −3 + 4i(z − 2i) + (z − 2i)2 . (4.4)
>> syms w
>> h ( w ) = exp ( w ) ; dth ( w ) = taylor (h ,w ,0 , ’ Order ’ ,5)
dth ( w ) =
w ^4/24 + w ^3/6 + w ^2/2 + w + 1
que corresponden a los primeros 5 términos del desarrollo de Laurent de f2 (z) centrado en
z0 = 2i,
1 1 1 1
+ 2 + 3 + +1 (4.5)
z − 2 i 2 (z − 2 i) 6 (z − 2 i) 24 (z − 2 i)4
Multiplicando las ecuaciones (4.4) y (4.5), obtenemos que el coeficiente b1 del desarrollo
de Laurent de f (z) centrado en z0 = 2i, está dado por:
4i 1 17
b1 = −3 + + = − + 2i,
2 6 6
en otras palabras, el residuo de f (z) en z0 = 2i es
17
Res [f (z)] = − + 2i.
z=2i 6
Para finalizar esta sección, mostramos la función ExpFracParc (ver Listado 4.3) que
halla la expansión en fracciones parciales de una función racional propia del tipo
b0 + b1 z + · · · + bM z M
f (z) = ,
a0 + a1 z + · · · + z N
donde M < N .
En el siguiente ejemplo, se utiliza la función ExpFracParc para hallar la expansión en
fracciones parciales de
z+1
f (z) = .
(z + 3)(z + 2)3 (z − 1)2
CAPÍTULO 4. INTEGRACIÓN COMPLEJA 123
>> syms z
>> f ( z ) =( z +1) /(( z +3) *( z +2) ^3*( z -1) ^2) ; p =[ -3 -2 1]; v =[1 3 2];
>> efp = E x p F r a c P a r c (f ,p , v )
efp =
1 / ( 5 4 * (z - 1) ^2) - 1 / ( 7 2 * (z - 1) ) - 1/(9*( z + 2) ) +
4 / ( 2 7 * (z + 2) ^2) + 1/(8*( z + 3) ) - 1/(9*( z + 2) ^3)
o, equivalentemente,
1/9 4/27 1/9 1/8 1/72 1/54
f (z) = − + 2 − 3 + − +
(z + 2) (z + 2) (z + 2) (z + 3) (z − 1) (z − 1)2
f u n c t i o n [ efp ]= E x p F r a c P a r c(f ,p , v )
% ExpFracParc Halla la e x p a n s i ó n en f r a c c i o n e s p a r c i a l e s
% de la f u n c i ó n r a c i o n a l propia f ( z )
% p : vector de polos
% v : vector de m u t i p l i c i d a d e s de los polos
% efp : e x p a n s i ó n en f r a c c i o n e s p a r c i a l e s
syms z ;
np = length ( p ) ;
efp = 0;
for t =1: np
if v ( t ) == 1
dphi ( z ) = (z - p ( t ) ) * f ( z ) ;
r = dphi ( p ( t ) ) ;
efp = efp + sym ( r /(z - p ( t ) ) ) ;
else
for s =1: v ( t )
dphi ( z ) = diff (( z - p ( t ) ) ^( v ( t ) ) * f ( z ) ,z , v ( t ) -s ) ;
r = dphi ( p ( t ) ) / f a c t o r i a l( v ( t ) -s ) ;
efp = efp + sym ( r /( z - p ( t ) ) ^ s ) ;
end
end
end
end
Problema 4.2. Sea f (z) una función de variable compleja definida por
(
1, si Re z < 0;
f (z) =
4 Re z, si Re z ≥ 0.
Problema 4.3. Sea f (z) = Log (z + 1). Determine la expresión analı́tica de la primitiva
Rz
F (z) = f (s) ds + c, en el dominio simplemente conexo
i−1
Problema 4.4. Sea f (z) una función de variable compleja definida como
f (z) = z i ≡ ei log z ,
donde log z = ln |z| + i arg z, con |z| > 0 y −π/2 < arg z < 3π/2. Determine el valor de la
integral Z
f (z) dz,
C
donde C es el arco de la circunferencia |z| = 1 desde z = −1 hasta z = 1 situado en el
semiplano superior.
Solución. La función f (z) = ei log z es analı́tica en el dominio simplemente conexo
D = {z ∈ C : |z| > 0, −π/2 < arg z < 3π/2}.
Por lo tanto, f (z) tiene primitiva en D dada por:
z i+1 e(i+1) log z
F (z) = ≡ + c,
i+1 i+1
además, la curva C está dentro de D. Luego, podemos escribir:
Z Z 1 " #1 0
i e(i+1) log z
e − e−π eπi 1 + e−π
f (z) dz = z dz = +c = = (1 − i).
C −1 i+1 i+1 2
−1
♦
Problema 4.5. Determine el valor de la integral
Z
sen (z + 1)
34
dz.
|z+1|=1 (z + 1)
Solución. Sea f (z) = sen (z + 1). Es claro que f (z) es analı́tica en la circunferencia
|z + 1| = 1 y en el interior de la misma; además, z0 = −1 es un punto interior de esta
circunferencia. Ası́, usando la extensión de la fórmula integral de Cauchy, podemos escribir:
Z
sen (z + 1) 2πi f (33) (−1)
34
dz = .
|z+1|=1 (z + 1) 33!
Hallemos f (33) (−1). Para ello utilizaremos el desarrollo de Taylor de f (z) centrado en
z0 = −1. Tal desarrollo tiene la forma
∞
X
f (z) = an (z + 1)n , |z + 1| < ∞,
n=0
f (n) (−1)
donde an = n! . Ası́, f (33) (−1) = 33! a33 . Como
∞
X (−1)n
f (z) = sen (z + 1) = (z + 1)2n+1 , |z + 1| < ∞,
n=0
(2n + 1)!
entonces f (33) (−1) = 1. Por lo tanto, el valor de la integral dada está dado por:
Z
sen (z + 1) 2πi
34
dz = .
|z+1|=1 (z + 1) 33!
♦
126 4.8. PROBLEMAS RESUELTOS
luego,
1
Res [f1 (z)] = − i.
z=z0 2
Por otra parte, como z1 = 1 es un polo simple de f2 (z), se tiene que
z3 1
Res [f2 (z)] = 2
= .
z=z1 (z + 2) z=1 9
Ası́,
f (z) = f1 (z) + f2 (z) + f3 (z),
luego,
Z Z Z Z
2 1/z 1/(z−1) 1/(z−2)
(1 + z + z )(e + e +e ) dz = f1 (z)dz + f2 (z)dz + f3 (z)dz.
C C C C
Ahora, los únicos puntos singulares de f1 (z), f2 (z) y f3 (z) son, respectivamente, 0, 1 y 2,
además, cada uno de ellos son puntos singulares esenciales, lo cual nos obliga a determinar
el desarrollo de Laurent para calcular el residuo en tales puntos. Se tiene que el desarrollo
de Laurent de f1 (z) centrado 0, válido en el anillo |z| > 0, está dado por
∞ ∞ ∞ ∞
X 1 −n X 1 −n X 1 1−n X 1 2−n
f1 (z) = (1 + z + z 2 ) z = z + z + z ,
n=0
n! n=0
n! n=0
n! n=0
n!
128 4.9. PROBLEMAS PROPUESTOS
que implica
3 1 28
+ =
Res [f2 (z)] = 3 + .
z=1 2 6 6
Ahora, el desarrollo de Laurent de f3 (z) centrado 2, válido en el anillo |z − 2| > 0, está
dado por
∞
X 1
f3 (z) = (7 + 5(z − 2) + (z − 2)2 ) (z − 2)−n
n!
n=0
∞ ∞ ∞
X 7 X 5 X 1
= (z − 2)−n + (z − 2)1−n + (z − 2)2−n ,
n! n! n!
n=0 n=0 n=0
luego
5 1 58
+ =
Res [f2 (z)] = 7 + .
z=1 2 6 6
De esta forma, utilizando el Teorema de los Residuos obtenemos
Z Z Z
10 28 58
f1 (z)dz = πi, f2 (z)dz = πi, f3 (z)dz = πi,
C 3 C 3 C 3
en consecuencia, el valor de la integral dada es
Z
10 28 58
(1 + z + z 2 )(e1/z + e1/(z−1) + e1/(z−2) ) dz = πi + πi + πi = 32πi.
C 3 3 3
♦
Z
4.3. Evalúe la integral Log z dz, donde C es la curva z(t) = eit , 0 ≤ t ≤ π.
C
Z
Entonces, evalúe la integral f (z) dz.
|z|=2
4.5. Para cada una de las siguientes integrales, diga las caracterı́sticas del contorno
cerrado simple C para que el valor de la integral sea cero, según el Teorema de
Cauchy-Goursat. (Justifique su respuesta.)
Z Z
cos z 1
a) dz c) z
dz
C z+2 C 1+e
Z Z
1
b) Log z dz d) dz
C C 1 − ez
1
4.6. Para la función de f (z) = determine la primitiva de f (z) tal que:
(z − i)2
4.7. Evalúe las siguientes integrales. Use la fórmula integral de Cauchy (o su extensión)
o bien el Teorema de Cauchy-Goursat donde sea necesario. En cada una de las
integrales, C es un contorno cerrado simple.
Z
z
a) 4
dz, a > 1.
|z−a|=a z − 1
Z
1 x2 y2
b) dz, donde C es el contorno + = 1.
C ez (z − 2) 9 16
Z
sen (ez + cos z) x2
c) dz, donde C es el contorno + y 2 = 1.
C (z − 1)2 (z + 3) 2
Z
cos(2z)
d) dz.
|z|=1 z 21
Z
ez
e) dz, donde C es un contorno cerrado simple que contiene a |z| ≤ a.
C z2 + a2
Z
z ez
f) dz, donde a es un punto interior del contorno cerrado simple C.
C (z − a)3
130 4.9. PROBLEMAS PROPUESTOS
4.8. En cada una de las siguientes integrales determine el mı́nimo valor de la constante
positiva M que satisface Z
f (z) dz ≤ M,
C
para cada función f (z) considerada en cada caso.
Z
1
a) dz
z
|z|=1
Z
1
b) dz
z2 +1
|z−i|=1/2
Z
1
c) dz
ez −1
|z|=2
Z
z
d) dz, donde C está conformado con los puntos del cuadrado con vértices
C z−i
en los puntos −1 + 2i, 1 + 2i, 1 y −1.
4.9. Determine los residuos de las siguientes funciones complejas en cada polo o en el
polo indicado:
2z + 1 z+1
a) f (z) = . e)
z2 − z − 2 (z − 1)2 (z + 3)
cos z
1 f) , z0 = 0
b) f (z) = 2 . z
z (1 − z)
z4 − 1 π
g) 4
, z0 = e 4 i
3z 2 + 2 z +1
c) f (z) = .
(z − 1)(z 2 + 9) sen z π
h) , z0 = e 3 i
z4 2
+z +1
z3 − z2 + z + 1 z
d) i) , z0 = π
z 2 + 4z sen z
4.10. Use el Teorema de los Residuos para evaluar las siguientes integrales:
Z
z
a) 2
dz, para
z +1
C (
i) |z| = 1/2
C=
ii) |z| = 2
Z
z 2 + 3i − 2
b) dz, para
z 3 + 9z
C (
i) |z| = 1
C=
ii) |z| = 4
CAPÍTULO 4. INTEGRACIÓN COMPLEJA 131
Z
(z 2 + 2)(z 2 + 4)
c) dz para
(z 2 + 1)(z 2 + 6)
C
i) |z| = 2
C= ii) |z − i| = 1
iii) |z| = 4
Z
1
d) dz, para
z 2 (1 + z 2 )2
C
(
i) |z| = 1/2
C=
ii) |z| = 2
Z
3z 2 + 2
e) dz, para
(z 2 + 4)(z − 1)
C
(
i) |z − 2| = 2
C=
ii) |z| = 4
Z
z 2 − 2z
f) dz, para
(z 2 + 4)(z − 1)2
C
(
i) |z| = 3
C=
ii) |z + i| = 2
Z
1
g) dz, para
(z + 1)3 (z − 1)(z − 2)
C
(
i) |z + 1| = 1
C=
ii) es el rectángulo de vértices en ± i, 3 ± i
Parte II
Cálculo Operacional
132
5
Funciones de Dominios Continuo y
Discreto
En este capı́tulo se da una introducción al cálculo operacional, comenzando con la
caracterización de las funciones de dominio continuo, pasando luego por la convolución
en el dominio continuo. En esta parte se estudia con interés la función delta de Dirac o
impulso unitario. Finalmente, se describen las funciones de dominio discreto y se define la
convolución de funciones de dominio discreto.
Note que normalmente en los libros textos de señales y sistemas, a una función de
dominio continuo se le denomina señal continua o señal en tiempo continuo. Tal deno-
minación corresponde a que la representación matemática de una señal continua es una
función definida en el conjunto de números reales. Recuerde que una señal es cualquier
fenómeno o magnitud fı́sica que puede ser representado de manera cuantitativa. En este
material no se usa el concepto de señal, pero el lector debe tener presente la relación que
existe entre señal en tiempo continuo y funciones de dominio continuo. A continuación se
describen las funciones de dominio continuo de uso más frecuente.
133
134 5.1. FUNCIONES DE DOMINIO CONTINUO
δ(t)
Aδ(t − t0 )
t0 t
i) lı́m δ(t) = ∞.
t→0
iv) δ(t) = δ(−t), en otras palabras, δ(−t) se comporta como δ(t) (si δ(t) se toma
como una función ordinaria, esta propiedad indicarı́a que el impulso unitario es
una función par).
1
δ(at) = δ(t).
|a|
Ejercicio 5.1. Demostrar cada una de las propiedades del impulso unitario.
a) 3t4 δ(t − 1)
b) tδ(3 − 2t)
Z ∞
c) sen (t2 ) δ′ (t − π 1/2 ) dt
−∞
136 5.1. FUNCIONES DE DOMINIO CONTINUO
b) Se tiene que
3
tδ(3 − 2t) = t δ −2 t − .
2
Ahora, aplicando la propiedad de escalado seguido de la propiedad de muestreo, se obtiene:
1 3 t 3 3 3
tδ(3 − 2t) = t δ t− = δ t− = δ t− .
| − 2| 2 2 t=3/2 2 4 2
u(t)
1
u(t − t0 )
t0 t
(
1, |t| < T,
pT (t) =
0, |t| > T.
pT (t)
1
−T T t
|t|
1− , |t| ≤ T,
qT (t) = T
0, |t| > T,
qT (t)
1
−T T t
(
−1, t < 0,
sgn(t) =
1, t > 0,
sgn(t)
t
−1
Exp (t)
A
r(t)
2
1 2 t
Sa (t)
t
−6π −5π −4π −3π −2π −π π 2π 3π 4π 5π 6π
sinc (t)
t
−6 −5 −4 −3 −2 −1 1 2 3 4 5 6
se comporta como el escalón unitario, en otras palabras, verificaremos que h(t) = 0, para
t < 0 y h(t) = 1, para t > 0. Como δ(t) = 0 para todo t 6= 0, entonces es claro que
h(t) = 0, para t < 0. Ahora, para t > 0 podemos escribir:
Z t Z 0− Z 0+ Z t
h(t) = δ(τ ) dτ = δ(τ ) dτ + δ(τ ) dτ + δ(τ ) dτ,
−∞ −∞ 0− 0+
luego
h(t) = 1, para t > 0.
Por lo tanto, todo lo anterior nos indica que h(t) se comporta como el escalón unitario, es
decir,
Z t
u(t) = δ(τ ) dτ.
−∞
d
δ(t) = u(t).
dt
Seguidamente mostramos otras relaciones entre funciones comunes. La demostración
de tales relaciones se deja como ejercicio para el lector.
1 1
sgn(t) = −1 + 2 u(t) y u(t) = + sgn(t).
2 2
n
′ d X
f (t) = f (t)
+ hk δ(t − tk ). (5.1)
dt t6=tk k=1
La expresión de la derivada generalizada de f (t) dada arriba, se deduce del hecho que
toda función diferenciable a trozos o, en general, toda función continua a trozos puede
expresarse como una suma de funciones que involucran escalones unitarios. Observe el
siguiente ejemplo.
f (t)
3
-1 1 2 3 4 t
donde
−t
−e , −1 < t < 0,
g(t) = 1, 0 < t < 2,
0, en otros casos.
Determinemos la derivada generalizada de f (t) por otro procedimiento, usando el hecho
que f (t) puede expresarse equivalentemente como:
Recuerde la
regla de la f (t) = e−t (u(t + 1) − u(t)) + t(u(t) − u(t − 2)) + 3(u(t − 2) − u(t − 3)) + u(t − 3).
derivada del
producto
Ahora, derivando la expresión anterior se tiene:
d −t d d
f ′ (t) = e (u(t + 1) − u(t)) + [t(u(t) − u(t − 2))] + [3(u(t − 2) − u(t − 3))]
dt dt dt
d
+ u(t − 3)
dt
= − e−t (u(t + 1) − u(t)) + e−t (δ(t + 1) − δ(t)) + (u(t) − u(t − 2))
+ t(δ(t) − δ(t − 2)) + 3(δ(t − 2) − δ(t − 3)) + δ(t − 3)
= − e−t (u(t + 1) − u(t)) + (u(t) − u(t − 2)) + e δ(t + 1) − δ(t) + δ(t − 2) − 2δ(t − 3),
Propiedades
Para funciones f , g y h de dominio continuo se satisface:
i) Conmutatividad
f ∗g =g∗f
f ∗ (g + h) = f ∗ g + f ∗ h
u(t − τ )
-2 t -1 1 2 τ
u(t − τ )
-2 -1 t 1 2 τ
Z t
Entonces, por (5.3), h(t) = dτ = t + 1, para −1 < t < 1.
−1+
u(t − τ )
1
-2 -1 1 t 2 τ
Z 1−
Entonces, por (5.3), h(t) = dτ = 2, para t ≥ 1.
−1+
h(t)
2
-2 -1 1 2 t ♦
1 1
-1 1 t -1 1 t
g(t − τ )
1
t−1 t 1 2 τ
g(t − τ )
1
t−1 t 1 2 τ
CAPÍTULO 5. FUNCIONES DE DOMINIOS CONTINUO Y DISCRETO 145
g(t − τ )
1
t−1 1 t 2 τ
g(t − τ )
1
1t−12 t τ
De esta forma, la expresión analı́tica de la convolución h(t) = f ∗ g (t) está dada por:
0, t ≤ 0,
1 1 − (t − 1)2 , 0 < t < 1,
2
h(t) =
1 2
2 (t − 2) , 1 ≤ t < 2,
0, t ≥ 2,
cuya representación gráfica se aprecia en la siguiente figura.
h(t)
1
1/2
-3 -2 -1 1 2 3 t
♦
146 5.2. FUNCIONES DE DOMINIO DISCRETO
δ(n)
b
1
b b b b b b b b
-4 -3 -2 -1 1 2 3 4 n
u(n)
1 b b b b b
···
b b b b
-4 -3 -2 -1 1 2 3 4 n
2 b
···
1 b
b b b b b
-4 -3 -2 -1 1 2 3 n
sgn(n)
1 b b b b
···
−4 −3 −2 −1 b
1 2 3 4 n
···
b b b b
-1
148 5.2. FUNCIONES DE DOMINIO DISCRETO
se comporta como el escalón unitario discreto. Como δ(k) = 0 para todo k 6= 0, entonces
es claro que f (n) = 0 para n < 0. Para n = 0, se tiene que
−1
X
f (0) = δ(k) + δ(0) = 1.
k=−∞
Por lo tanto, todo lo anterior nos dice que f (n) se comporta como el escalón unitario
discreto, es decir,
Xn
u(n) = δ(k).
k=−∞
Veamos ahora que la función g(n) = u(n) − u(n − 1) se comporta como el impulso
unitario discreto. Como u(n) = 0 para todo n < 0, entonces es claro que g(n) = 0 para
n < 0. Para n = 0, se tiene que
Ahora, se tiene que u(n) = 1 y u(n − 1) = 1 para n ≥ 1; luego, g(n) = 0 para n ≥ 1. Por lo
tanto, todo lo anterior nos indica que g(n) se comporta como el impulso unitario discreto,
es decir,
δ(n) = u(n) − u(n − 1).
∞
X
f ∗ g (n) = f (k) g(n − k), para n ∈ Z.
k=−∞
Propiedades
Para funciones f , g y h de dominio discreto se satisface:
i) Conmutatividad
f ∗g =g∗f
f ∗ (g + h) = f ∗ g + f ∗ h
Solución. Considerando las expresiones de f (n) y g(n), se tiene que la expresión analı́tica
de la convolución h(n) = f ∗ g (n), para cada n ∈ Z, es:
∞
X
h(n) = f ∗ g (n) = f (k)g(n − k)
k=−∞
3
X
= g(n − k)
k=2
= g(n − 2) + g(n − 3)
= u(n − 2) + u(n − 3),
2 b b b b
b
1
··· ···
b b b b b
-3 -2 -1 1 2 3 4 5 6 n
♦
150 5.2. FUNCIONES DE DOMINIO DISCRETO
Ejemplo 5.8. Si f (n) = (1/3)n u(n) y g(n) = 2n u(−n), entonces determine la convolución
f ∗ g.
u(k − n)
1
b b b b b
··· ···
b b b
n−1 n −1 1 2 k
u(k − n)
1 b b b
··· ···
b b b b b b
De esta forma, la expresión analı́tica de la convolución h(n) = f ∗ g (n) está dada por
h(n) = (6/5)2n u(−1 − n) + (6/5)3−n u(n),
cuya representación gráfica se aprecia en la siguiente figura.
h(n)
b
1
b
b
b
b b
b b b b b
-5 -4 -3 -2 -1 1 2 3 4 5 n
>> syms t
>> f ( t ) = t ^2*( sin ( t ) + exp ( - t ) )
f(t) =
t ^2*( exp ( - t ) + sin ( t ) )
2
t (sin(t)+exp(−t))
1000
800
600
400
200
−5 −4 −3 −2 −1 0 1 2 3 4 5
t
152 5.3. INTRODUCCIÓN AL CÁLCULO OPERACIONAL CON MATLAB
Las funciones simbólicas creadas en MATLAB se pueden operar entre sı́, diferenciar,
integrar, etc. Ahora, las funciones elementales ya están predefinidas y las funciones como
el impulso unitario y el escalón unitario, también están predefinidas. El impulso unitario,
δ(t), es dirac(t), y el escalón unitario, u(t), es heaviside(t). La gráfica de u(t) en el
intervalo [−10, 10], se obtiene de la siguiente forma:
>> syms t
>> ezplot ( h e a v i s i d e ( t ) ,[ -10 ,10])
heaviside(t)
0.8
0.6
0.4
0.2
Utilice el 0
comando
ezplot para −10 −8 −6 −4 −2 0 2 4 6 8 10
obtener las t
gráficas de las
funciones
sgn(t), r(t), Las otras funciones de dominio continuo de uso común, por ejemplo, la función signo
p1 (t) y q1 (t), en sgn(t), la rampa r(t), el pulso rectangular pT (t) o el pulso triangular qT (t), se pueden
el intervalo
[−10, 10]
definir a través de su relación con el escalón unitario. Observe el siguiente ejemplo donde
se definen las funciones sgn(t), r(t), p1 (t) y q1 (t).
>> syms t
>> u ( t ) = h e a v i s i d e ( t ) ;
>> sgn ( t ) = u ( t ) -u ( - t )
sgn ( t ) =
h e a v i s i d e ( t ) - h e a v i s i d e( - t )
>> r ( t ) = t * u ( t )
r(t) =
t * h e a v i s i d e( t )
>> p1 ( t ) = u ( t +1) -u (t -1)
p1 ( t ) =
h e a v i s i d e ( t + 1) - h e a v i s i d e( t - 1)
>> q1 ( t ) = (1+ t ) *( u ( t +1) -u ( t ) ) +(1 - t ) *( u ( t ) -u (t -1) )
q1 ( t ) =
( t - 1) *( h e a v i s i d e( t - 1) - h e a v i s i d e( t ) ) +
( t + 1) *( h e a v i s i d e( t + 1) - h e a v i s i d e( t ) )
>> syms t
>> f ( t ) = h e a v i s i d e ( cos ( t ) ) ;
>> ezplot (f ,[ -6* pi ,6* pi ])
heaviside(cos(t))
0.8
0.6
0.4
0.2
−15 −10 −5 0 5 10 15
t
Con el comando int se pueden calcular integrales que involucren a la función impulso
unitario. Seguidamente se calculan los valores de las integrales
Z 3 Z 3
(t−2) 1 1
e δ(2t − 4) dt = , e(t−2) δ′′ (2t − 4) dt = .
0 2 0 2
>> syms t
>> f1 ( t ) = exp (t -2) * dirac (2*t -4) ;
>> int ( f1 ,t ,0 ,3)
ans ( t ) =
1/2
>> f2 ( t ) = exp (t -2) * diff ( dirac (2*t -4) ,t ,2) ;
>> int ( f2 ,t ,0 ,3)
ans ( t ) =
1/2
f u n c t i o n [ h0 ]= C o n v C o n t(f ,g , t0 )
% ConvCont Halla el valor de la c o n v o l u c i ó n h ( t ) = f * g ( t )
% en t0
% h0 = h ( t0 )
syms t tau ;
h0 = int ( f ( tau ) * g ( t0 - tau ) , tau , - Inf , Inf ) ;
end
>> syms t
>> f ( t ) = h e a v i s i d e ( t +1) - h e a v i s i d e(t -1) ;
>> g ( t ) = heaviside(t);
>> h0 = C o n v C o n t(f ,g ,1/2)
h0 =
3/2
>> t0 = -4:.1:4;
>> h0 = C o n v C o n t(f ,g , t0 ) ;
>> plot ( t0 , h0 , ’ LineWidth ’ ,2)
>> axis ([ -4 ,4 , -1 ,3]); grid on
>> title ( ’ h ( t ) = f * g ( t ) ’)
h(t) = f*g(t)
3
2.5
1.5
0.5
−0.5
−1
−4 −3 −2 −1 0 1 2 3 4
En este caso, t0 es un vector, por ello ConvCont arroja como salida un vector h0, cuyas
componentes son los valores que alcanza h(t) en cada componente del vector t0.
CAPÍTULO 5. FUNCIONES DE DOMINIOS CONTINUO Y DISCRETO 155
En cuanto a las funciones de dominio discreto, con MATLAB también se pueden ma-
nejar, pero en forma limitada. Una función f (n) de dominio discreto es en sı́ una sucesión
de números reales. Este hecho nos permite definir funciones de dominio discreto como
un arreglo de números reales. Por ejemplo, la función sen (πn/2) puede definirse en MA-
TLAB como una discretización de la función sen (πt), primero creando n como un arreglo
de números enteros y después evaluando la función sen (πt/2) en ese arreglo. Esto puede
observarlo a continuación.
>> n = -10:10;
>> f = sin ( pi * n /2) ;
sen(nπ/2)
1
0.8
0.6
0.4
0.2
−0.2
−0.4
−0.6
−0.8
−1
−10 −8 −6 −4 −2 0 2 4 6 8 10
El proceso anterior se puede aplicar para generar todas las funciones elementales de
dominio discreto (cos(n), en , tan(n), etc.). En cuanto a las funciones impulso unitario
discreto δ(n) y escalón unitario discreto u(n), hemos definido las funciones impulso (ver
Listado 5.2) y escalon (ver Listado 5.3), que hallan los valores que alcanzan δ(n) y u(n)
en las componentes del arreglo n, respectivamente. A continuación usamos las funciones
impulso y escalon, para generar la gráfica de δ(n) y la de u(n) en el intervalo [−10, 10].
>> n = -10:10;
>> s u b p l o t (2 ,1 ,1) , stem (n , i m p u l s o( n ) ,’ fill ’ , ’ - - ’)
>> title ( ’\ delta ( n ) ’)
>> s u b p l o t (2 ,1 ,2) , stem (n , e s c a l o n( n ) ,’ fill ’ , ’ - - ’)
>> title ( ’ u ( n ) ’)
156 5.3. INTRODUCCIÓN AL CÁLCULO OPERACIONAL CON MATLAB
f u n c t i o n [ d ] = i m p u l s o( n )
% impulso Halla los v a l o r e s de delta ( n )
% d : a r r e g l o de v a l o r e s de delta ( n )
d = zeros ( size ( n ) ) ;
d ( n ==0) = 1;
end
f u n c t i o n [ u ] = e s c a l o n( n )
% escalon Halla los v a l o r e s de u ( n )
% u : a r r e g l o de v a l o r e s de u ( n )
u = zeros ( size ( n ) ) ;
u (n >=0) = 1;
end
δ(n)
1
0.8
0.6
0.4
0.2
0
−10 −8 −6 −4 −2 0 2 4 6 8 10
u(n)
1
0.8
0.6
0.4
0.2
0
−10 −8 −6 −4 −2 0 2 4 6 8 10
Los Listados 5.4, 5.5 y 5.6, muestran las funciones signo, rampa y pulsorectangular,
que definen respectivamente las funciones signo discreto,
−1, n < 0,
sgn(n) = 0, n = 0,
1, n > 0,
rampa discreta,
r(n) = n u(n),
CAPÍTULO 5. FUNCIONES DE DOMINIOS CONTINUO Y DISCRETO 157
Utilice el
y pulso rectangular discreto comando stem
para obtener las
gráficas de las
0,
n < n1 ,
funciones
pn1 ,n2 (n) = 1, n1 ≤ n ≤ n2 , sgn(n), r(n) y
p−7,−1 (n), en el
0, n > n2 , intervalo
[−10, 10]
donde n1 , n2 ∈ Z tales que n1 < n2 .
f u n c t i o n [ s ] = signo ( n )
% signo Halla los v a l o r e s de sgn ( n )
% s : a r r e g l o de v a l o r e s de sgn ( n )
s = zeros ( size ( n ) ) ;
s (n <0) = -1;
s (n >0) = 1;
end
f u n c t i o n [ r ] = rampa ( n )
% rampa Halla los v a l o r e s de r ( n )
% r : a r r e g l o de v a l o r e s de r ( n )
r = zeros ( size ( n ) ) ;
r (n >=0) = n (n >=0) ;
end
f u n c t i o n [ p ] = p u l s o r e c t a n g u l a r (n , n1 , n2 )
% pulsorectangular Halla los v a l o r e s de un pulso
% r e c t a n g u l a r entre n1 y n2
% p : a r r e g l o de v a l o r e s del pulso r e c t a n g u l a r
p = ones ( size ( n ) ) ;
p (n < n1 ) = 0;
p (n > n2 ) = 0;
end
>> n = -10:10;
>> f1 = e s c a l o n(n -4) ;
>> f2 = n .*( e s c a l o n( n +2) - e s c a l o n(n -3) ) ;
>> f3 = e s c a l o n( - n ) - i m p u l s o(n -1) - i m p u l s o(n -2) ;
>> f4 = p u l s o r e c t a n g u l a r (n ,0 ,4) - p u l s o r e c t a n g u l a r (n , -5 , -1);
f1(n) f2(n)
1 2
0.8
1
0.6
0
0.4
−1
0.2
0 −2
−10 −5 0 5 10 −10 −5 0 5 10
f (n) f (n)
3 4
1 1
0.5 0.5
0 0
−0.5 −0.5
−1 −1
−10 −5 0 5 10 −10 −5 0 5 10
Utilicemos el comando h = conv(f,g) para hallar los valores de h(n), por supuesto, en el
intervalo [−2, 13].
f(n) g(n)
1 1
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0 1 2 3 4 5 −2 0 2 4 6 8
h(n)
2
1.5
0.5
0
−2 0 2 4 6 8 10 12 14
Ahora bien, considerando las expresiones de f (t) y g(t) dadas en el enunciado del problema,
obtenemos:
Z ∞ Z 2a−
h(t) = pa (τ − a)u(t − τ ) dτ = u(t − τ ) dτ, para t ∈ R. (5.6)
−∞ 0+
u(t − τ )
τ
t 2a
u(t − τ )
1
τ
t 2a
u(t − τ )
1
τ
2a t
Problema 5.2. Sean f (t) y g(t) funciones de dominio continuo definidas respectivamente
por:
f (t) = r(t) y g(t) = sgn(t) + u(−t − 1).
Determinar la convolución h(t) = f ∗ g(t).
Calculemos cada una de las convoluciones r(t) ∗ sgn(t) y r(t) ∗ u(−t − 1).
Usando la definición de convolución obtenemos:
Z ∞
r(t) ∗ sgn(t) = r(τ ) sgn(t − τ ) dτ, para t ∈ R. (5.7)
−∞
r(τ ) sgn(t − τ )
−τ
r(τ ) sgn(t − τ )
t τ
−t
−τ
162 5.4. PROBLEMAS RESUELTOS
Problema 5.3. Sean f (t) y g(t) funciones de dominio continuo definidas respectivamente
por:
f (t) = e−t u(t) y g(t) = e−2t u(t).
Determinar la convolución h(t) = f ∗ g(t).
eτ u(t − τ )
et
τ
t
eτ u(t − τ )
et
τ
t
t
♦
f (t) g(t)
1
−1 1
−1 1 t t
−1
(τ + 1)(t − τ )
t−1 t t+1
τ t−1 t t+1 τ
−1 −1
−1
(τ + 1)(τ − t)
(τ + 1)(τ − t)
t−1 t t+1 t t+1
−1 1 τ t−1 1 τ
−1 −1 (τ + 1)(t − τ )
(τ + 1)(t − τ )
CAPÍTULO 5. FUNCIONES DE DOMINIOS CONTINUO Y DISCRETO 165
−1
(τ + 1)(τ − t)
(τ + 1)(t − τ )
−1
(τ + 1)(τ − t)
(τ + 1)(t − τ )
6 2 2 0 < t ≤ 1,
t (t−2) ,
1 < t ≤ 2,
2
0, en otros casos.
166 5.4. PROBLEMAS RESUELTOS
h(t)
−2 2
t
−2/3
Problema 5.5. Sean f (n) y g(n) funciones de dominio discreto definidas respectivamente
por:
f (n) = u(n + 1) − u(n − 2) y g(n) = u(n + 1) − u(n − 2).
Determinar la convolución h(n) = f ∗ g(n).
Calculemos los valores que toma h(n) para cada n ∈ Z. Seguidamente construiremos
gráficas de g(n − k) para cada n. Aquı́ solo se mostrará el gráfico de
g(n − k)
b b b
1
g(n − k)
b b b
1
-5 -4 -3 -2 -1 1 2 3 4 5 k
g(n − k)
b b
1 b
-5 -4 -3 -2 -1 1 2 3 4 5 k
g(n − k)
b
1 b b
-5 -4 -3 -2 -1 1 2 3 4 5 k
g(n − k)
1 b b b
-5 -4 -3 -2 -1 1 2 3 4 5 k
g(n − k)
1 b b b
-5 -4 -3 -2 -1 1 2 3 4 5 k
g(n − k)
1 b b b
De todo lo anterior se tiene que la convolución h(n) = f ∗ g(n) está dada por:
h(n)
3 b
b b
2
b
1 b
b b b b b b
-5 -4 -3 -2 -1 1 2 3 4 5 n
Problema 5.6. Sean f (n) y g(n) funciones de dominio discreto definidas respectivamente
por:
Solución. Considerando las expresiones de f (n) y g(n) dadas en el enunciado del problema
y la definición de la convolución de funciones de dominio discreto, podemos escribir para
n ∈ Z:
∞
X 3
X
h(n) = f ∗ g(n) = f (k)g(n − k) = (n − k) u(n − k).
k=−∞ k=−2
Calculemos los valores que toma h(n) para cada n ∈ Z. Para ello construiremos la gráfica
de ku(k) y usaremos la ecuación (5.11) variando los valores de n.
La gráfica de k u(k) es:
CAPÍTULO 5. FUNCIONES DE DOMINIOS CONTINUO Y DISCRETO 169
ku(k)
6 b
5 b
4 b
b
3
2 b
1 b
b b b b b b
-5 -4 -3 -2 -1 1 2 3 4 5 6 k
..
.
De todo lo anterior se tiene que la convolución h(n) = f ∗ g(n) está dada por:
0, n ≤ −2,
(n+2)(n+3)
h(n) = 2 , −1 ≤ n ≤ 4,
(n+2)(n+3) (n−4)(n−3)
2 − 2 , n ≥ 5.
170 5.4. PROBLEMAS RESUELTOS
h(n) b
···
b
b
b b b b
-5 -4 -3 -2 -1 1 2 3 4 5 6 n
♦
4
X (−1)m
Problema 5.7. Determinar la convolución f ∗ g, donde f (n) = δ(n − m) y
m
m=1
g(n) es una función cuya gráfica se muestra en la siguiente figura.
g(n)
b b
1
1 2 3 4 ···
b b b b b
· · · −4 −3 −2 −1 n
-1 b b
Verifique que Solución. Se tiene que la función g(n) se puede expresar como g(n) = g1 (n) − g2 (n),
g(n) =
g1 (n) − g2 (n),
donde
para las
funciones g1 (n) g1 (n) = u(−1 − n) − u(−3 − n) y g2 (n) = u(n − 1) − u(n − 3).
y g2 (n) dadas
· · · −3 −2 −1 b 4 6
b b b b
b
b
1 2 3 5 b
7 8··· n
b
-1
f (t) g(t)
1 1
a) −1
−1 1 t 1 t
−1
f (t) g(t)
1 1
b) 1
−1 1 t −1 t
−1
f (t) g(t)
2
c) 1
1
−1 2 t −1 1 t
f (t) g(t)
d) 1 1
−1 1 t −1 1 t
CAPÍTULO 5. FUNCIONES DE DOMINIOS CONTINUO Y DISCRETO 173
5.5. Sea f (n) la función de dominio discreto cuya gráfica se muestra en la siguiente
figura.
f (n)
2 b
b
1
−2 1 3
−1 2 n
-1 b
b
-2 b
6.1 Definición
Sea f (t) una función de dominio continuo definida de R en C. La transformada de
Fourier de f (t) es una transformada integral lineal definida por:
Z ∞
F (ω) = F{f (t)} = f (t) e−jωt dt,
−∞
174
CAPÍTULO 6. TRANSFORMADA DE FOURIER 175
para todo ω0 ∈ R. En otras palabras, la multiplicación de la función f (t) por ejω0 t , resulta
en un desplazamiento en la transformada de Fourier por ω0 .
CAPÍTULO 6. TRANSFORMADA DE FOURIER 177
F 1
f (at) ←→ F (ω/a),
|a|
Para a < 0, haciendo el cambio de variable r = at, la integral de (6.1) adquiere la forma:
Z −∞ Z ∞
1 −j( ω )r 1 ω
F{f (at)} = f (r) e a dr = − f (r) e−j( a )r dr.
a ∞ a −∞
F
F (t) ←→ 2π f (−ω).
Demostración. Como Z ∞
1
f (t) = F (ω) ejωt dω,
2π −∞
178 6.2. PROPIEDADES DE LA TRANSFORMADA DE FOURIER
entonces Z ∞
2πf (−t) = F (ω) e−jωt dω.
−∞
Intercambiando t y ω, se obtiene
Z ∞
2πf (−ω) = F (t) e−jωt dt = F{F (t)}.
−∞
F
f (t) ←→ F (−ω).
F
f ∗ g (t) ←→ F (ω) G(ω).
F 1
f (t) g(t) ←→ F ∗ G (ω).
2π
Es decir, la multiplicación en el dominio del tiempo, corresponde a la convolución en el
dominio de la frecuencia multiplicada por la constante 1/2π.
Demostración. Se tiene que
Z ∞
F{f (t) g(t)} = (f (t) g(t)) e−jωt dt
−∞
Z ∞ Z ∞
1
= F (σ) e dσ g(t) e−jωt dt
jσt
−∞ 2π −∞
Z ∞ Z ∞
1
= F (σ) g(t) e−j(ω−σ)t dt dσ
2π −∞ −∞
Z ∞
1
= F (σ)G(ω − σ) dσ
2π −∞
1
= F ∗ G (ω).
2π
d F
f (t) ←→ jω F (ω).
dt
Demostración. Se tiene que
Z ∞
1
f (t) = F (ω) ejωt dω.
2π −∞
F d
t f (t) ←→ j F (ω).
dω
180 6.2. PROPIEDADES DE LA TRANSFORMADA DE FOURIER
F 1
cos(ω0 t) f (t) ←→ [F (ω + ω0 ) + F (ω − ω0 )].
2
y
F j
sen (ω0 t) f (t) ←→ [F (ω + ω0 ) − F (ω − ω0 )].
2
CAPÍTULO 6. TRANSFORMADA DE FOURIER 181
Demostración. Para la prueba usaremos las transformadas de Fourier de cos(ω0 t) y sen (ω0 t),
a saber:
F{cos(ω0 t)} = π[δ(ω − ω0 ) + δ(ω + ω0 )]
y
π
F{sen (ω0 t)} = [δ(ω − ω0 ) − δ(ω + ω0 )],
j
las cuales la hallaremos más adelante. Aplicando la propiedad de multiplicación
1
F{cos(ω0 t) f (t)} = F{cos(ω0 t)} ∗ F{f (t)}.
2π
Ahora, usando la expresión de la transformada de Fourier de cos(ω0 t) y la propiedad
distributiva de la convolución, se obtiene
1
F{cos(ω0 t) f (t)} = [δ(ω − ω0 ) ∗ F (ω) + δ(ω + ω0 ) ∗ F (ω)]
2
1
= [F (ω − ω0 ) + F (ω + ω0 )] .
2
Por un razonamiento similar se demuestra que
j
F{sen (ω0 t) f (t)} = [F (ω + ω0 ) − F (ω − ω0 )].
2
(Se deja al lector la prueba.)
Ahora, tomando f1 (t) = f (t) y f2 (t) = f (t), y luego aplicando la propiedad de conjugación
F
f (t) ←→ F (−ω), de la ecuación (6.4) se deduce
Z ∞ Z ∞ Z ∞
2 1 1
|f (t)| dt = F (ω)F (ω) dω = |F (ω)|2 dω,
−∞ 2π −∞ 2π −∞
que era lo que deseábamos demostrar.
1
Escalamiento en tiempo F {f (at)} = F (ω/a)
|a|
Dualidad F {F (t)} = 2πf (−ω)
n o
Conjugación F f (t) = F (−ω)
1
Multiplicación F {f (t)g(t)} = F ∗ G (ω)
2π
dn
Diferenciación en tiempo F f (t) = (jω)n F (ω)
dtn
dn
Diferenciación en frecuencia F {tn f (t)} = j n F (ω)
dω n
t
Z 1
Integración F f (λ) dλ = F (ω) + πF (0)δ(ω)
jω
−∞
1
Modulación F {cos(ω0 t) f (t)} = [F (ω + ω0 ) + F (ω − ω0 )]
2
j
F {sen (ω0 t) f (t)} = [F (ω + ω0 ) − F (ω − ω0 )]
2
Z ∞ Z ∞
1
Teorema de Parseval |f (t)|2 dt = |F (ω)|2 dω
−∞ 2π −∞
CAPÍTULO 6. TRANSFORMADA DE FOURIER 183
1.
F
δ(t) ←→ 1
Se tiene que
Z ∞
F{δ(t)} = δ(t) e−jωt dt
Z−∞
∞
= δ(t)(cos(−ωt) + jsen (−ωt)) dt
Z−∞
∞ Z ∞
= δ(t) cos(ωt) dt − j δ(t) sin(ωt) dt
−∞ −∞
= cos(0) − jsen (0) = 1.
F
En la Figura 6.1 se aprecia gráficamente el par de transformadas δ(t) ←→ 1.
δ(t) F (ω)
1
1 F
←→
t ω
2.
F
1 ←→ 2π δ(ω)
F
En la Figura 6.2 se aprecia gráficamente el par de transformadas 1 ←→ 2πδ(ω).
f (t) = 1 F (ω)
1 2πδ(ω)
F
←→
t ω
3.
F
ejω0 t ←→ 2π δ(ω − ω0 )
4.
F 2
sgn(t) ←→
jω
5.
F 1
u(t) ←→ + πδ(ω)
jω
sgn(t) = 2u(t) − 1,
de donde se deduce:
sgn(t) 1
u(t) = + .
2 2
Ası́, empleando la última ecuación, se obtiene:
1 1
F{u(t)} = F{sgn(t)} + F{1}
2 2
1 2 1 1
= + (2πδ(ω)) = + πδ(ω).
2 jω 2 jω
CAPÍTULO 6. TRANSFORMADA DE FOURIER 185
6.
F 1
e−αt u(t) ←→
α + jω
Se tiene que
Z ∞ Z ∞
−αt −αt −jωt 1
F{e u(t)} = e u(t) e dt = e−(α+jω)t dt = .
−∞ 0 α + jω
7.
F π
sen (ω0 t) ←→ [δ(ω − ω0 ) − δ(ω + ω0 )]
j
F
←→ ω0
t −ω0 ω
−π
8.
F
cos(ω0 t) ←→ π [δ(ω − ω0 ) + δ(ω + ω0 )]
cos(ω0 t) F (ω)
π π
F
←→
t −ω0 ω0 ω
9.
F sen (ωT )
pT (t) ←→ 2T
ωT
Como
pT (t) = u(t + T ) − u(t − T ),
entonces
pT (t) F (ω)
1 1
F
←→
−T T t ω
10.
F 2a
e−a|t| ←→ , para a > 0
ω 2 + a2
Como
e−a|t| = eat u(−t) + e−at u(t),
CAPÍTULO 6. TRANSFORMADA DE FOURIER 187
entonces
F{e−a|t| } = F eat u(−t) + e−at u(t)
Z ∞ Z ∞
at −jωt
= e u(−t) e dt + e−at u(t) e−jωt dt
−∞ −∞
Z 0 Z ∞
= eat e−jωt dt + e−at e−jωt dt
−∞ 0
Z 0 Z ∞
= e(a−jωt) dt + e−(a+jωt) dt
−∞ 0
1 1 2a
= + = 2 .
a − jω a + jω ω + a2
e−a|t| F (ω)
1 2/a
F
←→
t ω
11.
1 F π −a|ω|
←→ e , a>0
a2 +t 2 a
Usando las propiedades de linealidad y dualidad, se tiene
1 1 2a 1 −a|ω|
π
F = F = 2πe = e−a|ω| .
a2 + t 2 2a a2 + t2 2a a
1
En la Figura 6.7 se aprecia gráficamente la transformada de Fourier de .
a2 + t2
F (ω)
1
a2 +t2
π/a
F
1/a2 ←→
t ω
1
Figura 6.7. Transformada de Fourier de
a2 + t2
δ(t) 1
1
u(t) πδ(ω) +
jω
δ(t − t0 ) e−jωt0
2
sgn(t)
jω
ejω0 t 2πδ(ω − ω0 )
π
sen (ω0 t) [δ(ω − ω0 ) − δ(ω + ω0 )]
j
π jω
cos (ω0 t)u(t) [δ(ω − ω0 ) + δ(ω + ω0 )] + 2
2 ω0 − ω 2
π ω0
sen (ω0 t)u(t) [δ(ω − ω0 ) − δ(ω + ω0 )] + 2
2j ω0 − ω 2
1
e−at u(t), Re {a} > 0
a + jω
1
te−at u(t), Re {a} > 0
(a + jω)2
tn−1 −at 1
e u(t), Re {a} > 0
(n − 1)! (a + jω)n
2a
e−a|t| , a > 0
a2 + ω2
4ajω
|t|e−a|t| , Re {a} > 0
a2 + ω2
1 π −a|ω|
, Re {a} > 0 e
a2 + t2 a
t −jπωe−a|ω|
, Re {a} > 0
a + t2
2 2a
r
2 π −ω2 /4a
e−at , a > 0 e
a
CAPÍTULO 6. TRANSFORMADA DE FOURIER 189
Observación 6.2. Para calcular la fase es necesario fijar una determinación del valor del
arg (F (ω)). En este escrito la fase se calculará como el argumento principal, es decir,
φ(ω) = Arg (F (ω)).
La siguiente proposición es un resultado de mucha utilidad para el cálculo de la mag-
nitud y la fase de una función. Esta proposición establece que si f (t) es una función real,
entonces la magnitud es una función par y la fase es una función impar.
Proposición 6.1. Si f (t) es una función real, entonces el espectro de magnitud A(ω) es
una función par, y el espectro de fase φ(ω) es una función impar.
donde Z ∞
Re {F (ω)} = f (t) cos(ωt) dt (6.5)
−∞
e Z ∞
Im {F (ω)} = − f (t) sen (ωt) dt. (6.6)
−∞
Ası́, los espectros de magnitud y fase de f (t) están dados respectivamente por:
p
A(ω) = |F (ω)| = (Re {F (ω)})2 + (Im {F (ω)})2 (6.7)
190 6.4. MAGNITUD Y FASE
y
0, Re {F (ω)} > 0, Im {F (ω)} = 0,
Im {F (ω)}
arctan , Re {F (ω)} > 0, Im {F (ω)} > 0,
Re {F (ω)}
π/2, Re {F (ω)} = 0, Im {F (ω)} > 0,
Im {F (ω)}
arctan
+ π, Re {F (ω)} < 0, Im {F (ω)} > 0,
Re {F (ω)}
φ(ω) = (6.8)
π, Re {F (ω)} < 0, Im {F (ω)} = 0,
Im {F (ω)}
arctan − π, Re {F (ω)} < 0, Im {F (ω)} < 0,
Re {F (ω)}
−π/2, Re {F (ω)} = 0, Im {F (ω)} < 0,
Im {F (ω)}
arctan
, Re {F (ω)} > 0, Im {F (ω)} < 0.
Re {F (ω)}
Veamos que A(ω) es par. Como la función real | · | es una función par y, por (6.5) y (6.6),
Re {F (ω)} y Im {F (ω)} son funciones reales (por ser f (t) una función real), entonces por
(6.7) la función A(ω) es par.
Veamos ahora que φ(ω) es impar. Por (6.5) y (6.6) se tiene que
= −φ(ω),
Luego,
A(ω)
1 1
A(ω) = =√
|2 + jω| 4 + ω2
ω
y
φ(ω)
π/2
1
φ(ω) = Arg = − arctan(ω/2).
2 + jω ω
−π/2
f (t) = u(t).
1
F (ω) = + π δ(ω).
jω
Ası́,
A(ω)
1 1
A(ω) =
+ π δ(ω) = , para ω 6= 0
jω |ω|
ω
y
φ(ω)
1 π/2
φ(ω) = Arg + π δ(ω)
jω
1
= Arg , ω 6= 0
jω
( ω
π/2, ω < 0,
=
−π/2, ω > 0.
−π/2
♦
192 6.4. MAGNITUD Y FASE
f (t) = sgn(t).
2
F (ω) = .
jω
Ası́,
A(ω)
2 2
A(ω) = = , para ω 6= 0
jω |ω|
ω
y
φ(ω)
π/2
2
φ(ω) = Arg , ω 6= 0
jω
(
π/2, ω < 0, ω
=
−π/2, ω > 0.
−π/2
f (t) = δ(t).
F (ω) = 1.
Ası́,
A(ω)
1
A(ω) = 1, para ω ∈ R
ω
y
φ(ω)
ω
♦
CAPÍTULO 6. TRANSFORMADA DE FOURIER 193
ω
y
φ(ω)
2
φ(ω) = Arg = 0, para ω ∈ R.
1 + ω2
ω
♦
>> syms t w
>> f1 ( t ) = exp ( -2* t ) * h e a v i s i d e( t ) ; F1 ( w ) = f o u r i e r( f1 )
F1 ( w ) =
1/( w * i + 2)
>> f2 ( t ) = exp ( - abs ( t ) ) ; F2 ( w ) = f o u r i e r( f1 )
F2 ( w ) =
1/( w * i + 2)
>> f3 ( t ) = h e a v i s i d e( t +1) - h e a v i s i d e(t -1) ; F2 ( w ) = f o u r i e r( f3 )
F2 ( w ) =
- ( cos ( w ) * i - sin ( w ) ) / w + ( cos ( w ) * i + sin ( w ) ) / w
f u n c t i o n [A , phi ] = m a g n i t u d f a s e (f , w0 )
% magnitudfase Halla el valor de la m a g n i t u d A ( w ) y de la
% fase phi ( w ) , en w = w0
% p = A ( w0 ) , phi = phi ( w0 )
syms t w ;
F ( w ) = f o u r i e r( f ) ;
A = abs ( double ( F ( w0 ) ) ) ;
phi = angle ( double ( F ( w0 ) ) ) ;
end
>> syms t
>> w0 = 1/2;
>> f1 ( t ) = exp ( -2* t ) * h e a v i s i d e( t ) ; [ A1 , phi1 ] = m a g n i t u d f a s e ( f1 , w0 )
A1 =
0.4851
phi1 =
-0.2450
>> f2 ( t ) = exp ( - abs ( t ) ) ; [ A2 , phi2 ] = m a g n i t u d f a s e ( f2 , w0 )
A2 =
1.6000
phi2 =
0
>> f3 ( t ) = h e a v i s i d e( t +1) - h e a v i s i d e(t -1) ; [ A3 , phi3 ] = m a g n i t u d f a s e
( f3 , w0 )
A3 =
1.9177
phi3 =
0
>> syms t
>> w0 = 0 . 0 1 : 0 . 1 : 1 0 ;
>> f1 ( t ) = exp ( -2* t ) * h e a v i s i d e( t ) ; [ A1 , phi1 ]= m a g n i t u d f a s e ( f1 , w0 ) ;
>> f2 ( t ) = exp ( - abs ( t ) ) ; [ A2 , phi2 ]= m a g n i t u d f a s e( f2 , w0 ) ;
>> f3 ( t ) = h e a v i s i d e( t +1) - h e a v i s i d e(t -1) ; [ A3 , phi3 ] = m a g n i t u d f a s e
( f3 , w0 ) ;
>> syms t
>> w0 = 0 . 0 1 : 0 . 1 : 1 0 ;
>> f1 ( t ) = exp ( -2* t ) * h e a v i s i d e( t ) ; [ A1 , phi1 ]= m a g n i t u d f a s e ( f1 , w0 ) ;
CAPÍTULO 6. TRANSFORMADA DE FOURIER 195
A (ω) φ (ω)
1 1
0.8 2
0.6 1
0.4 0
0.2 −1
0 −2
−10 −5 0 5 10 −10 −5 0 5 10
A (ω) φ (ω)
2 2
2 1
1.5 0.5
1 0
0.5 −0.5
0 −1
−10 −5 0 5 10 −10 −5 0 5 10
A3(ω) φ3(ω)
2 4
1.5 2
1 0
0.5 −2
0 −4
−10 −5 0 5 10 −10 −5 0 5 10
MATLAB también cuenta con el comando ifourier(F) que permite calcular la trans-
formada inversa de Fourier de F (ω). Observe el siguiente ejemplo donde se calcula la
transformada inversa de Fourier de las siguientes transformadas:
>> syms w t
>> u ( w ) = h e a v i s i d e ( w ) ;
>> F1 = ( 1 + 1 0 *w ) *( u (10* w +1) -u (10* w ) ) +(1 -10*w ) *( u (10* w ) -u (10* w -1) ) ;
>> f1 = simple ( i f o u r i e r( F1 ) )
f1 =
(20* sin ( x /20) ^2) /( pi * x ^2)
>> F2 = dirac ( w ) +2*( u ( -w +1) -u ( -w -1) ) ;
>> f2 = simple ( i f o u r i e r( F2 ) )
f2 =
((2* sin ( x ) ) / x + 1/2) / pi
>> F3 = i *( u (2* w ) -u ( -2* w ) ) ;
>> f3 = simple ( i f o u r i e r( F3 ) )
f3 =
-1/( pi * x )
ii) Como p1 (t) = u(t + 1) − u(t − 1), entonces usando convenientemente las propiedades de
linealidad y desplazamiento en tiempo, obtenemos:
♦
CAPÍTULO 6. TRANSFORMADA DE FOURIER 197
♦
198 6.6. PROBLEMAS RESUELTOS
Solución. Primero calculemos por definición la transformada de Fourier de q1/2 (t). Apli-
cando la definición de la trasformada de Fourier y considerando la identidad trigonométrica
1 − cos(2α) = 2 sen 2 α, podemos escribir:
Z ∞
F q1/2 (t) = q1/2 (t)e−jωt dt
−∞
Z 0 Z 1/2
−jωt
= (2t + 1) e dt + (1 − 2t) e−jωt dt
−1/2 0
jω jω
2 − 2e 2 + jω 2e− 2 − 2 + jω
= −
ω2 ω2
1 jω jω
= 4 − 2e 2 − 2e− 2
ω2
4 8 sen 2 (ω/4)
= (1 − cos(ω/2)) = .
ω2 ω2
Ahora, usando la propiedad de desplazamiento en tiempo obtenemos:
8 e−jω sen 2 (ω/4)
F (ω) = F q1/2 (t − 1) = e−jω F q1/2 (t) = .
ω2
♦
Ahora, calculemos la transformada F e2t u(−t) . Sea g(t) la función definida como
g(t) = e−2t u(t).
Por definición, la transformada de Fourier de g(t) está dada por:
Z ∞ Z ∞ " #r
e−(2+jω)t 1
G(ω) = e−2t u(t)e−jωt dt = e−(2+jω)t dt = lı́m − = .
−∞ 0+ r→∞ 2 + jω 2 + jω
0
Problema 6.7. Sea f (t) una función de dominio continuo cuya transformada de Fourier
está dada por:
F (ω) = p1 ((ω − 1)/2).
Determine la transformada de Fourier de g(t) = f (2t − 1) e−2jt .
= e−jω/2 F {f (2t)}
−jω/2 1
= e F (ω/2)
2
1 −jω/2
= e p1 ((ω − 2)/4).
2
De esta forma, la transformada de Fourier de g(t) es
1 −j(ω+2)/2
G(ω) = e p1 (ω/4)
2
♦
200 6.6. PROBLEMAS RESUELTOS
Problema 6.12. Determine los espectros de magnitud y fase del pulso rectangular p1 (t).
A(ω)
2
ω
−6π −5π −4π −3π −2π −π π 2π 3π 4π 5π 6π
202 6.6. PROBLEMAS RESUELTOS
φ(ω)
π
···
ω
−6π −5π −4π −3π −2π −π π 2π 3π 4π 5π 6π
···
−π
Problema 6.13. Determine los espectros de magnitud y fase del pulso triangular q1 (t).
2(1 − cos ω)
F (ω) = .
ω2
Luego, el espectro de magnitud de q1 (t) es:
2 (1 − cos ω)
A(ω) = , para ω 6= 0.
ω2
A(ω)
1
ω
−6π −5π −4π −3π −2π −π π 2π 3π 4π 5π 6π
Como 2(1 − cos ω) > 0 y ω 2 > 0, para todo ω ∈ R, entonces el espectro de fase de q1 (t) es:
2(1 − cos ω)
φ(ω) = Arg (F (ω)) = Arg = 0, para ω = 6 0.
ω2
♦
CAPÍTULO 6. TRANSFORMADA DE FOURIER 203
Problema 6.14. Determine los espectros de magnitud y fase de f (t) = u(t) − u(t − 5).
Solución. La transformada de Fourier de f (t) = u(t) − u(t − 5) (ver Problema 6.3) es:
1 − e−5jω 1 − cos(5ω) + j sen (5ω)
F (ω) = = .
jω jω
Luego, el espectro de magnitud de f (t) es:
p p
(1 − cos(5ω))2 + sen 2 (5ω) 2 − 2 cos(5ω)
A(ω) = = , para ω 6= 0.
|ω| |ω|
A(ω)
5
ω
−6π −5π −4π −3π −2π −π π 2π 3π 4π 5π 6π
···
ω
···
−π
♦
Problema 6.15. Determine los espectros de magnitud y fase de f (t) = q1/2 (t − 1).
Solución. La transformada de Fourier de f (t) = q1/2 (t − 1) (ver Problema 6.4) es:
8 e−jω sen 2 (ω/4)
F (ω) = .
ω2
Luego, el espectro de magnitud de f (t) es:
8 sen 2 (ω/4)
A(ω) = , para ω 6= 0.
ω2
204 6.6. PROBLEMAS RESUELTOS
A(ω)
1/2
ω
−3π −2π −π π 2π 3π
8 sen 2 (ω/4)
Como ω2
> 0, para ω 6= 0, podemos escribir:
−jω 8 sen 2 (ω/4)
arg (F (ω)) = arg (e ) + Arg = −ω, para ω 6= 0.
ω2
ası́, el espectro de fase de f (t) es, para ω 6= 0:
∞
X ∞
X
φ(ω) = (2kπ − ω)pπ (ω − 2kπ) − (ω + 2kπ)pπ (ω + 2kπ).
k=0 k=1
φ(ω)
π
··· ···
ω
−5π −4π −3π −2π −π π 2π 3π 4π 5π
−π
Problema 6.16. Determine los espectros de magnitud y fase de f (t) = e2t u(−t) + u(t).
Solución. La transformada de Fourier de f (t) = e2t u(−t) + u(t) (ver Problema 6.6) es:
2
F (ω) = + πδ(ω).
jω(2 − jω)
Luego, el espectro de magnitud de f (t) es:
2
A(ω) = √ , para ω 6= 0.
|ω| 4 + ω 2
A(ω)
ω
−3π −2π −π π 2π 3π
CAPÍTULO 6. TRANSFORMADA DE FOURIER 205
φ(ω)
π
2
− π2
6.2. Sea f (t) una función de dominio continuo cuya transformada de Fourier es F (ω).
Pruebe que la transformada de Fourier de la función g(t) = f (t) cos ω0 t, con
ω0 ∈ R, está dada por
1
F{g(t)} = {f (t) cos ω0 t} = [F (ω − ω0 ) + F (ω + ω0 )].
2
6.3. Sea f (t) una función de dominio continuo cuya transformada de Fourier es
F (ω) = p1 ((ω − 1)/2).
Determine la transformada de Fourier de las siguientes funciones de dominio con-
tinuo, utilizando las propiedades de la transformada de Fourier:
206 6.7. PROBLEMAS PROPUESTOS
6.6. Calcular los espectros de magnitud y fase de las siguientes funciones de dominio
continuo:
7.1 Definición
Sea f (t) una función de dominio continuo definida de R en C. La transformada de
Laplace bilateral de f (t) se define como
Z ∞
F (s) = L[f (t)] = f (t) e−st dt,
−∞
Esta ecuación nos indica que las funciones componentes de F (s) = u(x, y) + jv(x, y), son:
Z ∞ Z ∞
−xt
u(x, y) = f (t) e cos(yt) dt y v(x, y) = − f (t) e−xt sen (yt) dt.
−∞ −∞
207
208 7.1. DEFINICIÓN
Por lo tanto, los valores de σ = Re s que satisfacen la ecuación (7.2), determinan explı́cita-
mente la región de convergencia de F (s). Note que, por la Proposición 7.1, la transformada
de Laplace es analı́tica en su región de convergencia.
L
f (t) ←→ F (s), s∈D
CAPÍTULO 7. TRANSFORMADA DE LAPLACE 209
que se lee par de transformadas, denota que F (s) es la transformada de Laplace de f (t)
con región de convergencia el conjunto D ⊂ C, y que f (t) es la transformada inversa de
Laplace de F (s).
L
af (t) + bg(t) ←→ aF (s) + bG(s), s ∈ Df ∩ Dg ,
Ahora, la región de convergencia de L{af (t) + bg(t)} son todos los s = σ + jω tales que
Z ∞
|af (t) + bg(t)| e−σt dt < ∞,
−∞
pero
Z ∞ Z ∞ Z ∞
|af (t) + bg(t)| e−σt dt, ≤ |a| |f (t)| e−σt dt, +|b| |g(t)| e−σt dt. (7.4)
−∞ −∞ −∞
Ası́, por (7.4), (7.5) y (7.6), la región de convergencia de L{af (t) + bg(t)} son todos los
s = σ + jω tales que s ∈ Df ∩ Dg .
L
f (t − t0 ) ←→ e−t0 s F (s), s ∈ D,
para todo t0 ∈ R.
210 7.2. PROPIEDADES DE LA TRANSFORMADA DE LAPLACE
L
es0 t f (t) ←→ F (s − s0 ), (s − s0 ) ∈ D,
para todo s0 ∈ C.
Ahora, la región de convergencia de L{es0 t f (t)} son todos los s = σ + jω tales que
Z ∞
st
e 0 f (t) e−σt dt < ∞.
−∞
de lo cual se deduce que la región de convergencia de L{es0 t f (t)} son todos los s = σ + jω
tales que (s − s0 ) ∈ D.
CAPÍTULO 7. TRANSFORMADA DE LAPLACE 211
L 1
f (at) ←→ F (s/a), (s/a) ∈ D,
|a|
para todo a ∈ R.
Demostración. Por la definición de la transformada de Laplace,
Z ∞
L{f (at)} = f (at) e−st dt. (7.7)
−∞
Ahora, la región de convergencia de L{f (at)} son todos los s = σ + jω tales que
Z ∞
|f (at)| e−σt dt < ∞.
−∞
de lo cual se deduce que la región de convergencia de L{f (at)} son todos los s = σ + jω
tales que (s/a) ∈ D.
L
f (t) ←→ F (s), s ∈ D.
L
f ∗ g (t) ←→ F (s) G(s), s ∈ Df ∩ Dg .
Ahora, la región de convergencia de L{f ∗ g (t)} son todos los s = σ + jω tales que
Z ∞ Z ∞ Z ∞
−σt
−σt
|f ∗ g (t)| e dt =
f (τ )g(t − τ ) dτ e
dt < ∞.
−∞ −∞ −∞
en otras palabras, la región de convergencia de L{f ∗ g (t)} son todos los s = σ + jω tales
que s ∈ Df ∩ Dg .
Demostración. Asumamos que f (t) es de orden exponencial, esto es, existen constantes
K > 0 y a ∈ R tales que
Además, si f (t) es de orden exponencial, entonces (se deja como ejercicio para el lector la
prueba de esta ecuación):
α
lı́m f (t)e−st −α = 0, para todo s ∈ D. (7.8)
α→∞
CAPÍTULO 7. TRANSFORMADA DE LAPLACE 213
Se tiene Z ∞
d d
L f (t) = f (t) e−st dt.
dt −∞ dt
Ası́, usando la ecuación (7.8) y la integración por partes con u = e−st y dv = f ′ (t), se
obtiene Z ∞
d
−st α
L f (t) = lı́m f (t)e −α
+s f (t) e−st dt = sF (s).
dt α→∞ −∞
L d
t f (t) ←→ − F (s), s ∈ D.
ds
d
de donde se deduce que L {t f (t)} = − F (s). Además, como F (s) es analı́tica en todo
ds
d
s ∈ D, entonces F (s) también es analı́tica en todo s ∈ D, en otras palabras, D es la
ds
d
región de convergencia de F (s).
ds
para toda función f (t) de dominio continuo. De esta forma, usando la propiedad de con-
volución en la ecuación (7.9)
y como
lı́m e−st = 0,
s→∞
la expresión de arriba se reduce a
o, equivalentemente,
f (0+ ) = lı́m s F (s).
s→∞
CAPÍTULO 7. TRANSFORMADA DE LAPLACE 215
y como
lı́m e−st = 1,
s→0
o, equivalentemente,
lı́m f (t) = lı́m s F (s).
t→∞ s→∞
Sean f (t) y g(t) dos funciones de dominio continuo con transformada de Laplace F (s),
para s ∈ Df , y G(s), para s ∈ Dg , respectivamente.
1
Escalamiento en tiempo L {f (at)} = F (s/a), s ∈ {s ∈ C : (s/a) ∈ Df }
|a|
n o
Conjugación L f (t) = F (s), s ∈ Df
1.
L
δ(t) ←→ 1, s∈C
Se tiene que
Z ∞ Z ∞
−st
L{δ(t)} = δ(t) e dt = δ(t) e−σt e−jωt dt
Z−∞
∞ Z ∞ −∞
−σt
= δ(t) e cos(ωt) dt − j δ(t) e−σt sen (ωt) dt
−∞ −∞
= e−σ·0 cos(ω · 0) − je−σ·0 sen (ω · 0) = 1.
pero Z ∞
δ(t)(t) e−σt dt = 1;
−∞
2.
L 1
u(t) ←→ , Re s > 0
s
Se tiene que
Z ∞ Z ∞ ∞
−st −st e−st 1
L{u(t)} = u(t) e dt = e dt = − = .
−∞ 0+ s
0+ s
Ahora, como
Z ∞ Z ∞
|u(t)| e−σt dt = e−σt dt < ∞ ⇔ σ > 0,
−∞ 0+
3.
L 1
−u(−t) ←→ , Re s < 0
s
Hallemos la transformada de Laplace L{−u(−t)} mediante dos procedimientos:
(a) por definición, y (b) usando las propiedades de la transformada de Laplace.
(a) Se tiene que
Z ∞ Z 0− 0−
−st −st e−st 1
L{−u(−t)} = [−u(−t)] e dt = − e dt = = .
−∞ ∞ s
∞ s
218 7.3. ALGUNOS PARES DE TRANSFORMADAS
Como Z Z
∞ 0−
−σt
| − u(−t)| e dt = e−σt dt < ∞ ⇔ σ < 0,
−∞ ∞
entonces la región de convergencia de L{−u(−t)} son todos los s = σ + jω ∈ C tales
que Re s = σ < 0.
(b) Usando las propiedades de linealidad y escalamiento en tiempo, se obtiene
1 1
L{−u(−t)} = −L{u(−t)} = − − L{u(t)}(−s) = , Re (−s) > 0,
| − 1| s
o, equivalentemente,
1
L{−u(−t)} = , Re s < 0.
s
Observación 7.3. Note que las transformadas de Laplace de u(t) y −u(−t) son
algebraicamente iguales, pero tienen regiones de convergencia diferentes, más aún,
son complementarias.
4.
L 1
e−αt u(t) ←→ , Re s > −Re α
s+α
Usando la propiedad de desplazamiento en el dominio de s, se obtiene
1
L{e−αt u(t)} = L{u(t)}(s − (−α)) = , Re (s + α) > 0,
s+α
o, equivalentemente,
1
L{e−αt u(t)} = , Re s > −α.
s+α
La transformada de Laplace de e−αt u(t) también se puede hallar por definición, lo
cual se deja como ejercicio para el lector.
5.
L 1
t u(t) ←→ , Re s > 0
s2
Se tiene que Z Z
∞ ∞
−st
L{t u(t)} = [t u(t)] e dt = te−st dt.
−∞ 0+
Usando integración por partes, se obtiene
∞ Z
te−st 1 ∞ −st
L{t u(t)} = − + e dt
s 0+ s 0+
−st ∞
1 e = 1.
= 0+ −
s s 0+ s2
Como
Z ∞ Z ∞
−σt
|t u(t)| e dt = t e−σt dt
−∞ 0+
Z ∞
1
= e−σt dt < ∞ ⇔ σ > 0,
σ 0+
CAPÍTULO 7. TRANSFORMADA DE LAPLACE 219
6.
L 1
sen t u(t) ←→ , Re s > 0
s2 +1
Se tiene que
Z ∞
L{sen t u(t)} = [sen t u(t)] e−st dt
−∞
Z ∞ jt
e − e−jt
= e−st dt
0+ 2j
Z ∞ Z ∞
1 (j−s)t −(j+s)t
= e dt − e dt
2j 0+ +
"" #∞ "0 #∞ #
1 e(j−s)t e−(j+s)t
= − −
2j j−s + j+s +
0
0
1 1 1 1 2j 1
= − = = 2 .
2j s − j s+j 2j s2 + 1 s +1
Como Z ∞ Z ∞
−σt
|sen t u(t)| e dt ≤ e−σt dt < ∞ ⇔ σ > 0,
−∞ 0+
7.
L s
cos t u(t) ←→ , Re s > 0
s2 +1
Usando las propiedades de linealidad y desplazamiento en el dominio de s, se obtiene
jt
e + e−jt
L{cos t u(t)} = L u(t)
2
1
= L{ejt u(t)} + L{e−jt u(t)}
2
1
= [L{u(t)}(s + j) + L{u(t)}(s − j)]
2
1 1 1
= +
2 s+j s−j
s
= , {s ∈ C : Re (s + j) > 0} ∩ {s ∈ C : Re (s − j) > 0}
s2 + 1
o, equivalentemente,
s
L{cos t u(t)} = , Re s > 0.
s2 +1
220 7.3. ALGUNOS PARES DE TRANSFORMADAS
8.
L 2jT senh(T s)
pT (t) ←→ , Re s 6= 0
Ts
Se tiene que
Z ∞
L{pT (t)} = pT (t) e−st dt
−∞
Z T−
= e−st dt
−T +
T −
e−st
= −
s −T +
esT − e−sT
=
s
2jT sen (sT )
= , s 6= 0.
sT
Como
Z ∞ Z T−
−σt eσT − e−σT
|pT (t)| e dt = e−σt dt = <∞ ⇔ σ 6= 0,
−∞ −T + σ
2jT senh(T s)
pT (t) , Re s 6= 0
Ts
222 7.4. CÁLCULO DE LA TRANSFORMADA INVERSA DE LAPLACE
Proposición 7.2. Sea F (s) la transformada de Laplace de f (t), con región de convergencia
dada por Re s > −a, donde a > 0. Si |F (s)| < b|s|−m para b > 0, un entero m > 0, y
R0 > 0 tal que |s| > R0 , entonces f (t) = 0, para t < 0.
Demostración. Se deja como ejercicio para el lector. Ayuda: utilice un contorno cerrado
simple C apropiado y el hecho que |F (s)| < 1/|s|.
1
F (s) = , Re s > −2.
s+2
z(t) = j t, −R ≤ t ≤ R,
ω Re s > −2
C1
C2
−2
σ
Se tiene que
Z Z Z
1 st 1 st 1
e F (s) ds = e F (s) ds + est F (s) ds
2πj C 2πj C1 2πj C2
Z Z jR st
1 est 1 e
= ds + ds. (7.11)
2πj C1 s + 2 2πj −jR s + 2
Se deja como ejercicio para el lector verificar que
Z
est
lı́m ds = 0, para t ≥ 0.
R→∞ C1 s + 2
t
♦
b
pI1 r1 < Re s < r2
b
pD
1
r1 r2 σ
b
pI3
b
pD
2
b
pI2
5s − 1
F (s) = , −1 < Re s < 2.
s3 − 3s − 2
Ahora,
est (5s − 1)
Res est F (s) = Res
s=−1 s=−1 s3 − 3s − 2
d est (5s − 1)
= = (2t − 1) e−t (7.13)
dz (s − 2) s=−1
y
st
st e (5s − 1)
Res e F (s) = Res 3
s=2 s=2 s − 3s − 2
st
e (5s − 1)
= = e2t . (7.14)
(s + 1)2 s=2
♦
CAPÍTULO 7. TRANSFORMADA DE LAPLACE 225
donde F1 (s), F2 (s), . . . , FR (s) son funciones tales que se les conoce su transformada inversa
de Laplace f1 (t), f2 (t), . . . , fR (t). Entonces, la transformada inversa de Laplace de F (s)
está dada por
f (t) = f1 (t) + f2 (t) + · · · + fR (t).
donde
A1 = Res [F (s)] = 1,
s=−1
A2,2 = (s − 2)2 F (s) s=2 = 336,
d 2
A2,1 = [(s − 2) F (s)] = 800,
ds s=2
A3,2 = (s − 3)2 F (s) s=3 = 468,
d 2
A3,1 = [(s − 3) F (s)] = −801.
ds s=3
Ası́,
1 1 1
f (t) = A1 L−1 + A2,1 L−1 + A2,2 L−1
s+1 s−2 (s − 2)2
1 1
+A3,1 L−1 + A3,2 L−1 .
s−3 (s − 3)2
226 7.5. TRANSFORMADA DE LAPLACE CON MATLAB
>> syms t
>> f1 ( t ) = exp ( -2*t ) * h e a v i s i d e( t ) ;
>> F1 = l a p l a c e( f1 )
F1 =
1/( s + 2)
>> f2 ( t ) = (2* t -1) * exp ( - t ) * h e a v i s i d e ( t ) - exp (2* t ) * h e a v i s i d e( - t ) ;
>> F2 = l a p l a c e( f2 )
F2 =
2/( s + 1) ^2 - 1/( s + 1)
f u n c t i o n [ F ] = T r a n s L a p l a c e ( fn , fc )
% T r a n s L a p l a c e Halla la t r a n s f o r m a d a b i l a t e r a l de L a p l a c e
% de una f u n c i ó n f ( t ) = fn ( t ) + fc ( t )
% F : T r a n s f o r m a d a b i l a t e r a l de L a p l a c e
syms t s ;
Fc ( s ) = l a p l a c e( fc ) ;
Fn ( s ) = l a p l a c e( fn ( -t ) ) ;
F ( s ) = Fn ( - s ) + Fc ( s ) ;
end
donde fn (t) es una función no causal y fc (t) es una función causal. Recuerde que una
función g(t) se dice causal si g(t) = 0 para t < 0, y ella es no causal, si g(t) = 0 para t > 0.
Calculemos la transformada bilateral de Laplace de las funciones f1 (t) y f2 (t) dadas
arriba. Observe que f1 (t) es una función causal y f2 (t) = f2n (t) + f2c (t), donde
>> syms t
>> f1c ( t ) = exp ( -2* t ) * h e a v i s i d e ( t ) ; f1n ( t ) = sym (0) ;
>> F1 = T r a n s L a p l a c e ( f1n , f1c )
F1 ( s ) =
1/( s + 2)
>> f2c ( t ) =(2* t -1) * exp ( - t ) * h e a v i s i d e( t ) ;
>> f2n ( t ) = - exp (2* t ) * h e a v i s i d e ( -t ) ;
>> F2 = T r a n s L a p l a c e ( f2n , f2c )
F2 ( s ) =
2/( s + 1) ^2 - 1/( s + 1) + 1/( s - 2)
Por otra parte, MATLAB también cuenta con el comando ilaplace que halla la trans-
formada inversa de Laplace. Observe el siguiente ejemplo donde se halla la transformada
inversa de Laplace de las transformadas F1 (s) y F2 (s) dadas arriba.
>> syms s
>> F1 ( s ) = 1/( s + 2) ; f1 ( t ) = i l a p l a c e( F1 )
f1 ( t ) =
exp ( -2* t )
>> F2 ( s ) = 2/( s + 1) ^2 - 1/( s + 1) + 1/( s - 2) ;
>> f2 ( t ) = i l a p l a c e( F2 )
f2 ( t ) =
exp (2* t ) - exp ( - t ) + 2* t * exp ( - t )
228 7.6. PROBLEMAS RESUELTOS
Problema 7.3. Determinar la transformada de Laplace de f (t) = e−3t (u(t) − u(t − 4)).
donde
f1 (t) = e−3t u(t) y f2 (t) = e−3(t+4) u(t).
Ası́, aplicando la propiedad de linealidad tenemos:
donde D1 , D2 ⊂ C son las regiones de convergencia de L {f1 (t)} y L {f2 (t − 4)}, respecti-
vamente. Calculemos las transformadas L {f1 (t)} y L {f2 (t − 4)}.
De esta forma, sustituyendo convenientemente las expresiones de L {f1 (t)} y L {f2 (t − 4)}
en (7.16), se tiene que la transformada de Laplace de f (t) es:
1 e−4(s+3) 1 − e−4(s+3)
F (s) = − = , Re s > −3.
s+3 s+3 s+3
Seguidamente se muestra el diagrama de polos (rojo) y ceros (azul) de F (s).
−3 σ
b
.. b
♦
230 7.6. PROBLEMAS RESUELTOS
1 n (4j+2)t o 1 n o
= L e u(t) − L e−(4j−2)t u(t)
2j 2j
1 1
= L {u(t)} (s − (4j + 2)) − L {u(t)} (s + (4j − 2))
2j 2j
1 1
= − .
2j(s − (4j + 2)) 2j(s + (4j − 2))
| {z } | {z }
Re s>2 Re s>2
4
F (s) = , Re s > 2.
(s − (4j + 2))(s + (4j − 2))
2 σ
2 − 4j b
♦
CAPÍTULO 7. TRANSFORMADA DE LAPLACE 231
f (t)
−1 1 t
−1
es + e−s − 2
F (s) = − (es − e−s )
s
2 cosh s − 2
= − 2 senh s, Re s > 0.
s
♦
Problema 7.6. Sea f (t) una función de dominio continuo tal que f (t) = 0 para t < 0,
cuya transformada de Laplace es:
s e−2s
F (s) = , Re s > 2.
(s + 1)(s − 2)
donde
3s4 − 3s3 + s2 − 2s − 1
F1 (s) = , 0 < Re s < 1.
s2 (s − 1)3
CAPÍTULO 7. TRANSFORMADA DE LAPLACE 233
Problema 7.9. Sea f (t) una función de dominio continuo cuya transformada de Laplace
F (s), −1 < Re s < 0, satisface:
• s1 = 2 y s2 = 3, son los únicos ceros de F (s);
• p1 = −1, p2 = 0 y p3 = 1, son los únicos puntos singulares de F (s), además, p1 y p3
son polos simples y p2 es un polo de orden 2;
1
• F (4) = .
24
Determinar la transformada inversa de F (s).
Solución. Por las condiciones dadas en el enunciado del problema, una expresión ma-
temática de F (s) puede ser (existen otras considerando ceros con multiplicidad mayor que
1):
k(s − 2)(s − 3)
F (s) = ,
(s + 1)s2 (s − 1)
donde k ∈ C. Ahora bien, como F (4) = 1/24, entonces podemos escribir:
k(4 − 2)(4 − 3) 1
= ⇔ k = 5.
(4 + 1)42 (4 − 1) 24
Ası́,
5(s − 2)(s − 3)
F (s) = , −1 < Re s < 0.
(s + 1)s2 (s − 1)
Se tiene que el polo a la izquierda de la región de convergencia es: p1 = −1 y los polos a
la derecha de la región de convergencia son: p2 = 0 y p3 = 1. Por la Proposición 7.3, la
transformada inversa de Laplace de F (s) está dada por:
f (t) = Res est F (s) u(t) − Res est F (s) + Res est F (s) u(−t).
s=−1 s=0 s=1
Ahora,
st 5est (s − 2)(s − 3)
Res e F (s) = = −30 e−t ,
s=−1 s2 (s − 1)
s=−1
st
d 5e (s − 2)(s − 3)
Res est F (s) = = 25 − 30t,
s=0 ds (s + 1)(s − 1) s=0
st
5e (s − 2)(s − 3)
Res est F (s) = = 5 et .
s=1 (s + 1)s2
s=1
7.3. La transformada de Laplace de una función causal f (t) (f (t) = 0 para t < 0) es:
s e−2s
F (s) = , Re s > −1.
s2 + 2s + 1
Determinar la transformada de Laplace de las siguientes funciones:
t d
a) g(t) = 3f e) g(t) = t f (t)
3 dt
f) g(t) = δ(t − a) f (t), a ∈ R
b) g(t) = t f (t)
d
c) g(t) = t f (t − 1) g) g(t) = (t − 1) f (t − 1) + f (t)
dt
Z t
d h) g(t) = f (τ ) τ
d) g(t) = f (t)
dt 0
7.4. Determinar los valores inicial y final de las funciones causales f (t), cuya transfor-
mada de Laplace se indica a continuación, sin calcular la transformada inversa de
Laplace. Si no existe valor final, indicar el motivo.
1
a) F (s) = , Re s > −2
s+2
236 7.7. PROBLEMAS PROPUESTOS
1
b) F (s) = , Re s > −2
(s + 2)25
6
c) F (s) = , Re s > 0
s(s2 + 25)
s+2
d) F (s) = , Re s > −3
s+3
s2 + s + 3
e) F (s) = , Re s > 2
s4 − 9s2 + 4s + 12
s
f) F (s) = 2 , Re s > 3
s − 2s − 3
7.5. Determinar la transformada inversa de Laplace de cada una de las siguientes
transformadas, mediante el procedimiento por residuos.
1
a) F (s) = , Re s < −1
s2 +1
1
b) F (s) = , Re s > −a, con a < 0
s2 + a2
s−2
c) F (s) = , −1 < Re s < 3
(s + 1)(s − 3)
s2 + 2s − 2
d) F (s) = , Re s > 1
s(s − 1)
s3 + 4s2 + s − 1
e) F (s) = , −1 < Re s < 0
s3 (s − 1)(s + 1)2
2s3 + 3s2 + 6s + 4
f) F (s) = , −1 < Re s < 0
(s2 + 4)(s2 + 2s + 2)
7.6. Sea f (t) una función de dominio continuo cuya transformada de Laplace F (s)
satisface:
8.1 Definición
Sea f (n) una función de dominio discreto definida de Z en C. La transformada z
bilateral de f (n) se define como
∞
X
F (z) = Z[f (n)] = f (n)z −n ,
n=−∞
Observación 8.1. De la definición de la transformada z se infiere que f (n) son los coefi-
cientes del desarrollo
P en serie de−nLaurent de F (z) centrado en z0 = 0. Además, F (z) es la
suma de la serie ∞ n=−∞ f (n)z .
En la siguiente proposición se establece que si la transformada z existe, entonces ella
es una función analı́tica.
F (z) para cada z tal que |F (z)| < ∞. Por lo tanto, F (z) es una función analı́tica para z
tal que |F (z)| < ∞.
237
238 8.1. DEFINICIÓN
|a|
♦
En el resto del capı́tulo, llamaremos transformada z de una función a su transformada
z bilateral.
CAPÍTULO 8. TRANSFORMADA Z 239
que se lee par de transformadas, denota que F (z) es la transformada z de f (n) con región
de convergencia el conjunto D ⊂ C, y que f (n) es la transformada z inversa de F (z).
Z
f (−n) ←→ F (z −1 ), z −1 ∈ D.
Z
z0n f (n) ←→ F (z/z0 ) , z ∈ D,
Z
ejω0 n f (n) ←→ F e−jω0 z , z ∈ D,
Z
an f (n) ←→ F a−1 z , z ∈ D,
para z0 ∈ C, ω0 ∈ R y a ∈ R.
CAPÍTULO 8. TRANSFORMADA Z 241
Demostración. Como ejω0 n f (n) y an f (n) son casos particulares de z0n f (n) (tome respec-
tivamente z0 = ejω0 y z0 = a), entonces solo demostraremos que Z{z0n f (n)} = F (z/z0 ),
z ∈ D. Se tiene que
∞
X ∞
X
Z{z0n f (n)} = z0n f (n)z −n = f (n) (z/z0 )−n = F (z/z0 ), z ∈ D.
n=−∞ n=−∞
Z
f (n) ←→ F (z), z ∈ D.
Z
f ∗ g (n) ←→ F (z) G(z), z ∈ Df ∩ Dg .
Z{f (n) ∗ u(n)} = Z{f (n)} Z{u(n)}, z ∈ D ∩ {z ∈ C : |z| > 1}. (8.3)
Z d
n f (n) ←→ −z F (z), z ∈ D.
dz
Demostración. Se tiene que, para todo z ∈ D,
" ∞ #
d d X
−n
−z F (z) = −z f (n)z
dz dz n=−∞
∞
X d −n
= −z f (n) z
n=−∞
dz
∞
X ∞
X
−(n+1)
= z (n f (n)) z = (n f (n))z −n ,
n=−∞ n=−∞
d
de donde se deduce que Z{n f (n)} = −z F (z), z ∈ D.
dz
Teorema 8.10 (Teorema del Valor Inicial).
Sea f (n) una función tal que f (n) = 0 para n < 0. Si la transformada z de f (n) es F (z)
para todo z ∈ C tal que |z| > a con a > 0, entonces
De esta forma, tomando lı́mite en ambos lados de (8.4), cuando z → 1, entonces por (8.5),
se tiene:
N
z−1 X
lı́m F (z) = lı́m [f (n) − f (n − 1)] = lı́m f (N ).
z→1 z N →∞
n=0
N →∞
Sean f (n), f1 (n) y f2 (n), funciones de dominio discreto con transformadas z dadas por:
F (z) para z ∈ D, F1 (z) para z ∈ D1 y F2 (z) para z ∈ D2 , respectivamente.
(
D, si n0 ≤ 0,
Desplazamiento en tiempo Z {f (n − n0 )} = z −n0 F (z), z∈
D − {0}, si n0 > 0
n o
Conjugación Z f (n) = F (z), z∈D
z−1
Primera diferencia Z {f (n) − f (n − 1)} = z F (z), z ∈ D ∩ {|z| > 0}
( n
)
X
z
Acumulación Z f (k) = z−1 F (z), z ∈ D ∩ {|z| > 1}
k=−∞
d
Diferenciación en z Z {n f (n)} = −z F (z), z∈D
dz
Teorema del valor inicial f (0) = lı́m F (z)
z→∞
z−1
Teorema del valor final lı́m f (n) = lı́m F (z)
n→∞ z→1 z
CAPÍTULO 8. TRANSFORMADA Z 245
2.
Z z
u(n) ←→ , |z| > 1
z−1
Se tiene que
∞
X ∞
X
Z{u(n)} = u(n)z −n = (z −1 )n
n=−∞ n=0
1 z
= −1
= , |z| > 1.
1−z z−1
La región de convergencia de Z{u(n)} se muestra en la Figura 8.1. Aquı́ también
se aprecia la ubicación de los polos y ceros de Z{u(n)}. Esta gráfica se denomina
diagrama de polos y ceros.
y
|z| > 1
1
b b
0 1 x
3.
Z z
−u(−n − 1) ←→ , |z| < 1
z−1
Se tiene que
∞
X −1
X
Z{−u(−n − 1)} = [−u(−n − 1)]z −n = − z −n .
n=−∞ n=−∞
Haciendo el cambio de variable k = −n, se obtiene
1
" ∞
#
X X
k k
Z{−u(−n − 1)} = − z =− z −1
k=∞ n=0
1
= − −1 , |z| < 1
1−z
246 8.3. ALGUNOS PARES DE TRANSFORMADAS
o, equivalentemente,
z
Z{−u(−n − 1)} = , |z| < 1.
z−1
En la Figura 8.2 se muestra el diagrama de polos y ceros de Z{−u(−n − 1)}.
y
|z| < 1
1
b b
0 1 x
4.
Z z
αn u(n) ←→ , |z| > |α|
z−α
Se tiene que
∞
X ∞
X
n n −n
Z{α u(n)} = [α u(n)]z = (αz −1 )n
n=−∞ n=0
1
= , |αz −1 | < 1
1 − αz −1
o, equivalentemente,
z
Z{αn u(n)} = , |z| > |α|.
z−α
5.
Z z
−αn u(−n − 1) ←→ , |z| < |α|
z−α
Se tiene que
∞
X −1
X
n n −n
Z{−α u(−n − 1)} = [−α u(−n − 1)]z =− (α−1 z)−n .
n=−∞ n=−∞
6.
Z αz
n αn u(n) ←→ , |z| > |α|
(z − α)2
Se tiene que
z
Z{αn u(n)} = , |z| > |α|.
z−α
Luego, aplicando la propiedad de diferenciación en z, se obtiene
d
Z{n αn u(n)} = −z [Z{αn u(n)}]
dz
d z
= −z , |z| > |α|
dz z − α
o, equivalentemente,
αz
Z{n αn u(n)} = , |z| > |α|.
(z − α)2
7.
Z αz
−n αn u(−n − 1) ←→ , |z| < |α|
(z − α)2
Se tiene que
z
Z{u(−n − 1)} = − , |z| < 1. (8.6)
z−1
Aplicando las propiedades de linealidad y diferenciación en z, se tiene
d
Z{−n αn u(−n − 1)} = z Z{αn u(−n − 1)} (8.7)
dz
o, equivalentemente,
αz
Z{−n αn u(−n − 1)} = , |z| < |α|.
(z − α)2
248 8.3. ALGUNOS PARES DE TRANSFORMADAS
8.
Z z
e−aT n u(n) ←→ , |z| > e−aT
z − e−aT
Se tiene que
∞
X
Z{e−aT n u(n)} = e−aT n u(n)z −n
n=−∞
∞
X n 1 −aT −1
= e−aT z −1 = , e z <1
n=0
1 − e−aT z −1
o, equivalentemente,
z
Z{e−aT n } = , |z| > e−aT .
z − e−aT
9.
Z z(z − cos(ω))
cos(ωn)u(n) ←→ , |z| > 1
z2 − 2z cos(ω) + 1
Se tiene que
∞
X
Z{cos(ωn)u(n)} = cos(ωn)u(n)z −n
n=−∞
∞ jωn
X e + e−jωn
= z −n
2
n=0
"∞ ∞
#
1 X jω −1 n X −jω −1 n
= e z + e z
2 n=0 n=0
1 1 1
= + ,
2 1 − ejω z −1 1 − e−jω z −1
z(z − cos(ω))
Z{cos(ωn)u(n)} = , |z| > 1.
z2 − 2z cos(ω) + 1
Función Transformada z
δ(n) 1, z∈C
z
u(n) , |z| > 1
z−1
z
nu(n) , |z| > 1
(z − 1)2
z(z + 1)
n2 u(n) , |z| > 1
(z − 1)3
z
−u(−n − 1) , |z| < 1
z−1
z
αn u(n) , |z| > |α|
z−α
z
−αn u(−n − 1) , |z| < |α|
z−α
αz
n αn u(n) , |z| > |α|
(z − α)2
αz
−n αn u(−n − 1) , |z| < |α|
(z − α)2
z
e−aT n u(n) −aT
, |z| > e−aT
z−e
z(z − cos(ω))
cos(ωn)u(n) , |z| > 1
z 2 − 2z cos(ω) + 1
z sen (ω)
sen (ωn)u(n) , |z| > 1
z2 − 2z cos(ω) + 1
z(z − α cos(ω))
αn cos(ωn)u(n) , |z| > |α|
z 2 − 2zα cos(ω) + α2
zα sen (ω)
αn sen (ωn)u(n) , |z| > 1
z2 − 2zα cos(ω) + α2
250 8.4. CÁLCULO DE LA TRANSFORMADA Z INVERSA
z
F (z) = , |z| > |a|.
z−a
Solución. Sea C la circunferencia |z| = r > |a|. Ahora, considerando la expresión de F (z)
se tiene
Z
1
f (n) = F (z) z n−1 dz
2πj C
Z
1 zn
= dz, para cada n ∈ Z. (8.9)
2πj C z − a
f (n) = an u(n).
♦
CAPÍTULO 8. TRANSFORMADA Z 251
de donde se deduce
∞
z X
= nu(n) z −n , |z| > 1. (8.13)
(z − 1)2 n=−∞
Sustituyendo convenientemente (8.11), (8.12) y (8.13) en (8.10), se obtiene:
∞ ∞ ∞
1 X 3 X 1 X
F (z) = (−1)n u(n) z −n + u(n) z −n + nu(n) z −n
4 n=−∞ 4 n=−∞ 2 n=−∞
∞
X 1 3 1
= (−1) u(n) + u(n) + nu(n) z −n , |z| > 1.
n
n=−∞
4 4 2
>> syms n z
>> f1 ( n ) = 2^( -n ) ; F1 ( z ) = ztrans ( f1 )
F1 ( z ) =
z /( z - 1/2)
>> f2 ( n ) = n ; F2 ( z ) = ztrans ( f2 )
F2 ( z ) =
z /( z - 1) ^2
>> syms n z
>> f3 ( n ) = n * h e a v i s i d e( n ) +3^ n * h e a v i s i d e ( -n ) ;
>> F3 ( z ) = ztrans ( f3 )
F3 ( z ) =
z /( z - 1) ^2
f u n c t i o n [ F ] = TransZ ( fl , fr )
% TransZ Halla la t r a n s f o r m a d a z b i l a t e r a l de una
% f u n c i ó n f ( n )
% F : Transformada z bilateral
syms n z ;
Fl ( z ) = ztrans ( fl ( - n ) ) ;
Fr ( z ) = ztrans ( fr ( n ) ) ;
F ( z ) = Fl ( z ^( -1) ) + Fr ( z ) ;
end
254 8.5. TRANSFORMADA Z CON MATLAB
La función TransZ se basa en el hecho que toda función f (n) se puede expresar como:
donde fl (t) es una función no causal y fr (n) es una función causal. A continuación se
calcula la transformada z bilateral de la función f3 (n) dada arriba.
>> syms n z
>> f3l ( n ) = 3^ n * h e a v i s i d e( - n ) ; f3r ( n ) = n * h e a v i s i d e( n ) ;
>> F3 ( z ) = TransZ ( f3l , f3r )
F3 ( z ) =
z /( z - 1) ^2 + 1/(3/ z - 1) + 1/2
z 1 1 z z 1
F3 = + + = 2 − z − 3 + 2.
(z − 1)2 3
z −1 2 (z − 1)
Por otra parte, MATLAB también cuenta con el comando iztrans que halla la trans-
formada z inversa. Observe el siguiente ejemplo donde se halla la expresión matemática
de la transformada z inversa de las transformadas F1 (z), F2 (z) y F3 (z) dadas arriba.
>> syms z
>> F1 ( z ) = z /( z - 1/2) ; f1 ( n ) = i z t r a n s( F1 )
f1 ( n ) =
(1/2) ^ n
>> F2 ( z ) = z /( z - 1) ^2; f2 ( n ) = i z t r a n s( F2 )
f2 ( n ) =
n
>> F3 ( z ) = z /( z - 1) ^2 + 1/(3/ z - 1) + 1/2;
>> f3 ( n ) = i z t r a n s( F3 )
f3 ( n ) =
n - 3^ n + k r o n e c k e r D e l t a (n , 0) /2
Tenga presente que MATLAB utiliza el sı́mbolo kroneckerDelta(n, 0) para denotar δ(n).
De esta forma, las transformadas z inversas obtenidas son:
1
f1 (n) = (1/2)n , f2 (n) = n, f3 (n) = n − 3n + δ(n).
2
Note que iztrans halla la transformada z inversa de F (z), asumiendo que F (z) es la
transformada z unilateral de f (n).
CAPÍTULO 8. TRANSFORMADA Z 255
z z z2
F (z) = + = , |z| > 1.
z − 1 (z − 1)2 (z − 1)2
b b
1 x
♦
256 8.6. PROBLEMAS RESUELTOS
∞
X 0
X
F (z) = n 2n cos(2n) u(−n) z −n = n 2n cos(2n) z −n
n=−∞ n=−∞
0
X e2nj + e−2nj
= n 2n z −n
n=−∞
2
0 0
1 X 1 X
= n 2n e2nj z −n + n 2n e−2nj z −n .
2 n=−∞ 2 n=−∞
z e2j z e2j
= − − , |z| < 2.
(z − 2 e2j )2 (z e2j − 2)2
♦
CAPÍTULO 8. TRANSFORMADA Z 257
Luego, calculando a que converge cada una de las series de potencias anteriores, obtenemos:
1 1 1 1
F (z) = −1 + −1 + −
1 − 4−1 z 1 − 2−1 z 1 − 4−1 z −1 1 − 2−1 z −1
| {z } | {z } | {z } | {z }
|z|<4 |z|<2 |z|>1/4 |z|>1/2
b b b b b
1/2 2 x
b
♦
258 8.6. PROBLEMAS RESUELTOS
1 z
Hallemos los desarrollos de Laurent de 1 y , centrados en z0 = 0 y válidos en
z− 2 z − 12
|z| > 12 . Se tiene que:
∞ ∞
1 −1 1 −1
X
−n −n
X 1
1 = z =z 2 z = 2−n u(n)z −(n+1) , |z| > ,
z− 2
1 − 2−1 z −1 n=−∞
2
n=0
∞ ∞
z 1 X X 1
1 = = 2−n z −n = 2−n u(n)z −n , |z| > .
z− 2
1 − 2−1 z −1 n=−∞
2
n=0
z z
Hallemos los desarrollos de Laurent de y , centrados en z0 = 0 y válidos en
z−4 z − 14
1
4 < |z| < 4. Se tiene que:
∞
z −4−1 X
= z = − 4−(n+1) z (n+1)
z−4 1 − 4−1 z n=0
−1
X ∞
X
n −n
= − 4 z =− 4n u(−1 − n)z −n , |z| < 4,
n=−∞ n=−∞
∞ ∞
z 1 X
−n −n
X 1
1 = = 4 z = 4−n u(n)z −n , |z| > .
z− 4
1 − 4−1 z −1 n=−∞
4
n=0
1
De esta forma, el desarrollo de Laurent de F (z) centrado en z0 = 0 y válido en 4 < |z| < 4,
es:
∞ ∞
14 X n 1 X −n
F (z) = − 4 u(−1 − n)z −n + 4 u(n)z −n
15 n=−∞ 15 n=−∞
∞
X 1 −n 14 n 1
= 4 u(n) − 4 u(−1 − n) z −n , < |z| < 4.
n=−∞
15 15 4
Determinemos Z −1 {F1 (z)}. Para ello, hallemos primero el desarrollo de Laurent de F1 (z)
centrado en z0 = 0 y válido en 13 < |z| < 12 . La expansión en fracciones parciales de F1 (z)
está dada por:
73
21 − 17
10
73
15 1 1
F1 (z) = 1 +
+ , < |z| < .
z− 3 z − 12 (z + 2) 3 2
Ahora, se tiene que
∞ ∞
1 z −1 X
−n −(n+1)
X 1
= = 3 z = 3(1−n) u(n)z −n , |z| > ,
z − 31 1 − 3−1 z −1
n=0 n=−∞
3
260 8.7. PROBLEMAS PROPUESTOS
∞ ∞
1 −2 X X 1
= =− 2(n+1) z n = − 2(1−n) u(−n)z −n , |z| < ,
z − 21 1 − 2z n=0 n=−∞
2
∞ ∞
1 2−1 X
−(n+1) n
X
= = 2 z = 2(n−1) u(−n)z −n , |z| < 2.
(z + 2) 1 + 2−1 z n=−∞
n=0
Ası́, usando los series de potencias anteriores obtenemos que el desarrollo de Laurent de
F1 (z) centrado en z0 = 0 y válido en 13 < |z| < 21 es:
∞ ∞
73 X (1−n) 17 X (1−n)
F1 (z) = 3 u(n)z −n + 2 u(−n)z −n
21 n=−∞ 10 n=−∞
∞
73 X (n−1)
+ 2 u(−n)z −n .
15 n=−∞
8.2. La transformada z de una función causal f (n) (f (n) = 0 para n < 0) es:
1 1
F (z) = , |z| > .
(z − 21 )3 2
a) g(n) = 2n f (n − 2)
b) g(n) = nf (n)
c) g(n) = f (−n − 1) − f (−n)
d) g(n) = (n + 1)(1/3)n f (n − 2)
262
RESPUESTA A LOS PROBLEMAS PROPUESTOS 263
2
b) D = {z = x + i y : x = −3/2}, f ′ (z) = 2 eLog(z) Log (z) /z
f ′ (z) = 2x + 2
e) f (z) = ecos z Log (i) ,
c) D = {− 21 , − 21 ′
+ i}, f (z) = 2y + i f ′ (z) = −Log (i) eLog (i) cos(z) sen (z)
2.11 a) D = C, f ′ (z) = 2(z + 1) f) f (z) = esen z Log (z) ,
sen(z)
1 f ′ (z) = eLog(z) sin(z) cos(z) Log (z) +
b) D = C − {0}, f ′ (z) = 1 − z2
z
8 z
c) D = C − {3}, f ′ (z) = (z−3)2 2.22 f ′ (z) = log(i) ez ee log(i) ,
3 2 f ′ (πi/2) = 3πe3π/2 /2
d) D = C − {0, ±i}, f ′ (z) = − 4zz2 (z
+3 z +1
2 +1)2 (
√ 1 1 v ≤ −1/2
e) D = C − ± 2 2 ± 2 i , 2.27 f (D) :
u+v ≥1
z3
f ′ (z) = (z44+1) 2
2.28 a) f (D) = M1√− M2 , donde
f) D = C \ {z√∈ C : Im z = 0 y Re z ≤ M1 : |w| ≥ 2/2,
(
−4 ó z = ± 22 (1 − i)}, u−v ≥0
1
f ′ (z) = (z2 +i)(z+4) − 2zLog (z+4) M2 :
(z 2 +i)2 u+v ≤0
z
g) D = C, f ′ (z) = ez+e b) f (D) = M1 − M2 , donde
′ z z M1 : |w
( − 1 + 2i| ≥ 2,
h) D = C, f (z) = e cos(e )
|w| ≤ 1
2.13 D = {z = x + iy : xy < 1} M2 :
|w + 1 + i| ≤ 1
2.14 a) D = R2 , c) f (D) = M1 − M √2 , donde
v(x, y) = y 2 − x2 − 3x + c M1 : |w + 3| ≥ 8,
( √
b) D = R2 − {(1, 0)}, |w − 1| ≤ 2
−y M2 : √
v(x, y) = x2 −2x+y 2 +1 + c |w + 1| ≤ 2
c) D = R2 − {(0, 0)}, d) f (D) = M1 − M2 , donde
√
2 2
u(x, y) = x(xx2+y +1)
+c M1 : |w
( + 1 − 2i| ≤ 2,
+y 2
|w| ≥ 1
d) D = R2 − {(−1, −2)}, M2 :
u(x, y) = (x+1)x+1 |w − 1 − i| ≥ 1
2 +(y+2)2 + c
264 RESPUESTA A LOS PROBLEMAS PROPUESTOS
√
d) 6 5 d)
4.9 a) Res [f (z)] = 1/3, Res [f (z)] = 5/3 f4 (t)
z=−1 z=2
1
b) Res [f (z)] = 1, Res [f (z)] = −1
z=0 z=1
c) Res [f (z)] = 1/2, Res [f (z)] = 45 (1− 3i ), t
z=1 z=3i a
Res [f (z)] = 45 (1 + 3i )
z=−3i
d) Res [f (z)] = 1/4, Res [f (z)] = 83/4 e)
z=0 z=−4 f5 (t)
e) Res [f (z)] = −1/8, Res [f (z)] = 0 1
z=−3 z=0
f) Res [f (z)] = 1 ··· ···
z=0
√
2 -5π -4π -3π -2π -π π 2π 3π 4π 5π t
g) Res [f (z)] = 4 (1 + i)
z=eπi/4
√ √
1+ 3i 3+ 3
h) Res [f (z)] = −sen 2 12 f)
z=eπi/3
f6 (t)
i) Res [f (z)] = −π
z=π
1
4.10 a) i) 0, ii) 2πi
b) i) π3 (−2 − 4i3 ), ii) 2πi
√
2
√
2 t
−( 2
+ 1) ( 2
− 1)
6π √
c) i) 0, ii) (2√6−1)( 6+1)
, iii) 0
d) i) 0, ii) 0 g)
e) i) 2πi, ii) 6πi f7 (t)
π
f) i) 0, ii) 25 (−14 + 2i) 4
19πi
g) i) 108 , ii) − 19πi
108 3
Capı́tulo 5 2
5.1 a) 1
f1 (t) 1/2
4 1 1 3 2 t
2 2
3
2 h)
f8 (t)
1
4
1 2 3 t 3
b) 2
f2 (t) 1
1
-5 - 9 -4 -3 -2 - 3 -1 1 3 2 3 t
2 2 2
1 2 3 t
c)
f3 (t) 5.2 a) 0
1 b) 12
c) 1/2
a t
d) −9
266 RESPUESTA A LOS PROBLEMAS PROPUESTOS
g) h(t) = 12 3 − e−2t u(t)
h) h(t) = e−2(t−1) u(t − 1) + d)
e−2t (et − 1) u(t) g(n)
b
2
5.4 a) h(t) = −(t + 2)u(t + 2) + 1
2(t + 1) u(t + 1) −
2(t − 1)u(t − 1) +
n
(t − 2) u(t − 2) 1 2 3
b
-1
b) h(t) = −(t + 1)u(t + 2) +
b
2(t + 1) u(t + 1) − 2 u(t)− -2
2(t − 1)u(t − 1) +
(t − 1) u(t − 2)
e)
c) h(t) = 2(t + 2)u(t + 2) − g(n)
2(t + 1) u(t + 1) − 2t u(t)+ 2
2(t − 1)u(t − 1) + p1 (t − 2) b b b b b b b b b
1
··· ···
0, t ≤ −1, b b b b
6n
2
-6 -5 -4 -3 -2 -1 1 2 3 4 5
− (t+1)6(t−2) ,
−1 < t ≤ 0,
2
d) h(t) = (t−1) (t+2)−3t(t−2)
6 , 0 < t ≤ 1, f)
2
g(n)
(t−2) 2
2 , 1 < t < 2,
b b b b b b b b
0, t≥2 1
··· ···
b b b b b
5.5 a)
-6 -5 -4 -3 -2 -1 1 2 3 4 5 6n
g(n)
b
2
5
P 4
P
1 b 5.6 a) h(n) = u(n − k) − u(n + k)
k=1 k=0
-2 b b
c) h(n) = (n + 3) (u(n + 2) − u(n − 1)) +
(3 − n) (u(n − 1) − u(n − 3))
d) No existe
b)
3
g(n) e) h(n) =
P
(n − k)u(n − k)
b
2 k=−2
b
1 f) h(n) = (n − 1) (u(n + 1) − u(n − 4)) +
(n − 2) (u(n) − u(n − 5)) +
2
n P
1 2 k u(k − n)
k=−2
RESPUESTA A LOS PROBLEMAS PROPUESTOS 267
2
P (−1)k 1 1
g) h(n) = 2k+2
[δ(n − k)− 6.5 a) f (t) = 2π πδ(t) − 2 −
k=−2 jt
δ(n − (k + 1)) + δ(n − (k + 2)) b) f (t) = 1 −2
t (1 − cos(2t))
4π
−δ(n − (k + 3))]
1 4 sen t
5
4 P c) f (t) = 2π 1+
P (−1)ℓ t
h) h(n) = ℓ+1 δ(n − (k + ℓ))
k=0 ℓ=0 1
d) f (t) =
πt
Capı́tulo 6 √
6.6 a) A(ω) = 2/ 1 + 4ω 2 , para ω ∈ R
1 φ(ω) = arctan(2ω), para ω ∈ R
6.1 a) F (ω) = 1
2 − jω b) A(ω) = 2|sen ω|/|ω|, para ω 6= 0
∞
2 ejω sen ω φ(ω) =
P
(ω−kπ)pπ/2 ω − (2k+1)π ,
b) F (ω) = 2
ω k=0
16 sen 2 (ω/8) para ω > 0
c) F (ω) = 16 sen 2 (ω/8)
ω2 c) A(ω) = , para ω 6= 0
8 ω2
d) F (ω) = + 8j πδ ′ (ω) φ(ω) = 0
jω
1
1
e) F (ω) = πj δ ′ (ω) − d) A(ω) = 2 , para ω 6= 0
ω2 ω
f) F (ω) = 2π e−jωt0 φ(ω) = π, para ω 6= 0
g) F (ω) = 2π cos(πω) e) A(ω) = 2|sen ω|/|ω|, para ω 6= 0
jω
jω
−2ω,
0 < ω ≤ π2 ,
e 2 2 − 2 e 2 + jω
h) F (ω) = φ(ω) = 2π − 2ω, π2 < ω ≤ π
2 ω2
3π − 2ω, π < ω ≤ 2π
2a
i) F (ω) = 2 (periódica, con periodo 2π, para ω > 0)
a + ω2
1
j) F (ω) = 2π e−a|ω| f) A(ω) = 2 , para ω 6= 0
ω
1 φ(ω) = π, para ω 6= 0
k) F (ω) = + 2πδ(ω)
ω 8
g) A(ω) = , para ω 6= 0
2 sen ω 8j e−jω sen 2 (ω/4) |ω|
l) F (ω) = + φ(ω) = − 2 u(−ω) + π2 u(ω)
π
ω ω2
6.3 a) G1 (ω) = p1 ((−ω − 1)/2)
Capı́tulo 7
b) G2 (ω) = j [δ(ω + 1) − δ(ω − 3)]
c) G3 (ω) = ejω p1 ((ω − 1)/2) (s + 1)e−s + s
7.1 a) F (s) = , Re s < −1
1 s+1
d) G4 (ω) = e−3jω/5 p1 (−(ω + 5)/10)
5 1 − e−5(s−a)
e) G5 (ω) = −2 ejω p1 ((ω − 1)/2) + b) F (s) = , Re s > a
s−a
j(δ(ω + 1) − δ(ω − 3)) 4
f) G6 (ω) = jω p1 ((ω − 1)/2) c) F (s) = 2 , Re s > 4
s − 8s + 32
g) G7 (ω) = −p1 ((ω − 1)/2) + 2(s − 1)
d) F (s) = , −3 < Re s < 5
δ(ω + 1) + 3δ(ω − 3) (s + 3)(s − 3)
h) G8 (ω) = 1
e−j(ω+2)/2 p1 (ω/4) e) No existe
2
i) G9 (ω) = p1 ((ω + 1)/2) f) No existe
j) G10 (ω) = j(δ(ω + 3) − δ(ω − 1)) s−1
g) F (s) = 2 , Re s > 0
s
k) G11 (ω) = j(ej δ(ω + 3) − e−3j δ(ω − 1)) n!
p1 ((ω − 1)/2) h) F (s) = − , Re s < −1
l) G12 (ω) = + πδ(ω) (s + 1)n+1
jω s
i) F (s) = 2 , Re s > 0
m) G13 (ω) = 2j p1 ω+π−1 − p1 ω−π−1 s + a2
2 2
a
n) G14 (ω) = e−jω p1 ω−1 j) F (s) = 2 , Re s < −a
2 a − s2
268 RESPUESTA A LOS PROBLEMAS PROPUESTOS
1 − (1 + 3s)e−3s Capı́tulo 8
7.2 a) F1 (s) = , Re s > 0
s2 z
1 − (1 + 3s)e−3s 1 8.1 a) F (z) = , |z| > 21
b) F2 (s) = + , Re s > 0 z + 12
s2 2 2
−2s
1−e 2z z − a 2a +1
c) F3 (s) = , Re s > 0
s b) F (z) = , |z| > |a|
2 (z − a)(z − a1 )
2e−3s/2 es/2 − 1
d) F4 (s) = , Re s > 0 3 senh(1) z (z 2 − 9)
s2 c) F (z) = , |z| > 1
18e−6s (z 2 − 6 cos(1) z + 9)2
7.3 a) G(s) = , Re s > − 13
(s + 13 )2 3
d) F (z) = , |z| < 3
4e−2s (s + 2) z+3
b) G(s) = , Re s > −1
(s + 1)3 z 11 − 1
e) F (z) = , |z| > 1
2e−3s (3s + 5) z 5 (z − 1)
c) G(s) = , Re s > −1
(s + 1)3 − 53 z 1
2s e−2s f) F (z) = , < |z| < 2
d) G(s) = , Re s > −1 (z − 13 )(z − 2) 3
(s + 1)2
e−2s (4s2 + 6s − 2) − 74 z
e) G(s) = , Re s > −1 g) F (z) = , |z| > 2
(s + 1)3 (z − 14 )(z − 2)
(
0, a < 0, −5 z 1
f) G(s) = s∈C h) F (z) = , |z| < 5
−as
f (a)e , a ≥ 0, 25z 2 + 1
2e−3s s2 es + s(2 + es ) + 4
eα senh(β) z
i) F (z) = ,
g) G(s) = , z 2 − 2eα cosh(β)z + e2α
(s + 1)3
Re s > −1 |z| > eα+β
2e−2s 16
h) G(s) = , Re s > 0 8.2 a) G(z) = , |z| > 1
s(s + 1)3 z 2 (z − 1)
7.4 a) f (0+ ) = 1, lı́m f (t) = 0 2z
t→∞ 1
b) G(z) = , |z| >
+
b) f (0 ) = 0, lı́m f (t) = 0 (z − 12 )3 2
t→∞
c) f (0+ ) = 0, lı́m f (t) = 6 4z 2 (z − 1)
t→∞ 25 c) G(z) = , |z| < 2
(z − 2)2
d) f (0+ ) = ∞, lı́m f (t) = 0
t→∞
10z − 1 1
e) f (0+ ) = 0, lı́m f (t) = 0 d) G(z) = , |z| >
t→∞ 162 z 2(z − 61 ) 6
[2] David, A. (1997). Variable Compleja con Aplicaciones 2da. ed. Addison-Wesley
Iberoamericana, Wilmington, Delaware.
[3] Derrick, W. (1987). Variable Compleja con Aplicaciones. Grupo Editorial Ibe-
roamérica, México.
[4] D.J. Higham and N.J. Higham (2005): Matlab Guide. SIAM, Philadelphia.
[5] James, G. (2002). Matemáticas Avanzadas para Ingenierı́a 2da. ed. Prentice Hall,
México.
[6] Kamen, E.W. y Heck, B.S. (2008). Fundamentos de Señales y Sistemas Usando la
Web y Matlab 3era. ed. Prentice Hall, México.
[7] Lindner, D.K. (1999). Introduction to Signals and Systems. McGraw-Hill, Boston.
[8] Markushévich, A.I. (1967). Teorı́a de las Funciones Analı́ticas, vol. I. Editorial Mir,
Moscú.
[9] Soliman, S.S. y Srinath, M.D. (1999). Señales y Sistemas Continuos y Discretos 2da.
ed. Prentice Hall, México.
269