INSTITUTO TECNOLÓGICO DE MÉRIDA
Ingeniería Electrónica
Análisis de la Transformada de Fourier
Br. Cesar Ernesto Meneses Lazo
Dr. Jesús Sandoval Gío
2/12/09
Obtener la señal yvalores de un archivo de audio “.wav” y efectuar el análisis de Fourier para obtener sus graficas con repescto a las muestras y frecuencias, asi como también su grafica con shift etc. Para obtener losvalores del archivo de audio aplicamos la función “wavread”. >> Fs = 44000; >> [x,Fs]=wavread(‘jungle4’); A continuación se grafica para observar su comportamiento
>> Ts= 1/500; >> N=length(x); >> n =0:N-1; >> stem(n,x);
Luego de adquirir los datos dentro de la variable “x” podemos aplicar la transformada rápida de Fourier >> X = abs(fft(x)); En este caso para graficar asignamos al ejehorizontal la variable “k”con los mismos valores de “n” esto nos servirá para indicar que es una grafica de la transformada de Fourier y hablamos de las muestras en ese plano. >> K=n; >> figure; plot(k,X); Lo anterior nos sirve para observar la transformada de Fourier con respecto a las muestras, para poder observarlo con respecto a la frecuencia debemos primero obtener la resolución de frecuencia, lacual nos servirá para ver los diferentes valores a diferentes frecuencias (ejemplo: 1Hz,2Hz,etc.); y ya para poder graficarlo con respecto a las frecuencias por el numero de muestras. >> df = Fs/N;>> f=k*df; >> figure;plot(f,X);
Ahora para ver la grafica en forma de espejo, creamos “ks” para poder graficar desde -N/2 hasta N/2-1, para luego tener nuestra nueva “fs” mediante la (ks*df) parapoder después aplicar la función “fftshift”. >> ks = -N/2:N/2-1; >> fs = ks*df; >> Xs=fftshift(X); >> figure;plot(fs,Xs);
Ahora para reducir las amplitudes aumentadas por la transformada rápida deFourier lo dividimos entre el número de muestras y lo multiplicamos por dos debido a que son dos partes.
>> Xsc=(Xs/N)*2; >> figure;plot(fs,Xsc);
En conclusión podemos observar que los picos de…