Articles of fft

Transformación de Fourier 3D comparativamente lenta de Python numpy

Para mi trabajo, necesito realizar transformaciones de Fourier discretas (DFT) en imágenes grandes. En el ejemplo actual, necesito un FT 3D para una imagen de 1921 x 512 x 512 (junto con FFT 2D de 512 x 512 imágenes). En este momento, estoy usando el paquete numpy y la función asociada np.fft.fftn () . El […]

¿Cuál es la forma más rápida / eficiente de recorrer una gran colección de archivos y guardar una gráfica de los datos?

Así que tengo este progtwig que recorre aproximadamente 2000 archivos de datos, realiza una transformación de Fourier, traza la transformación y luego guarda la figura. Parece que cuanto más se ejecuta el progtwig, más lento parece ser. ¿Hay alguna forma de hacerlo funcionar más rápido o más limpio con un simple cambio en el código […]

Espectro de potencia y autocorrelación de datos en Numpy

Estoy interesado en calcular el espectro de potencia de un sistema de partículas (~ 100,000) en el espacio 3D con Python. Lo que he encontrado hasta ahora es un grupo de funciones en Numpy ( fft , fftn , ..) que calculan la transformada de Fourier discreta, de la cual el cuadrado del valor absoluto […]

Python numpy equivalente de bandpower () de MATLAB

Estoy utilizando el método ftt.ftt () de Python numpy para generar la transformada de Fourier de una señal. Sin embargo, quiero calcular la potencia de banda en un rango de frecuencias. MATLAB tiene el método bandpower (x, fs, freqrange), estoy tratando de simular específicamente esta syntax de la función. Fuente: https://www.mathworks.com/help/signal/ref/bandpower.html No parece que numpy […]

Cómo usar numpy con portaudio para extraer bajos, agudos medios

Como en este ejemplo, ¿Cómo extraer información de frecuencia de un flujo de audio de entrada (usando PortAudio)? Tengo curiosidad por portaudio y numpy … No estoy 100% seguro de fft, ¿cómo puedo pasar Numpy un trozo y recuperar tres valores de -1.0 a 1.0 para bajos, medios y agudos? No me importa si esto […]

¿Cómo escribir fftshift y ifftshift en R?

En numpy , tenemos las siguientes funciones: import numpy from numpy.fft import fft2, ifft2, fftshift, ifftshift Me gustaría reescribir estas funciones en R. fft en R funciona igual que fft o fft2 en python. También para ifft2 , tenemos que hacer fft(,inverse=T) Ahora me gustaría saber cómo reescribir las funciones fftshift y ifftshift (para matrices), […]

Filtrado del espacio de Fourier

Tengo una serie temporal vectorial real x de longitud T y un filtro h de longitud t << T. h es un filtro en el espacio de Fourier, real y simétrico. Es aproximadamente 1 / f. Me gustaría filtrar x con h para obtener y. Supongamos que t == T y FFT de longitud T […]

¿Cuál es la diferencia entre numpy.fft y scipy.fftpack?

¿Es el último solo un sinónimo del primero, o son dos implementaciones diferentes de FFT? ¿Cuál es mejor?

Cálculo del desplazamiento cuadrático medio utilizando python y FFT

Dada una matriz bidimensional, donde cada fila representa el vector de posición de una partícula, ¿cómo calculo el desplazamiento cuadrado medio de manera eficiente (usando FFT)? El desplazamiento cuadrado medio se define como donde r (m) es el vector de posición de la fila m, y N es el número de filas.

Python: diseño de un filtro de series de tiempo después del análisis de Fourier

Tengo una serie temporal de datos de temperatura de 3 horas que he analizado y he encontrado el espectro de potencia para usar el análisis de Fourier. data = np.genfromtxt(‘H:/RData/3hr_obs.txt’, skip_header=3) step = data[:,0] t = data[:,1] y = data[:,2] freq = 0.125 yps = np.abs(np.fft.fft(y))**2 yfreqs = np.fft.fftfreq(y.size, freq) y_idx = np.argsort(yfreqs) fig = […]