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.

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

How to use a size N complex DFT to compute a size N/2 real DCT? [duplicate]

Can someone (maybe @hotpaw?) Explain How to use a size N complex DFT to compute a size N/2 real DCT? Also, I noticed that the DCT matrix (which can be obtained in Matlab by dctmtx(n)) contains dct ...
0
votes
1answer
45 views

How can I “smooth” square wave/pulse data to remove/reduce harmonics?

First, I have reviewed close-match questions, but can't find something that seems to match what I'm looking for. Also, it's been 30 years since college, so forgive me for forgetting stuff.. :) I am ...
0
votes
1answer
47 views

Practical question about DFT in Matlab

Please consider this piece of code: ...
1
vote
2answers
84 views

DCT vs DFT why do we need/want phase?

Why do we use complex exponentials in the Fourier transform, why do we want the phase part? As opposed to in DCT where there is no phase and only magnitude? Moreover, what does this phase concretely ...
2
votes
6answers
577 views

Clarification on defining FFT bin sizes

When reconstructing the frequency domain for an FFT, what is the most self-consistent way to do this -- i.e. how is it best to define the bin width, $\Delta f$? For example, previously I thought it a ...
2
votes
1answer
66 views

MATLAB FFT Amplitude

I am trying to understand FFT, DFT through Matlab and I am fairly inexperienced in all of these. I have a time vector t and a corresponding simple sine wave x = 0.25*sin(t). Now I would like to view ...
1
vote
1answer
12 views

Matrix form of 2D-DFT for a vectorized image

I want to apply 2d DFT to a N by N image. However, image is vectorized such that it is NxN by 1. How can i find the matrix form of 2d DFT such that resulting vector from multiplication is the ...
0
votes
1answer
33 views

Proof that DFT does not require more than N points

I'm trying to show how the discrete Fourier transform (DFT) arises from the equation for the continuous-time Fourier Transform. I've run into an interesting caveat which I can't seem to find an ...
0
votes
2answers
69 views

Have you seen these kind of DFT artifacts?

we've been working with the Spiral DFT Implementation on an FPGA and managed to get it working. Unfortunately we see these strange artifacts and have no idea where they come from. Our input is a ...
0
votes
1answer
42 views

Using fourier coefficients to reconstruct data in matlab

When doing a discrete fourier transform on some data using matlab's fft function, its output is a set of fourier coefficients but I was wondering how do I go about converting these into an and bn so I ...
4
votes
3answers
345 views

Absolute convergence of periodic sinc interpolation

An $N$-periodic complex discrete-time sequence $[x_0, \dots, x_{N-1}]$ can be resampled to an $M$-periodic sequence $[y_0, \dots, y_{M-1}]$ with $M>N$, using sinc interpolation: $$\begin{align}y_m ...
-2
votes
0answers
38 views

Can I find FFT for signals with odd number of element

I am said to divide signal into two parts even and odd numbered {index/value of n} .What if number of elements in x[n] is odd out will I divide into two networks?
0
votes
0answers
41 views

DFT truncation of signals

How can I calculate 8 point DFT of signals whose length is less then 8 ( say 2,4) then what will I assume other members in formulae "0" or rotation or keep repeating the same two numbers Consider I ...
2
votes
1answer
49 views

Goertzel derivation

Referring to this link, about Goertzel,I am confused about the final equation after N iterations Where does real = (q1 - q2 * cosine) and ...
4
votes
4answers
433 views

How to get Fourier coefficients to draw any shape using DFT?

I'm teaching myself about Fourier Series and the DFT and trying to draw a stylised $\pi$ symbol by fourier epicycles as detailed by Mathologer on youtube (from 18:39 onwards), and the excellent ...
0
votes
1answer
44 views

FFT of a signal with 0 added between samples

I'm having a rough time doing the following question from a HW. Given a sequence of your choice, what happens to the DFT if we add a 0 between every sample? I.e. if we have x[n] = [A,B,C], we turn it ...
2
votes
2answers
104 views

Why is the Fourier (or cosine) transform decorrelating?

The discrete Fourier transform (DFT) and the discrete cosine transform (DCT) both decompose a signal into its frequency-domain spectrum. One property that I have seen praised across various domains ...
2
votes
1answer
139 views

Comparing arithmetic complexity of FFT radix-2 and convolution

Let's assume we have a discrete linear time invariant system and we have a real signal $x[n]$ with length N=50 as input for the system. The impulse response $h[n]$ of the system is considered to be ...
0
votes
0answers
45 views

Discrete Fourier Transform by longhand using MATLAB - what is the lower limit for N?

I'm seeking some guidance/reassurance on my understanding so far of the DFT, as demonstrated by my MATLAB script below. I deliberately compute the DFT longhand. My first question: Are my ...
0
votes
1answer
46 views

inverse discrete fourier transfor with plain python

I am trying to calculate inverse discrete fourier transform for an array of signals. I am using the following formula: And my python code looks as follow. ...
1
vote
0answers
28 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
32 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
38 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 ...
0
votes
1answer
47 views

STFT Spectrum Analyzer: power drops over 15 kHz with songs, not with pure signals

I'm developing an audio file spectrum analyzer for a University Project. My main goal is to have an application that plots the Db Spectrum of a 16 Bit WAV PCM audio file (at this time only mono files) ...
0
votes
2answers
51 views

