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 […]
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 […]
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 […]
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 […]
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 […]
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), […]
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 […]
¿Es el último solo un sinónimo del primero, o son dos implementaciones diferentes de FFT? ¿Cuál es mejor?
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.
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 = […]