PDF de programación - Práctico 1: Análisis Espectral usando MATLAB

PPRRÁÁCCTTIICCOO  11::  AAnnáálliissiiss  EEssppeeccttrraall  uussaannddoo  MMAATTLLAABB                                                            



OOBBJJEETTIIVVOOSS DDEELL TTRRAABBAAJJOO
Al finalizar el trabajo, el estudiante será capaz de:
► Ingresar archivos de sonido (en particular de
señales de voz) en formato .wap (formato de
sonido WAVE de Microsoft) al espacio de
trabajo de MATLAB, mediante el comando
wavread.


► Usar MATLAB para representar señales de
voz en el dominio del tiempo (formas de onda)
y en el dominio de la frecuencia (espectro),
mediante el comando plot.



► Usar

la

los
componentes frecuenciales de las muestras.

fft para obtener

función


GGEENNEERRAACCIIÓÓNN DDEE AARRCCHHIIVVOOSS

Crear un archivo de audio .wav
1. Abra el programa grabador de sonidos de



Windows (Figura 1).


2. En Archivo/Propiedades, seleccione el
formato PCM con los siguientes atributos: 8
kHz, 16 bits, Mono.


3. Delante de un micrófono, pronuncie su

nombre.


4. Repita el anterior punto para una nota
musical de cualquier instrumento, durante el
tiempo que tardó en pronunciar su nombre.


5. Guarde

los archivos generados con

la
extensión .wav ,a la que se llamará ‘dirección
del archivo de audio’, ejemplos: ‘denise.wav’
o ‘guitarra.wav’



Agosto 2010 – www.coimbraweb.com  



Figura 1. Programa grabador de sonido.



Cargar el archivo de audio a MATLAB.

1. Mediante el comando x=wavread(‘dirección
del archivo’); ingrese al espacio de trabajo de
MATLAB, cargando el archivo de audio al
vector x, tanto para su nombre como para la
nota musical.


2. Puede comprobar la gráfica en el dominio del
tiempo con plot(x); parecida a la de la Figura
2 para su nombre como para
la nota
musical:



Figura 2. Pronunciación de su nombre.


AANNÁÁLLIISSIISS EESSPPEECCTTRRAALL DDEE LLAASS MMUUEESSTTRRAASS
Utilice los siguientes comandos para el archivo
correspondiente, tanto para su nombre como
para la nota musical.

x=wavread(‘dirección del archivo’);
% archivo de audio para su nombre.
plot(x) % gráfica en el dominio del tiempo.

Y=fft(x); % transformada rápida de Fourier.
A=Y.*conj(Y); % potencia de la señal.

 11  

TTIICCOO  11::  AA
PPRRÁÁCCTT
AAnnáálliissiiss

EEssppeeccttrr

rraall  uussaann

nddoo  MMAATT

TTLLAABB          

                           

                   

uce mayor

100:3000);
% espectro
t(f,A(1:2901
)); % gráfica

la frecue
obtendrá u
n lo cual se
la Figura 3.
.

o de frecuen
a en el dom
encia.
na gráfica p

f=(1
plo


Co
de


ncia.
inio de

parecida a la
a


su nombre.

ponentes frec

cuenciales de

ación para la

a nota music

cal

pita la opera



Fig
gura 3. Comp

Re

IINNFF
Org
pre
mín


FFOORRMMEE
s
resultad
los
ganice
n un inform
esentarlos e
nimo, lo sigu
uiente:

dos
obten
me que cont

para
idos
a
tenga, como
o

ntación en e
s que resu
de tocar una
de duración d
entación e
de cada señ

el dominio de
ulten de pro
a nota musi
da cada señ
en el dom
ñal.

el tiempo de
e
onunciar su
u
cal.
ñal.
minio de

la
a

La represen
las señale
nombre y d

El tiempo d
La
repres
frecuencia



E
El ancho de
La frecuenc
L
p
potencia en

ada señal.
banda de c
cia a la cu
ual se prod

cada señal.


MMAATT

Se e
MAT

TTEERRIIAALL DDEE

RREEFFEERREENNCC

CCIIAA

entregará a
TLAB Tools.

los estudian


ntes el mater

rial

CCRROO

Pres
instr

AA
OONNOOGGRRAAMMAA
sentación
del
rucciones en

traba
ajo: de a
n el sitio web
b.

acuerdo a



sto 2010 – www

w.coimbraweb.c

om  

Agos


 22