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

learn more… | top users | synonyms

2
votes
0answers
49 views

FourierTransform of a gaussian at -1 and 2

Let ft[w_] := FourierTransform[Exp[-t^2], t, w] Then ft[#] & /@ Range[-2, 3] Evaluates to ...
1
vote
1answer
196 views
7
votes
1answer
461 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 ...
1
vote
0answers
59 views

Discrete Fourier Transform baseline subtraction

I am trying to remove sinusoidal variability from a set of evenly-spaced intensity vs. time data. These data contain periodic events that are separated by a normalized (but sinusoidally variable) ...
0
votes
1answer
101 views

Fourier series of the odd extension of a function [duplicate]

Using Mathematica to make a plot, I noticed that the interval it's defined over must be changed in order to do the Fourier series. How do I change the interval to be on an interval of $[-{\cal l}, ...
5
votes
2answers
1k 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, ...
1
vote
1answer
167 views

DFT: why are frequencies mirrored around the frequency axis?

I want to apply Mathematica's DFT procedure, Fourier[], to a Gaussian. Why a Gaussian? Because I know the analytical result. So I define ...
10
votes
1answer
218 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 ...
8
votes
0answers
156 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 ...
1
vote
1answer
98 views

Problem of Discrete Fourier?

I tried Fourier Transform for a range of figures. Firstly, I tried in this way. ...
0
votes
0answers
66 views

3-adic FFT algorithm

I am trying to implement the FFT algorithm when $n=3^k$, but maybe my code is wrong. Are there any suggestions for improving my code? ...
1
vote
1answer
127 views

Finding the discrete Fourier transform for a simple oscillator?

Here's what I want to do, initially at least: I have an oscillator, it's in some force field, I've integrated it using the Euler method, so I have a list of points of the form {t,x}, where t is the ...
1
vote
1answer
182 views

Precision while calculating Fourier transform

I'm trying to understand how precision works in Mathematica. Particularly I'm calculating discrete Fourier transform using the Fourier function and calculating it ...
0
votes
0answers
62 views

Monitor Progress of a function (Fourier[] in my case)

I am using PSD = Abs[Fourier[normedAutoValues, FourierParameters -> {-1, 1}]]; normedAutoValues contains around 200 *10^7 ...
3
votes
1answer
282 views

Approximate $h$ in $F(\theta)=\sin \theta \int_{-L}^{+L}h(z)e^{-ikz\cos \theta} \,dz$

Consider $$F(\theta)=\sin \theta \int_{-L}^{+L}h(z)e^{-ikz\cos \theta} \,dz$$ $$|z|\le L$$ $$0 \le \theta \le \pi$$ By having knowledge of $F(\theta)$, how can one approximate $h(z)$? In ...
6
votes
2answers
250 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 ...
0
votes
0answers
74 views

Using Fourier to return exponential function

I am trying to use Fourier to numerically demonstrate the following identity: $$ \frac{1}{2\pi}\int_{-\infty}^{\infty}\frac{e^{i\,s\,y}}{1 + a\,s}ds=e^{-y/a} $$ I'm getting the correct shapes, but my ...
2
votes
2answers
311 views

Using Periodogram to investigate the seasonality of the data

After asking how to get rid of the seasonality of this set of data, I'd like to ask how to investigate the seasonality in data using Mathematica. Which tool would you use to study seasonality? I.e. ...
0
votes
0answers
40 views

InverseFourier of nested lists

last week I've got an Mathematica package that is using several inverse Fourier transforms on a data cube ($n\times n\times n$-array). I'm implementing the same functionality in C++ and tried to ...
3
votes
0answers
100 views

Different FourierParameters in each dimension for discrete Fourier transform

I have a query about the Fourier utility, that does discrete Fourier transforms. I am unable to find any information on whether different FourierParameters can be used in each dimension. This is a ...
2
votes
0answers
140 views

Calling CUDA cuFFT

Following ybeltukov's excellent answers (see here and here) about how to access the CUDA library functions and interact with Mathematica through CUDAMemory, I tried to compile the code given in ...
1
vote
2answers
102 views

Why does InverseFourierTransform gives different result to doing it by hand?

When I can calculate an inverse fourier transform by the built-in function InverserFourierTransform, the result is different from what I calculate by definition ...
3
votes
1answer
547 views

Animated Wave Propagation using Fourier & InverseFourier

This is a continuation off of previous help on the first part of my project: fourier issue arising from input miscommunication Now I want to go one step further in the current code. Here's the code ...
4
votes
0answers
217 views

Help with Fourier analysis and filtering of an image

I am starting with the following image: For the curious, it is an STM scanned image of graphite. I am hoping to find some kind of periodic structure in the image related to the honeycomb structure ...
1
vote
0answers
56 views

PeriodogramArray or PowerSpectralDensity?

could anyone explain what is the difference between PowerSpectralDensity[data,omega] and PeriodogramArray[data]? When I plot them (starting from the same set of data of course) I get two different ...
1
vote
2answers
268 views

Get Y-Axis of FFT in decibels

I am trying to make a plot of a sound file similar to Audacity's FFT plot which looks like this when the x-axis is on a log scale: My plot of the same sound file (below) has the problem that the ...
0
votes
0answers
41 views

FourierSinCoefficients extremely slow

I'm just playing with Fourier decomposition of periodic functions. We had a similar example with standing waves and cosines in class, so I tried to expand a triangle wave in sines. This is what I came ...
13
votes
3answers
562 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. ...
4
votes
1answer
343 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 ...
1
vote
2answers
439 views

How to calculate the phase spectrum

I have the following frequency characteristics in the Fourier domain: $$H(\omega)=\frac{-\omega^2}{63170s^{-2}-\omega^2+355.1s^{-1}i\omega}$$ How do I find the phase spectrum from this? I should ...
1
vote
1answer
93 views

Changing Function Domain of Fourier Series

If I write FourierSinSeries[x, x, 10] Mathematica will output the first 10 terms of the Fourier Sine Series f(x) = x defined for 0 < x < Pi. How do I output the the Fourier Sine Series for ...
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, ...
2
votes
0answers
100 views

Incorrect output from InverseFourierTransform

Bug introduced in 8.0.4 or earlier and fixed in 9.0.1 or earlier I am trying to compute the inverse Fourier transform of this expression ...
0
votes
0answers
49 views

Multidimensional Fourier PDE solution and and exact solution are scaled by $(2 \pi)^{2}$

I'm solving a partial differential equation (Stokes equation) in a periodic domain in three dimensions with fourier transforms. I'm using Fourier[] mathematica routines to perform the discrete ...
1
vote
1answer
196 views

Extracting out the Fourier modes of a series, and using these to reconstruct data series

I would like to extract the most important frequency modes from a data set which exhibits strong annual periodicity, as well as some (less) important shorter-term frequency components. The following ...
2
votes
1answer
262 views

Fourier Transform on a Sound

I don't really know much above Fourier Transform , or how it's used in Mathematica but I want to know how I can find out, out of which Sin-waves a (static) sound is made of. So for example I take "G" ...
0
votes
1answer
217 views

Discrete Fourier Transform of a vector defined on a 3D lattice

In my research I need to compute the Discrete Fourier transform of a vector defined on a 3D lattice (a cube) to the "reciprocal" lattice. This is quite new to me so before proceeding i wanted to be ...
2
votes
0answers
50 views

Precision loss in FourierDCT

I have high precision data that I want to do Fourier transform with, but it looks like the precision is lost dramatically. ...
3
votes
2answers
145 views

Compare two shifted signals

I have one set of original data (signal A) and a second set of derived data (signal B) which was computed from A by position rotation (circular shift). Here are the data for the signals (A and B): ...
1
vote
0answers
326 views

Using 2D Fourier transform of an image to detect typical wavelengths

I am new to Mathematica, and I have read this post to understand how to perform Fourier transform on an image. My mission is to extract information on the typical distance between the black patches in ...
1
vote
0answers
123 views

The inverse Fourier transform of a rather complex function

I am using Mathematica 7.0 and trying to find the Fourier transform of a rather complex function g. The function g is defined by ...
5
votes
0answers
190 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 ...
1
vote
0answers
135 views

Inverse Fourier Transform - wolfram alpha vs Mathematica 9.0

I'm facing an interesting problem. I'm trying to do this simple task: InverseFourierTransform[0.031622*Exp[0.4995ω^2], ω, t] Mathematica can't handle it for some ...
9
votes
2answers
244 views

InverseFourierTransform of simple function takes forever

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

Finding periodicity of multidimensional data?

I have a data set consisting of measurements of an $n$-dimensional vector $\mathbf{x}(t)$, at $N$ distinct time points $t_i$. The time points $t_i$ are equally spaced. I need to find the dominant ...
2
votes
1answer
244 views
1
vote
1answer
720 views

Discrete Fourier Transform with x axis data?

Usually, Fourier Transform in mathematica is supplied by y data only, like the famous Sin[x] transform example. What if my data is a list of ...
1
vote
0answers
31 views

Inverse Fourier transform to get real function [duplicate]

I have got a sampled complex function, $H(\omega)$ and know that its inverse Fourier transform has to be real function $h(t)$, which means that $H(-\omega)=H^{\ast}(\omega)$ has to be fulfilled. If we ...
1
vote
1answer
223 views

Inverse Fourier transform

I am dealing with an inverse Fourier transform of a sampled complex function, $H(\omega)$. I want to obtain a real transformed function $h(t)$, which means that $H(-\omega)=H^{\ast}(\omega)$ has to be ...
6
votes
4answers
854 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]]