Stochastic geometry and communication networks

4 downloads 10789 Views 2MB Size Report
Oct 7, 2005 ... I Panorama of some stochastic-geometry models in action, .... Stoyan & Kendal & Mekke (1995) Stochastic Geometry and its Applications. Wiley ...
Stochastic geometry and communication networks

OUTLINE I Panorama of some stochastic-geometry models in action,

Bartek Błaszczyszyn

II Basic geometric models,

tutorial lecture

III Signal-to-Interference-and-Noise ratio (SINR) coverage model,

Performance conference

IV Modeling ad-hoc networks,

Juan-les-Pins, France, October 3-7, 2005

V Power control in CDMA: from static to dynamic modeling — a spatial Erlang formula. INRIA-ENS

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

2

PANORAMA...

I PANORAMA

• The mathematical principle of the Voronoi tessellation is widely used as a simple idealization of many complex “real” partitions of the plane (cells in cellular

• Boolean model is a first model of coverage of wireless network (it does not take into account interference). As the underlying model for the study of continuum

communication, Gupta & Kumar protocol model of ad-hoc networks; it takes into

percolation it can be used to address the questions of connectivity of ad-hoc

account only locations of antennas and ignores all other physical aspects of the

networks in the absence of interference.

communication technology as a.g. additive interference)

• The dual Delaunay graph can be used as a “protocol model” of neighbuorhood in networks ⇒ topology for routing

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

3

• Mathematical representation of interferences based on (Poisson) shot noise processes ⇒ a variety of results on coverage, connectivity and capacity of large interference-limited networks.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

4

BASIC MODELS... Poisson Point Process Planar Poisson point process (p.p.)

II BASIC GEOMETRIC MODELS

• Number of Points Φ(B) of Φ in subset B of the plane is Poisson random variable with parameter λ|B|, where | · | is the Lebesgue measure on the plane; i.e.,

• Poisson point process, • Voronoi tessellation and Delaunay graph,

P{ Φ(B) = k } = e−λ|B|

• Boolean model, • Shot-Noise model,

(λ|B|)k , k!

• Numbers of points of Φ in disjoint sets are independent.

• References.

Laplace transform of the Poisson p.p.

5

BASIC MODELS/Poisson p.p. ...

R

h(x) Φ(dx)

R

