Questions on the discrete and continuous Fourier analysis functions of Mathematica, as well as the FourierSeries` package.

learn more… | top users | synonyms

20
votes
2answers
3k views

How to make Fourier behave like FourierTransform?

I'm not very experienced with Fourier Transforms, so there may be something inherently wrong with attempting to do this, but how can I make the discrete Fourier behave like the continuous ...
20
votes
3answers
2k views

Can one find the beat of a tune with Fourier analysis?

I'm trying to find out if it's possible to find the beat of a tune by Fourier analysis with Mathematica. I'm taking a 44.1 kHz sample sound and hoping that I might get a nice peak for a frequency ...
18
votes
4answers
4k views

Numerical Fourier transform of a complicated function

Say I have a function $f(x)$ that is given explicitly in its functional form, and I want to find its Fourier transform[1]. If $f$ is too complicated to have an analytic expression for $\hat f(k)$, how ...
17
votes
4answers
1k views

Digital filter of image in Fourier space

Consider the contour plot below, of a function f(x,y) - which is calculated numerically by solving some equations. If one looks carefully, one can observe two types of waves, to the left and right of ...
16
votes
4answers
4k views

Plotting Fourier spectrum versus frequency of a signal

I have looked around here, and i am sure this has been answered, but i don't understand it. The thing is, I have taken a introductory signal processing course, and we had to use Mathematica, and i had ...
15
votes
1answer
2k views

What's the correct way to shift zero frequency to the center of a Fourier Transform?

I'm trying to apply a Fourier transform of a one dimensional list of a time history of some quantity using the Fourier function. I'm interested in the frequency ...
13
votes
3answers
388 views

Disentangling the data

An experiment yields a functional dependence on two variables $F(x,y)$, can be imagined as a 2D map. It is known that the function can be factored as follows $F(x,y)=fx(x) gy(y) +gx(x) fy(y)$. Given ...
13
votes
1answer
4k views

Coulomb potential as a Fourier transform

It is well known from theory that the Coulomb potential can be obtained as a Fourier transform in the following way: $$ \int \frac{\mathrm{d}^3p}{\left( 2 \pi \right)^3} \frac{e^{\mathrm{i} ...
13
votes
3answers
509 views

Next highly composite number?

R language has this function 'nextn' (link) which computes the next highly composite number greater than a given one, which is used to find the optimal padding size for the subsequent FFT operation. ...
12
votes
1answer
575 views

How can I obtain a zero-filled Fourier transform, as produced by MATLAB's fft?

The MATLAB FFT function, fft(X,n), can be used to return the n-point DFT. This effectively extends the original signal, X, to ...
12
votes
2answers
297 views

Parallelization of distinct array write access from subkernels

I'm working on an implementation of a multivariate FFT, which is (or at least should be) highly parallelizable due to the row-column-algorithm. However, i can't figure out how to implement that. The ...
11
votes
2answers
2k views

Finding fourier transform of data and hence frequencies

I looked at all the other questions related to mine before posting this, and they didn't solve my problem. I have a large data set which can be downloaded from here. I'm using the following code to ...
11
votes
2answers
319 views

Rotational Averaging of Image Periodograms (FFTS)

I'm writing a code to do rotational averages of image FFT's. I have a working version but it is pretty slow. Wondering if anyone has any thoughts on a way to do this more efficiently. Data looks ...
11
votes
3answers
2k views

How to approximate a given WAV file with trigonometric series?

I'd like (together with a few people) to prepare a presentation about Fourier series for middle/high school students. I thought it might be quite cool to play a violin sound from, say, a WAV file, ...
11
votes
2answers
8k views

Plotting the frequency spectrum of a data series using Fourier

testData = Table[N@Sin[500 x], {x, 0, 100}]; ListLinePlot[Abs[Fourier[testData]], PlotRange -> Full] Gives me Which I do not expect because the Fourier ...
11
votes
0answers
369 views

Proving (self) similarity with Mathematica - Reccurrence Plots, Similarity Plots etc

I posted this question in math.se but given the sheer tumultuous number of questions that keep appearing on math.se and also given that I am trying to accomplish this in Mathematica, I thought I'd ...
10
votes
1answer
480 views

Why does Mathematica return a Fourier transform for a function for which it is not defined?

The following function $$g(x) = (1 + x^{1/a} )^a $$ should NOT have a Fourier transform, as far as I am aware, for any real values of $a$ since $g(x)$ is not nice in the sense of decays quickly ...
10
votes
1answer
443 views

FourierSeries for rational function looks wrong

I have the following code: ser[x_] = FourierSeries[(\[Pi]^2 + a)/(3 x^2 + a), x, 10] // N // Chop It gives me some series, which I then try to plot. And ...
10
votes
1answer
395 views

Inconsistent results of second derivative of inverse fourier transform

I am trying to get the Green's function of a toy diffusion equation $$\frac{\partial^2 u(x,t)}{\partial x^2} = \frac{1}{\alpha^2}\frac{\partial u(x,t)}{\partial t}$$ with Mathematica 9. Then solve ...
10
votes
1answer
204 views

Fourier transform over a custom dimension

I have a multidimensional array A. Fourier[A] finds the discrete Fourier transform over all dimensions. How to find discrete ...
9
votes
2answers
219 views

InverseFourierTransform of simple function takes forever

i'm quite new to Mathematica. I'm trying to do something along the lines of: ...
9
votes
1answer
782 views

FourierTransform and Partial Derivatives?

I have a function of three variables $f(x,y,z)$, and it obeys a linear partial differential equation. I'm checking my by-hand calculations with Mathematica. I want to convert the PDE into the ...
9
votes
1answer
2k views

Computation of Hankel Transform using FFT (Fourier)

To address circular symmetric cases of 2D Fourier Transformations the so called Hankel Transform can be applied (for a detailed derivation of the relation between the 2D Fourier transform and the 1D ...
9
votes
0answers
161 views

Possible bug with FourierTransform linearity

Each component is easily transformed but the sum is not: FourierTransform[f''[x], x, k] FourierTransform[f'[x], x, k] ...
8
votes
3answers
5k views

Plotting Partial Sums of Fourier Series

I'd like to plot some partial sums for a Fourier Series problem, but I am not sure if the output I am getting is correct. I want to be able to plot the partial sums and the function on the same graph. ...
8
votes
2answers
574 views

Why Fourier doesn't show me the peaks?

I'm trying to identify the frequencies in my time history samples, and I can see a frequency in the time history, but can't see it in its Fourier transform. Here it is : the sample data: ...
7
votes
2answers
5k views

Does Mathematica implement the fast Fourier transform?

Is there a fast Fourier transform in Mathematica? Although looking in the help I could not find one. I am looking to implement the equivalent of fft in MATLAB.
7
votes
2answers
2k views

Fourier Transform of a Step Function

I'm trying to obtain the form of a sinc function that I know I'm supposed to get in Mathematica. I'm doing this because I intend to do a lot with Fourier transforms (FT) and I'd like to know I'm not ...
7
votes
2answers
203 views

Speed up Fourier for Booleans

I'm calculating a bit of the power spectrum of a two color cellular automaton using: ...
7
votes
1answer
663 views

Speeding up numerical Fourier Transform

I wrote this function NFourierTransform, which takes a function $f(k)$ and numerically calculates the fourier transform integral for discrete values of $k \in ...
6
votes
4answers
824 views

Why this gives 0?

fixed in 10.0.2. Why this gives 0? I am in shock, it should be Sinh[x]. FullSimplify[I InverseFourierTransform[FourierTransform[Cosh[t],t,w]/w,w,x]]
6
votes
2answers
486 views

Correct fourier scaling and high-resolution frequency identification

I have a dataset of amplitude versus time $(t,A(t))$ and I need to extract the dominant frequency and amplitude, and also get the amplitude at one other specific frequency. My data looks like this: ...
6
votes
2answers
851 views

Band-pass filter an image

Ok so i have some HR-TEM images. In this technique fringes appear whereever a crystalline structure is present. This fringes are separated by a certain distance depending on the material, this is our ...
6
votes
2answers
167 views

Filtering of an unwanted frequency through a Fourier Transform

I have a data set with evenly spaced data points. The plot is frequency vs. intensity. The overall shape of the plot is an upwards curve into a plateau, this cannot be seen in the data as this is an ...
6
votes
1answer
355 views

calling CUDA library functions

Mathematica has some predefined CUDA functions, like CUDAMap, CUDAFourier, etc. My question is how can one call functions which are part of CUDA library but not part of predefined Mathematica ...
6
votes
1answer
86 views

Calculation of 2D FFT for an image

I am trying to obtain 2D FFT of an image. From the 2D FFT, I wish to calculate the amplitude of the FFT averaged over a circle vs the distance in pixels by varying the radius of the pixels which is ...
6
votes
1answer
219 views

How to plot stability region of Fourier analysis

Stability limits for Runge–Kutta methods in the complex ΩΔt-plane are needed to be plotted like this: For example, let P=1+z where ...
6
votes
0answers
147 views

The same analytical expression gets inconsistent FourierTransform results

I have an expression which is just a linear combination of plane waves, and I'd like to calculate FT of it. I know what I will get should be a bunch of delta functions, but it turns out Mathematica ...
5
votes
2answers
844 views

Frequencies represented in a DFT list (Fourier[])

I have a sound sample (sample rate: 44100). When I have a list l with n successive elements from the sample, ...
5
votes
1answer
299 views

Extracting frequencies with Fouriertransformation

I know this question came up so many times, but after browsing through many of the threads, I'm still not sure, if I'm doing it correctly with my dataset, so I'll give it a try. If have some sample ...
5
votes
1answer
489 views

Fit data with linear combinations of Cos and Sin

Recall that if $f$ is a piecewise continuous function on the interval $[-\pi, \pi]$, then the Fourier series of $f$ is $$f(x) = \sum_{n=0}^{\infty} (a_n \cos(nx) + b_n \sin(nx)).$$ See for example, ...
5
votes
1answer
256 views

Prove an identity in quantum harmonic oscillator

Problem: In the context of quantum harmonic oscillator the eigenfunctions are given by: $ u_n(x) = (N_n/\sqrt{b}) H_n(x/b) \exp\left[-x^2/(2b^2)\right] $, where $N_n$ is the normalization factor: $ ...
5
votes
2answers
818 views

How to draw a graph about the FourierTransform of a 2D Function?

Here is a 2D function, and it cannot be calculated by the FourierTransfrom ...
5
votes
0answers
154 views

Periodogram for irregularly spaced data

I am trying to replicate some existing work on the length of the cycles in certain economic data. Specificially, I am using the Schularick and Taylor (2012) data set to replicate the periodogram shown ...
4
votes
4answers
622 views

Why does Fourier give a shifted frequency?

I have a signal that I want to identify the frequencies in it, I used the Fourier function but I can't get the frequency correctly. Here is a simplified example: ...
4
votes
1answer
212 views

Reconciling results from Fourier with those form FourierTransform

I need to use the discrete Fourier transform for function that represented as list of values. I started with an easy task to check my understanding. I tried to get the amplitude values for ...
4
votes
1answer
119 views

What do the X and Y axis stand for in the Fourier Transformation in wolfram mathematica?

In wolfram mathematica the function Fourier has the following declaration Fourier[list] And after a list is given to the function, a simple Plot gives the ...
4
votes
1answer
184 views

Finding the dominant frequency in a large dataset with Fourier + periodogramm

I know there are dozens of questions about how to extract the correct frequency out of data, however I probably don't interpret my data correctly, or do not apply those answers correctly. Problem I ...
4
votes
1answer
121 views

Workarounds for a possible bug in the linearity of FourierTransform

This is somewhat similar to this question, except the problem I am encountering is to do with Fourier transforms of scalar multiples of functions and their derivatives. I wish to input ...
3
votes
1answer
212 views

Programming a recursive FFT with Mathematica

I'm trying to program a FFT using a recursive function in Mathematica, though my program is not getting me anywhere at the moment. Could you say what's wrong with it and what I can do to make it work? ...