Solution

ID: past-exam-of-the-mathematics-course-of-the-university-of-cambridge/2015/iii/paper-39/2/solution

For the Erlang loss formula, assume a Poisson process of arrivals with rate , independent holding times with exponential distribution of mean , unit capacity usage per accepted call, and rejection without waiting when all units are occupied. Here is a positive integer and the offered load is . The call population on is a birth-death process; detailed balance gives
By Poisson arrivals see time averages, the fraction of arrivals rejected is the time probability of full occupancy. Therefore
The paper's is the offered load, rather than necessarily an arrival rate. Exponential holding times suffice for this derivation; no unproved insensitivity of loss networks assumption is needed.
A loss network with fixed routing has capacities , call types , and a link-route incidence matrix with exactly when a type- call requires link . An accepted call holds one unit on every link of its route for its holding time; if any required link is full, the call is rejected everywhere. Take independent Poisson processes of calls, independent holding times, and route offered loads . For the finite-variable formulation below, take ; a zero-capacity link and all routes using it may be removed first.
The Erlang fixed point approximation treats blocking at different links as independent and thins the traffic offered to link by acceptance at the other links. Its reduced offered load is
The equality uses : only routes with contribute. Treating this reduced stream as Poisson leads to . Independence and Poisson thinning here are the approximation, not exact properties of a general loss network.
To construct the convex potential for the Erlang fixed point, write and , the carried load of an Erlang loss resource, equivalently the mean occupancy in the single-link stationary distribution. Differentiating its finite sums gives, for ,
The variance is positive because the finite occupancy law is nondegenerate for . Furthermore maps increasingly onto , while increases from zero to . Set
In particular , is continuous and strictly increasing, and as .
For
the first sum is convex, and each integrated strictly increasing function is strictly convex. Thus is a strictly convex function on . Since tends to a positive limit on every coordinate, is a coercive function there; it has a unique minimizer. Its coordinate derivative is
The Erlang fixed point equation is exactly . For an unused link, and the optimum is ; for a used link with positive offered load, the derivative at zero is negative and the optimum is positive. Hence every fixed point gives the unique minimizer, and the minimizer gives a fixed point. The Erlang fixed point is unique under these fixed-routing assumptions.
A convergent cyclic substitution for the Erlang fixed point updates one coordinate at a time, repeatedly cycling through all links:
Always use the most recently available values of the other coordinates. Start, for example, with . In coordinates this is exact coordinate descent for : with other coordinates fixed, the displayed update is the unique coordinate minimizer.
Here is a convergence proof. All iterates stay in the initial compact sublevel set of the coercive function , and their potential values decrease to a limit. Each single-coordinate update is continuous, since its reduced load is continuous and finite and . Thus a whole-sweep map is continuous. If a full-sweep subsequence converges to , continuity gives , because both consecutive potential values converge to the same limit. Strict coordinate convexity means equality can occur only when no coordinate changes during that sweep. Thus minimizes every coordinate, satisfies the nonnegative-orthant first-order conditions, and is the unique global minimizer of . Every subsequential limit is therefore the same point, proving convergence of the full-sweep iterates and, by continuity of the updates, all intermediate iterates. This argument applies to sequential substitution; it does not assume simultaneous updates converge.

New to topics? Read the docs here!