Articles of fft

Alineación de memoria para FFT rápida en Python usando arreglos compartidos

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 […]

Análisis del espectro de Python

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 […]

Convolución circular 2D Vs convolución FFT

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 = […]

scipy / numpy FFT en datos del archivo

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 […]

Recrear datos de series de tiempo usando resultados FFT sin usar ifft

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, […]

Teorema de Parseval en Python

Estoy tratando de controlar la funcionalidad fft de Python, y una de las cosas extrañas con las que me he topado es que el teorema de Parseval no parece aplicarse, ya que da una diferencia de aproximadamente 50 ahora, mientras que debería serlo. 0. import numpy as np import matplotlib.pyplot as plt import scipy.fftpack as […]

FFT cambio de frecuencia inesperado después de la aplicación de la función de ventana

Obtuve este código de Python para el cálculo FFT de una señal de sonido: from math import * from cmath import exp, pi def fft(x): N = len(x) if N 0: print 20 * log10(sqr) # …in decibels else: print “-INF” Obtuve los siguientes resultados: SIN función de ventana (comentada): -20.3017238269 -16.9192604108 -12.5089302395 -8.97999530657 -5.96033201086 […]

La FFT inversa devuelve valores negativos cuando no debería

Tengo varios puntos (coordenadas x, y, z) en un cuadro 3D con masas asociadas. Quiero dibujar un histogtwig de la densidad de masa que se encuentra en las esferas de un radio R dado. He escrito un código que, siempre que no haya cometido ningún error que pueda tener, funciona de la siguiente manera: Mi […]

El filtro en el espacio de Fourier no se comporta como se supone que debe

Este es un seguimiento de una pregunta respondida que hice y que se puede encontrar aquí . Tengo varios puntos (coordenadas x, y, z) en un cuadro 3D con masas asociadas. Quiero dibujar un histogtwig de la densidad de masa que se encuentra en las esferas de un radio R dado. La idea es calcular […]

¿Cómo creo una matriz numpy de cadena?

Tengo un lector de archivos que lee n bytes de un archivo y devuelve una cadena de caracteres que representan esos datos (binarios). Quiero leer n bytes en una gran cantidad de números y ejecutar una FFT en ella, pero tengo problemas para crear una matriz a partir de una cadena. Un par de líneas […]