       Re: Problem with the zero-term of Fourier[]

• To: mathgroup at smc.vnet.net
• Subject: [mg19634] Re: Problem with the zero-term of Fourier[]
• From: Paul Abbott <paul at physics.uwa.edu.au>
• Date: Mon, 6 Sep 1999 04:20:47 -0400
• Organization: University of Western Australia
• References: <7qsftj\$17c@smc.vnet.net>
• Sender: owner-wri-mathgroup at wolfram.com

```Hans Steffani wrote:

> If have a long list of measured values in ixr.
> (Plus@@ixr)/Sqrt[Length[ixr]]
> gives
> -146.263
>
> but
> Fourier[ixr][]
> gives
> -56.1317-3.03437 I
>
> which is complex although it should be real (ixr is a list of
> real numbers) and it is different from the first term althugh
> it should be the same.

Assuming that ixr is real then I can only guess that the problem
arises from the bug with packed arrays of complex numbers?

I don't see this behaviour though. For the list

In:= ixr=Table[Random[], {20000}];

here is the DC term,

In:= Plus @@ ixr/Length[ixr]
Out= 0.498249

which agrees with the Fourier result:

In:= Chop[Fourier[ixr][]/Sqrt[Length[ixr]]]
Out= 0.498249

> I also tried
> 1/Sqrt[Length[ixr]] *
>   Sum[ixr[[tau]] Exp[2 Pi I(tau-1)(1-1)/Length[ixr]],{tau,Length[ixr]}]
> which also gives -146.263

In:= Sum[ixr[[tau]]*E^((2*Pi*I[tau - 1]*(1 - 1))/Length[ixr]),
{tau, 1, Length[ixr]}]/Length[ixr]
Out= 0.498249

> What is the problem?

You need to use the appropriate normalization in each case.

____________________________________________________________________
Paul Abbott                                   Phone: +61-8-9380-2734
Department of Physics                           Fax: +61-8-9380-1014
The University of Western Australia
Nedlands WA  6907                     mailto:paul at physics.uwa.edu.au
AUSTRALIA                            http://physics.uwa.edu.au/~paul

God IS a weakly left-handed dice player
____________________________________________________________________

```

• Prev by Date: Re: FindMinimum to minimize e . w . e
• Next by Date: Re: Help: Outer of a list of lists
• Previous by thread: Problem with the zero-term of Fourier[]
• Next by thread: Re: Problem with the zero-term of Fourier[]