Questions tagged [dft]

The Discrete Fourier Transform (DFT) is a mapping between a finite set of discrete points in a (primal) domain (time, space) and the dual frequency domain. DFT requires an input sequence which is discrete, such as a sampling from an analogue audio signal.

93 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
6
votes
1answer
260 views

Analytical Expression for Convolution of Two 2D DFTs

I am trying to do an image analysis problem on some images, and I want to calculate some things in the frequency domain which are giving me problems. Say I have an (discrete) image $I(x,y)$ of size $...
3
votes
1answer
116 views

Use samples of Fourier transform as DFT?

Consider the LTI system given by: $H(z) = 1 - \frac{1}{2}z^{-1}+\frac{3}{4}z^{-2}$ Let $x[n] = (\frac{1}{2})^nu[n]$ be the input to the system. We want to find the output for $n = 0,1,...,N_a$, using ...
3
votes
0answers
611 views

Incorrect Frequency results when using Multiple Signal Classification (MUSIC)

I am using MUSIC (Multiple Signal Classification) to determine Direction Of Arrival (DOA) and frequency of signals impinging on an Antenna Array. I am using a function ...
3
votes
0answers
5k views

1/3 octave spectra from fft

I have got a signal in frequency domain. This is a frequency response function from software, so I can do nothing about it and have to leave it in frequency domain. I want to transfer the data to 1/3 ...
2
votes
1answer
59 views

DFT from Fourier transform

I always studied the DFT starting from his formula, but for some reasons I need to do comparison between the FT and the DFT. I found the pdf in this link very useful http://www.robots.ox.ac.uk/~sjrob/...
2
votes
0answers
41 views

Understanding convolution of Chirp z algorithm

I dont´t understand how works the convolution part of the Chirp z. I understand how the DFT is transformed \begin{align*} x(k) = \sum_{n=0}^{N-1} x(n) W_N^{kn} \end{align*} to this expresion: \...
2
votes
3answers
757 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)...
2
votes
1answer
4k views

How to calculate the power of a discrete signal? / Clarification on PSD estimates

I want to calculate the channel power $ P_\mathrm{x}$ of a given discrete and complex signal $x[n]$ (with a length of N) in a given bandwidth $B$. I'm aware, that I could probably apply a sharp ...
2
votes
0answers
505 views

discrete fourier transform circular symmetry

