Questions tagged [fft]

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

Filter by
Sorted by
Tagged with
0
votes
0answers
14 views

How to relate power spectral density with the number of resource blocks in LTE?

I am running a simulation with 100 Resource blocks for LTE OFDM and I plotted a PSD for it. Now I am able to bring down the out of band emissions when I perform an algorithm on it. Based on the ...
3
votes
2answers
613 views

Understanding a simple case of fft2 in python

I have 2D arrays representing surfaces. I have created them by scanning some theoretically flat objects and I ended up with what can be called "imperfections", or else, deviation from the perfectly ...
0
votes
0answers
22 views

Strange phase for gaussian beam 2D

I have a 2D Gaussian beam : After doing FFT2 and angle I get strange results: This is the code I'm using (Python 3) ...
0
votes
1answer
935 views

How to plot Full spectrum with negative frequencies

I'm trying to do frequency analysis on signal of a rotating shaft. The signal is generated in two measurement planes i,e X-position and Y-Position. I was able to use fft and transform into frequency ...
0
votes
1answer
40 views

Phase differences between DFT and FFT (in MatLab)

I have implemented a (direct) DFT in MatLab (following script) and compared it to the built-in FFT routine. The magnitude response seems to be identical (excluding some possibly round-off errors), but ...
2
votes
1answer
110 views

Purpose of multiplying fft output by Hann window?

I'm working with some acoustic propagation models at the moment. A feature I've noticed in 2 of these models is that the output of the model is often multiplied by a Hann window in the frequency ...
2
votes
1answer
7k views

Simple fft to Gaussian pulse with MATLAB

trying simply to create femtosecond pulse in MATLAB exactly like in the attached image, the carrier frequency is around 374THz, and my sampling frequency is 10 times the carrier. The results of the ...
0
votes
1answer
42 views

DTFT frequency range

$$X(e^{j\omega}) = \sum_{n=-\infty}^\infty x[n] e^{-j\omega n} $$ The frequency term $\omega$ in DTFT is normalized as $\omega = \frac{\Omega}{f_\mathrm{s}}$ $\Omega= 2 \pi f$ is the angular ...
0
votes
0answers
56 views

Restoration problem with numpy fft function and self-made high-pass filter

Hi, I am using numpy.fft functions to do CV assignments in univ. But while implementing the high-pass and low-pass filter, I found something wrong during the process. The difference should converge ...
4
votes
2answers
3k views

Removing drift from integration of accelerometer data

I am trying to get a positional data from the accelerometer data using the following steps: Re-zero the accelerometer value Removing mean from accelerometer value First integration of accelerometer ...
2
votes
1answer
40 views

Electret microphone RC filter design contradiction among information sources

I'm building a circuit for an electret microphone and I want to build a bandpass filter around the op-amp. I was using this source until I have found that most sources ( e.g. here and here) indicate ...
1
vote
2answers
438 views

Color Artifacts in Fourier Transformed Image

I am working with images and Fourier transforms. I am trying to understand what might be causing some artifacts in my output image. I am starting with a 512x512, RGB image of Lenna. I FFT the image,...
2
votes
2answers
238 views

Verifying Channel Response Identity

If $x[n]$ is input signal, $y[m]$ is output signal, $h[m]$ is impulse response of the channel, then the following relation holds $\lvert Y(\omega)\rvert = \lvert H(\omega)\rvert^2 . \lvert X(\omega)\...
0
votes
1answer
635 views

2D FFT of Image and Sharpness

I'm trying to make some simply autofocus algorithm for my device and pretty confused with 2D FFT results. As I know for now - number of high-frequency's components in 2D FFT will say me that image is ...
0
votes
0answers
125 views

Correctly scaling Power Spectrum and Correlation Function of Galaxies from Mock Catalogues

I'm currently writing a code to compute the two-point correlation function of a set of galaxies. As I work with a loooot of galaxies, it was suggested to me not to do it in real space, but to compute ...
2
votes
1answer
117 views

Why Does 2D FFT of Gaussian Looks More Sharper than Gaussian Itself?

I am trying to understand why 2D FFT is done on a Gaussian process in a particular code. From my understanding from these posts: https://www.researchgate.net/post/...
1
vote
1answer
69 views

Amplitude and Phase of FFT gives messy results - Matlab [closed]

I'm trying to calculate the spectrum magnitude and phase of a signal, but when I plot them I get pretty messy results: The signal is displayed on the right, then I have the magnitude and phase of ...
2
votes
1answer
1k views