h(x)

] = e−λ (1−e ),dx , R P where h(·) is a real function on the plane and h(x) Φ(dx) = Xi ∈Φ h(Xi ). LΦ (h) = E[e

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

Φ of intensity λ:

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

6

BASIC MODELS ... Voronoi Tessellation (VT) and Delaunay graph

Poisson p.p. is the basis of the stochastic-geometry modeling of communication

= {Xi } on the plane and a given point x, we define the Voronoi cell of this point Cx = Cx (Φ) as the subset of the plane of all locations that are closer to x than to any point of Φ; i.e., Given a collection of points Φ

networks. This modeling consist in treating the given architecture of the network as a snapshot of a (homogeneous) random model, and analyzing it in a statistical way. In this

Cx (Φ) = {y ∈ R2 : |y − x| ≤ |y − Xi | ∀Xi ∈ Φ} .

approach the physical meaning of the network elements is preserved and reflected in the model, but their geographical locations are no longer fixed but modeled by

= {Xi } is a Poisson p.p. we call the (random) collection of cells {CXi (Φ)} When Φ

random points of, typically, homogeneous planar Poisson point processes. Consequently, any particular detailed pattern of locations is no longer of interest.

the Poisson-Voronoi tessellation (PVT).

Instead, the method allows for catching the essential spatial characteristics of the network performance basically through the densities of these point processes (i.e.,

Edges of the Delaunay graph connect nu-

the densities of the network devices). Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

Borders of Voronoi Cells

clei of the adjacent cells. 7

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

8

BASIC MODELS/VT ...

BASIC MODELS... Boolean Model (BM)

˜ Let Φ

= {(Xi , Gi )} be a marked Poisson p.p., where {Xi } are points and {Gi }

are iid random closed stets (grains). We define the Boolean Model (BM) as the union

Ξ=

VT is a frequently used generic model of tessellation of the plane.

[ i

Points denote locations of various structural elements (devices) of the network (base station antennas and/or network controllers in cellular networks, concentrators in

Xi ⊕ G i

where x ⊕ G

Known: • Poisson distribution of the number of grains intersecting

fixed telephony, access nodes in ad hoc networks, etc.).

any given set.

Cells denote mutually disjoint regions of the plane served in some sense by these

• Asymptotic results (λ → ∞)

devices.

for the probability of complete

covering of a given set.

= {x + y : y ∈ G}.

10

10

8

8

6

6

4

4

2

2

0

2

4

6

8

10

0

2

4

6

8

10

BM with spherical grains of random radii Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

9

BASIC MODELS/BM ...

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

10

BASIC MODELS... Shot-Noise (SN) model

˜ Let Φ

= {(Xi , Si )} be a marked p.p., where {Xi } are points and {Si } are iid random variables. Given a real response function L(·) of the distance on the plane BM is a generic coverage model.

we define the Shot-Noise field

Points denote locations of various structural elements (devices) of the network.

IΦ˜ (y) =

Granis denote independent regions of the plane served these devices .

X i

In wireless networks it is a simplified model (it does not take into account

Si L(y − Xi ) .

˜ is a marked Poisson p.p. then we call I ˜ the Poisson SN. When Φ Φ

interference) for the study of coverage and connectivity.

For the Poisson SN, the Laplace transform of the vector (I Φ (y1 ), . . . , IΦ (yn )) is known for any y1 , . . . , yn

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

11

∈ R2 (via Laplace transform of the Poisson p.p.).

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

12

BASIC MODELS/SN ...

BASIC MODELS ...

References SN is a good model for interference in wireless networks.

1. Stoyan & Kendal & Mekke (1995) Stochastic Geometry and its Applications. Wiley, Chichester

Marks Si correspond to emitted powers.

2. Okabe & Boots & Sugihara (2001) Spatial Tessellations: Concepts and

Response function correspond to attenuation function.

Applications of Voronoi Diagrams. Willey, Chichester.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

13

Φ = {Xi , (Si , Ti )} marked point process (Poisson)

• Motivations,

{Xi } points of the p.p. on R2 — antenna locations,

(Si , Ti ) ∈ (R+ )2 possibly random mark of point Xi — (power,threshold) ¾ ½ Si l(y − Xi ) y: ≥ Ti cell attached to point Xi : Ci (Φ, W ) = W + κIΦ (y) P where Iφ (y) = i6=0 Si l(y − Xi ) shot noise process, κ interference factor, W ≥ 0 external noise, l(·) attenuation (response) function. Ci is the region where the SINR from Xi is bigger than the threshold Ti . Coverage PROCESS:

Ξ(Φ; W ) =

• Snapshots and qualitative results, • Typical cell study

(coverage probability for the point is some distance to the antenna, simultaneous coverage of several points, mean area of the cell),

• Handoff study

(overlapping of cells, coverage probability for a typical point, distance to different

handoff states),

• Macroeconomic optimization example,

Ci (Φ, W ).

i∈N Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

14

SINR COVERAGE MODEL ...

III SINR COVERAGE MODEL

[

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

• References. 15

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

16

SINR COVERAGE MODEL ...

SINR COVERAGE MODEL / CDMA motivation ... Parameter values Intensity of Poisson process of base stations λBS ∼ 0.2 BS/km2 .

Motivation I: CDMA handoff cells

x0 — a point in R2 (location of an antenna), s0 ≥ 0 and t0 ≥ 0 — (pilot signal power of

the antenna and SINR threshold (bit energy-tonoise spectral power density

Pilot signal power s0 (x 1 , (s1 , t1 ))

pilot signal),

y

(x 3 , (s3 , t3 ))

(x 0 , (s0 , t0 )) (x 4 , (s4 , t4 ))

s0 l(y − x0 ) ≥ t0 w + Iφ (y)

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

∼ 30 mW

SINR threshold (bit energy-to-noise spectral power density

w

Eb /NO ) for the

φ = {xi , (si , ti )} — pattern of antennas, w ≥ 0 — external noise, 0 ≥ κ ≤ 1 — orthogonality factor, l(·) — attenuation function

(x 2 , (s2 , t 2))

−14 dB

Eb /NO ) for the pilot t0

External noise w



∼ −105 dB

Interference factor for pilots from different BS’s

κ=1 Attenuation function

l(x) = A max(|x|, r0 )−α or l(x) = (1 + A|x|)−α with α ∼ 3 − 6. 17

SINR COVERAGE MODEL / CDMA motivation ...

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

18

SINR COVERAGE MODEL / Motivations II

This is a relatively simple model, which takes into account only locations of the Base Stations, their pilot signal powers and SIR’s for the pilots.

Gupta & Kummar physical model for ad-hoc networks

In particular there is no any pattern of mobiles assumed yet and it does not take into

(see Modeling of ad-hoc networks)

account power control issues.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

19

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

20

SINR COVERAGE MODEL ...

SINR COVERAGE MODEL / Snapshots ...

Snapshots and qualitative results 500

20000

20000

20000

20000

18000

18000

18000

18000

16000

16000

16000

16000

14000

14000

14000

14000

12000

12000

12000

12000

400

10000

10000

10000

10000

8000

8000

8000

8000

6000

6000

6000

6000

4000

4000

4000

4000

2000

2000

2000

2000

300

200

100

20000

20000

18000

18000

16000

16000

14000

14000

12000

12000

10000

10000

8000

8000

6000

6000

4000

4000

2000 0

2000

4000

6000

8000 10000 12000 14000 16000 18000 20000

0

2000

κ = 0.5

4000

6000

8000 10000 12000 14000 16000 18000 20000

0

2000

κ = 0.1

4000

6000

8000 10000 12000 14000 16000 18000 20000

0

κ = 0.01

2000

4000

6000

2000

8000 10000 12000 14000 16000 18000 20000

0

κ=0

≡ 1, T = 0.4 and interference factor κ → 0.

Constant emitted powers Si , ei

100

200

300

400

500

0

2000

4000

6000

8000 10000 12000 14000 16000 18000 20000

0

2000

4000

6000

8000 10000 12000 14000 16000 18000 20000

≡ 1, T = 0.4, W = 0, l(r) = (Ar)−β and attenuation exponent β → ∞.

Constant emitted powers Si , ei

Small interference factor allows one to approximate SINR cells by a Boolean model SIR cells tend to Voronoi cells whenever attenuation is stronger, e.g. in urban areas.

(via perturbation methods) Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

21

SINR COVERAGE MODEL / Typical cell study ...

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

22

SINR COVERAGE MODEL / Typical cell study ...

Probability for a typical cell to cover a point

Φ — marked Poisson point process representing antennas in R2 , (0, (S, T )) — additional antenna located at fixed point 0 with random (S, T ) distributed as any mark of Φ, independent of it (thus Φ ∪ {(0, (S, T )} has 2 Poisson Palm distribution), y — location (of a mobile) in R . Given:

IΦ (y) that is

− 1), W

and the Laplace transform of

· Z ³ ´ i 1 − LS (ξl(y − z)) µ(dz) , E[exp(−ξIΦ (y))] = exp −

where LS (ξ)

´

pR = P y ∈ C 0 ³ ´ = P S(1/T − 1)l(R) − W − IΦ (y) > 0 . Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

be given via Laplace transforms of S(1/T

Rd

Probability for C0 to cover a given point y located at the distance R to the origin:

³

Res. For M/G case (general distribution of (S, T )) the coverage probability p R can

= E[e−ξS ] is the Laplace transform of S .

Cor. Fourier transform of the Poisson shot-noise variable I φ (y) → Rieman Boundary Problem

23

→ probability of coverage by the typical cell.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

24

SINR COVERAGE MODEL / Typical cell study ...

SINR COVERAGE MODEL / Typical cell study ... Special M/M case

Example Fourier transform FIΦ (ξ) of the homogeneous Poisson (intensity λ) shot noise with

exponential S (parameter m) and attenuation l(x)

£

−iξIΦ

FIΦ (ξ) = E e " r

for ξ

= A max(|x|, r 0 )

−4

Res. [Baccelli&BB&Muhlethaler (2004)]: Assume that {Si } are exponential r.vs. with

par.

µ and Ti = T are constant. Then the probability for C0 to cover a given point

located at the distance R: is equal

¤

pR = exp

r ¶ µ r λ iAξ m iAξ − π2 arctan r02 = exp λπ m iAξ 2 m ¸ 4 4 r − iAξ − r m 0 + λπr02 0 , 4 iAξ + r0 m

− 2πλ

Z

∞ 0

¾ u du . 1 + l(R)/(T l(u))

proof: Say the emitter is at the origin and consider the corresp. Palm distribution P;

pR = P(S ≥ T IΦ1 /l(R)) Z ∞ e−µsT /l(R) dP(IΦ ≤ s) = 0

= ψIΦ (µT /l(R)) ,

∈ R, where the branch of the complex square root function is chosen with

where LIΦ (·) is the Laplace transform of the value of the hom. Poisson SN I Φ .

positive real part.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

½

25

SINR COVERAGE MODEL / Typical cell study / M/M case ...

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

26

SINR COVERAGE MODEL / Typical cell study / M/M case ...

Some optimizations One can study the following optimization problems for the expected effective transmission range r Cor. For the attenuation function l(u)

= (Au)−β

pR (λ) = e−λR where C

2 T 2/β C

• given the density of stations λ find the targeted range r that optimizes the expected effective transmission range

,

³ ´ = C(β) = 2πΓ(2/β)Γ(1 − 2/β) /β .

1 √ r≥0 T 1/β 2λpC 1 √ = rmax (p) = argmaxr≥0 {rpr (λ)} = T 1/β 2λC ρ = ρ(p) = max{rpr (p)} =

rmax Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

× pr :

27

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

28

SINR COVERAGE MODEL / Typical cell study / M/M case optimization

SINR COVERAGE MODEL / Typical cell study ...

Probability for a typical cell to cover two points

(y1 , y2 ) — two point to be covered by a given cell C0 (Φ, W ) under Palm distribution of Φ ∪ {(0, (S, T )}

• given the targeted range R find the density of emitters λ that optimize the expected effective transmission range R × pR :

We need the joint Laplace transform of

(IΦ (y1 ), IΦ (y2 )) that is given by h ³ ´¤ E exp −ξ1 IΦ (y1 ) − ξ2 IΦ (y2 ) · Z ³ ´ i 1 − LS (ξ1 l(y1 − z) + ξ2 l(y2 − z)) µ(dz) . = exp −

1 R2 T 2/β C 1 max{RpR (λ)} = 2 2/β λ≥0 R T eC

λmax = λmax (R) =argmaxλ≥0 {RpR (λ)} =

Rd

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

29

SINR COVERAGE MODEL / Typical cell study ...

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

30

SINR COVERAGE MODEL / Typical cell study / Perturbation of Boolean model ...

Coverage probability via perturbation of Boolean model valid for small interference factor κ

Res.

(κ) (κ) Denote pR = P(x ∈ C0 ), where |x| = R and n o (κ) 2 C0 = y ∈ R : Sl(y) ≥ κIΦ (y) + W .

Assume F∗ (u)

= P((Sl(x) − W ) ≤ u) admits Taylor approximation at 0:

F∗ (u) = F∗ (0) + and R∗ (u)

(κ)

pR =

h (k) X F∗ (0) k=1

k!

k

u + R∗ (u)



}| ´{ P Sl(x) ≥ W provided E[(IΦ (x))

h

= o(u ) u & 0.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

correcting terms value for the Boolean model

31

2h



z h X k=1

κ

}|

(k) k F∗ (0)

k!

] < ∞.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

E

£

(IΦ (y))

{ ¤ k

error

z }| { + o(κh ) ,

32

SINR COVERAGE MODEL / Typical cell study ...

SINR COVERAGE MODEL / Typical cell study ... Numerical examples

p(κ) (1)

p(0.1) (y)

1

Mean cell area formula

v (κ) (1)

1

2 1.8

0.8

= E[|C0 |]. Recall that pR is the coverage probability for location at distance R. Denote the mean area of the cell of the BS located at 0 by v 0

0.8

1.6 1.4

0.6

0.6

1.2

0.4

0.4

0.8

0.2

0.2

0.4

1

0.6

Then

v0 =

Z

0.2

Rd

p|y| dy.

0

0.2

0.4

0.6

1

0.8

κ

0 1

1.2

1.4

1.6

1.8

2

function of the distance.

coefficient Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

33

SINR COVERAGE MODEL / Handoff study ...

0.02

0.04

0.06

0.08

0.1

κ

y

Probability of coverage as a Probability of coverage as a function of the interference

0

Mean cell area as a function of the interference coefficient.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

34

SINR COVERAGE MODEL / Handoff study / Overlapping of cells...

Overlapping of cells Deterministic scenario: given n cells C(xi , si , ti ; φ, w), i = 1, . . . , n Pn Res. The inequality i=1 ti /(1 + ti ) < 1 is a necessary condition for Tn i=1 C(xi , si , ti ; φ, w) 6= ∅ Random scenario:

Cor. If the distribution of the ratio T is such that T

≥τ

for some τ

> 0, then the number Ky of cells of the coverage process Ξ covering any given point y is a.s. bounded

Example: For the maximal pilot’s bit energy-to-noise spectral power density

τ = Eb /NO = −14 dB the pole handoff number (theoretical maximal handoff number) K ≤ 26.

1+τ . τ (Given point cannot be covered by (1 + τ )/τ or more cells, no matter how close Ky
1/(β/2−δ)

< C/|x − y| β for some

0, β > 0 and if the distribution of S has finite moment E[S ] < ∞ for some δ ∈ [1, β/2], then one can show that for any R, ε, α > 0³, there exists R0 > 0 such that ´ P P sup|y|R0 Si l(y, Xi ) < ε > 1 − α . Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

˜t ) is obtained using The exact stationary distribution of the Markov process ( Z backward simulation (coupling from the past) similar to that proposed by Kendall. Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

41

SINR COVERAGE MODEL ...

42

SINR COVERAGE MODEL ...

Macroeconomic optimization example:

densification / magnification

Increase the mean power m of existing antennas or increase the density λ of antennas?

References 2

C total budget of an operator per km , Cλ cost of one antenna, Cm cost of increasing the power of one antenna by 1W. constraint: λCλ + Cm λ/m = C.

1. Baccelli & Błaszczyszyn (2001), On a coverage process ranging from the Boolean model to the Poisson Voronoi tessellation, with applications to wireless communications Adv. Appl. Probab. 33 mean handoff level EN0

2. Tournois (2002), Perfect Simulation of a Stochastic Model for CDMA coverage

10

INRIA report 4348,

8

Plots of mean handoff as a functions of mean antenna power

1/m 6

under budget constraint with top: Cm

3. Baccelli, Błaszczyszyn & Tournois (2002), Spatial averages of downlink coverage

C = 1000, Cλ = 500 and from the

= 1, 2, 5.

Solution: Plot maximum = Optimal configuration.

4

0

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

characteristics in CDMA networks (INFOCOM ).

2

2

4

6

8

10

12

mean antenna power 1/m

43

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

44

AD-HOC NETWORKS ...

• A random set of users distributed in

IV MODELING AD-HOC NETWORKS

space and sharing a common Hertzian

• Ad-hoc networks,

medium.

• A few sg “interference aware” models for ad-hoc networks,

• Users constitute ad-hoc network that is

• Some optimization problem in capacity and medium access control,

emitter receiver

in charge of transmitting information far away via several hops.

• References.

• Users switch between emitter and receiver modes.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

45

AD-HOC NETWORKS ...

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

46

AD-HOC NETWORKS ...

Interference and successful reception

Multi-hop transmissions

• Emitter sends a packet in some given di-

• Emitter sends a packet emitting some

• The packet is received by some number

• Transmission distance attenuates the

power.

rection far away via several hops.

emitted power.

(possibly 0) of neigbouring receivers.

• An optimal receiver among them is in

emitter

charge of forwarding this packet in (one receiver

towards destination optimal receiver

emitter used sign

receivers.

of) his next emission time-slots.

al receiver

interferen

ce

• Reception is successful if the Signal to In-

• In the case of no reception, emitter re-

terference Ratio (SIR) at the receiver is large enough.

emits the packet next authorized time.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

• Emitted power causes interference at all

47

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

48

AD-HOC NETWORKS ...

AD-HOC NETWORKS / A few models / General settings ...

A few sg interference aware models

General settings

A way of choosing marks ei for nodes is called medium access control. Two

• Nodes distributed in (a subset of) the plane according to a Poisson p.p. Φ = {Xi },

possibilities can be considered:

• Nodes Xi are marked in some way (not necessarily independently) by ei = 1 when node is emitting or 0 when not (it is a potential receiver); • Call – –

• Some “central authority” assigns marks ei (in some optimal way) for each given configuration of nodes

⇒ {ei } are dependent in some way. • Each node independently switches between modes ei = 1 and ei = 0 ⇒ {ei } are independent, given configuration of nodes.

Φ1 = {Xi ∈ Φ : ei = 1} — emitters and

Φ0 = {Xi ∈ Φ : ei = 0} — potential receivers.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

49

AD-HOC NETWORKS / A few models ...

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

50

AD-HOC NETWORKS / A few models / Protocol model ... Examples of medium access centralized

independent

Voronoi tessellation principle: Protocol Model [Gupta & Kumar (2000)]

• When the node Xi ∈ Φ1 transmits it can be successfully received by node Xj ∈ Φ0 if |Xj − Xi | ≤ (1 + ∆)|Xj − Xk | ∀Xk ∈ Φ1 , where ∆

≥ 0 is some constant. ⇔ Xj ∈ modified Voronoi Cell CXi (Φ1 ) The Protocol Model with centralized medium access was studied asymptotically

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

(when number of nodes goes to ∞) by G & K (2000); to be commented. 51

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

52

AD-HOC NETWORKS / A few models ...

AD-HOC NETWORKS / A few models / Carrier-Sense model ... Examples of medium access

Hard-core principle: Carrier-Sense Model

´ hard-core model: • Matern ⇒ Nodes are first independently marked by some auxiliary marks in [0, 1]. A node Xi is selected as emitter (ei = 1) if its auxiliary mark is larger than all auxiliary marks within its neighbourhood of range R cs . Otherwise the node is

Close to some currently used protocols implemented e.g. in IEEE 802.11

• When the node Xi ∈ Φ1 transmits it can be successfully received by node Xj ∈ Φ0 if

marked receiver.

|Xj − Xi | ≤ R and |Xk − Xi | > Rcs ∀Xk ∈ Φ1 , where Rcs

> R ≥ 0 are some constants.

⇔ Each communication within the transmission range R is protected by the exclusion disc centered at the transmitter with radius R cs > R. (One can also

• Gibbs hard-core model: ⇒ The steady state distribution of a spatial birth-and-death process of nodes. When a node is born it is marked as emitter if in its neighbourhood of range R cs there is no any other emitter. Otherwise it is marked receiver.

think of the models with exclusion discs centered at the receiver and at both

The Carrier Sense Model, has not been studied yet by means of stochastic geometry

transmitter and receiver.)

tools (to the best of our knowledge).

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

53

AD-HOC NETWORKS / A few models ...

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

54

AD-HOC NETWORKS ...

SINR principle: Physical model

• Points Xi are independently identically marked by some Si ≥ 0 representing powers used when emitting signals.

• l(r) = (Ar)−β signal attenuation function, • T signal-to-interference ratio threshold, Connectivity for the physical model

• W > 0 constant external noise, κ > 0 interference factor.

• When the node Xi ∈ Φ1 transmits it can be successfully received by node Si l(|Xj − Xi |) Xj ∈ Φ0 if ≥T, W + κIΦ1 \Xi (Xj ) P where IΦ1 is the shot-noise process of Φ1 IΦ1 (y) = Xk ∈Φ1 Sk l(|y − Xk |). ⇔ Communication is successful if the signal-to-interference ratio at the receiver is bigger than the threshold T . Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

55

See Performance’05 tutorial by P. Thiran.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

56

AD-HOC NETWORKS ...

AD-HOC NETWORKS / Some optimization problem ...

Some optimization problem in capacity and medium access for the physical model Suppose emitter X0

= 0 has to send information to infinity along the real axis. Define also a modified progress

Define effective distance traversed in one hop, called also progress effective distance

effective distance

probab. of z ³ }| µsuccessful ´+{¶ z }| rec.{ ˜ = max p|X | (λp) |Xj | cos(arg(Xj )) . D j 0

rec. ind.z }| µsuccessful ³ ´+{¶ z }| { , D = max0 δ(Xj , 0, Φ1 )|Xj | cos(arg(Xj ))

Xj ∈Φ

Xj ∈Φ

Denote d(λ, p)

effective distance emitter optimal receiver

˜ p) = E[D] ˜. = E[D], d(λ,

Search for receivers that realize the progresses, D and

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

57

AD-HOC NETWORKS / Some optimization problem ...

˜ can be implemented. D

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

58

AD-HOC NETWORKS / Some optimization problem ...

Fact The mean total distance traversed in one hop by all transmissions initialized in some unit area (called density of (modified) progress) is equal to λp

d(λ, p) (resp.

˜ p)). λp d(λ,

³ zet ´ G(z) = 2 arccos √ dt . √ 2et {t:et / 2et≤1/z}

⊂ R of unit area, by the Campbell’s formula · X ¸ Z E Di = λp 1(x ∈ B)E[D0 ] dx Xi ∈Φ1 ∩B

Fact For all λ, p

∈ [0, 1], define an auxiliary function Z

2

proof: For B

For z

R2

Res. For the M/M model with W

= 0 The distribution function of the modified

progress is given by

= λp d(λ, p) .

˜ ≤ z) = e−λ(1−p)(rmax (λp))2 G(z/ρ(λp)) . FD˜ (z) = P(D

˜ p). > 0 we have d(λ, p) ≥ d(λ,

proof follows from Jensen’s inequality.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

59

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

60

AD-HOC NETWORKS / Some optimization problem ...

˜ is an extremal shot-noise maxX ∈Φ0 proof: D i

AD-HOC NETWORKS / Some optimization problem ...

g(Xi ) with the response function

g(x) = |x|p|x| (cos(arg(x)))+ . Its distribution function can be expressed by

the Laplace transform of the (additive) shot noise

P( max

Xi ∈Φ0

·

g(Xi ) ≤ z) = E exp

·X Xi

∈Φ0

ln(1(g(Xi ) ≤ z))

Cor. For the M/M model with W

¸¸

˜ is equal to = 0 the expectation of D

˜ p) = E[D] ˜ d(λ, Z 1 h³ 1 1 ´ G(z) i = 1/β √ 1 − exp 1 − dz , p 2T 2/β C T 2λpeC 0

Φ0 with intensity λ(1 − p) · ¸ Z ˜ ≤ z) = exp −λ(1 − p) P( D 1(g(x) > z) dx .

and thus, for Poisson p.p.

R2

R

Passing to polar coordinates in the integral R2

. . . dx completes the proof.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

61

AD-HOC NETWORKS / Some optimization problem ...

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

62

AD-HOC NETWORKS / Some optimization problem ...

˜ p) is attained for p Res. The maximal density of modified progress λp d(λ, satisfyingZ 1 ³

0

The optimization

h³ 1 ´ G(z) i G(z) ´ exp 1 − 1+ dz = 1 . p 2T 2/β C pT 2/β C ˜ p) λpd(λ,

Goal: optimize the density of progress λp d(λ, p) in MAP p, for fixed λ.

0.006

Note that replacing d by d˜ in the optimization problem will result in conservative

0.005

0.004

bound on the density of progress. 0.003

Example: Modified transport capacities for 0.002

the special case with T

= {13, 15, 17}dB

0.001

(curves from top to bottom).

0

0.02

0.04

0.06

0.08

0.1

0.12

0.14

p Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

63

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

64

AD-HOC NETWORKS / Some optimization problem ...

AD-HOC NETWORKS ...

Relation to Gupta & Kumar’s (2000) result

• The transport capacity of a bounded network with Protocol √and SIR (Physical) Model with centralized medium access is proportional to

λ (as λ → ∞).

This law holds true even if the network architecture and operation is optimally organized in a centralized manner.

• Gupta & Kumar (2000), The Capacity of Wireless Networks, IEEE Inf. Theory, • Dousse & Baccelli & Thiran (2003) Impact of Interferences on the Connectivity of

• Our density of progress can be seen as an instantaneous transport capacity.

Ad Hoc Networks, IEEE INFOCOM,

• We show that the independent medium access gives an instantaneous transport √ capacity of K(p)

References

λ.

• Closed form of K(p) allows for optimization of K(p) in p.

• Baccelli & Blaszczyszyn & Muhlethaler (2004), An Aloha protocol for multihop

mobile ad-hoc wireless networks, Proc of 16th ITC Specialist Seminar to appear

in IEEE Inf. Theory

• Question: Is this performance implementable and stable in the strong sense (existence of positive recurrent Markov process describing the system)?

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

65

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

66

POWER CONTROL ... CDMA: Interference limited radio channel

V POWER CONTROL IN CDMA: FROM STATIC TO DYNAMIC MODELING

• CDMA Power allocation algebra,

(Y2 , S 2 )

(Y5 , S5 )

• Decentralized Power Allocation Principle DPAP,

W

• Network architecture models, • Maximal load estimates,

X

(Y1 ,S1 )

• Feasibility probabilites, (Y3 ,S3 )

• Blocking rates via a spatial Erlang formula, • References. Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

(Y4 ,S4 )

S1 l(X − Y1 ) ≥C W + I(X) 67

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

68

POWER CONTROL ...

POWER CONTROL / Power allocation algebra...

CDMA Power allocation algebra NBS = {Yj }j : locations of base-stations (BS’s) in R2 , NM = {Xij }i : locations of mobiles served by BS No. j , Cij : SINR required for mobile Xij , Wij : total non-traffic noise at Xij (from common overhead channels, thermal noise), κj , γ : orthogonality factors,

Local and global problem

l(x, y): path-loss from y to x.

Power allocation feasible

Power allocation feasible if exist antenna powers 0



Sij

< ∞ such that

Sij l(Yj − Xij ) j X X X ≥ C i Wij + κj l(Yj , Xij ) Sij0 + γ l(Yk , Xij ) Sik0 |

{z

i0 6=i

own-cell interference

|

all

i, j .

for each BS

{z

}

other-cell interference

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

69

POWER CONTROL / Power allocation algebra...

Global problem

Fix total powers emitted on traffic channels by other BS’s:

Sk =

Power allocation is locally feasible in cell j if exist powers 0

P

Suppose for each BS local power allocation is feasible.

k i0 Si0 (k

≤ Sij < ∞ s.t.

6= j ).

Sij l(Yj − Xij ) X X ≥ Cij j j P j j Wi + κj l(Yj , Xi ) i0 6=i Si0 + γ l(Yk , Xi ) Sik0

κj

X i

P

{z

j k6=j l(Yk ,Xi )Sk

Cij

1+

κj Cij

bj =

i

i

l(Yj , Xij )

X HjW j i i i

}

X H j l(Yk , X j ) i

i

all .

l(Yj , Xij )

,

for j

6= k

and

j where Hi

ajj =

X

κj Hij ,

i

=

Cij

1 + κj Cij

.

Denote the matrix (ajk )

= A, the vector (bj ) = b and (Sj ) = S. Global power allocation is feasible if exist antenna powers 0 ≤ S j < ∞ (total

< 1.

powers emitted on traffic channels) such that

S ≥ b + AS

(local) pole capacity condition

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

Define:

ajk = γ

i0

k6=j

fixed=γ

70

POWER CONTROL / Power allocation algebra...

Fix BS j .

Res. Local power allocation is feasible iff

Global power allocation feasible

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

Local problem

|

and

i0

k6=j

}

Local power allocation feasible

m

71

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

72

POWER CONTROL / Power allocation algebra...

POWER CONTROL ...

Res.

Decentralized Power Allocation Principle (DPAP)

• Global power allocation feasible iff the spectral radius of the matrix A is less

j Each BS j verifies for the pattern NM of the mobiles it controls if

than 1.

• The minimal solution S is equal to

X

n

κj

A b.

n

|

• The minimal solution can be obtained as the limit of the iteration

equivalently

• A sufficient condition for the spectral radius to be less than one is that A is Principle

⇒ Decentralized Power Allocation

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

73

POWER CONTROL ...

}

|

{z

< 1,

}

other-cell-interference correcting term

total path loss of user i Hij < 1. own-BS path loss of user i j | {z } Xij ∈NM

X

user’s i weight

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

74

{Yj } located according to Poisson p.p, with intensity λBS .

• BS’s

located according to hexagonal

• All antenna parameters are i.i.d. marks.

grid, p.p, with spatial density λBS .

• All mobiles form independent Poisson p.p.

• All antenna parameters are i.i.d. marks.

NM

• All mobiles form independent Poisson p.p.

NM

{z

k6=j X j ∈N j i M

j j l(Yk , Xi ) Hi l(Yj , Xij )

Poisson-Voronoi (P-V) model (“too random”)

Hexagonal (Hex) model (“too regular”)

{Yj }

j Xij ∈NM



X X

POWER CONTROL / Network architectures ...

Network architecture models

• BS’s

Hij

pole capacity term

Ab, A2 b, . . .

substochastic (has row-sums less then 1).

X

• Each mobile is served by the nearest BS.

with spatial density λM .

j (Equivalently: Each BS j serves mobiles NM

• Each mobile is served by the nearest BS.

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

with intensity λM .

in its Voronoi cell.)

75

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

76

POWER CONTROL ...

POWER CONTROL / Maximal load estimations ... Explicit mean formulas

(Wrong) Idea: Given density of BS’s λBS find maximal density of mobiles λM ,

¯ – mean number of users per cell, R ∼ – mean distance between BS’s; M λBS = 1/(πR2 ), L(r) = (Kr)α – path-loss, α – path-loss exponent, ¯ ∼ – bit rate (SINR threshold) κ – (downlink) orthogonality factor, H Pmax – power limit

such that power allocation is feasible with probability 1.

downlink:

Maximal load estimations of P-V and Hex model

Res. Given density of BS’s λBS , for any λM > 0 in both P-V and Hex model, the spectral radius of A is equal ∞ with probability 1, and thus power allocation in not feasible!

1 − Pcch /Pmax ¯ ≤ M ¯ + f¯ + L(R)g(M ¯ , α)/Pmax ) H(κ f¯H = 1/(α − 2), f¯P V = 2/(α − 2),

g¯H () ≈ . . . , g¯P V () = . . .

uplink:

Conclusion: A reduction of mobiles (admission control) is necessary for any λM > 0. Calculate blocking probabilities.

1 ¯ ≤ M ¯ + f¯ + L(R)h(M ¯ , α)/Pmax ) H(1 ¯ P V () = . . . ¯ H () ≈ . . . , h h

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

77

POWER CONTROL / Maximal load estimations ...

Stochastic geometry and communication networks, B. Błaszczyszyn; tutorial lecture, Performance’05, Juan-les-Pins, France, October 3-7, 2005

78

POWER CONTROL ...

¯ M 35

30

Feasibility probabilities for DPAP

25 20 15

P

10 5 0

1

2

3

4

5

µ

j Xij ∈NM

|

R

Hexagonal model Poisson-Voronoi model

X

Hij

total path loss of user i own-BS path loss of user i