Random multiplicative processes

Introduction

Examples of random multiplicative processes include the distributions of incomes, rainfall, and fragment sizes in rock crushing processes. Consider the latter for which we begin with a rock of size w. We strike the rock with a hammer and generate two fragments whose sizes are pw and qw, where q = 1 - p. In the next step the possible sizes of the fragments are p²w, pqw, qpw, and q²w. What is the distribution of the fragments after N blows of the hammer?

    To answer this question, consider a binary sequence in which the numbers x₁ and x₂ appear independently with probabilities p and q respectively. If there are N elements in the product Π, we can ask what is <Π>, the mean value of Π? To compute <Π>, we define P(n) as the probability that the product of N independent factors of x₁ and x₂ has the value x₁ⁿ x₂^N-n. This probability is given by the number of sequences where x₁ appears n times multiplied by the probability of choosing a specific sequence with x₁ appearing n times:

P(n) = N!/(n! (N - n)!) pⁿ q^N-n.

The mean value of the product is given by

<Π>_mp = Σ_n=0 P(n) x₁ⁿx₂^N-n = (px₁ + qx₂)ⁿ.

    The most probable event is one in which the product contains Np factors of x₁ and Nq factors of x₂. Hence, the most probable value of the product is

Π = (x₁^px₂^q)^N.

Problems

1. The average value of the sum of random variables is a good approximation to the most probable value of the sum. Is there is a similar relation for a random multiplicative process? First consider x₁ = 2, x₂ = 1/2, and p = q = 1/2. Determine <Π> and <Π>_mp s.
2. Use the program to estimate <Π> and <Π>_mp for the same parameters as used to calculate the results analytically. Do your estimated values converge more or less uniformly to the exact values as the number of measurements becomes large? Do a similar simulation for N = 20. Compare your results with a similar simulation of a random walk and discuss the importance of extreme events for random multiplicative processes.

References

• S. Redner, Random multiplicative processes: An elementary tutorial, Am. J. Phys. 58, 267–273 (1990).

Java Classes

• ProductProcessApp

Updated 21 February 2007.