Pds-fourier

Solo disponible en BuenasTareas
  • Páginas : 2 (348 palabras )
  • Descarga(s) : 4
  • Publicado : 9 de diciembre de 2009
Leer documento completo
Vista previa del texto
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...
tracking img