Stationary distributions and reversibility of jump processes

[27/89%]
Show
LaTeX

Find the stationary distributions when they exist for the processes in Problem 15.3.

?

Show that for any finite state space irreducible CTMC {X(t):t≥0}\left\{ X(t): t \geq 0\right\} with all λi∈(0,∞)\lambda_{i} \in (0, \infty ), there is a unique stationary distribution.

?

This is a birth and death process such that αn≡α,βn=nβ,n≥0,0<α,β<∞\alpha_{n} \equiv \alpha , \beta_{n} = n \beta , n \geq 0,0 < \alpha , \beta < \infty. Show that this process has a stationary distribution that is Poisson with mean ρ=λμ\rho = \frac{\lambda }{\mu }.

?

Consider the Markov process for the M/M/1 queue, as given in Figure 7.4.

?
(a)

Find the steady-state process probabilities (as a function of ρ=λ/μ\rho = \lambda /\mu) from (7.15) and also as the solution to (7.23). Verify that the two solutions are the same.

(b)

For the remaining parts of the exercise, assume that ρ=0.01\rho = 0.01, thus ensuring (for aiding intuition) that states 00 and 11 are much more probable than the other states. Assume that the process has been running for a very long time and is in steady state. Explain in your own words the difference between π1\pi_{1} (the steady-state probability of state 11 in the embedded chain) and p1p_{1} (the steady-state probability that the process is in state 11). More explicitly, what experiments could you perform (repeatedly) on the process to measure π1\pi_{1} and p1p_{1}.

(c)

Now suppose you want to start the process in steady state. Show that it is impossible to choose initial probabilities so that both the process and the embedded chain start in steady state. Which version of steady state is closest to your intuitive view? (There is no correct answer here, but it is important to realize that the notion of steady state is not quite as simple as you might imagine.)

(d)

Let M(t)M(t) be the number of transitions (counting both arrivals and departures) that take place by time tt in this Markov process and assume that the embedded Markov chain starts in steady state at time 00. Let U1,U2,…U_{1}, U_{2}, \ldots be the sequence of holding intervals between transitions (with U1U_{1} being the time to the first transition). Show that these rv's are identically distributed. Show by example that they are not independent (i.e., M(t)M(t) is not a renewal process).

