Matem Atica Aplicada: Universidad Nacional de Rosario
Matem Atica Aplicada: Universidad Nacional de Rosario
Matem Atica Aplicada: Universidad Nacional de Rosario
MATEMÁTICA APLICADA
Ingenierı́a Eléctrica - Ingenierı́a Electrónica
2 Funciones complejas 23
2.1 Funciones complejas de variable compleja . . . . . . . . . . . . . . . . . . . . . . 23
2.2 Funciones analı́ticas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.2.1 Lı́mites y continuidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.2.2 Derivada . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.2.3 Ecuaciones de Cauchy-Riemann . . . . . . . . . . . . . . . . . . . . . . . 31
2.3 Funciones elementales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.3.1 La función exponencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.3.2 La función logaritmo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.3.3 Las funciones trigonométricas seno y coseno . . . . . . . . . . . . . . . . 37
2.4 Transformaciones o mapeos en el plano complejo . . . . . . . . . . . . . . . . . . 38
2.4.1 Algunas transformaciones elementales . . . . . . . . . . . . . . . . . . . . 40
2.4.2 Transformación o mapeo lineal . . . . . . . . . . . . . . . . . . . . . . . . 40
2.4.3 Transformación de inversión . . . . . . . . . . . . . . . . . . . . . . . . . 41
2.4.4 Transformación bilineal, racional o de Möebius . . . . . . . . . . . . . . . 41
2.5 Ejercicios propuestos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3 Integración compleja 49
3.1 Funciones complejas de variable real . . . . . . . . . . . . . . . . . . . . . . . . 49
3.2 La integral compleja . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.2.1 El Teorema de Cauchy-Goursat . . . . . . . . . . . . . . . . . . . . . . . 58
3.2.2 La fórmula integral de Cauchy . . . . . . . . . . . . . . . . . . . . . . . . 64
Capı́tulo 0 2022 2
ÍNDICE GENERAL
Capı́tulo 0 2022 3
ÍNDICE GENERAL
Bibliografı́a 261
Capı́tulo 0 2022 4
Capı́tulo 1
Adoptamos las siguientes convenciones con el fin de simplificar el uso de los números com-
plejos:
Capı́tulo 1 2022 5
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
Capı́tulo 1 2022 6
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
z n = zz...z, (n veces).
Si z 6= 0 entonces
1
z −n =
zn
y z 0 = 1.
1
Ejemplo 1.1.2. Para calcular (1 + i)−2 , notemos primero que (1 + i)−2 = y como
(1 + i)2
(1 + i)2 = (1 + i)(1 + i) = 1 + 2i + i2 = 1 + 2i − 1 = 2i
1 −2i −2i 1
entonces se tiene que (1 + i)−2 = = = = − i.▲
2i 2i(−2i) 4 2
Demostración. Ejercicio.
Capı́tulo 1 2022 7
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
Según esta definición se establece una correspondencia biunı́voca entre los pares (x, 0) y los
números reales x. De esta manera el conjunto de los números reales R se identifica con un
subconjunto C′ de C y esta correspondencia se establece también, entre sumas en R y sumas
en C′ y entre productos en R y productos en C′ .
La unidad imaginaria se define por i = (0, 1), que verifica
Por otra parte (x, y) = (x, 0) + (0, y) = (x, 0) + (0, 1) (y, 0) de lo que se tienen que:
(x, y) = x + iy.
Esta última igualdad muestra la equivalencia entre dos formas de expresión de un mismo
números complejo, como par ordenado z = (x, y) y en forma binómica z = x + iy.
P ↔ z = (x, y) = x + iy.
En virtud de esta correspondencia biunı́voca entre puntos del plano y números complejos,
se identifica a C con el plano que por este motivo se llama plano complejo.
En este plano, el eje x de las abscisas se llama eje real y el eje y de las ordenadas, se llama
eje imaginario.
−→
A cada número complejo z = (x, y) = x + iy le corresponde el vector OP , denominado
vector asociado a z, de origen O y extremo el punto P llamado afijo de z, como se observa
en la Figura 1.1.
Capı́tulo 1 2022 8
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
A partir de este concepto, se tienen las siguientes propiedades de números complejos conjugados:
1. z = z ⇔ z ∈ R 5. z1 + z2 = z1 + z2
2. z + z = 2 Re(z) 6. z1 z2 = z1 z2
7. −z = −z
3. z − z = 2i Im(z)
1 1
4. zz = Re2 (z) + Im2 (z) 8. = si z 6= 0
z z
Demostración. Ejercicio.
1. z1 − z2 = z1 − z2
z1 z1
2. = si z2 6= 0
z2 z2
Capı́tulo 1 2022 9
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
Demostración. Probaremos sólo (1) mientras que (2) queda como ejercicio.
Dados los complejos z1 y z2 , entonces
z1 − z2 = z1 + (−z2 ) = z1 + −z2 = z1 − z2 .
Capı́tulo 1 2022 10
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
1. |z| = | − z| = |z| = | − z|
2. zz = |z|2
Demostración. Probaremos aquı́ sólo la propiedad (6) conocida con el nombre de desigualdad
triangular. Las demás quedan como ejercicio.
|z1 + z2 |2 = (z1 + z2 ) z1 + z2 = z1 z1 + z1 z2 + z2 z1 + z2 z2 = |z1 |2 + z1 z2 + z1 z2 + |z2 |2 =
= |z1 |2 + 2 Re (z1 z2 ) + |z2 |2 ≤ |z1 |2 + 2 |z1 z2 | + |z2 |2 = |z1 |2 + 2 |z1 | |z2 | + |z2 |2 = (|z1 | + |z2 |)2
Luego, se tiene que |z1 + z2 | ≤ |z1 | + |z2 |.
−π < arg(z) ≤ π.
1. x > 0, y = 0 ⇒ arg(z) = 0
2. x < 0, y = 0 ⇒ arg(z) = π
Capı́tulo 1 2022 11
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
π
3. x = 0, y > 0 ⇒ arg(z) =
2
π
4. x = 0, y < 0 ⇒ arg(z) = −
2
arctan(y/x) si z ∈ Ic ó z ∈ IV c
5. x 6= 0, y =
6 0 ⇒ arg(z) = arctan(y/x) + π si z ∈ IIc
arctan(y/x) − π si z ∈ IIIc
π 3
Ejemplo 1.2.4. arg(−1 − i) = arctan(1) − π = − π = − π. ▲
4 4
Observación 1.2.4. Las fórmulas anteriores no son válidas para argumentos principales. En
π
efecto, sean z1 = −1 y z2 = i. Se tiene que arg(z1 z2 ) = arg(−i) = − mientras que arg(z1 ) +
2
π 3
arg(z2 ) = π + = π.
2 2
Capı́tulo 1 2022 12
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
z = r(cos θ + isen θ)
denominada forma trigonométrica de z, y
z = rθ
llamada forma polar de z. Notar que r = |z| y θ ∈ Arg(z).
√ π
Ejemplo 1.3.1. Sea z = 1 − i. Se tiene que |z| = 2 y arg(z) = − . Entonces las formas
4
trigonométrica y polar de z son:
√ π π
z = 2 cos − + isen −
4 4
y √
z = 2− π4
respectivamente. ▲
Observación 1.3.1 (Igualdad en forma polar). Sean z1 = rθ1 y z2 = rθ∗2 . Entonces z1 = z2 si
y sólo si r = r∗ y θ1 − θ2 = 2kπ, k ∈ Z.
Demostración. Ejercicio.
Capı́tulo 1 2022 13
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
donde k = 0, 1, 2, ..., n − 1.
√
Ejemplo 1.3.2. Para obtener 3 −8i, comencemos expresando z = −8i en forma polar como
z = 8− π2 . Entonces las raı́ces cúbicas de z están dadas por
√
3
wk = 8 − π2 +2kπ
3
Capı́tulo 1 2022 14
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
z = reiθ
1. z1 z2 = r1 r2 ei(θ1 +θ2 )
z1 r1
2. = ei(θ1 −θ2 )
z2 r2
3. n ∈ Z ⇒ z1n = r1n einθ1
√ √ θ1 +2kπ
4. n z1 = n r1 ei n con k = 0, 1, 2, ..., n − 1.
C∗ = C ∪ {∞}
∞ + z = z + ∞ = ∞, ∀z ∈ C∗
∞ − z = z − ∞ = ∞, ∀z ∈ C
z
= 0, ∀z ∈ C
∞
∞
= ∞, ∀z ∈ C
z
z
= ∞, ∀z 6= 0, z ∈ C∗
0
Capı́tulo 1 2022 15
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
El módulo del infinito se define como |∞| = +∞ pero carece de sentido hablar del argumento
del punto del infinito, ası́ como de las operaciones:
∞ 0
∞ − ∞, 0.∞, , .
∞ 0
Fue Riemann quien dio una respuesta satisfactoria a la representación del punto del infinito,
con la que se conoce como la proyección estereográfica y que no es otra cosa que uno de
los métodos utilizados en la confección de mapas, es decir, en la representación de superficies
esféricas.
De esta forma podemos representar a todo punto del plano complejo sobre la esfera sin
el punto correspondiente al polo norte, punto este último donde viene representado el punto
del infinito ∞. Nótese que los números complejos de módulo uno quedan representados en el
ecuador, en el mismo punto sobre el que se encontraban. Los del interior de la circunferencia
(de módulo menor que uno), quedan representados en el hemisferio sur y los de módulo mayor
Capı́tulo 1 2022 16
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
que uno (los de fuera de la circunferencia) quedan representados en el hemisferio norte. Ası́
pues, conforme nos alejamos del origen en el plano complejo, nos acercamos al polo norte en
su representación en la esfera. Luego, hemos conseguido representar cada punto de C∗ por un
único punto de la esfera, y recı́procramente. Ası́ representado, el plano complejo ampliado suele
llamarse esfera de Riemann.
Capı́tulo 1 2022 17
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
7. La evolución en el tiempo de una tensión alterna senoidal de valor eficaz 220 V (Voltios)
y frecuencia de 50 Hz (Hertz), como la que existe en la instalación eléctrica de nuestros
hogares, puede representarse de la siguiente forma
√
220 cos 2π |{z}
v(t) = 2 |{z} 50 t + ϕ [V ] , (1.2)
valor ef icaz f recuencia
Capı́tulo 1 2022 18
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
8. Calcule:
9 √
(a) (1 + i)14 (c) 3 π3 (e) i
π −25 √ 1
(b) 2e− 4 i (d) 4 −1 (f) (−1 + i) 3
(a) Halle las raı́ces de p(λ) para los valores m = 0, 46 10−3 kg, b = 0, 0075 Ns/m y
k = 18, 6 N/m.
(b) El movimiento analizado será estable si todas las raı́ces de p(λ) poseen parte real
negativa. Evalúe si para los valores de m, b y k dados se tiene estabilidad.
(c) Demuestre que siempre que m, b y k sean positivos, se tendrá estabilidad.
(a) S1 = {z ∈ C : |z| = 4}
n π πo
(b) S2 = z ∈ C : − < arg(z) ≤
2 4
n πo
(c) S3 = z ∈ C : 1 < |z| ≤ 3, −π < arg(z) ≤
2
(d) S4 = {z ∈ C : |Re(z)| ≤ 3, −2 ≤ Im(z) ≤ 1}
z − 3
(e) S5 = z ∈ C : <2
z + 3
(f) S6 = {z ∈ C : Re (z 2 ) = 1}
(g) S7 = {z ∈ C : |z + 2 − 3i| + |z − 2 + 3i| ≤ 10}
Capı́tulo 1 2022 19
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
Determine si los conjuntos Si dados en el ejercicio (11) son abierto, cerrado, acotado,
conexo, simplemente conexo y/o un dominio.
Capı́tulo 1 2022 20
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
Capı́tulo 1 2022 21
CAPÍTULO 1. EL SISTEMA DE LOS NÚMEROS COMPLEJOS
Capı́tulo 1 2022 22
Capı́tulo 2
Funciones complejas
f : S ⊆ C → C, w = f (z)
Ejemplo 2.1.1. Sea la función f : C → C tal que f (z) = z 2 . Entonces f (1 + 2i) = (1 + 2i)2 =
−3 + 4i, por lo que se dice que −3 + 4i es la imagen de 1 + 2i por la función f . ▲
Capı́tulo 2 2022 23
CAPÍTULO 2. FUNCIONES COMPLEJAS
Observación 2.1.1. Si una función f está dada sólo por su ley, sobreentenderemos que se
toma por dominio el mayor subconjunto de C donde f (z) existe, es decir
z
Figura 2.1: Representación del dominio de f (z) = z 2 +1
f (x + iy) = u + iv.
Cada número real u y v depende de las variables reales x e y, luego f (z) puede ser expresado
en términos de un par de funciones reales de las variables reales x e y
Capı́tulo 2 2022 24
CAPÍTULO 2. FUNCIONES COMPLEJAS
Por el contrario, sea por ejemplo la función dada por f (z) = 7x + 3iy. Si queremos expresar
a f (z) en términos de la variable compleja z entonces dado que Re(z) = x y Im(z) = y tenemos
que
f (z) = 7Re(z) + 3iIm(z).
donde p y q son polinomios. Este tipo de funciones tienen por dominio C − {z ∈ C : q(z) = 0}.
iz 1
Ejemplo 2.2.1. Sea f (z) = . Probaremos que lim f (z) = − .
2 z→i 2
Sea ε > 0. Nótese que:
1 iz 1 1 1
f (z) + = + = |i(z − i)| = |z − i| .
2 2 2 2 2
Capı́tulo 2 2022 25
CAPÍTULO 2. FUNCIONES COMPLEJAS
Teorema 2.2.1. Sean la función dada por f (z) = u(x, y) + iv(x, y) y z0 = x0 + iy0 .
Entonces,
lim f (z) = u0 + iv0
z→z0
si y sólo si
lim u(x, y) = u0 y lim v(x, y) = v0 .
(x,y)→(x0 ,y0 ) (x,y)→(x0 ,y0 )
luego
u(x, y) = x2 − y 2 , v(x, y) = 2xy.
Dado que
lim u(x, y) = lim x2 − y 2 = −1
(x,y)→(0,1) (x,y)→(0,1)
y
lim v(x, y) = lim 2xy = 0,
(x,y)→(0,1) (x,y)→(0,1)
Capı́tulo 2 2022 26
CAPÍTULO 2. FUNCIONES COMPLEJAS
en el caso en que alguno o ambos, z0 ó w0 , sean sustituidos por el punto del infinito ∞.
Definición (Lı́mite infinito). La afirmación lim f (z) = ∞ significa que para cada
z→z0
número positivo M , existe un número positivo δ tal que si 0 < |z − z0 | < δ entonces
|f (z)| > M .
1 1
Puesto que la definición anterior puede escribirse como 0 < |z − z0 | < δ ⇒ − 0 < ,
f (z) M
entonces tenemos que:
1
lim f (z) = ∞ ⇔ lim = 0.
z→z0 z→z0 f (z)
Capı́tulo 2 2022 27
CAPÍTULO 2. FUNCIONES COMPLEJAS
iz + 3
Ejemplo 2.2.3. Probemos que lim = ∞.
z→−1 z + 1
iz + 3 1 z+1
Para ello, sea f (z) = . Entonces dado que lim = lim = 0 se tiene por
z+1 z→−1 f (z) z→−1 iz + 3
lo tanto que:
iz + 3
lim = ∞.▲
z→−1 z + 1
Ejemplo 2.2.4. La función dada por f (z) = |z|2 = x2 + y 2 es continua en C pues las funciones
componentes de f , u(x, y) = x2 + y 2 y v(x, y) = 0 son continuas en R2 . ▲
Observación 2.2.1. Las funciones polinómicas y racionales son continuas en sus dominios.
Capı́tulo 2 2022 28
CAPÍTULO 2. FUNCIONES COMPLEJAS
2.2.2 Derivada
Definición (Derivada en un punto). Sea f una función cuyo dominio contiene un
entorno de z0 . Se dice que f es derivable en z0 si existe
f (z) − f (z0 )
lim . (2.1)
z→z0 z − z0
En tal caso, al resultado del lı́mite se lo denota por f ′ (z0 ) y se llama derivada de f en
z0 .
4w = f (z + 4z) − f (z),
que denota el cambio en el valor de f correspondiente a un cambio 4z en el punto en que
evaluamos f . Entonces (2.1) se puede expresar como:
dw 4w
= lim .
dz △z→0 4z
f′ : Ω → C / w = f ′ (z)
Ejemplo 2.2.5. Determinemos dónde es derivable la función definida por f (z) = z 2 . Como
Ejemplo 2.2.6. Estudiemos ahora dónde es derivable la función dada por f (z) = z. Para ello,
consideremos z = x + iy ∈ C y 4z = 4x + i4y. Entonces,
Capı́tulo 2 2022 29
CAPÍTULO 2. FUNCIONES COMPLEJAS
f (z + 4z) − f (z)
Luego, no existe lim para ningún z ∈ C por lo que f no es derivable en
△z→0 4z
todo C.
Demostración. Ejercicio.
Capı́tulo 2 2022 30
CAPÍTULO 2. FUNCIONES COMPLEJAS
Demostración. La prueba de este teorema es similar al caso de funciones reales de una variable
real.
′ 0z −n − 1(−n)z −n−1
f (z) = −2n
= nz n−1 , z 6= 0.▲
z
Además,
f ′ (z0 ) = ux (x0 , y0 ) + ivx (x0 , y0 ) = vy (x0 , y0 ) − iuy (x0 , y0 ).
Capı́tulo 2 2022 31
CAPÍTULO 2. FUNCIONES COMPLEJAS
Si en particular, 4z → 0 a lo largo del eje real (4y = 0) entonces (2.3) se puede expresar
como
Con una hipótesis adicional, se obtiene el recı́proco de este teorema, que enunciamos sin
demostración:
Ejemplo 2.2.9. Usaremos las condiciones de Cauchy-Riemann para determinar dónde es de-
rivable la función dada por f (z) = zz.
Capı́tulo 2 2022 32
CAPÍTULO 2. FUNCIONES COMPLEJAS
Sean z = x + iy, entonces la función f puede expresarse como f (z) = x2 + y 2 por lo que
u(x, y) = x2 + y 2 y v(x, y) = 0.
Las funciones ux = 2x, uy = 2y, vx = vy = 0 son continuas para todo (x, y) ∈ R2 y como
ux = vy 2x = 0 x=0
⇔ ⇔
uy = −vx 2y = 0 y=0
Capı́tulo 2 2022 33
CAPÍTULO 2. FUNCIONES COMPLEJAS
f : C → C / f (z) = ez .
Si z = x+iy entonces f (z) = ex+iy = ex eiy = ex cos y +iex sen y, por lo que u(x, y) = ex cos y
y v(x, y) = ex sen y. Como las ecuaciones de Cauchy-Riemann
se satisfacen para todo (x, y) ∈ R2 , y las derivadas parciales de primer orden de las funciones
u y v son continuas en R2 , entonces se tiene que f es una función entera. Además,
para todo z ∈ C.
p
Para cualquier z ∈ C se tiene que |ez | = |ex (cos y + isen y)| = ex cos2 y + sen 2 y = ex > 0
por lo que se tiene que el módulo de ez es ex y ez 6= 0 para todo z ∈ C.
Si escribimos la expresión ez = ex eiy como ez = reiθ observamos que r = ex y θ = y por lo
que un argumento de ez es y. Se tiene ası́ que
Arg(ez ) = y + 2kπ
donde k ∈ Z.
También, se tiene que para cualesquiera z1 y z2 ∈ C,
ez1
z1 z2
e e =e z1 +z2
, z
= ez1 −z2 .
e 2
Por otro lado, como e2πi = 1 entonces para cualquier z ∈ C se tiene que ez+2πi = ez e2πi = ez ,
por lo que la función f es periódica de perı́odo fundamental 2πi.
f (z) = ez = w (2.10)
ex eiy = reiθ .
Capı́tulo 2 2022 34
CAPÍTULO 2. FUNCIONES COMPLEJAS
Luego, por la igualdad de números complejos se tiene que ex = r y eiy = eiθ o equivalente-
mente
x = ln r , y = θ + 2kπ (k ∈ Z).
Por lo tanto, la ecuación (2.10) se cumple si
z = ln r + i(θ + 2kπ)
iz = Log (1 + i)
iz = ln |1 + i| + iArg(1 + i)
√ π
iz = ln 2 + i + 2kπ , k ∈ Z
π 4 √
z = + 2kπ − i ln 2, k ∈ Z.▲
4
u(r, θ) = ln r , v(r, θ) = θ
Capı́tulo 2 2022 35
CAPÍTULO 2. FUNCIONES COMPLEJAS
Si (r, θ) es arbitrariamente cercano al punto (r, π), entonces existirán valores de v próximos
a π y también otros en los que los valores de v están próximos a −π, por lo que se puede
concluir que h no es continua en el semieje real negativo.
Podemos concluir hasta el momento que como la función h no es continua en z = 0 y en los
puntos de la forma (r, π) con r ≥ 0 y θ = π, entonces tampoco será derivable allı́.
1
Dado que las primeras derivadas parciales de u y v son ur = , vθ = 1 y vr = uθ = 0, son
r
continuas en (r, θ) con r > 0 y −π < θ < π y se satisfacen las ecuaciones de Cauchy-Riemann
en polares (ver ejercicio 3 de la Sección 2.5)
1 1
ur = vθ , uθ = −vr
r r
entonces podemos concluir que h es continua y analı́tica en:
Ω = C − {z ∈ C : Re(z) ≤ 0, Im(z) = 0}
para todo z ∈ Ω.
Ω = C − {z ∈ C : Re(z − i + 1) ≤ 0, Im(z − i + 1) = 0} .
Capı́tulo 2 2022 36
CAPÍTULO 2. FUNCIONES COMPLEJAS
2. eLog z1 = z1
Observación 2.3.1. Las propiedades enunciadas en el Teorema 2.3.1 no son válidas para
logaritmos principales. Por ejemplo, observemos que:
π
log(−3i) = ln |−3i| + i arg(−3i) = ln 3 − i
2
mientras que
3
log(−3) + log(i) = ln 3 + i arg(−3) + ln 1 + i arg(i) = ln 3 + πi.
2
eiz − e−iz
sen z =
2i
y la función trigonométrica coseno de z como
eiz + e−iz
cos z = .
2
Capı́tulo 2 2022 37
CAPÍTULO 2. FUNCIONES COMPLEJAS
Las funciones sen z y cos z son enteras pues son suma y composición de funciones enteras.
Además,
d 1 d iz 1 eiz + e−iz
(sen z) = e − e−iz = ieiz + ie−iz = = cos z
dz 2i dz 2i 2
para cualquier z ∈ C. Análogamente se prueba que para todo complejo z se tiene que:
d
(cos z) = −sen z.
dz
3. sen z = 0 ⇔ z = nπ con n ∈ Z
π
4. cos z = 0 ⇔ z = + nπ con n ∈ Z
2
5. sen z = sen x cosh y + i cos xsenhy
Demostración. Ejercicio.
Capı́tulo 2 2022 38
CAPÍTULO 2. FUNCIONES COMPLEJAS
Entonces,
w = (1 + i)(x + iy) = (x − y) + i(x + y)
Capı́tulo 2 2022 39
CAPÍTULO 2. FUNCIONES COMPLEJAS
para cualquier x ∈ R.
Estas últimas son las ecuaciones paramétricas de una recta en el plano uv que contiene al
punto (−1, 1) y tiene la dirección del vector (1, 1) como se muestra en la Figura 2.4. ▲
Traslación
Sea w = z + β con β ∈ C. Entonces los puntos del plano z se trasladan en la dirección del
vector asociado al número complejo β en el plano w.
Rotación
Sea w = eiθ0 z con θ0 ∈ R. Entonces los puntos del plano z rotan un ángulo θ0 respecto al origen
en el plano w.
Si θ0 > 0, la rotación es en sentido positivo (en contra de las agujas del reloj), mientras que
si θ0 < 0, la rotación es en sentido negativo (a favor de las agujas del reloj).
Dilatación-contracción
Sea w = az con a > 0. Entonces si a > 1, puntos del plano z se dilatan un factor a en el plano
w en la dirección del vector asociado a z. Si por el contrario, 0 < a < 1, puntos del plano z se
1
contraen un factor en el plano w en la dirección del vector asociado a z.
a
Ejemplo 2.4.2. La tranformación del ejemplo 2.4.1, w = (1 + i)z se puede expresar como:
√ π
w = 2e 4 i z
w = αz + β
Capı́tulo 2 2022 40
CAPÍTULO 2. FUNCIONES COMPLEJAS
w1 = αz = aeiθ0 z
con a ∈ R+ , θ0 ∈ R es decir se trata de una dilatación/contracción y rotación.
Luego, la transformación
w = w1 + β
es una traslación en la dirección del vector asociado al complejo β.
Capı́tulo 2 2022 41
CAPÍTULO 2. FUNCIONES COMPLEJAS
Si pedimos que,
T (∞) = ∞ si c = 0
y
a d
T (∞) = y T − = ∞ si c 6= 0
c c
esto hace entonces que T sea continua en todo el plano complejo extendido C∗ .
Capı́tulo 2 2022 42
CAPÍTULO 2. FUNCIONES COMPLEJAS
∂u ∂v ∂u ∂v
= y =− .
∂x ∂y ∂y ∂x
Capı́tulo 2 2022 43
CAPÍTULO 2. FUNCIONES COMPLEJAS
∂u ∂u ∂x ∂u ∂y
= + = ux cos θ + uy sen θ (2.11)
∂r ∂x ∂r ∂y ∂r
∂v ∂v ∂x ∂v ∂y
= + = vx cos θ + vy sen θ (2.12)
∂r ∂x ∂r ∂y ∂r
∂u ∂u ∂x ∂u ∂y
= + = −ux rsen θ + uy r cos θ (2.13)
∂θ ∂x ∂θ ∂y ∂θ
∂v ∂v ∂x ∂v ∂y
= + = −vx rsen θ + vy r cos θ. (2.14)
∂θ ∂x ∂θ ∂y ∂θ
∂u 1 ∂v ∂v 1 ∂u
= y =− .
∂r r ∂θ ∂r r ∂θ
Capı́tulo 2 2022 44
CAPÍTULO 2. FUNCIONES COMPLEJAS
4. Funciones armónicas
Una función real de dos variables reales ϕ = ϕ(x, y) se dice armónica en un dominio D
del plano real si tiene derivadas parciales continuas de primer y segundo orden en dicho
dominio y si satisface la ecuación de Laplace
Consideremos una función compleja de variable compleja f (z) = u(x, y)+iv(x, y) analı́tica
en el dominio D del plano complejo.
Si dos funciones reales de dos variables reales u = u(x, y) y v = v(x, y) son armónicas
en un dominio del plano real D y sus derivadas parciales de primer orden satisfacen las
ecuaciones de Cauchy-Riemann en D, se dice que v es armónica conjugada de u.
Es evidente que si una función f (z) = u(x, y) + iv(x, y) es analı́tica en un dominio del
plano complejo D, entonces v es armónica conjugada de u. Recı́procamente, si v es una
armónica conjugada de u en un dominio D, la función f (z) = u(x, y)+iv(x, y) es analı́tica
en D.
De esto último se tiene el siguiente teorema: Una función compleja de variable compleja
f (z) = u(x, y) + iv(x, y) es analı́tica en un dominio D del plano complejo si y sólo si v
es una armónica conjugada de u.
Capı́tulo 2 2022 45
CAPÍTULO 2. FUNCIONES COMPLEJAS
(a) Verifique que la función dada es armónica en su dominio. Halle una función u o v,
según corresponda, de manera que f (z) = u(x, y)+iv(x, y) resulte analı́tica. Exprese
a f (z) en términos de la variable compleja z.
i. v(x, y) = 2x − 2xy con la condición u(0, 0) = 0
ii. u(x, y) = e−x (xsen y − y cos y) con la condición v(0, 0) = 0
(b) Sea f (z) = u(x, y) + iv(x, y) una función analı́tica en un dominio D. Demuestre
que las familias de curvas de ecuaciones u(x, y) = c y v(x, y) = k, donde c y k son
constantes reales arbitrarias, son ortogonales (es decir, cada miembro de una familia
tienen tangentes perpendiculares en su punto de intersección con cada miembro de
la otra familia).
(c) Para cada una de las familias de curvas de ecuaciones dadas, halle las trayectorias
ortogonales:
i. x3 y − xy 3 = c
ii. e−x (xsen y − y cos y) = c
6. Pruebe que:
π
(a) log (−ei) = 1 − i (c) Log (e) = 1 + 2kπi, k ∈ Z
2 √
1 π (d) Log −1 + 3i = ln 2 + 2 n + 13 πi,
(b) log (1 − i) = ln 2 − i k∈Z
2 4
√
7. Sea la función f : C → C / f (z) = z. Pruebe que f es analı́tica en:
C − {z ∈ C : Re(z) ≤ 0, Im(z) = 0} .
1
Demuestre además que la derivada está dada por f ′ (z) = √ .
2 z
8. Demuestre el Teorema 2.3.2.
(a) Explique por qué las funciones senhz y cosh z son enteras. Pruebe que:
d d
(senhz) = cosh z y (cosh z) = senhz.
dz dz
(b) Demuestre las siguientes proposiciones:
i. senhz y cosh z son periódicas de perı́odo fundamental 2πi.
ii. cosh2 z − senh2 z = 1, cosh(−z) = cosh z, senh(−z) = −senhz para cualquier
z ∈ C.
iii. senhz = 0 ⇔ z = nπi
π con n∈ Z
iv. cosh z = 0 ⇔ z = + nπ i con n ∈ Z
2
v. senhz = −isen (iz), cosh z = cos (iz) para cualquier z ∈ C.
10. Determine la región de analiticidad de las siguientes funciones. Obtenga una expresión
para f ′ .
3z + i √
(a) f (z) = z 3 − 2iz 2 + i (c) f (z) = 3 (e) f (z) = z − 1 + 3i
z −1
3 1
log(z + 4)
(b) f (z) = (z 2 + iz) (d) f (z) = 3e z − i (f) f (z) =
z2 + 1
Calcule, si es posible:
1 − cos z 1 − ez (c) lim
sen z
(a) lim (b) lim z→π z − π
z→0 sen z z→0 z2
12. Halle las imágenes por la función f de las curvas de ecuaciones Re(z) = k, Im(z) = c,
|z| = r y arg(z) = α, donde k, c, r y α son constantes reales, siendo:
13. Sea el mapeo lineal w = αz + β donde α y β son constantes complejas, el cual transforma
el punto z = 1 + i en el punto w = i, y el punto z = 1 − i en el punto w = −1.
Capı́tulo 2 2022 47
CAPÍTULO 2. FUNCIONES COMPLEJAS
√
14. Halle la imagen de la transformación w = z de las curvas de ecuaciones arg(z) = c y
|z| = k con c, k ∈ R.
15. Determine las imágenes de las curvas de ecuaciones Re(z) = k y Im(z) = c, por la
transformación w = ez .
16. En cada uno de los siguientes casos, determine la imagen de la región S por la transfor-
mación w dada. Represente geométricamente la región S y su imagen.
Capı́tulo 2 2022 48
Capı́tulo 3
Integración compleja
Observación 3.1.1. Algunos resultados que serán de utilidad son las siguientes fórmulas de
derivación y se deja su demostración como ejercicio:
• Si w1 = w1 (t) y w2 = w2 (t) son funciones complejas derivables en la variable real t, y c
es una constante compleja, entonces
Capı́tulo 3 2022 49
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
′
– (w1 + w2 )′ (t) = w1′ (t) + w2′ (t) w1 w1′ (t)w2 (t) − w1 (t)w2′ (t)
– (t) =
w2 w22 (t)
– (w1 w2 )′ (t) = w1′ (t)w2 (t) + w1 (t)w2′ (t) – (cw1 )′ (t) = cw1′ (t)
d ct
• Si c es una constante compleja entonces, [e ] = cect .
dt
Z Z Z
1 1 1
2
Ejemplo 3.1.1. (1 + it) dt = 2
1−t 2
dt + i 2t dt = + i. ▲
0 0 0 3
Demostración. Sólo probaremos aquı́ la propiedad (6). Supongamos que w(t) = x(t) + iy(t) y
W (t) = X(t) + iY (t) son continuas en [a, b].
Capı́tulo 3 2022 50
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
Si W ′ (t) = w(t) para t ∈ [a, b], entonces X ′ (t) = x(t) y Y ′ (t) = y(t). Luego,
Z b Z b Z b
w(t) dt = x(t) dt + i y(t) dt = [X(b) − X(a)] + i [Y (b) − Y (a)] =
a a a
Z π
π
1 2it 1
Ejemplo 3.1.2. e 2it
dt = − ie = − i e2πi − 1 = 0. ▲
0 2 0 2
Capı́tulo 3 2022 51
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
Ejemplo 3.1.4. Para representar gráficamente la curva C de ecuación z(θ) = eiθ , θ ∈ [0, 2π),
expresémosla en la forma z(θ) = cos θ + isen θ donde θ ∈ [0, 2π) por lo que la curva C tiene
por ecuaciones paramétricas:
x = cos θ
, θ ∈ [0, 2π)
y = sen θ
Capı́tulo 3 2022 52
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
3. C es una curva suave si z ′ es continua en [a, b] y z ′ (t) 6= 0 para todo t ∈ (a, b).
4. C es una curva suave a trozos o seccionalmente suave si es una curva que consiste
en un número finito de arcos suaves unidos por sus extremos.
Capı́tulo 3 2022 53
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
Observación 3.1.2. Sea la función continua z(t) = x(t) + iy(t) con t ∈ [a, b] cuya representa-
ción gráfica es una curva C en el plano complejo. Se puede probar que:
1. La representación paramétrica usada para la curva C no es única.
p
2. Si C es seccionalmente suave, entonces |z ′ (t)| = [x′ (t)]2 + [y ′ (t)]2 es integrable en [a, b].
Luego, Z Z b b
|z ′ (t)| dt = ds = long(C).
a a
3. La long(C) dada por la integral del ı́tem (2) es invariante bajo ciertos cambios de repre-
sentación de dicha curva C.
Teorema 3.2.1. Sean C una curva seccionalmente suave de ecuación z = z(t) con
t ∈ [a, b] y f , g funciones complejas de variable compleja continuas sobre C. Entonces
se verifican las siguientes propiedades:
Z Z
1. k ∈ C ⇒ kf (z) dz = k f (z) dz.
C C
Z Z Z
2. [f (z) + g(z)] dz = f (z) dz + g(z) dz.
C C C
Z Z
3. f (z) dz = − f (z) dz donde −C es una curva con el mismo conjunto de
−C C
puntos de C pero recorrido en sentido contrario al de C.
Capı́tulo 3 2022 54
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
es evidente que Z Z Z
f (z) dz = f (z) dz + f (z) dz.
C C1 C2
5. Como por hipótesis f es continua sobre C y dado que C es una curva seccionalmente
suave, entonces z(t) es continua y su derivada z ′ (t) es seccionalmente continua. Luego,
se tiene que la función real |f (z(t))| es continua en el intervalo cerrado [a, b], por lo que
existe M > 0 tal que |f (z(t))| ≤ M para cualquier t ∈ [a, b]. Entonces,
Z Z b Z b
f (z) dz = f (z(t)) z ′
(t) dt≤ |f (z(t))| |z ′ (t)| dt
C a a
Z b
≤M |z ′ (t)| dt = M long(C).
a
Z
Ejemplo 3.2.2. Calculemos f (z) dz donde f (z) = y − x − i3x2 con z = x + iy siendo
C
Capı́tulo 3 2022 55
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
Entonces, Z Z Z
f (z) dz = f (z) dz + f (z) dz = I1 + I2 .
C C1 C2
Para calcular I1 , se tiene que f (z(t)) = t y z ′ (t) = i con t ∈ [0, 1], entonces
Z Z 1
1
I1 = f (z) dz = ti dt = i.
C1 0 2
Capı́tulo 3 2022 56
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
Observación 3.2.2. Nótese que una primitiva es, necesariamente, una función analı́tica.
Capı́tulo 3 2022 57
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
w(t) = F (z(t))
es derivable en t. Además, sean F (z) = u(x, y) + iv(x, y) y z(t) = x(t) + iy(t). Entonces,
o equivalentemente,
w′ (t) = F ′ (z(t)) z ′ (t) = f (z(t)) z ′ (t)
donde a ≤ t ≤ b. Por la propiedad 6 del Teorema 3.1.1, se tiene que
Z Z b
f (z) dz = f (z(t)) z ′ (t) dt = F (z(t))|ba = F (z(b)) − F (z(a)) .
C a
Como z(b) = z2 y z(a) = z1 , el valor de esta integral de contorno es F (z2 ) − F (z1 ); y este
valor es ciertamente independiente del contorno C, en tanto esté contenido en Ω y una z1 con
z2 . Es decir, Z Z z2
f (z) dz = f (z) dz = F (z)|zz21 = F (z2 ) − F (z1 ) .
C z1
Z
Ejemplo 3.2.3. Si queremos calcular z 2 dz para toda curva seccionalmente suave C desde
C
el punto z = 0 hasta el punto z = 1 + i, observemos que la función
Z f (z) = z 2 admite como
1
primitiva a F (z) = z 3 (que es continua), entonces la integral z 2 dz es independiente de la
3 C
trayectoria. Luego,
Z Z 1+i 1+i
1 3 2 2
2
z dz = z dz = z = − + i.▲
2
C 0 3 0 3 3
Capı́tulo 3 2022 58
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
que garantizan que el valor de la integral de f a lo largo de cualquier contorno cerrado simple
seccionalmente suave es cero.
Antes de demostrar este teorema, recordemos lo que establece el Teorema de Green:
Sea C una curva simple cerrada seccionalmente suave orientada positivamente del plano,
y sea D la región interior a C. Si las funciones reales de dos variables reales P = P (x, y),
Q = Q(x, y), Py y Qx son continuas en D y sobre la frontera C, entonces
I ZZ
P (x, y) dx + Q(x, y) dy = (Qx − Py ) dA.
C D
Ahora estamos en condiciones de enunciar y probar uno de los teoremas más importantes
en la teorı́a de funciones de variable compleja.
Demostración. Sea C una curva simple cerrada seccionalmente suave y orientada positivamente
de ecuación,
z = z(t) = x(t) + iy(t)
donde t ∈ [a, b]. Sea además, f (z) = u(x, y) + iv(x, y) una función analı́tica dentro y sobre C
con derivada f ′ continua dentro y sobre C donde al interior de C lo denotamos por D.
Entonces dado que
Además, dado que f ′ es continua, las funciones reales u = u(x, y), v = v(x, y), ux y vx
son continuas dentro y sobre C por lo que aplicando el Teorema de Green y las condiciones de
Cauchy-Riemman se tiene que,
ZZ ZZ
J1 = (−vx − uy ) dA = 0 dA = 0,
D D
y
ZZ ZZ
J2 = (ux − vy ) dA = 0 dA = 0.
D D
Capı́tulo 3 2022 59
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
Observación 3.2.3. Nótese que una vez establecido que el valor de la integral es cero, la
orientación de C es irrelevante. Es decir, el Teorema 3.2.3 es válido también si se toma a C
con orientación negativa pues
I I
f (z) dz = − f (z) dz = 0
C −C
Una consecuencia inmediata del Teorema de Cauchy-Goursat y del Teorema 3.2.2. es el siguiente
resultado:
para cualquier curva C simple cerrada y seccionalmente suave contenida en Ω. Por el Teorema
3.2.2 de equivalencias se tiene que (3) ⇒ (1) de lo que resulta que f tiene una primitiva en
Ω.
Ejemplo 3.2.4. Dado que f (z) = sen z es entera, entonces tiene primitiva en C. Luego,
F (z) = − cos z es una primitiva de f en C. ▲
Capı́tulo 3 2022 60
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
Teorema 3.2.5. Sean C1 una curva simple cerrada seccionalmente suave orientada
positivamente y C2 otra curva simple cerrada seccionalmente suave orientada negativa-
mente interior a C1 . Si f es una función analı́tica en la región limitada entre y sobre
ambas curvas, entonces I I
f (z) dz + f (z) dz = 0.
C1 C2
Demostración. Sea C la curva cerrada simple seccionalmente suave constituida por C1 , seguida
de Γ1 , luego de C2 y por último de Γ2 , como se observa en la Figura 3.6.
Dado que f es analı́tica sobreI C y en la región delimitada por C, entonces por el Teorema
de Cauchy-Goursat se tiene que f (z) dz = 0. Pero además,
I I C
I I I
f (z) dz = f (z) dz + f (z) dz + f (z) dz + f (z) dz = 0. (3.1)
C C1 Γ1 Γ2 C2
Capı́tulo 3 2022 61
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
Capı́tulo 3 2022 62
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
o equivalentemente I I I
f (z) dz = − f (z) dz = f (z) dz.
C1 −C2 C2
Capı́tulo 3 2022 63
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
para cualquier curva simple cerrada orientada positivamente C que contiene en su interior al
punto z0 . Queda como ejercico para el lector probar este resultado.
Teorema 3.2.6 (Fórmula integral de Cauchy). Sea f una función analı́tica dentro y
sobre una curva C simple cerrada seccionalmente suave orientada positivamente. Si z0
es un punto interior a C entonces,
I
f (z)
dz = 2πif (z0 ) . (3.2)
C z − z0
Demostración. Dado que f es analı́tica en z0 entonces resulta continua en dicho punto, por lo
que dado ϵ > 0, existe un δ > 0 tal que
Capı́tulo 3 2022 64
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
es decir, I
f (z)
dz − 2πif (z0 ) = 0,
C z − z0
como querı́amos demostrar.
Capı́tulo 3 2022 65
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
I
ez
Ejemplo 3.2.7. Calculemos dz siendo C la circunferencia de ecuación |z| = 4 reco-
C z − πi
rrida en sentido antihorario.
La función definida por f (z) = ez es entera y además z0 = πi es un punto interior a C,
entonces se tiene que:
I
ez
dz = 2πif (πi) = 2πieπi = −2πi.▲
C z − πi
I
z
Ejemplo 3.2.8. Para obtener el valor de la integral 2
dz siendo C la elipse de ecuación
C 4+z
y2
cartesiana x2 + = 1 recorrida en sentido positivo, consideremos la función dada por
9
z z
f (z) = = .
4+z 2 (z + 2i)(z − 2i)
Como f es analı́tica en C − {2i, −2i}, consideremos las curvas simples cerradas seccional-
mente suaves recorridas positivamente C1 y C2 tal que C = C1 ∪C2 , como se indica en la Figura
3.11.
Luego, I I I
z z
z z+2i z−2i
dz = dz + dz = J1 + J2 .
C1 z − 2i
2
C 4+z C2 z + 2i
z
Dado que f1 (z) = es analı́tica dentro y sobre C1 , entonces se tiene que
z + 2i
J1 = 2πif1 (2i) = πi.
Capı́tulo 3 2022 66
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
z
De igual forma, llamando f2 (z) = se tiene que es analı́tica dentro y sobre C2 , por lo
z − 2i
que obtenemos que
J2 = 2πif2 (−2i) = πi.
Por lo tanto, I
z
dz = 2πi.▲
C 4 + z2
Estableceremos ahora un resultado fundamental que indica que si una función es analı́tica
en un punto, admite derivada de todos los órdenes en ese punto, es decir, es infinitamente
derivable en dicho punto, y todas las derivadas son analı́ticas en él. Tal resultado se conoce
como la Fórmula integral de Cauchy para derivadas.
Observación 3.2.6. Notemos que si en la fórmula integral de Cauchy para derivadas (3.8) se
toma n = 0 se tiene la fórmula integral de Cauchy (3.2).
I
zez
Ejemplo 3.2.9. Para calcular dz siendo C cualquier curva simple cerrada orien-
C (z − 1)
2
tada positivamente y seccionalemente suave que contiene en su interior al punto z0 = 1, consi-
deremos f (z) = zez la cual es entera, entonces
I
zez 2πi ′
dz = f (1) = 4πei.▲
C (z − 1)
2 1!
Capı́tulo 3 2022 67
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
Ejemplo 3.2.10. Obtengamos el máximo del módulo de f (z) = sen z sobre la región
R = {z = x + iy : 0 ≤ x ≤ π, 0 ≤ y ≤ 1} .
Como f es una función entera y no constante se sigue del teorema del módulo máximo que
siendo C la curva simple cerrada seccionalmente suave frontera de la región R como se observa
en la Figura 3.12.
Además si escribimos,
Capı́tulo 3 2022 68
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
Calculemos primero el máximo valor que asume |f (z)| para los z sobre la curva C1 . Para
ello, consideremos la parametrización,
x(t) = t
C1 : , 0 ≤ t ≤ π.
y(t) = 0
Entonces, √
|f (z)| = sen 2 t = |sen t| = sen t con 0 ≤ t ≤ π
π
asume su máximo en t = y es M1 = 1 por lo que |f (z)| toma su valor máximo sobre C1 en
2
π
z = y es M1 = 1.
2
Busquemos ahora el máximo de |f (z)| para los z sobre la curva C2 . Sea la parametrización,
x(t) = π
C2 : , 0 ≤ t ≤ 1.
y(t) = t
Entonces, p
|f (z)| = senh2 t = |senht| = senht con 0 ≤ t ≤ 1
e − e−1
asume su máximo en t = 1 y es M2 = senh1 = por lo que |f (z)| toma su valor máximo
2
e − e−1
sobre C2 en z = π + i y es M2 = .
2
Hallemos ahora el máximo valor que asume |f (z)| para los z sobre la curva C3 . Considere-
mos la parametrización,
x(t) = t
−C3 : , 0 ≤ t ≤ π.
y(t) = 1
Entonces, p
|f (z)| = sen 2 t + senh2 1 con 0 ≤ t ≤ π.
Consideremos la función,
π p p q
Como g ′ (t) = 0 si y sólo si t = ; g(0) = senh1, g(π) = senh1 y g π2 =
p p 2
π p
2
senh 1 + 1, entonces g(t) asume su máximo en t = y es M3 = senh2 1 + 1 por lo
π
2 p
que |f (z)| toma su valor máximo sobre C3 en z = + i y es M3 = senh2 1 + 1.
2
Por último, busquemos el máximo de |f (z)| para los z sobre la curva C4 la cual la podemos
parametrizar por,
x(t) = 0
−C4 : , 0 ≤ t ≤ 1.
y(t) = t
Entonces, p
|f (z)| = senh2 t = |senht| = senht con 0 ≤ t ≤ 1
Capı́tulo 3 2022 69
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
e − e−1
asume su máximo en t = 1 y es M4 = senh1 = por lo que |f (z)| toma su valor máximo
2
−1
e−e
sobre C4 en z = i y es M4 = .
2
Por lo tanto,
p
max |f (z)| = max |f (z)| = max {M1 , M2 , M3 , M4 } = senh2 1 + 1
z∈R z∈C
π
y se obtiene en el punto z = + i. ▲
2
Demostración. De la fórmula integral de Cauchy para las derivadas (3.8) se tiene que,
I I
(n) |f (z)|
f (z0 ) = n! f (z) n!
2πi n+1 dz ≤ |dz| ≤
C (z − z0 ) 2π C |z − z0 |n+1
I
n! M n! M n! M
≤ n+1
|dz| = n+1
2πr = n
2π r C 2π r r
para todo n ∈ N0 .
Como consecuencia de la desigualdad de Cauchy, se tiene el siguiente corolario:
Capı́tulo 3 2022 70
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
4. Utilizando las consecuencias del Teorema de Cauchy, evalúe las siguientes integrales:
Z
π
(a) sen z dz, siendo C la poligonal abierta cuyo origen es el punto z1 = , vértices
C 2
π π π
z2 = + i, z3 = − + i, y extremo z4 = − .
2 2 2
Z
(b) z 3 − iz 2 − 5z + 2i dz, siendo C el arco de elipse de ecuación |z − 3i|+|z + 3i| =
C
10, Re(z) ≥ 0, recorrido en el sentido en el cual Im(z) es creciente.
Z πi
(c) sen (3z) cos(3z) dz
0
Z π+i
(d) (z + 2) eiz dz
0
Capı́tulo 3 2022 71
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
I
ez
(a) dz, siendo C la circunferencia de centro 2 y radio 1, recorrida en sentido
C z −2
antihorario.
I
1
(b) dz, siendo C el rectángulo de vértices z1 = 4 + i, z2 = −4 + i, z3 = −4 − i
C z − 16
4
y z4 = 4 − i, recorrida en sentido horario.
I
5z
(c) dz, siendo C la curva de ecuación |z| = 5, recorrida en
C (z + 1)(z − 2)(z + 4i)
sentido antihorario.
1
6. Sea la función definida por g(z) = . Determine los valores de las integrales de la
z2
+9
función g para todas las posibles curvas simples cerradas seccionalmente suaves que no
contengan puntos de singularidad de g.
Capı́tulo 3 2022 72
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
y represéntelos gráficamente.
2. El diagrama de Smith es un método gráfico que se usa en ingenierı́a eléctrica para analizar
lı́neas de transmisión de alta frecuencia. Establece la relación entre dos variables com-
plejas: la impedancia normalizada z = r + ix, (r ≥ 0, x ∈ R) y el coeficiente de reflexión
Γ = a + ib, (a, b ∈ R).
La transformación
z−1
Γ(z) = ,
z+1
se aplica a una malla formada por rectas verticales infinitas y rectas horizontales semiin-
finitas en el semiplano derecho. La imagen de esta malla en el plano Γ es el diagrama de
Smith.
(a) Pruebe que el diagrama de Smith cuando se aplica la malla Re(z) ≥ 0 bajo la
transformación dada es |Γ| ≤ 1. Grafique ambas regiones e indique si se tratan de
conjuntos abiertos, cerrados, conexos, simplemente conexos y acotados.
(b) Calcule, si es posible, Γ (3e−πi ).
1
(c) Demuestre que z (Γ) = , es decir, que los valores de z correspondientes a
z (−Γ)
valores de Γ que se encuentran en posiciones diametralmente opuestas respecto al
origen son recı́procos uno del otro.
Capı́tulo 3 2022 73
CAPÍTULO 3. INTEGRACIÓN COMPLEJA
(a) Halle la función entera f = f (z) sabiendo que Re {f (z)} = u(x, y) y f (0) = 1.
(b) Pruebe que la función f obtenida se puede expresar como f (z) = e2iz .
(c) Halle todas las soluciones de la ecuación f (z) = i.
(d) Obtenga las trayectorias ortogonales de la curva de ecuación u(x, y) = c donde c una
constante real.
y2
siendo C la semi-elipse en el plano xy de ecuación cartesiana x + = 1 con x ≥ 0,
2
4
recorrida en sentido positivo.
6. Calcule I
2
4z − 2
3
f (z) dz
C z (z + i)
siendo C la curva de ecuación z = 2eiθ , 0 ≤ θ ≤ 2π orientada positivamente y f es una
función analı́tica en |z| < 3 tal que f ′ (0) = f (−i) = 1 y f (0) = i.
Capı́tulo 3 2022 74
Capı́tulo 4
X
∞
f (x) = an xn = a0 + a1 x + a2 x2 + .... + ar xr + ...
n=0
Las series de potencias también son de gran importancia en el comportamiento de las fun-
ciones complejas. En esta sección consideramos algunas de las propiedades de la expansión en
serie de potencias de funciones complejas haciendo, siempre que sea posible, una analogı́a con
la expansión en serie de potencias de las funciones reales correspondientes.
X
∞
an (z − z0 )n = a0 + a1 (z − z0 ) + a2 (z − z0 )2 + .... + an (z − z0 )n + ... (4.1)
n=0
Capı́tulo 4 2022 75
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
llamada serie geométrica de razón z y primer término a. Por el criterio del cociente o de
D’Alembert vemos que la serie geométrica es convergente si
n+1
az
lim = lim |z| = |z| < 1.
n→+∞ az n n→+∞
por lo que el radio de convergencia es R = 1 y el cı́rculo de convergencia es |z| < 1. Además,
la serie geométrica es divergente si |z| ≥ 1.
Teorema 4.1.1 (Series de Taylor). Sea f una función analı́tica en el disco abierto
|z − z0 | < R. Entonces la función f admite una representación en serie de potencias,
llamada serie de Taylor de f centrada en z0 ,
X
∞
f (z) = an (z − z0 )n
n=0
f (n) (z0 )
si |z − z0 | < R donde an = para n = 0, 1, 2, 3....
n!
Demostración. Por una cuestión de simplicidad, presentamos la prueba del teorema sólo para
el caso especial donde z0 = 0. Dado que z0 es arbitrario, la generalización del mismo es
consecuencia inmediata de considerar z0 = 0.
Sea f una función analı́tica en el disco abierto |z| < R. Consideremos una circunferencia C
orientada positivamente de ecuación |z| = r tal que r < R de modo tal que esté ı́ntegramente
contenida en el disco |z| < R como se observa en la Figura 4.1.
Entonces por las fórmulas integrales de Cauchy se tiene que para cualquier z interior a C,
I
1 f (w)
f (z) = dw, (4.2)
2πi C w − z
Capı́tulo 4 2022 76
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
y I
(n) n! f (w)
f (z) = dw (4.3)
2πi C (w − z)n+1
para todo n = 0, 1, 2, ...
Usando el hecho de que para cualquier número complejo a 6= 1 se verifica la identidad
1 am
= 1 + a + a2 + ... + am−1 + , m = 1, 2, 3, ... (4.4)
1−a 1−a
entonces
z n−1 (z/w)m
1 1 1 1 z z 2
= = 1+ + + ... + + , (4.5)
w−z w1− z
w
w w w w 1 − wz
o equivalentemente
1 1 1 1 1 1
= + 2 z + 3 z 2 + ... + m z m−1 + zm. (4.6)
w−z w w w w (w − z)wm
I
1 1 1 1 2 1 m−1 1 m
f (z) = + z + 3 z + ... + m z + z f (w) dw. (4.7)
2πi C w w2 w w (w − z)wm
Supongamos que |z| = r0 . Dado que z es un punto interior de la curva C entonces r0 < r,
y además como w ∈ C entonces
r0m M M r r0 m
|φm (z)| ≤ 2πr = . (4.12)
2π (r − r0 )rm r − r0 r
Capı́tulo 4 2022 77
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
Figura 4.1: Disco abierto |z| < R y circunferencia C : |z| = r postivamente orientada
r0
Dado que < 1 entonces φm (z) → 0 si m → +∞, por lo que de (4.8) se sigue
r
f ′ (0) f ′′ (0) 2 f (m−1) (0) m−1
f (z) = f (0) + z+ z + ... + z + .... (4.13)
1! 2! (m − 1)!
en el disco abierto |z| < R.
Capı́tulo 4 2022 78
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
1
4. Sea la función dada por f (z) = . El desarrollo en serie de Maclaurin será válido
1−z
para |z| < 1. Además, las derivadas de f vienen dadas por
n!
f (n) (z) =
(1 − z)n+1
para todo n = 0, 1, 2, 3, ... por lo que se tiene que
El teorema anterior nos permite obtener el desarrollo en serie de Taylor de una función de
un modo más eficiente, haciendo uso de los desarrollos por serie de Taylor ya conocidos para
algunas funciones elementales. Veamos algunos ejemplos:
f (z) = z 2 e3z
X
∞
3n z n+2
f (z) = z e 2 3z
= si |z| < +∞.
n=0
n!
Capı́tulo 4 2022 79
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
X∞
(−1)n
sen (iz) = (iz)2n+1 si |iz| < +∞
n=0
(2n + 1)!
X
∞
z 2n+1
senh(z) = si |z| < +∞.
n=0
(2n + 1)!
cosh z = cos(iz)
X∞
z 2n
cosh z = si |z| < +∞.
n=0
(2n)!
1 1 1 z
Para ello, notemos que = , por lo que sustituyendo z por − en (4.17)
2+z 2 1 − (−z/2) 2
se tiene entonces que
n
1X
∞
1 −z −z
= si <1
2+z 2 n=0 2 2
o equivalentemente,
1 X (−1)n z n
∞
= si |z| < 2.
2+z n=0
2n+1
Capı́tulo 4 2022 80
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
1 X
∞
= (−1)n (z − 1)n
z n=0
Teorema 4.1.4. Si la serie de potencias (4.18) converge a f (z), entonces puede in-
tegrarse término a término en toda curva simple seccionalmente suave C interior al
cı́rculo de convergencia de esa serie.
Además,
Z Z X
∞ X∞
an
f (z) dz = an (z − z0 )n dz = (z − z0 )n+1 .
C C n=0 n=0
n + 1
−1 X∞
=− wn si |w| < 1,
1−w n=0
Capı́tulo 4 2022 81
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
Observación 4.1.2. Si la curva C del Teorema 4.1.4 es además cerrada entonces tenemos que
I X
∞
an (z − z0 )n dz = 0.
C n=0
Teorema 4.1.5. Si la serie de potencias (4.18) converge a f (z), entonces puede de-
rivarse término a término en todo punto interior al cı́rculo de convergencia de dicha
serie. Además,
!
d X∞ X∞
′
f (z) = an (z − z0 ) n
= nan (z − z0 )n−1 .
dz n=0 n=1
z
Ejemplo 4.1.3. Por un lado, la función dada por f (z) = puede expresarse
(1 + z)2
d −1
f (z) = z ,
dz 1+z
y dado que
−1 X∞ X∞
=− n
(−z) = (−1)n+1 z n
1+z n=0 n=0
X
∞
f (z) = (−1)n+1 nz n si |z| < 1.▲
n=1
Capı́tulo 4 2022 82
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
f (z) = e1/z .
1
Si sustituı́mos z por en
z
X
∞
zn
z
e = si |z| < +∞
n=0
n!
obtenemos el desarrollo en serie de potencias
X∞
1 1
e 1/z
= si < +∞
z n n! z
n=0
y en definitiva
1 1 1
e1/z = 1 + + 2 + 3 + ... si 0 < |z| < +∞.
z 2z 6z
Claramente, la función f no es analı́tica en z0 = 0, por lo que no podemos hallar una serie
de Maclaurin en dicho punto. No obstante, obtuvimos una representación en serie de potencias
para f que contiene potencias negativas de z. Esto es posible, gracias al teorema de Laurent.
Teorema 4.1.6 (Serie de Laurent). Sean w = f (z) una función analı́tica en una
corona abierta R1 < |z − z0 | < R2 y C una curva simple cerrada seccionalmente suave
positivamente orientada en torno a z0 y contenida en dicha corona circular abierta
como se observa en la Figura 4.2.
Entonces, la función f admite la representación en serie de potencias, llamada serie
de Laurent siguiente
X
∞ X
∞
bn
f (z) = an (z − z0 ) + n
n=0 n=1
(z − z0 )n
Capı́tulo 4 2022 83
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
X
∞ X
∞
bn
f (z) = an (z − z0 ) +
n
si R1 < |z − z0 | < R2
n=0 n=1
(z − z0 )n
X
∞
f (z) = cn (z − z0 )n
n=−∞
X
∞ X
∞
bn
f (z) = an (z − z0 ) +
n
si R1 < |z − z0 | < R2
n=0 n=1
(z − z0 )n
Capı́tulo 4 2022 84
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
f (n) (z0 )
an = para todo n = 0, 1, 2, 3, ...
n!
Por lo tanto, el desarrollo en serie de Laurent para f se reduce al desarrollo en serie de
Taylor de f centrada en z0 .
De manera similar, si f es analı́tica en todo punto del plano complejo que verifique R1 <
|z − z0 |, el radio R2 se puede como R2 = +∞, por lo que la condición de validez de la
serie es
R1 < |z − z0 | < +∞.
converge a f (z) en todo punto interior de una corona R1 < |z − z0 | < R2 , entonces
I
1 f (z)
cn = dz, n ∈ Z
2πi C (z − z0 )n+1
A partir del teorema anterior, que establece la unicidad de representación en serie de Laurent
de una función, los coeficientes de una serie de Laurent para una función f no se suelan hallar
recurriendo al uso directo de sus representaciones integrales, sino por otros métodos de manera
similar a lo realizado para el desarrollo por series de Taylor.
Capı́tulo 4 2022 85
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
1 X ∞
= wn si |w| < 1. (4.19)
1 − w n=0
1 1 X∞
=− =− z n para |z| < 1 (4.20)
z−1 1−z n=0
z 1
y por otro lado como |z| < 1 entonces < < 1 entonces usando (4.19) se tiene que:
2 2
1 X z n
∞
1 1
− = = . (4.21)
z−2 2 (1 − z/2) 2 n=0 2
Capı́tulo 4 2022 86
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
1
Ejemplo 4.2.1. La función dada por f (z) = tiene un punto singular en z0 = −2i.
z + 2i
Dicho punto singular es una singularidad aislada dado que la función f resulta analı́tica en
cualquier entorno reducido de z0 = −2i. ▲
Capı́tulo 4 2022 87
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
Ejemplo 4.2.2. Consideremos la función logaritmo principal f (z) = log z. Esta función tiene
claramente una singularidad en z0 = 0.
Sin embargo, no se trata de una singularidad aislada de f pues no existe ningún entorno
reducido de z0 = 0 que no contenga puntos del semieje real negativo donde la función es no
analı́tica. ▲
Capı́tulo 4 2022 88
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
1 X (−1)n 2n+1
∞
sen z z2 z4
f (z) = = z =1− + + ...
z z n=0 (2n + 1)! 6 120
si 0 < |z| < +∞, entonces la función f tiene una singularidad evitable en z0 = 0. ▲
Ejemplo 4.2.6. La función dada por (4.27) tiene una singularidad evitable en z0 = 0 dado que
ez sen z
lim = lim ez (sen z + cos z) = 1.▲
z→0 z z→0
Capı́tulo 4 2022 89
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
z − zk −z 2 1
lim = lim 1/z = − 2 2
z→zk e1/z − 1 z→zk e 4k π
si k ∈ Z − {0}. Por lo tanto, las singularidades aisladas zk con k ∈ Z − {0} son todos polos
simples o de primer orden de f . ▲
Singularidad en el infinito
Si trabajamos en el plano complejo extendido C∗ = C ∪ {∞}, podemos estudiar el comporta-
miento de una función f (z) en el punto z = ∞ y determinar qué tipo de singularidad presenta
la función f en dicho punto.
1
Para ello, consideremos la transformación z = . El comportamiento de f (z) en z = ∞
w
corresponde al comportamiento de g(w) = f (1/w) en w = 0, por lo que si g es analı́tica en
w = 0 entonces la función f es analı́tica en z = ∞ y si la función g tiene una singularidad en
w = 0 entonces la función f tiene una singularidad en z = ∞.
Además, si la función g tiene un polo o una singularidad evitable o una singularidad esencial
en w = 0 entonces la función f tiene un polo o una singularidad evitable o una singularidad
esencial respectivamente en z = ∞.
Por ejemplo, la función dada por
f (z) = z 2
Capı́tulo 4 2022 90
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
1
tiene un polo de orden 2 en z = ∞ ya que la función definida por g(w) = 2 tiene un polo de
w
segundo orden en w = 0.
Si ahora consideramos por ejemplo la función definida por
f (z) = ez
vemos que posee una singularidad esencial en z = ∞. En efecto, basta con obtener el desarrollo
en serie de Laurent de la función dada por g(w) = e1/w en 0 < |w| < +∞ para observar que g
tiene una singularidad esencial en w = 0.
Res(f, z0 ) = b1 .
Capı́tulo 4 2022 91
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
Res(f, 2) = 3.▲
1 X z 2n+1
∞
senhz 1 1 z
f (z) = 4
= 4 = 3+ + + ...
z z n=0 (2n + 1)! z 6z 120
f (z) = e−1/z .
Res(f, 0) = −1.▲
Observación 4.2.1. Si una función w = f (z) tiene una singularidad evitable en z0 entonces
el desarrollo en serie de Laurent (4.29) de f en la corona 0 < |z − z0 | < R no contiene ningún
término en la parte principal, es decir
bn = 0, n = 1, 2, 3, ...
Res(f, z0 ) = 0.
Capı́tulo 4 2022 92
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
Cuando una función w = f (z) tiene una singularidad esencial en z0 , la única manera en que
podemos calcular el residuo en dicha singularidad aislada es obteniendo el desarrollo en serie
1
de Laurent de f en la corona 0 < |z − z0 | < R y eligiendo el coeficiente b1 de .
z − z0
Sin embargo, si la función f tiene un polo en z0 existen otros métodos alternativos para
calcular el residuo sin recurrir a su desarrollo en serie de Laurent. A continuación, enunciaremos
dos de ellos:
ϕ(z)
f (z) =
(z − z0 )m
ϕ(m−1) (z0 )
Res(f, z0 ) = .
(m − 1)!
Veamos con un ejemplo que el Teorema 4.2.2 no es aplicable en cualquier caso que la función
presente una singularidad aislada de tipo polo. Consideremos la función definida por
3
f (z) = .
ez − 1
Capı́tulo 4 2022 93
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
Teorema 4.2.3 (Cálculo del residuo en un polo). Sea f una función que admite la
representación
p(z)
f (z) =
q(z)
donde p y q son funciones analı́ticas en z0 y p(z0 ) 6= 0.
Entonces f tiene un polo de orden m en z0 si y sólo si q(z0 ) = q ′ (z0 ) = ... = q (m−1) (z0 ) =
0 y q (m) (z0 ) 6= 0.
Además,
1 dm−1
Res(f, z0 ) = lim m−1 {(z − z0 )m f (z)}.
(m − 1)! z→z0 dz
Ejemplo 4.2.13. Retomemos el ejemplo anterior donde vimos que la función dada por
3
f (z) =
ez − 1
tiene un polo simple en z0 = 0. Llamemos
p(z) = 3
y
q(z) = ez − 1
donde concluimos que p y q son analı́ticas en z0 = 0 y p(0) = 3 6= 0.
Entonces, dado que q(0) = 0 y q ′ (0) = 1 6= 0 podemos concluir (como ya lo sabı́amos) del
Teorema 4.2.3 que z0 = 0 es un polo simple de f .
Además, el residuo de f en z0 = 0 está dado por
1 3z
Res(f, 0) = lim {zf (z)} = lim z = 3.▲
0! z→0 z→0 e − 1
Capı́tulo 4 2022 94
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
2i
Ejemplo 4.2.14. Calculemos el residuo de f (z) = en la singularidad aislada z0 = 0.
z(ez− 1)
Para ello, llamemos p(z) = 2i y q(z) = z(ez − 1). Notemos que p y q son analı́ticas en
z0 = 0 y además p(0) = 2i 6= 0.
Entonces como q(0) = q ′ (0) = 0 pero q ′′ (0) = 2 6= 0 entonces se tiene que f tiene un polo
doble en z0 = 0.
Luego, el residuo de f en z0 = 0 es
1 d 2 2i(ez − 1 − zez )
Res(f, 0) = lim z f (z) = lim = −i.▲
1! z→0 dz z→0 (ez − 1)2
Capı́tulo 4 2022 95
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
Teorema 4.3.1 (El teorema de los residuos). Sea C una curva seccionalmente suave
simple cerrada orientada positivamente, dentro y sobre la cual la función f es analı́tica
excepto en un número finito z1 ,z2 ,...,zn de singularidades aisladas interiores a C. En-
tonces I Xn
f (z) dz = 2πi Res(f, zk ).
C k=1
Demostración. Sean z1 ,z2 ,...,zn singularidades aisladas de f interiores a la curva C simple ce-
rrada seccionalmente suave y orientada positivamente.
Consideremos además, curvas simples cerradas seccionalmente suaves y positivamente orien-
tadas Ck con k = 1, 2, ..., n interiores a C y disjuntas entre sı́ que contienen en su interior a los
puntos zk con k = 1, 2, ..., n respectivamente, como se observa en la Figura 4.3.
Las curvas Ck con k = 1, 2, ..., n junto con la curva C forman la frontera de una región
cerrada dentro y sobre la cual la función f es analı́tica y cuyo interior es un conjunto conexo.
Entonces, por una consecuencia del Teorema de Cauchy-Goursat se tiene que
I n I
X
f (z) dz − f (z) dz = 0. (4.32)
C k=1 Ck
Ahora bien, el residuo en cada una de las singularidades aisladas zk donde k = 1, 2, ..., n de
f está dado por I
1
Res(f, zk ) = f (z) dz (4.33)
2πi Ck
para todo k = 1, 2, ..., n.
Luego, sustituyendo (4.33) en (4.32) se tiene entonces que
I X
n
f (z) dz − 2πi Res(f, zk ) = 0
C k=1
o equivalentemente
I X
n
f (z) dz = 2πi Res(f, zk ).
C k=1
Ejemplo 4.3.1. Haciendo uso del teorema de los residuos obtengamos el valor de la integral
I z
e −1
2
dz
C z +z
Capı́tulo 4 2022 96
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
ez − 1
Las singularidades de f (z) = son z1 = 0 y z2 = −1, ambas aisladas e interiores a la
z2 + z
circunferencia positivamente orientada C de centro z = 0 y radio 4 por lo que
I z
e −1
2
dz = 2πi [Res(f, 0) + Res(f, −1)] .
C z +z
Capı́tulo 4 2022 97
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
cos z ez
(a) f (z) = , z0 = 0 (c) f (z) = , z0 = 1, z1 = 2i
1 + z2 1 + ez
ez 1
(b) f (z) = , z0 = 4i, z1 = 1 + 4i (d) f (z) = , z0 = 1, z1 = −1
z2 − z cos(πz)
4. En cada uno de los siguientes casos obtenga la serie de potencias (de Taylor o Laurent),
que represente a f en las regiones indicadas.
1
(a) f (z) = ; S1 = {z ∈ C : |z| < 1}, S2 = {z ∈ C : |z| > 1}.
z−1
1
(b) f (z) = 2 ; S1 = {z ∈ C : |z| < 1}, S2 = {z ∈ C : 1 < |z| < 3}, S3 =
z + 4z + 3
{z ∈ C : |z| > 3}, S4 = {z ∈ C : 1 < |z + 4| < 3}.
z
(c) f (z) = ; S1 = {z ∈ C : |z − 1| > 1}, S2 = {z ∈ C : 0 < |z − 2| < 1}.
(z − 1)(z − 2)
z
(d) f (z) = ; S1 = {z ∈ C : |z − 1| > 2}.
z2 −1
5. En cada uno de los casos siguientes obtenga todos los posibles desarrollos en serie de
Laurent de la función f alrededor del punto z0 indicado. Determine en cada caso la
región de convergencia.
1
(a) f (z) = , z0 = i
z2 + 1
1
(b) f (z) = z 2 sen z−1 , z0 = 1
1
(c) f (z) = z 2 e z , z0 = 0
7. Obtenga la función suma f (z) de cada una de las siguientes series para |z| < 1:
X
∞ X
∞
zn
n
(a) nz (b)
n=1 n=1
n
8. Una función se dice meromorfa si sus únicas singularidades son polos. Se puede demostrar
que el cociente de dos funciones enteras es una función entera o una función meromorfa.
Analice en cuál de los siguientes casos la función f es entera o meromorfa.
Capı́tulo 4 2022 98
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
p(x)
donde f (x) = es par con p y q polinomios reales sin factores comunes, y donde q no
q(x)
tiene raı́ces reales.
Capı́tulo 4 2022 99
CAPÍTULO 4. SERIES COMPLEJAS Y RESIDUOS
donde C es el contorno cerrado ilustrado en la Figura 4.4, que consta del eje real entre
−R y R y la semicircunferencia Γ, de radio R, en el semiplano superior z.
es decir, grado(q) ≥ 2 + grado(p). Ası́, suponiendo que esta condición es satisfecha, esta
aproximación puede utlizarse para calcular
Z +∞
f (x) dx.
−∞
+∞ Z
dx π
(a) Pruebe que 2 2
= .
−∞ (x + 4) 16
Z +∞
2x − 1
2
(b) Calcule dx.
0 x4 + 5x2 + 4
Sugerencia: Use el hecho de que si f es una función par entonces se tiene que
Z +∞ Z +∞
f (x) dx = 2 f (x) dx.
−∞ 0
Z 2π
√
dθ 2 3
(a) Pruebe que = π.
0 2 + sen θ 3
Z 2π
dθ
(b) Calcule .
0 3 − 2 cos θ + sen θ
La transformación de Laplace
Si queremos obtener la transformada de Laplace de una función f = f (t) la cual está definida
también para valores de t < 0, haremos uso de la función escalón unitario o función de
Heaviside
0 si t < 0
H(t) = .
1 si t ≥ 0
De esta manera,
0 si t < 0
f (t)H(t) =
f (t) si t ≥ 0
es una función causal y podemos ası́ aplicar (5.1) para obtener su transformada de Lapla-
ce. Sin embargo, la transformada de Laplace de f (t)H(t) contiene sólo información sobre el
comportamiento de f (t) para t ≥ 0.
La transformada de Laplace de t
Sea f (t) = t, t ≥ 0. Es claro que:
Z T2
T si s = 0
te−st dt = 2 (5.4)
0
1
+ e −sT
−
T
−
1
si s 6= 0
s2 s s2
Si en (5.5) hacemos T → +∞, la integral diverge si Re(s) ≤ 0 mientras que si Re(s) > 0 la
integral converge resultando que:
s
L {cos(ωt)} = si Re(s) > 0.
s2 + ω2
|f (t)| ≤ M eτ t
para toda t ≥ 0.
2
Ejemplo 5.1.2. La función dada por f (t) = et no es de orden exponencial en [0, +∞).
En efecto,
2
et
lim τ t = lim et −τ t = +∞
2
t→+∞ e t→+∞
para cualquier t ≥ 0. ▲
Observación 5.1.1. El teorema establece condiciones suficientes para asegurar que f es trans-
formable por Laplace, pero no necesarias. En efecto, la función dada por
1
f (t) = √
t
tiene una discontinuidad infinita en t = 0, sin embargo es transformable por Laplace pues si
consideramos el cambio de variable
x2
t=
s
se tiene que
Z +∞ Z +∞ √ r
1 1 2 2 π π
e−st √ dt = √ e−x dx = √
2
L √ = =
t 0 t s 0 s s s
si Re(s) > 0.
El Teorema 5.1.1 indica además de condiciones suficientes para la existencia de transfor-
mada de Laplace de una función f = f (t), que su función transformada
Z +∞
F (s) = L {f (t)} = e−st f (t) dt
0
L : FR → FC
tal que Z +∞
L {f (t)} = F (s) = e−st f (t) dt
0
donde para cualesquiera α y β contantes y f y g funciones de R se tiene que:
Teorema 5.2.1 (Primer teorema de traslación). Sea f = f (t) una función cuya trans-
formada de Laplace es F = F (s) con Re(s) > τ . Entonces la función dada por eat f (t)
es transformable por Laplace y su transformada está dada por
L eat f (t) = F (s − a) = L {f (t)}|s−a
1
Ejemplo 5.2.1. L {te−2t } = L {t}|s+2 = si Re(s) > 0 + (−2) = −2.▲
(s + 2)2
2
Ejemplo 5.2.2. L {e3t sen (2t)} = L {sen (2t)}|s−3 = si Re(s) > 0 + 3 = 3.▲
(s − 3)2 + 4
Z Z +∞ n
(n) dn +∞
−st ∂ −st
F (s) = n e f (t) dt = e f (t) dt =
ds 0 0 ∂sn
Z +∞
= (−1) n
tn e−st f (t) dt = (−1)n L {tn f (t)}
0
si Re(s) > τ .
d d 3 6s
Ejemplo 5.2.3. L {tsen (3t)} = (−1) (L {sen (3t)}) = −
1
2
= . ▲
ds ds s +9 (s2 + 9)2
Ejemplo 5.2.4. Si queremos obtener L {tn } donde n ∈ N, notemos que:
d d 1 1
n = 1 ⇒ L {t} = − (L {1}) = − = 2,
ds ds s s
d d 1 2
n = 2 ⇒ L t2 = − (L {t}) = − 2
= 3,
ds ds s s
d d 2 6
n = 3 ⇒ L t3 = − L t2 = − = ,
ds ds s3 s4
y ası́ sucesivamente se puede probar que para cualquier n ∈ N,
n!
L {tn } = .▲
sn+1
f (t) = g(t)
para todo t ≥ 0.
Por lo tanto, esta solución y(t) de (5.7) es llamada antitransformada de Laplace o trans-
formada inversa de Laplace de F , se simboliza por L−1 y se caracteriza por
Si las funciones f y g del Teorema 5.3.1 son seccionalmente continuas en [0, +∞), el teorema
de Lerch establece que si sus transformadas son iguales entonces las funciones f y g son “casi”
iguales. Sin embargo, no identificaremos esta situación y trabajaremos como si fuesen continuas
pues su desarrollo no es sencillo y escapa de los intereses de este curso.
También, tenemos que tener en cuenta que si la función f = f (t) no es causal, deberemos
escribir
L−1 {F (s)} = f (t)H(t).
Es sencillo verificar que el operador antitransformada de Laplace L−1 es lineal, es decir si
F y G son funciones cuyas antitransformadas existen, entonces
−1 1 1 1 1 1 −1 −3t 1 2t
L = − L−1 + L−1 = e + e .▲
(s + 3)(s − 2) 5 s+3 5 s−2 5 5
F (s) = L {f (t)} .
Entonces, Z +∞
′
L {f (t)} = e−st f ′ (t) dt (5.8)
0
pues lim e−st f (t) = 0 debido a que la función dada por e−st f (t) es de orden exponencial.
t→+∞
y ası́, se puede probar, mediante inducción matemática, que para cualquier n ∈ N, la transfor-
mada de Laplace de la derivada n−ésima de f está dada por:
L f (n) (t) = −f (n−1) (0) − sf (n−2) (0) − ... − sn−1 f (0) + sn F (s) (5.14)
con las condiciones iniciales x(0) = 1 y x′ (0) = 0, apliquemos a ambos miembros de (5.15) la
transformada de Laplace
L {x′′ + 5x′ + 6x} = L 2e−t . (5.16)
Por las propiedades de la transformada de Laplace y usando (5.13) se tiene entonces que
2
s2 X(s) − sx(0) − x′ (0) + 5 (sX(s) − x(0)) + 6X(s) = (5.17)
s+1
donde X(s) = L {x(t)}. Usando las condiciones iniciales, entonces (5.17) se expresa como
2
s2 + 5s + 6 X(s) = +s+5 (5.18)
s+1
o equivalentemente
1 1 1
X(s) = + − . (5.19)
s+1 s+2 s+3
Entonces, antitransformando (5.19) se tiene la solución del problema dado:
Ejemplo 5.4.2. Supongamos que queremos hallar las funciones x = x(t) e y = y(t) que
satisfacen el siguiente sistema de ecuaciones diferenciales
′
x + y ′ + 5x + 3y = e−t
2x′ + y ′ + x + y = 3
Ahora bien, para resolver el sistema (5.20), podemos utilizar cualquiera de los métodos
resolución de sistemas de ecuaciones. Aquı́, usaremos la regla de Cramer. Para ello, calculemos
primero el determinante de la matriz de coeficientes de (5.20),
s+5 s+3
4 = = (s + 5)(s + 1) − (s + 3)(2s + 1) = −(s + 2)(s − 1)
(5.21)
2s + 1 s + 1
Entonces, para obtener X(s) usamos (5.21) y (5.20) donde se tiene que
4X 2s2 + 14s + 9 9 1 11 1 25 1
X(s) = = =− − + (5.23)
4 s(s + 2)(s − 1) 2 s 16 s + 2 3 s−1
L {g ′ (t)} = L {f (t)}
o equivalentemente
sL {g(t)} − g(0) = L {f (t)} . (5.28)
Como g(0) = 0 entonces (5.28) se reduce a
1
L {g(t)} = L {f (t)}
s
o equivalentemente
Z t
1
L f (x) dx = L {f (t)} . (5.29)
0 s
Z
t 1 6 2
Ejemplo 5.4.3. L 3
x + sen (2x) dx = L {t3 + sen (2t)} = 5 + 2
.▲
0 s s s(s + 4)
También, la función escalón unitario es especialmente útil para dar expresiones compac-
tas de ciertas funciones que tienen una ley definida por tramos. Sea por ejemplo,
f1 (t) si 0 < t < t1
f (t) = f2 (t) si t1 < t < t2 (5.30)
f3 (t) si t > t2
cuya gráfica se muestra en la Figura 5.4.
Para obtener una expresión compacta de (5.30), podemos hacer uso de las “operaciones de
encendido”
1. encender f1 (t) en t = 0,
f (t) = f1 (t)H(t) + [f2 (t) − f1 (t)] H(t − t1 ) + [f3 (t) − f2 (t)] H(t − t2 ). (5.31)
que se conoce como la construcción compacta de (5.30) haciendo uso de la función “sombrero
de copa” y cuya expresión es,
1 si a < t < b
H(t − a) − H(t − b) =
0 si t < a ó t > b
e−as
L {H(t − a)} = (5.34)
s
si Re(s) > 0.
En particular, y como es de esperar, si a = 0 entonces (5.34) se reduce a:
1
L {H(t)} = (5.35)
s
si Re(s) > 0.
Teorema 5.5.1 (Segundo teorema de traslación). Sean f = f (t) una función cuya
transformada de Laplace es F = F (s) y a > 0. Entonces
Z +∞ Z +∞ Z +∞
−st −s(τ +a) −as
f (t − a)e dt = f (τ )e dτ = e f (τ )e−sτ dτ = e−as L {f (τ )} (5.37)
a 0 0
Ejemplo 5.5.3. Para calcular la transformada de la Laplace de la función del Ejemplo 5.5.1
definida por 2
2t si 0 ≤ t < 3
f (t) = t+4 si 3 ≤ t < 5
9 si t ≥ 5
primero expresamos la misma en una forma compacta haciendo uso de la función escalón de
Heaviside. Del Ejemplo 5.5.1 se tiene que:
f (t) = 2t2 H(t) + t + 4 − 2t2 H(t − 3) + (5 − t) H(t − 5)
ó equivalentemente
f (t) = 2t2 H(t) + t + 4 − 2t2 H(t − 3) − (t − 5) H(t − 5). (5.38)
Ahora bien, para obtener la transformada de Laplace de f (t), haremos uso del segundo
teorema de traslación donde primero debemos escribir la expresión que multiplica a H(t − 3)
en términos de t − 3. Para ello, sea
z =t−3
luego,
2 −3s
2 −3s e−3s
L {f (t)} = 2 3
− 2e L t − 11e L {t} − 11 − e−5s L {t}
s s
o equivalentemente
4 4 11 11 e−5s
L {f (t)} = 3 − e−3s 3
+ 2 + − .▲
s s s s s2
Dado que
s+3 A Bs + C (A + B)s2 + Cs + A
= + = ⇔ A = 3, B = −3, C = 1
s(s2 + 1) s s2 + 1 s(s2 + 1)
se sigue que,
−1 s+3 −1 3 s 1
L =L −3 2 + 2 = 3 − 3 cos t + sen t (5.40)
s(s2 + 1) s s +1 s +1
es convergente si −st
e = e−xT < 1
entonces x > 0, siendo s = x + iy. Luego se sigue que
X
+∞
1
e−snT = (5.43)
n=0
1 − e−sT
o equivalentemente,
Z T
1
L {f (t)} = f (t) e−st dt si Re(s) > 0.
1 − e−sT 0
Ejemplo 5.6.1. Para obtener la transformada de Laplace de la onda periódica diente de sierra
f (t) = t, t ∈ [0, 1]
periódica de perı́odo T = 1 como se observa en la Figura 5.6
calculamos,
Z −s
1 1
−st 1 e e−s 1
L {f (t)} = te dt = − − 2 + 2 .▲
1 − e−s 0 1 − e−s s s s
Z +∞
δ(t) dt = 1.
−∞
La llamada función impulso unitario es una función idealizada. Ciertamente no puede ser
una función pues es “imposible” que una función sea nula excepto en un único punto, y su
integral sea no nula.
Esta función se fundamenta a través de la Teorı́a de Distribuciones. Puede ser interpre-
tada fı́sicamente como un shock o impulso que transmite energı́a unitaria instantáneamente
(idealización).
Si se desea considerar la función impulso unitario en el instante de tiempo t = a, la deno-
tamos δ(t − a), y está definida como el caso lı́mite del pulso unitario:
0 si t < a − T2
Φ(t) = 1
si a − T2 ≤ t < a + T2
T
0 si t ≥ a + T2
conocida con el nombre de propiedad de filtrado. Además, si α < a < β entonces se tiene
que Z β
f (t)δ(t − a) dt = f (a).
α
Haciendo uso de la propiedad de filtrado se tiene entonces que la transformada de Laplace
de δ = δ(t − a) con a > 0 es
Z +∞
L {δ(t − a)} = δ(t − a)e−st dt = e−as .
0
H ′ (t − a) = δ(t − a).
y en particular tomando a = 0,
L δ (n) (t) = sn .
5.8 Convolución
En Matemática y, en particular, en Análisis Funcional, una convolución es un operador ma-
temático que transforma dos funciones f y g en una tercera función que en cierto sentido
representa la magnitud en la que se superoponen f y una versión trasladada e invertida de g.
La convolución y las operaciones relacionadas se encuentran en diversas aplicaciones de
Ingenierı́a y Matemática. Por ejemplo, en Eléctrica y Electrónica, la convolución permite
determinar la respuesta de un sistema lineal invariante en el tiempo ante cualquier entrada,
a partir del conocimiento de la respuesta del sistema ante una única entrada particular, el
impulso.
También, la convolución, es particularmente útil en la solución de ecuaciones diferencia-
les. La transformada de Laplace de la solución (inicialmente desconocida) de una ecuación
diferencial algunas veces se reconoce como el producto de las transformadas de dos funciones
desconocidas. Por ejemplo, cuando transformamos el problema con condiciones iniciales
x′′ + x = cos t,
Esto nos sugiere que deberı́a existir una forma de combinar las dos funciones sen t y cos t
para obtener la función x(t) cuya transformada sea el producto de las transformadas de aquellas.
Pero, evidentemente, x(t) no es sólo el producto de sen t y cos t porque
1 1 s
L {cos tsen t} = L sen (2t) = 2 6= 2 .
2 s +4 (s + 1)2
Ası́ es que
L {cos tsen t} 6= L {cos t} L {sen t} .
Por tales motivos, estamos ante la necesidad de definir la convolución entre las funciones
f = f (t) y g = g(t), denotada f ∗ g, de tal modo que la transformada de f ∗ g resulte ser el
producto de las transformada de f y g.
f (τ ) = g (τ ) = 0
si τ < 0 y
g (t − τ ) = 0
si τ > t por lo que (5.44) se puede expresar como
Z t
(f ∗ g) (t) = f (τ ) g (t − τ ) dτ. (5.45)
0
Demostración. Ejercicio.
El siguiente teorema nos permite hallar la transformada de Laplace de la convolución de
dos funciones:
Demostración. Sean Z +∞
F (s) = L {f (x)} = f (x)e−sx dx
0
y Z +∞
G(s) = L {g(y)} = g(y)e−sy dy.
0
Entonces por el teorema de Fubini,
Z +∞ Z +∞ ZZ
−s(x+y)
F (s)G(s) = f (x)g(y)e dx dy = f (x)g(y)e−s(x+y) dx dy (5.46)
0 0 R
t(x, y) = x + y, τ (x, y) = y
donde |J(x, y)| = |tx τy − τx ty | = 1 se tiene que (5.46) se puede expresar como
ZZ
F (s)G(s) = f (t − τ )g(τ )e−st dτ dt (5.47)
D
donde D es la región del plano representada en la Figura 5.10 que se obtiene de la transformación
dada, por lo que (5.47) se puede expresar como
Z +∞ Z t Z +∞
−st
F (s)G(s) = e f (t − τ ) g(τ ) dτ dt = e−st f (t) ∗ g(t) dt
0 0 0
o equivalentemente
F (s)G(s) = L {f (t) ∗ g(t)} .
x′′ + x = cos t,
con las condiciones iniciales x(0) = x′ (0) = 0, planteada en la introducción de esta Sección.
Como
s s 1
X(s) = 2 2
= 2 2
(s + 1) s +1s +1
entonces
−1 −1 s −1 1
x(t) = L {X(s)} = L ∗L = cos t ∗ sen t.
s2 + 1 s2 + 1
Entonces,
Z t
t
x(t) = cos(t − τ )sen (τ ) dτ = sen t.▲
0 2
9. Use la transformada de Laplace para resolver los siguientes sistemas de ecuaciones dife-
renciales:
′′
x −y =0
(a) tal que x(0) = y(0) = 1, x′ (0) = y ′ (0) = −1.
y ′′ − x = 0
′′
x = y − 2x
(b) tal que x(0) = 4, y(0) = 2, x′ (0) = y ′ (0) = 0.
y ′′ = x − 2y
′′
2x − y ′′ − x′ − y ′ = 3y − 9x
(c) tal que x(0) = x′ (0) = 1, y(0) = y ′ (0) = 0.
2x′′ − y ′′ + x′ + y ′ = 5y − 7x
10. Vibraciones mecánicas
Los sistemas mecánicos de traslación pueden ser usados para modelar muchas situaciones
e involucran tres elementos básicos: masas (con masa M , medida en Kg), resortes (con
rigidez del resorte k, medida en Nm−1 ) y amortiguadores (con coeficiente de amortigua-
miento B, medido en Nsm−1 ).
Las variables asociadas son el desplazamiento x(t) (medido en m) y la fuerza F (t) (me-
dida en N). Convencionalmente, los elementos básicos son representados simbólicamente
de la siguiente manera:
Suponiendo que estamos tratando con resortes y amortiguadores ideales (esto es, supo-
niendo que se comportan linealmente), las relaciones entre las fuerzas y los desplazamien-
tos en el tiempo t son:
Dado que M = 1, F (t) = 4sen (ωt), F1 (t) = kx(t) = 25x(t) y F2 (t) = Bx′ (t) =
6x′ (t), se tiene que:
x′′ (t) + 6x′ (t) + 25x(t) = 4sen (ωt).
El flujo de corriente en el circuito está relacionado con la carga q(t) (medida en coulombs
C) mediante la relación:
i(t) = q ′ (t).
Las relaciones entre el flujo de corriente i(t) y la caı́da de voltaje v(t) a través de estos
elementos en el tiempo t son:
La interacción entre los elementos individuales que forman un circuito eléctrico está de-
terminada por las leyes de Kirchhoff. En particular, la ley de las tensiones de Kirchhoff
establece que:
12. Utilice la función escalón de Heaviside para expresar en forma compacta los pulsos si-
guientes y obtenga sus transformadas de Laplace.
6
(a) 100 −2 sen t si π ≤ t ≤ 2π
(b) f (t) =
0 si t < π ó t > 2π
- t
0 b
14. La entrada θi (t) y la salida θ0 (t) de un servomecanismo están relacionados por la ecuación
diferencial
θ0′′ + 8θ0′ + 16θ0 = θi , t ≥ 0
′ 1 − t si 0 < t < 1
e inicialmente θ0 (0) = θ0 (0) = 0. Para θi = f (t), donde f (t) =
0 si t > 1
s − 1 e−s
muestre que L {θi } = + 2 y luego obtenga una expresión para la respuesta del
s2 s
sistema en el tiempo t.
15. Durante el intervalo de tiempo [t1 , t2 ], una fuerza electromotriz constante e0 actúa en
un circuito serie RC. Supongamos que el circuito está inicialmente en estado de reposo,
verifique que la corriente en el circuito en el tiempo t es:
e0 −(t−t1 )/RC
i(t) = e H(t − t1 ) − e(t−t2 )/RC H(t − t2 ) .
R
16. Obtenga la trasnformada de Laplace F (s) de cada una de las siguientes funciones f (t)
cuyas gráficas son las siguientes:
(a) Onda rectangular periódica
A6
- t
0 T 2T 3T 4T
−A
6
A Z Z
Z Z
Z Z
Z Z
Z Z
ZZ
Z Z
Z
Z - t
0 T 2T 3T 4T
donde, para hacer que el sistema sea fı́sicamente realizable, los grados m y n de los
polinomios P (s) y Q(s) deben ser tales que n ≥ m.
Es importante destacar que, en general, una función de transferencia sólo se usa para
caracterizar un sistema lineal invariante en el tiempo. Es una propiedad del propio sistema
y es independiente tanto de la entrada como de la salida del sistema.
(a) La respuesta x(t) de un sistema invariante en el tiempo a una función de fuerza u(t)
está determinada por la ecuación diferencial
En el ejercicio anterior encontramos que para un sistema invariante en el tiempo que tiene
función de transferencia G(s), la respuesta x(t) del sistema inicialmente en reposo, a una
entrada u(t) está determinada por la relación en las transformadas
Si la entrada u(t) es la función impulso δ(t) entonces la respuesta del sistema estará
determinada por
X(s) = G(s)L {δ(t)} = G(s).
Es importante destacar el que teorema del valor inicial no dá el valor inicial de f (0), sino
más bien el valor de f (t) cuando t → 0+ .
Teorema del valor final : Si las funciones f (t) y f ′ (t) son transformables por Laplace y
existe lim f (t), entonces lim f (t) = lim sF (s).
t→+∞ t→+∞ s→0
(c) Use los teoremas del valor inicial y del valor final para encontrar el salto en t = 0 y
el valor lı́mite cuando t → +∞ para la solución del problema de valor inicial
7y ′ (t) + 5y(t) = 4 + e−3t + 2δ(t)
con y(0− ) = −1.
22. Obtenga la expresión del producto de convolución en cada uno de los siguientes casos:
23. Utilice el teorema de convolución para obtener las transformadas inversas de Laplace de
las funciones siguientes:
1 s
(a) F (s) = (b) F (s) = 2
s(s2 + 1) (s + 1)2
(a) La respuesta θ0 (t) de un sistema a una fuerza de transmisión θi (t) está dada por la
ecuación diferencial lineal:
d2 θ0 dθ0
2
+2 + 5θ0 = θi .
dt dt
i. Determine la respuesta al impulso del sistema.
ii. Use la integral de convolución para hallar la respuesta del sistema a una entrada
escalón unitario en el tiempo t = 0, suponiendo que está inicialmente en estado
de reposo.
iii. Confirme el resultado obtenido en (ii) por cálculo directo.
(b) Halle la respuesta del sistema a la entrada de pulso u(t) = A [H(t) − H(t − T )],
suponiendo que está inicialmente en reposo, caracterizado por la ecuación diferencial:
x′′ (t) + 7x′ (t) + 12x(t) = u(t).
c s
• L {c} = • L {cos (ωt)} =
s s2 + ω2
n! e−as
• L {tn } = • L {H(t − a)} = , a≥0
sn+1 s
1
• L ekt = • L {δ(t − a)} = e−as , a ≥ 0
s−k
Algunos resultados importantes
• Derivada de la transformada:
dn
L {tn f (t)} = (−1)n (L {f (t)}) = (−1)n F (n) (s)
dsn
• Primer teorema de traslación:
L eat f (t) = L {f (t)}|s−a = F (s − a)
• Transformada de derivadas:
L f (n) (t) = sn F (s) − sn−1 f (0) − sn−2 f ′ (0) − ... − f (n−1) (0)
• Transformada de integrales:
Z t
1 1
L f (x) dx = L {f (t)} = F (s)
0 s s
• Transformada de la convolución:
L {f (t) ∗ g(t)} = L {f (t)} L {g(t)} = F (s)G(s)
z+i
3. Dada la función definida por f (z) = .
z2 + 1
(a) Halle y clasifique las singularidades aisladas de f en el plano complejo extendido.
(b) Determine el cı́rculo de convergencia de la serie de Taylor de f centrada en 1 + i.
está dada por la ecuación diferencial de segundo orden θ0′′ (t) + 6θ0′ (t) + 10θ0 = θi (t), donde
t ≥ 0. Determine la respuesta del sistema en el tiempo t suponiendo que el sistema está
inicialmente en estado de reposo.
Z t
6. Halle la función y = y(t) sabiendo que y(t) = δ (t − 1) + 2 y (τ ) cos (t − τ ) dτ.
0
7. Obtenga:
n R √ o
t
(a) L e−t 0 τ sen 3τ dτ −1 e−3s
(b) L
s(s2 + 1)
Métodos numéricos
eso, muchas veces las respuestas a los modelos deben validarse con resultados experimentales y
observables en los mismos fenómenos.
Muchas veces sucede que no podemos hallar analı́ticamente las soluciones exactas a nuestros
problemas matemáticos, por ende tenemos que recurrir a una solución aproximada del problema.
Las aproximaciones numéricas introducidas al modelo matemático ahora reemplazarán el
problema matemático por un problema numérico.
Entonces, al tener que reemplazar nuestro problema matemático por un problema numérico,
una representación alternativa al esquema original de la Fig.6.1 es la siguiente:
accesibles para ser procesados por una computadora. Éste es uno de los objetivos principales
del area de estudio de los Métodos Numéricos.
Históricamente el descubrimiento y desarrollo de muchos métodos numéricos ha sido de
carácter empı́rico. Su estudio se ha consolidado como ciencia especı́fica a partir del desarrollo
de las computadoras y la tecnologı́a. El area de la Matemática que se encarga de estudiar los
fundamentos de los métodos numéricos y analizar sus errores se conoce como Análisis Numérico.
Esta area de la matemática moderna comprende a la vez de muchas sub-areas y se encuentra
en pleno desarrollo dado que existen problemas cientı́ficos que aún no se ha encontrado un
algoritmo satisfactorio.
Entre las distintas areas de las matemática, los métodos numéricos se aplican en:
Los métodos numéricos que desarrollaremos en este capı́tulo permitirán simplificar procedi-
mientos matemáticos de manera que podamos auxiliarnos con una computadora en problemas
como: resolver sistemas de ecuaciones lineales, hallar raı́ces de ecuaciones no lineales y de siste-
mas no lineales, calcular derivadas numéricamente, operar con matrices y aproximar soluciones
de ecuaciones diferenciales ordinarias. Ese será nuestro objetivo.
Las principales caracterı́sticas que tiene que tener un algoritmo son la eficiencia (rapidez
de cómputos) y la precisión (exactitud de los resultados). Sin embargo, hay más caracterı́sticas
deseables en un algoritmo numérico:
• Preciso y exacto: Las salidas del algoritmo deben aproximarse tanto a la solución del
problema como aproximarse entre sı́.
• Estable: Pequeños cambios en los datos de entrada (input) deben generar pequeños
cambios en los resultados de salida (output).
ENTRADA: x0 (aproximación
inicial), tol (tolerancia)
PASO 1: x1 ←− 21 x0 + x20
PASO 2: Si |x1 − x0 | ≤ tol, entonces TERMINAR.
En caso contrario, actualizar x0 ←− x1 y volver al PASO 1
PASO 3: Asignar la salida a x1
PASO 4: TERMINAR √
SALIDA: x1 aproximación de 2
x4 = 1.414213562374690
√
la aproximación de 2 = 1.414213562373095... con 12 dı́gitos exactos.
Es necesario aclarar que nuestro problema numérico en este caso consiste en calcular la raı́z
aproximada de dos con cierta tolerancia fijada de antemano mediante un esquema iterativo.
Este tolerancia se fija como un valor que acotación entre la distancia entre aproximaciones
sucesivas construidas por el esquema, esto es,
d(xk+1 , xk ) = |xk+1 − xk |.
Ahora bien, necesitamos poder implementar este como otros algoritmos en una computado-
ra. Para ello, introduciremos brevemente los programas (software) que más se usan para realizar
cálculos numéricos actualmente. En la Bibliografı́a o en los Links de interés puede encontrarse
más información sobre estos programas.
Congresos:
Asociaciones:
• Asociación Argentina de Matemática Aplicada, Computacional e Industrial (ASAMACI)
• Asociación Argentina de Mecánica Computacional (MECOM)
Revistas internacionales:
6.2.3 Expresiones
Para trabajar las expresiones tenemos varias formas:
FORMA 1: Con el operador asignación (=) asignamos valor(es) a variable(s)
variable = expresión
expresión
>> (1+sqrt(5))/2;
• Coma fija escalada con 4 cifras decimales (por defecto): format short
>> pi
ans = 3.1416 (3 cifras decimales exactas)
• Con el formato format short g, MATLAB elige la mejor notación entre coma fija y
punto flotante con 4 decimales, y con format long g realiza lo mismo con 15 decimales.
−1 0 1 1 + eps
A la vez, MATLAB tiene un formato especial para números muy grandes (que no puede
representar) considerados infinito:
>> 1.0/0.0
Warning: Divide by zero
ans = Inf
También tiene una representación especial para resultados que no están definidos como
números reales como N aN (Not a Number) o bien Inf (Infinity):
>> 0/0
Warning: Divide by zero
ans = NaN
>> Inf/Inf
ans = NaN
Con el comando complex creamos un número del campo complejo a partir de dos argu-
mentos, su parte real y su parte imaginaria:
>> z = complex(3,2)
z = 3.0000 + 2.0000i
La función disp imprime en pantalla un mensaje de texto o el valor de una variable pero
sin imprimir su nombre:
>> disp(’Fin del programa’)
Fin del programa
>> A=rand(2,5); (defino una matriz random 2×5)
>> disp(A)
0.9340 0.7577 0.3922 0.8235 0.3171
0.6787 0.7431 0.6555 0.6948 0.9502
La función fprintf imprime en pantalla mensajes permitiendo controlar los distintos valores
de formatos de salida. Su estructura es:
fprintf(’texto’)
La ’cadena de control’ de fprintf contiene los distintos formatos de escritura. Por ejemplo:
Para imprimir números con decimales como salida, se puede hacer de la siguiente manera:
>> raiz = 1.414213562374690;
>> fprintf(’Valor final aproximado: %12.5f \n’, raiz)
Valor final aproximado: 1.41421
Para imprimir un entero empleando el mismo formato se debe poner un cero después del
punto decimal:
>> iter = 9;
>> fprintf(’Cantidad de iteraciones: %9.0f \n’, iter)
Cantidad de iteraciones: 9
El comando fprintf se utiliza mucho también para escribir datos sobre un archivo de texto.
6.2.6 Operaciones
MATLAB puede operar distintas variables por medio de operadores. Los operadores aritméticos
son los siguientes:
+adición o suma
- sustracción o resta
* multiplicación
división
/
^ potenciación
\ división inversa
>> s= 2^.5
c = 1.4142 (raı́z cuadrada de 2)
>> c= 3 \ 1
c = 0.3333 (divide de derecha a izquierda)
>> 1==1
ans = 1 (verdadero)
>> 3.1415926==pi
ans = 0 (falso)
>> 3.1415926∼=pi
ans = 1
& y (AND)
| o no excluyente (OR)
&& y (AND con “cortocircuito”)
|| o no excluyente (OR con “cortocircuito”)
xor o excluyente
∼ negación lógica (not)
>> 1&1
ans = 1
>> 1&0
ans = 0
>> 1|0
ans = 1
>> xor(1,1)
ans = 0
Las matrices y vectores se definen como variables con nombres distinguiendo masyúsculas
de minúsculas. A diferencia de otros lenguajes de programación, en MATLAB no hace falta
preestablecer su dimensión dado que según el número de elementos definidos determina su
dimensión. Incluso su dimensión se puede modificar posteriormente.
Podemos definir un vector fila x separando sus componentes por espacios en blancos o
bien por comas (,):
>> x = [1 2 4] (o bien x = [1, 2, 4])
x = 1 2 4 (vector fila)
Si los números están separados por punto y coma (;) obtenemos el vector columna y:
>> y = [1; 3; 9]
y =
1
3
9 (vector columna)
Definimos las matrices por filas separando los elementos de una misma fila por espacios
en blancos o comas (,) y separando las filas por punto y coma (;):
>> A=[8 1 6; 3 5 7; 4 9 2]
A =
8 1 6
3 5 7
4 9 2
Internamente las matrices se almacenan por columnas y en una columna (aunque se in-
troduzcan por filas). Existen una gran variedad de funciones predefinidas en MATLAB que
generan de matrices predefinidas (zeros, ones, eye, rand, magic, ...). Usaremos luego.
>> Z=zeros(3,4)
Z =
0 0 0 0
0 0 0 0
0 0 0 0
>> U=ones(1,10)
Z =
1 1 1 1 1 1 1 1 1 1
También se puede crear una nueva matriz o vector a partir de otras matrices o vectores
definidos previamente. Sin embargo, debemos tener cuidado en que las dimensiones de conca-
tenación tiene que ser compatibles. En caso de no ser compatibles las dimensiones, MATLAB
dará un error advirtiendo esto. Veamos los siguientes ejemplos:
>> z = [x 8 16 32]
z = 1 2 4 8 16 32
>> B = [A; x]
B =
8 1 6
3 5 7
4 9 2
1 2 4
>> C = [B; z]
??? Error using ==> vertcat
CAT arguments dimensions are not consistent.
Accedemos a los elementos de un vector con el ı́ndice entre paréntesis, siempre y cuando el
mismo esté definido:
>> z(5)
ans = 16
>> z(7)
??? Attempted to access z(7); index out of bounds because numel(z)=6.
Accedemos a los elementos de una matriz con los ı́ndices entre paréntesis separados por una
coma:
>> B(2,3)
ans = 7
También podemos crear vectores o matrices muy rápidamente dando valores de algún/os
de sus elementos. Automáticamente se asignará cero a todos los demás elementos del vector o
matriz:
>> u(3) = .1, u(10) = .5
u = 0 0 .1 0 0 0 0 0 0 .5
>> C(3,5) = 1
C =
0 0 0 0 0
0 0 0 0 0
0 0 0 0 1
Es importante recordar que MATLAB diferencia entre vectores filas y vectores columnas.
Si x es el vector fila e y es el vector columna definidos anteriormente, entonces resultará un
error si queremos sumarlos (o restarlos):
>> x+y
??? Error using ==> +
Matrix dimensions must agree.
Lo correcto en dicho caso serı́a trasponer uno de los vectores dependiendo si queremos como
solución un vector fila o columna:
>> x+y’ (vector fila)
>> x’+y (vector columna)
Recordemos del Algebra Lineal que podemos definir el producto escalar interno entre
vectores del mismo espacio vectorial:
X
n
< x, y >= x.yT = xi yi
i=1
Recordemos también del Algebra Lineal que el producto matricial entre las matrices
M = (mij ) de dimensión r × s y la matriz N = (nij ) de dimensión s × t, está definido como:
X
s
pij = mik nkj i = 1, ..., r, j = 1, ..., t
k=1
En cambio, si hacemos:
>> M*N’
??? Error using ==> mtimes
Inner matrix dimensions must agree.
Utilizamos el operador dos puntos (:) para construir vectores como rangos de la forma:
Este operador también es muy potente en el uso de los ı́ndices de matrices. Por ejemplo:
>> A
A =
8 1 6
3 5 7
4 9 2
>> A(3, :)
ans = 4 9 2 (devuelve la 3er. fila de A)
>> A(:,2)
ans =
1
5
9 (devuelve la 2da. columna de A)
>> A(2:3, 2:3)
5 7
9 2 (devuelve una submatriz 2x2 de A desde (2,2) a (3,3))
Podemos usar este operador (:) para invertir el orden de elementos de un vector o permutar
filas/columnas de una matriz:
>> x = x(3:-1:1)
x = 4 2 1
>> D = A(:,3:-1:1)
D =
6 1 8
7 5 3
2 9 4 (permutamos las columnas de A)
Si en cambio usamos corchetes [ ] dentro del paréntesis, hacemos referencia a los ı́ndices
del vector o matriz:
>> A([1 3])
ans = 8 4
Las funciones intrı́nsecas anteriores de las Tabla 6.2 y Tabla 6.3 se aplican tanto a escalares,
vectores y matrices. Este acto de transformar un código desarrollado para una data de tipo
escalar y usarlo en una data de tipo no escalar, se llama vectorización.
Este proceso tiene la ventajas de impactar enormemente en la velocidad de cálculos (10 o
incluso 1000 veces más rápido), permite escribir en forma sintética y general códigos (funciones
y scripts .m). Además, MATLAB internamente tiene automatizado paralelizar operaciones
lógicas y aritméticas sobre datos, lo que se conoce como (paralelización automática implı́cita).
Las mayorı́a de las funciones intrı́nsecas están paralelizadas, por ende, son más rápidas y
eficientes. Si quisiéramos hacer esto en otros lenguajes habrı́a que programar estructuras de
bucles.
Ejemplos de vectorización: Operamos sobre los elementos de un vector o matriz usando
funciones intrı́nsecas de MATLAB:
>> sqrt([1,4; 9,16])
ans =
1 2
3 4 (raı́z cuadrada componente a componente)
>> abs([0,1,2,-5,-6,-7])
ans =
0 1 2 5 6 7 (valor absoluto componente a componente)
6.2.10 Programación
En MATLAB podemos programar y ejecutar nuestros programas muy fácilmente, sin embargo
como lenguaje de programación en sı́ mismo MATLAB no tiene tantas posibilidades como otros
lenguajes (Fortran, C++, etc).
Un programa se escribe en un archivo de texto *.m. Estos archivos se crean y modifican con
cualquier Editor de textos aunque es recomendable usar el editor propio Editor de MATLAB.
Existen dos grandes tipos de archivos *.m:
if [condición]
sentencias
end
Ejemplos de bifurcaciones tipo if:
if p=inf
disp(’norma infinito’)
end
La sentencia de bifurcación switch realiza una función análoga a un conjunto de if ... elseif
concatenados:
switch [expresión]
case caso 1
bloque 1
case {caso 2,..., caso n}
bloque 2
otherwise
bloque 3
end
Ejemplos de bifurcaciones tipo switch:
switch p
case inf
disp(’norma infinito’)
case 2
disp(’norma euclidea’)
otherwise
disp(’norma p-ésima’)
end
Primero se evalúa [expresión] que debe ser escalar o string. Se compara con caso 1 y si es
igual se ejecuta el bloque 1 de sentencias (es posible agrupar varios casos entre llaves). Ası́ se
sigue pero si ninguno es igual, se ejecutan las sentencias del caso otherwise.
Bucles: estos permiten repetir las mismas operaciones o comandos sobre datos distintos
(for, while). La sentencia de bucle for repite un conjunto de sentencias un número predeter-
minado de veces:
for variable=[rango]
sentencias
end
Ejemplos de bifurcaciones tipo for:
for i=1:10
a = a + 1
end
for k=20:-2:1
y(k) = k*sin(k*pi/4)
end
t = [-1:1e-2:1];
for p=[ 10 5.5 1 -3.5 -8 ]
y=@(t) p*cos(t)
plot(t,y(t)), hold on
end
while 1
iter = iter + 1;
if iter == 50, break
end
end
La sentencia break termina la ejecución del bucle interno que comprenden dicha sentencia.
6.2.11 Almacenamiento
En el trabajo con Métodos Numéricos aplicados a problemas reales o a problemas cientı́ficos,
muchas veces encontramos que nuestra data (input) viene en formato tipo: *.dat, *.txt, etc
dado que manejamos grandes cantidades de datos. Es importante poder leer dichos datos,
almacenarlos en alguna variable de nuestro entorno de trabajo en MATLAB y luego poder
trabajar con estos valores almacenados. También es importante que una vez que hayamos
hecho nuestros cálculos, la salida de nuestros resultados (output) a los cálculos numéricos rea-
lizados, se guarde en algún espacio de trabajo para que luego puede ser leı́do nuevamente y
analizado. Recordemos que cuando cerramos nuestra sesión, el contenido de las variables se
borra automáticamente. Para esto usaremos principalmente el comando load y el comando save.
El comando load: Para recuperar los datos almacenados en un archivo de texto (prefe-
rentemente plano) utilizamos el comando load que lee todo el archivo de datos y lo almacena
en una variable cuyo nombre coincide con el nombre del archivo.
De la siguiente manera guardamos todo el contenido del archivo entrada.dat en una va-
riable llamada datos sin importar cual sea la extensión del archivo:
>> load entrada.dat
Si en cambio deseamos que la variable tenga distinto nombre, por ejemplo queremos guardar
el contenido del archivo en una variable llamada A:
>> A = load(’matriz.dat’)
Si tipeamos:
>> load
recuperamos todas variables definidas en la sesión anterior de MATLAB guardadas previamen-
te con el comando save (ver abajo).
El comando save: Para grabar datos que en un archivo se utiliza el comando save.
De esta última manera, cada elemento de A se almacena con 8 cifras decimales. En cambio,
si deseamos almacenar 16 cifras decimales:
>> save salida.dat A -ascii -double
En formato ASCII sólo se guardan valores, no información como nombres de matrices, etc.
6.2.12 Scritps
Los scripts (archivos de comandos) son archivos *.m con una secuencia de comandos seme-
jante a la que se ingresan por tecleado. Los comandos se ejecutan secuencialmente cuando se
teclea el nombre del script (sin extensión) en la Ventana de Comandos o bien se puede ejecutar
con algún ı́cono de MATLAB.
El script debe almacenarse en el directorio actual donde estamos trabajando. Dentro del
mismo conviene poner punto y coma (;) al final de cada sentencia para evitar una salida
cuantiosa de las variables que usamos.
Un archivo (script o función) *.m puede llamar otros archivos *.m, sean funciones o scripts.
Esto hace que la programación sea de tipo modular.
Es importante remarcar que las variables definidas en scripts son variables internas, es decir
variables del espacio de trabajo desde el que se ejecuta el archivo. Por ejemplo, si el archivo se
ejecuta desde la Ventana de Comandos, se crean interactivamente en el Workspace de MATLAB
variables con el mismo nombre del script. Al terminar la ejecución del script, dichas variables
permanecen en memoria.
nombre función: Es el nombre de la función que tiene que coincidir con el nombre del archivo,
es decir, nombre.m. No debe incluir espacios ni palabras reservadas por MATLAB.
argumento(s) de entrada: Parámetro(s) que le paso a la función. Si hay más de uno, van
entre paréntesis separados por comas. Si no hay, no hace falta poner paréntesis.
argumento(s) de retorno: Resultado(s) devuelto por la función. Si hay más de uno van
entre corchetes separados por comas. Si sólo hay un valor no hace faltan corchetes. Si no hay
se omiten corchetes y el signo igual (=).
Los argumentos de entrada o valores de salida, dependerá de cómo esté armada nuestra
función. En general, usaremos:
function nombre funcion(variables entrada): cuando la función no devuelve ningún valor.
function variables salida = nombre funcion: cuando la función no precisa variables de entrada.
function [y] = nombre funcion(a,b,c): cuando la función devuelve un solo resultado.
function y = nombre funcion(a,b,c): se puede escribir de esta manera omitiendo corchetes.
function [y1,y2] = nombre funcion(a,b,c): Esta será la forma general que usaremos.
Una función no modifica los argumentos que recibe, estos pueden ser: strings, escalares,
vectores, matrices o incluso otras funciones. Los resultados de una función se obtienen a través
de los valores de retorno.
Los valores de salida deben calcularse dentro de la función. Las variables definidas dentro
de la función son variables locales. Dichas variables pertenecen a dicho espacio de trabajo y
no son vistas desde otros espacios de trabajo.
Opcionalmente (pero recomendable) desde la segunda lı́nea se ponen comentarios que deta-
llen la función (ver Figura 6.5). Esto genera información asociada a mi función que cree como
usuario. Para acceder a esta información podemos usar el comando Help. Cuando se teclea en
la Ventana de Comandos help nom func, se muestran los comentarios de las primeras lı́neas
del archivo nom func.m que comienzan con %. Esto genera un Help de mi función.
Opcionalmente (pero recomendable) se termina la función con end.
• Tipo inline: Son funciones de una lı́nea de comandos definidas por una cadena de
caracteres (string).
6.2.14 Ploteos
El comando plot: MATLAB dispone cuatro funciones básicas para crear gráficos 2-D:
• plot() crea un gráfico a partir de vectores y/o columnas de matrices con escalas lineales
sobre ambos ejes.
• semilogx() Idem con escala lineal en el eje de ordenadas y logarı́tmica en el de abscisas.
• semilogy() Idem con escala lineal en el eje de abscisas y logarı́tmica en el de ordenadas.
• loglog() Idem con escala logarı́tmica en ambos ejes.
En lo sucesivo usaremos la función plot, pero puede consultarse el Help las restantes en
modo similar. Veamos un ejemplo de un ploteo en 2D:
>> x = [1 3 2 4 5 3];
>> plot(x), grid on, (grid on activa una grilla en la gráfica)
>> x = [1 6 5 2 1]; y = [1 0 4 3 1];
>> plot(x,y), grid on
5 4
4.5 3.5
4 3
3.5
Figura 6.6: Ploteos sencillos
2.5
3 2
>>2.5
x = -1:.01:1; 1.5
>> 2
plot(x,exp(x),’g’), grid on 1
1.5 0.5
>> x = -2*pi:pi/100:2*pi;
1 0
>> 1
y = sin(x); z = cos(x);
1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6
8 1
0.8
7
0.6
6
0.4
5
Figura 6.7: Ploteos de funciones
0.2
4 0
−0.2
3
−0.4
2
Capı́tulo 6
1
2022
−0.6
172
−0.8
0 −1
−2 −1.5 −1 −0.5 0 0.5 1 1.5 2 −8 −6 −4 −2 0 2 4 6 8
CAPÍTULO 6. MÉTODOS NUMÉRICOS
Color Sı́mbolo
Azul b Tipo de marca Sı́mbolo
Tipo de linea Sı́mbolo
Verde g Puntos .
Continua -
Rojo r Cı́rculos o
A trazos –
Cyan c Cruces x
Punteada :
violet m Sumas +
A trazos y puntos -.
Amarillo y Asteriscos *
Negro k
>> x = 0:0.4:10;
>> y1 = sin(x).*exp(-0.4*x);
>> y2 = cos(x).*exp(-0.5*x);
>> plot(x,y1,’b+’,x,y2,’ro’), grid on
0.8
0.6
El comando subplot:
0.2
Una misma ventana de figura de MATLAB se puede dividir en m
particiones horizontales y n verticales. Para esto usamos la función subplot(m,n,i) donde:
0
• m: número de −0.2
subdivisiones de filas.
• n: número de subdivisiones
−0.4
0 1 2 3
de4 columnas.
5 6 7 8 9 10
1 1
0.5 0.5
0 0
−0.5 Figura
−0.5 6.9: Subplot 2×2
−1 −1
0 5 10 0 5 10
• X,Y: coordenadas
0
x e y en arreglos0
bidimensionales.
• Z: arreglo−0.5bidimensional con valores de la función.
• nivel: vector
−1
que contiene los niveles
−0.5
de contorno ó entero que representa el número de
0 5 10 0 5 10
niveles de contorno (se divide el valor mı́nimo y máximo de Z en m-1 niveles).
>> x = -2:0.2:2;
>> y = x;
>> [X,Y] = meshgrid(x,y); (uso meshgrid para hacer el "mallado")
>> Z = X.*exp(-X.^2-Y.^2);
>> mesh(X,Y,Z)
>> contour(X,Y,Z,10), grid on
1.5
1
Figura 6.10: Curvas de nivel de z = xe−(x
2 +y 2 )
0.5
0
Capı́tulo 6 2022 174
−0.5
CAPÍTULO 6. MÉTODOS NUMÉRICOS
También podemos utilizarlo para graficar funciones implı́citas f (x, y) = 0 en 2D. Por ejem-
plo, la intersección de dos cónicas, una parábola y una elipse:
>> x = -3:0.1:3;
>> y = x;
>> [X,Y] = meshgrid(x,y);
>> Z1 = X.^2-2*X-Y+0.5;
>> Z2 = X.^2+4*Y.^2-4;
>> contour(X,Y,Z1,[0,0],’b’), hold on,
>> contour(X,Y,Z2,[0,0],’r’), grid on,
El vector [0,0] en los argumentos de contour especifica el nivel de contorno. Solo nos interesa
el contorno de nivel z = 0. Los niveles de contorno se ingresan en forma de vector.
1
Figura 6.11: Intersección de la parábola y la elipse
El comando
0 mesh: Pasos para graficar el campo en dos variables
−1 z = f (x, y)
sobre un dominio
−2
rectangular [a, b] × [c, d].
2. Creamos dos matrices X (con filas copias de x) e Y (con columnas copias de y). Estas
matrices se crean con meshgrid y representan las coordenadas (x, y) de todos los puntos
de la retı́cula.
3. La matriz de valores Z del campo z se calcula con las matrices de coordenadas X e Y.
4. Ejecutamos la función mesh con las matrices X, Y y Z como argumentos. Se plotea una
superficie cuyos puntos tienen como cordenadas (X(i,j),Y(i,j),Z(i,j)).
>> x = -2:0.05:2;
>> y = -1:0.1:1;
>> [X,Y] = meshgrid(x,y);
>> Z = X.*exp(-X.^2-Y.^2);
>> mesh(X,Y,Z)
0.5
El comando surf: grafica en perspectiva la grafica de una función faceteada sobre una
retı́cula (el color depende del valor z). El color de las facetas depende del valor de la función y
−0.5
podemos agregar 1 una barra de color con escalas con colorbar.
0.5 2
0 1
>> surf(X,Y,Z), colorbar
−0.5
0
−1
−1 −2
0.4
0.3
0.5
0.2
−0.1
Capı́tulo 6 2022 176
−0.5 −0.2
1
0.5 2 −0.3
CAPÍTULO 6. MÉTODOS NUMÉRICOS
• Errores del modelo (model errors): son las simplificaciones que hacemos en la cons-
trucción del modelo que describa el fenómeno fı́sico para representarlo como un modelo
matemático y poder resolverlo.
• Errores numéricos (numerical errors): son los que introducimos al usar un algorit-
mo en una computadora. Están asociados al “algoritmo” (tenemos un tiempo finito de
ejecución) y al “hardware” (tenemos un espacio finito de almacenamiento).
Por otro lado, una vez planteado el método numérico para resolver dicho problema e im-
plementado el algoritmo en algún lenguaje, el proceso de cálculo introducirá errores dado que
la aritmética y la representación de números en la máquina no puede hacerse en forma exacta.
Esto introducirá otro tipo de errores que llamaremos errores de redondeo.
En la gráfica podemos suponer que los puntos son calculados por alguna técnica numérica,
mientras que el centro de las circunferencias concéntricas es el valor exacto. Si tomamos las
aproximaciones del item (a) y (c) vemos que ambas son inexactas, sin embargo, la última pre-
senta un esquema de cálculo más preciso ya que la distancia entre puntos es mucho menor. Es
decir, el item (a) es el resultado inexacto e impreciso, mientras que el item (c) representa un
esquema inexacto, aunque preciso. Por otro lado, si comparamos los items (a) y (b), vemos que
en este último el esquema se vuelve más exacto, aunque impreciso. Nuestra situación ideal es
el item (d) que proviene de un método exacto y preciso dado que los puntos están centradas
respecto al centro de las circunferencias y entre sı́. A continuación definiremos estos conceptos.
∆x = xe − xa .
∆x xe − xa
Ex = =
xe xe
siempre que xe 6= 0.
∆x xe − xa
Ex% = Ex ∗ 100 = ∗ 100 = ∗ 100
xe xe
siempre que xe 6= 0.
|xe − xa | 1
< × 10−d
|xe | 2
|∆y | = |ye − ya | = 3
|ye − ya |
|Ey | = = 0.000003 −
7 → Ey% = 0.0003%
|ye |
1
|Ey | < 0.000005 = × 10−5 − 7 → d = 5▲
2
Ejemplo 6.3.3. Sean los valores exactos y aproximados ze = 0.000012 y za = 0.000009,
donde ξ está entre 0 y x. El último término se conoce como fórmula del error en la forma
de Lagrange (consultar [Apos1984]). La diferencia entre la función exacta y el polinomio de
Taylor de orden n asociado en el punto x será el error de truncamiento de Taylor.
Xn
xk
En (x, a) = e −
x
.
k=0
k!
donde podemos obtener una estimación de la función en x con un número finito de términos.
Es claro que en la medida que n crece, el factorial en el denominador le gana a la potencia en
el numerador.
El script ExpTaylor.m explora esto ploteando el error relativo de la suma parcial como fun-
ción de la cantidad de sumandos n.
Código en MATLAB:
% Script: ExpTaylor
% Plotea, como funcion de n, el error relativo en la aproximacion de Taylor de
la funcion exponencial
% 1 + x + x ^2/2! +...+ x^n/n! a exp(x).
En la Figura 6.16 observamos que cuando aumentamos la cantidad de términos del polinomio
de Taylor n, el error relativo de la suma parcial de la serie disminuye. Esto coincide con
lo esperable matemáticamente puesto que si observamos la fórmula del error de Lagrange, el
factorial de n + 1 de denominador tiende a cero más rápidamente comparado con el numerador.
Sin embargo, cuando ploteamos las gráficas para x = 1, 5, 10, −1, −5, −10 como se muestra
en la Figura 6.17, observamos que el ploteo muestra que la “convergencia matemática” no
siempre es válida en máquina dado que los errores no convergen a cero a medida que la serie
aumenta la cantidad de términos. En algunos casos, el error se acuesta en un valor pequeño,
e incorporando más términos, no hacemos diferencia. Notamos también que el piso del error
relativo depende del valor de x: Errx=−10 > 1 × 10−10
Errx=10 < 1 × 10−15
Es decir, podemos notar que además del error analı́tico de truncamiento existe “otro error”
que aparece en este cálculo numérico de la exponencial.
x = 1.00
0
10
−2
10
−4
10
Error Relativo Suma Parcial.
−6
10
−8
10
−10
Figura 6.16: Error de Taylor de ex en función de n
10
−12
10
−14
10
−16
10
0 2 4 6 8 10 12 14 16
Orden Suma Parcial
x = 1.00 x = 5.00 x = 10.00
10 10 10
10 10 10
Error Relativo Suma Parcial
0 0 0
10 10 10
0 0 0
10 10 10
Figura 6.17: Error de Taylor de e en función de n para x = 1, 5, 10, −1, −5, −10.
x
donde
d1 = {1, 2, ..., 8, 9}
dj = {0, 1, ..., 8, 9} j > 1.
Entonces
f ltrunc (x) = ±0.d1 d2 ...dk × 10n
es su representación de coma flotante por truncamiento
d1 = {1, 2, ..., 8, 9}
dj = {0, 1, ..., 8, 9} 1 < j ≤ k,
y d∗k se calcula redondeando el número d1 d2 ...dk .dk+1 dk+2 ... al entero más próximo.
p(x) = (x − 1)6
Código en MATLAB:
% Script: ZoomPoli
% Plotea (x-1)^6 alrededor de x=1 con escala creciente
% pero evaluado via x^6 - 6x^5 + 15x^4 - 20x^3 + 15x^2 - 6x +1
clear, clc, close all;
k=0;
n=100;
for delta = [.1 .01 .008 .007 .005 .003 ]
x = linspace(1-delta,1+delta,n)’;
y = x.^6 - 6*x.^5 + 15*x.^4 - 20*x.^3 + 15*x.^2 - 6*x + 1;
k = k+1;
subplot(2,3,k), plot(x,y,x,zeros(1,n)), grid on
axis([1-delta 1+delta -max(abs(y)) max(abs(y))])
end
−7 −12 −13
x 10 x 10 x 10
1
2
5 0.5
1
0 0 0
−1
−5 −0.5
−2
−1
0.9 1 1.1 0.99 1 1.01 0.995 1 1.005
y
x
g(x) = √ √ , x≥0
x+1+ x
algebraicamente equivalentes. Vamos a comparar los resultados de calcular f (500) y g(500)
usando 4 cifras significativas con redondeo simétrico. ¿Deberı́amos obtener el mismo resultado
en nuestra máquina imaginaria? Normalizando antes de cada operación, obtenemos:
√ √
f (500) = 0.5000 × 10 3
0.5000 × 10 + 0.1000 × 10 − 0.5000 × 10
3 1 3
√ √
= 0.5000 × 103 0.5000 × 103 + 0.0010 × 103 − 0.5000 × 103
√ √
= 0.5000 × 103 0.5010 × 103 − 0.5000 × 103
= 0.5000 × 103 (22.3830292855... − 22.3606797749...)
= 0.5000 × 103 0.2238 × 102 − 0.2236 × 102
= 0.5000 × 103 0.0002 × 102
= 0.5000 × 103 0.2000 × 10−1
= 0.1000 × 102 = 0.1000 × 102 = 10.
0.5000 × 103
g(500) = √ √
0.5000 × 103 + 0.1000 × 101 + 0.5000 × 103
0.5000 × 103
= √ √
0.5000 × 103 + 0.0010 × 103 + 0.5000 × 103
0.5000 × 103
= √ √
0.5010 × 103 + 0.5000 × 103
0.5000 × 103
=
22.3830292855... + 22.3606797749...
0.5000 × 103 0.5000 × 103
= =
0.2238 × 102 + 0.2236 × 102 0.4474 × 102
= 1.11747575636... × 10 1
f (500) = 10,
g(500) = 11.17.
>> f(500)
ans = 11.174755300746853
>> g(500)
ans = 11.174755300747199
Notemos que incluso en MATLAB obtenemos una precisión de 16 cifras significativas y tam-
bién observamos la diferencia entre los cálculos con ambas fórmulas si calculamos los errores
absolutos y relativos anteriores.
>> exac = 11.174755300747198;
Figura 6.19: Dependencia del Error total del tamaño del paso h (step size)
Código en MATLAB:
function ssum = SerieSeno(x,tol,n)
% SerieSeno.m: Evalúa la representación en serie de la función SENO
% ssum = SerieSeno(x)
% ssum = SerieSeno(x,tol)
% ssum = SerieSeno(x,tol,n)
% ENTRADA:
% x: argumento de la función SENO
% tol: (optional) tolerancia para sumatoria (por defecto: tol = 5e-9)
% La serie termina cuando abs(T k/S k) < tol, donde:
% T k: k-ésimo término y S k suma incluyendo el k-ésimo término
% n: (optional) máximo número de terminos (por defecto: n = 15)
% RESULTADO:
% ssum: valor de la serie al cabo de n terminos o al satisfacer la tolerancia
if nargin < 2, tol = 5e-9;
end
if nargin < 3, n = 15;
end
term = x; % Inicializa la serie
ssum = term;
fprintf(’Valor aproximado de la función seno (%f)\n\n n término serie\ n’,x)
fprintf(’%3d %11.3e %12.8f \ n’,1,term,ssum)
for k=3:2:(2*n-1)
term = -term .* x.*x./(k.*(k-1)); % Próximo término en la serie
ssum = ssum + term;
fprintf(’%3d %11.3e %12.8f\n’,k,term,ssum)
if abs(term./ssum) < tol % Verifica convergencia
break
end
end
fprintf(’\n Error Total con %d términos es %g \ n\n’,(k+1)/2,abs(ssum-sin(x)))
end
Ventana de comandos en MATLAB:
>> sum = SerieSeno(pi/4,5e-9,100)
Valor aproximado de la función seno (0.785398)
n termino serie
1 7.854e-01 0.78539816
3 -8.075e-02 0.70465265
5 2.490e-03 0.70714305
7 -3.658e-05 0.70710647
9 3.134e-07 0.70710678
11 -1.757e-09 0.70710678
Error Total incluyendo 6 términos es 6.92801e-12
sum = 0.707106781179619
−6
x 10
4
3.5
En
2.5 caso de utilizar 15 términos no nulos en la sumatoria:
Error Total
|x|31 |π/2|31
2 |E31 (x)| ≤ ≤ ≈ 1.4611 × 10−28 .
31! 31!
1.5
0.5
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
x
−16
x 10
3.5
3
Figura 6.21: Error total con 15 términos.
2.5
1.5
CAPÍTULO 6. MÉTODOS NUMÉRICOS
• Propagación del error en la suma y resta. Sean los números exactos x e y con
valores aproximados x̃, ỹ, y sean los errores absolutos δx , δy respectivamente. Esto es:
x = x̃ + δx , y = ỹ + δy .
Por lo tanto, el error relativo del producto de dos (o más magnitudes) es igual a la suma
de los errores relativos de dichas magnitudes.
Análogamente, esto también vale para el cociente. Consultar [MaFi2000].
x ' ±q × β n
x ' ±q × 2n
donde 1
2
≤ |q| < 1.
Sin embargo, en un computadora sólo tenemos una cantidad finita de espacio de números
binarios para el almacenamiento de la mantisa q y el exponente n. Por ende, no podemos
representar todos los números reales, sino una cantidad finita de los mismos.
Para comprender esto, veamos a continuación una máquina virtual “más chica” que las
computadoras actuales.
(0.d1 d2 d3 d4 )2 × 2n
s →
7 d1 d1 ∈ {0, 1},
e → 7 d1 d2 . . . d11 di ∈ {0, 1} i = 1, . . . , 11,
q → 7 d1 d2 . . . d53 d1 = 1, di ∈ {0, 1} i = 2, . . . , 53.
0 ≤ (q)2 < 1
0 ≤ (q)10 ≤ 252 − 1.
En realidad, como siempre vale que d1 = 1, el formato IEEE permite almacenar información
de una palabra de 65 de bits en una de 64 bits.
Teniendo 11 bits para el exponente, estos representarı́an un rango decimal de:
Pero como (e)10 = 0 y (e)10 = 2047 están reservados para números especiales, el rango exacto
del exponente en base 10 es:
1 ≤ (e)10 ≤ 211 − 2 = 2046.
Para que (e)10 tenga signo negativo, corremos el rango restando 1023 a la representación, y se
reacomoda para que el exponente sea un entero en el intervalo:
y obtenemos
−1022 ≤ (e′ )10 ≤ 1023.
Estos serán los posibles valores del exponente en base 10 en la representación de punto flotante.
Existen exponentes y números especiales reservados en MATLAB, por ejemplo los valores
reservados del exponente:
(e)10 = 0 → todos ceros
(e)10 = 0 → todos unos
En la siguiente Tabla 6.4, vemos las demás combinaciones:
Binario Decimal
eps 2 ˆ(-52) 2.220446049250313e-16
realmin 2 ˆ(-1022) 2.225073858507201e-308
realmax (2-eps)*2ˆ(1023) 1.797693134862316e+308
Cualquier número por encima de realmax generará un overflow , y cualquier número debajo
de realmin generará un underflow . Sin embargo, algunas máquinas permiten números de
punto flotante subnormales en el intervalo [eps*realmin,realmin]. Esto permite llenar (en
parte) el gap (salto) entre 0 y realmin. El menor número subnormal positivo de punto flotante
es: 4.940656458412465e-324. Cualquier valor menor es cero.
A modo de interés general y para resaltar la importancia del estudio de los métodos
numéricos en problemas reales, se linkean algunos desastres históricos causados por errores
en los cálculos de los métodos:
• 25 de febrero de 1991 (Arabia Saudita): Falla del misil Patriot en sistema de defensa.
Causa: acumulación de errores de redondeo.
• 4 de junio de 1996 (Holanda): Explosión del cohete Ariane 5.
Causa del desastre: error de “overflow”.
• Conversión aritmética del EURO.
• Enero de 1982 - Noviembre de 1983 (Canadá): Bolsa de valores de Vancouver.
Causa: Rendodeo truncado en vez de simétrico.
• 5 de abril de 1992 (Alemania): Error de redondeo cambia al Parlamento durante la
elección de Schleswig-Holstein.
• 23 de agosto de 1991 (Holanda): El hundimiento la plataforma offshore Sleipner A.
Causa del desastre: Aproximación inexacta del método de elementos finitos.
Matricialmente
a11 a12 . . . a1n x1 b1
a21 a22 . . . a2n x2 b2
.. .. .. . . .. = .. .
. . . .. . .
an1 an2 . . . ann xn bn
Nuestro problema numérico reside en que, dada una matriz A ∈ Rn×n y un vector (columna)
lado derecho b ∈ Rn , queremos encontrar el vector (columna)x ∈ Rn tal que:
Ax = b.
Una condición necesaria y suficiente para que el SEL tenga solución es:
det(A) 6= 0
pues existe A−1 y la solución de dicho sistema viene dada por:
x = A−1 b.
Calcular esta solución puede ser complicado desde el punto de vista numérico aún en casos
“simples”. Ahora bien, ¿cómo resolvemos Ax = b numéricamente? Existen dos grandes tipos
de métodos para resolver este problema:
• Métodos Directos: Son métodos numéricos que se ejecutan en un número finito de
pasos y generarı́an la solución x̃ exacta si no fuera por los errores de redondeo (no hay
truncamiento). Por ejemplo: Eliminación Gaussiana, Factorización LU, Choslesky y
otros.
• Métodos Iterativos: Son métodos numéricos que generan iterativamente una sucesión
de vectores {xk } que converge a una solución aproximada x∗ con una tolerancia deseada
y donde es necesario considerar el error de truncamiento. Por ejemplo: Jacobi, Gauss-
Seidel, Sobre-relajación Sucesiva (SOR), Métodos de tipo Krylov (GMRES, etc) y otros.
Eliminación Gaussiana: La idea consiste en realizar una cantidad finita de pasos donde
en cada paso realizamos operaciones elementales entre matrices para obtener sistema de ecua-
ciones lineales que tengan el mismo conjunto solución e ir triangularizando superiormente la
matriz de coeficientes original. Si suponemos que la matriz de coeficientes es cuadrada con n
ecuaciones y n incógnitas, luego de (n − 1) pasos obtenemos una matriz triangular superior
equivalente a la original y con el método de la sustitucion regresiva obtenemos la solución del
SEL. Todo este proceso se conoce como eliminación Gaussiana y en general es muy eficiente
para resolver sistemas de ecuaciones lineales sobretodo cuando las matrices son “llenas”.
• Sustitución: Reemplazar una fila por la suma de esa fila más un múltiplo escalar de
cualquier otra fila
rk + mrj 7−→ rk .
El siguiente Teorema afirma que realizar operaciones elementales entre matrices, se generan
sistemas de ecuaciones que son equivalentes.
PASO 1
• La 1era. fila se usa para eliminar los elementos de la 1era. columna debajo de la diagonal
principal.
• El elemento a11 = 10 se usa como pivote.
• Los valores mk1 con k = 2, 3 son los multiplicadores:
a11 = 10 (pivote) 10 −7 0 7
m21 = a21 /a11 = −0.3 −3 2 6 4
m31 = a31 /a11 = 0.5 5 −1 5 6
10 −7 0 7
r2 − m21 r1 7−→ 0 −0.1 6 6.1
r3 − m31 r1 7−→ 0 2.5 5 2.5
En rojo son los coeficientes de la matriz que se modifican en cada paso del algoritmo.
PASO 2
• Usamos la 2da. fila para eliminar los elementos de la 2da. columna debajo de la diagonal
principal.
• El elemento a22 = −0.1 se usa como pivote.
• Los valores mk2 con k = 3 es el multiplicador:
10 −7 0 7
a22 = −0.1 (pivote) 0 −0.1 6 6.1
m32 = a32 /a22 = −25 0 2.5 5 2.5
10 −7 0 7
0 −0.1 6 6.1
r3 − m32 r2 7−→ 0 0 155 155
Usando sustitución regresiva:
155x3 = 155 =⇒ x3 = 1.
Sustituyendo en la 2da. ecuación:
−0.1x2 + 6(1) = 6.1 =⇒ x2 = −1.
Sustituyendo en la 1er. ecuación:
10x1 + (−7)(−1) = 7 =⇒ x1 = 0.
Finalmente, la solución al sistema es:
0
x = −1 .
1
arq r = q + 1, q + 2, . . . , n
• Los números
arq
mrq = r = q + 1, q + 2, . . . , n
aqq
con los que se multiplica la fila pivote para restarle las correspondientes filas
posteriores se llaman multiplicadores.
El siguiente Teorema asegura que si la matriz es inversible, entonces puedo aplicar el método
de eliminación Gaussiana y obtener una matriz triangular superior inversible.
Ax = b
aqq = 0
éste aqq no puede usarse como pivote. Entonces, buscamos en la columna k-ésima el coeficiente
akq 6= 0 para algún k entre q + 1 ≤ k ≤ n. Una vez hayado dicho valor, intercambiamos la
fila q-ésima con la k-ésima. Este proceso se conoce como pivoteo trivial y aparece en nuestra
implementación en las siguientes sentencias.
% Verificación
if (Aum(q,q) == 0)
disp(’A no es inversible. No hay solución o no es única’);
break
end
|aiq |
(multiplicadores) 0 ≤ |miq | = ,≤ 1 i = q, q + 1, ..., n.
aq̃q
Notemos que los multiplicadores caen en el intervalo [0, 1] donde hay mayor densidad de repre-
sentación en punto flotante. Esto permite reducir el error de redondeo en la representación de
los números.
Para evitar más aún la propagación de errores de redondeo usamos otra estrategia conocida
como pivoteo parcial escalado.
• Error de redondeo: Puede ser serio. Se recomienda el uso ”Pivoteo parcial” o ”Pivoteo
parcial escalado” para minimizarlo.
2 3 2 3
F lops 3
n + o(n2 ) n2 + o(n) 3
n + o(n2 )
Definición (Factorización LU). Diremos que una matriz A invertible admite una
factorización LU si podemos expresarla como el producto
A = LU
siendo:
Suponiendo que A admite una factorización LU, ¿cómo resuelvo Ax = b? Notemos que
Ax = (LU)x = L(Ux) = Ly = b
o
donde y = Ux. Por lo tanto 1 ) Ly = b (hallo y)
2o ) Ux = y (hallo x)
Si la eliminación Gaussiana para resolver el SEL Ax = b pueda llevarse a cabo hasta el
final sin intercambio de filas entonces es posible realizar la factorización LU. Es importante
porque una vez hallada la factorización LU podemos resolver el SEL para distintos vectores b.
Veamos cómo determinamos las matrices triangulares L y U.
Ejemplo 6.4.5. Volvamos nuevamente al ejemplo 3 × 3 original:
10 −7 0 x1 7
−3 2 6 . x2 = 4 .
5 −1 5 x3 6
Matricialmente hicimos:
10 −7 0 1 0 0 10 −7 0
0 −0.1 6 = −0.3 1 0 −3 2 6 = M1 A
0 2.5 5 0.5 0 1 5 −1 5
10 −7 0 1 0 0 10 −7 0
0 −0.1 6 = 0 1 0 0 −0.1 6 = M2 (M1 A)
0 0 155 0 −25 1 0 2.5 5
Por ende
U = M2 M1 A.
Suponiendo que las matrices triangulares inferiores M1 y M2 son inversibles (lo son), entonces
M−1 −1
1 M2 U = A.
Llamando L1 = M−1 1 y L2 = M−1 2 a las inversas de las matrices triangulares inferiores (que
también serán triangulares inferiores)
L1 L2 U = A =⇒ LU = A
siendo L = L1 L2 matriz triangular inferior por ser producto de matrices triangulares inferiores.
Por ende, los elementos de L y U son calculados en el proceso de eliminación Gaussiana.
El siguiente Teorema nos garantiza este hecho:
Teorema 6.4.3. Consideremos el SEL Ax = b donde A de orden n × n. Supongamos
que podemos realizar el proceso de eliminación Gaussiana sin intercambio de filas.
Entonces A puede factorizarse como el producto de una matriz triangular inferior L y
una matriz superior U. Esto es, A = LU.
Una vez construidas L y U, usamos la sustitución regresiva y progresiva para resolver este
último SEL equivalente al original.
Resolviendo por eliminación Gaussiana con pivoteo parcial, tuvimos que permutar en el PASO
2 la 2da. fila con la 3era. dado que el pivote lo encontramos de calcular:
max{|−0.001| , |2.5|}.
Entonces
10 −7 0 7
a22 = 2.5 (pivote) 0 2.5 5 2.5 .
m32 = a32 /a22 = −0.0004 0 −0.001 6 6.001
Matricialmente
10 −7 0 1 0 0 10 −7 0
0 −0.001 6 = −0.3 1 0 −3 2.099 6 = M1 A,
0 2.5 5 0.5 0 1 5 −1 5
10 −7 0 1 0 0 10 −7 0
0 2.5 5 = 0 0 1 0 −0.001 6 = PM1 A.
0 −0.001 6 0 1 0 0 2.5 5
Notemos que es este caso estamos introduciendo una nueva matriz en el proceso matricial.
Esta matriz P se conoce como matriz de permutación
1 0 0
P = 0 0 1 .
0 1 0
Cada vez que intercambio filas en la eliminación Gaussiana matricialmente estoy multiplicando
por una matriz de permutación.
Definición (Matriz de permutación). Una matriz de permutación P de orden
n × n es tal que cada fila y cada columna tiene solo un elemento igual a 1. Es decir,
es la matriz identidad con filas y columnas intercambiadas.
P−1 = P .
T
0 0 1 0
Si hacemos PA permutamos la 1era. con 2da. fila de A y también la 3era. con la 4ta.
• Sean Mk con k = 1, . . . , n−1 las matrices unitarias triangulares inferiores que se obtienen
con los multiplicadores debajo de la k-ésima diagonal en el paso k-ésimo.
U = Mn−1 Pn−1 . . . M2 P2 M1 P1 A
• Llamando: −1
L = M′n−1 . . . M′2 M′1 = L1 L2 . . . Ln−1
a la matriz triangular inferior que contiene los multiplicadores donde Lk = M′−1
k
• Y llamando
P = Pn−1 . . . P2 P1
la matriz de permutación que contiene todos los intercambios.
• Obtenemos que
U = L−1 PA
o equivalentemente
PA = LU
que se conoce como la factorización indirecta LU de A.
PA = LU.
1o ) Ly = Pb
2o ) Ux = y
MATLAB tiene incorporada una función intrı́nseca que usando el método de eliminación
Gaussiana con pivoteo parcial devuelve la matriz triangular inferior L, la matriz triangular
superior U y la matriz de permutación P.
El comando lu de MATLAB:
>> [L,U] = lu(A)
Devuelve una triangular superior U y una matriz L producto de la matriz triangular inferior y
la matriz de permutación.
>> [L,U,P] = lu(A)
Devuelve la triangular inferior L con diagonal de unos, la triangular superior U y la matriz de
permutación P de modo que puedo reconstruir P*A = L*U.
kxk > 0 ∀x 6= 0.
kxk = 0 sı́ y solo sı́ x = 0.
kλxk = |λ|kxk λ ∈ R.
kx + yk ≤ kxk + kyk. (desigualdad triangular)
Observaciones:
• La noción de norma generaliza el concepto de longitud.
• Trabajaremos con los espacios vectoriales V = Rn (vectores) o V = Rm×n (matrices).
Sobre estos espacios definiremos algunas normas muy usadas en el cálculo numérico.
B = {x ∈ V : kxkp ≤ 1},
También pueden definirse normas matriciales, o sea, normas en el espacio de matrices Rm×n .
X
n
p = 1 7−→ kAk1 = max |aij | (Norma 1)
j=1,...,n
i=1
!1/2
X
n
p = 2 7−→ kAk2 = |aij |2 (Norma Euclı́dea)
i,j=1
X
n
p = ∞ 7−→ kAk∞ = max |aij | (Norma infinito)
i=1,...,n
j=1
La noción de iteración:
Una técnica fundamental de resolución en los Métodos Numéricos es la idea de iteración.
Esta consiste en repetir un proceso determinado para obtener una sucesión de valores hasta
un resultado suficientemente aproximado (que idealmente converge a la solución). Se necesita
un esquema iterativo (regla, fórmula, función), una condición inicial y un tolerancia de
corte. Estos esquemas se utilizan en hallar raı́ces de ecuaciones, soluciones de sistemas de
ecuaciones lineales y no lineales, entre otros. Nosotros estudiaremos el proceso de iteración que
consiste en sustituir repetidamente un valor obtenido en una misma fórmula.
x = G (x) .
p(k) −→ r
que vemos converge al vector x = [2, 4, 3]T que es la solución del SEL inicial.
Ejemplo 6.4.11. Usando la misma aproximación inicial x0 = [1, 2, 2]T obtenemos la sucesión:
h iT
(0) (0) (0)
x0 = x1 , x 2 , x 3 = [1.0000, 2.0000, 2.0000]T
h iT
(1) (1) (1)
x1 = x1 , x 2 , x 3 = [1.7500, 3.3750, 2.9200]T
h iT
(2) (2) (2)
x2 = x1 , x 2 , x 3 = [1.9500, 3.9688, 2.9862]T
h iT
(3) (3) (3)
x3 = x1 , x 2 , x 3 = [1.9956, 3.9961, 2.9990]T
h iT
(4) (4) (4)
x4 = x1 , x 2 , x 3 = [1.9993, 3.9995, 2.9998]T
h iT
(5) (5) (5)
x5 = x1 , x 2 , x 3 = [1.9999, 3.9999, 3.0000]T
h iT
(6) (6) (6)
x6 = x1 , x 2 , x 3 = [2.0000, 4.0000, 3.0000]T
que converge al vector solución x = [2, 4, 3]T . Notemos que lo hace en menos iteraciones.
Nos preguntamos ahora, ¿cuando convergen estos esquemas iterativos? Es decir, bajo qué
condiciones podemos tener asegurarada la convergencia de estos esquemas para cierta tolerancia
y cierto punto de arranque. Para responder esto, introduciremos un concepto muy importante
del Análisis Numérico.
• Las iteraciones de Jacobi y Gauss-Seidel aplicadas al mismo SEL pueden ser ambas con-
vergentes, ambas divergentes, o bien una cualquiera de ellas convergente y la otra no.
• Si bien estos algoritmos son fáciles de entender e implementar, son poco usados en apli-
caciones. En general, los métodos iterativos más usados para resolver grandes SEL son
los de tipo Krylov. Consultar [Saad2003].
El vector residual:
Al resolver en una máquina un sistema de ecuaciones lineales Ax = b con precisión finita,
obtenemos una aproximación x̃ del vector solución x:
x̃ = x + ∆x
que se verá afectada por errores de redondeo (en los métodos directos), y por errores de redon-
deo y truncamiento (en los métodos iterativos). Ahora queremos conocer cuánto difiere nuestra
solución x̃ de la solución exacta x. Para esto, introduciremos los siguientes conceptos.
r = b − Ax̃
κ(A) = kAk.kA−1 k
El siguiente Teorema establece una cota para el error relativo de la aproximación de un SEL
que depende del número de condición.
Teorema 6.4.6. El error relativo al resolver Ax = b está acotado por:
f (r) = 0
para cierta función f definida en (a, b) tal que f, f ′ , f ′′ continuas en un intervalo que contenga
un entorno (r − δ, r + δ). En la Figura 6.25 vemos el ejemplo de f (x) = ex − x − 2 donde
x ∈ [−3, 3].
16
exp(x) − x − 2
14
2
CAPÍTULO 6. MÉTODOS NUMÉRICOS
f (x)
r =x+h≈x− .
f ′ (x)
Esta expresión nos sugiere un esquema de iteración. La iteración de N-R será de la forma:
f (xn )
xn+1 = xn − , n ≥ 0,
f ′ (xn )
siendo x0 una aproximación inicial.
16
f(x) = exp(x) − x − 2
14 tangente en (x0,f(x0))
y=f(x)
12 Figura 6.26: Interpretación geométrica de N-R
10
Geométricamente, si considaremos la primera iteración:
8
f (x0 ) 0 − f (x0 )
x1 = x0 − =⇒ f ′ (x0 ) =
6 ′
f (x0 ) x 1 − x0
f(x )
es la4 pendiente de la recta tangente a la gráfica de la 0función que pasa por los puntos: (x0 , f (x0 ))
y (x21 , 0). Es decir, calculamos la pendiente de la recta tangente a la gráfica en el punto x0 , donde
esta recta corta al eje de las ordenadas tenemos la primer aproximación x1 . Ası́ sucesivamente.
0
x1 x0
Capı́tulo
−2 6 2022 226
−4
−4 −3 −2 −1 0 1 2 3 4
CAPÍTULO 6. MÉTODOS NUMÉRICOS
√
Ejemplo 6.5.1. Calcular la raı́z cuadrada M siendo M > 0.
Definamos la función
f (x) = x2 − M.
Como f infinitamente derivable y f ′ (x) = 2x,
f (xn ) x2 − M 1 M
xn+1 = xn − ′ = xn − n = xn + .
f (xn ) 2xn 2 xn
Para el caso M = 2,
1 2
xn+1 = xn + , n ≥ 0.
2 xn
Notemos que este es nuestro ejemplo rcuad2.m visto en la sección introductoria.
Notemos que pasamos como argumento del N-R las funciones f y su derivada f ’ definidas
en forma de funciones anónimas.
Ejemplo 6.5.3. Aproximar la raı́z de g(x) = ex − x − 2.
>> g=inline(’exp(x)-x-2’);
>> dg=inline(’exp(x)-1’);
>> x=NewtonRaphson(g,dg,1,1e-5,30)
Notemos que ahora pasamos como argumento de N-R las funciones g y su derivada g’ defi-
nidas ”en lı́nea“.
Ahora bien, un problema que aparece en el método de N-R es que necesitamos conocer
explı́citamente una expresión para la derivada de la función f que me define la ecuación no
lineal. Pero, ¿cómo aproximamos numéricamente f ′ (x)? Una forma natural (aunque mala) de
hacerlo es mediante el cociente incremental
f (x + h) − f (x)
f ′ (x) ≈ .
h
Usando Taylor, el error que cometemos con esta fórmula es
1
f (x + h) = f (x) + hf ′ (x) + h2 f ′′ (ξ).
2
Reacomodando
f (x + h) − f (x) 1 ′′
f ′ (x) = − hf (ξ).
h 2
Por ende, el error de truncamiento es
1
− hf ′′ (ξ) = o(h).
2
Es importante obsvervar que, si bien bajo ciertas hipótesis podemos construirnos un esquema
numérico, no siempre hay garantı́as de que el método numérico converja a la raı́z de nuestra
ecuación. Por ejemplo, puede suceder que haya una oscilación divergente o un ciclo como
muestran las Figuras 6.27 y 6.28 siguientes.
1.5
0.5
0
Figura 6.27: Ejemplo de divergencia de N-R.
−0.5
−1
−1.5
−2
−20 −15 −10 −5 0 5 10 15 20
Los siguientes Teoremas nos brindan condiciones suficientes para asegurar la convergencia
del método de N-R tanto local como globalmente.
Teorema 6.5.1. Carácter local. Supongamos que f ∈ C 2 (a, b) y que existe una raı́z
r de f en (a, b). Si f ′ (r) 6= 0 entonces existe un δ > 0 tal que la sucesión {xk }k≥0
definida por el proceso iterativo de Newton-Raphson:
f (xn )
xn+1 = xn − n≥0
f ′ (xn )
|xn+1 − r| |en+1 |
α = 'C cuando n→∞
|xn − r| |en |α
Ejemplo 6.5.5. Aproximar las raices de h(x) = x3 − 3x + 2. Se sabe que r1 = −2 raı́z simple
y r2 = 1 raı́z doble.
Es importante que grafiquemos la función para poder observar en qué lugar aproximadamente
estarán la o las raı́ces.
>> x = -3:1e-2:3;
>> h = inline(’x.^3-3*x+2’);
>> plot(x,h(x)), grid on
20
15
10
0
Figura 6.29: Raı́z simple y doble de h(x) = x3 − 3x + 2
>> x = NewtonRaphsonDeriv(h,-2.5,1e-5,30);
−10
Metodo iterativo de Newton-Raphson con derivada numerica
−15
Iter x f(x) Error
1 -2.111111 -1.075446 0.388889
−20
2 -2.007407
−3 −2 -0.066996
−1 0.103704
0 1 2 3
Resultado final:
Raiz aproximada -2.000000
Resultado final:
Raiz aproximada 1.000000
Notemos que para alcanzar la misma tolerancia e inicializando con una aproximación inicial
tan cercana como en el caso anterior, tuvimos que iterar un mayor número de veces. Esto se
debe a que la raı́z r2 = 1 es doble, es decir, h(r1 ) = h′ (r1 ) = 0 y h′′ (r1 ) 6= 0.
El método de N-R presenta una velocidad de convergencia lineal para raı́ces múltiples y
convergencia cuadrática para raı́ces simples, como es el primer caso.
Veamos algunas consideraciones sobre el método de N-R:
VENTAJAS:
• La convergencia de N-R cerca de una raı́z simple es cuadrática.
• Modificando el esquema iterativo de N-R para una raı́z múltiple puede recuperarse la
convergencia cuadrática.
DESVENTAJAS:
• Puede darse la división por cero (o casi cero) si f ′ (xn ) ' 0.
Sabemos de Algebra y Geometrı́a Analı́tica que estas ecuaciones son ecuaciones generales de
segundo grado que definen implı́citamente curvas en el plano z = 0 como se muestra en la Figura
6.30. Una raı́z del SENL es un punto del plano (p∗ , q ∗ ) ∈ R2 que satisface simultáneamente
ambas ecuaciones. Geométricamente son los puntos donde se intersecan ambas curvas
f1 (p∗ , q ∗ ) = 0,
f2 (p∗ , q ∗ ) = 0.
Nuestro objetivo en esta Sección es desarrollar un método numérico para aproximar estas raı́ces
iterativamente. Es decir, contruiremos una sucesión de puntos del plano que converja a la o las
raices del SENL.
parabola
2
elipse
1.5
0.5
−0.5
x2 − y + 0.5
f1 (x, y) = 0 ⇔ x = ,
2
−x2 − 4y 2 + 8y + 4
f2 (x, y) = 0 ⇔ y = .
8
Definimos a continuación las funciones llamadas generatrices:
x2 − y + 0.5
g1 (x, y) := ,
2
−x2 − 4y 2 + 8y + 4
g2 (x, y) := .
8
Tomando una aproximación inicial (p0 , q0 ), generamos la sucesión {(pk , qk )}k≥0 tal que
p2k − qk + 0.5
pk+1 = g1 (pk , qk ) = k ≥ 0,
2
−p2k − 4qk2 + 8qk + 4
qk+1 = g2 (pk , qk ) = k≥0
8
Deseamos que esta sucesión de pares ordenados converjan a alguna de las raı́ces del SENL dado
inicialmente.
Observación: Notemos además que hay otras posibles generatrices que podemos definir a
partir de los campos escalares inciales, como por ejemplo
−x2 + 4x + y − 0.5
ge1 (x, y) := ,
2
−x2 − 4y 2 + 11y + 4
ge2 (x, y) := .
11
Debido a que nuestro esquema iterativo de Punto Fijo necesita de la definición de las gene-
ratrices obtenidas, veamos primero cómo definir nuestras generatrices en un m-file:
Este m-file define las generatrices en un archivo separado que luego pasaremos como argu-
mento a la función PuntoFijo.m como cadena de caracteres (string).
Iter Pk Qk
0 2.0000000 0.0000000
1 2.2500000 0.0000000
2 2.7812500 -0.1328125
3 4.1840820 -0.6085510
4 9.3075467 -2.4820360
5 44.8062311 -15.8910907
6 1011.9947186 -392.6042650
7 512263.2073904 -205477.8225378
.
. .
.
. .
13 Inf -Inf
14 Inf -Inf ¡DIVERGE!
15 Inf -Inf
En este caso observamos que las iteraciones sucesivas no converjen a la raı́z que se observa
en la gráfica, aún cuando el punto de arranque lo tomemos cerca de dicha raı́z. Incluso, si
probamos con otros puntos de arranque cerca de dicha raı́z, veremos que el esquema diverge
también. Será necesario entonces tomar otras generatrices.
Consideremos las generatrices ge1 y ge2 dadas anteriormente:
−x2 + 4x + y − 0.5
f1 (x, y) = 0 ⇔ x = ge1 (x, y) = ,
2
−x2 − 4y 2 + 11y + 4
f2 (x, y) = 0 ⇔ y = ge2 (x, y) = .
11
Consideremos la sucesión obtenida a partir de estas:
−p2k + 4pk + qk − 0.5
pk+1 = ge1 (pk , qk ) = , k ≥ 0,
2
−p2k − 4qk2 + 11qk + 4
qk+1 = ge2 (pk , qk ) = , k ≥ 0,
11
tomando alguna aproximación inicial (p0 , q0 ).
En este caso, observamos que si bien el esquema corta por número de iteraciones, las apro-
ximaciones sucesivas generadas por G2.m converjen a la otra raı́z partiendo del mismo punto
de arranque que en el caso con G1.m. Podemos verificarlo evaluando nuestros campos (f1 , f2 )
en el valor aproximado (p∗ , q ∗ ).
p∗ = g1 (p∗ , q ∗ )
q ∗ = g2 (p∗ , q ∗ )
x∗ = G(x∗ )
para todo punto (x, y) en un entorno de (p∗ , q ∗ ). Entonces la iteración de punto fijo
genera una sucesión {(pk , qk )}k≥0 convergente al punto fijo (p∗ , q ∗ ).
|x − 2| < 0.5
|x|
+ |y − 1| < 1
4
Inicializando en (p0 , q0 ) = (2, 0) no podemos garantizar la convergencia, pero aún ası́ converge
al punto fijo (1.9006767, 0.3112186).
ϕ(t) = 6e−t/2 + t − 2.
Usando técnicas de Cálculo Matemático conocidas, la solución general del PVI será de la forma
ϕ(t) = ce−t/2 + t − 2,
Notemos que como la solución depende de un parámetro, en la medida que varı́amos dicho
2
parámetro se van “moviendo” las curvas soluciones del PVI. Cuando traducimos esto a la
solución de un modelo matemático de un problema fı́sico, observamos que la condición inicial
0
nos fija la curva que relaciona una variable independiente (como el tiempo t), y otra dependiente
(como la posición o una temperatura). Esta condición inicial medida al comienzo de nuestro
experimento−2y que serı́a parte de la data de nuestro problema, determina lac =solución
−2 final que
c = −1
es la descripción final de nuestro modelo fı́sico. c=0
−4 c=1
c=2
6.6.2 Interpretación geométrica: el campo de direcciones c=3
c=4
−6
Suele ocurrir−1que a priori0 1
no conozcamos 2 3
la solución 4
de nuestro 5
problema 6
o bien no podamos
encontrarla por métodos clásicos conocidos. En ese caso, podemos ayudarnos a tener una
descripción de la solución usando el campo de direcciones.
Un campo de direcciones es una gráfica del plano cartesiano tal que en cada punto (t, y)
del rectángulo R = {(t, y) : a ≤ t ≤ b, c ≤ y ≤ d} representamos la pendiente de la solución
del PVI y = y(t) dado por la fórmula implı́tica m = f (t, y(t)). Es decir, para una cantidad
finita de valores mi,j = f (ti , yj ) representamos un vector pendiente a la solución que pasa por
el punto (ti , yj ) para i = 1, .., m y j = 1, ..., m.
4.5
campo de direcciones
3.5
2.5
1.5
En la Figura 6.32 se observa la gráfica del campo de direcciones del PVI original realizado
con 1la función CampoDirec.m en MATLAB.
0.5
Capı́tulo 6 2022 244
0
CAPÍTULO 6. MÉTODOS NUMÉRICOS
R = {(t, y) ∈ R2 : a ≤ t ≤ b, c ≤ y ≤ d}
Demostración. Sean y1 < y2 en el intervalo (c, d). Fijando t en el intervalo (a, b) y usando el
Teorema del Valor Medio, existe ξ con y1 < ξ < y2 tal que
∂f ∂f
|f (t, y1 ) − f (t, y2 )| = (t, ξ)(y1 − y2 ) = (t, ξ) |(y1 − y2 )| ≤ L|y1 − y2 |.
∂y ∂y
R = {(t, y) : a ≤ t ≤ b, c ≤ y ≤ d}.
Por tanto, podemos asegurar que nuestro PVI original tiene solución única. Ahora bien, aún
sabiendo que existe y es única, ¿cómo hacemos para hallarla en los casos en que no podamos
obtener una “expresión cerrada” la por métodos clásicos y conocidos del Análisis Matemático?
En ese caso, recurriremos a “aproximar” la solución.
{(tk , yk )} k = 1, . . . , n
que sean aproximaciones de la solución y = ϕ(t) evaluada en los puntos tk , esto es,
yk ≈ ϕ(tk ) k = 1, . . . , n.
Para construir dicho conjunto de puntos que verifiquen “aproximadamente” la EDO, primero
particionamos el intervalo [a, b] en n subintervalos iguales:
tk = a + kh k = 0, 1, . . . , n
y ′′ (ξ)
y(t) = y(t0 ) + y ′ (t0 )(t − t0 ) + (t − t0 )2
2!
donde ξ está entre t y t0 .
Evaluando esta expresión en t = t1 , siendo h = t1 − t0 el tamaño del paso:
y ′′ (ξ1 ) 2
y(t1 ) = y(t0 ) + hy ′ (t0 ) + h
2!
donde ξ1 está entre t1 y t0 .
y ′′ (ξ1 ) 2
y(t1 ) = y(t0 ) + hf (t0 , y(t0 )) + h.
2!
Entonces, si el tamaño del paso h es “suficientemente pequeño”, podemos despreciar el
último término del desarrollo de Taylor que tiene h2 , obteniendo ası́
y1 := y0 + hf (t0 , y0 ) ≈ y(t1 )
y2 := y1 + hf (t1 , y1 ) ≈ y(t2 )
..
.
yn := yn−1 + hf (tn−1 , yn−1 ) ≈ y(tn )
yk+1 = yk + hf (tk , yk ), k = 0, 1, . . . , n − 1,
1.8 n=6
n=12
n=24
1.6
exacta
1.4
1.2
0.8
ek = ϕ(tk ) − yk k = 0, 1, ..., n
Teorema 6.6.3. Sea y = ϕ(t) la solución del PVI. Si ϕ(t) ∈ C 2 [a, b] y {(tk , yk )}nk=0
son las aproximaciones generadas con el método de Euler, entonces:
Definición. El error final del intervalo [a, b] se conoce como error global final:
E(ϕ(b), h) ≈ Ch
h 1 1
E(ϕ(b), h/2) ≈ C ≈ Ch ≈ E(ϕ(b), h)
2 2 2
1
E(ϕ(b), h/2) ≈ E(ϕ(b), h)
2
Si reducimos a la mitad el tamaño del paso en el método de Euler, entoces cabe esperar que
el error final global se reduzca a la mitad.
Observamos la relación entre el tamaño del paso h definido por h = (b − a)/n y el error
final global |ϕ(b) − yn | para las aproximaciones del método de Euler
h n |ϕ(b) − yn | o(h) ≈ Ch
1.000000 3 0.294390 0.256000
0.500000 6 0.135455 0.128000
0.250000 12 0.065139 0.064000
0.125000 24 0.031961 0.032000
0.062500 48 0.015833 0.016000
0.031250 96 0.007880 0.008000
0.015625 192 0.003931 0.004000
2. Métodos multipasos
• Método de Adams-Bashforth-Moulton
• Método de Milne-Simpson
• Método de Hamming
Muchos de estos métodos pueden ser extendidos para aproximar soluciones de Sistemas
de EDO. También existen métodos especı́ficos para resolver Problemas de Contorno de
EDO. El campo es un campo de investigación, de desarrollo y de aplicaciones nuevas totalmente
abierto con nuevos problemas y desafı́os todo el tiempo.
Sin usar MATLAB, indique qué resultados se obtienen al ejecutar los siguientes comandos:
a = [ 1, 5, -1]
b = 3:-2:-2
c = [4 6 -1 0 -6 7 9 1 5]
d = [-1 1 5; 2 3 0; 1 5 -2]
Escriba un único comando para resolver las siguientes instrucciones:
(i) Obtener el máximo de cada fila de la matriz d.
(ii) Obtener la suma de los valores absolutos de todos los elementos de d.
(iii) Obtener la suma de los elementos de posición par del vector c.
(iv) Calcular el mı́nimo valor entre los elementos de posición 3,4,5 y 6 de c.
(v) Calcular la primer columna del producto externo de a por b.
(vi) Eliminar el cuarto elemento del vector c.
(vii) Generar un vector llamado z con los elementos de posición impar del vector c.
(viii) Agregar una cuarta columna a la matriz d cuyos elementos sean -2,3 y 0.
(ix) Armar una matriz diagonal con la diagonal de la matriz d.
(x) Calcular la norma del vector a-b.
3. (a) Escriba la función suma.m que reciba como argumentos dos vectores fila y devuelva
como resultado el vector columna suma.
(b) Modifique la función anterior de modo que los vectores de entrada sean filas o co-
lumna pero siga devolviendo un vector columna.
4. (a) Escriba una función llamada cuad.m que calcule las raı́ces de la ecuación cuadrática
ax2 + bx + c = 0 con a, b, c ∈ R, a 6= 0, recibiendo como argumentos de entrada los
coeficientes y devolviendo como salida las raı́ces:
√ √
−b + b2 − 4ac −b − b2 − 4ac
x1 = x2 = .
2a 2a
(b) Modifique dicha función para que se imprima por pantalla indicando si las raı́ces son
reales, reales iguales (en cuyo caso solo imprima una) o complejas conjugadas.
5. (a) Genere el script res_sel.m que lea una matriz cuadrada nxn A almacenada en el
archivo de datos matriz.dat y un vector nx1 b desde el rhs.dat.
(b) Resuelva el sistema de ecuaciones lineales Ax=b usando el operador aritmético barra
invertida \, y guarde el vector solución x en un archivo de datos de salida sol.dat.
6. (a) Escriba la función multi.m que reciba como argumentos dos funciones y un vector
de abscisas, y devuelva como resultado un vector con la multiplicación de ambas
funciones evaluadas en dichas abscisas.
(b) Calcule el producto de f1 y f2 del Ej. 7 en el vector x = [−1, 0, 1] usando multi.m.
Ploteos en MATLAB:
7. (a) Escriba m-archivos para evaluar cada una de las siguientes funciones:
(b) Grafı́quelas en el intervalo [−3, 3], todas en la misma gráfica y con distintos colores.
(c) Repita el item (b) pero graficando cada función en distintas gráficas de la misma
figura usando el comando subplot.
2. (a) Modifique el script ZoomPoli.m para que plotee el polinomio de sexto grado de la
siguiente manera algebraicamente equivalente a la ya implementada:
p(x) = (x − 1)6
(b) Modifique el script ZoomPoli.m para que plotee el polinomio cúbico de las siguientes
dos maneras algebraicamente equivalentes:
(a) Halle a mano las raı́ces de la ecuación con la resolvente utilizando aritmética de 4
dı́gitos y redondeo por truncamiento.
(b) Demuestre que las raı́ces de una ecuación general de 2do. grado ax2 + bx + c = 0
pueden obtenerse a través de las siguientes fórmulas:
−2c −2c
x1 = √ , x2 = √ .
b + b2 − 4ac b − b2 − 4ac
(c) Halle a mano las raı́ces de la ecuación con las fórmulas presentadas en el apartado
anterior, utilizando aritmética de 4 dı́gitos y redondeo por truncamiento.
(d) Compare errores relativos porcentuales de los resultados en los apartados a) y c).
¿Por qué se producen estas diferencias? Calcule valores exactos usando MATLAB.
(a) Utilice la función cuad.m del ejercicio 4(a) de la sección 6.7.1 para hallar sus raı́ces.
(b) Usando las fórmulas mejoradas del Ejercicio 2, implemente la función cuad_mejor.m
para calcular raı́ces.
(c) Compare los resultados obtenidos en a) y b) y obtenga conclusiones.
x x2 x3 xn
ex ≈ 1 + + + + ... + ,
1! 2! 3! n!
recibiendo como parámetros de entrada x y n (número de términos de la suma).
(b) Modifique la función anterior para que imprima en columnas el valor de las sumas
parciales, el término que se está sumando y el error absoluto cometido. Use la
función intrı́nseca exp de MATLAB para calcular el error.
(c) Corra los casos x = 1, 10, 100 cada uno para n = 10, 15, 20. Analice el error cometido.
6. (a) Calcule en forma exacta sen (π/2 + 2π10j ), con j entero positivo.
(b) Calcule en MATLAB la misma expresión, para j = 1, 10, 20, 50, 100, 1000.
(c) Intente dar una explicación a los resultados obtenidos.
7. Analice la función SerieSeno y las salidas por Ventana de Comandos en respuesta de:
(a) SerieSeno(pi/4,5e-9),
(b) SerieSeno(pi,5e-9),
(c) SerieSeno(5*pi,5e-9).
8. (a) Usando la función SerieSeno.m, grafique el error total cometido para x ∈ [0, π/2]
usando 6 y 15 términos no nulos del polinomio de Taylor del sen(x) respectivamente.
Considere que la función intrı́nseca sin de MATLAB devuelve el valor exacto del
seno. ¿Qué controla el error total en cada caso? Sugerencia: Vea las transparencias.
(b) Repita el item anterior para la función exponencial desarrollada en el Ejercicio 5.
9. Usando la máquina virtual que aparece en la Tabla de las transparencias con mantisa
binaria de 4 bits, 8 posibles exponentes y redondeo simétrico (si x1 y x2 son los números
de la Tabla más próximos al número x que quiero representar, entonces x redondea a x1
sólo si x < x1 +x
2
2
, caso contrario se redondea a x2 ), evalúe las siguientes expresiones:
(a) 1.1
(b) eps
(c) 20
(d) 8 + 5 ∗ eps
(e) (8 + 3 ∗ eps) + 3 ∗ eps
redondeo
Ejemplo: (5 + 3 ∗ eps) − 4 ∗ eps = (5 + 0.375) − 0.5 = 5,375 − 0.5 −−−−−→ 5.5 − 0.5 = 5.
¿Cuanto vale el epsilon máquina en nuestro caso?
10. Anticipe los resultados de evaluar las siguientes expresiones. Luego verifique en MATLAB.
(a) Resuelva a mano con el método de Gauss con pivoteo parcial escalado.
(b) Compruebe sus resultados utilizando la función Gauss.m en MATLAB.
(a) Resuelva a mano dicho sistema con eliminación Gaussiana. ¿Cuántas soluciones hay?
(b) Modifique la función Gauss.m para que imprima por pantalla la matriz aumentada
en cada paso y resuelva en máquina nuevamente el sistema.
(c) En la salida de pantalla, ¿resulta nula la última fila de la matriz ampliada? Explique
el error numérico cometido.
3. (a) Desarrolle en MATLAB un m-file que tome como argumento una matriz cuadrada
triangular superior A, el vector lado derecho b, y devuelva una solución x mediante
una sustitución regresiva (backward sustitutition).
(b) Repita lo anterior para una sustitución progresiva (forward sustitutition).
Sugerencia: Ayúdese con la función DescompLU.m.
(a) Genere un vector lado derecho bi en cada caso usando el comando ones de MATLAB
y resuelva el SEL Ai x = bi mediante algún método directo.
(b) Calcule el vector residual r = b − Ax̃ en los SEL anteriores y la norma del mismo
krkp para p = 1, 2, inf.
6. Para las matrices del ejercicio anterior, calcule con ayuda de MATLAB:
(a) El número de condición κ(A) = kAkp .kA−1 kp utilizando el comando inv(A) para
calcular la inversa y el comando norm(A,p) para la norma matricial con p = 1, 2, inf.
(b) Compare los resultados anteriores usando el comando cond(A) de MATLAB.
¿Están bien condicionadas las matrices?
(c) ¿Observa alguna relación entre el vector residual r y el número de condición κ(A)
de las matrices obtenido?
Métodos Iterativos:
7. Considere el SEL
3x1 − x2 + x3 = 1
3x1 + 6x2 + 2x3 = 0
3x1 + 3x2 + 7x3 = 4
(a) Calcule a mano las dos primeras iteraciones del método de Jacobi.
(b) Resuelva en MATLAB ulilizando la función Jacobi.m, tomando como aproximación
inicial x(0) = 0 y una tolerancia de 1 × 10−3 .
(c) ¿Puede asegurar la convergencia del método de Jacobi en este caso? Justifique.
(d) Repita los 3 items anteriores utilizando el método de Gauss-Seidel.
8. (a) Compare las implementaciones en MATLAB de Jacobi.m y GaussSeidel.m
(b) Modifique ambas para que se imprima por pantalla la iteración k-ésima, la aproxi-
mación, el vector error (iter, x, err).
9. Compare la velocidad de convergencia en la resolución de los siguientes SELs,
(
−x + 3y = 1
x + z = 2
5x − y + z = 10
(a)
6x − 2y = 2 (b) −x + y = 0 (c) 2x + 8y − z = 11
x + 2y − 3z = 0 −x + y + 4z = 3
(a) Plotee dicho polinomio en el intervalo [1, 2]. ¿Cuantas raı́ces observa?
(b) Plotee dicho polinomio en el intervalo [1.4, 1.7]. ¿Y ahora?
(c) Aproxime dichas raı́ces utilizando NewtonRapshon.m y compare los errores absolutos.
Sug. Tome como valor exacto las raı́ces obtenidas usando la función roots.
p
3. Considere la función s(x) = sign(x − 2) |x − 2| para x ∈ [0.4].
(a) Modifique NewtonRaphson.m para que imprima otra columna con el cociente |e|en+1 n|
|
6. Determine analı́ticamente los puntos fijos de cada una de las siguientes generatrices:
g1 (x, y) = x − y 2 g1 (x, y) = sen(y)
a) b)
g2 (x, y) = −x + 6y g2 (x, y) = −6x + y
(a) Usando la aproximación inicial (p0 , q0 ) = (1.1, 2.0), calcule dos iteraciones a mano
mediante el esquema de Punto Fijo.
(b) Realice lo mismo pero utilizando el esquema iterativo de PuntoFijo.m.
(c) Grafique las curvas involucradas en el plano z = 0.
2. Usando la función Euler.m, resuelva los siguientes PVI. En una misma figura, grafique
las aproximaciones para distintos pasos y las soluciones exactas:
′
y = −ty, t ∈ [0, 5]
Solución Exacta: ϕ(t) = e−t /2 .
2
(a)
y(0) = 1
′
2 − 2y , t ∈ [0, 10]
1 2
y = 1+t t
(b) Solución Exacta: ϕ(t) = 1+t2.
y(0) = 0
′ 1
y = 4 y 1 − 201
y , t ∈ [0, 20] 20
(c) Solución Exacta: ϕ(t) = 1+19e −t/4 .
y(0) = 1
[Alav2006] J.Alavez Ramı́rez, ”Métodos Numéricos I”, Universidad Juárez Autónoma de Tabasco,
México, 2006.
[Apos1984] T.M. Apostol, “Calculus”, Vol. 1, Ed. Reverté, Buenos Aires, 1984.
[Atki1985] K. Atkinson, “Elementary Numerical Analysis”, John Wiley & Sons, New York, 1985.
[BrCh2004] J.W. Brown, R.V. Churchill, “Variable compleja y aplicaciones”, McGraw-Hill, Madrid,
2004.
[Codd1968] E.A. Coddintgton, “Introducción a las Ecuaciones Diferenciales Ordinarias”, Compañia
Editorial Continental, S.A., México, 1968.
[FuGi1995] R. Fuster, I. Giménez, “Variable Compleja y Ecuaciones Diferenciales”, Reverté, Barce-
lona, 1995.
[GoVa2013] G.H. Golub, C.F. Van Loan, “Matrix Computations”, 4th Ed., The Johns Hopkins Uni-
versity Press, Baltimore, 2013.
[Jame2002] G. James, “Matemáticas avanzadas para Ingenierı́a”, Pearson Educación, México, 2002.
[KiCh1994] D. Kincaid, W. Cheney, “Análisis Numérico. Las matemáticas del cálculo cientı́fico”,
Addison-Wesley Iberoamericana, Delaware, 1994.
[Krey2003] E. Kreyszig, “Matemáticas avanzadas para Ingenierı́a”, Limusa Wiley, México, 2003.
[MaFi2000] J.M. Mathews, K.D. Fink, “Métodos Numéricos con MATLAB”, 3era Ed., Pearson Pren-
tice Hall, Madrid, 2000.
[Mole2015] C.B. Moler, ”Numerical computing with MATLAB”, SIAM, Philadelphia, 2004.
[Saad2000] Y. Saad, ”Iterative Methods for Sparse Linear Systems”, SIAM, Philadelphia, 2003.
[SiSo2010] J.W. Signorelli, J. Sorribas, ”Notas de Clases de Informática Aplicada y Métodos Compu-
tacionales”, FCEIA, UNR, Rosario, 2010.
[SLSS2011] M. Spiegel, M. Lipschutz, J. Schiller, D. Spellman, “Variable Compleja y Ecuaciones
Diferenciales”, Mc Graw Hill, México, 2011.
[TrBa1997] Ll.N. Trefethen, D. Bau, ”Numerical Linear Algebra”, SIAM, Philadelphia, 1997.
[Wuns1997] A.D. Wunsch, “Variable Compleja con Aplicaciones”, Addison Wesley, México, 1997.
• https://la.mathworks.com/products/matlab.html
• https://www.gnu.org/software/octave/
• https://www.scilab.org/
• www.fortran.com/