3.3.

Transformada de Fourier de un Escalón Unitario

La Transformada de Fourier de un Escalón Unitario se define como
Z ∞
F[u(t)] = u(t)e−j2πf t dt
−∞

Debido a que la función se anula en el intervalo [∞, 0[ y toma valor 1 en el intervalo

restante, la expresión queda
Z ∞
F[u(t)] = e−j2πf t dt
0

Mediante integración directa resulta

" #∞
−1 −j2πf t 1 δ(f )
F[u(t)] = e = +
j2πf 0
j2πf 2

3.4. Limitación en ancho de banda

Para demostrar la igualdad
sen(at)
x(t) = x(t) , considerando convenientemente que a = 2πfa > wc = 2πfc
πt
Se procede a calcular la Transformada de Fourier a ambos lados de la ecuación, que-
dando !
f
X(f ) ∗ rect = X(f )
2fa
Debido a que
sen(at) sen(2πfa t) 2fa
= = 2fa sinc(2fa t)
πt πt 2fa
Por otro lado, como la función rect posee un largo τ = 2fa , donde fa > fc , ésta poseerá
valor unitario, por lo que se estarı́a cumpliendo la igualdad

4. Simulaciones en Matlab
4.1. Cálculo matemático de la Serie de Fourier
Sea la señal y(t) definida como
∞
X
y(t) = p(t − k)
k=−∞

Los coeficientes asociados a la Serie de Fourier Trigonométrica están dados por

Z 1
2
Z 1
2
Z 0 √ Z 1
2 √ 2
A0 = y(t)dt = (p(t + 1) + p(t))dt = t + 1 dt + t dt =
−1
2
−1
2
−1
2
0 3

12
 Z 1 Z 1
2 2
An = 2 −1
y(t)cos(2πnt)dt = 2 −1
(p(t + 1) + p(t))cos(2πnt)dt
2 2
Z 0 √ Z 1
2 √
An = 2 −1
t + 1 cos(2πnt)dt + 2 t cos(2πnt)dt
2
0

√  √  √   √ 
" ! !
1 
An = 2 C 2 n − C 2n sen(2πn) + S 2n − S 2 n cos(2πn)
4πn3/2
√ √  !#
√ 2 sen(πn) S 2n
!# "
1
+ 2n sen(πn) +2 −
4π n n3/2
Z 1 Z 1
2 2
Bn = 2 −1
y(t)sen(2πnt)dt = 2 −1
(p(t + 1) + p(t))sen(2πnt)dt
2 2
Z 0 √ Z 1
2 √
Bn = 2 −1
t + 1 sen(2πnt)dt + 2 t sen(2πnt)dt
2
0

√  √   √  √ 
" ! !
1 
Bn = 2 C 2 n − C 2n cos(2πn) + S 2 n − S 2n sen(2πn)
4πn3/2
√ √
√ √
!!# " #
C( 2n) − 2n cos(πn)
+ n 2 cos(πn) − 2 +2
4πn3/2
Donde Z x
C(x) = cos(t2 )dt
0
Z x
S(x) = sen(t2 )dt
0
Son las Integrales de Fresnel que se obtienen al computar An y Bn en Wolfram Alpha

4.2. Caracterı́sticas de la Señal y(t)

La descripción de la señal y(t) a través el tiempo se muestra a continuación

13
Figura 1. Señal y(t) de 0 a 10 segundos

14
Figura 2. Código de Matlab asociado a la descripción de la señal y(t)

Donde
X k - Iterador
X M - Número de Muestras Discretas
X T M - Tiempo de Muestreo
X F M - Frecuencia de Muestreo
X T - Vector de Tiempo

15
 X f - Vector de Función y(t)

Figura 3. Código de Matlab asociado a la descripción de la fase del espectro de

la señal y(t)

Donde

X F - Transformada de Fourier de la función y(t)

Figura 4. Fase del espectro de la señal y(t)

16
4.3. Aproximación de la Señal y(t) a través de 3 armónicos

√ √
S( 2) − S(2) S( 2)
A1 = − = 0,0589759 − 0,113632 = −0,0546461
2π 2π
√
S(2) − S(2 2) S(2)
A2 = √ − √ = −0,002507 − 0,0193239 = −0,0218309
4 2π 4 2π
√ √ √
S( 6) − S(2 3) S( 6)
A3 = √ − √ = 0,00675802 − 0,019264 = −0,01250598
6 3π 6 3π
√ √ √ √
−C(2) + C( 2) + 2 + 2 C( 2) + 2
B1 = + = −0,549857 + 0,309255 = −0,240602
2π 2π
√ √ √ √
− 2C(2) + 2C( 2) − 4 + 2 2 C(2) − 2
B2 = + √ = −0,0462009−0,0850656 = −0,1312665
8π 4 2π
√ √ √ √ √ √ √ √
− 3C(2 3) + 3C( 6) + 6 + 3 2 3C( 6) + 3 2
B3 = − + = −0,181407+0,0905445 = −0,09
18π 18π

Figura 5. Código de Matlab asociado a la descripción de la Serie de Fourier

asociada a la señal y(t), incluyendo solamente los tres primeros armónicos

Donde

X sf - Aproximación mediante la Serie de Fourier Trigonométrica de la función y(t),

con 3 armónicos, los cuales fueron obtenidos mediante los resultados anteriores

17
 Figura 6. Comparativa entre la señal y(t) original, y la aproximada mediante
Serie de Foyrier

4.4. Aproximación más precisa de la Señal y(t)

A4 = −0,000647497−0,0077184 = −0,00836637 , B4 = −0,0232519−0,0464097 = −0,0696616

A5 = 0,00244113−0,00854631 = −0,00610518 , B5 = −0,108739+0,0521784 = −0,0565606

A6 = −0,000289766−0,00442155 = −0,004711316 , B6 = −0,0155216−0,0321247 = −0,0776463
A7 = 0,0012467−0,00502674 = −0,00378004 , B7 = −0,0776494−0,0364674 = −0,1141168
A8 = −0,000163401−0,00295779 = −0,003121191 , B8 = −0,0116467−0,0246291 = −0,0362758
A9 = 0,000754489−0,00338916 = −0,002634671 , B9 = −0,0603871+0,0279639 = −0,0324232
A10 = −0,0001047−0,00215851 = −0,00226321 , B10 = −0,00931941−0,0199971 = −0,02931651
A11 = −0,001971923 , B11 = −0,0267576

18
 A12 = −0,0017384549 , B12 = −0,02461261
A13 = −0,001547891 , B13 = −0,0227886
A14 = −0,0013899134 , B14 = −0,021218
A15 = −0,001257213 , B15 = −0,0198515

Figura 7. Código de Matlab asociado a la descripción de la función sf (t) que

realiza una aproximación de la señal y(t) a través de la Serie de Fourier Trigo-
nométrica

Donde

X sf 2 - Aproximación mediante la Serie de Fourier Trigonométrica de la función y(t),

con 6 armónicos, los cuales fueron obtenidos mediante los resultados anteriores

X sf 3 - Aproximación mediante la Serie de Fourier Trigonométrica de la función y(t),

con 10 armónicos, los cuales fueron obtenidos mediante los resultados anteriores

X sf 4 - Aproximación mediante la Serie de Fourier Trigonométrica de la función y(t),

con 15 armónicos, los cuales fueron obtenidos mediante los resultados anteriores

19
 Figura 8. Gráfico comparativo entre la señal y(t) y su aproximación, sf (t)

Con el uso de 10 armónicos la aproximación es bastante precisa. Sin embargo, con los
resultados de 15 armónicos la aproximación no nota cambios significativos, por lo que
ese serı́a el criterio para escoger la cantidad de armónicos a utilizar.

20