Filter - Spatial Domain Versus Frequency Domain

Note: This is a followup question of How can I create a custom band-pass filter?. Can you explain the following python source code?. Bandpass filter bank. Suppose, I have an equation ...
0
votes
0answers
19 views

Wavelet Transform Implementation Method

I have been watching the wonderful youtube videos here: https://www.youtube.com/watch?v=4TTpwIZrUAo&list=PLn0OLiymPak2G__qvavn3T8k7R8ssKxVr&index=3 detailing morlet wavelets and fit ...
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 ...
0
votes
2answers
71 views

What should my reference value be when converting FFT bin amplitudes to dB?

I want to transform my FFT output values into a dB scale, but I'm struggling to determine the function I should run each bin amplitude through. My understanding of the decibel scale is that a value ...
4
votes
2answers
2k views

Implementing Convolution in Frequency Domain?

Suppose, we have a bitmap image represented as a 2D integer array, int [,] image2D; whose FFT is Complex[,] fftImage2D; ...
4
votes
1answer
6k views

Difference Between Linear Convolution and Circular Convolution

If I understood correctly (and this page should confirm: http://www.cs.ioc.ee/~khoros2/linear/convolution-teo/front-page.html) if I convolve linearly (the usual point-to-point multiplication and ...
1
vote
1answer
122 views

Applying Image Filtering (Circular Convolution) in Frequency Domain

To filter an image we can: Use a 3x3, 5x5, 7x7, etc. filter, that is convolve the image and the filter in the space domain. Use a FFT on both the image and the filter, multiply them together in the ...
2
votes
2answers
1k views

Using MATLAB `xcorr()` to Find a Signal Hidden Behind Noise

I am trying to develop some code to find a transient signal in a data set I have. Before I get there however I am trying to write code using simulated data. I first am creating a signal with a ...
4
votes
1answer
47 views

Multidimensional FFT with mixed directions

Let's say I have three-dimensional complex array $X[m,n,p]$, $m \in \{0, 1, \ldots, M-1\}$, $n \in \{0, 1, \ldots, N-1\}$, $p \in \{0, 1, \ldots, P-1\}$; i.e. $X \in \mathbb C^{M \times N \times P}$. ...
0
votes
0answers
162 views

How to Extract frequency bands from Short-time Fourier Transform?

I am reproducing the results reported by the author in a research paper. In which he takes the STFT of EEG signals. The input data has a shape of (1 X500) and I used Scipy library built-in function to ...
0
votes
0answers
36 views

Zero padding frequency domain data for improving time domain resolution before inverse FFT

I am working in MATLAB with measured data from a VNA. The bandwidth is 1 GHz from 5-6 GHz with 801 samples in the frequency domain. Everywhere I have read that I need to split the data into 2 parts at ...
0
votes
0answers
36 views

Image from Inverse FFT operation is sparse

I'm trying to train a deep nerual network that takes in input, an signal in the frequency domain, and attempts to learn a mapping to another signal in the frequency domain. Basically, the input to ...
0
votes
1answer
239 views

Matrix cross correlation in python

I'm currently performing matrix cross correlation in python using : C = scipy.signal.correlate2d(A,A) where A is a 2D matrix, typically a picture. As you can ...
0
votes
1answer
79 views

FFT: do I need to include gaps in data in my FFT domain?

I am currently doing an FFT on a very large amount of vibration data, which I am going to sample with no pattern which I want to PSD eventually My question is in essence can I go from my data ...
2
votes
1answer
82 views

How to set STFT-parameters to visualize woodpecker / Sawtooth Signal? ​

I always wanted to work with FFTs and Spectrograms to characterise sounds and their frequencies using python. Lately I recorded a woodpecker at his work. Woodpecker Wav File To know where to look ...
1
vote
1answer
48 views

Frequency analysis : relation between FFT size and sampling

I am currently trying to perform an experiment given an audio signal. I am therefore sampling the audio frequency range (0 to 22kHz) with a 48k sampling rate. When performing a FFT on my signal, I ...
0
votes
1answer
63 views

What's the formula for choosing FFT sample rate and window length to make sure the frequency I want is in the center of a bin?

I don't know why I am having such a hard time with this math, but I am. I am doing phase-sensitive detection of a few sensors, which are triggered by a modulated light source. The modulation frequency ...
1
vote
1answer
40 views

Quantized error in scatter plots for a hamming windowed LTE signal

I have a time domain signal based on LTE standards. Carrier spacing = 15KHz, Symbols per slot = 6, extended Cyclic Prefix = 512 samples FFT size = 2048, No. of sub carriers = 1200, channel BW = 20MHz, ...
0
votes
1answer
44 views

How to classify the movement type of a target into two categories(direct motion or else) by using a doppler radar?

How to classify the movement type of a target just into two categories(one is direct motion which is almost straight line movement and the other is anything else) by using a doppler radar with a ...
0
votes
2answers
117 views

Understanding FFT size in OFDM modulation

I know that the FFT size is the the number of subcarriers OFDM system have. But I also know that the FFT has to be $ 2^N $ big (in order to process N samples from time domain to frequency domain). ...
2
votes
0answers
43 views

N-dimensional Cooley-Tukey FFT

I am trying to formulate an n-dimension algorithm for the Cooley-Tukey FFT. In two dimensions, it is very clear: $$ Y[k_1 + k_2 n_1] = \sum_{l_2 = 0}^{n_2} \Biggl[ \Biggl( \sum_{l_1 = 0}^{n_1 - 1} X ...
0
votes
1answer
73 views

Incorrect dominant frequency from FFT

DATA source: Bearing vibration data FEMTO IEEE PHM 2012 challenge. DATA context: Sample rate: 25.6 kHz Sample period: 0.1 Seconds (2560 points) Motor runs at 1800 rpm Steps taken: Low pass butter ...
0
votes
0answers
99 views

Spatial smoothing - 2D MUSIC [FMCW MIMO radar processing]

I'm using a 2D Music algorithm for the estimation of Range-Azimuth info in a ULA FMCW MIMO radar. The algorithm procedure is pretty simple as per Belfiori F., Application of 2D MUSIC Algorithm to ...
0
votes
3answers
56 views

In case of Complex DFT spectrum, why the x axis range from N/2 to N point mean a negative frequency?

When we transform a complex signal into frequency signal by using a complex DFT, The range from N/2 to N point on the X axis of the spectrum mean a negative frequency... But i cant understand why ...
2
votes
3answers
617 views

Generating a timeseries with an arbitrary power spectrum

I need to generate random numbers with a given PSD. To do so, I found this recipe. In case it's behind a paywall, here's now it works: Given your power spectrum $P_k$ with $0\leq k<N$ points ...
0
votes
0answers
54 views

OFDM modulator: Transmission time calculation

I have a problem with a certain task. The frequency of the signal sampling at the OFDM modulator output is 20 MHz. The length of the protection period is 2 microseconds and it is 20% of the ...
0
votes
0answers
45 views

Method for the Localisation of Impulsive Noise in Speech Signal

I'm sorry if this question has already been answered, but I couldn't find anything. I'm working on a program for noise detection and need to locate the position of impulsive noise. The following image ...
0
votes
0answers
112 views

retrieving complete 2d autocorrelation function from 2d power spectral density function (numpy.fft.ifft2)

My question is at the boundary between signal processing and python. I have a 2d power spectral density function (PSDF) constructed on Nx/2 + 1 and ...
8
votes
2answers
44k views

How to analyze audio input/spectrum correctly

I am quite new to Java and Signal Processing, but I was assigned a project dealing in audio processing. My topic given to me is a game that has a component of analyzing a song (any given song) and ...
1
vote
1answer
339 views

Two voice pronunciation comparison similarity MFCC + DTW

Calculate each MFCC to compare wave file A and wave file B, and then use FastDTW to measure the distance after two sets of MFCCs. We compared the four wave files and obtained the Euclidean distance ...
0
votes
1answer
82 views

Does dividing the magnitude spectrum of white noise by sqrt(2) give an RMS magnitude spectrum?

I understand that the RMS Amplitude of a sinusoidal signal is around 0.707 ($\frac1{\sqrt2}$) times the Peak Value, but this is not true for noise. However, an FFT of a noise signal indicates ...
0
votes
1answer
28 views

Clarification on the application of window pre-FFT

Good morning, I have a simple question about the application of a Hanning window before performing an FFT. My application is pretty textbook, i'm decoding a 4-FSK signal. Question: why does the ...
1
vote
1answer
36 views

DFT of a series of RC exponentials

Context: I'm trying use matlab to apply a single-pole filter to a time-domain ramp waveform that is generated by a sequence of time-shifted "RC steps" that are added together. The time domain voltage ...