Mean, variance and correlation function of Langevin equation

Click For Summary
SUMMARY

The discussion focuses on calculating the mean, variance, and conditional correlation function of the Langevin equation using numerical simulations in Matlab. The user has successfully simulated the Langevin equation and obtained the solution X(t) but struggles with calculating , - ()^2, and . Key insights include the importance of conditioning on specific values and the impact of numerical methods, such as the Euler-Maruyama approximation, on variance calculations. The reference to Kloeden and Platen's work highlights the significance of numerical methods in stochastic processes.

PREREQUISITES
  • Understanding of Langevin equations and their applications in stochastic processes.
  • Proficiency in numerical methods, particularly the Euler-Maruyama approximation.
  • Familiarity with conditional expectations and variance calculations.
  • Basic knowledge of Matlab for numerical simulations.
NEXT STEPS
  • Research how to compute conditional expectations in stochastic processes.
  • Study the Euler-Maruyama method in detail to improve numerical accuracy.
  • Explore the implications of conditioning on specific values in statistical mechanics.
  • Read Kloeden and Platen's book on stochastic differential equations for deeper insights.
USEFUL FOR

Students and researchers in physics and applied mathematics, particularly those working with stochastic processes and numerical simulations of Langevin equations.

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   Reactions: 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)>##...
 

Similar threads

Replies
8
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
7
Views
2K
  • · Replies 19 ·
Replies
19
Views
2K
  • · Replies 19 ·
Replies
19
Views
3K
Replies
1
Views
2K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K