Questions about stochastic calculus in Mathematica, for example how to use `ItoProcess` and `RandomFunction`.

learn more… | top users | synonyms

4
votes
1answer
50 views

How to use Tandem Queueing Network Process?

Mathematica provides QueueingNetworkProcess and QueueingProcess. However, I can't seem to figure out how to create a tandem ...
0
votes
1answer
53 views

Why does RandomFunction return variable number of data points?

I am using the following code snippet to generate Compound Poisson process random numbers: ...
1
vote
1answer
61 views

ItoProcess for coupled SDEs

I am trying to create an ItoProcess from the following system of SDEs: $\begin{bmatrix} \mathrm d x\\\mathrm d y\end{bmatrix} = \begin{bmatrix} 0 & 1\\ 0 & \theta\end{bmatrix} \begin{bmatrix} ...
2
votes
3answers
162 views

Hypergeometric function with a matrix argument

I am looking for the evaluation of a Hypergeometric function with a matrix argument as for example in Ploev and Edelman or as showcased in this Wikipedia article. From what I understand from ...
0
votes
0answers
49 views

How to deal with matrices involved in system of SDEs?

This question is in continuation of the the previous posts Solving Stochastic differential equation and Fast Simulations with Compile. What I want to do is numerically solving the epidemic model which ...
1
vote
0answers
57 views

Use of Ito's lemma in ItoProcess

In the documentation for the ItoProcess it says: Converting an ItoProcess to standard form automatically makes use of ...
0
votes
0answers
116 views

Stochastic Simulation using the Gillespie algorithm

I'm trying to reproduce the simulation with demographic stochasticity in Figure 1 from the paper entitled "Dynamical Resonance Can Account For Seasonality of Influenza Epidemics" ...
0
votes
1answer
98 views

How to define a new stochastic process which is function of another process?

I need to define a new process from for example Wiener process like $U(t)=f(W(t))$, (for example $f(x)=1+x^2$ ) and then calculate the average like $\langle U(t)U(s)\rangle$. Is it possible?
2
votes
1answer
112 views

Plot A Function Of A Stochastic Process

I am trying to do something very simple in Mathematica 9. I want to play around with option pricing and for that I thought it best to use the new stochastic process functionality. So, first of all I ...
2
votes
1answer
64 views

How to obtain SliceDistribution or StationaryDistribution for an ItoProcess when it is known to exist?

According to this reference page StationaryDistribution[proc] represents the stationary distribution of the process proc, when it exists. When I define the OrnsteinUhlenbeckProcess by the ...
2
votes
1answer
97 views

How to use some other driving process than the WienerProcess?

According to the following reference page http://reference.wolfram.com/language/ref/ItoProcess.html The driving process dproc can be any process that can be converted to a standard Ito process ...
2
votes
1answer
108 views

Expectation of GeometricBrownianMotionProcess

