Stack Exchange Network

Stack Exchange network consists of 174 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [fft]

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

1
vote
1answer
80 views

Calculating THD where peaks are not exact harmonics

I have captured an FFT of a 1 kHz signal that I am generating using a WM8904 codec running at 48k samples/second connected via I2S to a ARM microcontroller (Atmel SAME70). It appears to have some ...
-1
votes
1answer
102 views

Spectral Subtraction python implementation error

I am trying to implement Spectral subtraction in python. I am new to DSP. I used this article as reference and below code is what I can come up with. Of course I did not get the expected output. I ...
0
votes
0answers
52 views

Comparing Octave/Matlab FFT results against CPU FFT

I wrote a program which runs FFT from Intel MKL and am trying to compare its results with that of Octave, but things don't quite match. The program does: Populates a 1D array of, say, 1024 elements, ...
1
vote
1answer
376 views

Amplitude of frequency bin in FFT doesn't match time-domain amplitude

I have a sinusoidal current that I am sampling at about 357k SPS. The current signal is about 3A pk-pk @ 750hz (top figure). I am wondering why when I take the fft using numpy/scipy's fft function, ...
0
votes
0answers
22 views

Is there any relationship between the frequency band and row number of spectrogram [duplicate]

Can we derive some equation or relationship between the frequency band of the spectrogram and the row number of the spectrogram matrix. For example, if we have a spectrogram ...
1
vote
3answers
308 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)...
0
votes
2answers
109 views

Software to extract frequency spectrum data from an audio file

I've been trying to learn how computers recognize musical notes, and I found apps like spectrogram which shows the frequency spectrum of a sound, but what I need is the data which is plotted or the ...
1
vote
3answers
183 views

Denoising a signal

I'm starting hydraulic experiments, where I'd have to measure velocity in an unsteady flow with a device called Acoustic Doppler Velocimeter. In DSP terms, I'd have a nonstationary signal in a shape ...
2
votes
1answer
66 views

Small shift in frequency in FFT

I'm trying to compare the estimation of the power spectral density (in this case of a fake cosine signal) by two different methods: 1st, MATLAB's fft function, and 2nd, a direct calculation from the ...
0
votes
1answer
58 views

methods of crest factor minimization

i am studying the methodes CF minimization. and have a doubt. in each methods there are steps with DFT and IDFT or like this FFT and IFFT. Everyone provides the fact, that they use DFT and IDFT for ...
0
votes
1answer
165 views

How can I show DTFT result in MATLAB? [closed]

I wanna show DTFT result and convolution result in f-domain are same. Additionally, I wanna show sampled function and inverse fourier transform result are same too. How can I show this? If I use ...
0
votes
1answer
42 views

Why are these results different about 'frequency shifting'?

