Mean, variance and correlation function of Langevin equation

PhoenixWright
Messages
20
Reaction score
1

Homework Statement



I have simulated Langevin equation (numerically in Matlab) for some specific conditions, so I have obtained the solution ##X(t)##.

But now, with the solution I have obtained, I have to calculate ## <X(t)|x_0>, <X^2(t)|x_0>-(<X(t)|x_o>)^2 ## and the conditional correlation function ##<X(t)X(0)>##.

##X(0)=1##

The Attempt at a Solution



I know how to calculate the mean and variance, even for several variables. Nevertheless, I don't know how to calculate the mean and variance of ##X## conditioned to a single value ##x_0##, nor conditional correlation function...
I would appreciate if someone could help me.
Thank you so much.
 
Physics news on Phys.org
I assume that you have a certain volume filled with particles governed by the Langevin equation. The state of each particle is given by X(t). If you want to know the unconditional mean in the volume, you compute it using all particles in your volume. For a conditional mean you simply discard all particles that do not fulfill your condition. So if you want to know the mean velocity of all particles that have a mass of 1, your condition is ##< U | m = 1>##

I'm not sure if this is what you are doing, so please give more details if I'm not describing your problem correctly.

Also note that the error in computing the variance depends much on the numerical method used for the stochastic part of the Langevin equation, see e.g. the book of Kloeden and Platen.
 
  • Like
Likes PhoenixWright
bigfooted said:
I assume that you have a certain volume filled with particles governed by the Langevin equation. The state of each particle is given by X(t). If you want to know the unconditional mean in the volume, you compute it using all particles in your volume. For a conditional mean you simply discard all particles that do not fulfill your condition. So if you want to know the mean velocity of all particles that have a mass of 1, your condition is ##< U | m = 1>##

I'm not sure if this is what you are doing, so please give more details if I'm not describing your problem correctly.

Also note that the error in computing the variance depends much on the numerical method used for the stochastic part of the Langevin equation, see e.g. the book of Kloeden and Platen.

Thank you for your answer.
Actually, the problem doesn't explain the physical situation we are solving. It only asks us to solve a certain Langevin equation with the initial condition ##X(0)=1##, and then to calculate ##
<X(t)|x_0>, <X^2(t)|x_0>-(<X(t)|x_o>)^2, <X(t)X(0)>##. So, I guess we can think a situation like the one you have explained. And that's more or less what I can't understand, as I have a very specific solution for ##X(t)## (I used Euler-Maruyama aproximation for solving it, by the way), but it's not specified what ##x_0## is, so I don't know which condition I should impose in order to calculate ##<X(t)|x_0>, <X^2(t)|x_0>-(<X(t)|x_o>)^2, <X(t)X(0)>##...
 
Thread 'Need help understanding this figure on energy levels'
This figure is from "Introduction to Quantum Mechanics" by Griffiths (3rd edition). It is available to download. It is from page 142. I am hoping the usual people on this site will give me a hand understanding what is going on in the figure. After the equation (4.50) it says "It is customary to introduce the principal quantum number, ##n##, which simply orders the allowed energies, starting with 1 for the ground state. (see the figure)" I still don't understand the figure :( Here is...
Thread 'Understanding how to "tack on" the time wiggle factor'
The last problem I posted on QM made it into advanced homework help, that is why I am putting it here. I am sorry for any hassle imposed on the moderators by myself. Part (a) is quite easy. We get $$\sigma_1 = 2\lambda, \mathbf{v}_1 = \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix} \sigma_2 = \lambda, \mathbf{v}_2 = \begin{pmatrix} 1/\sqrt{2} \\ 1/\sqrt{2} \\ 0 \end{pmatrix} \sigma_3 = -\lambda, \mathbf{v}_3 = \begin{pmatrix} 1/\sqrt{2} \\ -1/\sqrt{2} \\ 0 \end{pmatrix} $$ There are two ways...
Back
Top