I was reading digital signal processing. In topic Discrete Fourier Transform, circular symmetry it is said that circular advance is obtained by shifting x(n) in clockwise direction. i.e. to obtain x(n+...
2
votes
0answers
55 views

DFT of signal with frequency-dependent time shift

I have a signal $s[n]$, and I want to compute the DFT of it with a frequency selective shift (computing all $k$ here but I really only care about $k > 0$) $$S[k] = \sum_{n=0}^{N-1} s[n-\dfrac{N}{...
2
votes
0answers
50 views

Separating out a sequence of DFTs with overlapping windows

Suppose I have a time series partitioned into $N$ equally-sized sections. Choosing some $k<N$, for each of the sections $i \in[0, N-k]$, I concatenate that section with the $k-1$ sections to the ...
2
votes
0answers
81 views

Why does filtering the same discrete dataset while choosing different sampling rates yield different results?

If I have a two-dimensional discrete dataset (one space, one time), and I create subsets of this dataset by "sampling" (although the original dataset isn't continuous) it at different rates, should I ...
2
votes
0answers
470 views

DFT sinusoid's fundamental frequency, intuitively

I am trying to understand, intuitively, the mapping of a sinusoidal signal's fundamental frequency from time domain to frequency domain using $\textrm{DFT}$ formula. Given the time domain samples $\{x[...
2
votes
0answers
7k views

Horizontal banding (flickering) due to electronic rolling shutters

It is a well-known artifact that in CMOS cameras with electronic rolling shutter, horizontal banding (flickering), i.e. brightness intensity variations, are observed when the image is recorded under ...
2
votes
0answers
250 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 ...
1
vote
0answers
101 views

How to express STFT and ISTFT as a 1d convolution and 1d deconvolution in tensorflow/keras

I'm trying to implement this paper in tensorflow and keras. At the end of section 3 it says. ...
1
vote
0answers
39 views

DFT equivalent circular convolution weight matrix with a symmetric filter of length 2K+1

$\DeclareMathOperator{\diag}{diag}$In a research paper, I read that: For optimization, the $n\times n$ weight matrix of DFT can be equivalent to circular convolution with a symmetric filter of length ...
1
vote
0answers
42 views

How to select the sign of the square root of each element of a DFT in obtaining the square root of a polynomial?

I want to find the square root of a polynomial by the following process: Compute the N-element DFT of its coefficients, maybe padded with zeros. Compute the complex square root of each of the N ...
1
vote
0answers
58 views

Equivalence of the Power Spectral Density definitions

I am trying to show the equivalence of the following Power Spectral Density definitions in Matlab: Definition 1: $$ P(\omega) = \sum_{k=-\infty}^{\infty} r(k)e^{-j\omega k} $$ Definition 2: $$ P(\...
1
vote
0answers
78 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
129 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
825 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
1answer
131 views

Spatial spectrum of EEG data: non-uniform DFT?

I want to compute the spatial spectrum of EEG data collected using a non-uniformly sampled grid of sensors (as in the figure below). One way to do this would be to interpolate the data on a ...
1
vote
1answer
289 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
0answers
69 views

Normalised magnitudes in sliding dft

I am new to signal processing. I am trying to use the output of a sliding DFT to analyse frequency peaks. Take for instance this implementation: https://github.com/bronsonp/SlidingDFT/blob/master/...
1
vote
0answers
649 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 ...
1
vote
0answers
63 views

Determining periodic pattern when there is pattern which period is multiple of former

I have periodic signal, say, it has year period, yearly pattern. I want to check if also has quarterly pattern. Event if there is no visible quarterly pattern, Fourier Transform of this signal has ...
1
vote
0answers
681 views

How does the size of a spectrogram relate to window length and frame rate?

I'm trying to understand the methodology from section 4 in this report about audio classification but I have some trouble understanding how they pre-process the data (dsp is not my area). $30$ ...
1
vote
0answers
53 views

What is going wrong with my DFT program?

Please let me ask your help on my problem. My following DFT program outputs: [Testing frequency's amplitude: $56.9708$, and its phase: $1.32018 \pi \textrm{ rad}.$] Testing frequency's amplitude ...
1
vote
0answers
163 views

Zero padded-OFDM signal?

Why do we need 2N-FFT at the receiver when we want to demodulate ZP-OFDM signal ?
1
vote
0answers
97 views

Proving DFT equation

How can I prove the equation $(3)$? I can't understand why there is a $2/N$ in $(3)$. Why he just get the second term in DFT? Consider a sinusoidal input signal of frequency $\omega$ given by ...
1
vote
0answers
471 views

How to convert a spatial frequency in a 2D-DFT into the units radians per pixel?

Let's say I have a 2D image, and I take the discrete Fourier Transform (via FFT) of that image. In the frequency domain, I get the following image: In this image, let's just assume all the spatial ...
1
vote
0answers
376 views

Discrete Fourier Transform for text analysis?

I would like to determine the number of text-blocks verifying a roughly similar pattern. I have the intuition that I could do it using some Discrete Fourier Transform (DFT)-like methodology. Example: ...
1
vote
0answers
398 views

Short VHF pulse detection from SDR

I need to detect a short (10ms) VHF pulse that occurs every 1.5s at a specific frequency (e.g., 150.20 MHz). I have an SDR that support IQ sample rates up to 10MSPS but due to computing constraints 2....
1
vote
0answers
76 views

Is negative frequency important in calculating Spectral Moments

A follow on question form this one (What is meant by "spectral moment"?) In answers to the above referenced question spectral moment is defined as: $$ m_k=\int_{-\infty}^{\infty} f^k|X(f)|^...
1
vote
0answers
38 views

Variation in Signal Power in Mimo system

I am exploring methods for reducing computations required by MUSIC (Multiple Signal Classification) Algorithm. Most of the computational power is spent on calculating SVD. My approach is to take ...
1
vote
0answers
931 views

Instantaneous power estimation by discrete hilbert transform - how far does it smooth?

In my research area, instantaneous power in a specific frequency band is commonly estimated by the following procedure: Apply a bandpass filter on the raw signal (e.g. 80-90Hz bandpass). Estimate ...
1
vote
0answers
149 views

Multiply filter kernel by sine

1)Suppose that we have low-pass filter with cf equal to 0.1. How frequency response would look like if I multiply filter kernel by sine with frequency say 0.3? How this is related to converting low-...
1
vote
0answers
488 views

Using Fourier Transform on Gyroscope

The original idea is to calculate distance from accelerometer input. However, accelerometer reading also contains the gravitational values, thus to remove gravity, I tried using Gyroscope. The idea ...
1
vote
0answers
692 views

Extraction of fundamental signal information-Fourier full cycle algorithm

After filtering my noisy input signal using an anti-aliasing and FIR filter, I now wish to get the basic signal information (peak voltage and impedance; $R$ and $X$) from the pre-filtered as well as ...
1
vote
0answers
1k views

DFT, cepstrum, LPC for feature extraction

I am a newbie in speech processing and experimenting to get a feel. I have extracted some speech segments using a window function and I want to find distance between a pair of segments. First, I took ...
1
vote
1answer
783 views

FFT analysis of longer data sets

I have a problem with a FFT analysis in MATLAB, which is probably related to my limited understanding of the fundamentals of Fourier analysis. I use MATLAB’s built-in fft analysis function. The ...
0
votes
1answer
56 views

Anyone explain to me this video?

I was watching a video in time 24:48 I would like to know where you got the value .9 (1.14z + .941) and 1.0232 + .757 Does anyone explain how he got those numbers?
0
votes
0answers
38 views

Where to apply zero-padding for convolution dealiasing and appropriate scale

After hours of browsing the DSP posts and resources online, I still struggle to understand why my code diverges when I activate zero-padding dealiasing. When I deactivate it, everything 'works well', ...
0
votes
0answers
17 views

DFT of a function and array convolution

I saw some questions (and answers) on this subject, but they were all about a specific example and I'm not sure I understood. I'm trying to understand the meaning of computing the DFT of an array ...
0
votes
0answers
58 views

Convolution of a non-symmetrical window function by a cosine signal in the frequency domain

A have a time signal: The associated DFT spectrum of this signal: The time signal can be considered as a non-symmetrical rectangular window function multiplied by a cosine signal with a frequency $...
0
votes
1answer
27 views

Bin sizes for non-uniform discrete Fourier transforms

For a non-uniform discrete Fourier transforms, do the specified frequencies – i.e., $f_k$ in – refer to the midpoint of the bin or the lower bound? I read the answer here, but that stated that ...
0
votes
1answer
22 views

What is the form of the spectral derivative in the all-positive-frequency notation in DFT?

The Discrete Fourier Transform (DFT) of a function $u:[0,2\pi] \to \mathbb R$ sampled over $N$ equidistant points $\theta_j = 2\pi j/N,\, j = 0, \dots, N-1,$ is defined by $$ \tilde U_k = \frac1N \...
0
votes
0answers
35 views

Find repeating sequences in audio file

I have a long audio file (12+ hours). I know that there are some unknown small (2 minutes each or thereabout, it varies) repeating (not bit by bit: it is recorded with a mic) chunks in it, repeating ...
0
votes
0answers
28 views

Do twiddle factors of the fixed-point DFT have to be scaled to the input signal?

I have sensor data that is digitized by a 12-bit ADC that has analog range -1.65 to 1.65 V. The digital sensor data has thus scaling $3.3\cdot 2^{-12}$. I need to perform fixed-point DFT using ...