Consider the Markov process illustrated below. The transitions are labeled by the rate qijq_{ij} at which those transitions occur. The process can be viewed as a single-server queue where arrivals become increasingly discouraged as the queue lengthens. The word time average below refers to the limiting time average over each sample path of the process, except for a set of sample paths of probability 00.  {#fig-1 width="70%"}

?
(a)

Find the time-average fraction of time pip_{i} spent in each state i>0i > 0 in terms of p0p_{0} and then solve for p0p_{0}.

(b)

Find a closed form solution to ∑ipiνi\sum_{i}p_{i}\nu_{i}, where νi\nu_{i} is the rate at which transitions out of state ii occur. Show that the embedded chain is positive recurrent for all choices of λ>0\lambda > 0 and μ>0\mu > 0 and explain intuitively why this must be so.

(c)

For the embedded Markov chain corresponding to this process, find the steady-state probabilities πi\pi_{i} for each i≥0i \geq 0 and the transition probabilities PijP_{ij} for each i,ji, j.

(d)

For each ii, find both the time-average interval and the time-average number of overall state transitions between successive visits to ii.

Remark.
?

Hint for (a): First find an equation relating pip_{i} to pi+1p_{i+1} for each ii. It also may help to recall the power series expansion of exe^{x}.

Consider the three-state Markov process below; the number given on edge (i,j)(i,j) is qijq_{ij}, the transition rate from ii to jj. Assume that the process is in steady state.  {#fig-1 width="45%"}

?
(a)

Is this process reversible?

(b)

Find pip_{i}, the time-average fraction of time spent in state ii for each ii.

(c)

Given that the process is in state ii at time tt, find the mean delay from tt until the process leaves state ii.

(d)

Find πi\pi_{i}, the time-average fraction of all transitions that go into state ii for each ii.

(e)

Suppose the process is in steady state at time tt. Find the steady-state probability that the next state to be entered is state 11.

(f)

Given that the process is in state 11 at time tt, find the mean delay until the process first returns to state 11.

(g)

Consider an arbitrary irreducible finite-state Markov process in which qij=qjiq_{ij} = q_{ji} for all i,ji, j. Either show that such a process is reversible or find a counter example.

Consider a preemptive resume LCFS queueing system with two classes of customers. Type AA customer arrivals are Poisson with rate λA\lambda_{A} and Type BB customer arrivals are Poisson with rate λB\lambda_{B}. The service time for type AA customers is exponential with rate μA\mu_{A} and that for type BB is exponential with rate μB\mu_{B}. Each service time is independent of all other service times and of all arrival epochs.

Define the `state' of the system at time tt by the string of customer types in the system at tt, in order of arrival. Thus state ABAB means that the system contains two customers, one of type AA and the other of type BB; the type BB customer arrived later, so is in service. The set of possible states arising from transitions out of ABAB is as follows:

ABAABA if another type AA arrives.

ABBABB if another type BB arrives.

AA if the customer in service (BB) departs.

Note that whenever a customer completes service, the next most recently arrived resumes service, so the state changes by eliminating the final element in the string.

?
(a)

Draw a graph for the states of the process, showing all states with two or fewer customers and a couple of states with three customers (label the empty state as EE). Draw an arrow, labeled by the rate, for each state transition. Explain why these are states in a Markov process.

(b)

Is this process reversible? Explain. Assume positive recurrence.

(c)

Characterize the process of type AA departures from the system (i.e., are they Poisson, do they form a renewal process, at what rate do they depart, etc.?)

(d)

Express the steady-state probability P{A}\mathbb {P}\left\{ A\right\} of state AA in terms of the probability of the empty state P{E}\mathbb {P}\left\{ E\right\}. Find the probability P{AB}\mathbb {P}\left\{ AB\right\} and the probability P{ABBA}\mathbb {P}\left\{ ABBA\right\} in terms of P{E}\mathbb {P}\left\{ E\right\}. Use the notation ρA=λA/μA\rho_{A} = \lambda_{A}/\mu_{A} and ρB=λB/μB\rho_{B} = \lambda_{B}/\mu_{B}.

(e)

Let QnQ_{n} be the probability of nn customers in the system, as a function of Q0=P{E}Q_{0} = \mathbb {P}\left\{ E\right\}. Show that Qn=(1−ρ)ρnQ_{n} = (1-\rho )\rho^{n} where ρ=ρA+ρB\rho = \rho_{A}+\rho_{B}.

Remark.
?

Hint for (b): If there is a transition from one state SS to another state S′S', how is the number of transitions from SS to S′S' related to the number from S′S' to SS?

?
(a)

Generalize Exercise 7.22 to the case in which there are mm types of customers, each with independent Poisson arrivals and each with independent exponential service times. Let λi\lambda_{i} and μi\mu_{i} be the arrival rate and service rate respectively of the iith user. Let ρi=λi/μi\rho_{i} = \lambda_{i}/\mu_{i} and assume that ρ=ρ1+ρ2+⋯+ρm<1\rho = \rho_{1}+\rho_{2}+\cdots +\rho_{m} < 1. In particular, show, as before that the probability of nn customers in the system is Qn=(1−ρ)ρnQ_{n} = (1-\rho )\rho^{n} for 0≤n<∞0 \leq n < \infty.

(b)

View the customers in (a) as a single type of customer with Poisson arrivals of rate λ=∑iλi\lambda = \sum_{i}\lambda_{i} and with a service density ∑i(λi/λ)μiexp⁡(−μix)\sum_{i}(\lambda_{i}/\lambda )\mu_{i}\exp (-\mu_{i}x). Show that the expected service time is ρ/λ\rho /\lambda. Note that you have shown that, if a service distribution can be represented as a weighted sum of exponentials, then the distribution of customers in the system for LCFS service is the same as for the M/M/1 queue with equal mean service time.

Suppose a Markov process with the set of states AA is reversible and has steady-state probabilities pip_{i}, i∈Ai \in A. Let BB be a subset of AA and assume that the process is changed by setting qij=0q_{ij} = 0 for all i∈B,j∉Bi \in B, j \notin B. Assuming that the new process (starting in BB and remaining in BB) is irreducible, show that the new process is reversible and find its steady-state probabilities.

?

Consider a queueing system with two classes of customers. Type AA customer arrivals are Poisson with rate λA\lambda_{A} and type BB customer arrivals are Poisson with rate λB\lambda_{B}. The service time for type AA customers is exponential with rate μA\mu_{A} and that for type BB is exponential with rate μB\mu_{B}. Each service time is independent of all other service times and of all arrival epochs.

?
(a)

First assume there are infinitely many identical servers, and each new arrival immediately enters an idle server and begins service. Let the state of the system be (i,j)(i,j), where ii and jj are the numbers of type AA and type BB customers respectively in service. Draw a graph of the state transitions for i≤2,j≤2i \leq 2, j \leq 2. Find the steady-state PMF, {p(i,j);i,j≥0}\left\{ p(i,j); i,j \geq 0\right\}, for the Markov process.

(b)

Assume for the rest of the exercise that there is some finite number mm of servers. Customers who arrive when all servers are occupied are turned away. Find the steady-state PMF, {p(i,j);i,j≥0,i+j≤m}\left\{ p(i,j); i,j \geq 0, i+j \leq m\right\}, in terms of p(0,0)p(0,0) for this Markov process.

(c)

Let QnQ_{n} be the probability that there are nn customers in service at some given time in steady state. Show that Qn=p(0,0)ρn/n!Q_{n} = p(0,0)\rho_{n}/n! for 0≤n≤m0 \leq n \leq m, where ρ=ρA+ρB\rho = \rho_{A}+\rho_{B}, ρA=λA/μA\rho_{A} = \lambda_{A}/\mu_{A}, and ρB=λB/μB\rho_{B} = \lambda_{B}/\mu_{B}. Solve for p(0,0)p(0,0).

Remark.
?

Hint for (a): Note that the type AA and type BB customers do not interact.

Hint for (b): Combine (a) with the result of Exercise 7.25.

?
(a)

Generalize Exercise 7.26 to the case in which there are KK types of customers, each with independent Poisson arrivals and each with independent exponential service times. Let λk\lambda_{k} and μk\mu_{k} be the arrival rate and service rate respectively for the kkth user type, 1≤k≤K1 \leq k \leq K. Let ρk=λk/μk\rho_{k} = \lambda_{k}/\mu_{k} and ρ=ρ1+ρ2+⋯+ρK\rho = \rho_{1}+\rho_{2}+\cdots +\rho_{K}. In particular, show, as before, that the probability of nn customers in the system is Qn=p(0,…,0)ρn/n!Q_{n} = p(0,\ldots ,0)\rho^{n}/n! for 0≤n≤m0 \leq n \leq m.

(b)

View the customers in (a) as a single type of customer with Poisson arrivals of rate λ=∑kλk\lambda = \sum_{k}\lambda_{k} and with a service density ∑k(λk/λ)μkexp⁡(−μkx)\sum_{k}(\lambda_{k}/\lambda )\mu_{k}\exp (-\mu_{k}x). Show that the expected service time is ρ/λ\rho /\lambda. Note that what you have shown is that, if a service distribution can be represented as a weighted sum of exponentials, then the distribution of customers in the system is the same as for the M/M/m/m queue with equal mean service time.

If X(t)X(t) is an irreducible positive recurrent Markov chain, and u(⋅)u(\cdot ) is a bounded function on the integers, show that

1t∫0tu(X(s))ds→ a.s. ∑i∈Sπiu(i) \frac{1}{t} \int _{0}^{t} u(X(s)) d s \xrightarrow {\text{ a.s. }} \sum _{i \in S} \pi _{i} u(i)

where π\mathbf{\pi } is the stationary distribution of X(t)X(t).

?

Let XX be an irreducible, continuous-time Markov chain with generator G\mathbf{G} on the state space T=S∪{∞}T = S \cup \left\{ \infty \right\}, where SS is countable and non-empty. Show that a distribution π\mathbf{\pi } on TT satisfies πG=0\mathbf{\pi } \mathbf{G} = \mathbf{0} if and only if, for j∈S,∑i∈Tπigij=0j \in S, \sum_{i \in T} \pi_{i} g_{i j} = 0.

?

Let YY be an irreducible discrete-time Markov chain on a countably infinite state space SS, having transition matrix Y=(yij)\mathbf{Y} = \left(y_{i j}\right) satisfying yii=0y_{i i} = 0 for all states ii, and with stationary distribution ν\mathbf{\nu }. Construct a continuous-time process XX on SS for which YY is the jump chain, such that XX has no stationary distribution.

?

Consider an immigration-death process XX, being a birth-death process with birth rates λn=λ\lambda_{n} = \lambda and death rates μn=nμ\mu_{n} = n \mu. Find the transition matrix of the jump chain YY, and show that it has as stationary distribution

πn=12(n!)(1+nρ)ρne−ρ \pi _{n} = \frac{1}{2(n!)}\left(1+\frac{n}{\rho }\right) \rho ^{n} e^{-\rho }

where ρ=λ/μ\rho = \lambda / \mu. Explain why this differs from the stationary distribution of XX.

?

In a model for a biofilm population, we assume there are nn colonizable 'niches' (or 'food sources'). Let X(t)X(t) be the number of occupied niches at time tt, and assume XX is a Markov chain that evolves as follows. The lifetime of any colony is exponentially distributed with parameter μ\mu; if X(t)=iX(t) = i, the rate of establishment of a new colony in an empty niche is λi(n−i)\lambda i(n-i). The usual independence may be assumed.

For X(0)=1,2,…X(0) = 1,2, \ldots, find the mean time until the population is extinct, which is to say that no niches are occupied. Discuss the implications when nn is large.

?
?
(a)

Let T>0T > 0 and let X={X(t):0≤t≤T}X = \left\{ X(t): 0 \leq t \leq T\right\} be an irreducible, non-explosive Markov chain with stationary distribution π\pi, and suppose that X(0)X(0) has distribution π\pi. Let Y(t)=X(T−t)Y(t) = X(T-t) for 0≤t≤T0 \leq t \leq T. We call XX reversible (in equilibrium) if XX and YY have the same joint distributions.

(i) Show that YY is a (left-continuous) Markov chain with transition probabilities p^ij(t)=(πj/πi)pji(t)\widehat{p}_{i j}(t) = \left(\pi_{j} / \pi_{i}\right) p_{j i}(t) and generator G^\widehat{\mathbf{G}} satisfying πjg^ji=πigij\pi_{j} \hat{g}_{j i} = \pi_{i} g_{i j}, where the pji(t)p_{j i}(t) and G=(gij)\mathbf{G} = \left(g_{i j}\right) are those of XX. Show that YY is irreducible and non-explosive with stationary distribution π\pi.

(ii) Show that XX is reversible in equilibrium if and only if the detailed balance equations πigij=πjgji\pi_{i} g_{i j} = \pi_{j} g_{j i} (for all ii and jj) hold.

(iii) Show that a measure ν\mathbf{\nu } satisfies νG=0\mathbf{\nu } \mathbf{G} = \mathbf{0} if it satisfies the detailed balance equations.

(b)

Let XX be irreducible and non-explosive with stationary distribution π\pi, and assume X(0)X(0) has distribution π\pi.

(i) Kolmogorov’s criterion. Show that XX is reversible if and only if, for all nn and all finite sequences k1,k2,…,knk_{1}, k_{2}, \ldots , k_{n} of states,

gk1,k2gk2,k3⋯gkn−1,kngkn,k1=gk1,kngkn,kn−1⋯gk2,k1 g_{k_{1}, k_{2}} g_{k_{2}, k_{3}} \cdots g_{k_{n-1}, k_{n}} g_{k_{n}, k_{1}} = g_{k_{1}, k_{n}} g_{k_{n}, k_{n-1}} \cdots g_{k_{2}, k_{1}}

(ii) Kelly’s criterion. Show that XX is reversible if, for all distinct triples i,j,k∈Si, j, k \in S, we have gijgjkgki=gikgkjgjig_{i j} g_{j k} g_{k i} = g_{i k} g_{k j} g_{j i}, and in addition there exists c∈Sc \in S such that gic>0g_{i c} > 0 for all i≠ci \neq c.

(c)

Show that every irreducible chain XX with exactly two states is reversible in equilibrium.

(d)

Show that every non-explosive birth-death process XX having a stationary distribution is reversible in equilibrium.

Let XX be an irreducible continuous-time Markov chain on the state space SS with transition probabilities pjk(t)p_{j k}(t) and unique stationary distribution π\mathbf{\pi }, and write P(X(t)=j)=aj(t)\mathbb {P}\left(X(t) = j\right) = a_{j}(t). If c(x)c(x) is a concave function, show that the function d(t)=∑j∈Sπjc(aj(t)/πj)d(t) = \sum_{j \in S} \pi_{j} c\left(a_{j}(t) / \pi_{j}\right) increases to c(1)c(1) as t→∞t \rightarrow \infty.

The relative entropy (or Kullback-Leibler divergence) of two strictly positive probability mass functions f,gf, g on a subset SS of integers is defined as

D(f;g)=∑i∈Sf(i)log⁡(f(i)/g(i)) D(f ; g) = \sum _{i \in S} f(i) \log (f(i) / g(i))

Prove that, if XX has a finite state space and stationary distribution π\mathbf{\pi }, the relative entropy D(a(t);π)D(a(t) ; \mathbf{\pi }) decreases monotonely to 0 as t→∞t \rightarrow \infty.

?

With the notation of the preceding problem, let uk(t)=P(X(t)=k∣X(0)=0)u_{k}(t) = \mathbb {P}\left(X(t) = k \mid X(0) = 0\right), and suppose the chain is reversible in equilibrium (see Problem (6.15.16)). Show that u0(2t)=∑j(π0/πj)uj(t)2u_{0}(2 t) = \sum_{j}\left(\pi_{0} / \pi_{j}\right) u_{j}(t)^{2}, and deduce that u0(t)u_{0}(t) decreases to π0\pi_{0} as t→∞t \rightarrow \infty.

?

Let XX and YY be independent recurrent birth-death processes with the same parameters (and no explosions). It is not assumed that X0=Y0X_{0} = Y_{0}. Show that:

?
(a)

for any A⊆R,∣P(Xt∈A)−P(Yt∈A)∣→0A \subseteq \mathbb {R},\left|\mathbb {P}\left(X_{t} \in A\right)-\mathbb {P}\left(Y_{t} \in A\right)\right| \rightarrow 0 as t→∞t \rightarrow \infty,

(b)

if P(X0≤Y0)=1\mathbb {P}\left(X_{0} \leq Y_{0}\right) = 1, then E[g(Xt)]≤E[g(Yt)]\mathbb {E}\left[g\left(X_{t}\right)\right] \leq \mathbb {E}\left[g\left(Y_{t}\right)\right] for any increasing function gg.

Let X={X(t):t≥0}X = \left\{ X(t): t \geq 0\right\} be a Markov chain having stationary distribution π\mathbf{\pi }. We may sample XX at the times of a Poisson process: let NN be a Poisson process with intensity λ\lambda, independent of XX, and define Yn=X(Tn)Y_{n} = X\left(T_{n}\right). Show that Y={Yn:n≥0}Y = \left\{ Y_{n}: n \geq 0\right\} is a discrete-time Markov chain with the same stationary distribution as XX. (This exemplifies the 'Pasta' property: Poisson arrivals see time averages.) [The full assumption of the independence of NN and XX is not necessary for the conclusion. It suffices that {N(s):s≥t}\left\{ N(s): s \geq t\right\} be independent of {X(s):s≤t}\left\{ X(s): s \leq t\right\}, a property known as 'lack of anticipation'. It is not even necessary that XX be Markov; the Pasta property holds for many suitable ergodic processes.]

?

Ergodic theorem in continuous time. Let XX be an irreducible recurrent continuous-time Markov chain with generator G\mathbf{G} and finite mean return times mjm_{j}.

?
(a)

Show that 1t∫0tI{X(s)=j}ds→ a.s. 1mjgj\frac{1}{t} \int_{0}^{t} I_{\left\{ X(s) = j\right\} } d s \xrightarrow {\text{ a.s. }} \frac{1}{m_{j} g_{j}} as t→∞t \rightarrow \infty;

(b)

deduce that the stationary distribution π\pi satisfies πj=1/(mjgj)\pi_{j} = 1 /\left(m_{j} g_{j}\right);

(c)

show that, if ff is a bounded function on SS,

1t∫0tf(X(s))ds→ a.s. ∑iπif(i) as t→∞. \frac{1}{t} \int _{0}^{t} f(X(s)) d s \xrightarrow {\text{ a.s. }} \sum _{i} \pi _{i} f(i) \quad \text{ as } t \rightarrow \infty .

The dynamics of a system is given by a jump Markov process with transition intensities

Q=(−3123−4142−6). Q = \begin{pmatrix} -3 & 1 & 2 \\ 3 & -4 & 1 \\ 4 & 2 & -6 \end{pmatrix}.

Show that all states of the process communicate, and find its stationary distribution.

?

Let ξ(t)\xi (t) be an integer-valued continuous-time random walk with a reflecting barrier at state 0, i.e., a jump Markov process with transition intensities qi,i+1=λ>0q_{i,i+1} = \lambda > 0 for i≥0i \geq 0, qi,i−1=μ>0q_{i,i-1} = \mu > 0 for i≥1i \geq 1, qii=−(λ+μ)q_{ii} = -(\lambda +\mu ) for i≥1i \geq 1, q00=−λq_{00} = -\lambda, qij=0q_{ij} = 0 for all other j≠ij \neq i. Show that a stationary distribution of the process ξ(t)\xi (t) exists if and only if λ<μ\lambda < \mu, and find it in this case.

?

A homogeneous jump Markov process ξ(t)\xi (t) with values 1, 2, 3, 4, 5, and 6 is given by the transition intensity matrix

Q=(−2101001−4110100−4040000−2020040−4000200−2). Q = \begin{pmatrix} -2 & 1 & 0 & 1 & 0 & 0 \\ 1 & -4 & 1 & 1 & 0 & 1 \\ 0 & 0 & -4 & 0 & 4 & 0 \\ 0 & 0 & 0 & -2 & 0 & 2 \\ 0 & 0 & 4 & 0 & -4 & 0 \\ 0 & 0 & 2 & 0 & 0 & -2 \end{pmatrix}.
?
(a)

Find the transition probability matrix of the embedded Markov chain ξn∗\xi_{n}^{*}. Use it to classify the states and find the absorption probabilities for each of the essential classes.

(b)

Find all stationary distributions of the process ξ(t)\xi (t) and of the embedded Markov chain ξn∗\xi_{n}^{*}. Do they coincide? Why?

(c)

Find the limiting distributions for the Markov process and for the Markov chain, if ξ(0)=2\xi (0) = 2.

(d)

Find the limiting distributions for the Markov process and for the Markov chain, if ξ(0)=3\xi (0) = 3.

The transition probabilities of the Markov chain (Xt,t⩾0)\left(X_{t}, t \geqslant 0\right) with phase space X={1,2,3}\mathcal{X} = \left\{ 1,2,3\right\} have the form:

p11(h)=1−λh+o(h),p12(h)=λh+o(h),p13(h)=o(h)p21(h)=o(h),p22(h)=1−μh+o(h),p23(h)=μh+o(h)p31(h)=νh+o(h),p32(h)=o(h),p33(h)=1−νh+o(h), \begin{align} p_{11}(h)& =1-\lambda h+o(h), & p_{12}(h)& =\lambda h+o(h), & p_{13}(h)& =o(h) \\ p_{21}(h)& =o(h), & p_{22}(h)& = 1 - \mu h + o(h), & p_{23}(h)& = \mu h + o(h) \\ p_{31}(h)& = \nu h + o(h), & p_{32}(h)& = o(h), & p_{33}(h)& =1 - \nu h + o(h), \end{align}

as h→0+h \to 0+. Prove that such a chain satisfies the condition of the ergodic theorem. Find its infinitesimal matrix and stationary distribution.

?

Find all invariant measures corresponding to the transition probability matrix

P(t)=(Pij(t))=[(2+e−3t)/3(1−e−3t)/3(2−2e−3t)/3(1+2e−3t)/3] P(t) = \left(P_{ij}(t)\right) =\left[\begin{smallmatrix} \left(2+e^{-3 t}\right) / 3 & \left(1-e^{-3 t}\right) / 3 \\ \left(2-2 e^{-3 t}\right) / 3 & \left(1+2 e^{-3 t}\right) / 3 \end{smallmatrix}\right]
?

References

[AtLa06] Athreya K.B., Lahiri S.N. Measure Theory and Probability Theory. Springer, 2006.
[Gal13] Gallager R.G. Stochastic Processes: Theory for Applications. Cambridge University Press, 2013.
[GrSt20] Grimmett G.R., Stirzaker D.R. Probability and Random Processes. 4th ed., Oxford University Press, 2020.
[КоФо03] Коршунов Д.А., Фосс С.Г. Сборник задач и упражнений по теории вероятностей. 2-е изд., Новосибирский государственный университет, 2003.
[ЖРШ16] Жуковский М.Е., Родионов И.В., Шабанов Д.А. Основы теории случайных процессов. МФТИ, 2016.
[ПУУ86] Прохоров А.В., Ушаков В.Г., Ушаков Н. Задачи по теории вероятностей. Наука. Главная редакция физико-математической литературы, 1986.