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

learn more… | top users | synonyms

6
votes
2answers
109 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 ...
0
votes
0answers
23 views

Discrete Fourier Transform and data manipulate

I solve a system od ODE and i save the discrete solution (i.e.: x,y,vx,vy, or a subset of it) in a file. Then i try to call the DFT on this file but the problem is that i must have also the time in ...
11
votes
4answers
824 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 ...
10
votes
1answer
482 views

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

I'm trying to Fourier transform a one dimensional list of a time history of some quantity using the Fourier function. I'm interested in the frequency spectrum, but ...
2
votes
0answers
46 views
3
votes
1answer
151 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: ...
4
votes
1answer
162 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 ...
1
vote
2answers
96 views

2D Fourier transform of a table - ArrayPlot issue

Recently, I have encountered an issue with ArrayPlot after performing a Fourier transform of a table. Picture presented above is an ArrayPlot of a 2D table. Using Fourier on this 2D table I obtain ...
3
votes
0answers
95 views

Fourier transformation of HeavisideTheta functions

I want to find 2D-Fourier transformation of the function given below f = HeavisideTheta[y1]*HeavisideTheta[y2 - y1] For the purpose, I use built-in function in ...
5
votes
1answer
147 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, ...
0
votes
2answers
89 views

How to obtain the form of sin(nx) and cos(nx) from the result of FourierSeries [closed]

How to obtain the form of sin(nx) and cos(nx) from the result of FourierSeries[]. For Example, $$\text{FourierSeries}[x,x,5]=i e^{-i x}-i e^{i x}-\frac{1}{2} i e^{-2 i x}+\frac{1}{2} i e^{2 i ...
1
vote
1answer
143 views

Plot Fourier transform of $\sin (2 t)$

I'm trying to plot the Fourier transform of $\sin (2 t)$. I tried using the FourierTransform Function (I'm expecting a peak at w = 2) but it gives me undefined at ...
0
votes
0answers
54 views

Recognizing a generalized function as zero

Consider the following expression: HeavisideTheta[x]HeavisideTheta[-x] Mathematica will not simplify it further. But $\Theta(x)\Theta(-x)$ is zero if considered ...
4
votes
1answer
192 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 ...
2
votes
1answer
148 views

How do you use fourier transforms to perform a deconvolution

I've been trying to teach myself about deconvolution through fourier transforms, but I seem to be missing something simple, as my results are garbage. I start by defining a test function and a window ...
0
votes
1answer
73 views

ListLinePlot “Cannot Take Positions” error with Fourier of WAV file

When I try to do a Fourier analysis of a WAV file with a sample rate of 22050 I get an error "Cannot take positions". What am doing wrong? ...
6
votes
0answers
110 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] ...
1
vote
1answer
115 views

Fourier Series of “Split” Defined Function

Not sure how to phrase this in a concise way, anyway it seems like the function FourierSeries assumes the interval for which to compute the Fourier coefficients of a given function is $[-\pi,\pi]$, ...
4
votes
1answer
148 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: $ ...
9
votes
2answers
435 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 ...
2
votes
1answer
127 views

Integrate message: can't prove Integration limits are real

I'm working on solving differential equations through Fourier series, I made a function to help me calculate the coefficients that looks like this: ...
1
vote
1answer
64 views

Two dimensional NIntegration failure - non-numerical values

I've got this function: E0[x0_, y0_, z_]:= A/w[z]*Exp[-(x0^2 + y0^2)/(w[z]*w[z])]*Exp[(I*2*Pi*(x0^2 + y0^2))/(λ*2*R[z])]*Exp[I*ϕ[z]]; Where the ...
0
votes
1answer
99 views

Fourier-Analysis How to display every period correctly? [duplicate]

I want to express a function as a sum of fourier series. I tried the code below. It works. But only get one period (from 0 to 1). How to display every period correctly? ...
2
votes
1answer
136 views

Why does FourierTransform[] converge while same integral manually written does not?

I'm studying the quantum mechanics of an infinite square well from a computational standpoint. My eigenfunctions are defined as $u_n(x)=\sqrt{2/L}\sin\left(n\pi x/L\right),\quad 0 \le x \le L, \quad ...
3
votes
3answers
92 views
0
votes
1answer
77 views

Is there a way to plot discrete functions in mathematica? [closed]

Suppose you have a discrete function Triangle[n/5]. How would you plot it in mathetica?
11
votes
1answer
428 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 ...
2
votes
1answer
179 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 ...
0
votes
2answers
81 views

fourier issue arising from input miscommunication

xdomain = Table[i , {i, -10, 10, .1}]; ListPlot[InverseFourier[ Fourier[E^-#^2 & /@ xdomain]*E^#^2 & /@ xdomain]] So I want to numerically fourier ...
4
votes
0answers
100 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 ...
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 ...
0
votes
1answer
149 views

Plotting partial sums of Fourier sine series

How do I plot this on Mathematica version 5.2? $\frac{4}{\pi} \sin{x} + \frac{4}{3 \pi} \sin{3 x} + \cdots + \frac{4}{(2 N+1) \pi} \sin{(2 N+1) x}$ over $x \in [-\pi,\pi]$ for $N= 3, 6, 12, 24$. I ...
3
votes
2answers
403 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, ...
6
votes
2answers
433 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 ...
7
votes
2answers
4k 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 ...
6
votes
2answers
4k 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.
10
votes
1answer
317 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 ...
0
votes
0answers
85 views

Cross-correlation in SDEs

Is it possible to derive a cross-correlation function between a stochastic variable and a state-variable in an SDE, such as for the simple model here, or better between two state variables in a two ...
5
votes
2answers
428 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 ...
6
votes
0answers
110 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 ...
0
votes
1answer
136 views

Error messages when using NInverseFourierTransform

I have two functions that I need to inverse Fourier transform and I was trying to get Mathematica to help me. I tried simply using theInverseFourierTransform ...
8
votes
2answers
355 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: ...
8
votes
1answer
775 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 ...
10
votes
0answers
238 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 ...
1
vote
1answer
251 views

Problem with Fourier and InverseFourier function

I have the following code Image[img, ImageSize -> 300] data = ImageData[img]; fdct = Fourier[data] tdata = InverseFourier[fdct]; Image[tdata] but it seems to ...
-5
votes
1answer
1k views

Determine the 2D-Fourier transform of an Image

I am new to Mathematica and using Version 8.0. I want to calculate the 2D-Fourier Transform of a digital image and see its Magnitude and Phase Spectrum, and reconstruct it.
8
votes
1answer
409 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 ...
17
votes
2answers
2k 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 ...
6
votes
2answers
274 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. ...
15
votes
4answers
736 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 ...