Clear Filters
Clear Filters

fft and random sequence

5 views (last 30 days)
yasser
yasser on 9 May 2014
Answered: yasser on 10 May 2014
generating 2 random sequences with randn with same variance then applying fft to them would result in 2 sequences of non = variances why is that

Answers (3)

the cyclist
the cyclist on 9 May 2014
Because the sample variance is not necessarily equal to the population variance.
  2 Comments
yasser
yasser on 9 May 2014
please explain more
the cyclist
the cyclist on 9 May 2014
Here is the simplest possible explanation I can think of. Imagine you use the command
randn(2,1)
to draw a sample with two values from a distribution that is N(0,1). (In other words, a normal distribution with population mean of 0 and standard deviation of 1.)
Suppose that draw gives you values of [0.7824 -0.5000]. The sample mean of that sample is 0.1412. The sample standard deviation is 0.9068.
Do you get upset that the sample mean is not exactly zero, and that the sample standard deviation is not exactly 1? Hopefully not. It's a random sample, with random variation.
Similarly in your example, there is random variation.

Sign in to comment.


Star Strider
Star Strider on 9 May 2014
The Fourier transform integrates over the original x variable, so the Fourier transform is a function of the (complex) frequency variable only. You are taking the mean and standard deviation (variance) of F(w) rather than f(x).
So if you have different lengths of the vectors in the original ( x or t ) domain, the means and variances should be close to the same in that domain. Without appropriate scaling for different lengths or sampling frequencies, the means and variances (taken over the frequency variable) of the Fourier transforms of those data would not be the same in the complex frequency domain.
I’m not certain what you did, so this is only theoretical. The analytical maths are not trivial, and neither the Symbolic Math Toolbox or my back-of-the-envelope analysis yielded a definitive answer.

yasser
yasser on 10 May 2014
thanks u all

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!