Compound and Non-Homogeneous Poisson Processes

Beyond Constant-Rate, Unit-Jump Arrivals

The basic Poisson process counts events that all look the same and arrive at a constant rate. Real systems break both assumptions. In a wireless network, each arriving packet carries a random payload size. During rush hour, the call arrival rate at a base station increases. Two natural extensions handle these situations: the compound Poisson process (random marks on arrivals) and the non-homogeneous Poisson process (time-varying rate).

Definition:

Compound Poisson Process

Let {N(t)}\{N(t)\} be a Poisson process with rate λ\lambda and let Y1,Y2,Y_1, Y_2, \ldots be i.i.d. random variables (called marks or jumps) independent of {N(t)}\{N(t)\}. The process X(t)=k=1N(t)YkX(t) = \sum_{k=1}^{N(t)} Y_k is called a compound Poisson process. Its mean and variance are: E[X(t)]=λtE[Y1],Var(X(t))=λtE[Y12].\mathbb{E}[X(t)] = \lambda t \, \mathbb{E}[Y_1], \qquad \text{Var}(X(t)) = \lambda t \, \mathbb{E}[Y_1^2].

The compound Poisson process has independent and stationary increments but is not a counting process (jumps can be arbitrary). It is the continuous-time analog of a random walk with random step sizes.

,

Example: Aggregate Data Traffic

Packets arrive at a link according to a Poisson process with rate λ=1000\lambda = 1000 packets/second. Each packet has a random size YkY_k uniformly distributed on [100,1500][100, 1500] bytes. Model the total data volume X(t)X(t) as a compound Poisson process and find the mean and standard deviation of the data volume in a 1-second window.

Definition:

Non-Homogeneous Poisson Process

A counting process {N(t):t0}\{N(t) : t \geq 0\} is a non-homogeneous (inhomogeneous) Poisson process with intensity function λ(t)0\lambda(t) \geq 0 if:

  1. N(0)=0N(0) = 0.
  2. It has independent increments.
  3. For all 0s<t0 \leq s < t: N(t)N(s)Poisson ⁣(stλ(u)du).N(t) - N(s) \sim \text{Poisson}\!\left(\int_s^t \lambda(u)\,du\right).

The cumulative intensity (or mean function) is Λ(t)=0tλ(u)du\Lambda(t) = \int_0^t \lambda(u)\,du, so E[N(t)]=Λ(t)\mathbb{E}[N(t)] = \Lambda(t).

When λ(t)=λ\lambda(t) = \lambda is constant, this reduces to the homogeneous Poisson process. The non-homogeneous version models time-varying arrival rates such as diurnal traffic patterns in cellular networks.

,

Theorem: Thinning and Superposition of Poisson Processes

(Thinning.) Let {N(t)}\{N(t)\} be a Poisson process with rate λ\lambda. If each arrival is independently retained with probability pp and deleted with probability 1p1-p, the retained arrivals form a Poisson process with rate pλp\lambda, independent of the deleted process (rate (1p)λ(1-p)\lambda).

(Superposition.) If {N1(t)},,{Nm(t)}\{N_1(t)\}, \ldots, \{N_m(t)\} are independent Poisson processes with rates λ1,,λm\lambda_1, \ldots, \lambda_m, then the merged process N(t)=N1(t)++Nm(t)N(t) = N_1(t) + \cdots + N_m(t) is Poisson with rate λ1++λm\lambda_1 + \cdots + \lambda_m.

Thinning is like coin-flipping each arrival: heads keeps it, tails removes it. Since independent coin flips on independent arrivals preserve independence, the sub-processes remain Poisson. Superposition works because independent increments and orderliness are preserved under merging.

,

Why This Matters: Spatial Poisson Process for Base Station Locations

In stochastic geometry models for cellular networks, the locations of base stations in R2\mathbb{R}^2 are modeled as a homogeneous Poisson point process (PPP) with spatial intensity λBS\lambda_{\text{BS}} (base stations per km2^2). This is the two-dimensional extension of the temporal Poisson process studied in this section.

Superposition explains why the PPP is such a good fit: when many independent operators deploy base stations according to their own (possibly non-Poisson) point processes, the merged set of all base stations converges to a PPP by the superposition property.

Thinning models heterogeneous networks (HetNets): from the PPP of all base stations, independently classify each as macro, micro, or pico cell. Each tier forms an independent PPP with a fraction of the total intensity.

The key results — coverage probability, rate distribution, handover rate — all build on the Poisson properties of this chapter. See the FSP/Telecom chapters on stochastic geometry for the full development.

,

Common Mistake: Inter-Arrival Times of a Non-Homogeneous Poisson Process

Mistake:

Assuming that the inter-arrival times of a non-homogeneous Poisson process are exponential. Students often write TkExp(λ(t))T_k \sim \text{Exp}(\lambda(t)), which is not well-defined since λ(t)\lambda(t) changes over time.

Correction:

For a non-homogeneous Poisson process, the inter-arrival times are not exponential (unless λ(t)=λ\lambda(t) = \lambda is constant). The survival function of the first arrival time is P(T1>t)=eΛ(t)\mathbb{P}(T_1 > t) = e^{-\Lambda(t)}, which is exponential only when Λ(t)=λt\Lambda(t) = \lambda t. The correct characterization uses the cumulative intensity Λ(t)=0tλ(u)du\Lambda(t) = \int_0^t \lambda(u)\,du.

Common Mistake: Variance of a Compound Poisson Process

Mistake:

Computing Var(X(t))=λtVar(Y)\text{Var}(X(t)) = \lambda t \cdot \text{Var}(Y). This is wrong because it ignores the randomness in the number of terms.

Correction:

By the law of total variance: Var(X(t))=E[N(t)]Var(Y)+Var(N(t))(E[Y])2=λt(Var(Y)+(E[Y])2)=λtE[Y2]\text{Var}(X(t)) = \mathbb{E}[N(t)]\text{Var}(Y) + \text{Var}(N(t))(\mathbb{E}[Y])^2 = \lambda t(\text{Var}(Y) + (\mathbb{E}[Y])^2) = \lambda t \, \mathbb{E}[Y^2]. The correct formula uses the second moment of YY, not just the variance.

Quick Check

Claims arrive at an insurance company as a Poisson process with rate λ=2\lambda = 2/day. Each claim amount is Exp(μ)\text{Exp}(\mu) with mean μ1=500\mu^{-1} = 500 euros. What is the expected total payout in a 30-day month?

30,00030{,}000 euros

60,00060{,}000 euros

1,0001{,}000 euros

15,00015{,}000 euros

Compound Poisson process

The process X(t)=k=1N(t)YkX(t) = \sum_{k=1}^{N(t)} Y_k where N(t)N(t) is a Poisson process and YkY_k are i.i.d. marks independent of N(t)N(t). Models random-sized arrivals.

Related: Poisson process

Non-homogeneous Poisson process

A Poisson process with time-varying intensity λ(t)\lambda(t). The count in (s,t](s,t] is Poisson with mean stλ(u)du\int_s^t \lambda(u)\,du. Models time-varying arrival rates.

Related: Poisson process