Episode #125 of the Stack Overflow podcast is here. We talk Tilde Club and mechanical keyboards. Listen now

Questions tagged [fft]

The fast Fourier transform is an efficient algorithm to compute the discrete Fourier transform (DFT) and its inverse.

400 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
2
votes
0answers
248 views

Is there any problem with my FFT method?

I am trying to analysis a system by Fourier Transform. my system has 3 dominant modes. (1,0.05) (0.51,0.01),(0.46,0.01) the pairs are defined as (freq,damping_ratio) I've got those result with ...
2
votes
0answers
111 views

Complex Cepstral Transform and Audio Watermarking

I am doing research in audio watermarking and have come across the paper that uses the complex transform to embed the watermark in order to improve the robustness and inaudibility of the embedded ...
2
votes
0answers
483 views

How to find the phase lag between a sinusoidal input force and the response?

I have a system that is forced by a sinusoidal function with specific amplitude, phase and frequency. I get a deformation of sinusoidal form as well, and I would like to find the phase lag. I use FFT ...
2
votes
1answer
254 views

Unexpected peaks in spectral analysis

I am doing the spectral analysis of some antenna recorded signals by computing PSD (logarithmic absolute value of the FFT) and the spectrogram of the signal. The problem is that I get some unexpected ...
1
vote
3answers
26 views

FFT error when number of samples growth - Matlab/octave

I see when I increase number of samples, the imaginary error part of the FFT growth due to the error of sin(n*pi)!=0 where n is the index of the sample in Matlab or octave. Is there anyway to force ...
1
vote
0answers
39 views

Why are low frequency peaks not sharp in scipy fft?

I am using numpy/scipy to plot graphs of sine waves. frequencies at 15 Hz or higher give nice, sharp peaks, but at lower frequencies, the peaks are smeared and the actual peak frequency can't be seen ...
1
vote
2answers
72 views

Multi Resolution STFT for audio

I am looking at trying to achieve a multi resolution STFT analysis on a log sine sweep audio measurement. I have an impulse, and I would like to perform two different windowed length FFTs on the ...
1
vote
0answers
39 views

Scaling issue after doing manipulations in FFT

I have two filters which must apply (in series) to same input data $x[n]$ e.g of length 256 (real time asio input buffers). One filter is e.g. an impulse response $h_1[n]$ of a speaker and other one ...
1
vote
0answers
76 views

SNR computation in frequency domain with scipy.fftpack.fft

I have a noisy time series (gaussian coloured noise) to which I add a signal. I call the sum of noise+signal a segment. In the picture you can see an instance: I want to compute the signal-to-noise ...
1
vote
0answers
29 views

how to cut off frequencies above 10 Hz with Welch method

This is the first time I have worked with any code so I am rather unfamiliar with all the lingo and processes involved. I am trying to convert a signal to the frequency domain and want to find the ...
1
vote
1answer
35 views

Measuring group delay through a linear phase filter

I'm having a difficult time wrapping my head around a phase/ group delay problem and it corresponding to where you start sampling so I built up a system to test some of these issues. The system passes ...
1
vote
1answer
28 views

Is fft2 in MATLAB unitary? Some differences happen

I meet a problem when implementing fft2 in MATLAB. The question is I try to simulate the realistic measurements $Y = |FCXF^H|^2$ - the intensity of Fourier domain of object $X$, where $F$ denotes ...
1
vote
0answers
34 views

When plotting the cosines of the phases of some sensors, I got an interesting plot

I have a laser that I modulate at some frequency, typically 1250 Hz. I have a sensor that tracks what happens inside a reactor when it's illuminated. I am doing Fourier analysis on both the laser ...
1
vote
2answers
110 views

Interpolation of estimated frequency between fft bins: what about phase?

I am developing an additive re-synthesier, which involves identifying partials from FFT frames. Partial identification requires amplitude peaks. The simplest way to do this is to just look at the ...
1
vote
0answers
64 views

Additive Synthesis - clean playback of a spectrum frame

First of all - I have only a limited of Fourier transforms and the mathematics involved, so please pardon any potentially bad terminology or huge gaps in conceptual knowledge. I'm experimenting with ...
1
vote
0answers
91 views

FFT of accelerometer signal: How to get a correct amplitude & noise reduction

I’m a student studying on vibration analysis & testing. I’ve performed low level sine sweep test of a simple rectangular PCB with clamped boundary conditions on two ends of it. That test was ...
1
vote
0answers
11 views

Complex Spectral Phase Evolution (CSPE) Performance depending on signal windowing?

