Procesamiento de señal - Matlab

El programa
Lo que haremos será lo siguiente:
1. Cargaremos el audio en MATLAB
2. Separaremos los canales izquierdo y derecho (Trabajamos con una señal en estéreo)
3. Aplicaremoslatransformada de Fourier a la señal.
4. Representaremos gráficamente los resultados
5. Identificaremos las frecuencias más destacadas
6. Hallaremos las notas musicales
 Empezamos por cargar laseñal en MATLAB, para ello usaremos dos variables, una que contendrá la señal en sí (wav), y otras que nos dará la frecuencia de muestreo (fs).
El segundo comando lo que hace estransponer lamatriz que contiene los datos, para más comodidad, pero no es esencial su comprensión.
[wav fs] = wavread(‘nota’);
wav=transpose(wav);
Una vez hecho esto separaremos los canales derechoeizquierdo en dos variables distintas, para que así podamos aplicar más tarde Fourier sin problemas.
x1 para el canal izquierdo y x2 para el derecho.x1=wav(1,1:length(wav));
x2=wav(2,1:length(wav));Ahora que tenemos cada canal por separado, nos interesa conocer la frecuencia de las ondas que componen el sonido; porque el sonido está formado por una suma de muchas ondas sinusoidalespuras. Paraello usaremos una operación matemática llamada Transformada de Fourier, que a grandes rasgos lo que hace es hallar que frecuencias son las más habituales en una onda.
Usaremos paraello el comandofreqz de MATLAB que nos toma nuestra señal y nos da los datos ya procesados en frecuencias. Peso de las frecuencias (X) y rango de frecuencias (f) para cada canal de audio (1izquierdo, 2 derecho)[X1 f1]=freqz(x1,1,length(wav),fs);
[X2 f2]=freqz(x2,1,length(wav),fs);
Ahora representaremos la señal en un par de ventanas. Dibujaremos arriba la gráfica Amplitud en funcióndel tiempo y abajo la dePeso en función de la frecuencia.
Es de especial mención el hecho de que el eje abscisas (eje OX) de la segunda gráfica no representa el tiempo, si no la frecuencia.