Libro Señales PDF
Libro Señales PDF
Libro Señales PDF
Seales y Sistemas
Alfredo Restrepo Palacios, Ph.D
Bogot, Agosto 13 2015
ndice general
0. Bases matemticas
0.1. INTRODUCCIN . . . . . . . . . . . . . . . . .
0.2. Relaciones y Funciones . . . . . . . . . . . . . .
0.2.1. Par ordenado . . . . . . . . . . . . . . .
0.2.2. Producto cartesiano . . . . . . . . . . .
0.2.3. Relacin . . . . . . . . . . . . . . . . . .
0.2.4. Particin . . . . . . . . . . . . . . . . . .
0.2.5. Relaciones de orden . . . . . . . . . . .
0.2.6. Suprmum, Infmum . . . . . . . . . . .
0.2.7. Mximo y mnimo . . . . . . . . . . . .
0.3. Funciones . . . . . . . . . . . . . . . . . . . . .
0.3.1. Sobreyeccin . . . . . . . . . . . . . . .
0.3.2. Inyeccin. (funcin uno-a-uno) . . . . .
0.3.3. Biyeccin. (correspondencia uno-a-uno)
0.4. Relaciones Ternarias . . . . . . . . . . . . . . .
0.4.1. Propiedades de relaciones ternarias . .
0.5. Estructuras algebricas . . . . . . . . . . . . . .
0.5.1. Grupos . . . . . . . . . . . . . . . . . . .
0.5.2. Anillos . . . . . . . . . . . . . . . . . . .
ii
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
4
4
6
6
6
7
11
12
12
13
13
15
15
15
18
19
20
22
25
NDICE GENERAL
iii
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
27
28
29
29
29
30
31
32
33
33
33
33
34
35
35
35
36
37
37
38
39
42
43
44
45
47
48
49
51
iv
NDICE GENERAL
51
51
52
53
54
54
56
57
58
58
59
59
59
59
60
64
65
65
66
69
69
71
72
74
74
75
77
78
NDICE GENERAL
0.15.1. Definicin . . . . . . . . . . . . . . . . . . . . . . . . . . .
78
79
80
81
. . . . . . . . . . . . . .
3
82
84
0.16. Sucesiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
85
0.16.1. Subsucesiones . . . . . . . . . . . . . . . . . . . . . . . . .
87
88
0.17.1. Bolas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
89
. . . . . . . . . . . . . . . . . . . . . . . .
89
90
90
91
91
92
94
0.18. Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
97
98
99
vi
NDICE GENERAL
108
109
109
110
110
112
113
114
114
116
119
119
119
120
121
122
129
129
131
131
132
133
133
134
1. Seales y sistemas
152
1.0. INTRODUCCIN . . . . . . . . . . . . . . . . . . . . . . . . . . . 152
NDICE GENERAL
vii
166
viii
NDICE GENERAL
185
185
190
193
197
197
199
200
202
203
204
205
205
209
211
211
213
213
213
213
215
216
217
NDICE GENERAL
3.1.
3.2.
3.3.
3.4.
3.5.
3.6.
3.7.
3.8.
ix
NDICE GENERAL
NDICE GENERAL
xi
334
336
336
337
357
357
361
361
364
367
368
373
375
375
376
378
380
380
382
384
384
385
387
390
391
NDICE GENERAL
xii
5.11.
5.12.
5.13.
5.14.
5.15.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
392
395
398
400
402
404
406
408
415
415
417
419
420
421
421
422
422
423
423
424
424
425
426
427
428
430
432
NDICE GENERAL
xiii
432
432
433
433
435
438
438
438
438
441
443
444
444
445
xiv
NDICE GENERAL
Notacin
A continuacin se da el significado de parte de la terminologa matemtica
que se usa. Se asume que se tiene cierta familiaridad con lgica matemtica
as como con la teora de conjuntos: en particular, se asume que conoce el
significado de la conjuncin , de la disyuncin , de la implicacin y de
la doble implicacin . Las letras son smbolos de elementos, conjuntos, etc.
x A : x pertenece a A (o, x es miembro de A).
P(x) : x cumple con la propiedad P.
x A : para algn (o existe un) elemento x en el conjunto A
x A : para cada (o para todo) elemento x del conjunto A
a = b : quiere decir que los smbolos a y b representan la misma cosa.
: tal que
B = {x A : Q(x)} : B es el conjunto de los elementos de A que cumplen con la
propiedad Q.
A B : el conjunto A contiene al conjunto B.
A B = {x : x A x B}.
A B = {x : x A x B}.
S
Ba = {x : a A, x Ba }.
aA
T
Ba = {x : a A, x Ba }.
aA
Conjunto vaco ; x, x 3 A A = .
NDICE GENERAL
NDICE GENERAL
, : xi
o, O : omicron
, $, : pi
, P : ro
, , : sigma
, T : tau
, : upsilon
, , : phi f i
, X : ji
, : psi
, : omega
Captulo 0
Bases matemticas
0.1.
INTRODUCCIN
0.1. INTRODUCCIN
i {a, b} = {b, a}
ii {1, 2, 2} = {1, 2}
Entre las personas que en tiempos modernos relacion la lgica y las
matemticas est el lgico italiano Giuseppe Peano (1858-1932). Peano dio
una lista de propiedades que cumplen los nmeros naturales. En Pars en 1900
durante la segunda conferencia internacional de matemticas Pean entreg a
B Rusell una copia de su Formulario donde trataba la axiomatizacin de las
matemticas. La teora de conjuntos, originada por Georg Cantor1 , se usa como
base para las matemticas hoy en da. La idea de basar las matemticas en la
teora de conjuntos y, a la larga en lgica, tom as fuerza desde finales del
siglo XIX. Los resultados de Gdel (1906-1978) indicando la imposibilidad de
obtener sistemas formales completos en los cuales basar las matemticas, quizs
le restaron algo de mpetu a esta tendencia racionalista en las matemticas de
comienzos de siglo. Poincar sostuvo la importancia de la intuicin en las
matemticas a comienzos del siglo XX.
1 Georg
0.2.
Relaciones y Funciones
De los varios conceptos que se definen en este captulo, el lector seguramente tendr una idea intuitiva clara. Sin embargo, en aras de anclar la
intuicin con el mtodo axiomtico, los definimos en todo caso. As el caso de
par ordenado, que definimos de dos maneras dentro de la teora de conjuntos.
0.2.1.
Par ordenado
0.2.2.
Producto cartesiano
Podemos generalizar la nocin de tupla y definir las n-tuplas para cada entero
positivo n.
Ejercicio. Defina 4-tuplas.
0.2.3.
Relacin
A = {0, 1, 2}
A A = {(0, 0), (0, 1), (0, 2), (1, 0), (1, 1), (1, 2), (2, 0), (2, 1), (2, 2)}
Sea R = {(1, 0), (2, 1)} decimos entonces que 1 R 0 y 2 R 1. (Se lee: 1
est relacionado con 0 y 2 est relacionado con 1, o R relaciona a uno con cero,
R relaciona a dos con uno; tambin, uno erre cero, dos erre uno y uno no
erre dos.)
Ejemplo.
Propiedades de relaciones
Sea R una relacin entre los elementos del conjunto A.
i) R es una relacin reflexiva si
a A, (a, a) R
Figura 4: Transitividad. Si est relacionada con y con entonces est relacionada con
Con respecto a la figura 4, note que (a, b) y (b, c) determinan dos vrtices
opuestos de un rectngulo y que uno de los vrtices restantes es (b, b)
que est sobre la diagonal, el vrtice restante, (a, c) es el requerido por
la propiedad de transitividad.
Ejemplo. es transitiva y reflexiva. A B, B C A C, A A
10
Ejemplo. no es transitiva. A B, B C ; A C
n
P
m=0
11
(nm ) = 2
conjunto de n elementos?
Relaciones de equivalencia y clases de equivalencia
La relacin de equivalencia en muy til en matemticas. Decimos que una
relacin es una relacin de equivalencia si es reflexiva, simtrica y transitiva. Sea E una relacin de equivalencia entre los elementos de un conjunto
A. Para cada elemento a de A definimos su clase de equivalencia como el
conjunto de los elementos de A que son equivalentes a a y se denota [a] :
[a] = x a : (x, a) E
Ejercicio. Muestre que b [a] si y slo si la [a] = [b].
Ejercicio. Muestre que, dados (a, b) A, [b] = [a], o, [a] [b] =
Ejercicio. Muestre que A =
0.2.4.
[a]
aA
Particin
12
lencia entre los elementos de A, tal que las clases de equivalencia resultantes
sean los elementos de la particin dada.
0.2.5.
Relaciones de orden
R es un orden parcial si es una relacin transitiva, antisimtrica y reflexiva. R es un orden lineal si es un orden parcial y es una relacin conectada.
El orden parcial tpico es el dado por ser subconjunto de entre el conjunto
de subconjuntos de un conjunto dado, por ejemplo en la figura 6 se muestra el
orden , para los subconjuntos de {0, 1, 2}
0.2.6.
Suprmum, Infmum
Sea A un conjunto linealmente ordenado por la relacin . Sea B un subconjunto de A y sea q un elemento de A. Si b B, se tiene b q, se dice que q
es una cota superior de B. q es una cota superior mnima (suprmum) de B si:
0.3. FUNCIONES
13
0.2.7.
Mximo y mnimo
Sea A un conjunto linealmente ordenado por la relacin . Sea B un subconjunto de A. Si B tiene suprmum r que pertenece a B, decimos que r es el
mximo de B. Similarmente, si B tiene infmum t, y t pertenece a B, decimos
que t es el mnimo de B.
Ejemplo. En los reales, en el intervalo abierto (0, 1) no tiene ni mximo ni
0.3.
Funciones
14
0.3. FUNCIONES
15
i) f (C D) = f (C) f (D)
ii) f (C D) = f (C) f (D)
Ejercicio. Si C y D son subconjuntos del rango de ( f, A, B), demuestre o d un
contraejemplo:
i) f 1 (C D) = f 1 C) f 1 (D)
ii) f 1 (C D) = f 1 (C) f 1 (D)
0.3.1.
Sobreyeccin
0.3.2.
0.3.3.
16
0.3. FUNCIONES
17
en cuantificar las variables que modela y en disear con base en stas. As, el
uso de nmeros se convierte en una segunda naturaleza para el ingeniero. En
las secciones siguientes se dan definiciones de los nmeros naturales, enteros,
racionales, reales, complejos y de los cuaternios. Dado que en ingeniera las
magnitudes se cuantifican, se podra argumentar que es importante tener una
idea clara de cada uno de estos tipos de nmeros. Los naturales se definen con
base en nociones de la teora de conjuntos, los enteros con base en los naturales,
los racionales con base en los enteros, los reales con base en los racionales. En
cada uno de estos casos tenemos una relacin de orden lineal y las operaciones
de suma y producto. Luego, tenemos RN , en donde tenemos la operacin de
suma. En R2 tenemos adems el producto complejo, y en R4 el producto de
cuaternios.
Definicin 3. Una biyeccin p : A A de un conjunto finito, a s mismo, se conoce
como una permutacin de ste.
Definicin 4. Una funcin t : A A de un conjunto finito A, con al menos dos
elementos, a s mismo, se conoce como una transposicin si A tiene dos elementos
a1 y a2 para los que t(a1 ) = a2 y t(a2 ) = a1 y, para el resto de los elementos de A,
f () = .
Ejercicio. Muestre que si una permutacin se descompone como una composi-
18
0.4.
Relaciones Ternarias
Una relacin ternaria entre los elementos de un conjunto A es un subconjunto del producto cartesiano A3 .
Ejemplo. Para el conjunto A = {0, 1, 2}, donde
AAA=
{(0, 0, 0), (0, 0, 1), (0, 0, 2), (0, 1, 0),
(0, 1, 1), (0, 1, 2), (0, 2, 0), (0, 2, 1),
(0, 2, 2), (1, 0, 0), (1, 0, 1), (1, 0, 2),
(1, 1, 0), (1, 1, 1), (1, 1, 2), (2, 0, 0),
(2, 0, 1), (2, 0, 2), (2, 1, 0), (2, 1, 1),
(2, 1, 2), (2, 2, 0), (2, 2, 1), (2, 2, 2)}
considere la relacin R = {(0, 1, 2), (1, 2, 0), (2, 0, 1)}
0.4.1.
19
20
0.5.
Estructuras algebricas
Dado que esta seccin es un resumen del tema que se trata, debe usarse
como una referencia o un diccionario, ms que como un texto sobre el tema.
El concepto de operacin, tal como se usa en aritmtica, lo podemos abstraer y generalizar, considerando simplemente la estructura consistente de un
conjunto en el que se define una operacin entre sus elementos; as, se hacen
explcitas ciertas propiedades generales que pueden tener estas estructuras.
Las operacin se asume binaria, es decir que a cada par de elementos
del conjunto le corresponde un elemento del conjunto que es el resultado de
ejecutar la operacin con los dos elementos dados; as, la operacin es tambin
clausurativa.
Las transformaciones lineales juegan una papel importante en ingeniera.
Los conceptos de superposicin y de respuesta proporcional son intuitivos y
muy tiles. Qu es una transformacin lineal? Para responder a esta pregunta
formalmente necesitamos el concepto de espacio vectorial y, aunque el lector
probablemente haya trabajado en repetidas ocasiones con espacios vectoriales, para una definicin formal, debemos considerar antes los conceptos ms
bsicos de grupo y anillo. Estas y otras estructuras resultan de considerar diversas posibilidades; por una parte, consideramos las posibilidad de que la
operacin sea asociativa o no, que haya un elemento identidad o no, que existan inversos o no, que la operacin sea conmmutativa o no, que haya ms de
una operacin, que una operacin distribuya con respecto a la otra o no, etc. etc.
Sean A un conjunto y o una operacin o: A A A entre los elementos
de A. Note que as una operacin es una funcin que asigna un elemento a
cada par ordenado de elementos del conjunto. Cuando el conjunto A es finito
y pequeo, puede ser conveniente hacer explcita la operacin tabulando los
posibles resultados en una tabla.
As, decimos que un grupoide (A, o) est dado por un conjunto no vaco
a
b
c
d
a
a
c
d
b
21
b
d
b
a
c
c
b
a
c
a
d
c
d
b
d
Tabla 1: La tabla de una posible operacin entre los elementos del conjunto {a, b, c, d}
A y una operacin o entre los elementos de ste. Note que en principio, no
requerimos que la operacin sea asociativa.
Ejemplo. Considere el conjunto {a, b, c, d} y la operacin en la Tabla 1
a (b c) = (a b) c
En cierto sentido, la asociatividad nos dice que la operacin, en principio binaria, pueda ser aplicada a un nmero arbitrario de operandos. En un
semigrupo no requerimos la existencia de un elemento identidad.
Ejemplo. (N, +) es un semigrupo.
Ejercicio. D un ejemplo de una operacin no asociativa.(Loas octoniones)
Ejercicio. Diga si la operacion del grupoide en la tabla 2 es asociativa.
Un monoide (mdulo) es un semigrupo en el cual hay un elemento identidad. Es decir, un elemento que muchas veces denotamos como e tal que:
a A a e = e a = a. Al elemento identidad tambin se le conoce como el
mdulo de la operacin. Por ejemplo, el cero es el mdulo de la suma y el uno
es el mdulo del producto.
22
Tabla 2: La tabla de una operacin entre los elementos del conjunto finito , , , ,
0.5.1.
Grupos
23
eg =
ee =
fe =
ff =
ge =
gg =
20
3
, f = ej
2
3
, g = ej
4
3
Galois (1811-1933) fue un matemtico francs reconocido, junto con Abel, por su
estudio de las races de los polinomios y la idea de cuerpo (field).
3 Niels Henrick Abel (1802-1829) matemtico Noruego, demostr que no existe una operacin
algebrica explcita con forma de races, productos, cocientes, sumas y restas, para encontrar las
races de un polinomio de orden 5.
4 Tanto Galois como Abel murieron bastante jvenes, Galois a causa de una herida recibida en
un duelo y Abel a causa de la tuberculosis.
24
Homorfismos de grupos
Dados dos grupos (G, ) y (F, ), una funcin : G F que respete las
operacines de los grupos, en el sentido que (g1 g2 ) = (g1 ) (g2 ), se
conoce como un homomorfismo. Un isomorfismo es un homomorfismo que es
tambin una biyeccin.
El espacio cociente G/H. Dado un subgrupo H del grupo G, el conjunto de
los subconjuntos de la forma gH := {gh : h H}, g G, de G se denota como
G/H. El espacio cociente es un grupo con respecto a la operacin g1 H g2 H :=
(g1 g2 )H si y slo si H es normal, es decir, si Hg = gH para cada g en G
(en donde Hg := {hg : h H}). La razn es que, dado que (g1 g2 )1 =
g1
g1
, queremos que (gH)1 = Hg1 y si g1 H , Hg1 = (gH)1 entonces
2
1
(gH) (gH)1 , (g g1 )H. Para un subgrupo normal H, gHg1 = H.
K es un conjunto de generadores del grupo G, si el subgrupo ms pequeo de G que contiene a K es G mismo.
1
1
i
i
j
j
k
k
1
1
1
i
i
j
j
k
k
1
1
1
i
i
j
j
k
k
i
i
i
1
1
k
k
j
j
i
i
i
1
1
k
k
j
j
25
j
j
j
k
k
1
1
i
i
j
j
j
k
k
1
1
i
i
k
k
k
j
j
i
i
1
1
k
k
k
j
j
i
i
1
1
0.5.2.
Anillos
26
Un anillo de divisin no conmutativo se conoce como un cuerpo de inclinacin (skew field). (Ej. los cuaternios.)
Ejercicio. Demuestre que todo dominio de enteros finito es un cuerpo.
Ideales
Cuando se tiene una relacin de equivalencia tipo mdulo-N (por ejemplo
los enteros, mdulo 2) entre los elementos de un anillo A, denotamos la clase
de equivalencia que contiene al cero (mdulo aditivo) como la clase NA (en el
ejemplo, 2Z = {0, 2, 4, ...}). Al escoger un elemento de cada una de las clases
de equivalencia obtenemos otro anillo, que denotamos A/NA (en el ejemplo,
Z/2Z = {2Z, 2Z + 1} o, abreviadamente, = {0, 1}).
27
0.5.3.
Espacio vectorial
28
0.5.4.
Matrices y Determinantes
11
21
.
.
.
n1
12
22
...
...
n2
...
1m
2m
nm
X
p
(1)(p)
m
Y
i=1
i.p(i)
29
0.5.5.
Homomorfismos
Cuando se tiene una funcin, cuyo dominio y rango cuentan con alguna
estructura algebrica y sta es respetada por la funcin, se dice que se tiene un
homomorfismo; si adems se tiene que la funcin es una biyeccin, se dice que se
trata de un isomorfismo y que los espacios dominio y rango, son algebraicamente
equivalentes o que son isomorfos.
0.5.6.
Algebra de conjuntos
Considere una coleccin de subconjuntos de un conjunto dado, clausurativa con respecto a la toma de uniones (operacin suma), intersecciones (operacin producto) y complementos. Note que la unin distribuye con respecto
a la interseccin y viceversa. Este tipo de estructura algebrica se conoce como
un lgebra booleana con complemento. Ac se tienen las leyes de De Morgan.
Alternativamente, considere las operaciones de suma y producto mdulo-2
sobre el conjunto {0, 1}.
0.6.
Bertrand Russell escribi una vez que el nfasis que puso Kant en la intuicin y en lo apriori se debi a que algunas de las demostraciones de Euclides
requeran de las figuras. Ms adelante, Hilbert dio una axiomatizacin de la
geometra en la que los axiomas son suficientes y no se requiere del uso de
figuras. El teorema de la incompletitud de Godel nuevamente puso en alguna
medida en tela de juicio al mtodo axiomtico. En la prctica, el matemtico
30
usa tanto su intuicin como el mtodo axiomtico. La corriente axiomatizadora del s.XIX produjo las definiciones de los nmeros naturales y reales que se
dan a continuacin, basadas en las teora de conjuntos de Cantor.
0.7.
Aunque todos tenemos una idea intuitiva de los nmeros naturales, ellos
se pueden definir axiomticamente, con base en la teora de conjuntos.
Los nmeros naturales se pueden definir en forma inductiva, as: Sea N el
conjunto que:
i) contiene al conjunto vaco como elemento: N.
ii) si x N, entonces x {x} N.
Nota. x {x} se conoce como el sucesor de x y lo denotamos S(x).
iii) x N S(x) , .
iv) El nico elemento de N que no es el sucesor de algn elemento de N es
.
v) x, y N
n , y S(x) , S(y).
de conjuntos.
A los elementos de N los llamamos nmeros naturales. As, los siguientes
4 conjuntos son nmeros naturales: {}, {{}}, {{}, {{}}}, {{}, {{}}, {{}, {{}}}} y los denotamos:
0 := = {}
1 := {} = {{}}
2 := {0, 1} = {{}, {{}}}
3 := {0, 1, 2} = {{}, {{}}, {{}, {{}}}} note que S(n) = {0, 1, 2, . . . , n}
31
0.7.1.
32
0.7.2.
Con base en el orden lineal que acabamos de definir, podemos dar un orden
circular para cada uno de los conjuntos n N. Para n = {0, 1, 2, . . . , n 1} El
orden consiste precisamente de aquellas triplas (a, b, c) para las que a < b < c,
c < a < b, b < c < a. Por ejemplo, para 4 = {0, 1, 2, 3}, el orden circular est
dado por
C = {(0, 1, 2), (1, 2, 0), (2, 0, 1); (0, 1, 3), (1, 3, 0), (3, 0, 1); (0, 2, 3), (2, 3, 0), (3, 0, 2);
(1, 2, 3), (2, 3, 1), (3, 1, 2)}.
Ejercicio. Muestre que el nmero de elementos del orden circular de n tiene
n
3
elementos.
0.7.3.
33
Divisor y mltiplo
0.7.4.
Nmeros Primos
0.7.5.
Factor
primos de cada nmero natural n, n = pk1 pk2 pkt , p1 < p2 < < pt es nica.
(Gauss)
[[Elaborar]]
Ejercicio. Muestre que hay infinitos nmeros primos.
0.7.6.
34
0.7.7.
Figura 12: Una inyeccin del conjunto de los nmeros pares al conjunto de los nmeros naturales
35
biyeccin de N a A.
0.7.8.
Conteos
0.7.9.
Cardinalidad
0.8.
36
El conjunto de los nmeros enteros se puede definir en la siguiente forma; inicialmente consideramos el conjunto {(x, y) : x {0, 1}, y N {0}} =
{0, 1} N {0}, en donde tenemos dos copias de los nmeros naturales; una de
los enteros de la forma (0, n) que llamamos positivos, y la otra de enteros no
positivos (1, n); luego, identificamos los elementos (0, 0) y (1, 0), de tal forma
que slo haya un cero entero que, sin embargo, no tendr signo (en la tradicin anglosajona, que seguimos aqu, el cero no es ni positivo ni negativo,
mientras que en la tradicin germano-francesa, el cero es positivo y negativo).
Formalmente, definimos el conjunto de los nmeros enteros as:
Z := {(0, y) : y N {0}} (1, y) : y N {0}} {}
donde = {(0, 0), (1, 0)}. Esta definicin y la notacin que comnmente usamos
para los nmeros enteros se relacionan as:
+n=n= (0, n), n > 0, (nmeros positivos)
n= (1, n), n < 0, (nmeros negativos)
0= .
Por convencin, el nmero entero cero no es ni positivo, ni negativo.
Esta convencin la mantendremos para el cero racional y para el cero real.
Similarmente, en el caso de los nmeros complejos resulta contraproducente
definirle ngulo al origen del plano complejo, o cero complejo, y para usar una
sola convencin en todos los casos hemos decidido no definir signo para el
cero en ningn caso. (Como dijimos ya, a diferencia de la convencin adoptada aqu, algunos
matemticos usan la convencin que el cero es negativo y positivo)
0.8.1.
37
{((x1 , x2 ), (y1 , y2 )) : x1 = 1 y1 = 0 ,, x1 = 1 y1 = 1 y2 x2 ,, x1 = 0 y1 =
0 x2 y2 }
Ejercicio. Muestre que aunque los enteros son contables, no se pueden contar
0.8.2.
Intervalos enteros
El intervalo entero /n, m/, donde n y m son nmeros enteros, es el conjunto de nmeros enteros mayores o iguales que n y menores o iguales que
m.
Ejemplo. /5, 9/ = {5, 6, 7, 8, 9}
Ejemplo. /9, 5/ =
Ejercicio. Muestre que (Z, +, ) es un grupo.
Ejercicio. Es (Z, +, ) un anillo?.
0.8.3.
Algoritmo de la divisin
Ejemplo. J5K3 = 2,
J1K3 = 2,
J1KN = N 1
[[algoritmo de Euclides]]
38
0.8.4.
Congruencias mdulo N
K.F. Gauss6 fue tal vez la persona que reconoci por primera vez la utilidad
y la generalidad de este tipo de relacin.
Considere el conjunto Z de los nmeros enteros; sean m y n nmeros enteros
y sea N un entero positivo; se dice que los nmeros m y n son congruentes
mdulo N, si su diferencia es un mltiplo entero de N, es decir, si r Z m
n = rN.
Ejercicio. Muestre que la relacin de congruencia mdulo 3 es una relacin de
0
r1 < 3 y n = 3c1 + 21 , 0 r2 < 3. Note que | r1 r2 |< 3 y muestre que: hay
una k entera tal que m n = 3k, si y solo si r1 = r2
Ejercicio. Muestre que los nmeros enteros m y n son congruentes mdulo N
si y slo si los residuos (segn el algoritmo de la divisin) al dividirlos por N
son iguales.
Ejercicio. Muestre que si m es congruente con n, mdulo N, y que si a es
0.9.
39
El trmino razn tiene entre otros sentidos los de proporcin y de explicacin as tambin los trminos griego oo y el latino ratio. La relacin entre
dos nmeros es un concepto al que los pitagricos dieron gran importancia y
lo que inicialmente tena como sentido razn de ser, pas a tener tambin en de
cociente. As, en nmero racional es una relacin entre dos enteros. El trmino racional tiene pues los dos sentidos mencionados anteriormente desde el
nacimiento de la cultura occidental, y no es raro que la implicacin del teorema
de Pitgoras de la existencia de longitudes irracionales, o inconmensurables,
fuera causa de incomodidad.
El conjunto de los nmeros racionales es el tercer conjunto de nmeros
que definimos. Como veremos, tambin es un conjunto contable (lo cual podra sorprendernos inicialmente). Desde ciertos puntos de vista, es similar al
conjunto de los nmeros reales que definiremos en la siguiente seccin; sin
embargo, los racionales con respecto a los reales adolecen de ciertas fallas,
geomtrica, topolgica, algebrica y analtica. Sin embargo, en su definicin,
son ms intuitivos que los reales.
Ejercicio. Encuentre un conteo del conjunto Z Z
40
no tiene suprmum.
p
CASO 1. Si ( q )2 < 2 entonces existe un nmero entero tal que,
2
p 2 p + ( 1 )
< 2
( ) <
q
q
entonces x =
p1
q
>
p
q
y x B. Como,
2
!2
p + ( 1 )
p
=
+
q
q
2p
+
q2
y queremos
(p + 1 )
q2
p2
p2
<2 2
2
q
q
1
2
41
1
2
q2
de donde
p2
q2
<2
2p
1
+ 2 < 2q2 p2
2p + 1 < (2q2 p2 )2
2
p + ( 1 )
p2
<
2 <
q
q2
entonces x =
p1
q
es el infmum de B. Como,
1 2
p2
= 2
q
2p
q2
y queremos
p2 p
q2 q
1 2
p2
< 2 2
q
q2
1
2
<
p2
2
q2
42
0.9.1.
Decimos que el nmero racional [(x, y)] es menor o igual que el nmero
racional [(m, n)], o que [(x, y)] [(m, n)] si,
1. 0 m y x 0
2. x 0, m 0, xn ym
3. x 0, m 0, ym xn
Note que este orden se puede interpretar con base en el dibujo de la figura
14. Si una lnea q est despus que otra p movindose en sentido horario, el
nmero q es menor que el nmero p. Note que dadas dos lneas p y q, se puede
encontrar un punto de Z (N {0}) entre las dos lneas (entre dos racionales
siempre hay otro racional) y por lo tanto, hay otra lnea que contiene el punto
mencionado y que est entre p y q segn el orden definido.
Ejercicio. Dados los nmeros racionales p y q, encuentre un nmero racional r
43
0.9.2.
44
0.9.3.
45
0.9.4.
Cortes de Dedekind
46
Dedekind. Este corte se usa para definir el nmero real raz de dos, en la seccin
siguiente.
9 Bernard Bolzano (Praga 1781, 1848). Matemtico, lgico, filsofo , telogo y sacerdote catlico,
con ancestros italianos. Uno de los primeros matemticos que en tiempos modernos abog por
el rigor en el anlisis matemtico, con sus 3 pobras Beytrage zu einer begrundeteren Darstellung
der Mathematik (1810), Der binomische Lehrsatz (1816) y Rein analytischer Beweis (1817).
10 Bertrand Russell (1872-1970) Matemtico y filsofo ingls se le conoce como el fundador de la
lgica moderna.
11 Georg Riemann (1826-1866) Matemtico alemn conocido por sus desarrollos en geometra no
eucldea, teora analtica de nmeros y por sentar las bases para la topologa moderna.
12 Karl Weierstrass (1815-1897) Matemtico alemn, tutor de Cantor y de Engel, entre otros.
Aunque no public ningn artculo, s escribi 3.
0.10. TOPOLOGA
47
0.10.
Topologa
48
Definicin 15. Un subconjunto X de un conjunto Y es denso si X = Y. Equivalentemente, todo conjunto abierto, no vaco, de Y contiene al menos un punto de
X.
Definicin 16. Un subconjunto X de un espacio C es denso en ninguna parte
= .
si su clausura tiene interior vaco: Int(A)
Definicin 17. Se dice que un espacio topolgico X es de Hausdorff si dados dos
puntos cualesquiera, hay dos conjuntos abiertos disyuntos conteniendo uno de los
puntos cada uno.
Definicin 18. Una base para un espacio topolgico es una coleccin de conjuntos
abiertos tal que cada conjunto abierto en el espacio es una unin de elementos de la
base.
Ejercicio. Suponga que |Bd(A)| 1 entonces A es enumerable y denso en ningu-
na parte.
Ejercicio. Demuestre o d un contraejemplo: Todo conjunto con frontera enu-
merable es enumerable.
Si en un espacio topolgico con base contable el conjunto de puntos de
acumulacin de un conjunto A no tiene puntos de acumulacin, A es contable.
Ejercicio. Demuestre, o d un contraejemplo: ABd(Bd(A)) = .
0.10.1.
0.11.
49
50
1
n
<q
k
n=l
l
El teorema de Arqumedes:
x R : x > 0, n N : n > 0 0 <
1
<x
n
Proposicin
Entre 2 reales s y r hay un racional q:
Por teorema de Arqumedes:
1
1
1
rn + 1
< s r(s > r) r < + r < s : q = + r =
n
n
n
n
0.11.1.
0.11.2.
0.11.3.
51
La relacin ser subconjunto de sobre los reales as definidos como subconjunto de racionales provee una relacin de orden lineal para los nmeros
reales.
Ejercicio. Muestre que la relacin ser subconjunto de, es un orden lineal para
R
Por ejemplo, note que {q Q : q 0, , q2 < 2} es subconjunto de {q
Q : q < 3}, correspondientemente decimos que el nmero real raz de dos es
menor o igual que el nmero real tres. Si se tiene un conjunto de nmeros
reales, la unin de los conjuntos racionales correspondientes es tambin un
nmero real, y corresponde al mximo; la interseccin al mnimo.
Ejercicio. Sea A un conjunto de nmeros reales. Muestre que su unin es un
52
reales tiene interseccin no vaca. (De hecho, que hay uno y slo un punto en
tal interseccin.)
0.11.4.
xa . yb mod
z.
Ejercicio. Muestre que para cada real > 0 hay enteros m y n tales que | m
n |<
.
Para terminar esta seccin, podramos decir que con los nmeros naturales contamos, con los enteros llevamos contabilidades y con los racionales
medimos magnitudes fsicas que asumimos que se pueden subdividir en un
nmero arbitrariamente grande de partes iguales. El uso de nmeros reales es
53
antetodo terico ya que ninguna medida tiene precisin definitiva. Sin embargo, su utilidad es grande ya que al modelar magnitudes como nmeros reales
podemos usar tcnicas del clculo para maximizar, obtener transformadas, etc.
Una de las enseanzas del positivismo es que no es conveniente ser positivista
a ultranza: ni Kronecker acept el infinito de los naturales ni Mach acept los
tomos; pero sin el infinito y sin los tomos, la ciencia y la tcnica estaran
innecesariamente limitadas hoy en da.
Los reales por su parte tienen fallas algebricas: hay ecuaciones de la forma
p(x) = 0, con puna funcin polinmica que no tienen soluciones en R. Esta falla
la resuelven los complejos.
Denotemos con bxc el mximo entero menor o igual que x. As, x bxc es
la parte decimal de x.
0.11.5.
El conjunto de Cantor
El conjunto de Cantor es de importancia en topologa y en sistemas dinmicos; tambin en teora de la medida ya que es un conjunto no contable que
mide cero. Lo definimos a continuacin. Considere la representacin en base
3 de los nmeros reales en el intervalo cerrado [0, 1]. Cada nmero all est
representado por (al menos) una sucesin {t0 , t1 , t2 , } donde cada ti {0, 1, 2}.
La idea es que cada nmero lo escribimos como 0.t0 t1 t2 y que el nmero es
P
ti 3i . El por qu todo nmero real tiene tal representacin lo explicamos ms
i=0
adelante. Defina ahora una sucesin de subconjuntos del intervalo [0, 1]. Sea
A0 = [0, 1], sea A1 el conjunto de los nmeros tales que t0 , 1, A2 el conjunto
de los nmeros tales que t0 , 1 y t1 , 1, etc. Ai es el conjunto de los nmeros
T
con t0 , 1, . . . , ti1 , 1. Ahora defina C :=
Ai . Tenemos que:
iRN
54
P 2n1
y en el lmite, C mide 1
3n = 0.
n=1
0.11.6.
0.12.
La circunferencia S1 = T1
En esta seccin definimos la esfera unidimensional S1 , que para dimensin 1 coincide con el toro unidimensional T1 y con la lnea proyectiva RP1 .
Considere las clases de equivalencia de R1 que resultan cuando se considera la relacin congruencia mdulo 2. Dado un nmero real x, su clase de
equivalencia [x] est dada por el conjunto de los nmeros reales de la forma
x + 2n, donde n es un nmero entero. Note que la clase de equivalencia [x]
de cada nmero real x es contable. Vea la figura 16.
Las clases de equivalencia determinan una particin de R1 . Adems, cada
clase de equivalencia [x] tiene un representante nico en el intervalo [0, 2) es
0.12. LA CIRCUNFERENCIA S1 = T1
55
56
Figura 18: El intervalo de S1 de radio con centro en [x], corresponde a una coleccin infinita
de intervalos de la lnea de los nmeros reales
0.12.1.
Circularidad de S1
(mn(JxK, JyK, JzK), med (JxK, JyK, JzK), max(JxK, JyK, JzK)). (Donde med quiere decir mediana) y sus permutaciones pares estn en T.
0.12. LA CIRCUNFERENCIA S1 = T1
57
0.12.2.
58
0.12.3.
El promedio Circular
zi
z = P1N
| 1 zi |
0.12.4.
Promedio y mediana en S1
N
X
e ji
i=1
Claramente la media no se encuetra definida para muestras cuya suma compleja sea cero, tales como las muestras equiespaciadas.
Para una muestra de angular se aplica el proceso de reduccin de forma
reiterativa hasta llegar a una muestra irreducible. Si la muestra resultante es
trivialmente equiespaciada, diremos que es la mediana angular de la muestra.
En caso contrario, diremos que la muestra no tiene mediana.
0.12. LA CIRCUNFERENCIA S1 = T1
0.12.5.
59
La lnea proyectiva P1
[a b]/2 , si a b , /2,
dp ([x], [y]) = mn{ds (a, b) : a x, b y} =
,
de lo contrario
2
Para n = 1, Sn = Pn , topolgicamente. Para dimensiones mayores, el espacio
proyectivo no es homeomorfo a la esfera.
0.12.6.
0.12.7.
Definimos Tn =
n
Q
k=1
0.12.8.
Definimos RPn como el conjunto de las lneas que pasan por el origen
en Rn+1 ; equivalentemente, RPn es la esfera de dimensin n Sn , donde se
han identificado los puntos antpodas, es decir, el espacio topolgico cociente
correspondiente.
60
0.13.
Una vez que se detectaron fallas algebricas en los reales, result la conveniencia de inventar un nuevo tipo de nmero. La frmula para la solucin
cuadrtica era conocida en el ao 2000 A.C. Desde los griegos hasta el siglo
XVI, la solucin de estas ecuaciones fue considerada como la solucin de un
problema geomtrico; as, la ecuacin x2 x 1 = 0 era la interseccin de la
parbola y = x2 y la recta y = x + 1. Cuando la recta y la parbola no se intersecan, se deca que la solucin de la ecuacin era imaginaria (porque no se
poda visualizar). A los reales se les llam as, en oposicin. Ms tarde, el nombre imaginario tom un sentido ms particular y a estos nmeros imaginarios
se les llama ahora nmeros complejos.
Basndose en resultados de del Ferro13 y Tartaglia14 , Cardano15 public
en 1545 su Ars Margna, donde consider races de nmeros negativos por
primera vez en una publicacin, considerando las soluciones respectivas como intiles por carecer de interpretacin geomtrica. El Ars Margna tambin
inclua una solucin de la ecuacin cbica x3 = 3px + 2q correspondiente
aqla interseccin q
de una recta con la curva cbica y = x3 . La solucin es
p
p
3
3
q + q2 + p3 + q q2 + p3 .
Unos 30 aos despus, Bombelli16 not que si la recta es tal que p3 > q2
13 Scipione
61
resulta la raz de un nmero negativo an en casos en los que hay una solucin
geomtrica al problema. Por ejemplo, la recta y = 15x+4, que interseca la cbica
en x = 4, tiene este problema. Bombelli tuvo la (gran) idea de suponer que para
3
3
la ecuacin x3 = 15x + 4, 2 + 11i es de la forma 2 + ni y 2 11i es de la forma
2 ni. Definiendo la multiplicacin y la adicin de nmeros complejos (como
la conocemos hoy en da) encontr que las cosas an cuadraban bien en la
frmula de Cardano.
62
63
porque implica contradicciones, definirle coordenada polar (de fase o ngulo) al nmero complejo , con coordenadas rectangulares (0, 0).
64
0.13.1.
Nomenclatura
En coordenadas rectangulares, la primera componente de un nmero complejo se denomina su parte real y la segunda su parte imaginaria. Las funciones
Re e Im estn definidas as:
Re(x, y) = x (parte real)
Im(x, y) = y (parte imaginaria)
Resulta conveniente denotar los nmeros con coordenadas rectangulares
65
(0, 0), (1, 0) y (0, 1) abreviadamente, as: := (0, 0), 1 := (1, 0) y j := (0, 1).
Por lo tanto, en coordenadas rectangulares, (a, b) = a(1, 0) + b(0, 1) = a1 + bj,
que escribimos simplemente como a + jb.
En coordenadas polares, la primera componente de un nmero complejo
se conoce
como su magnitud mientras que la segunda se conoce como su fase:
(r, ) = r (magnitud)
(r, ) = (fase).
Como dijimos, al origen del plano no le definimos fase.
0.13.2.
Equivalencias de ngulos
Ejercicio. e j2
0.13.3.
1/2
= e j ?
Estructura algebrica
Lo que diferencia al conjunto de los nmeros complejos del plano de coordenadas rectangulares y del cilindro de coordenadas polares, es la estructura
algebrica dada por las operaciones (producto) y + (suma) se definen
como sigue:
+ (suma): (x, y) + (a, b) = (x + a, y + b) y es una funcin asociativa, conmutativa, con mdulo aditivo y con inversos aditivos.
66
vas.
Sea z = (x, y), si (x, y) + (a, b) = , entonces (a, b) = (x, y), por lo tanto el
inverso aditivo de z est dado por z = (x, y).
Cada nmero complejo z, con excepcin de , tiene un inverso multiplicativo
que denotamos z1 . Para obtener una frmula en coordenadas cartesianas para
z1 , procedemos como sigue.
Si (x, y) (a, b) = (1, 0) entonces, resolviendo para a y b encontramos los siguientes valores, siempre y cuando (x, y) , ,
(x, y)
y
x
=: (a, b) = 2
, 2
2
x + y x + y2
0.13.4.
67
Figura 24: Equivalencia geomtrica de las representaciones polar y rectangular de los nmeros
complejos
p
x2 + y2
arctan(y/x),
2,
=
2 ,
arctan(y/x) + ,
si x > 0,
si x = 0, y > 0
si x = 0, y < 0
si x < 0
68
de las magnitudes)
Ejercicio. Muestre que (zw) = z + w. (La fase de un producto es la suma de
las fases)
Ejercicio. Muestre que |z + w| |z| |w|. (La magnitud es subaditiva)
Ejercicio. Muestre que ||z| |w|| |z w| (Vea la figura 27)
Ejercicio. Sea g : R1 C dada por g(t) = cos t. Grafique la magnitud y la fase
69
Figura 26: Una forma de definir la funcin arctan: con rango (/2, /2)
e j
1+e j
sen()
0.13.5.
0.13.6.
70
71
Figura 29: Representacin geomtrica del inverso multiplicativo, el inverso aditivo y el conjugado
nadas polares?
0.13.7.
Complejo conjugado
Ejercicio. Muestre que |z|2 = z2 .
Ejercicio. Muestre que |z | = |z|, z = z.
Ejercicio. Muestre que, si z , 0, z1 =
z
.
|z|2
72
tan()+tan()
1tan() tan() , muestre que zw
= zw.
Re[z] Im[z]
Im[z] Re[z]
muestre que las operaciones de suma y producto de complejos coinciden con
las correspondientes de matrices.
0.13.8.
La proyeccin estereogrfica
73
74
=
z/ = , si z ,
z/ = , si z ,
0.13.9.
Funciones lineales
0.13.10.
sobreyectiva. tambin, dibuje las imgenes con respecto a f de los ejes real e
imaginario, de las lineas Re[z] = Im[z], Re[z] = Im[z], de las lineas Re[z] = 1,
Re[z] = 1, Im[z] = 1, Im[z] = 1 y de las curvas |z| = 1, |z| = 0,5, |z| = 2.
Ejercicio. Diga si la funcin f : C C, dada por f (z) = ez es inyectiva o
sobreyectiva. tambin, dibuje las imgenes con respecto a f de los ejes real e
75
imaginario, de las lineas Re[z] = Im[z], Re[z] = Im[z], de las lineas Re[z] = 1,
Re[z] = 1, Im[z] = 1, Im[z] = 1 y de las curvas |z| = 1, |z| = 0,5, |z| = 2.
Ejercicio. Diga si la funcin f : C C, dada por f (z) = ez tiene inversa
0.13.11.
Transformaciones de Mbius
76
invariantes, y que las partes positiva y negativa del eje imaginario se intercambian. Finalmente, veamos que, con el procedimiento indicado, la magnitud de
cada punto resulta invertida.
Si mostramos que los ngulos y son iguales, el resultado buscado se
sigue, ya que entonces los tringulos 0NB y CNA son similares por tener 2
1
.
ngulos iguales, entonces N0 = NC , entonces 1 = NC
1 osea 0B =
0B
CA
0B
0A
1z
1+z
en la forma ya indicada.
la forma indicada.
Ejercicio. Muestre que la suma de dos funciones de Mbius es una funcin de
Mbius.
Ejercicio. Muestre que la composicin de dos funciones de Mbius es una
funcin de Mbius.
0.14. RN Y CN
77
1z
1+z .
0.14. RN y CN
El producto cartesiano RN , N N, tiene una estructura natural de espacio
vectorial, dada por la base {ei : i /1, n/} con ei = (0, ..,0, 1, 0, ..,0) (con un
solo 1, en la posicin i de de la n-tupla). A esta base se le conoce como la
base (ortonormal) estndar de RN , y con cuerpo de escalares dado por los
reales. Toda transformacin lineal tiene una matriz asociada y, en el caso de
una transformacin lineal RN RM resulta que los vectores columna son las
imagenes de los vectores de la base estndar de RN .
La misma base, con cuerpo de escalares C, sirve para el espacio CN .
Sinembargo, con cuerpo de escalares R, la base tiene 2N elementos: {ei : i
/1, n/} {di : i /1, n/} con ei = (0, ..,0, 1, 0, ..,0) y di = (0, ..,0, j, 0, ..,0) (con una
sola j, en la posicin i de de la n-tupla). En un caso la dimensin del espacio es
N y, en el otro, 2N.
78
0.15.
0.15.1.
Definicin
Si x es un punto en R4 = R1 R3 (como explicamos abajo, resulta conveniente pensar que las ltimas 3 componentes de un cuaternio son las coorde-
79
0.15.2.
Suma y producto
La operacin suma se define en forma estndar, con la ley del paralelogramo, como (a1+bi+cj+dk)+(1+i+j+k) = (a+)1+(b+)i+(c+)j+(d+)k,
mientras que el producto queda determinado por la ley distributiva y los productos elementales:
11 = 1
ii = jj = kk = 1
ij = ji = k
ki = ik = j
jk = kj = i
Note que el producto resultante no es conmutativo. El resultado del producto se puede recordar con la ayuda de la figura 34; el producto de i, j, o k
con s misma, es -1. El producto de elemento i, j, o k con otro diferente, es la
que no se usa, con signo positivo cuando seguimos el orden indicado en la
80
0.15.3.
0.15.4.
81
1
1+i+ j+k .
Ejercicio. Diga por qu, el producto en R3 dado por [x, y, z][a, b, c] = [yc
zb, xc az, xb ay] junto con la operacin normal suma, no permite tener un
anillo de divisin sin divisores de cero.
Un cuaternio es unitario si su magnitud vale 1. Por lo tanto, el anillo de los
cuaternios unitarios es la esfera tridimensional S3 (el conjunto de los puntos
de R4 a distancia unitaria del origen). En forma anloga a como los nmeros
complejos de magnitud 1 viven en la esfera unidimensional, los cuaternios de
magnitud 1 viven en la esfera tridimensional.
Si u es un cuaternio puro, u 2 = [0, u1 , u2 , u3 ][0, u1 , u2 , u3 ] = [u1 u22
2
u3 , 0, 0, 0] osea, su cuadrado es un escalar (de valor la norma). En general, el
e = u e u e. El cuadrado
producto de dos cuaternios puros est dado por u
de un cuaternio unitario puro es -1. Esto nos permite escribir, usando las series
82
0.15.5.
83
nante vale 1.
Alternativamente, considere el caso especial en que q = cos + i sen . con
e = e1 i + e2 j + e3 k,
Lq (e) = q e1 i + e2 j + e3 k q
= qq ie1 + q jq e2 + qkq e3
2
= ie1 + q je2 + q e3
= ie1 + (cos 2 i sen 2)(e2 + e3 i) j
y tenemos que el eje x (el subespacio generado por i) permanece invariante
mientras que el plano yz (el subespacio generado por {j, k}) se rota en un
ngulo de valor 2. (Lo cual es interesante ya que entonces [0, ], como en
el caso de la recta proyectiva).
en el caso general, q = cos + e sen , descomponemos a e en una compo y otra paralela. Usando e = ek +e , donde con = e u,
ek = |u|u
nente normal a u,
y e = e |u|u (la componente perpendicular a u),
84
queda:
qeq = (cos + u sen )(ek + e )(cos u sen )
= ek cos2 + e cos2 ek u cos sen e u cos sen
k sen cos + ue
sen cos ue
k u sen2 ue
u sen2
+ ue
= ek cos2 +ek sen2 + 0 + (e cos2 e sen2 + 2e u sen cos )
= ek + (e cos2 + 2e u sen cos )
donde vemos que el sumando paralelo permanece invariante mientras que el
sumando ortogonal se rota un ngulo 2.
0.15.6.
0.16. SUCESIONES
85
a j
aa +
1
aa bb)
((a + b)(a + b)
2
1
(ab ba)
2
La norma satisface |xy| = |x||y|. La base 1, i, j, i j es ortonormal. El ordenamiento da una orientacin para H.
=
0.16.
Sucesiones
86
0.16. SUCESIONES
87
0.16.1.
Subsucesiones
Ak
n=0 k=n
\
n=0 k=n
Ak
88
0.17.
Algunos espacios topolgicos son espacios mtricos y algunos de los espacios mtricos son espacios normados. Una norma determina una mtrica y
una mtrica determina una topologa.
Un espacio mtrico (X, d) consiste en un conjunto X y una funcin d :
X X [0, ), que se conoce como mtrica, tal que:
i. x0 inX, d(x, x) = 0
ii. x, y X, si d(x, y) = 0, entonces x = y (definitivamente positiva)
iii. x, y X, d(x, y) = d(y, x), (simetra)
iv. x, y, z X, d(x, y) + d(y, z) d(x, z) (desigualdad del tringulo)
En un espacio mtrico, la bola con centro c X y radio r (0, ) est dada
por el conjunto {x X : d(x, c) < r}, que denotamos Bc (r).
Ejercicio. Muestre que la funcin d : R2 R2 [0, ), dada por d(x, y) =
p
(x + y)2 + (x y)2 para cada x, y R2 , con x = (x1 , x2 ), y = (y1 , y2 ), determina
una mtrica para R2 . Esta mtrica se conoce como la mtrica eucldea para R2 .
Ejercicio. Para la mtrica del ejemplo anterior, dibuje la bola con centro en el
mtrica dada por d(x, y) = |x1 y1 | + |x2 y2 |. Esta mtrica se conoce como la
mtrica del taxista.
Ejercicio. En R2 , para cada x, y R2 , con x = (x1 , x2 ), y = (y1 , y2 ), considere la
mtrica dada por d(x, y) = max{|x1 y1 |, |x2 y2 |}. Esta mtrica se conoce como
la mtrica del producto.
0.17.1.
89
Bolas
0.17.2.
Mtrica para S1
,
de lo contrario
as, estamos tomando la longitud mnima de los recorridos entre x y y.
90
0.17.3.
Espacios normados
0.17.4.
Espacios topolgicos
puntos de acumulacin.
91
Equivalentemente, dado un subconjunto A de un espacio mtrico X, decimos que el punto x X es un punto de acumulacin de A, si cada bola que
contenga a x, contiene un punto de A diferente de x.
A partir de una norma se puede definir una mtrica d(x, y) = |x y| pero no
viceversa: se necesita un espacio vectorial y la propiedad de la homogeneidad
positiva. Similarmente, con una mtrica se puede definir una topologa, pero
no viceversa: la desigualdad del tringulo no es una propiedad topolgica.
0.17.5.
Cauchy supo de la idea de Dedekind para definir los nmeros reales pero
de alguna manera, no le satisfaca. Cauchy dio otra definicin de los reales,
basada en sucesiones de racionales en vez de subconjuntos de racionales. La
desventaja es que un nmero real, con la definicin de Cauchy, es una clase de
equivalencia de sucesiones (algo ms complejo que un corte de Dedekind). La
ventaja es que la idea se generaliza y se aplica a espacios mtricos en general.
0.17.6.
92
para cada positiva hay un nmero natural N tal que, siempre que n y m sean
mayores que N, d(sn , sm ) < Esto lo denotamos como d(sn , sm ) 0. Si una
sucesin es convergente, claramente es una sucesin de Cauchy: esto se sigue
de la desigualdad del tringulo: d(sn , sm ) d(sn , l) + d(l, sm ).
Ejercicio. Diga si R2 es un espacio mtrico completo. Diga si (0,1) lo es. Diga si
0.17.7.
93
a. an = (1)n
b. an = (0,9)n
c. an = ( j)n
d. an = ( n2 )
n
e. an = ( n+1
)
1
)
f. an = sen( n+1
g. an = (n)n
h. an = sen n1
PRUEBA M de Weierstrass.
Sea { fn } una sucesin de funciones definidas en A, y sea {mn } una sucesin de
P
nmeros tales que x A, | fn (x)| mn . Suponga que
mn converge. Entonces
x A
a
n=1
n=1
fn (x).
n=1
94
Lema 3. Cualquier sucesin {an } contiene una subsucesin que es o bien no decreciente
o bien no creciente.
Demostracin. Definimos un punto cumbre a el nmero n tal que am < an m <
n. Un punto valle se define de manera anloga.
CASO 1. La sucesin tiene infinitos puntos cumbre. Entonces estos conforman
una sucesin no creciente.
CASO 2. La sucesin tiene un nmero finito de puntos cumbre. Encuentre el
ltimo de los puntos cumbre, el siguiente punto no es punto cumbre, tome ste.
Como no es cumbre, hay puntos a nivel o ms abajo, tome uno de estos puntos
y repita el procedimiento. Lo anterior resulta en una sucesin no creciente.
Para la sucesin no decreciente usamos los puntos valle.
Corolario 2. Toda sucesin acotada tiene una subsucesin convergente.
Teorema 3. Una sucesin {an } converge si y slo si es un sucesin de Cauchy.
Demostracin. n, m > N |an am | < n > N |an aN | <
0.17.8.
95
ahora proseguimos con la idea de Cauchy para definir los nmeros reales.
suponga que inicialmente tenemos solamente los nmeros racionales (y los
enteros y los naturales). Note que Q es un espacio mtrico con la mtrica
racional : Q Q dada por d(p, q) = |p q|. Se puede tener una sucesin
de Cauchy {sn } de racionales, (es decir (sn , sm ) 0) pero que no converja a
ningn racional, indicando que los racionales son completos.
Ejercicio. D un ejemplo de una sucesin de Cauchy de racionales que no
96
n
xn+1 xn + 2yn
x0 7
ecuacin de recurrencia,
=
con, = los primeros valores
yn+1
xn + yn
y0
5
son 7/5 = 1,4, 17/12 = 1,416 . . . , 41/29 = 1,413 . . . , . . . .
ejemplo anterior.
Ejercicio. Muestre que la sucesin {qn } del ejemplo anterior es de Cauchy.
lentes.
Ejercicio. Diga si Q es completo.
Ejercicio. Diga si el intervalo de reales (0, 1) es completo.
0.18. SERIES
0.18.
97
Series
Definicin 20. El trmino serie (series en ingls) quiere decir aqu el lmite de
una sucesin de sumas parciales. En otros contextos, por ejemplo una serie de tiempos
(time series) es una sucesin de datos tomados a medida que avanza el tiempo. En
cierto sentido, lo que un profesional de la estadstica llama una serie de tiempos es lo
que un ingeniero llama una seal discreta.
si la sucesin {an } converge y, el tal caso, lmite. Tal lmite lo denotamos as:
n1
P
1 + (1/2)i y lo llamamos el valor de la serie. Si ha ledo sobre la paradoja de
i=1
Zenn, y sabe el resultado de la serie, quizs vea aqu una solucin matemtica
a la paradoja, es decir, tenemos que el resultado de una suma infinita es finito:
n1
P
1 + (1/2)i = 2.
i=1
n
P
i=0
1
i+1
P
n=1
1
n
dola con las reas indicadas en la figura 38, sabiendo que la integral de
[0, ) es infinita.
1
x
sobre
98
0.18.1.
P
i=0
i=0
Es posible que
i=0
i=0
P
(1)n+1
converge a ln 2, lo que podramos
Por ejemplo, resulta que la serie
n
n=1
+
2 4 3 6 8 5 10 12
1
1
1 1
1
1
1
1
= (1 ) + ( ) + ( )
+
2
4
3 6
8
5 10
12
1 1 1 1
1
1
= + +
+
2 4 6 8 10 12
1
1 1 1 1 1
= (1 + + + )
2
2 3 4 5 6
1
= ln 2
2
P
P
Decimos que la serie
zi converge en magnitud, si la serie | zi | converge.
i=0
i=0
0.18. SERIES
99
incondicionalmente.
0.18.2.
R
P
P
1
1
1
que
>
n
x dx = , (vea la figura 39) concluimos que la serie
n diverge.
n=1
n=1
P
n=1
N
P
n=1
1
n2
<1+
RN
1
dx
x2
1
n2
converge ya que
< 1 + 1 = 2.
P
n=1
1
,
n2
100
0.18. SERIES
0.18.3.
101
La serie geomtrica
zn ,
n=1
n=0
N
P
zn
n=0
tambin es til y tambin tiene una frmula sencilla, que deducimos a continm
1zm+1
m
uacin, y que, 1z
1z + z = 1z . Esto permite demostrar la frmula de la suma
N
P
de potencias de z. Para evaluar
zn procedemos inductivamente. el truco a
n=0
1z2
1z ,
N
P
n=0
zn =
1zN +1
1z ,
si z , 1, y, claro est,
zn = N + 1, si z = 1.
N
P
n=0
n=0
P
1
; si |z| 1, zN = |z|N e jN diverge. Por lo tanto, la frmula para
cero y
zn es 1z
n=0
P
n=0
Ejercicio. Calcule
P
n=0
2n .
zn =
1
1z
102
Ejercicio. Calcule
2n .
n=1
0.18.4.
La transformada zeta
P
zeta como la funcin W : D C dada por la serie de Laurent W(z) =
wn zn ,
P
n=0
zn =
1
1z ,
con regin de
0.18. SERIES
103
1/2
1
=
2 z 1 z/2
1X z n
=
( )
2 n=0 2
1 X z n
( )
2 n=0 2
1 X n n
2 (z)
2 n=0
Ejemplo.
1
(1 z)(2 z)
1
1
=
1z 2z
1
1
=
1
z(1 z ) 2 z
1
1
=
+
1
z(1 z ) 2z(1/2 z1 )
W(z) =
Corresponde a diferentes sucesiones en cada una de las regiones |z| < 1, 1 <
|z| < 2 y 2 < |z|.
104
0.19.
0.20. INTERVALOS DE RN Y DE ZN
105
0.20.
Intervalos de Rn y de Zn
en la figura (44
Sean r = [r1 , r2 , . . . , rN ] y s = [sa , s2 , . . . , sN ] elementos de ZN . El intervalo
entero, de dimensin N, /r, s/ es el conjunto de elementos k = [k1 , k2 , . . . , kN ]
tales que, para cada i en /1, N/, r j ki si ,
Las uniones de intervalos son los dominios naturales de las seales. Los intervalos multidimensionales son los dominios de las seales multidimensionales; por ejemplo, las imgenes tienen como dominio intervalos de dimensin
dos y se consideran seales de dimensin dos.
Por otra parte, los intervalos proveen bases para la topologa estndar de
n
R . Estos intervalos se conocen como las bolas con respecto a la mtrica del
producto.
106
0.21.
Clculo
0.21.1.
por fn (x) = xn , converge a la funcin c(x) dada por c(0) = 0 y c(x) = 1 si x (0, 1].
0.21. CLCULO
107
Observe la figura siguiente. note que las fn son continuas y que el lmite c no
lo es.
0.21.2.
Convergencia uniforme
Wierstrass.
108
0.21.3.
Figura 46: En la definicin de lmite, los puntos del intervalo J, con la posible excepcin de a,
tienen imgenes dentro del intervalo I; en el caso de la definicin de continuidad, la imagen de J
es un subconjunto de I
Se puede comprobar que, cuando el rango y el dominio son espacios mtricos,
con mtricas d y D, la definicin anterior es equivalente a: > 0, > 0
0 < |x a| < d( f (x), b) < .
Cuando el rango y el dominio son R o C, tenemos la definicin usual:
> 0, > 0 0 < |x a| < | f (x), b| < .
0.21. CLCULO
0.21.4.
109
Continuidad
0.21.5.
0
Continuidad uniforme
intervalo (0,1).
Ejercicio. Muestre que la funcin sen t es uniformemente continua.
110
0.21.6.
f : R1 S1 dada por f () = e j .
f : R1 S1 dada por f () = (e j 1).
0.21.7.
Derivada
f (x + h) f (x)
h
0.21. CLCULO
111
112
Figura 47: La ecuacin de la tangente est basada en la funcin lineal (x) = f 0 (x)(t x)
0.21.8.
0.21. CLCULO
113
Figura 48: El ngulo con vrtice z0 y el ngulo con vrtice f (z0 ) son iguales
usando f 0 (1 + j).
Si f tiene derivada compleja en cada uno de los puntos de una regin
abierta, se dice que f es holomorfa en la regin.
0.21.9.
Aproximacin lineal
114
J(x) =
fM (x)
x1
0.21.10.
f1 (x)
xN
fM (x)
xN
Diferencial
h0
para RN y RM .
Se puede mostrar que T es nica y se conoce como el diferencial T f (c), y
tambin como la funcin tangente a f en c. T f (c) es la matriz Jacobiana
M N de f .
Ejercicio. Sea f : R2 R1 dada por f (x, y) = x2 + y2 . Encuentre la ecuacin del
0.21.11.
L =
y
b a y
0.21. CLCULO
115
x u(x, y)
f =
y
v(x, y)
y si es diferenciable su Jacobiano est dado por:
x u
J =
x v
y u
y v
LL = J
y
y
y si f tiene derivada compleja w,
x u
v
x
y u a b
=
y v b a
f () =
x4 y4
x3 +y3
116
0.21.12.
n=1
z0 )n .
an xn .
n=0
n=1
2. g(z) = b0 +
bn (z
n=1
bn (z z0 )n + h(z) = b0 +
como
P
n=
n=1
bn (z z0 )n .
P
nan xn1 (Por qu?) y, por lo tanto f 00 (0) = (2)(1)a2 , en general, procedin=1
los coeficientes: an = n! .
Por ejemplo, las funciones seno y coseno son funciones analticas mientras
que la funcin |x| no lo es. Las series de Taylor son usadas en los programas
de computadores y calculadoras que usan aproximaciones a funciones tales
como la raz cuadrada y las trigonomtricas.
Una diferencia importante entre funciones de la variable real y funciones
de la variable compleja, es que en el caso complejo, tener derivada compleja es
0.21. CLCULO
117
condicin suficiente para que una funcin de variable compleja tenga derivada
(continua) de cualquier orden.
Ejercicio. Que dice el teorema de Morera?
de cero?
Ejercicio. Cul es la de la funcin coseno?
Ejercicio. Si una serie converge en valor absoluto o en magnitud, se puede
P
n=0
(x)n
n!
P
n=0
(y)n
n!
P
n=0
(x+y)n
n!
P
n=0
1
n!
P
n=0
xn
n!
y dado que f (x + y) =
P
n=0
reales x y y.
Sugerencias:
(x+y)n
n!
P
n=0
(x)n
n!
P
n=0
(y)n
n!
118
P
C (con dominio C) utilizando la misma serie de Taylor: g(z) =
cn zn por
ejemplo, si cn =
n=0
1
n! ,
e j
1+e j
1
2
sen
+ j 2(1+cos
)
0.21. CLCULO
0.21.13.
119
Series de Laurent
P
i=
si =
si +
i=0
si es la suma
i=
i=1
n=
cn (z z0 )n .
El residuo de una funcin con expansin en series de Laurent es el coeficiente del trmino que incluye a z1 .
Teorema 4. si se sabe que la serie
P
n=
concluir que la serie converge para cada z con |z| < |z0 |.
Ejercicio. Demuestre el teorema anterior.
0.21.14.
Integral de Riemann-Stieltjes
0.21.15.
Particiones de intervalos
120
0.21.16.
0.21. CLCULO
121
0.21.17.
La integral de Riemann
f (x)dx :=
a p
N1
X
i=0
f (x)dx :=
a p
N1
X
i=0
f (x)dx = lm
a p
f (x)dx = lm
a
a p
122
si se cumple la igualdad
b
Z
f (x)dx =
f (x)dx
a
decimos que f tiene integral de Riemann y que esta est dada por este valor
comn
Z b
Z b
Z b
f (x)dx =
f (x)dx =
f (x)dx
a
1
Ejercicio. Sea f (x) =
0.21.18.
xQ
xRQ
La integral de Riemann-Stieltjes
0.21. CLCULO
123
Nomenclatura
Una integral de Riemann-Stieltjes se expresa as:
b
Z
f d o
f (x)d(x)
a
R2
0
cos x d(x).
124
Z
f (x)d(x) +
Calcule
R3
0
f (x)d(x) +
R3
0
(x)d f (x).
Rb
a
f (x)d(x) =
0.21. CLCULO
125
t1
Z
f (x)d(x) =
f (x)0 (x)dx +
XZ
ti
f (x)0 (x)dx
ti1
f (x)0 (x)dx +
tn
(b)
Z
f (x)d(x) =
(a)
f (1 (x))dx
f (x)dx =
a
dG(t)
= f (t)
dt
de
1 t+h
h t
f (x)dx es f (t).
126
f (z)d
Z
Re[ f ()(t))]dRe[)(t)]
Im[ f ()(t))]dIm[)(t)]
a
Z
+j
Z
Re[ f ()(t))]dIm[)(t)]
!
Im[ f ()(t))]dRe[)(t)]
0.21. CLCULO
127
f (z)
wz dw
= f (z)
1
wz dw,
donde
Rx
a
fn (t)dt =
Rx
a
lm fn (t)dt.
128
n=0
tonces
lm
m Z
X
n=0
Z
fn (t)dt =
lm
m
X
fn (t)dt
n=0
1
ba
f (t)dt
a
Note que f es la altura del rectngulo con base [a, b] y la misma rea que la
0.21. CLCULO
0.21.19.
129
Integral de lnea
f d donde
f d =
f ((t))d(t)
a
si es derivable, tenemos
b
f ((t))0 (t)dt
a
Z
zn d =
S1
e jnt ( je jt )dt
0
2
Z
=j
e j(n+1)t dt
0
2 j si n = 1
=
0
si n Z {1}
Ejercicio. Qu dice el teorema del residuo?
0.21.20.
P
{xn }, entonces X(z) =
xn zn .
n=
130
Z
zn d =
= jrn+1
e j(n+1)t dt
0
jr 2 n = 1
=
0
n , 1
2 j n = 1
=
0
n , 1
(El resultado no cambia.)
Ejercicio. Muestre que el resultado no cambia tampoco si es una curva simple
I
X(z)d =
Z
=
xn zn d
n=
2 X
0
xn zn d
n=
Z
X
n=
xn zn d
0
= 2 jx1
Ejercicio. D la frmula para xn en general.
0.21. CLCULO
0.21.21.
131
El funcional Dirac
El funcional : CR C dado por (g) = g(0), tiene representacin de Riesz
R
de la forma (g) = g(x)du(x), donde u es la funcin escaln.
0.21.22.
Transformadas
132
g(x)H(x)dx =
2. g(t) =
R
1
2
g(t) f (t)dt.
f (t)dg(t).
0.21.23.
La distribucin de Dirac
de D es la funcin constante 1.
0.21. CLCULO
0.21.24.
133
La distribucin constante
R
Sea dada por ( f ) = f (t)dt, que, como vemos, tiene ncleo la funcin
constante 1. Por definicin, su transformada est dada por,
K( f ) = (F)
Z
=
F()d
Z Z
=
f (t)e jt dtd
Z
Z
?
=
f (t)
e jt ddt
0.21.25.
u(t) f (t)dt =
R
0
f (t)dt.
M( f ) = (F)
Z
=
F()d
Z0 Z
=
f (t)e jt dtd
0
Z
Z
?
=
f (t)
e jt ddt
R
0
134
0.21.26.
R
||>0
f (t)
t dt. Por
( f ) = (F)
Z
F()
d
= lm
0 ||>0
Z
Z
1
= lm
f (t)e jt dtd
0 ||>0
Z
Z
ejt
?
ddt
=
f (t) lm
0 ||>0
R
jt
resultando un distribucin con ncleo ||>0 e d = jsgnt, donde sgn es la
funcin sgnum. (sgn vale 1 para argumento positivo, -1 para argumento negR jt
ativo y 0 para argumento nulo). Sugerencia: evale e dt usando clculo de
0
residuos.
As, en cierto sentido, 1/x y jsgn(x) son un par de Fourier.
Variacin de seales discretas
Sea X : Z C. La variacin total de la seal se define como
|xi
i=
xi1 |.
Variacin promedio
N
P
1
|xi
2N+1
N
i=N
xi1 |.
1
N
P
i=
|xi xi1 |
0.21. CLCULO
135
Problemas
1. Muestre o d un contraejemplo en relacin con la afirmacin: Si una
relacin es simtrica y transitiva, entonces es reflexiva.
2. C y D son dos subconjuntos del rango de f . Demuestre o d un contraejemplo:
a) f 1 (C D) = f 1 (C) f 1 (D)
b) f 1 (C D) = f 1 (C) f 1 (D)
3. C y D son dos subconjuntos del dominio de f . Demuestre o d un contraejemplo:
a) f (C D) = f (C) f (D)
b) f (C D) = f (C) f (D)
4. Muestre o d un contraejemplo: Para ningn conjunto hay una biyeccin de
ste a su conjunto potencia.
5. D ejemplos de conjunto finito, contable infinito y no contable.
6. Muestre que hay infinitas ternas pitagricas. (Nmeros naturales a, b, c
tales que a2 + b2 = c2 )
136
9.
0.21. CLCULO
137
j
1+ j .
1
1+ j .
2024
P
e j2/7 .
n=0
138
29. Sea f : Z C dada por f (n) = e j(2/5)n , indique en el plano complejo los
puntos f (n) para n / 10, 10/.
30. Grafique la magnitud y la fase de f contra n, para n /0, 10/, donde f es
la funcin definida en el problema anterior.
31. Encuentre k //0, 4 tal que e j(2/5)1024 = e j(2/5)k .
32. Sea z = e j2(1/)(12/15) . D el mnimo valor entero de n, positivo, tal que
zn = 1.
33. Encuentre un polinomio equivalente a a(1 z9 )/(1 z), para z , 1.
34. Sea f : C C dada por f (z) = (1 + 2j)z. Sea T el tringulo mostrado.
dibuje la imagen f (T) de T.
0.21. CLCULO
139
8
P
n=0
(1 + j)n .
140
a)
b)
47. Encuentre la extensin en series de Taylor de ln (1 + x) alrededor de x = 0.
48. Encuentre la extensin en series de Taylor de
1
1+x
alrededor de x = 0.
P
n=1
1
n
0.21. CLCULO
b)
P
n=1
51. Exprese
141
1
n2
9
P
n=0
P
n=1
5
P
n=1
j
( 2 )n
2
P
n=1
( 2 )n y
en el plano complejo.
54. Sea g : (1, 1] R dada por g(x) = ln (1 + x). Encuentre los coeficientes
P
de la serie de Taylor a0 +
an xn de g al rededor de cero. Diga si la serie
n=1
P
por la izquierda en 1). Diga si f 0 (1 ) =
nan xn1 con x = 1.
n=1
55. Encuentre la transformada zeta U(z) con su respectiva regin de convergencia (RC) del escaln discreto un (un = 1, si n 0, un = 0, si n < 0). A
continuacin, encuentre otra seal discreta con transformada zeta U(z)
pero con diferente regin de convergencia.
56. Demuestre que si
n=0
10
P
n=10
P
n=0
1
zn = z10 1z
1z .
an tambin.
142
sen t
t ,
t R1 .
R
1 h
h 0
cos xdx.
P
71. Sean f, : R R dadas por (x) =
u(x n 0,5) y, g(x) = x2 , para
n=0
R
x , 0 y g(0) = 0. Calcule g(x)d(x).
R
R
72. Sea (t) = u(t 0,25) u(t 0,75). Calcule 0 sen(t)d(t), 0 sen2 (t)d(t)
R
y 0 cos(2t)d(t).
R 2
73. Calcule la integral de lnea ez d donde es el camino dado por
(t) = 1 + jt, t[0, 1], (t) = j + 1 t, t[1, 2].
R 2
74. Calcule et ejt dt.
75. Calcule
R1
1
0.21. CLCULO
143
R 2
76. Calcule 0 f ((t))d(t), donde f : C C est dada por f (z) = z2 y
: R C est dada por (t) = e j(t+u(t1)) (u(t) es el escaln).
77. Sea f : [0, ) R1 dada por f (x) = 0 si x = 0 y f (x) = sen(1/x) si x > 0.
Diga si f es de variacin acotada en [0, 1]. Explique.
78.
79. Demuestre que la funcin f : C {0} C dada por f (z) = 1z puede ser
implementada en la siguiente forma:
Coloque el polo sur de una esfera de radio 12 sobre el origen del plano.
Suba los puntos del plano a la esfera, usando proyeccin estereogrfica.
Rote la esfera manteniendola tangente al plano, y manteniendo invariantes las imgenes del eje real, de tal forma que al final, con proyeccin
estereogrfica inversa, la imagen del punto que parti del punto j del
plano inicialmente, aterrice ahora sobre el punto j.
80. Demuestre o d un contraejemplo: para cada nmero complejo z , 0,
para todos nmeros complejos w1 , w2 y w3 , no colineales en el plano, los
tringulos w1 w2 w3 y zw1 zw2 zw3 son similares.
81. Represente los productos de complejos y de cuaternios con matrices; es
decir, diga cmo obtener una matriz a partir de cada complejo (o cuaternio) de tal forma que al multiplicar las matrices correspondientes, resulte
la matriz correspondiente al producto de los complejos (o cuaternios).
144
1
n
1
n=1 n2 .
1
x4 +y4
x4 y4
x3 +y3
x4 +y4
+ j x3 +y3 derivable.
1+z
,
z3 +2z2 +5 6
0.21. CLCULO
145
z1
z+1 .
95. Diga si la siguiente funcin tiene derivada compleja en el origen. Tambin, evale las ecuaciones de Cauchy-Riemann all.
f (x + jy) =
x4 y4
x4 + y4
+
j
x3 + y3
x3 + y3
146
1n
|n|+1
0.21. CLCULO
147
2 3
5 6
8 7
f (x, y, z) = det 4
x
y
z
3
6
7
148
5
f (w, x, y, z) = det
9
2
x
8
2
3
7
y
3
4
8
0
z
Referencias
[1] Mathematical Analysis. Tom M. Apostol. Adyson-Wesley,
Reading MA, 1974. (Integral de Riemann-Stieltjes)
[2] Essay on Theory of Numbers. Richard Dedekind. Dover, N.Y.,
1963.
[3] Science and Hypothesis. Henri Poincar. Dover, N.Y., 1963.
[4] Contributions to the Founding of the Theory of Transfinite
Numbers. Georg Cantor. Dover, N.Y., 1955.
[5] Geometry of Complex Numbers. Hans Schwerdteger. Dover,
N.Y.,1979.
[6] Elementary Theory of Analytic Functions on One or Several
Complex Variables. Henri Cartan. Dover, N.Y., 1995.
[7] Principles of Mathematics. Bertrand Russell. Norton, N.Y.,
1902. (Relaciones)
[8] The Large, the small and the Human Mind. Roger Penrose.
[9] Theory of Complex Functions. Reinhold Remmert. Springer,
N.Y.,1991.
149
150
REFERENCIAS
[10] The Real Line and Calculus. William Eaton. Class notes for
665A. The University of Texas at Austin, 1985. (Huecos de Q)
[11] Visual Complex Analysis. Tristan Needham. Clarendon, Oxford, 1997. (Visualizaci del producto de complejos)
[12] Quaternion-Fourier Transforms for Analysis of Two Dimensional Linear Time-Invariant Partial Diferential Systems. Proc.
32nd Conf. Decision Contr., Dec. 1993, pp. 1830-1841.
[13] Calculus. Clculo Infinitesimal. Michael Spivak. Revert,
Mxico, 1992. (Sucesiones, definicin de reales)
[14] Georg Cantor, His Mathematics and Philosophy of the Infinite. J.W. Dauben. Princeton University Press, Princeton, 1979.
[15] Hypercomplex Fourier Tansform of Color Images. IEEE International Conference on Image Procesing (ICIP 2001). S.J. Sangwine. Thessaloniki, Grece, Oct. 2001, vol. I, pp. 137-140.
[16] Quaternions and Rotation Sequences. J.B. Kuipers. Princeton
University Press, Princeton, 1999.
[17] Computational Conformal Mapping. P.K. Kythe. Birkhuser,
Boston, 1998. (ecs- Cauchy-Riemann y diferenciabilidad)
[18] On Quaternions and Octonions. J.H. Conway and D.A. Smith.
AK Peters, Natick, 2003.
[19] Elementary Geometry in Hyperbolic Space. Werner Frenchel
- de Gruyter, Berln, 1989.
[20] Circular Processing of the hue variable: a Particular Trait
of Colour Image Procesing. A. Restrepo, C.Rodriguez and
REFERENCIAS
151
C.Vejarano. VISSAP-Second International Conference on Computer Vision Theory and Applications. Barcelona, Spain, March
2007
[21] Colour spaces.A. Restrepo, C.Rodriguez and C.Vejarano.
VISSAP-First International Conference on Computer Vision
Theory and Applications. Setubal, Portugal, 2006.
[22] Directional Statistics. K.V. Mardia and P.E. Jupp. Wiley,
Chichester, 2000.
[23] Color Appearance Models. Mark D. Fairchild. Wiley, 2005.
Captulo 1
Seales y sistemas
1.0.
INTRODUCCIN
152
1.0. INTRODUCCIN
153
154
Definicin 24. Una seal discreta es una funcin cuyo dominio es un intervalo
de ZN . El rango de una seal discreta es un conjunto de nmeros; por ejemplo, el
conjunto de los nmeros complejos, o el intervalo /0, 255/.
Definicin 25. Un subconjunto W de RN es denso en ninguna parte si para
cada punto de W hay in intervalo de RN que contiene al punto y a ningn otro
elemento de W.
Definicin 26. Una seal continua en sentido laxo es una seal continua que se ha
dejado sin definir en un subconjunto denso en ninguna parte de su dominio.
La ingeniera es en parte el arte de asignar nmeros a magnitudes fsicas con
el objeto de manejar los diferentes aspectos en un diseo. As como los fsicos,
los ingenieros cuantificamos las variables con que trabajamos, lo que nos permite
usar herramientas matemticas y poder hablar objetivamente de los resultados
esperados de un diseo. Esta actividad, que podramos llamar aritmetizar, va
acompaada de aquella de matematizar, donde modelamos matematicamente
los fenomenos fisicos.
Ejemplo. En la figura 1.1 indicamos la funcin de voltaje v definida en el
intervalo (t1 , t2 ), los valores del voltaje son nmeros reales: v : (t1 , t2 ) R1
Seguramente, no cualquier funcin es una seal que provenga de algn
fenmeno fsico, pero la funcin matemtica es un poderoso modelo de seal.
Muchas de las herramientas usadas en ingeniera para el anlisis de seales,
estn basadas en el clculo y en otras disciplinas matemticas.
En un principio, en ingeniera de comunicaciones, las seales consideradas
eran generalmente continuas: seales de voltaje o corriente, en el dominio del
tiempo; por lo general eran seales de radio. La idea de poder transmitir palabras a distancia era muy atractiva. As como el estudio de la electricidad se
impuls con los experimentos de Volta con patas de ranas, probablemente la
idea de usar un cable coaxial para la transmisin de informacin a distancia se
1.0. INTRODUCCIN
155
inspir en los axones largos y mielinados del sistema nervioso de los mamferos. El problema que an persiste con un cable coaxial (o de fibra ptica), es
la atenuacin que sufre la seal con la distancia; as, resulta muy importante
poder amplificar una seal. La tecnologa de tubos de vaco con electrodos
conectados a terminales externos permiti efectuar tal amplificacin. Esta misma tecnologa permiti el desarrollo del tubo de rayos catdicos, que es con lo
que se construan las pantallas de TV y los osciloscopios anteriormente. A finales de la dcada de los 40 del s.XX, se desarroll un dispositivo amplificador
semiconductor que era pequeo, confiable y permita integracin (es decir, la
fabricacin de un circuito electrnico muy denso: (muy denso: miles de componentes en un milmetro cuadrado, as como la posibilidad de producirlos en
gran cantidad): este dispositivo es el transistor bipolar, inicialmente construido con germanio o silicio y despues con arseniuro de galio. Probablemente
en un futuro se usen las seales luminosas tanto como las electrnicas en los
circuitos integrados.
La tecnologa de las comunicaciones y del control para los vuelos Apolo (en
la dcada de los 60)fue anloga en gran medida (aunque usaron computadores
rudimentarios al final). En los sesentas, el anlisis de seales discretas o series
de tiempo era ms del dominio de los estadsticos que de los ingenieros.
A partir de los setentas, la tecnologa digital entr en auge; se desarrollaron
varias familias de circuitos integrados que implementaban funciones lgicas,
se invent el microprocesador, aparecieron las calculadoras electrnicas de
bolsillo (de 4 operaciones), las memorias bajaron de precio y aumentaron en
capacidad, etc. (La industria electrnica japonesa que hasta el momento era de
calidad regular lleg a sobrepasar el desempeo de la industria estadounidense
en algunos aspectos). Paralelamente, la teora matemtica bsica que sera usada para el anlisis y el tratamiento lineal de seales digitales, se sofistic
en algn grado interpretando algunos resultados aplicables y desarrollando
otros como por ejemplo la transformada rpida de Fourier. Hoy en da, para
156
tasas de muestreo por debajo de, aproximadamente, 109 Hz, la tecnologa digital
supera en muchos aspectos a la tecnologa anloga; el caso ms patente es el
tratamiento de seales de audio donde el uso de tcnicas digitales permite
una gran versatilidad en la implementacin de sistemas para el tratamiento
de las seales. Desde los ochentas, compaas como Texas Instruments producen procesadores (DSP o digital signal processor) con arquitectura orientada
al tratamiento de seales.Desde la dcada de los noventas ha habido un boom
de las comunicaciones, que ha sido posible en parte por la tecnologa para el
tratamiento digital de seales.
1.1.
157
158
dos funciones pueden tener el mismo grfico y distinto rango; como acabamos
de decir, se asume que el rango es el conjunto de valores que la seal puede
tomar, fsicamente. por ejemplo, la seal almacenada en un disco compacto
es digital ya que cada componente est codificado con 16 bits. La seal del
ejemplo anterior es una seal digital, as como la del ejemplo siguiente:
Ejemplo. f : {1, 2, 3, 4, 5} /0, 15/, dada por f (1) = 10, f (2) = 0, f (3) = 1,
159
1 si t 0
u(t) =
0 si t < 0
La seal escaln discreta u : Z C est dada por
1 si n 0
un =
0 si n < 0
Tambin, la seal impulso discreta, : Z C est dada por
1 si n = 0
n =
0 si n , 0
Por razones que explicamos ms adelante, no definiremos la seal impulso
continua como una funcin, sino como un funcional.
La seal escaln continua se conoce tambin, especialmente en matemticas, como funcin de Heaviside, en honor a Oliver Heaviside quien not su
importancia en la teora de sistemas lineales.
Ejercicio. Puede ser digital una seal continua?
Ejemplo. Sea v : R1 N donde v(x) = parte entera(x)
160
Tamao
Definimos el tamao de un conjunto a continuacin. Si el conjunto es contable, su tamao est dado por su cardinalidad. Si el conjunto es un subconjunto
de RN , su tamao est dado por la correspondiente medida (de Lebesgue): su
longitud, su rea o su volumen, etc.
Ejemplo. El tamao del intervalo /4, 9/ es 6.
El tamao de Z es infinito.
El tamao del producto cartesiano /0, 1023//0, 1023//subsetZ/timesZ es 1.048.576
(o un mega en la jerga de los programadores).
El tamao del intervalo (0, 10)/subsetR es 10, as como tambin lo es el del
intervalo [0, 10].
El tamao de R1 es infinito.
El tamao del producto cartesiano [2, 5] [0, 2]/subsetR/timesR es 4.
Soporte
El soporte de una funcin es el subconjunto del dominio donde la funcin
toma valores diferentes de cero.
Ejemplo. El soporte de la seal escaln continua u : R C, es el intervalo
161
[0, ).
El soporte de la seal seno sen : R1 C, es R1 Z.
El soporte de la seal escaln discreta, u : Z C, es el conjunto de los nmeros
naturales.
El soporte de la seal pulso discreta p : Z C, pn = un un10 es /0, 9/.
Longitud y duracin
La longitud de una seal est dada por el tamao de su dominio. La
duracin de una seal est dada por el tamao del intervalo ms pequeo
que contenga su soporte.
Ejemplo. A continuacin se muestra una seal digital de longitud finita. Sea
s : /0, 3/ {0, 2} dada por s = [0, 2, 0, 2], s(0) = 0, s(1) = 2, s(0) = 0, s(1) = 2.
Figura 1.4: Una funcin del conjunto 4 al conjunto de los nmeros reales
162
Ejercicio. Grafique la seal g(t) del ejemplo anterior para t (1, 2).
1.2.
Funcionales
1.2. FUNCIONALES
163
1.2.1.
Funcionales de CN a C
164
1.2.2.
1.2.3.
1.2.4.
El espacio dual de CN
1.2. FUNCIONALES
1.2.5.
165
Funcionales de CZ a C
1.2.6.
1.2.7.
1.2.8.
166
El subconjunto de CN dado por las seales con promedio finito, lo denotamos como A. Dado que la suma {an +bn } de dos seales {an } y {bn } de promedio
finito tiene promedio finito, y que el producto {w an }porunescalarwtienepromedio f inito, seconcluyequeAesunsubespacio,
El funcional promedio : A C, dado por (s) = s, es un funcional en A.
1.2.9.
1.2.10.
Dado que la suma de dos seales con una cantidad finita de componentes
no nulas tiene a su vez una cantidad finita de componentes no nulas, tenemos
nuevamente un subespacio lineal de seales, que llamamos B Un ejemplos de
funcional definido en B es el funcional suma. Note que el promedio de toda
seal en B es cero.
1.2.11.
El dual de c0 es l1 .
1.2.12.
Los espacios lp
Para p [0, ), se tiene que (x+ y)p xp + yp . Para cada p [0, ), definimos
entonces el espacio lp como el conjunto de las sucesiones bilaterales complejas
P
p
{sn } tales que
n= |sn | R. Adems, dado que |x + y| |x| + |y| , ponemos
1.3. SISTEMAS
167
l = {s CZ : sup{|sn | : n Z} R}
Los espacios lp , as como sus espacios duales, son muy interesantes.
1.2.13.
p
p1
El espacio dual de l1 es l .
El espacio dual de l contiene a l1 como subconjunto propio.
1.2.14.
1.3.
Funcionales de CR a C
Sistemas
168
1.3. SISTEMAS
169
1.3.1.
Espacio vectorial
170
Ejercicio. 0 v = .
1.3.2.
Transformacin lineal
Para cada sistema lineal, continuo o discreto, se tiene por lo tanto que la respuesta
a la seal cero es la seal cero.
Ejemplo. Considere el espacio vectorial V de las funciones continuas h : [0, 1]
C, con las operaciones de suma de funciones, sobre el cuerpo de los complejos.
Considere el conjunto de los funcionales lineales, definido sobre este espacio.
Como se enuncia en el teorema de Riesz, para cada funcional lineal : V C,
de este tipo, hay una funcin g : R C, de variacin acotada, tal que se
puede representar como una integral de Riemann-Stieltjes:
1
Z
(h) =
h(t)dg(t)
0
Ejercicio. En relacin con el ejemplo anterior, si (h) = h(0), para cada h, Cmo
1.3. SISTEMAS
1.3.3.
171
1.3.4.
Normas Elep
172
P
P
cada una de las series
|sn | y
|sn | convergen (a nmeros reales). Sea l1 el
n=
n=0
1
X
|sn | +
n=
n=0
n=1
n=0
|sn | +
|sn |
|sn |
Ejercicio. Muestre que efectivamente se tiene una norma (que cumple con
1 + 2 = 3.
1
P
n=
2n +
P
n=0
2n =
X
X
2
|sn | +
|s2n |
ksk2 =
n=
n=0
P
1
n=1
1.3. SISTEMAS
P
n=1
P
n=1
173
1
n2
(1)n+1
n
Ms generalmente, definimos
1
X p p
kskp :=
|s|
n=
Eles maysculas
Para seales continuas, tambin se define la norma Ele-uno y se denota con
ele mayscula. Sea h : R1 C una seal continua. La norma L-1 de h est
dada por:
Z
khk1 =
|h(t)|dt
khk2 =
|h(t)|2 dt
Originalmente en ingeniera elctrica, las seales consideradas eran funciones de voltaje o de corriente casi siempre. Dado que la potencia instantnea
que disipa una resistencia es proporcional al cuadrado del voltaje entre sus
174
1
t
k f k1 y k f k2 . Diga si f L1 y/o f L2 .
Ejercicio. Sea f la funcin dada por
t , si t (, 1) (1, )
f (t) =
0, si t [1, 1]
Encuentre k f k1 y k f k2 . Diga si f L1 y/o f L1 .
1.3. SISTEMAS
175
1,
f (t) =
sen t ,
t
si t = 0
si t , 0
R1
1
t
1
t
1
t2
R
0
khk1 =
|h(t)|dt
176
1
t2
1.3. SISTEMAS
177
178
khk1 =
|h2 (t)|dt
1.4.
La norma ele-infinito
179
Por lo tanto, una seal es acotada si y slo si su norma ele-infinito es finita. Por
ejemplo, la norma de la seal continua cos t es 1, mientras que la seal discreta
2n no es acotada.
Ejercicio. Encuentre la norma L de la seal g(t) = 1 2|t| .
Ejercicio. Sea F : R C dada por F() = e j + a + j. Encuentre kFk .
180
1.5.
TA
La variacin es un funcional. Definimos la variacin de una seal discreta
P
{sn } como Var(s) = nZ ksn sn1 k. Es interesante anotar que el filtro mediana
(definido ms adelante) nunca aumenta la variacin de una seal, lo que puede
considerarse como una propiedad alisadora (smoothing) del filtro.
1.6.
TA
Fue deinida en el captulo anterior, al hablar de la integral de RiemannStieltjes.
1.7.
181
Variacin promedio
FALTA
1.8.
Problemas
R sen(t)
4. Muestre que ( t )dt = /2.
0
5. Muestre que
sen(t)
t |dt
182
1
.
1+t2
11. D un ejemplo de una seal con norma L1 finita, pero con norma L2
infinita.
12. D un ejemplo de una seal con norma L2 finita, pero con norma L1
infinita.
13. Calcule la norma L de s(t) = arctan (x) .
14. Con ejemplos, muestre que ni L1 es subconjunto de L2 , ni L2 de L1 .
15. Demuestre que L2 [0, 2) L1 [0, 2).
16. Demuestre que l1 l2 .
17. Encuentre una seal continua que est el L1 pero no en L2 , y viceversa.
18. Muestre que cada seal discreta s : Z C sumable en magnitud, es
sumable en magnitud al cuadrado.
1.8. PROBLEMAS
183
1
.
1+t2
Referencias
[1] Geometry of Complex Numbers. Hans Schwerdteger. Dover,
N.Y.,1979.
[2] Elementary Theory of Analytic Functions of One or Several
Complex Variables. Henri Cartan. Dover, N.Y., 1995.
[3] Mathematical Analysis. Tom M. Apostol. Adyson-Wesley,
Reading MA, 1974.
[4] Signals and Systems. A. V. Oppenheim, A. S. Willsky y I. T.
Young. Prentice Hall, Londres, 1983.
[5] Fourier and Wavelet Analysis. Bachman, Narici and Beckenstein. Springer, N.Y., 2000.
[6] Functional Analysis. F. riesz y B. Sz-Nagy. dover, N.Y., 1990.
[7] Fuctional Analysis. G. Bachman y L. Narici. Academic Press,
N.Y., 1966.
184
Captulo 2
Convolucin y sistemas de
convolucin
2.0.
INTRODUCCIN
La convolucin se puede ver como una operacin entre seales; as, dadas
dos seales, resulta otra seal. La convolucin es una herramienta de utilidad
en probabilidad, anlisis de sistemas, anlisis de Fourier y otras reas; su utilidad es ms que todo terica. La solucin de muchas ecuaciones diferenciales
(y de diferencia) lineales, no homogneas y con coeficientes constantes, en
funcin del trmino no homogneo 1 , es una convolucin continua o discreta,
respectivamente.
Dada la fcil representabilidad en el dominio de la frecuencia de Fourier de
los sistemas de convolucin, la convolucin es un concepto muy til y popular
1 Los
185
186
en tratamiento de seales.
La convolucin tambin es un concepto importante en teora de sistemas,
en ecuaciones diferenciales y en ecuaciones de diferencia. Consideremos esto
desde el punto de vista de la ingeniera elctrica. Un circuito puede verse como
un sistema que transforma una seal de entrada en una seal de salida. Si el
circuito est compuesto por resistencias, bobinas y condensadores lineales e
invariantes, adems de amplificadores lineales, el circuito se puede modelar
con una ecuacin diferencial, lineal y de coeficientes constantes. Es decir que
para cada corriente o voltaje en el circuito hay una ecuacin integro-diferencial
en esta variable, con un trmino no homogneo que depende de la seal de
entrada al circuito y, si se quiere, de las condiciones iniciales del circuito. En
realidad, en este curso no asumiremos condiciones iniciales, debido a que
asumimos que el circuito ha estado funcionando siempre, desde t = 1 .
Vale la pena resaltar que si derivamos una ecuacin integro-diferencial al
querer convertirla en una ecuacin diferencial, perdemos informacin; adems,
es posible que el trmino no homogneo no tenga derivadas del orden requerido.
Tanto el voltaje entre los terminales de un condensador, como la corriente
por una bobina son integrales, de la corriente y del voltaje, repectivamente,
y por esto asumimos que son funciones continuas. Es decir, que si para un
Rt
intervalo [a, b] se tiene que la integral a(t) := t b()d existe como un nmero
0
real, para toda t [a, b] entonces a(t) es continua en (a, b). Adems, asumimos
que la corriente por una inductancia y el voltaje sobre un condensador, por ser
integrales, adems de ser funciones continuas, cuentan con derivadas izquierda y derecha en cada punto t.
Adndum. Note que
1 Esto
R
1
d
0
est relacionado con el hecho de que las transformadas, de Laplace y de Fourier que
consideramos, son bilaterales.
2.0. INTRODUCCIN
187
(2.1)
s(t) =
e(t )dg()
(2.2)
188
R
s(t)
L
(2.3)
R
s(t)
L
(2.4)
en el primer caso, la solucin est dada por s(t) = keRt/L y como requerimos que
la solucin sea acotada, ponemos k = 0. As, tenemos s(t) = 0 y, en particular,
2.0. INTRODUCCIN
189
g(t) =
pasa altas?
CONVOLUCION DISCRETA Y ECUACIONES DE DIFERENCIA
Por otra parte, la convolucin discreta es la solucin de las ecuaciones de
diferencia lineales y con coeficientes constantes, en trminos del trmino no
homogneo de la ecuacin. As por ejemplo, para la ecuacin de diferencia, sn =
0,5sn1 + en , donde nuevamente asumimos que e es la entrada a algn sistema
discreto gobernado por esta ecuacin y que s es la salida correspondiente,
resulta que, asumiendo causalidad, la salida s para una entrada e arbitraria
P
est dada por,
enk hk . Donde h es la solucin de la ecuacin cuando la
k=
190
2.1.
Dado que a veces una serie es ms intuitiva que una integral, quizs la
convolucin discreta sea ms fcil de visualizar que la continua; por esto
comenzamos con ella. Sean f y g seales discretas con dominio Z; la convolucin de f y g, que se denota f g, es tambin una seal discreta. Sea h = f g;
P
para cada n, hn est dada por, hn =
fk gnk .
k=
Ejercicio. Diga si
P
k=
fk gnk =
k=
fk gNk , resulta
k=
n, n / 10, 10/
fn =
0, n < /1, 4/
gn =
1, n /1, 4/
por,
1/3, n {1, 0, 1}
hn =
0,
n < {1, 0, 1}
191
k
2k
g2k
fk
-2
4
1
-2
-1
3
1
-1
0
2
1
0
1
1
1
1
2
0
0
2
3
-1
0
3
4
-2
0
4
5
-3
0
5
192
y0 = 1/3(x1 + x0 + x1 ) = 2/3
y1 = 1/3(x0 + x1 + x2 ) = 1
y2 = 1
y2 = 0
...
y, en general,
0,
1/3,
yn =
2/3,
1,
n 2
n = 1
n=0
n1
Ejercicio. Para el ejemplo anterior, demuestre que, en general, para una entrada
x, yn = 1/3(xn1 + xn + xn+1 )
193
X
=
2k xnk
n=0
=xh
con hn = 2n un
2.2.
una seal escaln. sea s = e h. Las figuras 2.6 y 2.7 ilustran el clculo de s(1).
La funcin e(1 ), en el eje , es una versin reflejada y desplazada de e(),
194
195
0,
t 3
t + 3, 3 t 1
s(t) =
2,
1 t 1
3 t, 1 t 3
0,
t3
(vea la figura 2.8). As decimos que la convolucin de dos pulsos es un trapecio.
Similarmente, la convolucin de un pulso consigo mismo, es un triangulo.
Ejercicio. Calcule y grafique h h, para la h definida en el ejemplo anterior.
196
Figura 2.8: g = f h
197
2.2.1.
Curvgrafos (Splines)
Definimos el curvgrafo de orden cero, como s0 (t) = u(t + 0,5) u(t 0,5),
El spline de orden 1 como s1 (t) = [s0 s0 ](t) y en general, en forma inductiva,
sn (t) = [sn1 s0 ](t).
Ejercicio. Grafique s0 (t), s1 (t), s2 (t) y s3 (t).
Ejercicio. Muestre que el lmite de sn (t), cuando n tiende a infinito, es una
2.2.2.
Convolucin de Stieltjes
e(t )dg()
Z
e(t )dg() =
e(t )h()d
198
La convolucin de Stieltjes es til para caracterizar los sistemas de convolucin que tengan respuesta escaln que sea una funcin discontinua (seal con
saltos), ya que evita el uso de herramientas de anlisis funcional tales como
los impulsos de Dirac.
Ejemplo. Sea g(t) la seal mostrada en la figura 2.9 y suponga que quere-
199
Z
+
Z
+
Z
= 0 + sen(t) +
donde z = 2e j/4 e j .
Ejercicio. Calcule f dg si f (x) = x3 y g = esta dada por, g(x) = 1/2, x [1, 1]
(x dg) + (w dg)
2.3.
200
2.3.1.
Respuesta escaln
u(t )dg() +
= 0g() 1g()
= g()
201
u(t )dg() =
g()du(t ) g()
= g(t) g()
Y si g() = 0 (que es lo que asumiremos de ahora en adelante con respecto
a la funcin g (la pepa) de los sistemas de convolucin continuos, a menos
que se diga lo contrario) se tiene que g es la respuesta escaln del sistema.
El resultado obtenido muestra que g es la respuesta escaln de un sistema
que tenga relacin entrada/salida dada por la ecuacin 2.5. Tambin, dadas las
propiedades de superposicin y homogeneidad de la operacin de convolucin indicadas en los ejercicios finales de la seccin 2.2, podemos concluir que
los sistemas de convolucin son lineales. Una vez que definamos invarianza, veremos que los sistemas de convolucin tambin son invariantes. Sinemabrgo,
hay sistemas lineales e invariantes que no son de convolucin.
Muchos sistemas de inters en ingeniera electrnica son sistemas de convolucin; por ejemplo, los circuitos de resistencias, inductancias y capacitancias lineales, ya que stos se modelan con ecuaciones diferenciales lineales con
coeficientes constantes y, como mostraremos en la seccin 2.9, la solucin de
estas ecuaciones es una convolucin de Stieltjes.
Ejemplo. El circuito de la figura 2.11 es un sistema de convolucin con respues-
202
2.3.2.
203
2.3.3.
204
2.3.4.
=
e0 (t )dg()
0
e(t )dg()
s(t) =
= g()e() g()e()
g()de(t ),
g(t )de(),
= g()e() g()e()
Z
=
g(t )de()
g(t )de()] =
d
g(t )de()
dt
205
s0 (t) =
e(t )dg0 ()
ciones?
Ejercicio. Hay casos, como cuando e(t) = cos t, en que e() , 0, qu pasa en
estos casos?
2.3.5.
2.3.6.
206
u()h(t )d
g(t) = [u h](t) =
u()h(t )d
0
h(x)(dx) =
t
h(x)(dx)
1 1
n, si t [ 2n , 2n ]
dn (t) =
0, si t R1 [ 1 , 1 ]
2n 2n
207
0, t , 0
(t) =
, t = 0
R
R
Para cada n, dn = 1 y, por lo tanto, lm dn = 1; sin embargo. la funcin
(t) tiene integral (de Riemann o de Lebesque) nula, es decir,
Z
(t)d(t) = 0
R
R
Como se ve, lm dn (t)dt , lm dn (t)dt. La razn por la cual es comn
el uso de la notacin para integrales para denotar funcionales lineales es que
la integral es lineal, pero aparte de esta propiedad de linealidad, el funcional
lineal Dt con dominio el conjunto de funciones contnuas f : R1 R1 , y rango
el conjunto de los reales, dado por Dt ( f ) = f (t) es en realidad diferente de la
R
integral f ()(t )d.
Aunque nemotcnicamente sea conveniente el uso de distribuciones de
Dirac con notacin de integral, es probable que no sean convenientes en un
estudio inicial de los sistemas de convolucin continuos por parte de un futuro
208
Z
e(t )g ()d + 0 +
Z
e(t )g ()d = 0 + 0 +
e(t )e( )d
Note que aunque, eh = he, se tiene que, siempre y cuando g() = e() = 0,
Z
Z
e(t )dg() =
e()dg(t )
sabemos, la respuesta escaln de este circuito est dada por g(t) = et/RC u(t),
la cual tiene un salto en t = 0. Si queremos obtener la salida de este circuito
cuando la entrada es e(t) = sen t, debemos evaluar la siguiente integral:
Z
s(t) =
sen(t )dg()
209
e(t )dg().
Figura 2.13: La respuesta escaln de un sistema de convolucin. Note que as como un RC pasa
altas, la respuesta escaln de este sistema no es diferenciable.
2.3.7.
una de ellas tiene rea igual a uno. Este resultado se puede inferir del hecho que la funcin de
densidad de probabilidad de la suma de dos variables aleatorias independientes es la convolucin
de las densidades de las variables que se suman.
210
el producto de la carga de e y
Z
dg(), es decir,
"Z
s(t)dt =
# "Z
#
dg(t)
e(t)dt
s(t)dt =
e(t)dt
h(t)dt
tenemos
s(t)dt =
f (u)du
dg
h()d
211
forma u(t + 1/2) u(t 1/2), la salida tiene una integral igual a la integral de
la respuesta impulso del sistema.
2.4.
La transformada de Hilbert
t, t , 0
h(t) =
0, t = 0
Por lo tanto, es la respuesta de un sistema de convolucin con respuesta escaln
h(t) y con funcin caracterstica 1t .
2.5.
ek hnk
k=
0, si n , 0,
n =
1, si n = 0
212
La respuesta impulso de un sistema es la respuesta del sistema cuando la seal que se aplica a la entrada es una seal impulso. A continuacin
se mostrar que la respuesta impulso de un sistema de convolucin discreto es precisamente su funcin caracterstica h. Si la entrada al sistema es la
seal impulso , entonces la salida evaluada en la instancia n est dada por,
P
P
k hnk =
nk hk = hn ya que nk = 1 cuando k = n, y 0 en caso
k=
k=
contrario.
Ejercicio. Se tiene un sistema promedio mvil, el cual para cada entrada e y
P
una convolucin: rn =
ek hnk
k=
2.6.
213
2.6.1.
2.6.2.
2.7.
214
1
= (1 et/rC )u(t)
rsC + 1
215
tiene funcin de transferencia H(s) = 1, muestre que pasa altas = identidad pasa bajas.
Ejercicio. Encuentre la funcin de transferencia (la relacin [volt. salida]/[volt.
entrada]) para un circuito similar al de la figura 2.15, en el dominio j. Grafique
la magnitud y la fase. Considere los siguientes casos:
Inductancia, capacitancia.
Resistencia, paralelo (capacitancia, inductancia).
Paralelo (capacitancia, inductancia), resistencia.
Serie (capacitancia, inductancia), resistencia.
Resistencia, serie (capacitancia, inductancia).
2.8.
Correlacin determinstica
La correlacin determinstica 3 de dos seales se puede definir en trminos de una convolucin. La correlacin determinstica de las seales
3 En
contraste con la correlacin estocstica est dada por la esperanza E[XY] de dos variables
aleatorias o de dos procesos estocsticos.
216
X
rn =
sk tkn
k=
por lo tanto, rn = sn tn .
Ejercicio. Es cierto que corrdet(g, h) = corrdet(h, g)? Demuestre o d un con-
traejemplo.
La correlacin determinstica es un estimador usado en estimacin espectral. La densidad espectral de potencia de un proceso aleatorio estacionario est
dada por la transformada de Fourier de su funcin de autocorrelacin. Norbert
Wiener del MIT, con una amplia gama de intereses, hizo aportes importantes
en esta rea de a caracterizacin espectral de procesos aleatorios y, en general,
en anlisis de Fourier.
2.9.
(2.6)
2.10. PROBLEMAS
217
asumiendo g con derivadas hasta de orden n, veamos qu sucede si reemplazamos la funcin g = x dg en el lado izquierdo de la ecuacin 2.6:
Z
+ an1 (x dg)
(n)
(x dg)
(n1)
+ + a1 (x dg) + a1
(x dg)
dg
Z
= x [dg
(n)
+ an1 dg
(n1)
+ + a1 dg + a1
dg]
= x du
=x
por lo tanto, g es solucin de la ecuacin 2.6.
La diferenciabilidad o integrabilidad de g las garantizamos en este caso
por ser estas derivadas e integrales variables fsicas.
2.10.
Problemas
218
P
si se puede representar con la ecuacin de la forma, rn =
ek hnk .
k=
r(t)dt.
2.10. PROBLEMAS
219
220
2.10. PROBLEMAS
221
16. Muestre que la convolucin iterada del curvgrafo de orden creo, consigo
mismo, en el lmite, es una campana de gaus.
17. Encuentre la funcin de transferencia con respuesta escaln g(t) = (1
et )u(t). Exprese la salida como una suma de senusoides si la entrada es
cos2 (t).
18. Se tiene un sistema de convolucin dicreto con respuesta impulso dada
por hn = Un 2|n| .
a) Encuentre la funcin de transferencia H() y muestre que es real y
par.
20
222
X
n=
yn senc(
t nT
)
T
Referencias
[1] Mathematical Analysis. Tom M. Apostol. Addison-Wesley,
Reading MA, 1974.
[2] Signals and Systems. A. V. Oppenheim, A. S. Willsky y I. T.
Young. Prentice Hall, Londres, 1983.
[3] Ecuaciones Diferenciales. D. L. Kreider, R. G. Kuller y D. R.
Ostberg. Fondo Educativo Interamericano, Bogot , 1973.
223
Captulo 3
Representacin de seales
peridicas en el dominio de la
frecuencia (de Fourier)
3.0.
Seales peridicas
Algunas de las seales que se originan en procesos fsicos son aproximadamente peridicas durante algn intervalo de tiempo. Por ejemplo la seal de
voltaje tomada por un micrfono de un violn tocando un La. A los sonidos
graves les corresponden frecuencias ms bajas que a los sonidos ms agudos.
Parece ser que el sistema auditivo humano reconoce patrones que son aproximadamente peridicos durante periodos cortos de tiempo. Por otra parte,
las funciones peridicas (principalmente sinusoidales) han jugado un papel
importante en fsica y en matemticas, y su teora est bastante desarrollada.
Por estas y otras razones es de utilidad el estudio de las seales peridicas en
224
225
ingeniera.
Tambin, es importante diferenciar entre la frecuencia con que un proceso
peridico se repite y las frecuencias que un anlisis de Fourier produce. As,
definiremos frecuencia de Fourier como la frecuencia de una sinusoide o
una exponencial compleja. como veremos en el captulo 4, las exponenciales
complejas son las eigenfunciones, o funciones propias de los sistemas de convolucin. Aunque la frecuencia de Fourier ha jugado un papel fundamental en la
teora de sistemas lineales y en el anlisis de seales , nuevas tcnicas como
las onditas (ondelettes, wavelets), permiten la consideracin de otros conjuntos
de seales bsicas.
3.0.1.
Definiciones
226
1, si n es mltiplo de 3
rn =
0, si n no es mltiplo de 3
encuentre el periodo.
227
cuentre el periodo.
3.0.2.
Frecuencia
Si se tiene una seal peridica, continua con periodo T, o discreta con periodo N, mayor que cero, su frecuencia est dada por T1 o N1 , respectivamente.
La frecuencia de una seal continua est en [0, ], mientras que la frecuencia
de una seal discreta es un nmero en el intervalo [0, 1]. Por cierta razn que
explicamos a continuacin, para seales discretas, la frecuencia 0 es equivalente a la frecuencia 1, mientras que para frecuencias de seales contnuas
asumiremos que 0 es equivalente a .1
Por convencin, decimos que la frecuencia de una seal constante (discreta o continua) es cero. As, las seales discretas con periodo 1 (las discretas
constantes) tienen 2 frecuencias: uno y cero. Las seales constantes (y otras)
contnuas tambin tienen 2 frecuencias: 0 y .
Ejercicio. D un ejemplo de una seal continua peridica no constante, con
periodo cero.
Como se define ms adelante, para exponenciales complejas y para sinusoides, la frecuencia angular de una sinusoide o de una exponencial compleja
con frecuencia f , es = 2 f .
Ejercicio. Sea s una seal discreta con periodo N, con N > 1. Sea r un nmero
228
3.0.3.
Unidades
3.1.
3.1.1.
229
Caso discreto
3.1.2.
Notacin
1, n 3Z
xn =
0, de lo contrario
1, n 4Z
yn =
0, de lo contrario
{. . . 0, 1, 2, 1, 0, 3, 0, 1, 2, 1, 0, 3, 0, . . . }, y = {. . . 0, 0, 2, 0, 0, 2, 0, 0, 2 . . . } Cul
es el periodo de x + y?
230
Figura 3.3: Una seal con periodo 3, otra con periodo 4 y su suma
Ejercicio. Caracterizacin del periodo de la suma de seales peridicas. Sean x
3.1.3.
Caso continuo
231
232
cero?
Conjetura 1. La suma de dos seales peridicas contnuas es peridica si y solo si el
cociente de sus perodos es un nmero racional.
Aunque la propiedad de semi periodicidad de una seal ha sido definida
rigurosamente y es intuitiva, no ha sido explotada en ingeniera. Es posible
observar diagramas de tiempo correspondientes a sonidos naturales tales
como una nota musical siendo tocada por un violn o cantada por una persona
y se observa que la seal se repite pero slo aproximadamente. En teora de caos
se estudian sistemas dinmicos no lineales de orden tercero o mayor y se ha
observado que algunos de ellos producen seales que se repiten aproximadamente. Por otra parte, cuando se audifica una seal (estrictamente) peridica,
el sonido no es agradable, suena como un pito.
El estudio de los sistemas dinmicos ha recibido un impulso reciente dado el inters en la teora del caos. All, las seales semi peridicas son ms
comunes.
Ejercicio. Demuestre que la suma de dos sinusoides es peridica si y slo si el
3.2.
233
3.2.1.
Caso continuo
3.2.2.
Caso discreto
234
3.3.
Una clase importante de seales peridicas est dada por las sinusoides
y las exponenciales complejas. Estas son las funciones propias de los sistemas
lineales e invariantes.
3.3.1.
Caso continuo
3.3.2.
Caso discreto
235
236
237
238
3.4.
3.4.1.
Series de Fourier
P
se puede representar con una serie de Fourier: f (t) =
cn e jnt donde los
n=
f (t)e jnt dt
0
239
X
f (t) =
cn e jn0 t
Si f es continua en t0 , entonces,
n=
donde
Z
1 T
f (t)ej0 nt dt
T 0
0 se conoce como la frecuencia angular fundamental de f y est dada
por: 0 = 2/T (la frecuencia fundamental es 1/T). Los mltiplos de la frecuencia fundamental se conocen como frecuencias armnicas. La fundamental
es la primera armnica. Cuando un conjunto de frecuencias tienen mnimo
comn mltiplo, se dice que estn armnicamente relacionadas. El valor del
coeficiente c0 , el cual est dado por el promedio de la seal sobre un periodo,
se conoce tambin como el nivel DC o nivel promedio de la seal.
La representacin en series de Fourier de una seal peridica permite su
representacin en el dominio de la frecuencia de Fourier.
cn =
P
este caso) de Fourier
cn e jn0 t .
n=
Ejercicio. Sea s(t) una onda cuadrada, como se muestra en la figura 3.9. En-
240
241
Ejercicio. Suponga que se tiene una seal cuadrada de voltaje, positiva, con
uniformemente.
3.4.2.
Propiedades de simetra
3.4.3.
Frmula de Parseval
La siguiente identidad, conocida como frmula de Parseval, permite relacionar la potencia promedio de una seal peridica con la norma l 2 de sus
242
coeficientes de Fourier:
1
|cn | =
T
n=
| f (t)|2 dt
1
1
dt =
2
1 + t2
"
#
1 t2
1
+
dt
(1 + t2 )2 1 + t2
243
1
dt = arctan t|
=
2 2
(1 + t )
Z /2
t2
tan2 u
dt
=
sec2 udu =
2
2
2
2
2
(1 + t )
/2 (1 + tan u)
3.4.4.
Series de sinusoides
X
f (t) =
cn e jnt
n=
donde
cn =
1
2
f (t)e jnt dt
X
n=0
an cos nt +
bn sen nt
n=1
donde
R
R
1
1
a0 = c0 = 2
f
(t)dt
y
a
=
f (t) cos(nt)dt y bn =
n
P
n 1. Tambin como, f (t) =
e
n=0 n cos(nt + ).
R
1
f (t) sen(nt)dt,
sucesin impar.
244
3.4.5.
Suponga que f es una onda cuadrada impar. La fundamental y las siguientes 3 armnicas de f son entonces sen t, 1/3 sen 3t, 1/5 sen 5t, 1/7 sen 7t,
recuerde que una onda cuadrada no tiene armnicas pares. Las armnicas se
muestran en la figura 3.13 y las sumas parciales a continuacin en las figuras
3.14 y 3.15.
Sea f una seal continua peridica, con periodo 2. Suponga que f (t) =
R
P
1
jnt
f (t)e jnt dt. Para cada nmero natural positivo
donde, cn = 2
n= cn e
Pm
jnt
m, defina fm as: fm = n=m cn e , entonces,
fm (t) =
m
X
cn e jnt
n=m
=
=
=
=
m Z
1 X
f ()e jn de jnt
2 n=m
Z
m
X
1
f ()
e jn e jnt d
2
n=m
Z
m
X
1
f ()
e jn(t) d
2
n=m
Z
1
f ()SEm (t )d
2
SEm (x) =
Si x es mltiplo de 2 y
SEm = 2m + 1
245
Si x es un mltiplo de 2.
A continuacin se grafica la funcin senc enmascarada de orden 4; note
que es peridica con periodo 2.
Pm
n=m
zn =
zm (1z2m+1 )
1z
3.5.
Suponga que tanto f como g son funciones peridicas con periodo T y con
expansiones en series de Fourier dada por:
f (t) =
X
n=
cn e
jn0 t
g(t) =
X
n=
dn e jn0 t
246
Figura 3.13: Las seales sen t, 1/3 sen 3t, 1/5 sen 5t, 1/7 sen 7t
Figura 3.16: La suma parcial seales sen t + 1/3 sen 3t + 1/5 sen 5t
Figura 3.17: La suma parcial seales sen t + 1/3 sen 3t + 1/5 sen 5t + 1/7 sen 7t
247
248
donde
cn =
1
T
dn =
f (t)e jn0 t dt
0
1
T
g(t)e jn0 t dt
0
Suponga que uno est interesado en saber qu funcin h(t) tiene serie de
Fourier con coeficientes cn dn , es decir,
h(t) =
cn dn e jn0 t
n=
1
=
T
n=
T X
1
=
T
1
=
T
n=
T
g(t s) f (s)ds
0
cn e jnt
n=
si h(t) =
P
n=
g(t) =
X
n=
dn e jnt
249
si r(t) =
P
n=
dn e j2nt , grafique r.
250
3.6.
P
n=
cn e jnt con cn =
Sea Sm (s, t) =
(1)n+1
jn
m
P
n=m
si n , 0 y c0 = 0.
de Fourier.
Sm (s, t) =
m
X
(1)n+1
n=1
jn
m
X
2
e jnt e jnt =
(1)n+1 sen(nt)
n
n=1
1
f (+ ) f ( ) s(t)
2
f (+ )+ f ( )
.
2
y
lm+ g(t) = f (+ ) +
f (+ ) + f ( )
= g()
2
f (+ ) + f ( )
() = g()
2
tx
Sm ( f, x)
f (x+ ) f (x )
2
252
n
X
2
r=1
sen
r
2
n
Z
0
3.7.
253
254
de Fourier de tales seales; estas sumas se conocen tambin como seres discretas
de Fourier.
Las seales discretas peridicas se pueden representar como sumas finitas
de exponenciales complejas, a diferencia de las seales peridicas contnuas
las cuales se representan como series de exponenciales complejas. De hecho, si
el periodo de la seal discreta en cuestin es N, no ms de N trminos en la
suma son necesarios.
Como se ver en el captulo 5, tal expansin permite tambin la representacin en el dominio de la frecuencia de seales discretas de longitud finita, conocida como DFT. La representacin en el dominio de la frecuencia de
seales discretas de longitud infinita se trata en la seccin siguiente; la representacin en el dominio de la frecuencia de estas seales se obtiene por medio
de la Transformada de Fourier de tiempo discreto.
El camino para obtener la expansin de las seales peridicas discretas en
sumas de exponenciales discretas tomado aqu, es algo indirecto; inicialmente,
la idea es obtener cierta biyeccin lineal de CN a CN . Empezamos por mostrar
que el resultado de ciertas sumas de exponenciales complejas es cero.
3.7.1.
Sumas tipo I
N1
P
e j N k . Sea WN = e j N y considere la
k=0
suma
N1
P
k=0
k
WN
,
donde N > 1.
N1
P
k=0
k
WN
= 0 (Vea la figura 3.22.)
N1
P
k=0
k
WN
= 0.
k
Geomtricamente se puede argumentar como sigue. Los nmeros WN
,k
/0, N 1/, estn uniformemente distribuidos sobre la circunferencia de radio 1
255
N1
P
0
k
=
WN
+ WN
0
k+1
N
WN
= WN x, dado que WN
= WN
, se tiene que
k=0
N1
P
k=0
k
WN
; por lo tanto
N1
P
k=0
k=0
N1
P
k=0
k
WN
=
k+1
WN
=
k
WN
= WN x y x = xWN , entonces, x = 0,
ca.
3.7.2.
Sumas tipo II
N1
P
k=0
rk
WN
.
256
entero de N, o no.
Ejercicio. Demuestre que si r es un mltiplo de N (es decir, si r es congruente
N1
P
k=0
rk
WN
=N
N1
P
rk
WN
=
k=0
r(N1)
0
r
{WN
, WN
, . . . , WN
}
N1
P
k
WN
= 0 (Sugerencia:
k=0
(N1)
0
1
{WN
, WN
, . . . , WN }. En
0
r
{WN
, WN
, . . . , WN
r(M1)
0
r
} = {WN
, WN
, . . . , WN
0
r
, WN
, WN
,
r(M1)
. . . , WN
r(M1)
0
r
donde la lista WN
, WN
, . . . , WN
r(M1)
0
r
, . . . , WN
, WN
, . . . , WN
se repite m veces.)
N1
P
k=0
rk
WN
valen N, cuando
ca.
3.7.3.
Una biyeccin de CN a CN
257
(r1)(k1)
N
X
(k1)(i1)
WN
(j1)(k1)
WN
k=1
N
X
0
WN
k=1
=N
Adems, para i , j
Ci j =
N
X
(k1)(i1)
WN
(j1)(k1)
WN
k=1
N
X
(k1)(i j)
WN
k=1
=0
ya que [i j]N , 0, por lo tanto, C es una matriz diagonal con componentes en
la diagonal de valor N y N1 C es la matriz identidad.
(N1)k
0k
1k
Corolario 3. Los N vectores de la forma [WN
, WN
, . . . , WN
forman una base para CN .
258
1
NB
Corolario 4. Dada una seal x = [x0 , x1 , . . . , xN1 ] de duracin finita, sus componentes se pueden expresar as:
xn =
N1
X
Xk e j N nk
k=0
N1
1 X j 2 nk
xn e N
N n=0
N1
P
k=0
3.8.
Problemas
3.8. PROBLEMAS
259
260
8. Encuentre la representacin en series de Fourier (de exponenciales complejas) de la seal f (t) = sen(2t) + cos(3t) 0,5 sen(3t) + sen(5t).
9.
X
n=
cn e jn0 t =
an cos(n0 t) +
n=0
bn sen(n0 t)
n=1
3.8. PROBLEMAS
261
que | f (t)
N
P
ck e | < .
jtk
k=0
n= cn e
jn0 t
. Si s, d
262
sen[(m+1/2)x]
1
T
f (t ) f ()
0
3.8. PROBLEMAS
263
an cos(nt) =
hn e jnt
n=
n=0
a) Grafique f (t).
b) Encuentre a0 , a1 y a2 .
c) grafique g(t) dada por
P
n=
hn e jnt .
264
41. Se tiene una onda peridica de periodo 1, nivel promedio 0,7 y amplitud
P
jnt
pico a pico de 1. Si su serie de Fourier es
n= cn e , encuentre la
potencia promedio de la version truncada de su serie de Fourier, dada
P
por 3n=3 cn e jnt .
2
42. Para la seal continua g : R C dada por : g(t) = sen( 2
15 t) + cos( 12 t).
3.8. PROBLEMAS
265
Referencias
[1] Digital Signal Processing. A. V. Oppenheim y R.W. Shafer.
Prentice Hall, Englewood Cliffs, 1975.
[2] Mathematical Analysis. T.M. Apostol. 2nd Ed. AddisonWesley, Reading, 1974.
[3] Fourier Series. G.P. Tolstov. Dover, N.Y., 1962.
[4] Almost Periodic Functions. H. Bohr. chelsea, N.Y., 1947.
[5] Sobre la periodicidad de sumas de seales discretas peridicas. A. Restrepo y L. Chacn. Memorias del simposio de
Tratamiento de Seales, Imgenes y Visin Artificial, Universidad de los Andes, Bogot, 1995.
[6] On the period of sums of discrete periodic signals. A. Restrepo y L. Chacn. IEEE Signal Processing Letter, vol. 5, No. 7, pp.
164-166, 1998.
[7] Anlisis de Fourier. J. Duoandikoetxea. Addison- Wesley,
Wilmington, 1995.
[8] Fourier Analysis. T.W. Krner. Cambridge University Press,
Cambridge, 1998.
266
Captulo 4
La Transformada de Fourier
de Seales continuas
4.0.
Introduccin
267
268
4.0.1.
Motivacin
P
As, la seal s(t), tiene una serie de Fourier de la forma, s(t) =
cn e jnt ,
n=
R
1
s(t)e jnt dt.
con coeficientes dados por, cn = 2
f (t) =
dn e j T nt
(4.1)
n=
con
1
dn =
T
T/2
f (t)e j T nt dt
(4.2)
T/2
4.0. INTRODUCCIN
269
270
Sea
g(t) = lm f (t)
T
= lm
= lm
dn ej T nt
n=
1
F(n )e jn t
T
n=
X
1
= lm
F(n )ejn t
T
2
n=
g(t) =
1
2
R T/2
T T/2
G()e jt d
4.0. INTRODUCCIN
271
la figura 4.4.
1
dn =
T
=
1
T
T/2
f (t)e jn0 t dt
T/2
Z 1
jn0 t
1
1
2
1
e jn T t
T 1
n=0
T,
i
h 2n
=
2n
1 T e j T e j T , n , 0
T j2n
n=0
T,
=
sen( 2n
)
T 2nT , n , 0
2
2n
= senc (
)
T
T
2
= senc (n)
T
272
4.0.2.
Objetivo
Quizs tenga sentido decir que en la prctica las seales tienen duracin
finita; sin embargo, no es conveniente restringirnos a conjuntos de seales de
duracin finita, la generalidad tiene sus ventajas. En todo caso, muchas seales
continuas, de longitud infinita, de duracin infinita y no peridicas, se pueden
expresar como las integral de una funcin de la forma S()e jt :
Z
1
S()e jt d
(4.3)
s(t) =
2
En estos casos, la funcin S() provee la caracterizacin en el dominio de la
frecuencia de Fourier de la seal s(t). S() se conoce como la transformada de
Fourier de s(t).
Para seales s en L1 , su transformada S() se calcula as:
Z
S() =
s(t)e jt dt
(4.4)
273
que ser en alguna medida lento: hay que empezar por seales en L 1, luego
en L 2, luego L p, luego distribuciones funcionales.
Para seales en L2 se define una transformada de Fourier L 2 que es un
lmite en norma de funciones en L1 L2 - La transformada de una distribucin
est dada por la distribucin aplicada a la transformada de la funcin.
4.1.
integrable.
Si g : R1 C es una seal integrable en magnitud, la siguiente frmula
determina un nmero complejo G() para cada R1 :
Z
G() =
g(t)e jt dt
Z
|g(t)e
jt
|dt =
|g(t)|dt
Note que G() se define como la suma de los siguientes lmites, cuando T
tiende a infinito, de
Z
Z
g(t)e
jt
dt +
g(t)e jt dt
0
(4.5)
274
g(t)e jt dt
(4.6)
1,
s(t) =
0,
t [, ]
t < [, ]
275
e jt dt
1 jt
e | , siempre que , 0
j
e jt e jt
j
2 sen
=
si = 0
e jt dt = 2
resumiendo,
=0
2,
S() =
2 sen , , 0
= 2senc
Esta funcin aparece con relativa frecuencia en teora de seales; la funcin
senc x est definida como senc = senx x , para x , 0, y senc (0) = 1. Como se
puede comprobar, la funcin senc es continua. La funcin senc tambin es
conocida como el ncleo de Fourier continuo.
Ejercicio. Calcule la transformada de la seal continua f dada por f (t) = e|t| .
1
1+t2
est en L1 .
276
R 2
2
et
0
1
.
1+t2
cos(2tx)dt
peridicas tales como las sinusoides y las exponenciales complejas la obtenemos por medio de su representacin como series (o sumas) de Fourier, es
decir, con el conjunto de coeficientes {cn } de la serie, graficado contra = n0
(donde 0 = 2/T y T es un tiempo de repeticin) mejor que contra n. Las
sinusoides no nulas no son seales en L1 ni en L2 . Estas seales no tienen transformada (segn la definimos aqu); la integral correspondiente no converge
para ninguna frecuencia .
Ejercicio. Para qu valores de existe el lmite de
RT
0
cos te jt dt cuando T
tiende a infinito?
De aqu en adelante, hasta la seccin 4.5, solo se considerar la transformada
de Fourier de seales integrables en magnitud. En esta forma, se tendr la
seguridad de que la transformada es tambin una seal cuyo dominio es la
277
bilidad, puede preguntarse en este punto las razones por las cuales la funcin
impulso no se utiliza en la definicin de transformada de Fourier de las seales
sinusoidales; algunas de stas son:
i El autor no considera conveniente su uso en un curso introductorio tal
como Teora de Seales.
ii El impulso, en matemticas, no es ninguna funcin con dominio un
conjunto de nmeros; es una entidad ficticia que tal vez sera el ncleo del
funcional de Dirac, el funcional de Dirac est dado por el funcional que
asigna a cada funcin en un espacio vectorial de funciones con dominio
R, el valor de la funcin cero.
iii La teora matemtica en la cual se definen los impulsos es la teora de las
distribuciones. Esta teora es ms sofisticada que el clculo y cubrirla es
innecesario ya que:
a La funcin escaln, junto con la integral de Riemann-Stieltjes y las
series de Fourier, son suficientes para desarrollar la teora de seales
y sistemas en el contexto del clculo.
b Las series de Fourier permiten una representacin adecuada de
seales peridicas continuas en el dominio de la frecuencia.
4.2.
278
de convolucin y la transformada inversa. La propiedad de modulacin (lineal) de la transformada de Fourier se aprovecha para compartir un canal por
multiplexacin en frecuencia; por ejemplo, ste es el sistema usado par la transmisin de radio AM comercial. La propiedad de convolucin junto con la de
la transformada inversa, permiten el modelaje y el diseo de sistemas lineales
invariantes (e.g. filtros) en el dominio de la frecuencia.
Una tercera razn para la popularidad de la frecuencia de Fourier se encuentra en el campo del procesamiento de audio y de voz: hay una correspondencia entre la frecuencia percibida de un sonido y las frecuencias de Fourier
de sus componentes, especialmente de la fundamental; de hecho, hay quienes
aseguran que los nicos tonos puros son los de tipo sinusoidal. Sin embargo, el odo humano no es muy sensitivo a fase de Fourier: si se escuchan dos
seales peridicas y luego se escuchan nuevamente, esta vez retardando una
con respecto a la otra, es difcil escuchar una diferencia.
Otra razn es que las exponenciales complejas (en general, las exponenciales) son funciones propias de los sistemas lineales invariantes: la respuesta
de un sistema lineal e invariante a una exponencial compleja es una exponencial compleja de la misma frecuencia, multiplicada por una constante (el
valor propio que resulta estar dado por la transformada de Fourier de la
funcin caracterstica del sistema evaluada en la frecuencia de la exponencial
compleja). Cuando un sistema lineal e invariante con funcin caracterstica
real se excita con una sinusoide, la salida es tambin una sinusoide (real) de la
misma frecuencia, posiblemente desfasada y con amplitud diferente (o la seal
cero). En modelos dispersivos de medios de propagacin de ondas, las nicas
seales que no se distorsionan, en cuando a forma de onda, son la sinusoides.
A continuacin se listan algunas de las propiedades de la transformada de
Fourier. Se asume que la seal s en su dominio natural (en contraste con la
seal s en el dominio de la frecuencia) es una funcin integrable en magnitud
(i.e. s L1 ). Estrictamente hablando, se asume que la integral de Fourier se
279
define con una integral de Lebesgue, en vez de con una de Riemann; tambin
, que dos funciones son equivalentes si el subconjunto de su dominio donde
difieren (es decir, el soporte de su diferencia) tiene medida cero.
i) (Teorema de Riemann-Lebesgue.) Si la seals : R1 C es integrable en magnitud, su transformada de Fourier S() tiende a cero a medida que ||
R
g(t)e jt dt, despus
crece: lm S() = lm S() = 0. G() =
se toma G() (G()) = 2G() y entonces
Z
Z
= 0,5 g(t)e jt dt +
g(t)e j(t ) dt
haciendo t0 = t
G() = 0,5
Z
g(t) g(t
jt
) e
dt
|g(t + ) g(t)dt = 0
280
y
Z
G() =
Z
=
g(t)e jt dt
g(t /)e jt dt
y
Z
2G() =
g(t)e jt dt +
Z
g(t /)e jt dt
g(t) + g(t /) e jt dt
resulta
G() 0,5
g(t) + g(t /) e jt dt
dado que |e
281
S( + h) S() =
Z
j(+h)t
s(t)e
Z
s(t)e jt dt
dt
por lo tanto,
Z
|S( + h) S()|
n|to
=0
282
jt
j0 t
=
s(t)[e
e
]dt
as,
|s(t)[e jt e j0 t ]|dt
|S() S(0 )|
Z
Z
=
s()
e jt r(t )dt d
Teorema de Fubini
Z
Z
=
s()e j
e j e jt r(t )dt d
Z
Z
j
=
s()e
d
e ju r(u)du
con u = t
= [F(s)]()[F(r)]()
283
A R1
r(t)ejt dt,
r(t)est dt,
sBC
S()e jt d
Z r
1
r
|| jt
e S()d = s(t)
r
donde
Kr (t) = r
sen( 2t )
t
2
!2
lm
1
[ f Kr ](t) = f (t)
2
284
h sen(r/2) i
j
1 ||
d = r r/2 . Sugerencia: Un
r e
tringulo es la convolucin de un pulso con sigo mismo.
donde Kr = r
Rr
r
||
r
e jt S()d = s Kr ,
h sen(t/2) i2
t/2
1
2 [s
Kr ](t) = s(t).
285
ces.
vi) (Transformada del complejo conjugado.) Si la transformada de s es S, la
transformada del complejo conjugado s de s est dada por,
[F(s )] = S ()
R
Para ver esto, dado que para todo nmero complejo z, (z ) , que ( z) =
R
z , y que (e j ) = e j , se tiene que,
Z
Z
jt
s (t)e
dt =
[s(t)e jt ] dt
#
"Z
jt
s(t)e dt
=
= [S()]
vii) (Modulacin.) Si la transformada de s(t) es S(), entonces la transformada
de e j0 t s(t) es S( 0 ):
Z
F[e j0 t s(t)]() =
e j0 t s(t)e jt dt
Z
=
s(t)e j(0 )t dt
= S( 0 )
Esta propiedad es base para el anlisis de sistemas de modulacin vestigial (para TV). DSB y SSB (radio no comercial), AM y otros.
286
1
1
+
1 + (10 )2 1 + (10 + )2
S() =
Figura 4.9: Sistema bsico de modulacin y demodulacin DSB. (No se incluyen filtros ni
amplificadores.)
Z
=
s(x)e j(x+) dx
=e
S()
[u(t 2/5) u(t 8/5)] sen 5t, donde u(t) es la funcin escaln.
287
Z
Z
=
f (t))e jt dt +
g(t))e jt dt
idad y desplazamiento
S() = [F(s)]()
Z
r(t + h) r(t) jt
=
e
dt
h
Z
Z
r(t + h) jt
r(t) jt
=
e
dt
e
dt
h
h
1 jh
[e
1]R()
=
h
Ejercicio. Calcular el lmite cuando h tiende a cero de 1h [e jh 1] (R/ j).
Ejercicio. Sea g(t) = et . Grafique g0 (t) y g00 (t). Encuentre y grafique las
288
tomando t = 0 resulta,
Z
Z
1
s()r()d =
S()R()d
2
ahora, asuma que s(t) = r (t), entonces s(t) tiene transformada S() =
R () o, equivalentemente, s (t) = r(t) y R() = S (), y se tiene que,
Z
Z
1
S()S ()d
s()s ()d =
2
La frmula deseada.
La frmula de Parseval nos dice que la energa de una seal es proporcional a la energa de su transformada. Esto nos permite llamar a
|S()|2 el espectro de energa de la seal s(t). Por esta razn tambin,
289
1
s()r ()d =
2
S()R ()d
es integrable en magnitud.
Ejercicio. Encuentre una relacin entre la correlacin de dos seales
S()R (
290
Z
Z
s(t) cos(t)dt
=
s(t) sen(t)dt
Re[S()] =
s(t) cos(t)dt
Z
=
s(u)eju du
Z
=
s(u)eju
= S()
291
4.3.
r(t) =
e j0 (t) dg()
Z
e j0 t e j0 dg()
Z
j0 t
=e
e j0 dg()
292
R
dado que e j0 dg() es un nmero complejo, que depende de 0 , vemos
que, a menos que tal nmero sea cero, r(t) es una exponencial compleja de
frecuencia 0 . A este factor de amplificacin para exponenciales complejas,
en funcin de la frecuencia, lo llamamos la funcin de transferencia del
sistema, que denotamos H():
Z
H() =
e j dg()
Z Z
=
e jt s(t )dg() dt
Z
Z
j
=
e
e ju s(u)du dg()
Z
Z
j
=
e dg()
e ju s(u)du
= H()S()
donde H es la funcin de transferencia del sistema, y S la transformada de
Fourier de la entrada s(t).
293
294
,0
= lm g(t) lm g(t)
t
= lm g(t)
t
295
Las transformadas involucradas en los dos ejercicios siguientes tienen sentido en la teora de la transformada de Fourier L2 , que cubriremos ms adelante.
Las incluimos debido a su popularidad en el rea de tratamiento de seales.
296
1,
(, )
H() =
1/2, {, }
0,
< [, ]
Ejercicio. Sea
t=0
1,
e(t) =
senc t , t , 0
t
Suponga que esta seal se aplica a un filtro con la funcin de transferencia H()
dada en el ejercicio anterior. Grafique aproximadamente el espectro |Y()|2 de
la seal de salida.
Ejercicio. Encuentre la funcin de transferencia de un pasa bajas RC de primer
fundamental de una onda cuadrada. Suponga que se tiene una onda cuadrada
297
4.4.
298
por
H() = H(0 + j)
G() = G(0 + j)
Es decir, la transformada de Fourier corresponde a la transformada de Laplace
evaluada sobre el eje imaginario del plano complejo s.
Asumiendo que la funcin de transferencia H(s) es una funcin racional de
la variable s (lo cual ocurre siempre que el sistema tenga una ecuacin diferencial correspondiente), por el teorema fundamental del lgebra, tenemos que
tanto el numerador como el denominador se pueden factorizar. As, podemos
escribir
a(s z1 )(s z2 ) (s zn )
H(s) =
(s p1 )(s p2 ) . . . (s pm )
para una funcin de transferencia con n ceros y m polos. Al restringirse al
dominio de la frecuencia de Fourier, tenemos:
H() =
a( j z1 )( j z2 ) ( j zn )
( j p1 )(j p2 ) . . . ( j pm )
|a||(j z1 )||( j z2 )| |( j zn )|
|( j p1 )||( j p2 )| . . . |( j pm )|
299
1
s2 +s+1
H(s) =
s
s2 +4s+1
H(s) =
s2 +1
s2 +4s+s
H(s) =
s2
s2 +4s+1
300
4.5.
(t) =
cn e j0 nt
n=
301
h:
(t) = (t) h(t)
Z
=
h()(t )d
Z
h()
=
=
=
=
Z
X
n=
X
n=
X
n=
cn e j0 n(t) d
n=
h()cn e j0 n(t) d
Z
cn e
j0 nt
h()cn e j0 n d
cn e j0 nt H(n0 )
H(n0 )cn e j0 nt
n=
por lo tanto, la salida tiene representacin como serie de Fourier dada por:
(t) =
dn e j0 nt
n=
302
enteros de una cierta frecuencia fundamental), la suma es una seal casiperidica, y se tiene el mismo resultado: los coeficientes de las exponenciales de
la seal de salida estn dados por el producto del coeficiente de cada exponencial de entrada por la funcin de transferencia evaluada en la frecuencia de la
exponencial. en general, si la entrada est dada por,
(t) =
cn e j0 nt
n=
la salida ser,
(t) =
dn e j0 nt
n=
donde dn = cn H(n0 ).
4.5.1.
Respuesta a sinusoides
Suponga que la respuesta impulso h(t) de un sistema de convolucin continuo es real. Entonces, la magnitud de la funcin de transferencia H() del
sistema es par, y la fase es impar. Si a este sistema se le aplica la seal cos(0 t),
la respuesta r del sistema, por linealidad, estar dada por:
r(t) = 0,5|H(0 )|e jH(0 ) e j0 t + 0,5|H(0 )|e jH(0 ) e j0 t
= 0,5|H(0 )|e jH(0 )+j0 t + 0,5|H(0 )|e jH(0 ) j0 t
(Por qu?)
aplica a un filtro pasa bajas ideal con frecuencia de corte en 100Hz. cul es la
seal de salida del filtro? (Expresela como una serie de sinusoides).
303
Ejercicio. Una onda tringulo de frecuencia 1cps (ciclo por segundo) es pasada
Ejemplo. Suponga que se tiene una seal cuadrada g(t), como se muestra en
304
bn sen(nt)
n=0
donde a0 = 1/2, bn = 2/(n) con n impar y bn = 0 con n par. Suponga que g(t)
4.5.2.
305
306
1
1+ j
se
Figura 4.25: La parte real de la funcin de transferencia del pasa bajos RC con RC=1
307
Figura 4.26: La parte imaginaria de la funcin de transferencia del pasa bajos RC con RC=1
1
1+ j
308
4.5.3.
1
1+j
Comportamiento transiente
309
4.5.4.
R
Considere la respuesta cos(0 [t ])dg() a una sinusoide aplicada desde
R
menos infinito y la respuesta cos(0 [t ])u(t )dg() a una sinusoide
aplicada en cero:
Z
Z
Z
cos(0 [t])dg()
cos(0 [t])u(t)dg() =
cos(0 [t])dg()
R
t
cos(0 [t])dg() =
0.
Ejercicio. diga qu condiciones debe cumplir g para tener comportamiento
asinttico.
4.5.5.
representacin de un sistema con respuesta escaln discontinua en trminos de una suma de sistemas identidad
(linea de retardo o filtro pasa todas) y un sistema con respuesta escaln continua
310
caso, la respuesta escaln est dada por g(t) = u(t ). (Dado que este sistema
es claramente lineal e invariante, y dada la propiedad de desplazamiento en
tiempo de la transformada de Fourier, es ese sentido, diremos que la funcin
de transferencia de la linea de retardo, tambin conocida como un filtro pasa
todas, est dada por: H() = e = j.)
Suponga que se tiene un sistema de convolucin con respuesta escaln
g(t) con derivadas excepto en un conjunto finito de puntos, con una posible
cantidad finita de saltos, y acotada. Note que podemos expresar una respuesta
escaln as como la suma de: una seal derivable en sentido relajado, es decir
continua y con derivada excepto en un conjunto finito de puntos, por una
parte y una seal escalonada, es decir, una suma de seales escaln desplazadas
y escaladas, por la otra. La seal escalonada corresponder a la suma (es decir el
paralelo) de sistemas identidad y la seal sin saltos a un sistema de convolucin
con funcin caracterstica.
Ejemplo. La respuesta escaln de un RC pasa altas es g(t) = et u(t) = u(t)
(1 et )u(t). As, decimos que un pasa altas es un pasa todas menos un pasa
bajas. Vea las figuras 4.29, 4.30 y 4.31.
311
312
4.6.
4.6.1.
Filtro ideal
e j5 , cuya magnitud es constante e igual a uno y con fase lineal dada por
H() = 5. Si la entrada al filtro es e(t) = cos(at) + cos(bt), la salida est dada
por,
cos(at + H(a)) + cos(bt + H(b))
= cos(at + 5a) + cos(bt + 5b)
= e(t + 5)
Por lo tanto, el nico efecto de este filtro es adelantar la seal que se filtra.
4.6.2.
Distorsin
313
1,
HI () =
0,
HR C() =
|| 1
|| > 1
1
1 + j
314
Figura 4.32: a. Seal de entrada b. Seal de salida del filtro ideal c. Seal de salida del filtro RC
315
La distorsin de fase que sufre una seal que se filtra con un filtro de fase
no lineal, es similar a lo que le sucede a una onda que viaja a travs de un
medio dispersivo.
4.6.3.
Un filtro pasa todas ideal nos permite ver lo que sucede en la banda
de paso de un filtro que no distorsiona. Un filtro pasa todas es un filtro con
respuesta escaln de la forma: g(t) = u(t ), es decir, une escaln atrasado o
adelantado, dependiendo del signo de . Por lo tanto, la funcin de transferencia de un pasa todas est dada por:
Z
H() =
ej0 t dg(t)
: R1 S1
316
4.6.4.
Dcadas y octavas
Las dcadas y las octavas se usan para medir relaciones (es decir cocientes)
en forma logartmica. As, se dice que 20.000 est 3 dcadas por encima de 20 y
317
que 1000 est 3 octavas por debajo de 8000. Las frmulas para medir la relacin
entre f1 y f2 en dcadas y en octavas, respectivamente, estn dadas por,
D = log10 (
f1
)
f2
O = log2 (
f1
)
f2
Figura 4.33: Do=C (261.626Hz), Re=D (293.665Hz), Mi=E(329.628Hz), Fa=F (349.228Hz), Sol=G
(391.995Hz), La=A (440.0Hz), Si=B (493.883Hz)
4.6.5.
318
x1/x2
log(x1/x2)
log2 (x1/x2)
10 log(x1/x2)
20 log(x1/x2)
300/30
1
3.322
10
20
880/440
0.301
1
3.01
6.02
25/50
-3.01
-1
-3.01
-6.02
Relacin natural
en dcadas
en octavas
dB de potencia
dB de voltaje
319
1
1
]1/2 )
|) = 20 log([
1 + jRC
1 + (RC)2
= 10 log(1 + (RC)2 )
= 20 log(/c )
donde c = 1/RC.
Ejercicio. Para la funcin de transferencia H del ejercicio anterior, con c = 10,
grafique 20 log10 H() contra log10 . Si la grfica es aproximadamente lineal,
encuentre la pendiente.
Ejercicio. Si la entrada al circuito RC pasa bajas es una seal escaln, y el
320
el voltaje para que disipe una potencia de 30W? (Asuma que el cambio de
temperatura no afecta significativamente la resistencia del bombillo, lo cual es
falso.)
Ejercicio. Se tiene una resistencia de 1K, si se aplica un voltaje V=10V, Qu
potencia disipa? Sea x el voltaje que se debe aplicar a la resistencia para que
disipe 5W, Cunto vale la relacin r = x/10 Cunto vale 20 log(r)?
Ejercicio. Se tiene un pasa bajas RC con frecuencia (angular) de corte de 10 r/s.
321
pasa bajas de primer orden y un RC pasa altas de primer orden. a) Si las dos
etapas se aslan con un seguidor. b) Si la salida de una se conecta directamente
a la entrada del otro.
Ejercicio. Suponga que tiene un filtro pasa todas con funcin de transferencia
tal como se indica en la figura 4.36; si la entrada es e(t) = sen t + sen 2t, Cul
es la salida s(t)? Grafique e(t) y s(t) para t [0, 2].
Ejemplo. Si la seal cos t 0,2 cos 3t se filtra con el pasa todas de la figura 4.36,
322
2. Calcule los coeficientes de Fourier {cn } para una onda tringulo con periodo 1/60 seg., amplitud de dos voltios pico a pico, nivel DC de un voltio,
como se muestra en la figura siguiente.
323
324
4.7.
4.7.1.
325
Definicin
Sea s(t) una seal con magnitud al cuadrado integrable, es decir, una seal
en L2 . Sea N un nmero natural y considere las versiones truncadas de s dadas
por: s(t)[u(t + N) u(t N)]. Note que cada sn es integrable en magnitud, y en
magnitud al cuadrado, es decir est en L1 L2 . Por lo tanto, cada sN (t) tiene una
transformada de Fourier L 1 que denotaremos SN (). Se puede mostrar que
{SN } es una sucesin de Cauchy en L2 y, dado que L2 es un espacio completo,
existe una funcin S() tal que kS Sn k2 tiene a cero a medida que N crece. A
sa S la llamamos tambin el lmite en norma L 2 de {SN }. Este lmite nos da la
transformada de Fourier L 2 de la seal s.
R N
Ejercicio. Calcule el lmite, a medida que N crece, de N senc (t)e jt dt.
Ejercicio. Encuentre la transformada de Fourier L 2 de la seal senc .
4.8.
326
R}, determina una familia contable de seminormas para L con la que se define
su topologa: { fn } converge a 0 si y slo si el lmite cuando n tiende a infinito
de p, ( f ) es cero, para cada y cada naturales.
El espacio de funcionales continuos sobre L, que se denota L0 se llama
espacio de distribuciones temperadas.
Teorema 8. La transformada de Fourier es una aplicacin continua de L en L tal que
1.
g(x)H(x)dx =
g.
2. g(t) =
R
1
2
(4.7)
4.8.1.
Distribuciones
( f ) =
g(t) f (t)dt
( f ) =
f (t)dg(t)
327
La distribucin de Dirac
considere la distribucin dada por ( f ) = f (0), que no tiene ncleo pero
s semilla: el escaln u(t). Por definicin, su transformada es la distribucin
dada por
( f ) = (F)
= F(0)
donde F es la transformada de f . As, s tiene ncleo, ya que
( f ) = F(0)
Z
f (t)e j0t dt
=
Z
=
f (t)1dt
Z Z
=
f (t)e jt dtd
Z
Z
?
=
f (t)
e jt ddt
328
Z
f (t)dt
=
0
Z
Z
?
=
f (t)
e jt ddt
R
0
329
e jt
d.
||>0
e jt
d
||>0
4.9.
330
331
4.9.1.
Pares de Hilbert
z(t) tiene Z() = 0 para < 0 entonces Re(z) y Im(z) son pares de Hilbert.
Demostracin
Suponga f (t) := Rez(t) y g(t) := =z(t) pare de Hilbert. Entonces:
G() =
1
sgn()F()
j
Entonces
Z() = F() + jG()
= F() + sgn()F()
Entonces
<0
0,
Z() =
2F(), > 0
332
Suponga que no existe Z() = 0 para cada > 0 desde que F y G sean reales,
Z() = F() + jG() = F () + jG () = [F() + jG()]
entonces para > 0, F() = jG(), y por hiptesis tambin para < 0,
F() = jG(), tenemos:
F() = jsgn()G()
G() =
1
sgn()F()
j
4.10.
Apndice I
4.10. APNDICE I
333
1
K1
=
1/A + K + j/A0
1 + j/AK0
334
4.11.
4.11.1.
c(t) =
1
2X1
n
2nt
(E + F) + (E F)
sen
cos
2
n
2
T
(4.8)
n=1
2X
n
1
RC
(E + F) + (E F)
sen
cos(n0 t + atan (n0 RC))
p
2
2
2
n=1 n 1 + (n0 RC)
(4.9)
donde 0 = 2/T. Note que las seales tienen armnicas pares nulas. La
seal r(t) es tambin la respuesta en estado estable del circuito; es decir, es la
respuesta a la entrada e(t) = c(t)u(t) para t >> 5RC.
r(t) =
336
4.11.2.
Respuesta escaln
Suponga que en el circuito de la figura 4.44, el voltaje inicial del condensador, vc (0 ) es B y que a la entrada se aplica la seal i(t) = B + (A B)u(t); as,
para t < 0, vi (t) = B y para t > 0, vi (t) = A. Resolviendo la ecuacin diferencial
correspondiente con la condicin inicial dada encontramos que,
t
vc (t) = A[a e RC ] + Be RC
4.11.3.
(4.10)
4.12. PROBLEMAS
337
D = F(1 e 2 ) + Ce 2
de las que podemos deducir
T
C = (E + Fe 2 )
1 e 2
1e
= (E + Fe 2 )
1 e 2
D = (F + Ee
T
T
2
1 e 2
1e
= (F + Ee 2 )
1 e 2
C D = (E F)
1 e 2
T
1 + e 2
T
= (E F)tanh ( )
4
1e 2
T
1+e 2
T
= tanh ( 4
).
T
T
En la figura 4.48 se grafica tanh ( 4
) contra 2
. En las figuras 4.49 y 4.50 se
muestra la salida en estado estable para dos valores extremos de T, en un caso
para T << y en el otro para T >> .
4.12.
Problemas
338
Figura 4.48: Grfica de la relacin de los voltajes pico a pico de entrada y salida del circuito RC
4.12. PROBLEMAS
339
a) e|t|
b)
1
1+t2
c) cos(10t)e|t|
2
d) et
e) u(t) u(t 1)
2
f ) e5t cos(20t)
2. Diga si la seal senc (t) es integrable en magnitud, o no. Si s, d el valor
de la integral, si no, demuestre.
3. La respuesta impulso de un sistema de convolucin continuo es h(t) =
u(t)et . Encuentre la respuesta r(t) cuando la entrada al sistema es e(t)
a) cos(t) + cos(2t)
b) u(t) u(t 1). (Encuentre la transformada de Fourier de la salida.)
4. Se tiene un filtro pasa bajas ideal con frecuencia de corte en 70Hz. Encuentre la salida si la entrada es la seal mostrada en la figura 4.51.
340
6. Considere el circuito mostrado en la figura 4.52, el cual modela la impedancia de una punto de oscilodcopio X10 (atenuado por 10) y la impedancia
de entrada del amplificador del osciloscopio. Se dice que la punto est
compensada si C=9pF.
a) Encuentre la respuesta escaln para C=9pF, 5pF y 12 pF. (Ha visto
estas formas de onda en la pantalla de un osciloscopio?)
b) Encuentre la funcin de transferencia del circuito para C=9pF, 5pF
y 12pF.
4.12. PROBLEMAS
341
te 2 .
10. Se tiene un sistema de convolucin cuya respuesta impulso est dada
por
1/2, |t| 1
h(t) =
0,
|t| > 1
encuentre la respuesta del sistema a la entrada e(t) = sen(t) + sen(10t).
11. Se tiene un filtro con funcin de transferencia H() =
1
1+ j .
Si la entrada
R
est dad por e(t) = sen(t) + sen(2t), sea s(t) la salida. Encuentre e2 (t)dt
R
y s2 (t)dt. Sugerencia: utilice la frmula de Parseval para series de
Fourier (captulo 3).
j
1+ j .
342
15.
4.12. PROBLEMAS
343
Figura 4.55: Pasabajs RC con frecuencia de corte 1Hz y entrada con nivel DC
344
j
1+ j .
4.12. PROBLEMAS
345
22. Sea el funcional delta de Dirac, dado por ( f ) = f (0), encuentre g(t) para
que el funcional ( f ) dado por ( f ) = (F), donde F es la transformada
de Fourier L1 de f , se puede expresar como:
Z
( f ) =
f (t)g(t)dt
346
muestre que el funcional ( f ) dado por ( f ) = (F), donde F es la transformada de Fourier L1 de f , se puede expresar como:
Z
( f ) =
f (t)G(t)dt
4.12. PROBLEMAS
347
r(t)
.
e j0 t
348
4.12. PROBLEMAS
349
e jt
d
||>
350
38. Encuentre una funcin de fase no nula, que sea par e impar. (Sugerencia:
= .)
R
39. Muestre que, si g es integrable en magnitud, entonces lm0 |g(t +
) g(t)|dt = 0.
40. El circuito mostrado est en estado estable sinusoidal. Encuentre la frecuencia de oscilacin y grafique las seales a(t), b(t) y c(t).
41. A un RC pasa bajas de primer orden con frecuencia de corte de 3dB de
1, se aplica la seal e(t) = arctan (t). Encuentre la energa de la seal de
salida.
42. Se tiene un sistema de convolucin continuo con funcin caracterstica
dada por cos(10t)e|t| . Diga qu tipo de filtro es (pasa bajas, pasa altas,
etc.) y d las frecuencias de corte de 3dB.
43. Se tiene la seal x(t) con transformada X() = [u(+11)u(+9)+u(
9) u( 11)]e j . Grafique la transformada de la seal x(t) cos(10t).
44. Encuentre la respuesta del siguiente sistema a la seal 1 + cos(t) + cos(2t).
Exprsela como una suma de cosenos desfasados.
45. Con un filtro con funcin de transferencia H() = 1+1j se filtra la seal
| cos(5t)|, encuentre la potencia promedio de las primeras 3 armnicas de
la seal de salida.
46. Utilizando las propiedades de la transformada de Fourier, sabiendo que
la transformada de f (t) = u(t + 1) u(t 1) es F() = 2senc (), encuentre
la transformada de g (t), si g(t) = e j10t [u(t) 2u(t 1) + u(t 2)].
47. Utilizando las propiedades de la transformada de Fourier
, sabiendo que
( )2
(at)2
la transformada de Fourier de f (t) = e
es F() = a e 2a , encuentre
el ancho de banda de un filtro con respuesta impulso h(t) = f 0 (t).
4.12. PROBLEMAS
351
352
4.12. PROBLEMAS
353
354
|s(t)|2 dt =
|S()|2 d
Referencias
[1] Digital Signal Processing. A. V. Oppenheim y R.W. Shafer.
Prentice Hall, Englewood Cliffs, 1975.
[2] Signals and Systems. A. V. Oppenheim, A. S. Willsky y I. T.
Young. Prentice Hall, Londres, 1983.
[3] Mathematical Analysis. T.M. Apostol. 2nd Ed. AddisonWesley, Reading, 1974.
[4] Fourier Transforms. S. Bochner y K. Chandrasekhanran.
Princeton University Press, Princeton, 1949.
[5] Anlisis de Fourier. J. Duoandikoetxea. Addison- Wesley,
Wilmington, 1995.
[6] Sistemas Digitales y Analgicos, Transformadas de Fourier, Estimacin Espectral. A. Papoulis. Marcombo, Barcelona,
1986.
[7] Special Functions and their Applications. N.N. Lebedev.
Dover, N.Y. 1972.
[8] Exercices Rsolus de Mathmatiques. M. Carbon, P. Marry,
N. Point y D. Vial. Dunod, Paris, 1986.
355
356
REFERENCIAS
[9] Fourier and Wavelet Analysis. G. Bachman, L.Narici, E. Beckenstein. Springer, N.Y., 2000.
Captulo 5
Introduccin
357
358
sn e jn
n=
X
X
jn
jn
sn e
sn e
n=
n=
X
=
|sn ||e jn |
=
n=
|sn |
n=
= ksk1
Ejercicio. Muestre que cada seal discreta s : Z C en l2 , est en l1 .
Note que la transformada S() de una seal discreta es una funcin peridica, para la que 2 es un tiempo de repeticin: S( + 2) = S(), ya que 2
es un tiempo de repeticin de cada trmino sn e jn .
Tambin note que S() tiene serie de Fourier con coeficientes {sn }. Esta
propiedad nos permitir dar una frmula para la transformada inversa DTFT1
basada en los resultados sobre series de Fourier, vistos en el captulo 3. Antes,
veamos un ejemplo.
Ejemplo. Sea s : Z C la seal dada por
1, n {1, 0, 1}
sn =
0, n < {1, 0, 1}
5.0. INTRODUCCIN
359
entonces
S() =
sn e jn
n=
1
X
sn e jn
n=1
= e j + 1 + e j
= 1 + 2 cos
Tanto S como su magnitud y fase, se grafican a continuacin.
360
361
y, cada exponencial compleja sk e jk tiene periodo 2/k. en general, el periodo de S es el mnimo comn mltiplo de los periodos de las exponenciales
complejas no nulas en la serie, ya que stas son funciones linealmente independientes.
5.0.1.
Transformada inversa
X
1
|S()|2 d
|sn |2 =
2
0
n=
que relacione la norma l2 de la seal s con la potencia promedio de la seal S,
en forma similar al caso correspondiente a las series de Fourier.
5.1.
362
nZ
nZ
sm e j(m+k)
mZ
= e j(k)
sm e j(m)
mZ
= e j(k) S()
Ejercicio. Como se puede mostrar, la transformada S() de la seal s(t)
363
S( 0 )e jn d
= 0
Z 20
1
S()e j(+0 )n d
=
2 0
Z
e j0 n 20
=
S()e jn d
2 0
=
e j0 n
=
2
S()e jn d
0
= e j0 n sn
R 2
1
dado que 2
S()e jn d y que la integral es independiente de dnde
0
empiece y dnde termine, siempre y cuando se integre sobre un periodo.
Por lo tanto, un desplazamiento en frecuencia del espectro de la seal,
corresponde a la multiplicacin de la seal por una exponencial compleja
discreta en tiempo, es decir, a una modulacin de la seal.
364
f (t)dt =
0
T+
f (t)dt
5.2.
365
366
2
7 n
367
5.3.
Filtros discretos
368
5.4.
Y() =
yn e jn
n=
X
X
Y() =
xi hni ejn
n= i=
X
xi hni e jn
n= i=
Y() =
xi hni ejn
i= n=
X
n=
jn
xi hni e
X
k=
xi hk e j(k+i)
369
370
entonces
Y() =
=
xi hk e j(k+i)
i= k=
X
X
j(i)
xi e
hk e j(k)
i=
k=
= X()H()
por lo tanto, se tienen dos representaciones del sistema: una en tiempo y la
otra en frecuencia, como se muestra en la figura 5.15.
Ejemplo. Se tiene un promedio mvil con ancho de ventana igual a tres, donde
la salida yi est dada por el promedio 13 (xi1 + xi + xi+1 ) de los tres componentes
de la seal de entrada en una ventana centrada en i.
Ejercicio. Demuestre que el sistema promedio mvil es un sistema de convolu-
dada por
H() =
1 j
1
(e
+ e j0 + e j ) = (1 + 2 cos )
3
3
371
Figura 5.16: La respuesta impulso de un promedio mvil con ancho de ventana igual a tres
372
1/2, n {1, 1}
hn =
1,
n=0
0,
n Z {1, 0, 1}
por lo tanto, su funcin de transferencia est dada por,
H() = 1 + cos()
Ejercicio. Grafique la magnitud y la fase e a + cos()
Ejercicio. Suponga que el filtro del ejemplo anterior se coloca en cascada con
373
Y() =
5.5.
1
1
(1 2 cos )(1 + cos ) = (cos + cos(2))
3
3
a0 + a1 z + + an zn
bo + b1 z + + bn z m
por el teorema fundamental del lgebra, tanto el numerador como el denominador pueden ser factorizados.
Qn
(z ci )
H(z) = Qmi=1
i=1 (z pi )
a las constantes ci se les llama los ceros de H(z) y las pi los polos de H(z).
Ahora, dado que H() = H(ei ) tenmos que
Qn
|(ei ci )|
|H(z)| = Qmi=1 i
i=1 |(e pi )|
que interpretamos geomtricamente para cada como el producto de las
distancias de e j a los ceros, partido por el producto de las distancias a los
polos. Note que dado el diagrama de polos y ceros, partido por el producto de
la figura 5.19, podemos graficar aproximadamente |H()| (figura 5.20).
374
375
5.6.
5.6.1.
N1
X
Xk e j N nk
k=0
hk e j0 (nk)
k=
= e j0 n
hk e j0 k
k=
=e
j0 n
H(0 )
N1
X
k=0
Xk e j N nk
376
N1
X
Xk H(
k=0
2 j 2 k
k)e N
N
5.7.
A continuacin consideramos la representacin en el dominio de la frecuencia de Fourier de las seales discretas de longitud finita. Las seales discretas de
longitud finita son seales cuyo dominio es un conjunto finito , normalmente
un intervalo entero, y su rango es el conjunto de los nmeros complejos.
nosotros asumiremos normalmente que el dominio de las seales discretas de
longitud finita es un intervalo de la forma /0, N 1/.
Las seales digitales discretas de longitud finita son las seales que mejor
se pueden representar en un computador digital: son arreglos (arrays) unidimensionales; su representacin en el dominio de la frecuencia resulta ser
tambin una seal de longitud finita y por lo tanto, tambin un array. Por esta
y otras razones, la transformada correspondiente, conocida como la Transformada Discreta de Fourier (Discrete Fourier Transform, DFT) es una se las
377
N1
1 X j 2 nk
xn e N
N n=0
N1
1 X
nk
x n WN
N n=0
Con base en los resultados del captulo 3, decimos que la DFT es una biyeccin lineal (un homomorfismo) de CN CN . en particular, la DFT es invertible:
a partir de X podemos recuperar x como DFT1 (X); sta transformada inversa
est dada por: x = [x0 , x1 , . . . , xN1 ] = DFT1 (X).
Ejercicio. Para la seal mostrada en la figura 5.21, calcule su transformada
discreta de Fourier.
Ejemplo. Sea x la seal de longitud 3 dada por x = [1, 2, 3]. su transformada
378
1
(1 + 2e j2k/3 + 3e j4k/3 )
3
por lo tanto,
X0 =
1
(1 + 2e j20/3 + 3ej40/3 ) = 3
3
X1 =
1
(1 + 2e j2/3 + 3e j4/3 )
3
X2 =
1
(1 + 2e j4/3 + 3e j8/3 )
3
X = [X0 , X1 , X2 ]
5.7.1.
379
Yk =
=
N1
2
1 X
x[n+m]N e j N nk
N n=0
2
2
2
2
1
(xm e j N 0k + xm+1 e j N k + + xN1 e j N (N1m)k + x0 e j N (Nm)k
N
2
+ + xm1 e j N (N1)k )
2
2
2
2
2
2
1
(xm e j N mk e j N mk + xm+1 e j N mk e j N (m+1)k + + xN1 e j N (N1)k e j N mk
N
2
2
2
+ x0 e j N mk + + xm1 e j N mk e j N (N1+m)k )
2
2
2
1 j 2 mk
e N (xm e j N mk + xm+1 e j N (m+1)k + + xN1 ej N (N1)k + x0
N
2
+ + xm1 e j N (N1+m)k )
2
2
2
2
1
(x0 + + xm1 e j N (N1+m)k + xm e j N mk + xm+1 e j N (m+1)k
N
2
+ + xN1 e j N (N1)k )
= e j N mk
= e j N mk Xk
380
5.7.2.
5.7.3.
N1
N1
X
1 X
mk 1
nk
r m WN
sn WN
N m=0
N n=0
N1 N1
1 XX
mk
nk
rm WN
s n WN
N2 m=0 n=0
N1 N1
1 XX
(m+n)k
rm sn WN
N2 m=0 n=0
y, dado que,
ui = [DFT (RS)]i =
1
N1
X
k=0
ik
Rk Sk WN
381
se tiene que,
ui =
N1
X
ik
Rk Sk WN
k=0
N1 N1
N1
X
1 XX
(m+n)k ik
=
rm sn WN
WN
N2 m=0 n=0
k=0
N1
N1
N1
1 X X X (imn)k
rm
sn
WN
= 2
N m=0 n=0
k=0
(imn)k
=N
(imn)k
=0
WN
k=0
WN
k=0
por lo tanto
N1
1 X
ui =
rm sJimKN
N m=0
382
Figura 5.22: Los trminos rm sim a ser sumados, para computar ui . Note el eje m
5.7.4.
N1
X
nk
,
Rk WN
nZ
nk
,
Sk WN
nZ
k=0
sn =
N1
X
k=0
383
[3, 3, 3, 3, 3, 3, 2, 1, 0].
Figura 5.23: Las seales rn y sn , y las correspondientes extensiones peridicas rk y s2k , necesarias
para el clculo de u2
384
5.7.5.
n=0
k=0
N1
N1
X
1 X
xk r[k]N =
Xk Rk
N
k=0
rk = x[k]N , Rk = Xk
k=0
N1
N1
X
1 X
xk x[k]N =
Xk Xk
N
k=0
k=0
N1
N1
1 X 2 X
|xk | =
|Xk |2
N
k=0
k=0
5.8.
385
5.8.1.
simtrica, simultneamente.
386
G.
G () = [F(g)] () = [F(g )]()
Similarmente, si s es una seal discreta con transformada de Fourier S.
S () = [F(s)] () = [F(s )]()
387
5.8.2.
a b
2
388
Figura 5.25: La restriccin de una seal discreta a [0, ]. a. es suficiente para reconstruir
la transformada en [, ]. b. ya que sta es par y, para cada real. c. ya que es peridica con
periodo 2
389
R1 .
P
jn
Im (sn ) = 0 y S() =
implican,
n= sn e
Re [S()] =
sn cos(n)
funcin par
n=
Im [S()] =
sn sen(n)
funcin impar
n=
p
Im [S()]2 + Re [S()]2
S() =
2n e jn +
n=0
2n e jn 1
n=0
+
1
j
j
1 e2
1 e2
2 cos()
=
1
1,25 cos()
Ejercicio. Grafique la magnitud y la fase de la transformada de 2|n| .
Nota. Hemos visto 4 transformadas de Fourier: series de Fourier, trasformada
390
5.9.
sn e jn =
n=
N1
X
sn e jn
n=0
N1
1 X j 2 nk
sn e N
N n=0
1 2
S( k)
N N
391
5.10.
Una seal continua se puede procesar con un sistema continuo. Sin embargo, la versatilidad y resistencia al ruido de los sistemas digitales hacen
prctica otra alternativa para seales que no incluyen frecuencias por encima
de aproximadamente 300kHz ; por ejemplo, seales ssmicas, de audio, potenciales nerviosos, etc. As, la seal continua se convierte en una seal digital, se
procesa digitalmente y se reconvierte en seal anloga con un conversor D/A.
El primer paso consiste en muestrear y cuantificar (en nmero de bits dados, e.g.
16) la seal continua; esto se logra con un S/H (Sampler and Holder), y un
conversor A/D (anlogo/digital) y un filtro pasa bajas que suaviza la salida del
D/A. Esta es la tecnologa usada para el tratamiento digital de audio y la que
seguramente ser usada en TV comercial en el futuro.
392
Para que una seal continua pueda ser filtrada pueda ser filtrada efectivamente as, es necesario que no se pierda informacin al pasar de la seal
continua a la seal discreta, es decir, que la seal continua sea recuperable a
partir de sus muestras. Una condicin que garantiza este requisito, es que la
seal continua tenga transformada con soporte acotado y que las muestras se
tomen a una tasa lo suficientemente alta. Esto se explica a continuacin.
Es importante advertir que los efectos del proceso de cuantificacin no se
consideran aqu. La cuantificacin tiene dos efectos: introducir un ruido con
densidad uniforme en la seal que se muestrea y hacer los sistemas discretos
efectivamente no lineales. La no linealidad de los filtros puede hacer que se
produzcan oscilaciones inesperadas (limit cycles, por ejemplo) o incluso
seales caticas. Una referencia clsica para este tema es [1].
5.10.1.
Suponga que tiene una seal continua g(t) la cual se muestrea a una tasa
de muestras por unidad del parmetro t, obteniendo la seal discreta x dada
por xm = g(nT), donde T es un nmero real positivo y n Z. Asuma que tanto
x como g y G son seales en ele uno. Sean G() y X() las transformadas,
respectivamente, de las seales g(t) y xn . Por lo tanto,
1
T
X() =
xn e jn
(5.1)
X()e jn d
(5.2)
n=
1
xn =
2
Z
G() =
g(t)e jt dt
(5.3)
G()e jt d
(5.4)
g(t) =
393
1
2
G()e jnT d
Z (2r+1)/T
1 X
=
G()e jnT d
2 r= (2r1)/T
1
2r jnT
xn =
)e
d
G(/T +
2 /T r=
T
con el cambio de variable = T,
Z X
1
2r
xn =
G( +
) e jn d
2 T r= T
T
comparando con 5.2 tenemos que,
X() =
1 X
2r
G( +
)
T r= T
T
394
395
5.10.2.
1
2
/T
G()e jt d
/T
396
xn e jn , tenemos que,
Z /T X
1
T
g(t) =
xn e jTn e jt d
2 /T n=
Z /T X
T
g(nT)e jTn e jt d
=
2 /T n=
n=
X
T
j(tTn)
g(t) =
g(nT)
e
d
2 /T
n=
=
X
n=
g(nT)
sen[ T (t nT)]
T (t
nT)
g(nT)senc ( (t nT))
T
n=
(5.5)
esta ltima frmula expresa la seal g(t) en trminos de sus muestras {g(nT)}.
Desde cierto punto de vista esta es una formula de interpolacin; la interpolacin
se hace usando funciones senc .
Este resultado nos dice que no hay prdida de informacin al pasar de una
seal limitada en frecuencia a la sucesin de sus muestras, si stas se toman
con una rapidez de al menos la tasa de Nyquist.
Esta frmula es vlida tambin para seales g L2 . El lector interesado
puede consultar el artculo siguiente:
A.J. Jerry. The Shannon Sampling Theorem, its Various Extensions . . . Proc.
IEEE vol. 65, pp. 1565-1596, Nov. 77.
La frmula 5.5 se puede interpretar as: sobre la lnea de los nmeros
reales, marque las posiciones t = nT, coloque una seal sen( T (t nT) centrada
en t = nT, para cada marca; luego, sume todas las seales. Note que cada
senc se anula en las posiciones t = NT excepto en la que est centrada.
397
la figura 5.30 se muestran las seales g0 senc (t), g1 senc (t 1) y g2 senc (t 2). En
la figura 5.31 se muestran la seal f (t) = g0 senc (t)+ g1 senc (t1)+ g2 senc (t2).
Note que efectivamente, f (0) = g0 , f (1) = g1 y f (2) = g2 .
398
5.10.3.
Suponga que nos dan N muestras {s0 , s1 , . . . , sN1 } de la seal continua g(t),
g : [0, T] C, espaciadas regularmente en el intervalo [0, T]:
sn = g(nT/[N 1]), N {0, 1, . . . , N 1}
Es posible entonces encontrar una seal continua que es una suma de
exponenciales complejas y que concuerda con las muestras dadas de g. Esta
seal est dada por,
N1
X
2
g =
Xk e j A kt
k=0
donde
Xk =
N1
1 X j 2 nk
sn e N
N
k=0
indicadas.
399
Ejercicio. Interpole las muestras 0, 1, 0,8, 0,6, 0 usando la frmula dada. Grafique
400
5.11.
Resumen grfico
401
402
5.12.
Serie de Fourier:
cn =
1
T
f (t)e jn0 t dt
f (t) =
cn e jn0 t
t R1 periodo T = 2/0
n=
Suma de Fourier:
Xk =
N1
2
1X
f (t)e jn N nk dt
T n=0
xn =
N1
X
Xk e jn N nk
n Z, periodo N
k=0
sn e jnn dt
sn =
n=
1
2
S()e jnn d
0
xn =
N1
X
k=0
Xk e jn N nk
n N, k N N = /0, N 1/
403
t R1 , R1
404
5.13.
|cn |2 =
n=
1
T
| f (t)|2 dt
0
Suma de Fourier:
N1
N1
X
1 X
|xn |2 =
|Xk |2
N n=0
k=0
|sn |2 =
n=
1
2
|S()|2 d
0
N1
N1
X
1 X
|xn |2 =
|Xk |2
N n=0
n=0
1
| f (t)| dt =
2
|F()|2 d
405
406
5.14.
Resumen convolucin
cn dn e
k=
jn0 t
1
=
T
f (t )g()d
0
Suma de Fourier:
N1
X
nk
Rk Sk WN
=
k=0
N1
1 X
sm rJnmK
N m=0
S()R()d =
0
sk rnk
k=
Z
Z
1
jt
F(W)G(t W)dW
f (t)g(t)e d =
2
407
408
5.15.
Problemas
1, si n {0, 1, . . . , L 1}
rn =
0, si n < {0, 1, . . . , L 1}
a) Grafique rn contra n.
b) Muestre que la transformada de Fourier R() de r, est dada por
j L1
2
R() = e
sen( L )
2
sen
2
5.15. PROBLEMAS
409
0,
n<0
hn =
n
3
n0
Si la entrada al sistema es la seal xn = cos(n/2) + cos(n), exprese la salida
yn como una suma de cosenos (desfasados).
8. Se tiene un sistema de convolucin discreto con funcin de transferencia
H() = cos(). Encuentre la respuesta impulso del sistema.
9.
410
5.15. PROBLEMAS
411
yn = xn + 0,5yn1 .
La relacin se puede iterar hacia atrs, y es vlida en el lmite,
cuando n :
yn = xn + 0,5xn1 + 0,25yn2
= xn + 0,5xn1 + 0,25xn2 + 0,125yn3
=
Siempre que la entrada tienda a 0 cuando n , la salida tambin.
a) Diga si el sistema es de convolucin. Explique.
b) Exprese la salida como una suma de exponenciales complejas, si la
entrada est dada po JnK3 .
16. La seal g(t) tiene transformada G(), con G() = 0 para || > 210000.
Muestreando g, se obtiene la seal discreta x : Z C dada por xn =
g(nT), T = 105 . La seal x es filtrada con el sistema discreto con ecuacin
de diferencia yn = xn xn1 , produciendo la seal discreta y = {yn }.
Finalmente, a partir de la seal y se obtiene la seal continua s(t), con la
siguiente frmula de interpolacin:
s(t) =
X
n=
yn senc [
t nT
]
T
412
Note que sta es una forma verstil de implementar filtros anlogos con
base en filtros digitales. Diga cules inconvenientes tiene.
17. Diga qu tipo de filtro es el promedio mvil ponderado que para cada
seal de entrada x produce a la salida la seal y dada por yn = 0,25(xn1 +
2xn + xn+1 ). es decir,
a) Encuentre la respuesta impulso.
b) Compruebe que es de convolucin.
c) Encuentre la funcin de transferencia. Ahora, de qu tipo es?
d) Tambin, diga i el sistema es lineal, invariante, estable y causal. En
cada caso, explique.
18. Calcule la DFT de la seal [1,1,0,0,1]. Expresela como un vector real.
19. Se tiene un sistema de convolucin discreto con respuesta impulso dada
por hn = un 2n ,
a) Encuentre la funcin de transferencia H() y muestre que es real u
par.
b) Si la entrada es la seal sn = e j(20/45)n , exprese la seal de salida
como una suma de exponenciales complejas.
c) La seal g(t) tiene transformada G() con G() = 0 para || >
210000. Se obtiene luego la seal discreta x usando xn = g(nT)
con n Z, T = 105 . Luego, la seal discreta es filtrada con un
sistema de convolucin discreto con ecuacin de diferencia dada por
yn = 0,5(xn xn1 ). Finalmente, se reconstruye una seal continua
P
s(t) =
0 yn senc (/T[t nT]).
1) Encuentre la funcin de transferencia del sistema discreto. Diga
qu tipo de filtro es.
5.15. PROBLEMAS
413
1
X =
1
1
j
1
j
1
1
1
1
1
j
1
j
P3
k=0
Yk e j 4 nk .(D el
Referencias
[1] Digital Signal Processing. A. V. Oppenheim y R.W. Shafer.
Prentice Hall, Englewood Cliffs, 1975.
[2] Computer Based Exercises for Signal Processing. C-S. Burrus,
et.al., Prentice Hall, Englewood Cliffs, 1994.
[3] Dicrete-Time Signal Processing. A. V. Oppenheim y R.W.
Shafer. Prentice Hall, Englewood Cliffs, 1989.
[4] Digital Signal Processing. J.G Proakis y D.G. Manolakis. Pentice Hall, Upper Saddle River, 1996.
[5] Chaos in Digital Filters. L.O. Chua y T. Lin. IEEE trans. on
Circuits and Systems, vol. 35, no. 6, pp. 648-658,1988.
414
Captulo 6
Clasificacin
415
416
y(t) =
e()dgt ()
417
Los sistemas no lineales se han estudiado con las llamadas series de Volterra
(en honor a Vito Volterra, matemtico italiano, quien public varios resultados relacionados a finales de siglo pasado). Sin embargo, esta teora no es
suficientemente general, en particular, no modela el filtro mediana.
As, hay dos reas grandes en las que se clasifican los sistemas no lineales:
en el caso continuo, los sistemas pueden ser lisos, y si son de forma y0 =
f (y) + x, donde x es la entrada y y es la salida, cuando f tiene Jacobiano, y
los sistemas arrugados (osea no lisos) pero lineales a trozos (o PL: picewise
linear), como en el caso del circuito de Chua. Los sistemas discretos no lineales
que se estudian normalmente son de la forma yn = f (yn1 ) + xn , donde f tiene
Jacobiano o es lineal a trozos.
la seccin 6.1 pas al captulo
cero, revisar la introduccin.
6.1.
Sistemas lineales
Sean (V, , K, +, ) y (W, , K, +, ) espacios vectoriales sobre el mismo cuerpo. La funcin T : V W es una transformacin lineal si para cada
, K y cada x, y V se tiene que T( x y) = T(x) T(y). As, una
transformacin lineal es un homomorfismo entre espacios vectoriales.
En otras palabras, T es lineal si cumple con las propiedades de homogeneidad:
para cada K y cada u V, T(u) = T(u), y de superposicin: para cada
u, v V, T(u + v) = T(u) + T(v).
Ejercicio. Demuestre que si T es una transformacin lineal entre los espacios
418
Ejemplo. Considere el espacio vectorial V de las funciones contnuas con dominio [0, 1] y rango complejo, con la operacin de suma de funciones, sobre el cuerpo de los complejos. Considere el conjunto de funcionales lineales
definidos sobre este espacio. Como se puede mostrar, para cada funcional lineal f : V C, de este tipo hay una funcin g de variacin acotada tal que f se
puede representar como una integral de Riemann-Stieltjes:
1
Z
f (h) =
h(t)dg(t)
0
"Z
T(h) =
Z
h(t)e
j2t
dt,
h(t)e
j2t
#
dt
Rb
Teorema 10. Riesz. La integral de Riemann-Stieltjes a s()d() donde es una
funcin dada de variacin acotada, define un funcional lineal sobre el espacio de las funciones contnuas. Conversamente, cada funcional lineal definido sobre el espacio de las
funciones contnuas con soporte incluido en el intervalo [a, b], puede ser representado
en esta forma.
Suponga que se tiene un sistema lineal S y que las posibles seales de
entrada e son las seales contnuas con soporte acotado por el intervalo [a, b].
419
Z
[S(e)](t) =
s()dt ()
a
Z
[S(e)](t) =
s()d()
a
Z
[S(e)](t) =
s()h(t, )d
a
Para seales de entrada con soporte lo suficientemente centrado y localizado en el intervalo [a, b], con respecto a un posible corrimiento, tenemos que si
t () = ( t) entonces el sistema es invariante con respecto al corrimiento.
6.2.
420
6.2.1.
Permutaciones
421
6.2.2.
Muestra
6.2.3.
Estadsticas
422
6.2.4.
Estadsticas de orden
6.2.5.
6.2.6.
423
Algunas estadsticas L
6.2.7.
Resulta conveniente expresar las estadsticas L en notacin matricial, como el producto, xT donde el vector = [0 , 1 , . . . , N1 ] es un vector de
coeficientes, y T indica vector transpuesto. As, para N = 5, los vectores de
coeficientes para la mediana, el mnimo, el mximo, el rango y el midrango
estn dados por,
estadstica
mediana
mnimo
mximo
rango
midrango
promedio
vector de coeficientes
[0, 0, 1, 0, 0]
[1, 0, 0, 0, 0]
[0, 0, 0, 0, 1]
[1, 0, 0, 0, 1]
[0,5, 0, 0, 0, 0,5]
[0,2, 0,2, 0,2, 0,2, 0,2]
424
6.2.8.
6.2.9.
Filtros L
geneidad.
Ejercicio. D un contraejemplo que ilustre el hecho que el filtro mediana no
425
Figura 6.1: Un filtro mediana tiene tres bloques principales: uno que toma segmentos de la
seal de entrada, otro que ordena los datos en los segmentos y otro que toma el valor ordenado
intermedio, o mediana
6.2.10.
Seales raz
Los filtros no lineales son ms difciles de analizar, debido a que hay que
hacerlo caso por caso y no se tiene un volumen de conocimiento para ninguno
en particular, similar al que se ha acumulado para los filtros no lineales. La
caracterizacin de las seales que pasan sin sufrir alteracin por un filtro, o
seales raz del filtro, constituye una caracterizacin terica importante de
un filtro no lineal. (En el caso de un filtro de convolucin, las exponenciales
complejas con frecuencias para las que la funcin de transferencia del filtro
vale 1 + j0).
Tyan, Brandt y Longbotham han caracterizado las seales raz del filtro
mediana. Como dijimos, una seal x = {xn } es una seal raz (o seal invariante,
o seal idempotente) del filtro mediana si al filtrarla obtenemos nuevamente la
misma seal.
Las seales raz del filtro mediana se pueden clasificar en dos conjuntos.
Uno contiene seales localmente montonas con grado de monotonicidad local
igual al ancho de la ventana del filtro en cuestin. Por ejemplo, las seales
426
lomo-3 son seales raz del filtro mediana con ventana de ancho 3. el otro
conjunto, contiene seales binarias y peridicas, no necesariamente localmente
montonas. Tyan anot en 1981 que las seales localmente montonas eran
invariantes, posteriormente, Brandt y Longbotham, independientemente en
1987, indicaron la relevancia de los conjuntos de seales peridicas binarias. A
continuacin damos ejemplos de seales peridicas binarias. Por simplicidad,
las seales tienen rango {0, 1} en todos los casos.
Ejemplo. La seal {. . . , 0, 1, 2, 2, 1, 0, 0, 1, 2, 2, 1, 0, 0, . . . } es localmente montona
6.3.
Invarianza y Causalidad
Los nombres de invarianza y de causalidad resultan de considerar la evolucin de las seales en la variable tiempo. Hoy en da el dominio de una seal
no necesariamente se refiere a tiempo si no, por ejemplo a desplazamiento;
sin embargo, se sigue utilizando la misma terminologa. Antes de definir las
propiedades, considere el uso de cierta notacin para seales desplazadas.
Sea f : R1 C una seal continua. Sea t un nmero real y sea g : R1 C
la seal continua tal que, t R1 g(t) = f (t ), decimos que la seal g es
una versin desplazada de la seal f y a denotamos como f . Usando esta
convencin, 0 f = f y t R1 g(t) = f (t )
427
6.3.1.
Invarianza
428
6.3.2.
Causalidad
429
430
6.4.
Estabilidad
6.4. ESTABILIDAD
431
Quizs sea conveniente anotar que aunque normalmente un sistema inestable es indeseable, no siempre es as. Por ejemplo, la maniobrabilidad de un
avin puede incrementar notablemente bajo condiciones de inestabilidad. Lo
importante es poder predecir el resultado y que el sistema produzca la salida
deseada.
Ejercicio. Diga si el promedio mvil es estable en sentido EASA, o no.
Ejercicio. D un ejemplo de un sistema no estable.
Ejercicio. Se tiene un sistema discreto T el cual, para cada seal de entrada x,
produce la seal de salida T(x), x = {xn } y T(x) = {[T(x)n ]}. Si T est definida
como se indica en cada caso abajo, diga si T es:
Lineal
Causal
Invariante
Estable
1. [T(x)n ] = gn xn Cmo dependen sus respuestas de g?
2. [T(x)n ] =
Pn
3. [T(x)n ] =
Pn+n0
k=n0
xk
k=nn0
xk
4. [T(x)n ] = xnn0
5. [T(x)n ] = ex n
6. [T(x)n ] = a xn + b
432
6.5.
6.5.1.
Caso discreto
6.5.2.
Caso continuo
6.6.
433
Sistemas de convolucin
6.6.1.
434
es cero para tiempos menores que cero; como conclusin tenemos que la respuesta escaln de un sistema de convolucin causal tiene que ser cero para
tiempos menores a cero. Por lo tanto la respuesta impulso, si existe, es cero
para tiempos menores que cero (la derivada de una constante es cero).
Para sistemas discretos, la situacin es similar y el resultado es el mismo:
suponga que S es un sistema de convolucin discreto; como S es lineal, la
respuesta a la seal cero es la seal cero, como la seal impulso discreta es cero
para instancias menores que cero y el sistema es causal, la respuesta impulso
es cero para instancias menores que cero.
Por lo tanto, si un sistema de convolucin es causal, entonces su respuesta
escaln y su respuesta impulso son cero en la parte negativa de su dominio.
Si h(t) es la respuesta impulso de un sistema continuo de convolucin causal,
entonces para t < 0, h(t) = 0. Si hn es la respuesta impulso de un sistema
discreto de convolucin causal, para n < 0, hn = 0.
Conversamente, suponga que se tiene un sistema de convolucin continuo
y que su respuesta escaln se anula en la parte negativa de su dominio. Sea e
una seal de entrada; la salida est dada entonces por,
Z
r(t) =
e(t )dg()
Z =
=
e()dg(t )
=
=t
=
e()dg(t )
y, dados los lmites de la integral, concluimos que la salida r para cada entrada
e, para cada tiempo t, solo depende de los valores que e tome con anterioridad
a t. Es decir, r slo depende de la restriccin de e al intervalo(, t) de su
dominio. Por lo tanto, el sistema es causal.
435
yn =
n
X
hnk xk
k=
6.6.2.
436
|h()x(t )|d : t R1 }
sup{
Z
= sup{
|h()||x(t )|d : t R1 }
Z
sup{
|h()|kxk d : t R1 }
Z
= sup{kxk
|h()|d : t R1 }
Z
= kxk sup{
|h()|d : t R1 }
= kxk khk1
<
En forma similar se puede mostrar que si se tiene un sistema de convolucin
discreto, kyk < kxk khk1 , por lo tanto, si la respuesta impulso {hn } es sumable
en magnitud (i.e. si est en l1 ) y si la entrada es acotada (i.e. si est en l )
entonces la salida est acotada y su norma ele-infinito es menor que el producto
de las normas ele-uno de h y ele-infinito de x.
Concluimos que, si la respuesta impulso de un sistema de convolucin
tiene norma ele-uno finita, entonces el sistema es estable en sentido EASA.
Viceversa, suponga que se tiene un sistema de convolucin continuo cuya
respuesta impulso no est en L1 , es decir, su norma ele-uno es infinita. Queremos mostrar que entonces podemos encontrar una seal de entrada acotada a
la que le corresponda una seal de salida no acotada.
Suponga primero que la integral de la respuesta impulso es infinita. En este
caso, la respuesta escaln no es acotada y por lo tanto el sistema es inestable.
437
1,
x(t) =
1,
h(t) 0
h(t) < 0
y(0) =
h()x()d
Z
=
|h()|d
= khk1
=
Conjetura 2. Un sistema de convolucin con respuesta impulso no integrable en
magnitud no es estable EASA.
Ejercicio. Diga si es cierto que un sistema de convolucin continuo es estable
438
6.7.
6.7.1.
6.8.
6.8.1.
Sistemas discretos
439
x=
X
i=
xi .i
xn =
xi [i ]n
i=
440
Figura 6.7: La descomposicin de la seal en la figura 6.6 como una suma de seales impulso
pesadas y desplazadas.
seal x.
y = S(x) = S(
=
xi .i )
i=
xi S(i )
i=
y si S es invariante,
y = S(x) =
xi .i h
i=
xi hni = x h
i=
441
1/3, n {1, 0, 1}
hn =
0
n < {1, 0, 1}
muestre que la salida r, para una entrada s, est dada por rn = (1/3)(sn1 + sn +
sn+1 ).
6.8.2.
Sistemas contnuos
N
X
k=N
como se muestra en la figura 6.9. a este tipo de seales las llamaremos seales
simples. como el sistema es lineal e invariante, la respuesta correspondiente
442
N
X
k=N
443
es una funcin simple de orden N que aproxima a e(t). El lmite de eN (t) cuando
N tiende a infinito es e(t) (Por qu?). Por lo tanto, la respuesta del sistema a
e(t) est dada por
r(t) = S lm eN (t)
N
= lm
N1
X
k=N
"
#
(k + 1)T
kT
)
e( ) u(t kT ) u(t
N
N
N
Z
r(t) =
Z
r(t) =
=T
e(dg())
e(t )dg()
Es L1 un subconjunto de ?
6.8.3.
Z
f (g) =
g()d()
a
444
6.9.
Sistemas de Volterra
6.9.1.
L
X
i=0
1 Vito
h1i xni +
L X
L
X
i=0 j=0
h2i j xni xn j +
L X
L X
L
X
Volterra fue un matemtico italiano de comienzos del siglo 20 que hizo aportes importantes a la teora del anlisis funcional
6.10. PROBLEMAS
445
6.10.
Problemas
1. Se tiene un sistema discreto causal que cumple con la ecuacin de diferencia yn = (1/2)yn1 + xn , donde x es la entrada y y la salida. Encuentre
la funcin de transferencia del sistema.
2. Se tiene un sistema continuo el cual, para cada seal de entrada e(t)
Rt
produce la seal de salida s(t) = e()d. Diga si el sistema es estable.
Si s, demuestre, si no, d un contraejemplo.
3. Se tiene un sistema discreto rango mvil con ventana de ancho tres.
La relacin entrada-salida est dada por, yn = rango{xn1 , xn , xn+1 } =
max{xn1 , xn , xn+1 } mn{xn1 , xn , xn+1 }. Diga si el sistema es estable y si
es invariante. en cada caso, si s, demuestre, si no, d un contraejemplo.
4. Se tiene un sistema continuo que produce a la salida la seal s(t) dada
R
por s(t) = e|t2| r()d cuando la entrada es la seal r(t).
a) Encuentre la salida correspondiente a la entrada u(t).
b) Encuentre la salida correspondiente a la entrada u(t 10).
c) Diga si el sistema es estable. (Demuestre o d un contraejemplo.)
d) Diga si la seal es lineal. (Demuestre o d un contraejemplo.)
446
6.10. PROBLEMAS
447
448
Referencias
[1] Digital Signal Processing. A. V. Oppenheim y R.W. Shafer.
Prentice Hall, Englewood Cliffs, 1975.
[2] Mathematical Analysis. T.M. Apostol. 2nd Ed. AddisonWesley, Reading, 1974.
[3] Fourier Series. G.P. Tolstov. Dover, N.Y., 1962.
[4] Almost Periodic Functions. H. Bohr. Chelsea, N.Y., 1947.
[5] Sobre la periodicidad de sumas de seales discretas peridicas. A. Restrepo y L. Chacn. Memorias del simposio de
Tratamiento de Seales, Imgenes y Visin Artificial, Universidad de los Andes, Bogot, 1995.
[6] On the period of sums of discrete periodic signals. A. Restrepo y L. Chacn. IEEE Signal Processing Letter, vol. 5, No. 7, pp.
164-166, 1998.
[7] Anlisis de Fourier. J. Duoandikoetxea. Addison- Wesley,
Wilmington, 1995.
[8] Fourier Analysis. T.W. Krner. Cambridge University Press,
Cambridge, 1998.
449
450
REFERENCIAS
Nuevo
en cap V
FILTROS DISCRETOS
El enfoque tomado es matemtico, lo cual redunda en precisin y en generalidad. Se dan ejemplos resueltos as como ejercicios y problemas al final de
cada captulo. Un aspecto en el que este curso se diferencia de la mayora de
los cursos introductorios de seales es que se ha evitado el uso de funciones
delta de Dirac, para el caso continuo. Para esto, se usa la respuesta escaln de
los sistemas de convolucin continuos, la integral de Riemann-Stieltjes y las
series de Fourier.
El profesor Restrepo es Ingeniero Electrnico de la Pontificia Universidad
Javeriana, MSc. y Ph.D. de la Universidad de Texas en Austin y, actualmente, es
profesor investigador en el departamento de Ingeniera Elctrica y Electrnica
de la Universidad de los Andes, en Bogot. Sus reas de investigacin incluyen
el tratamiento no lineal y estadstico de seales.
ndice de figuras
1.
2.
3.
4.
. . . . . . . . . . . . . . . . . . . . . .
. . . . . .
Una relacin simtrica . . . . . . . . . . . . . . . . . . . . . . . . . . .
El producto cartesiano de A y A
5.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
10
12
14
15
16
16
18
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Clasificacin de conjuntos segn su cardinalidad . . . . . . . . . . . . . .
Cada lnea representa un nmero racional . . . . . . . . . . . . . . . . . .
34
35
40
. . . . . . . . . . . . .
Diagrama de una funcin con dominio A y rango B . . .
f es una sobreyeccin de A en B, pero no una inyeccin .
f es una inyeccin de A en B, pero no es sobreyectiva . .
Una biyeccin . . . . . . . . . . . . . . . . . . . . .
Un ordenamiento circular para {0, 1, 2} . . . . . . . . .
Un orden parcial, no conectado
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Una inyeccin del conjunto de los nmeros pares al conjunto de los nmeros
naturales
13.
14.
15.
16.
6.
7.
8.
9.
10.
11.
12.
7
8
9
. . . . . . . . . .
Los puntos indicados hacen parte de la clase de equivalencia de x . . . . . . .
453
43
55
17.
18.
kxk . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
El intervalo de
S1
. . . . .
La solucin grfica de las races reales de una cbica . . . . . .
El modelo de coordenadas rectangulares: el producto R1 R1 .
El modelo de coordenadas polares:(0, ) S1 {} . . . . . . .
El modelo geomtrico cono para las coordenadas polares . . . .
El modelo geomtrico cono para las coordenadas polares . . . .
infinita de intervalos de la lnea de los nmeros reales
19.
20.
21.
22.
23.
24.
.
.
.
.
.
.
. . . . . . . . . . . . . . . . . . . .
La fase de z es igual a la fase de z ms (o menos) . . . . .
Una forma de definir la funcin arctan: con rango (/2, /2) .
Ilustracin de la frmula ||z| |w|| |z w| . . . . . . . . . .
Representacin geomtrica del producto complejo . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
56
61
63
63
64
64
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
67
68
69
70
70
.
.
.
.
.
.
.
.
.
.
.
.
.
.
71
72
73
75
76
79
86
87
89
99
100
100
102
104
. . . . . . . . . . . . . . . . . . . . . . . . . . .
La proyeccin estereogrfica de x S1 (/2] es p(x) R1 . . . .
La proyeccin estereogrfica . . . . . . . . . . . . . . . . . .
Proyeccin estereogrfica de la circunferencia S1 a la lnea R1 . . .
conjugado
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.
42.
.
.
.
.
.
.
25.
26.
27.
28.
29.
55
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
s es una sucesin unilateral . . . . . . . .
.
La sucesin a . . . . . . . . . . . . . .
.
1
Funcin para definir una mtrica en S . .
.
No tiene epgrafe . . . . . . . . . . . .
.
El criterio de la integral para mostrar divergencia .
El criterio de la integral para mostrar convergencia .
No tiene epgrafe . . . . . . . . . . . . . . . .
No tiene epgrafe . . . . . . . . . . . . . . . .
Orden para productos elementales positivos
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
43.
44.
45.
46.
. . . . . . . . . . . . . . . . . . . . . . . . . . . 104
. . . . . . . . . . . . . . . . . . . . . . . 106
n
Las funciones x para n = 1, 2, 4, 6 y la funcin lmite . . . . . . . . . . . . . 107
Grfica de la funcin x
47.
48.
49.
50.
51.
52.
. . . . . . . . . . . .
La ecuacin de la tangente est basada en la funcin lineal (x) = f 0 (x)(t x) .
El ngulo con vrtice z0 y el ngulo con vrtice f (z0 ) son iguales . . . . . .
El conjunto {(n, k) N N; k n} . . . . . . . . . . . . . . . . . . . . .
Una funcin escaln desplazada . . . . . . . . . . . . . . . . . . . . .
Teorema fundamental del clculo . . . . . . . . . . . . . . . . . . . . .
Promedio de una funcin . . . . . . . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
108
112
113
118
124
126
128
1.1.
1.2.
1.3.
1.4.
1.5.
1.6.
1.7.
1.8.
1.9.
1.10.
1.11.
1.12.
1.13.
1.14.
. . . . . . . . . . . . . . . . . . . . .
Una seal discreta de longitud finita . . . . . . . . . . . . . . .
Una seal continua. (La funcin correspondiente no es continua) . .
Una funcin del conjunto 4 al conjunto de los nmeros reales . . . .
Una seal discreta con rango complejo . . . . . . . . . . . . . .
Un sistema es una funcin entre conjuntos de seales . . . . . . .
Representacin de un sistema como una caja negra . . . . . . .
La seal de voltaje h se aplica a una resistencia de un ohmio . . . .
Grfica 1t , 1 , t12 . . . . . . . . . . . . . . . . . . . . . . . . .
t
La seal cos 5t para t [0,1, 0,6] . . . . . . . . . . . . . . . .
La seal de la figura anterior, en valor absoluto . . . . . . . . . .
La seal dela figura anterior, al cuadrado . . . . . . . . . . . . .
Tubo de rayos catdicos . . . . . . . . . . . . . . . . . . . . .
La seal | sen(t)| y una onda tringulo F(t) . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
153
157
160
161
162
168
169
174
176
177
177
178
180
182
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
2.3.
2.4.
2.5.
2.6.
2.7.
2.8.
2.9.
2.10.
2.11.
2.12.
2.13.
. . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . .
Las seales e y h a ser convolucionadas . . . . . . . . . .
Para calcular s(1) es necesario multiplicar h() y e(1 ) . .
El valor de s(1) es el rea sombreada . . . . . . . . . . .
g= f h . . . . . . . . . . . . . . . . . . . . . . . . .
Una seal con salto en t = 0 . . . . . . . . . . . . . . .
Un sistema de convolucin . . . . . . . . . . . . . . . .
Circuito RC pasa bajas de primer orden . . . . . . . . . .
Circuito RC pasa altas . . . . . . . . . . . . . . . . . .
Las seales discretas f y g a ser convolucionadas
g2k
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
191
191
194
194
195
196
198
200
202
207
2.14.
2.15.
2.16.
2.17.
2.18.
. . . . . . . . . . . . .
Circuito RC pasa bajas e primer orden en el plano s .
Dos sistemas de convolucin en cascada . . . . . .
Un circuito RLC pasa banda de segundo orden . . .
Respuesta escaln . . . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
209
212
215
219
219
220
3.1.
3.2.
3.3.
3.4.
3.5.
3.6.
3.7.
3.8.
3.9.
3.10.
3.11.
3.12.
. . . . . . . . . . . . . . .
. . . . . . . . . . .
Una seal con periodo 3, otra con periodo 4 y su suma .
Las seales sen t, 1/2 sen 2t y su suma . . . . . . . . .
La seal g(t) = sen t + sen t . . . . . . . . . . . . .
La seal sen t . . . . . . . . . . . . . . . . . . . .
La seal cos t . . . . . . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
225
226
230
231
231
235
236
240
240
241
242
245
. . . . .
. . . . . .
Onda tringulo de frecuencia 60, con nivel DC de 6 voltios . .
La seal senc enmascarada, SE4 (x), de orden 4 . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
3.13.
3.14.
3.15.
3.16.
3.17.
3.18.
3.19.
3.20.
3.21.
3.22.
3.23.
4.1.
4.2.
4.3.
4.4.
4.5.
4.6.
4.7.
4.8.
4.9.
. . . . .
. . . . . . . . . . . .
La suma parcial seales sen t + 1/3 sen 3t . . . . . . . .
La suma parcial seales sen t + 1/3 sen 3t + 1/5 sen 5t . . .
.
.
.
.
La suma de los tres nmeros complejos mostrados es cero .
Una seal peridica discreta . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
. . . . . . . . . . . . . . . . . .
b. El resultado de agregar intervalos nulos entre periodos .
c. el resultado en el lmite. . . . . . . . . . . . . . . . .
Pulsos . . . . . . . . . . . . . . . . . . . . . . . . .
Una seal pulso . . . . . . . . . . . . . . . . . . . . .
La funcin senc x . . . . . . . . . . . . . . . . . . . .
La transformada de Fourier . . . . . . . . . . . . . . .
La funcin 1 ||
r para varios valores de r . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Las seales sen t, 1/3 sen 3t, 1/5 sen 5t, 1/7 sen 7t
La suma parcial seales sen t + 0
.
.
.
.
.
.
.
.
La suma parcial seales sen t + 1/3 sen 3t + 1/5 sen 5t + 1/7 sen 7t
. . . . . . . .
La onda tringulo par de periodo 2 . . . . . . . . .
Onda sierra . . . . . . . . . . . . . . . . . . .
Seal generada por la mquina de Michelson
Dos seales peridicas con periodo 2
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
246
246
246
247
247
249
249
250
253
255
259
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
269
269
269
271
274
276
281
284
.
.
.
.
.
.
.
.
.
286
291
293
294
295
295
296
297
299
. . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . .
Representacin en tiempo y en frecuencia de un sistema de convolucin . .
Descomposicin espectral de un rayo de luz . . . . . . . . . . . . . . .
La funcin de transferencia de un pasa bajas ideal . . . . . . . . . . . .
Esta seal se aplica al filtro pasa bajas ideal . . . . . . . . . . . . . . .
Transformada de Fourier de la seal en la figura anterior . . . . . . . . .
Falta epgrafe . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
NO TIENE EPGRAFE . . . . . . . . . . . . . . . . . . . . . . . . .
amplificadores.)
4.10.
4.11.
4.12.
4.13.
4.14.
4.15.
4.16.
4.17.
.
.
.
.
Un sistema de convolucin
.
.
.
.
.
.
.
.
.
303
303
304
304
305
305
305
306
307
307
308
310
311
311
4.32. a. Seal de entrada b. Seal de salida del filtro ideal c. Seal de salida del filtro RC 314
4.33. Do=C (261.626Hz), Re=D (293.665Hz), Mi=E(329.628Hz), Fa=F (349.228Hz),
Sol=G (391.995Hz), La=A (440.0Hz), Si=B (493.883Hz) . . . . . . . . . . . . 317
4.34. Representacin en decibeles de la funcin de transferencia . . . . . . . . . . 320
4.35. La magnitud y la fase de un pasa todo con fase lineal . . . . . . . . . . . . . 321
4.36. Un pasa todo con fase no lineal . . . . . . . . . . . . . . . . . . . . . . . 322
4.37. a) La seal de entrada. b) La seal de salida (alteracin de la forma de onda) . 323
4.38. Onda tringulo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 324
4.39. Espectro de la seal s(t) . . . . . . . . . . . . . . . . . . . . . . . 329
4.40. Espectro de la seal modulada f (t) . . . . . . . . . . . . . . . . . 330
4.41. Espectro de la seal h(t) . . . . . . . . . . . . . . . . . . . . . . . 330
4.42. Espectro de la seal resultante b(t) . . . . . . . . . . . . . . . . . 330
4.43. Modulacin SSB usando la transformada de Hilbert . . . . . . . 331
negativo de la respuesta escaln de un pasa bajas de primer orden, que no salta
4.44.
4.45.
4.46.
4.47.
4.48.
4.49.
4.50.
4.51.
4.52.
4.53.
4.54.
4.55.
4.56.
4.57.
4.58.
4.59.
4.60.
4.61.
4.62.
4.63.
4.64.
4.65.
. .
. . . . . . . . . . . . . . . . . . . . . .
Una onda cuadrada . . . . . . . . . . . . . . . . . . . . . . .
Respuesta del circuito RC a una onda cuadrada . . . . . . . . . .
.
.
.
.
.
.
.
.
333
334
335
336
Grfica de la relacin de los voltajes pico a pico de entrada y salida del circuito RC 338
Respuesta estable para T << . . . . . . . . . . . . . . . . . . . . . . . 338
Respuesta estable para T >> . . . . . . . . . . . . . . . . . . . . . . . 338
Seal peridica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 339
Circuito . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 340
Magnitud y fase . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 341
Rectifivador de media onda y pasabajas RC . . . . . . . . . . . . . . . . . 342
Pasabajs RC con frecuencia de corte 1Hz y entrada con nivel DC . . . . . . . 343
Seal triangular . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 343
La convolucin de dos pulsos es un tringulo . . . . . . . . . . . . . . . . 343
Pasa-altas = identidad - pasabajas . . . . . . . . . . . . . . . . . . . . . 344
Problema 20 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 345
La seal mostrada se filtra con el sistema de convolucin . . . . . . . . . . . 345
Sistema . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 348
Onda cuadrada . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 349
Circuito impedancias . . . . . . . . . . . . . . . . . . . . . . . . . . . 351
Seal de entrada . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 352
Sistema . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 352
5.1.
5.2.
5.3.
5.4.
5.5.
5.6.
5.7.
La funcin 1 + 2 cos
. . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . .
La fase de la funcin 1 + 2 cos . Como = , la funcin tambin es impar .
La respuesta impulso de un promedio mvil . . . . . . . . . . . . . . .
La respuesta impulso de un promedio mvil causal . . . . . . . . . . . .
En los dos casos se ha sombreado el rea sobre un periodo . . . . . . . . .
La circunferencia de Fourier para seales discretas . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
359
360
360
362
363
364
365
365
366
367
367
368
368
369
370
371
371
372
374
374
378
382
383
386
. . . . . . . . . . . . . . .
5.26. Las Xk s son muestras a escala de S() . . . . . . . .
5.27. La transformada de una seal continua . . . . . . . .
5.28. La transformada de la seal muestreada correspondiente
5.29. Una seal discreta de duracin finita . . . . . . . . .
5.30. Las seales g0 senc (t), g1 senc (t 1) y g2 senc (t 2) . . .
5.31. La seal f (t) = g0 senc (t) + g1 senc (t 1) + g2 senc (t 2)
5.32. Serie de Fourier . . . . . . . . . . . . . . . . . . .
peridica con periodo 2
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
388
391
394
394
397
397
398
400
5.33.
5.34.
5.35.
5.36.
5.37.
5.38.
. . .
. . . . . . . . . . . .
Suma de Fourier . . . . . . .
DFT . . . . . . . . . . . .
.
.
.
.
Dos seales discretas de longitud finita .
Un filtro de retardo . . . . . . . . . .
.
.
.
.
.
.
400
401
401
401
408
410
. . . . . . . . . . . . . . . . . . .
6.2. La seal s y una versin desplazada s de sta . . . . . . . . . . . . . . . .
6.3. Si la seal de entrada a un sistema invariante se desplaza, la salida se desplaza
en la misma forma. . . . . . . . . . . . . . . . . . . . . . . . . . . . .
6.4. Las tres seales mostradas x, y, z para t < , coinciden . . . . . . . . . . . .
6.5. La seal impulso discreta (o delta de Kronecker) . . . . . . . . . . . . . . .
6.6. Una seal discreta . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
6.7. La descomposicin de la seal en la figura 6.6 como una suma de seales impulso
pesadas y desplazadas. . . . . . . . . . . . . . . . . . . . . . . . . . .
6.8. La respuesta impulso de un promedio mvil . . . . . . . . . . . . . . . .
6.9. Una seal simple . . . . . . . . . . . . . . . . . . . . . . . . . . . .
425
427
Transformada de Fourier
DTFT
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
6.1. Un filtro mediana tiene tres bloques principales: uno que toma segmentos de la
seal de entrada, otro que ordena los datos en los segmentos y otro que toma el
valor ordenado intermedio, o mediana
428
429
439
439
440
441
442
ndice de cuadros
1.
2.
3.
4.
.
.
2
4
j 20
j
j
Complete la tabal mostrada e = e 3 , f = e 3 , g = e 3 . . . . . . . . . . .
La tabla de multiplicacin de los cuaternios. . . . . . . . . . . . . . . . .
La tabla de una posible operacin entre los elementos del conjunto {a, b, c, d}
La tabla de una operacin entre los elementos del conjunto finito , , , ,
.
.
.
.
21
22
23
25
462