I am look into CSPE. "Signal Analysis Using the Complex Spectral Phase Evolution (CSPE) Method" The method is simple. It compares the original signal's FFT and shifted signal FFT in phase domain so ...
1
vote
1answer
42 views

How to achieve a smoothed frequency phase response with no outliers

I have two sampled signals, which I then perform wavelet analysis on, by synthesising the Filterbank, perform the dft on each kernel, and the dft on the sampled signals, and multiply together. ...
1
vote
4answers
107 views

Detecting a three-tone signal in real-time

I'm new to digital signal processing, and I've been looking into FFT, filtering, Matlab etc. What I'd like is to build a program which will detect when a certain BEEP is transmitted through a trunked ...
1
vote
1answer
146 views

FFT Heart Rate Monitor on Arduino

Part of my project is to build heart rate monitor from Nano and a PulseSensor. I am having trouble working with the Fast Fourier Transform part of the code. This code here is from a library (...
1
vote
0answers
56 views

Zero-padding vs. nonzero-padding in computation of auto-correlation with FFT

Isn't the usual zero-padding in the computation of the auto-correlation function with FFT just one of many possible extrapolations of the original signal? If I have a measured signal which has good ...
1
vote
2answers
88 views

Multiple different length Windows for high resolution realtime FFT

I am looking at trying to achieve a 'realtime' (as quick as possible to acquire, process and present data) FFT application to analyse audio. I have setup my app to acquire a set of samples, apply a ...
1
vote
2answers
92 views

Checking Parseval's Theorem for Gaussian Signal by Using Scipy

I'm trying to check Parseval's theorm for Gaussian signal. It's well known that fourier transform of $\exp(-t^2)$ is $\sqrt{\pi}\exp(-\pi^2 k^2)$. So I implement it by using quad and simps. I think ...
1
vote
1answer
133 views

Best way to resample an OFDM signal without interpolating + filtering?

I am working on a custom transceiver with a sampling rate of 3.84 MHz (256 OFDM subcarriers at 15 kHz, although 56 subcarriers are null). I would like to resample the signal to 5 MHz, in order to use ...
1
vote
0answers
74 views

Applying DFT to OFDM signals above Nyquist rate

Consider an OFDM signal which has $N=8$ subcarriers and is band-limited to $[-B/2,B/2]$, $\sum_k x_k e^{j2\pi kf_0t}$ where $f_0=\frac{B}{N}$. For this case the Nyquist rate is $2\times B/2=B$, i.e., ...
1
vote
0answers
62 views

FIR Audio Equalizer Crackling Noise

I've implemented a few FIR filters into a song equalizer, but there seems to be a crackling that I'm unsure how to fix. Heres a video for reference: https://www.youtube.com/watch?v=IlzWRYpv0rk. Noise ...
1
vote
0answers
33 views

Determine 3 most dominant frequencies in a signal, PSD or just the absolute value of a Fourier transform?

I have a noisy ECG signal recorded for 5 minutes. My goal is to determine heart rate every 2 seconds. To find out 3 most dominant heart rates (Beats per minute) in a 2s-signal should I calculate its ...
1
vote
0answers
100 views

How to sample a transfer function (angular spectrum) in the frequency domain?

I am having this problem in Fourier optics, where I am using the Angular spectrum method (a lti filter) to calculate the electric field at the required plane given ...
1
vote
0answers
48 views

Frequency bins too powerful at the start and at the end of each sample computed DFT data

I've taken the following link http://coding-geek.com/how-shazam-works/#comments as a project to start practising and learning signal processing and python programming for me to get started in speech ...
1
vote
0answers
93 views

dsPIC33E ADC to FFT help

I am sampling an audio range signal with a bandwidth of 3100Hz and then applying a FFT using the DSP library example from Microchip to determine the most dominant frequency of the signal. On the last ...
1
vote
1answer
73 views

Break a signal into segments, add them together and then perform an FFT on the result

If you break a 4 second signal into 0.5 second segments, add them together and then perform an FFT on the resultant 0.5 second data. Will the FFT / frequency spectrum be the simlar to an FFT on the ...
1
vote
0answers
72 views

Can Reassigned-STFT be an invertible transform? Algorithm for “Signal => Reassigned-STFT => masking => Reassigned-ISTFT => Output signal”

It is well known in DSP that this is a successful method to modify audio signals: Signal => STFT => masking => ISTFT => Output signal For example, ...
1
vote
0answers
17 views

artifact identification in irregular-samped data

There is a dataset combined by different datasets a,b,c,d,e. The time intervals are irregular(i1, i2, i3:intervals between a and b, b and c...), but the sampling rate SR in every single dataset is the ...
1
vote
0answers
520 views

Remove a line in image with FFT

I would like to remove vertical line from an image (an example). I took a 2D FFT and try to apply a mask to suppress the line. Nonethelesse the approch is not very efficient, because i lose an ...
1
vote
1answer
297 views

Image zooming with Fourier transform

I want to go from this image into this one: So basically I need to scale the white square. The authors of the paper claim that this can be done in four steps: zero-padding in real space (image is ...
1
vote
0answers
674 views

Acceleration time series to velocity conversion using DFTs

I have 1000 Hz time series data for acceleration (512 data points), which I want to convert to velocity. I am trying to use the omega arithmetic method to achieve this. Following are the steps I am ...
1
vote
3answers
638 views

How do I calculate peak amplitude of the signal components after zero padding and FFT?

I am learning about DFT and trying to apply it to some audio processing. I am new to DSP but experienced in programming and have some background in math and physics. The FFT algorithm I use (lomontFFT)...
1
vote
0answers
53 views

which is the relationship between domains (x and k) in FFT?

I need to plot an escalar field f(x,y) in an interval xv= (xmin,xmax) in Nx steps and ...
1
vote
0answers
170 views

Frequency Domain Variant of Pearson's Sample Correlation Coefficient?

Suppose that you need to find (i.e. acquire) $x_1$, a short data packet with a known bit sequence preamble (i.e. a PRBS BPSK), in $x_2$, quantized samples recorded from the output of an ADC. You don'...
1
vote
0answers
86 views

1D as a 2D FFT - am I understand it properly

in one of my question about non $ 2^L $ points FFT I got answer with advice to read about that in following book: Rabiner, Lawrence R., and Bernard Gold. "Theory and application of digital signal ...
1
vote
4answers
211 views

How to see the amplitude of each frequency?

do you know software that exports spectrograms to images of 500x22000 pixels? all the spectrograms that I know show diffused frequencies and it does not look very good I just want to see the amplitude ...
1
vote
0answers
2k views

USRP + GNU Radio as 1 GHz bandwidth spectrum analyzer

I have maybe not trivial but some interesting question. I need to create 802.11ac channels spectrum analyzer with USRP. Hardware: USRP B210 Software: GNU Radio 3.7.9 PC: Intel Core i3 i3210m 16 GB ...
1
vote
0answers
110 views

When Is a Power of 2 FFT Slower than Smaller Sized Non Power of 2 FFT in MATLAB

Knowing that computing an FFT is faster if the amount of samples is a power of 2 I have always tried to pad the inputs to Matlab's FFT with zeros until the next power of 2 is achieved. Matlab's ...
1
vote
2answers
1k views

Extracting peak frequencies using FFT vs. time-domain peak finding

I am trying to compare two processing algorithms that I believe should be producing very similar results. I am processing acceleration signals over the course of the day and we expect to see a ...
1
vote
1answer
253 views

How to Combine 8 N/8 FFT's into one N FFT

I need to make in FPGA (using Verilog) an FFT. Input data is N=8192 points at 1 GSPS. However, the FPGA operates at 125 MHz, therefore the data is split into 8 channels (each one at 125 MHz). This ...
1
vote
3answers
351 views

Why is pitch detection necessary in pitch correction?

I am trying to understand the process of changing the pitch of a signal. After performing the FFT on a windowed signal, we have the frequency information for the signal. Why can we not simply change ...
1
vote
0answers
76 views

Memory requirements for FFT on STM32F103C8

I have a limited system and would like to implement FFT within STM32F103C8 without any extra memory buffers. So I want to know how many memories are needed to implement if I have 2592x1944x8bit size'...
1
vote
0answers
879 views

Vibration Acceleration Analysis with FFT and inverse FFT

I've been reading time series data off of an accelerometer, and converting it to one RMS acceleration in g. The sensor works in a range of 0-6.4 kHz approximately. Ideally i would like to restrict the ...
1
vote
0answers
155 views

Speaker equalization for white noise

I need a speaker to play broadband white noise for my scientific application. I would like to know whether my relative naive implementation is correct or not. The speakers I use have a large frequency ...
1
vote
0answers
160 views

Imfilter, FFT-bsed convolution speed

I want to convolve two large matrices and it is working fast - as long as one of those matrixes contains only separate points (point sources). If I replace the points by eg. Gaussians (line "HERE" ...