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)>##...
 
Hello everyone, I’m considering a point charge q that oscillates harmonically about the origin along the z-axis, e.g. $$z_{q}(t)= A\sin(wt)$$ In a strongly simplified / quasi-instantaneous approximation I ignore retardation and take the electric field at the position ##r=(x,y,z)## simply to be the “Coulomb field at the charge’s instantaneous position”: $$E(r,t)=\frac{q}{4\pi\varepsilon_{0}}\frac{r-r_{q}(t)}{||r-r_{q}(t)||^{3}}$$ with $$r_{q}(t)=(0,0,z_{q}(t))$$ (I’m aware this isn’t...
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
Back
Top