Acabo de empezar a escribir un progtwig para manipular un poco de audio en Python. Antes de escribir cualquier función de filtrado, quería hacer una prueba para comparar la señal de entrada con la señal de salida después de que la señal de entrada pasara por un rfft e irfft. Por alguna razón, el archivo […]
Tengo este fragmento de código y quiero trabajar de forma independiente con cada uno de los cuatro arreglos de micro que utilizo (Ps3 Eye Cam). Por ejemplo, quiero calcular la correlación cruzada entre el micrófono 1 y el micrófono 4 excluyendo 2 y 3. ¿Cómo puedo usar los canales que quiero? import audioop import alsaaudio […]
Estoy trabajando con algunos archivos de audio en Python usando los paquetes scipy.fftpack y scipy.io. Lo que esto implica es que tengo archivos wave que estoy importando, jugando con ellos utilizando las transformadas de Fourier y luego enviando a un nuevo archivo wave. Sin embargo, me estoy encontrando con problemas en los que después de […]
Me gustaría calcular un espectro de potencia utilizando Python3. De otro hilo sobre este tema obtuve los ingredientes básicos. Creo que debería ser algo como: ps = np.abs(np.fft.fft(x))**2 timeres = t[1]-t[0] freqs = np.fft.fftfreq(x.size, timeres) idx = np.argsort(freqs) plt.plot(freqs[idx], ps[idx]) plt.show() Aquí están los tiempos y x es el recuento de fotones. También he intentado: […]
Realizo el iFFT en un espectro simétrico (usando Python). ¿Por qué el resultado no es una señal de valor real sino que contiene valores complejos? # My symmetric spectrum spectrum = numpy.array( [1+1j,2+2j,3+3j,3-3j,2-2j] ) # Perform the iFFT print numpy.fft.ifft(spectrum) Salida: (2.2+0.2j) (-1.98979431354+0.2j) (0.59464641547+0.2j) (-0.74743281997+0.2j) (0.942580718037+0.2j)
Escribo una aplicación de procesamiento de imágenes que necesita hacer varias cosas y tiene que hacerlo lo más posible en tiempo real. La adquisición de los datos y su procesamiento se ejecuta en procesos separados (principalmente por motivos de rendimiento). Los datos en sí son bastante grandes (imágenes en escala de grises de 2MPix de […]
Estoy tratando de estimar la PSD de la variabilidad de la frecuencia cardíaca de una señal de ECG. Para probar mi código, extraje el intervalo RR de la base de datos de ECG de Fantasia . He extraído la señal que se puede acceder aquí . Para calcular el PSD, estoy usando el método de […]
Estoy tratando de entender la teoría del FTT y la convolución (correlación cruzada) y por esa razón he creado el siguiente código para entenderlo. El código es Matlab / Octave, sin embargo, también podría hacerlo en Python. En 1D: x = [5 6 8 2 5]; y = [6 -1 3 5 1]; x1 = […]
Miré en muchos ejemplos de scipy.fft y numpy.fft. Específicamente este ejemplo, el Análisis de Frecuencia de Scipy / Numpy FFT es muy similar a lo que quiero hacer. Por lo tanto, utilicé el mismo posicionamiento de subttwig y todo parece muy similar. Quiero importar datos de un archivo, que contiene una sola columna para que […]
Analicé los datos de sunspots.dat (abajo) usando fft, que es un ejemplo clásico en esta área. Obtuve resultados de FFT en partes reales e imaginarias. Luego intenté usar estos coeficientes (los primeros 20) para recrear los datos siguiendo la fórmula de la transformada de Fourier. Pensando partes reales corresponden a a_n e imaginarios a b_n, […]