LiZnB的博客

Random Processes

发布于 # 概率

1. Geometric random variable

P(X=k)=(1p)k1p,k=1,2,3,\begin{equation*} P(X=k)=(1-p)^{k-1} p, \quad k=1,2,3, \cdots \tag{1} \end{equation*}

Mean of X=E(X)=1pX=E(X)=\frac{1}{p}

Memoryless Property (or Markovian property)

P(X=m+nX>m)=P(X=n),m,n1\begin{equation*} P(X=m+n \mid X>m)=P(X=n), \quad m, n \geq 1 \tag{2} \end{equation*}

2. Exponential Random Variable

The cumulative distribution function (cds)

FX(x)=P(Xx)={1eλx,0x<0,x<0(3)F_{X}(x)=P(X \leq x)= \begin{cases}1-e^{-\lambda x}, & 0 \leq x<\infty \tag{3}\\ 0, & x<0\end{cases} P(X>x)=1P(Xx)=eλx,x0.P(X>x)=1-P(X \leq x)=e^{-\lambda x}, \quad x \geq 0 .

The probability density function (pdf)

fX(x)={λeλx,x>00,x0(4)f_{X}(x)= \begin{cases}\lambda e^{-\lambda x}, & x>0 \tag{4}\\ 0, & x \leq 0\end{cases}

Mean of X=E(X)=1λX=E(X)=\frac{1}{\lambda}

3. Stochastic Processes

A stochastic process is a collection of random variables {X(t):tT}\{X(t): t \in T\}. Note that X(t)X(t) is a r.v. for each tTt \in T. TT is called the index set of the process, a tTt \in T is called an index. Very often, the index tt is interpreted as time.

If TT is a countable set, then the stochastic process is called a discrete time process.

If TT is an interval of R\mathbb{R}, then the stochastic process is called a continuous time process.

The set of all values that X(t)X(t) may assume is called the state space, denoted by SS. SS can be a discrete state space (countable). In this case, the stochastic process is called a chain. Eg, a machine has only 2 states, working or failed. SS can be a continuous state space. Eg, X(t)X(t) represents the distance covered by a robot arm in time tt.

There are 4 different types of stochastic processes:

  1. discrete time, discrete state space processes (discrete time chain)
  2. discrete time, continuous state space processes
  3. continuous time, discrete state space processes (continuous time chain)
  4. continuous time, continuous state space processes

4. Poisson Process

Poisson random variable

P(X(t)=k)=eλt(λt)kk!,k=0,1,2,\begin{equation*} P(X(t)=k)=\frac{e^{-\lambda t}(\lambda t)^{k}}{k!}, \quad k=0,1,2, \cdots \tag{5} \end{equation*}

A collection of Poisson random variables {X(t)\{X(t) : t0}t \geq 0\} is a Poisson process. λ\lambda is a parameter and is called the rate of the Poisson process.

Superposition of Poisson Processes

λ=Σi=1nλi\lambda = \Sigma_{i=1}^{n} \lambda_{i}

Decomposition of a Poisson Process

5. Discrete Time Markov Chain (DTMC)

Markov property

P(Xn=jXn1=i,Xn2=i2,Xn3=i3,,X0=in)P\left(X_{n}=j \mid X_{n-1}=i, X_{n-2}=i_{2}, X_{n-3}=i_{3}, \cdots, X_{0}=i_{n}\right)

=P(Xn=jXn1=i)\begin{equation*} =P\left(X_{n}=j \mid X_{n-1}=i\right) \tag{6} \end{equation*}

One-Step Transition Probability

pijpij(1)=P(Xn=jXn1=i),n1\begin{equation*} p_{i j} \equiv p_{i j}(1)=P\left(X_{n}=j \mid X_{n-1}=i\right), \quad n \geq 1 \tag{7} \end{equation*}

Transition Probability Matrix (TPM)

P=[pij]=[p00p01p02p10p11p12p20p21p22]0pij1,i,jNjpij=1,iN\begin{align*} P=\left[p_{i j}\right] & =\left[\begin{array}{cccc} p_{00} & p_{01} & p_{02} & \cdots \\ p_{10} & p_{11} & p_{12} & \cdots \\ p_{20} & p_{21} & p_{22} & \cdots \\ \cdot & \cdot & \cdot & \cdots \end{array}\right] \tag{8}\\ 0 \leq p_{i j} \leq 1, & i, j \in N \\ \sum_{j} p_{i j} & =1, \quad i \in N \tag{9} \end{align*}

Sojourn Time

Given a state ii of a DTMC {XnS:nN}\left\{X_{n} \in S: n \in N\right\}, the sojourn time TiT_{i} of the state ii is the discrete r.v. that gives the no. of time steps for which the DTMC resides in state ii before transiting to a different state.

P(Ti>n)=P(Xn=iXn1=i)P(Xn1=iXn2=i)P(X2=iX1=i)P(X1=iX0=i)=piin\begin{aligned} & P\left(T_{i}>n\right) \\ & =P\left(X_{n}=i \mid X_{n-1}=i\right) P\left(X_{n-1}=i \mid X_{n-2}=i\right) \cdots \\ & \quad \cdots P\left(X_{2}=i \mid X_{1}=i\right) P\left(X_{1}=i \mid X_{0}=i\right) =p_{i i}^{n} \end{aligned} P(Ti=n)=k=nP(Ti=k)k=n+1P(Ti=k)=P(Ti>n1)P(Ti>n)=piin1piin=piin1(1pii).\begin{aligned} P\left(T_{i}=n\right) & =\sum_{k=n}^{\infty} P\left(T_{i}=k\right)-\sum_{k=n+1}^{\infty} P\left(T_{i}=k\right) \\ & =P\left(T_{i}>n-1\right)-P\left(T_{i}>n\right) \\ & =p_{i i}^{n-1}-p_{i i}^{n}=p_{i i}^{n-1}\left(1-p_{i i}\right) . \end{aligned}

Thus, TiT_{i} is a geometric r.v. with success probability (1pii)\left(1-p_{i i}\right).

E(Ti)=11piiE\left(T_{i}\right)=\frac{1}{1-p_{i i}}

6. Chapman-Kolmogorov Equation (C-K eqn)

pij(m+n)=kSpik(m)pkj(n)p_{i j}(m+n)=\sum_{k \in S} p_{i k}(m) p_{k j}(n) P(m+n)=P(m)P(n).\Longrightarrow P(m+n)=P(m) P(n) .

Let P(n)=[pij(n)]P(n)=\left[p_{i j}(n)\right] be the matrix of nn-step transition probabilities. Note that P(0)=IP(0)=I.

P(n)=P(n1)P,n1\begin{equation*} P(n)=P(n-1) \cdot P, \quad n \geq 1 \end{equation*}

State Probabilities

pj(n)=P(Xn=j),n=0,1,2,,j=0,1,2,\begin{equation*} p_{j}(n)=P\left(X_{n}=j\right), \quad n=0,1,2, \cdots, j=0,1,2, \cdots \end{equation*} pj(n)=P(Xn=j)=iSP(Xn=j,X0=i)=iSP(Xn=jX0=i)P(X0=i)=iSpi(0)pij(n),j=0,1,2,\begin{align*} p_{j}(n) & =P\left(X_{n}=j\right) \\ & =\sum_{i \in S} P\left(X_{n}=j, \quad X_{0}=i\right) \\ & =\sum_{i \in S} P\left(X_{n}=j \mid X_{0}=i\right) P\left(X_{0}=i\right) \\ & =\sum_{i \in S} p_{i}(0) p_{i j}(n), \quad j=0,1,2, \cdots \end{align*} π(n)=[p0(n)p1(n)p2(n)]\pi(n)=\left[\begin{array}{llll} p_{0}(n) & p_{1}(n) & p_{2}(n) & \cdots \end{array}\right] π(n)=π(0)P(n)=π(0)Pn,n=0,1,2,\begin{equation*} \pi(n)=\pi(0) P(n)=\pi(0) P^{n}, \quad n=0,1,2, \cdots \end{equation*}

7. Continuous Time Markov Chain (CTMC)

Sojourn Times in States

TiT_{i} is an exponential r.v., i.e., Ti=EXP(λi)T_{i}=\operatorname{EXP}\left(\lambda_{i}\right),

P(Ti>x)=eλix,x0P\left(T_{i}>x\right)=e^{-\lambda_{i} x}, \quad x \geq 0

If λi=0\lambda_{i}=0, the state ii is an absorbing state.

If λi=\lambda_{i}=\infty, the state ii is an instantaneous state. If λi(0,)\lambda_{i} \in(0, \infty), the state ii is a stable state.

8. Chapman-Kolmogorov Equation (C-K eqn)

H(s,t)=[pij(s,t)]H(s, t)=\left[p_{i j}(s, t)\right] H(s,t)=H(s,u)H(u,t),0sut\begin{equation*} H(s, t)=H(s, u) H(u, t), \quad 0 \leq s \leq u \leq t \end{equation*}

9. Kolmogorov Differential Equation

H(s,t)t=H(s,t)Q(t),0st\begin{equation*} \frac{\partial H(s, t)}{\partial t}=H(s, t) Q(t), \quad 0 \leq s \leq t \end{equation*} H(s,t)s=Q(s)H(s,t),0st\begin{equation*} \frac{\partial H(s, t)}{\partial s}=-Q(s) H(s, t), \quad 0 \leq s \leq t \end{equation*} qii(t)=limh0pii(t,t+h)1h\begin{equation*} q_{i i}(t)=\lim _{h \rightarrow 0} \frac{p_{i i}(t, t+h)-1}{h} \end{equation*} qij(t)=limh0pij(t,t+h)h,ij.\begin{equation*} q_{i j}(t)=\lim _{h \rightarrow 0} \frac{p_{i j}(t, t+h)}{h}, \quad i \neq j . \end{equation*} jqij(t)=0, for all i\begin{equation*} \sum_{j} q_{i j}(t)=0, \quad \text { for all } i \end{equation*}