sin(t)
+
sin(3* 1) / 3
+
sin(5* 1) / 5
+
sin(7* 1) / 7
+
sin(9* 1) / 9;
plot(l,
y)
,
,
,
,
,
,
,
,
rrl:n.,'
",.:r,
\7 :•.•. ; .... ; .... ;..
,
0.6
....
~
....:.. . . ..... ..:....
~
.... ; ... : .... : ....
~
.. .
"
"""
,
,
,
,
,
,
0.4 ....:.....:.....
····:····~·· · -1 - --t.---f----~-
--
,
,
,
,
,
,
,
,
,
0.2
____
~
____.. . .. ..: ___ ..: ____
~
____ ; ___ : ____ : ____
~
__ _
,
,
,
,
,
,
,
,
,
"
,,""
O ____ :__ • • _:_____: ____:____ , ____, ___ ; ____ ; __
._¡. __ .
,
,
,
,
,
,
,
,
,
-0 .2
,
""'"
--
..... .......
_-_
...
--_
..
_. _-,---_ . _._ . _--- "----,-
--
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
-- ..... . ..,-----, --- ...----,---- , ---'---- r---- r- --
-0.4
,
,,"'"
,
,
,
,
,
,
,
,
,
____ '- __ •• '... _-'___ -' ____ , ____ ,
___ • ____ L ____ L __ _
-0 .6
,
,
,
"
'"
,
,
,
"
'"
-0 .8
,
,
,
,
,
,
----
..
--.-
.....
-....
--_
...
_-- .-._- . --
,
,
,
,
,
,
,
,
,
,
,
,
- 1
~o
-"""--:-,
-:-3
-47-~5:--C6:-~7:--C8:--C9,---J1O
Adición de tres armónicos
t
=
0:.02:3.14;
y
=
zeros(lO, max(size(t))); % abriendo espacio para almacenar los
valores de los armónicos
x
=
zeros(size(t));
% programa para generar los armónicos de una señal cuadrada
for
k
=
1:2:19
end
x
=
x
+
sin(k* I) / k;
y((k+
1) / 2,:)
= x;
ANÁLISIS ESPECTRAL
Como vimos en el capítu lo IlI, una herramienta poderosa en el análisis de
señales es la transformada de Fourier.
A
continuación veamos un ejemplo
de la obtención de ésta utilizando el algoritmo
FFT:
1
=
0:1 / 1000:0.25;
O/o
vector con datos muestreados
al khz
x
=
sin(2* pi* 60' 1)
+
sin(2' pi ' 120' 1);% suma de sella les d e 60
y
120 hz
y l
=
x
+
2* randn(size(I));
% adicionando ruido con d esv iación
% es tándar d e 2
subplot(211)
plot(t(l:100),
yl(I :100)),
title('señal con ruido en e l dominio del tiempo' );
165
1...,155,156,157,158,159,160,161,162,163,164 166,167,168,169,170,171,172,173,174,175,...196