Questions tagged [fft]

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

382 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged 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
472 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
227 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
1answer
22 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
31 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
58 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
49 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
54 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
9 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
34 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
71 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
106 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
41 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
65 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
63 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
81 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
72 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
56 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
32 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
62 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
42 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
58 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
57 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
56 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
405 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
248 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
467 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
507 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
52 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
151 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
66 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
191 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
1k 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
101 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
968 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
217 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
290 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
71 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
789 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
143 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
138 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" ...
1
vote
0answers
184 views

Understanding why spectral flatness cannot be computed

I have a real signal recorded 512Hz for 1s. After resampling the signal at 256Hz, I would like to compute its spectral flatness. To do so, I start by taking its discrete Fourier transform with ...
1
vote
0answers
218 views

What is the algorithmic latency of STFT processing?

I was reading this paper, and in section 3 I got a bit confused. It claims that the algorithmic latency is equal to synthesis frame length. If the overlap between two consecutive frames is 50%, the ...
1
vote
0answers
74 views

Compute hilbert transform from DCT- FFTW3 REDFT00

My problem is similar to this but not the same. Edit: Here is a description of Neumann boundary conditions for differential equations: Neumann-type boundary conditions means derivatives are ...
1
vote
0answers
514 views

Symmetry and mirror of inverse FFT

I am trying to do an inverse FFT to get the real signal from a derived formula which describes the Fourier transformed signal. I am doing this with FFT. I am fairly new to the theory of FFTs, so my ...
1
vote
0answers
2k views

FMCW radar FFT processing

I'm working on the processing of raw FMCW radar data. In my case the radar works in the [6.014-7.76] GHz band (I have 187 frequency points starting from 6.014 GHz and ended with 7.76GHz) and the data ...
1
vote
0answers
54 views

Increase power for particular frequencies

I am taking an FFT of two signals ( A very clean signal, with little noise, and a very noisy signal). Signal A The signal that has very little noise has the following FFT output: Signal B The ...
1
vote
0answers
86 views

Ideal channel in OFDM simulation

I am trying to simulate an OFDM based system in MATLAB where the estimated channel in the reciever is fed back to transmitter. As a proof of concept I wanted to return the ideal channel. I used the ...
1
vote
0answers
490 views

Radix-8 butterfly with Winograd and Cooley-Tukey algorithm

I saw the Winograd radix-8 kernel algorithm below, shown in the image. Comparing to the mathematical formula of Cooley-Tukey, there is a multiplication by $\cos$ and $\sin(\pi/8)$, which can't be ...