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 = […]
Estoy tratando de volver a implementar una de las cajas de herramientas de matlab. Ellos usan FFT allí. cuando realizo la misma operación en los mismos datos obtengo resultados diferentes a los de matlab. Solo echa un vistazo: MATLAB : Msig = 0 0 0 0 0 0 0 0 0 0 0 0 0 […]
¿Cómo puede determinarse el ángulo de rotación mediante la correlación de fase (usando fft) de 2 imágenes? El algoritmo dado en http://en.wikipedia.org/wiki/Phase_correlation devuelve un cambio lineal, no angular. También menciona que las imágenes se deben convertir a coordenadas log-polares para calcular la rotación. ¿Cómo se logra esta conversión en python? ¿Y después de la conversión […]
¿Hay alguna manera de cambiar la fase de un archivo wav en python? Estoy tratando de lograr la reducción activa de ruido. Lo que planeo hacer es grabar el ruido ambiental y luego cambiar su fase en 180 grados fuera de fase. Luego, grabaré otro archivo wav con alguien que esté hablando con el micrófono […]
Lo que bash es filtrar mis datos con fft. Tengo una señal ruidosa grabada con 500Hz como una matriz 1d. Mi alta frecuencia debería cortarse con 20Hz y mi baja frecuencia con 10Hz. Lo que he intentado es: fft=scipy.fft(signal) bp=fft[:] for i in range(len(bp)): if not 10<i<20: bp[i]=0 ibp=scipy.ifft(bp) Lo que obtengo ahora son números […]
¿Existe una función de correlación cruzada o convolución 2D basada en FFT integrada en scipy (u otra biblioteca popular)? Hay funciones como estas: scipy.signal.correlate2d – “el método directo implementado por convolveND será lento para datos grandes” scipy.ndimage.correlate – “La matriz está correlacionada con el kernel dado usando un cálculo exacto (es decir, no FFT)”. scipy.fftpack.convolve.convolve […]