Resolution of Discrete Fourier Transform is 1/T - Mathematical proof?

In many articles I see that the frequency resolution of the Discrete Fourier Transform (DFT) equals Fs/N where Fs is the sampling rate and N is the total number of samples. Fs/N is equivalent to 1/T ...
1
vote
1answer
137 views

Cepstrum Calculation of Rational Function H(z)

I am trying to solve my first problems at cepstrum calculation. I want to calculate the complex cepstrum $\hat{h}[n]$ of a signal $h[n]$ with Z-Transform: $$H(z)=\frac{(1-0.5z^{-1})(1+4z^{-2})}{(1-0....
0
votes
2answers
68 views

What is the interpretation of Fourier Transform containing only imaginary part?

The FT of a unit step function is taken as: $$ X(\omega) = \int_0^\infty e^{-j\omega t}dt = \frac{-1}{jw}e^{-j\omega t} \Biggr |_{0}^{\infty} = \frac{j}{\omega} $$ The transform only has the ...
0
votes
0answers
18 views

Calculating N for DFT

We want to compute an N-point DFT of a one-second duration compact disc (CD) audio signal x[n], whose sample rate is $f_s = 44.1Khz$ with a DFT sampling of 1 Hz. (a) What is the number of necessary x[...
1
vote
1answer
37 views

How Could I extract the signal in the below Image

This is my first post so sorry if I make any mistake. The entire signal length in the screen shot below is 10 second. I ran a first order low pass filter (digital) at 10Hz cut-off frequency, and I ...
0
votes
2answers
95 views

Circular Convolution Matrix of $ {H}^{H} {H} $

We all know that Discrete Fourier Transform (DFT) corresponds to circular (not linear) convolution. That is to say, if $x(n),h(n)$ and $y(n)$ is the original signal, the filter and output signal in ...
0
votes
3answers
56 views

dft of sampled sine using python

I'm trying to write a python script to perform a 100-point DFT over a finite length sample of a sinewave at 1/8 the sampling frequency. I'm curious why my DFT magnitude plot has two spikes in it when ...
0
votes
1answer
68 views

PSD subtraction

I am interested seeing the difference between two power spectral densities (PSD) as a noise reduction exercise The blue line is the psd of my signal, and the the orange line is the psd of the ...
0
votes
0answers
13 views

Interpreting results of multiplying signal and wavelets in the frequency domain

I have followed a lot of stuff on mike cohens website and have managed to make a filterbank of Morlet wavelets and convolve with an incoming block of samples. This of course is a computationally heave ...
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 ...
0
votes
1answer
43 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 ...
1
vote
0answers
38 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(\...
0
votes
1answer
112 views

DFT of 2d real signal and Hermitian symmetry

Knowing that DFT of n-values real signal in 1d consists of n/2+1 different values where the second half of the spectrum is complex conjugate of the first one (Hermitian symmetry). However in the ...
0
votes
0answers
19 views

Criteria for Inverse STFT

So far I have studied the STFT and how it works. Following this question: Inverse Short Time Fourier Transform algorithm described in words , I got a grasp of how it works and how to implement it. ...
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 ...
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 ...
0
votes
1answer
29 views

Processing using Matlab

I'm new to Matlab.As a beginner, I want to work on few signals which are used in day today life and extract useful information through it. Can someone please help me from where I could start or share ...
4
votes
1answer
47 views

Amplitude Estimating Using a Windowed DFT

Let's say we want to estimate the amplitude A of a mono-frequent signal using a windowed DFT. The frequency of the signal is unknown, and the frequency resolution of the DFT is limited, thus it cannot ...
2
votes
2answers
58 views

Significance of modular arithmetic in DFT?

In what ways does modular arithmetic plays a part in DFT? Why is it a so integral part of DFT?
0
votes
1answer
46 views

Help me understand the stages involved in filtering a signal using Discrete Fourier Transform

I have a series of discrete values measured from a sensor. I want to filter the frequencies coming from this sequence of values. Then, if I understood the process correctly this is what I do: I ...
0
votes
1answer
49 views

How can I find the value of Inverse DFT

I have dt signal $x[n]={[6.29, 8.11,-7.46,8.26,2.64,-8.04,-4.43,0.93,-9.29]}$ And I need to give the function value of: 1) sum of $x[n] = \sum\limits_{k=0}^{N-1} X[k] \, e^{-j 3 \pi k/5}$ from k=0 ...
0
votes
1answer
39 views

Is there a name for the result of summing the bins of an FFT?

Is there a name for the result of summing the bins of a DFT? I don't mean to sum the squares of the bins, but to simply add the magnitude of the frequency bins together to get a single result. Is ...
0
votes
1answer
49 views

Finding n Amplitudes by DFT, what is correct normalization

Please forgive me if this has already been asked. Let us assume an example with $x(t) = \sum_{i=1}^N A_i \sin(2\pi f_i t) $ under a given sampling frequency $f_s$, frequencies $\omega_i$ and ...
0
votes
3answers
57 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
1answer
75 views

Proving that the IDTFT is the inverse of the DTFT?

The DTFT is given by: $$X(e^{j\omega}) = \sum_{n=-\infty}^{\infty}x[n]e^{-j\omega n}$$ The IDTFT is given by: $$x[n]=\frac{1}{2\pi}\int_{0}^{2\pi}X(e^{j\omega})e^{j\omega n}d\omega$$ I have been ...
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 ...