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

learn more… | top users | synonyms

3
votes
0answers
129 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
338 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, ...
1
vote
1answer
410 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
80 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
256 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 ...
3
votes
1answer
346 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
2answers
120 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 ...
0
votes
1answer
153 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? ...
8
votes
0answers
142 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
218 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]$, ...
5
votes
1answer
222 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: $ ...
10
votes
2answers
1k 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
266 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
79 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
128 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
174 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
95 views

Filling in explicit zeros in a (x,y) list for nonexisting y values for discrete Fourier transform in Mathematica

I have a list of (x,y) values in Mathematica for various discrete x values, as in ...
0
votes
1answer
86 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?
2
votes
1answer
250 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
96 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 ...
5
votes
0answers
117 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 ...
0
votes
1answer
260 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 ...
12
votes
1answer
1k 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 ...
6
votes
2answers
682 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 ...
4
votes
2answers
650 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, ...
10
votes
1answer
421 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 ...
5
votes
2answers
663 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 ...
7
votes
0answers
155 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
2answers
451 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: ...
1
vote
1answer
348 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 ...
-4
votes
1answer
3k 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.
9
votes
1answer
633 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 ...
8
votes
2answers
406 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. ...
16
votes
4answers
985 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 ...
9
votes
1answer
1k 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 ...
2
votes
1answer
105 views

Tough Inverse Fourier Integral: why does the sign matter?

I'm computing the Inverse Fourier Integral of ((a^2 + omega ^2) c^2) /((b^2 + omega^2) ((r^2 + omega^2 - omegaInt^2)^2 + (2 omegaInt r)^2)) ...
0
votes
1answer
301 views

How to manipulate a Fourier sine series?

I would like to plot the Fourier sine series for $f(x) = x$. I would like to manipulate the parameter n from $0$ to $10$ and plot the sine waves over $(0,1)$. My ...
9
votes
1answer
421 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 ...
2
votes
1answer
736 views

Fourier transform of sampled data

I have got some impulse response data that I would like to transform via Fourier to get the amplitude-frequency characteristics of the performing loudspeaker. The ...
0
votes
1answer
150 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 ...
1
vote
0answers
446 views

Doing local FFT on huge 3D vector data cell mesh and visualizing it spatially?

Simulation type: I'm running a simulation with the OOMMF micromagnetics package http://math.nist.gov/oommf/ where are magnet is represented by a mesh of 3 million cells, it gets excited by a ...
1
vote
2answers
1k views

Discrete FFT of non-periodic signal excited by short pulse

Say I have a signal F(t) which represents a excitation by a pulse (so it is not periodic and declining to zero amplitude). According to this answer, for instance, ...
0
votes
1answer
591 views

How can I find the Fourier series from discrete data?

I would like to define the function with Fourier series when I just have discrette data. (I don't have a specific function). I can draw the graph using the data, but I don't have a idea how can I get ...
8
votes
2answers
1k 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, ...
1
vote
1answer
471 views

Fourier transformation of solution of differential equation

I'm trying to plot a Fourier transform of solution of differential equation. I have tried with: ...
9
votes
2answers
6k 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 ...
12
votes
1answer
507 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
2answers
1k views

Obtaining the Fourier transform of an operator

We know that if we have a function $f(x)$, and we call $g(\omega)$ its Fourier transform, then the Fourier transform of $x f(x)$ is $$\imath \frac{\mathrm{d} g(\omega))}{\mathrm{d}\omega} $$ and ...
7
votes
1answer
551 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 ...
2
votes
0answers
302 views

Fourier Coefficients in Mathematica [duplicate]

I'm calculating Fourier Coefficients by hand and trying to verify them in Mathematica. However, in Mathematica I get them wrong by a factor of 2. Is there some part of the Mathematica functions I'm ...