I think 'spec_ori' & 'spec_shifted' spectrum amplitude and phase are same in below code. Because of 'spectrum shiting property' (What effect does a delay in the time domain have in the frequency ...
0
votes
0answers
31 views

How to calculate FFT when many samples are observed?

I was given a classwork(which I flunked) where a continuous time wave of the format below was given: $x(t) = A\sin(\omega_1 t) + A\sin(\omega_2 t)$. with different frequencies for $\omega_1$ and $\...
1
vote
1answer
300 views

Frequency analysis to determine low-pass cut off frequency

I collected some data for a practical application, where the signal represents force data obtained from an impact of a punch against a force plate attached to a quasi-rigid rig (it moves once the ...
0
votes
0answers
115 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
0answers
80 views

How to know whether numerical roundoff error or not in MATLAB?

This is result of real and imaginary parts of cos fft. I think e-13 number is numerical error but, I don't know how to know they are real numerical error. If I have to try change primitive types, ...
0
votes
0answers
35 views

Confused about obtaining impulse response from overlapped FFT'd signal

I have two audio signals which I have performed an FFT to plot the magnitude vs frequency response in realtime. To do this I have divided the incoming signal into small blocks of size 2048, ...
1
vote
1answer
283 views

Why is the second half of the FFT negative frequencies

There are lots of questions here about what the negative frequencies in fft mean, but I'm confused on why the second half of the fft calculation is the correct calculation for a negative frequency. ...
1
vote
1answer
203 views

Why MATLAB fft cos makes imaginary parts?

The cos fourier transform has no imaginary parts, but in this code it has imaginary parts that little big. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% figure(7); fs=128; t=-4*pi:1/fs:4*pi-1/fs;...
0
votes
0answers
118 views

FFT amplitude vs IQ amplitude

I'm very new to digital signal processing, so please let me know if my question doesn't make sense. I calculated the amplitudes (using the formula : $10 \log_{10}(i^2 + q^2)$) for some IQ data ...
0
votes
2answers
218 views

Ifft through Matrix multiplication

I am still new to MATLAB, so apologies if I sound lazy to you. I am attempting to model a transformation as a set of matrix operations. I start with a vector, up-sample it by $U$ (up-sampling rate), ...
1
vote
1answer
125 views

Understand Inverse DFT (FFT)

I was reading this article to understand the fft https://jackschaedler.github.io/circles-sines-signals/dft_walkthrough.html I understand that the fft takes a block (256 or 512) of the signal, ...
1
vote
2answers
84 views

Radix2 fft - zero padding output size

With zeros padding the FFT radix2 has different input size than output. How to deal with that? I see two solutions but not sure which one is better/appropriate. Maybe you know some other and better ...
0
votes
1answer
443 views

Numpy fft2 frequency shift

I don't understand how to make frequency shift in fft2 or higher dimensions. Could anyone explain it, please? import numpy as np import matplotlib.pyplot as plt from mpl_toolkits.mplot3d import ...
1
vote
3answers
151 views

Deconvolution Using Complex Division in The Frequency Domain

Consider these two signals: a = [1 1 0 0 0 0 0 0] b = [1 0 1 0 0 0 0 0] their convolution is c = a * b = [1 1 1 1 0 0 0 0] ...
0
votes
1answer
45 views

Meaning of a formula about transient frequency bins

I'm trying to understand the following formula in this paper: Sadly, the paper does not define most of the symbols it is using, so my guesses so far are: $T_{tr}$ refers to the sampling interval i....
0
votes
1answer
39 views

Complex FFT reversing

I've succeeded to compute complex FFT of a data array of interleaved values using the arm_cfft_f32() API of the CMSIS DSP library for Cortex-M4, and got in the same array the FFT results as frequency ...
0
votes
1answer
23 views

Non-decaying Harmonics in Vibration Signal

I am analyzing vibration data collected by an accelerometer positioned in a stamp press machine. The fundamental frequency is around 10 Hz. The problem here is that when I take the FFT of the data I ...
0
votes
0answers
19 views

What means transition band in bandpass filters? [duplicate]

im no a specialist in dsp, but i would like to know what means a transition band in band-pass filters, in this statement: ...
0
votes
1answer
38 views

FFT for distance-limited sum of pairwise products

I have a sequence of fixed length $N$ input arrays ($x_0$, $x_1$, ..., $x_{n-1}$), and my goal is to compute output arrays as such: $$d_{00}[l] = \sum_{i = 0}^{N - l} x_0[i]x_0[i + l]$$ $$d_{01}[l] = ...
0
votes
2answers
49 views

What to do for the right FFT in Matlab (two peaks and incorrect amplitude)?

This is the FFT I'm using. But two peaks (did not expect the one in the beginning) are occurring and the amplitude is not as specified (expecting a value of 2). Any help is much appreciated! ...
0
votes
1answer
44 views

Harmonics in Quantal Histogram

I'm studying what I believe to be a system that produces a (noisy) quantal response (i.e. there should be a minimum increment of response). To do this, I've plotted a histogram of some ~6000 amplitude ...
0
votes
1answer
92 views

THD: estimation using periodogram or fft?

I want to estimate the THD of a voltage signal (output from a frequency inverter). At the begining, I need to estimate the power of the signal at frequency domain. To doing so, I should use a ...
0
votes
2answers
125 views

How to assess perceptual volume from FFT

So far I could estimate signal volume from FFT using amplitude information, then adding all amplitudes together and divide this number by the number of frequencies. But I think this method does not ...
0
votes
1answer
405 views

Checking for VHF pulse with SDR and python

So I have a small NooElec Software Defined Radio and I've got it working well with my Raspberry Pi using GQRX as a controller and to visualize the SDR output. My goal though is to use a python script ...
1
vote
4answers
489 views

Why does the frequency in the DFT have to be an integer?

I don't know why the result of DFT and FFT in MATLAB like below images.. fs=128; t=0:1/fs:1-1/fs; x=cos(2*pi*3.5*t); X=fft(x); N=length(x); n=0:N-1; f=fs*n/N; plot(f,abs(X)/N); If I set the ...
0
votes
0answers
67 views

Making longer a time domain signal by adding values on its frequency domain. What am I doing wrong?

Background Here's the thing: using software for Finite Element Acoustic Simulation I got the dataset of frequency response from a room; software works by solving the wave equation in the interval $`[...
0
votes
0answers
142 views

Split Radix Cooley Tukey Code implementation

I've implemented the Radix-2 cooley tukey algorithm from wikipedia in java & c and I've also managed to do a Radix-4 Variant of this algorithm. I'd like to implement a split radix variant of the ...
0
votes
0answers
72 views

Simple tests for normalisation of power spectra

I am trying to compute power spectra of multidimensional (2d, 3d) arrays. In particular, I want to compute radially averaged power spectra $P(k) = P(|\mathbf{k}|)$ using FFTW, where $P(\mathbf{k}) = |\...
0
votes
0answers
97 views

Normalize FFT signal (DC part)

I'm struggling for few hours now with FFT normalization. I'm working in Matlab to prove my concept before moving to embedded application. I have 4 sine wave signals with: ...
0
votes
0answers
44 views

calculation of derivative using FFT

Here is my issue: I am trying an algorithm using FFT to fit and get a function's approximative in order and calculate it's derivative in an more easy way ( for any data), but I get some discontinuity....
0
votes
1answer
100 views

How to interpret the signal in frequency domain

I'm having difficulties understanding some basics regarding FFT. I was hoping you could answer a few questions I have. These questions will probably seem silly to you, sorry for that. I'm a complete ...
0
votes
0answers
79 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 ...
4
votes
3answers
1k views

Correct method for processing FFT in realtime

I am attempting to build a 'live' FFT for continuous analysis of an audio signal. I am relatively new to this, but have attempted the full thing, and then updating the graph at 20Hz refresh rate. My ...
0
votes
2answers
46 views

Why the angle of a real signal in frequency is continuous?

This original signal is shown below (time domain): It is a real signal. According to the signal processing theory, the form part of the angle is conjugate to that of the latter part. That is, $\...
3
votes
2answers
120 views

Auto-correlation function, an inverse problem

$x[n]$ is a complex function $n=0,1,2,\cdots,L-1 $ we assume $x[n]$ is periodic in its index: $x[n+L]=x[n]$ Its auto-correlation function $C[n]$ is uniquely defined as: $$ C[n]=\sum_{i=0}^{L-1} x[i+...
0
votes
0answers
57 views

non-negative FFT spectrum, subtlety

definations $x(t)$ is a function of time. Physically, it can be voltage, displacement, magnetization and so on. It can be real, complex, vectors or more fancy numbers. $C(t)$ is the auto correlation ...
0
votes
0answers
41 views

How to properly filter an image in the frequency domain using matlab

my problem is the next: I´m starting to understand the world of image processing and I´ve been kind of bugged with the idea of filtering images in frequency domain. I have made a test programm in ...
0
votes
2answers
67 views

How to detect a long signal inside an even longer one?

I'm trying to detect a long pattern inside a long recorded audio. In order to do so I am implementing a matched filter in time by performing a FIR filtering using the coefficients of my pattern ...
0
votes
1answer
252 views

extracting ultra accurate phase and frequency with FFT with real world data

I was doing some experiments in Matlab, and write this code to experiment with it. I want to know if it's possible to extract accurate within -+0.01° phase info from samples of an ADC, the matlab FFT ...