| |
 |
|
|
Science Forum Index » Statistics - Math Forum » Creating White Noise in Matlab...
Page 1 of 1
|
| Author |
Message |
| miki... |
 Posted: Mon May 05, 2008 6:51 am |
|
|
|
Guest
|
Hello All,
How can I generate a Gaussian white noise process in Matlab?
To be more specefic, I want to generate the process x(t)=w(t) where
E[w(t)] = 0
E[w(t)w(s)] = (Q*Q)*dirac(t-s)
Well, In matlab (in discrete form), If (t_(k+1) - t_k) = dt < 1
then one might think that x_t_k = Q*randn is the white noise.
But, I dont think so. something is missing here that should include
the discrete sampling dt.
For example, in order to solve the random walk equation
x_dot = w(t)
then x_t_(k+1) = x_t_k + Q*sqrt(dt)*randn is the correct solution.
BECAUSE THERE IS NO WHITE NOISE (as defined above) IN DISCRETE SPACE.
So what is the correct form of just white noise process in Matlab as I
asked in the beginning of the question.
Thanks in advance,
Miki |
|
|
| Back to top |
|
| Herman Rubin... |
| Posted: Mon May 05, 2008 9:06 pm |
|
|
|
Guest
|
In article <c420f084-23f9-4c3e-891f-0ce8a527c581 at (no spam) 24g2000hsh.googlegroups.com>,
miki <miki.livne at (no spam) gmail.com> wrote:
Quote: Hello All,
How can I generate a Gaussian white noise process in Matlab?
To be more specefic, I want to generate the process x(t)=w(t) where
E[w(t)] = 0
E[w(t)w(s)] = (Q*Q)*dirac(t-s)
Well, In matlab (in discrete form), If (t_(k+1) - t_k) = dt < 1
then one might think that x_t_k = Q*randn is the white noise.
But, I dont think so. something is missing here that should include
the discrete sampling dt.
For example, in order to solve the random walk equation
x_dot = w(t)
then x_t_(k+1) = x_t_k + Q*sqrt(dt)*randn is the correct solution.
BECAUSE THERE IS NO WHITE NOISE (as defined above) IN DISCRETE SPACE.
So what is the correct form of just white noise process in Matlab as I
asked in the beginning of the question.
Thanks in advance,
Miki
In a sense, one can get a discrete white noise process;
one cannot get a continuous process of this type.
It is in continuous time that there is no such process.
The Wiener process does not quite satisfy a Lipshitz
condition of order 1/2, let alone enough to have a
derivative.
--
This address is for information only. I do not claim that these views
are those of the Statistics Department or of Purdue University.
Herman Rubin, Department of Statistics, Purdue University
hrubin at (no spam) stat.purdue.edu Phone: (765)494-6054 FAX: (765)494-0558 |
|
|
| Back to top |
|
| |
|
Page 1 of 1
All times are GMT - 5 Hours
The time now is Mon May 12, 2008 4:35 am
|
|