I am trying to compute $$\mathbb E\left[\max\left(\frac{S_{1/2}+S_1}{2}-K,0\right)\right]$$ where $K=100$ and $S_t$ is a geometric brownian motion (with $S_0=100$, drift $r=0.05$ and volatility ...
0
votes
0answers
109 views

Differential equation with random variable

How can I derive analytically or compute numerically the solution to following differential equation $$ dy/dt = y\cdot X\cdot (y\cdot X - g(y,X))\cdot X $$ where X is a random variable (e.g. from a ...
0
votes
0answers
151 views

Time series (stochastic process) estimating parameters using characteristic function

I have a time series of assets ${A_1, A_2, ..., A_n}$, which is described by a sophisticated distribution having the following characteristic function: $\phi(u; t;\theta)$, where $\theta$ is a vector ...
0
votes
0answers
101 views

Ito process estimating in Mathematica

How can I estimate parameters of Ito process in Mathematica? I have some time-process data (for example assets) and want it to be described Heston model - vector Ito process.
1
vote
1answer
80 views

Evaluating an ItoProcesss function

The following is my Stochastic D.E.: ...
0
votes
0answers
141 views

Approximating a stochastic integral with a Wiener process

I would like some assistance to solve the following: $\sum\limits_{i=0}^{n-1} e^{-k(n\Delta t -i\Delta t)}\Delta z_i$ where $z$ is a Wiener process, $\Delta z_i = z((i+1)\Delta t)-z(i\Delta t)$. ...
0
votes
0answers
176 views

Solving stochastic master equation

Can someone help please me in solving the following stochastic master equation for the density matrix $\rho$? where $\rho$ is the density matrix, $\sigma_i$ are the pauli matrices, $dt$ is the ...
0
votes
0answers
48 views

Determine appropriate initial conditions for exit time

These are some constants: ...
4
votes
2answers
970 views

Solving Stochastic differential equation

How do I solve the following simple stochastic differential equation: $$ m x''[t] + \Gamma x'[t] + k x[t] = \sqrt[]{2 kb T/\Gamma)} \eta[t] $$ here $\eta[t]$ is brownian motion , i.e. Wiener ...
1
vote
1answer
229 views

Finding the time at which an Ito process satisfies a constraint

I want to run an Ito stochastic process. I have the following parameters ...
1
vote
2answers
209 views

Problem in simulating discrete time stochastic

I have been playing with some stochastic questions and specially the problem here.It seems no matter for the first time in bet, gambler is going to lose the first bet.Am I right?How can we correct the ...
2
votes
1answer
229 views

Boundary condition for stochastic differential equation

I have a simple stochastic differential equation (SDE) with white noise: ...
3
votes
0answers
226 views

Malliavin Derivative with Mathematica is it possible?

Is it possible to define a Malliavin calculus with Mathematica 9? Consider a random variables on the Wiener-space $\Omega=\mathcal{C}([0,1])$ of the form ...
1
vote
0answers
166 views

Simulations with MonteCarlo and Autoregressive methods

I am trying to find the best-fit trend for my data. Here I just generated it but let's say I don't know my data trend at all. ...
5
votes
1answer
162 views

Error in ARIMAProcess example

I am trying to compile the sample of ARIMAProcess of MMA here .It doesn't work.What is wrong?Could you please help?I am using exact same code. ...
6
votes
3answers
401 views

Share experiences, preferably the surprising ones, with using ItoProcess

Can people please share their experiences, preferably the surprising ones, with using ItoProcess? I am a big fan of ItoProcess and have already used it for several finance-related tasks, though I ...
1
vote
1answer
493 views

ItoProcess function

While looking in the help manual of Mathematica concerning the ItoProcess function I found the following: ItoProcess[{a,b,c},x,t]: represents an Ito Process y(t)=c(t,x(t)), where ...
1
vote
0answers
61 views

Changing GeometricBrownianMotionProcess function

Since the function GeometricBrownianMotionProcess is given by Mathematica I have some technical questions. If we consider the following example: ...
8
votes
1answer
810 views

Mathematica code for hidden Markov models (HMM)

I am looking for some simple Mathematica code to model an HMM with just a few states and an equal number of observable signals (emissions). I am hoping to generate sample paths and keep track of the ...
2
votes
2answers
179 views

Stochastic problem

I have to organize a small sports league and I am puzzled on how to create the game plan. We are 8 persons playing table soccer with 2 vs. 2 matches. The idea is that each person plays once with ...
4
votes
0answers
430 views

ItoProcess for stochastic reaction-diffusion equation

I am trying to simulate a stochastic differential equation in time and space, but I'm unsure if this can be done in Mathematica. The sde that I would like to study is: $$ ...
11
votes
2answers
544 views

Efficient GeometricBrownianMotionProcess Monte Carlo simulation

Following the answers in this post, I'm trying to implement something similar. If the GBM stays inside the corridor [L, U] between predefined times it should return ...
6
votes
1answer
770 views

Monte Carlo simulation using geometric Brownian motion

I'm relatively new to Mathematica programming, so forgive my rather unsophisticated question: I'm trying to do a Monte Carlo simulation using geometric Brownian motion (GBM). I want to write a ...
15
votes
3answers
617 views

Plotting the solution of a vector stochastic differential equation

I have a vector stochastic differential equation, $$\mathrm dq = p\,\mathrm dt\qquad q(0)=0$$ $$\mathrm dp = (-q -p)\mathrm dt+\mathrm dW\qquad p(0)=10$$ This can be entered to give me the process ...