MathGroup Archive 2001

[Date Index] [Thread Index] [Author Index]

Search the Archive

RE: Fourier Transform and convolution.

  • To: mathgroup at
  • Subject: [mg26617] RE: [mg26583] Fourier Transform and convolution.
  • From: Adil.Reghai at
  • Date: Sat, 13 Jan 2001 22:36:01 -0500 (EST)
  • Sender: owner-wri-mathgroup at


I was doing a Fourier transform using the Excel Tool Pack and I have not
zero padded before the computation.
Now with Mathematica, I have tried your codes and it seems to work very



-----Original Message-----
From: BobHanlon at [mailto:BobHanlon at]
To: mathgroup at
Subject: [mg26617] Re: [mg26583] Fourier Transform and convolution.


dist = NormalDistribution[mu, sigma];

p = PDF[dist, x]

1/(E^((x - mu)^2/(2*sigma^2))*(Sqrt[2*Pi]*sigma))

f = FourierTransform[p, x, w, FourierParameters -> {1, 1}]//PowerExpand

E^(I*mu*w - (sigma^2*w^2)/2)

f == CharacteristicFunction[dist, w]


InverseFourierTransform[f^2, w, x, 
    FourierParameters -> {1, 1}]//PowerExpand

1/(E^((x - 2*mu)^2/(4*sigma^2))*(2*Sqrt[Pi]*sigma))

This is N(2*m, Sqrt[2] * sigma)

% == PDF[NormalDistribution[2*mu, Sqrt[2]*sigma], x]


Bob Hanlon

In a message dated 2001/1/9 2:17:00 AM, Adil.Reghai at writes:

>We wanted to perform the distribution of a sum of two independent variables
>X and Y which are both N(0,1).
>The way We wanted to do it is to use the convolution property. In order
>compute this convolution we perform 
>a Fourier Transform of the pdf n(0,1) and then square each point. After
>we performed an inverse Fourier transform. 
>Instead of obtaining a N(0,sqrt(2)) we obtain a distribution which is
>shifted in the form of a U.

If you have received this e-mail in error or wish to read our e-mail 
disclaimer statement and monitoring policy, please refer to or contact the sender.

  • Prev by Date: Re: Equation
  • Next by Date: Re: Sum of Squares
  • Previous by thread: Re: Fourier Transform and convolution.
  • Next by thread: worldplot GIF