For questions regarding the numerical treatment of processes whose behaviors are determined by both deterministic (predictable) and non-deterministic (random) actions.

learn more… | top users | synonyms

2
votes
2answers
58 views

Extracting time scales information from empirical cumulative distribution function

I have a stochastic process (a Markov chain actually) that has two absorbing states. I am using a difference equation to calculate the first passage time to either of the absorbing states. There are ...
2
votes
2answers
177 views

visualization of 3D probability flow

I have a master equation for $P(N_A^+,N_B^+,N_C^+,t)$, with $N_A^+,N_B^+,N_C^+$ all discrete. The numerical integration is done by this Matlab program using Euler's method. Despite the crude Euler's ...
3
votes
1answer
168 views

Slight mistake in Stochastic Galerkin code

I'm following Paul Constantine's Primer on Stochastic Galerkin Method, Section 3.1 (2D Poisson Example). In this matlab code, the example attempts to solve the PDE $$\alpha(w)(u_{xx}+u_{yy})=1 \text{ ...
3
votes
0answers
49 views

Stochastic Collocation for time evolving ODE

For an Stochastic Differential Equation, e.g., $$ \frac{du}{dt} = \alpha*\sin(u*t) $$ where $\alpha$ is normally distributed with nonzero mean, I am trying to use a stochastic collocation approach ...
4
votes
1answer
231 views

Examples of numerical solution of stochastic differential equation(SDE)?

I want to simulate a nonlinear stochastic differential equation $$ {\rm d}X_t = f(X_t) {\rm d}t + g(X_t){\rm d}B_t $$ where $f,g \in C^{\infty}({\mathbb R}^n ,{\mathbb R})$ and $B_t$ is ...
10
votes
2answers
146 views

Numerical method for equation solving that works on stochastically computed functions

There are many well known numerical methods for solving equations of the type $$ f(x) = 0, \quad x \in \mathbb{R}^n,$$ e.g. bisection method, Newton's method, etc. In my application $f(x)$ is ...
4
votes
2answers
124 views

Convergence tests in Markov Chain Monte Carlo

For a relatively simple Markov chain Monte Carlo process, such as using Metropolis to find calculate thermal averages for an Ising model, how is it possible to determine whether quantities have ...
2
votes
1answer
149 views

How to solve this numerical technique problem?

Well, in a numerical technique test we were given the following problem: A physical phenomenon is modeled such that, $F(f,d) = A(f)/d^2 + L$; Where, $F$ is a function of frequency $f$ and distance ...
5
votes
0answers
215 views

Stochastic Galerkin Projection Approach for using Generalized Polynomial Chaos Expansion (GPCE) in solving PDE

I am not sure if this is very general question but I want to know Is there any way that I can define the Test and trial function in the way that I want and dont use the default functions. so if I want ...
3
votes
2answers
326 views

Differences between “least square”, “mean square” and “least mean square”?

I was wondering what differences are between the terminology: "least square (LS)" "mean square (MS)" and "least mean square (LMS)"? I get confused when reading in Spall's Introduction to Stochastic ...
0
votes
3answers
140 views

Translating from SBML to Petri Net (matrix representation)

Is there a piece of software that can easily take a representation of a biochemical network from a SBML file and translate it into a Petri Net? Specifically, I'm looking for something that can give ...
3
votes
1answer
94 views

Introduction for (numerical) linear algebra of random variables

I am in search of an introduction into numerical linear algebra - or, at least, pure linear algebra - that treats the case when the input data are random variables. A typical application would be to ...
2
votes
1answer
48 views

Optimal sample size for Stochastic Steepest Descent

Suppose $g(x_{1:n})$ is the estimate of a gradient, which is calculated at each step of a Stochastic Steepest Descent algorithm. A dataset $x_{1:n}$ is simulated at each step, so if $n$ is small the ...
5
votes
4answers
671 views

Simulated Annealing proof of convergence

I implemented downhill simplex simulated annealing algorithm. Algorithm is very hard to tune, w.r.t. parameters including cooling schedule, starting temperature... My first question is about ...
3
votes
1answer
185 views

Where can I find coded examples of stochastic collocation applied to an elliptical PDE using smolyak sampling?

I'm having some troubles implementing a collocation method to solve a stochastic partial differential equation of the form: $\nabla (a(x,w)\nabla u(x,w))=f(x,w)$ in $D$, $u=g$ in $\partial D$ where ...
3
votes
1answer
146 views

Sampling strategies to solve a stochastic partial differential equation

Suppose I had a stochastic partial differential equation of the form: $\nabla^2U=F(x,D)$, where $x\in\Omega\equiv [0,1]$ and $F(x,D)$ is a function which depends on position $x$ and a uniform random ...
3
votes
3answers
2k views

How to choose a good step size for stochastic gradient descent?

For the purpose of model fitting in a large time series dataset, I am using stochastic gradient descent of the negative log likelihood. The model is nonlinear and non-convex. Is there a thumb rule for ...
5
votes
2answers
282 views

What is the deterministic counterpart of Robbins-Monro algorithm?

From Wikipedia, assume that we have a function $M(x)$, and we want to solve the equation $M(x) = 0$. But we cannot directly observe the function $M(x)$, we can instead obtain measurements of the ...
4
votes
1answer
178 views

Convert ODE into discrete probabilistic model

how can I turn an ODE equation into a discrete probabilistic model? I take for example the Verhulst equation for the growth of a population. $$\frac{dP}{dt} = rP(1-P/K)$$ I was thinking to simulate ...
4
votes
0answers
125 views

Richardson extrapolation for strong rate of convergence of SDE

Is it possible to apply Richardson extrapolation with Euler-Maruyama scheme to improve strong rate of convergence of stochastic differential equations?