What is the solution to this ODE (and SDE)?

  • Thread starter Thread starter Only a Mirage
  • Start date Start date
  • Tags Tags
    Ode
Only a Mirage
Messages
57
Reaction score
0
I'm trying to analyze the following Ito stochastic differential equation:

$$dX_t = \|X_t\|dW_t$$

where X_t, dX_t, W_t, dW_t \in \mathbb{R}^n. Here, dW_t is the standard Wiener process and \|\bullet\| is the L^2 norm. I'm not sure if this has an analytical solution, but I am hoping to at least find an analytical expression for the expected value E[X_t].

In order to gain intuition for this problem, I'm considering the following ordinary differential equation:

$$\dot{z}(t) =\|z(t)\|b(t)$$

where z(t), b(t) \in \mathbb{R}^n and everything is completely deterministic. Does anyone know the analytical solution to this second equation, and under what conditions it exists?
 
Physics news on Phys.org
Only a Mirage said:
I'm trying to analyze the following Ito stochastic differential equation:

$$dX_t = \|X_t\|dW_t$$

where X_t, dX_t, W_t, dW_t \in \mathbb{R}^n. Here, dW_t is the standard Wiener process and \|\bullet\| is the L^2 norm. I'm not sure if this has an analytical solution, but I am hoping to at least find an analytical expression for the expected value E[X_t].

In order to gain intuition for this problem, I'm considering the following ordinary differential equation:

$$\dot{z}(t) =\|z(t)\|b(t)$$

where z(t), b(t) \in \mathbb{R}^n and everything is completely deterministic. Does anyone know the analytical solution to this second equation, and under what conditions it exists?

It is most easily solved in spherical polar coordinates where \|z\| = r so that z = r\mathbf{e}_r and \mathbf{b} = b_r\mathbf{e}_r + \mathbf{b}_n where \mathbf{b}_n \cdot \mathbf{e}_r = 0. We then have <br /> \dot r \mathbf{e}_r + r \dot{\mathbf{e}_r} = r(b_r\mathbf{e}_r + \mathbf{b}_n). Since \mathbf{e}_r and \dot{\mathbf{e}_r} are orthogonal we have <br /> \dot r = rb_r(t), \\<br /> \dot{\mathbf{e}_r} = \mathbf{b}_n.<br /> The radial and angular components thus decouple and the radial component has solution <br /> r(t) = r(0) \exp\left( \int_0^t b_r(s)\,ds\right).<br /> The angular component represents the motion of a point on the unit (n-1)-sphere. If \mathbf{b} is continuous then a solution should exist, but finding a co-ordinate representation of it valid for all time may be impossible, as one can see from the angular components in \mathbb{R}^3, <br /> \dot \theta = b_\theta(t), \\<br /> \dot \phi = \frac{b_\phi(t)}{\sin \theta(t)}<br /> with solution <br /> \theta(t) = \theta(0) + \int_0^t b_\theta(s)\,ds, \\<br /> \phi(t) = \phi(0) + \int_0^t \frac{b_\phi(s)}{\sin \theta(s)}\,ds, which becomes invalid if ever \theta(t) &lt; 0 or \theta(t) &gt; \pi.
 
  • Like
Likes 1 person
pasmith said:
It is most easily solved in spherical polar coordinates where \|z\| = r so that z = r\mathbf{e}_r and \mathbf{b} = b_r\mathbf{e}_r + \mathbf{b}_n where \mathbf{b}_n \cdot \mathbf{e}_r = 0. We then have <br /> \dot r \mathbf{e}_r + r \dot{\mathbf{e}_r} = r(b_r\mathbf{e}_r + \mathbf{b}_n). Since \mathbf{e}_r and \dot{\mathbf{e}_r} are orthogonal we have <br /> \dot r = rb_r(t), \\<br /> \dot{\mathbf{e}_r} = \mathbf{b}_n.<br /> The radial and angular components thus decouple and the radial component has solution <br /> r(t) = r(0) \exp\left( \int_0^t b_r(s)\,ds\right).<br /> The angular component represents the motion of a point on the unit (n-1)-sphere. If \mathbf{b} is continuous then a solution should exist, but finding a co-ordinate representation of it valid for all time may be impossible, as one can see from the angular components in \mathbb{R}^3, <br /> \dot \theta = b_\theta(t), \\<br /> \dot \phi = \frac{b_\phi(t)}{\sin \theta(t)}<br /> with solution <br /> \theta(t) = \theta(0) + \int_0^t b_\theta(s)\,ds, \\<br /> \phi(t) = \phi(0) + \int_0^t \frac{b_\phi(s)}{\sin \theta(s)}\,ds, which becomes invalid if ever \theta(t) &lt; 0 or \theta(t) &gt; \pi.

Thanks a lot for the answer. Can you explain why <br /> \theta(t) = \theta(0) + \int_0^t b_\theta(s)\,ds, \\<br /> \phi(t) = \phi(0) + \int_0^t \frac{b_\phi(s)}{\sin \theta(s)}\,ds, becomes invalid if
ever \theta(t) &lt; 0 or \theta(t) &gt; \pi?
 
Thread 'Direction Fields and Isoclines'
I sketched the isoclines for $$ m=-1,0,1,2 $$. Since both $$ \frac{dy}{dx} $$ and $$ D_{y} \frac{dy}{dx} $$ are continuous on the square region R defined by $$ -4\leq x \leq 4, -4 \leq y \leq 4 $$ the existence and uniqueness theorem guarantees that if we pick a point in the interior that lies on an isocline there will be a unique differentiable function (solution) passing through that point. I understand that a solution exists but I unsure how to actually sketch it. For example, consider a...
Back
Top