This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Fractional Power Control for
Decentralized Wireless Networks

Nihar Jindal, Steven Weber, Jeffrey G. Andrews The contact author N. Jindal (nihar@umn.edu) is with the University of Minnesota, S. Weber is with Drexel University, J. Andrews is with the University of Texas at Austin. This research was supported by NSF grant no. 0634763 (Jindal), no. 0635003 (Weber), nos. 0634979 and 0643508 (Andrews), and the DARPA IT-MANET program, grant no. W911NF-07-1-0028 (all authors). An early, shorter version of this work appeared at Allerton 2007 [1]. Manuscript date: September 22, 2025.
Abstract

We consider a new approach to power control in decentralized wireless networks, termed fractional power control (FPC). Transmission power is chosen as the current channel quality raised to an exponent s-s, where ss is a constant between 0 and 1. The choices s=1s=1 and s=0s=0 correspond to the familiar cases of channel inversion and constant power transmission, respectively. Choosing s(0,1)s\in(0,1) allows all intermediate policies between these two extremes to be evaluated, and we see that usually neither extreme is ideal. We derive closed-form approximations for the outage probability relative to a target SINR in a decentralized (ad hoc or unlicensed) network as well as for the resulting transmission capacity, which is the number of users/m2 that can achieve this SINR on average. Using these approximations, which are quite accurate over typical system parameter values, we prove that using an exponent of s=12s^{*}=\frac{1}{2} minimizes the outage probability, meaning that the inverse square root of the channel strength is a sensible transmit power scaling for networks with a relatively low density of interferers. We also show numerically that this choice of ss is robust to a wide range of variations in the network parameters. Intuitively, s=12s^{*}=\frac{1}{2} balances between helping disadvantaged users while making sure they do not flood the network with interference.

I Introduction

Power control is a fundamental adaptation mechanism in wireless networks, and is used to at least some extent in virtually all terrestrial wireless systems. For a single user fading channel in which the objective is to maximize expected rate, it is optimal to increase transmission power (and rate) as a function of the instantaneous channel quality according to the well-known waterfilling policy [2]. On the other hand, if the objective is to consistently achieve a target rate (or SNR), then the power should be adjusted so that this target level is exactly met. Such an objective is philosophically the opposite of waterfilling, since power is inversely related to the instantaneous channel quality: we call this channel inversion. Although suboptimal from an information theory point of view, some channel inversion is used in many modern wireless systems to adapt to the extreme dynamic range (often >50>50 dB due to path loss differences as well as multipath fading) that those systems experience, to provide a baseline user experience over a long-term time-scale.

I-A Background and Motivation for Fractional Power Control

In a multi-user network in which users mutually interfere, power control can be used to adjust transmit power levels so that all users simultaneously can achieve their target SINR levels. The Foschini-Miljanic algorithm is an iterative, distributed power control method that performs this task assuming that each receiver tracks its instantaneous SINR and feeds back power adjustments to its transmitter [3]. Considerable work has deeply explored the properties of these algorithms, including developing a framework that describes all power control problems of this type [4], as well as studying the feasibility and implementation of such algorithms [5, 6], including with varying channels [7]; see the recent monographs [8][9] for excellent surveys of the vast body of literature. This body of work, while in many respects quite general, has been primarily focused on the cellular wireless communications architecture, particularly in which all users have a common receiver (i.e., the uplink). More recently, there has been considerable interest in power control for decentralized wireless networks, such as unlicensed spectrum access and ad hoc networks [10, 11, 12, 13, 14, 15]. A key distinguishing trait of a decentralized network is that users transmit to distinct receivers in the same geographic area, which causes the power control properties to change considerably.

In this paper, we explore the optimal power control policy for a multi-user decentralized wireless network with mutually interfering users and a common target SINR. We do not consider iterative algorithms and their convergence. Rather, motivated by the poor performance of channel inversion in decentralized networks [16], we develop a new transmit power policy called fractional power control, which is neither channel inversion nor fixed transmit power, but rather a trade-off between them. Motivated by a recent Motorola proposal [17] for fairness in cellular networks, we consider a policy where if HH is the channel power between the transmitter and receiver, a transmission power of HsH^{-s} is used, where ss is chosen in [0,1][0,1]. Clearly, s=0s=0 implies constant transmit power, whereas s=1s=1 is channel inversion. The natural question then is: what is an appropriate choice of ss? We presume that ss is decided offline and that all users in the network utilize the same ss.

I-B Technical Approach

We consider a spatially distributed (decentralized) network, representing either a wireless ad hoc network or unlicensed spectrum usage by many nodes (e.g., Wi-Fi or spectrum sharing systems). We consider a network that has the following key characteristics.

  • Each transmitter communicates with a single receiver that is a distance dd meters away.

  • Channel attenuation is determined by path loss (with exponent α\alpha) and a (flat) fading value HH.

  • Each transmitter knows the channel power to its intended receiver, but has no knowledge about other transmissions.

  • All multi-user interference is treated as noise.

  • Transmitters do not schedule their transmissions based on their channel conditions or the activities of other nodes.

  • Transmitter node locations are modeled by a homogeneous spatial (2-D) Poisson process.

These modeling assumptions are made to simplify the analysis, but in general reasonably model a decentralized wireless network with random transmitter locations, and limited feedback mechanisms. In particular, the above assumptions refer to the situation where a connection has been established between a transmitter and receiver, in which case the channel power can be learned quickly either through reciprocity or a few bits of feedback. It is not however as easy to learn the interference level since it may change suddenly as interferers turn on and off or physically move (and reciprocity does not help). The fixed transmit distance assumption is admittedly somewhat artificial, but is significantly easier to handle analytically, and has been shown to preserve the integrity of conclusions even with random transmit distances. For example, [16, 18] prove that picking the source-destination distance dd from an arbitrary random distribution reduces the transmission capacity by a constant factor of E[d2]/(E[d])21E[d^{2}]/(E[d])^{2}\geq 1. Therefore, although fixed distance dd can be considered best-case as far as the numerical value of transmission capacity, this constant factor will not change fractional power control’s relative effect on the transmission capacity, which is the subject of this paper.

I-C Contributions and Organization

The contributions of the paper are the suggestion of fractional power control for wireless networks and the derivation of the optimum power control exponent s=12s^{*}=\frac{1}{2}. The exponent s=12s=\frac{1}{2} is shown to be optimal for an approximation to the outage probability/transmission that is valid for relatively low density networks that are primarily interference-limited (i.e., the effect of thermal noise is not overly large); if the relative density or the effect of noise is large, then our numerical results show that no power control (s=0s=0) is generally preferred. In the relatively large parameter space where our primary approximation is valid, fractional power control with the choice s=12s^{*}=\frac{1}{2} is shown to greatly increase the transmission capacity of a 1-hop ad hoc network for small path loss exponents (as α2\alpha\to 2), with more modest gains for higher attenuation channels. The results open a number of possible avenues for future work in the area of power control, and considering the prevalence of power control in practice, carry several design implications.

The remainder of the paper is organized as follows. Section II provides background material on the system model, and key prior results on transmission capacity that are utilized in this paper. Section III holds the main results, namely Theorem 3 which gives the outage probability and transmission capacity achieved by fractional power control, and Theorem 4 which determines the optimum power control exponent ss^{*} for the outage probability approximation. Section IV provides numerical plots that explore the numerically computed optimal ss^{*}, which provides insight on how to choose ss in a real wireless network. Section V suggests possible extensions and applications of fractional power control, while Section VI concludes the paper.

II Preliminaries

II-A System Model

We consider a set of transmitting nodes at an arbitrary snapshot in time with locations specified by a homogeneous Poisson point process (PPP), Π(λ)\Pi(\lambda), of intensity λ\lambda on the infinite two-dimensional plane, 2\mathbb{R}^{2}. We consider a reference transmitter-receiver pair, where the reference receiver, assigned index 0, is located without loss of generality, at the origin. Let XiX_{i} denote the distance of the ii-th transmitting node to the reference receiver. Each transmitter has an associated receiver that is assumed to be located a fixed distance dd meters away. Let Hi0H_{i0} denote the (random) distance–independent fading coefficient for the channel separating transmitter ii and the reference receiver at the origin; let HiiH_{ii} denote the (random) distance–independent fading coefficient for the channel separating transmitter ii from its intended receiver. We assume that all the HijH_{ij} are i.i.d. (including i=ji=j), which implies that no source-destination (S-D) pair has both a transmitter and receiver that are very close (less than a wavelength) to one another, which is reasonable. Received power is modelled by the product of transmission power, pathloss (with exponent α>2\alpha>2), and a fading coefficient. Therefore, the (random) SINR at the reference receiver is:

SINR0=P0H00dαiΠ(λ)PiHi0Xiα+η,\displaystyle{\rm SINR}_{0}=\frac{P_{0}H_{00}d^{-\alpha}}{\sum_{i\in\Pi(\lambda)}P_{i}H_{i0}X_{i}^{-\alpha}+\eta}, (1)

where η\eta is the noise power. Recall our assumption that transmitters have knowledge of the channel condition, HiiH_{ii}, connecting it with its intended receiver. By exploiting this knowledge, the transmission power, PiP_{i}, may depend upon the channel, HiiH_{ii}. If Gaussian signaling is used, the corresponding achievable rate (per unit bandwidth) is log2(1+SINR0)\log_{2}(1+{\rm SINR}_{0}). The Poisson model requires that nodes decide to transmit independently, which corresponds in the above model to slotted ALOHA [19]. A good scheduling algorithm by definition introduces correlation into the set of transmitting nodes, which is therefore not well modeled by a homogeneous PPP. We discuss the implications of scheduling later in the paper.

II-B Transmission Capacity

In the outage-based transmission capacity framework, an outage occurs whenever the SINR falls below a prescribed threshold β\beta, or equivalently whenever the instantaneous mutual information falls below log2(1+β)\log_{2}(1+\beta). Therefore, the system-wide outage probability is

q(λ)=(SINR0<β)q(\lambda)=\mathbb{P}({\rm SINR}_{0}<\beta) (2)

Because (2) is computed over the distribution of transmitter positions as well as the iid fading coefficients (and consequently transmission powers), it corresponds to fading that occurs on a time-scale that is comparable or slower than the packet duration (if (2) is to correspond roughly to the packet error rate). The outage probability is clearly a continuous increasing function of the intensity λ\lambda.

Define λ(ϵ)\lambda(\epsilon) as the maximum intensity of attempted transmissions such that the outage probability is no larger than ϵ\epsilon, i.e., λ(ϵ)\lambda(\epsilon) is the unique solution of q(λ)=ϵq(\lambda)=\epsilon. The transmission capacity is then defined as c(ϵ)=λ(ϵ)(1ϵ)bc(\epsilon)=\lambda(\epsilon)(1-\epsilon)b, which is the maximum density of successful transmissions times the spectral efficiency bb of each transmission. In other words, transmission capacity is area spectral efficiency subject to an outage constraint.

For the sake of clarity, we define the constants δ=2/α<1\delta=2/\alpha<1 and SNR=pdαη{\rm SNR}=\frac{pd^{-\alpha}}{\eta}. Now consider a path-loss only environment (Hi0=1H_{i0}=1 for all ii) with constant transmission power (Pi=pP_{i}=p for all ii). The main result of [18] is given in the following theorem.

Theorem 1 ([18])

Pure pathloss. Consider a network where the SINR at the reference receiver is given by (2) with Hi0=1H_{i0}=1 and Pi=pP_{i}=p for all ii. Then the following expressions give bounds on the outage probability and transmission attempt intensity for λ,ϵ\lambda,~\epsilon small:

qpl(λ)\displaystyle q^{\rm pl}(\lambda) \displaystyle\geq qlpl(λ)=1exp{λπd2(1β1SNR)δ},\displaystyle q^{\rm pl}_{l}(\lambda)=1-\exp\left\{-\lambda\pi d^{2}\left(\frac{1}{\beta}-\frac{1}{\rm SNR}\right)^{-\delta}\right\}, (3)
λpl(ϵ)\displaystyle\lambda^{\rm pl}(\epsilon) \displaystyle\leq λupl(ϵ)=log(1ϵ)1πd2(1β1SNR)δ.\displaystyle\lambda^{\rm pl}_{u}(\epsilon)=-\log(1-\epsilon)\frac{1}{\pi d^{2}}\left(\frac{1}{\beta}-\frac{1}{\rm SNR}\right)^{\delta}. (4)

Here pl denotes pathloss. The transmission attempt intensity upper bound, λupl(ϵ)\lambda^{\rm pl}_{u}(\epsilon), is obtained by solving qlpl(λ)=ϵq^{\rm pl}_{l}(\lambda)=\epsilon for λ\lambda. These bounds are shown to be approximations for small λ,ϵ\lambda,\epsilon respectively, which is the usual regime of interest. Note also that log(1ϵ)=ϵ+O(ϵ2)-\log(1-\epsilon)=\epsilon+O(\epsilon^{2}), which implies that transmission density is approximately linear with the desired outage level, ϵ\epsilon, for small outages. The following corollary illustrates the simplification of the above results when the noise may be ignored.

Corollary 1

When η=0\eta=0 the expressions in Theorem 1 simplify to:

qpl(λ)\displaystyle q^{\rm pl}(\lambda) \displaystyle\geq qlpl(λ)=1exp{λπd2βδ},\displaystyle q^{\rm pl}_{l}(\lambda)=1-\exp\left\{-\lambda\pi d^{2}\beta^{\delta}\right\}, (5)
λpl(ϵ)\displaystyle\lambda^{\rm pl}(\epsilon) \displaystyle\leq λupl(ϵ)=log(1ϵ)1πd2βδ.\displaystyle\lambda^{\rm pl}_{u}(\epsilon)=-\log(1-\epsilon)\frac{1}{\pi d^{2}\beta^{\delta}}. (6)

III Fractional Power Control

The goal of the paper is to determine the effect that fractional power control has on the outage probability lower bound in (3) and hence the transmission capacity upper bound in (4). We first review the key prior result that we will use, then derive the maximum transmission densities λ\lambda for different power control policies. We conclude the section by finding the optimal power control exponent ss.

III-A Transmission capacity under constant power and channel inversion

In this subsection we restrict our attention to two well-known power control strategies: constant transmit power (or no power control) and channel inversion. Under constant power, Pi=pP_{i}=p for all ii for some common power level pp. Under channel inversion, Pi=p𝔼[H1]Hii1P_{i}=\frac{p}{\mathbb{E}[H^{-1}]}H_{ii}^{-1} for all ii. This means that the received signal power is PiHiidα=p𝔼[H1]dαP_{i}H_{ii}d^{-\alpha}=\frac{p}{\mathbb{E}[H^{-1}]}d^{-\alpha}, which is constant for all ii. That is, channel inversion compensates for the random channel fluctuations between each transmitter and its intended receiver. Moreover, the expected transmission power is 𝔼[Pi]=p\mathbb{E}[P_{i}]=p, so that the constant power and channel inversion schemes use the same expected power. We would like to emphasize the distribution of HH is arbitrary and can be adapted in principle to any relevant fading or compound shadowing-fading model. For some possible distributions (such as Rayleigh fading, i.e. Hexp(1)H\sim\exp(1)), the value 𝔼[H1]\mathbb{E}[H^{-1}] may be undefined, strictly speaking. In practice, the transmit power is finite and so Pi=p𝔼[H1]Hii1P_{i}=\frac{p}{\mathbb{E}[H^{-1}]}H_{ii}^{-1} is finite. The value 𝔼[H1]\mathbb{E}[H^{-1}] is simply a normalizing factor and can be interpreted mathematically to mean that Hmin(H,δ)H\to\min(H,\delta) for an arbitrarily small δ\delta. Such a definition would not affect the results in the paper.

A main result of [16] extended to include thermal noise is given in the following theorem, with a general proof that will apply to all three cases of interest: constant power, channel inversion and fractional power control. Note that cp and ci{\rm ci} are used to denote constant power and channel inversion, respectively.

Theorem 2

Constant power. Consider a network where the SINR at the reference receiver is given by (2) with Pi=pP_{i}=p for all ii. Then the following expressions give good approximations of the outage probability and transmission attempt intensity for λ,ϵ\lambda,\epsilon small.

qcp(λ)\displaystyle q^{\rm cp}(\lambda) \displaystyle\geq qlcp(λ)=1(H00βSNR)𝔼[exp{λπd2𝔼[Hδ](H00β1SNR)δ}|H00βSNR]\displaystyle q^{\rm cp}_{l}(\lambda)=1-\mathbb{P}\left(H_{00}\geq\frac{\beta}{\rm SNR}\right)\mathbb{E}\left[\exp\left\{-\lambda\pi d^{2}\mathbb{E}[H^{\delta}]\left(\frac{H_{00}}{\beta}-\frac{1}{\rm SNR}\right)^{-\delta}\right\}\Big{|}H_{00}\geq\frac{\beta}{\rm SNR}\right]
\displaystyle\approx q~lcp(λ)=1(H00βSNR)exp{λπd2𝔼[Hδ]𝔼[(H00β1SNR)δ|H00βSNR]}\displaystyle\tilde{q}^{\rm cp}_{l}(\lambda)=1-\mathbb{P}\left(H_{00}\geq\frac{\beta}{\rm SNR}\right)\exp\left\{-\lambda\pi d^{2}\mathbb{E}[H^{\delta}]\mathbb{E}\left[\left(\frac{H_{00}}{\beta}-\frac{1}{\rm SNR}\right)^{-\delta}\Big{|}H_{00}\geq\frac{\beta}{\rm SNR}\right]\right\}
λcp(ϵ)\displaystyle\lambda^{\rm cp}(\epsilon) \displaystyle\approx λ~cp(ϵ)=log(1ϵ(H00βSNR))1πd21𝔼[Hδ]𝔼[(H00β1SNR)δ|H00βSNR]1.\displaystyle\tilde{\lambda}^{\rm cp}(\epsilon)=-\log\left(\frac{1-\epsilon}{\mathbb{P}\left(H_{00}\geq\frac{\beta}{\rm SNR}\right)}\right)\frac{1}{\pi d^{2}}\frac{1}{\mathbb{E}[H^{\delta}]}\mathbb{E}\left[\left(\frac{H_{00}}{\beta}-\frac{1}{\rm SNR}\right)^{-\delta}\Big{|}H_{00}\geq\frac{\beta}{\rm SNR}\right]^{-1}. (7)

Channel inversion. Consider the same network with Pi=p𝔼[H1]Hii1P_{i}=\frac{p}{\mathbb{E}[H^{-1}]}H_{ii}^{-1} for all ii. Then the following expressions give tight bounds on the outage probability and transmission attempt intensity for λ,ϵ\lambda,\epsilon small:

qci(λ)\displaystyle q^{\rm ci}(\lambda) \displaystyle\geq qlci(λ)=1exp{λπd2𝔼[Hδ]𝔼[Hδ](1β𝔼[H1]SNR)δ}\displaystyle q^{\rm ci}_{l}(\lambda)=1-\exp\left\{-\lambda\pi d^{2}\mathbb{E}[H^{\delta}]\mathbb{E}[H^{-\delta}]\left(\frac{1}{\beta}-\frac{\mathbb{E}[H^{-1}]}{\rm SNR}\right)^{-\delta}\right\} (8)
λci(ϵ)\displaystyle\lambda^{\rm ci}(\epsilon) \displaystyle\leq λuci(ϵ)=log(1ϵ)1πd21𝔼[Hδ]𝔼[Hδ](1β𝔼[H1]SNR)δ.\displaystyle\lambda^{\rm ci}_{u}(\epsilon)=-\log(1-\epsilon)\frac{1}{\pi d^{2}}\frac{1}{\mathbb{E}[H^{\delta}]\mathbb{E}[H^{-\delta}]}\left(\frac{1}{\beta}-\frac{\mathbb{E}[H^{-1}]}{\rm SNR}\right)^{\delta}. (9)
Proof:

The SINR at the reference receiver for a generic power vector {Pi}\{P_{i}\} is

SINR0=P0H00dαiΠ(λ)PiHi0Xiα+η,{\rm SINR}_{0}=\frac{P_{0}H_{00}d^{-\alpha}}{\sum_{i\in\Pi(\lambda)}P_{i}H_{i0}X_{i}^{-\alpha}+\eta}, (10)

and the corresponding outage probability is

q(λ)=(SINR0<β)=(P0H00dαiΠ(λ)PiHi0Xiα+η<β).q(\lambda)=\mathbb{P}({\rm SINR}_{0}<\beta)=\mathbb{P}\left(\frac{P_{0}H_{00}d^{-\alpha}}{\sum_{i\in\Pi(\lambda)}P_{i}H_{i0}X_{i}^{-\alpha}+\eta}<\beta\right). (11)

Rearranging yields:

q(λ)=(iΠ(λ)PiHi0XiαP0H00dαβη).q(\lambda)=\mathbb{P}\left(\sum_{i\in\Pi(\lambda)}P_{i}H_{i0}X_{i}^{-\alpha}\geq\frac{P_{0}H_{00}d^{-\alpha}}{\beta}-\eta\right). (12)

Note that outage is certain when P0H00<ηβdαP_{0}H_{00}<\eta\beta d^{\alpha}. Conditioning on P0H00P_{0}H_{00} and using f()f(\cdot) to denote the density of P0H00P_{0}H_{00} yields:

q(λ)=(P0H00ηβdα)+ηβdα(iΠ(λ)PiHi0Xiαp0h00βdαη|P0H00=p0h00)f(p0h00)d(p0h00).q(\lambda)=\mathbb{P}\left(P_{0}H_{00}\leq\eta\beta d^{\alpha}\right)+\int_{\eta\beta d^{\alpha}}^{\infty}\mathbb{P}\left(\sum_{i\in\Pi(\lambda)}P_{i}H_{i0}X_{i}^{-\alpha}\geq\frac{p_{0}h_{00}}{\beta d^{\alpha}}-\eta~\Big{|}~P_{0}H_{00}=p_{0}h_{00}\right)f(p_{0}h_{00}){\rm d}(p_{0}h_{00}). (13)

Recall the generic lower bound from [16]: if Π(λ)={(Xi,Zi)}\Pi(\lambda)=\{(X_{i},Z_{i})\} is a homogeneous marked Poisson point process with points {Xi}\{X_{i}\} of intensity λ\lambda and iid marks {Zi}\{Z_{i}\} independent of the {Xi}\{X_{i}\}, then

(iΠ(λ)ZiXiα>y)1exp{πλ𝔼[Zδ]yδ},\mathbb{P}\left(\sum_{i\in\Pi(\lambda)}Z_{i}X_{i}^{-\alpha}>y\right)\geq 1-\exp\left\{-\pi\lambda\mathbb{E}[Z^{\delta}]y^{-\delta}\right\}, (14)

Applying here with Zi=PiHi0Z_{i}=P_{i}H_{i0} and y=p0h00βdαηy=\frac{p_{0}h_{00}}{\beta d^{\alpha}}-\eta:

q(λ)\displaystyle q(\lambda) \displaystyle\geq (P0H00ηβdα)+ηβdα(1exp{πλ𝔼[(PiHi0)δ](p0h00βdαη)δ})f(p0h00)d(p0h00)\displaystyle\mathbb{P}\left(P_{0}H_{00}\leq\eta\beta d^{\alpha}\right)+\int_{\eta\beta d^{\alpha}}^{\infty}\left(1-\exp\left\{-\pi\lambda\mathbb{E}[(P_{i}H_{i0})^{\delta}]\left(\frac{p_{0}h_{00}}{\beta d^{\alpha}}-\eta\right)^{-\delta}\right\}\right)f(p_{0}h_{00}){\rm d}(p_{0}h_{00}) (15)
=\displaystyle= 1ηβdαexp{πλ𝔼[(PiHi0)δ](p0h00βdαη)δ}f(p0h00)d(p0h00)\displaystyle 1-\int_{\eta\beta d^{\alpha}}^{\infty}\exp\left\{-\pi\lambda\mathbb{E}[(P_{i}H_{i0})^{\delta}]\left(\frac{p_{0}h_{00}}{\beta d^{\alpha}}-\eta\right)^{-\delta}\right\}f(p_{0}h_{00}){\rm d}(p_{0}h_{00})
=\displaystyle= 1(P0H00ηβdα)𝔼[exp{λπd2𝔼[(PiHi0)δ](P0H00βηdα)δ}|P0H00ηβdα].\displaystyle 1-\mathbb{P}\left(P_{0}H_{00}\geq\eta\beta d^{\alpha}\right)\mathbb{E}\left[\exp\left\{-\lambda\pi d^{2}\mathbb{E}[(P_{i}H_{i0})^{\delta}]\left(\frac{P_{0}H_{00}}{\beta}-\frac{\eta}{d^{-\alpha}}\right)^{-\delta}\right\}\Big{|}P_{0}H_{00}\geq\eta\beta d^{\alpha}\right].

The Jensen approximation for this quantity is:

q(λ)1(P0H00ηβdα)exp{λπd2𝔼[(PiHi0)δ]𝔼[(P0H00βηdα)δ|P0H00ηβdα]}.q(\lambda)\approx 1-\mathbb{P}\left(P_{0}H_{00}\geq\eta\beta d^{\alpha}\right)\exp\left\{-\lambda\pi d^{2}\mathbb{E}[(P_{i}H_{i0})^{\delta}]\mathbb{E}\left[\left(\frac{P_{0}H_{00}}{\beta}-\frac{\eta}{d^{-\alpha}}\right)^{-\delta}\Big{|}P_{0}H_{00}\geq\eta\beta d^{\alpha}\right]\right\}. (16)

For constant power we substitute PiHi0=pHi0P_{i}H_{i0}=pH_{i0} (for all ii) into (15) and (16) and manipulate to get the expressions for qlcp(λ)q^{\rm cp}_{l}(\lambda) and q~lcp(λ)\tilde{q}^{\rm cp}_{l}(\lambda) in (7). To obtain λ~cp(ϵ)\tilde{\lambda}^{\rm cp}(\epsilon), we solve q~lcp(λ)=ϵ\tilde{q}^{\rm cp}_{l}(\lambda)=\epsilon for λ\lambda. For channel inversion, P0H00=p𝔼[H1]P_{0}H_{00}=\frac{p}{\mathbb{E}[H^{-1}]} while for i0i\neq 0 we have PiHi0=p𝔼[H1]Hi0HiiP_{i}H_{i0}=\frac{p}{\mathbb{E}[H^{-1}]}\frac{H_{i0}}{H_{ii}}. Plugging into (15) and using the fact that HiiH_{ii} and Hi0H_{i0} are i.i.d. yields (8), and (9) is simply the inverse of (8). ∎

Note that channel inversion only makes sense when SNR𝔼[H1]=pdαη𝔼[H1]\frac{\rm SNR}{\mathbb{E}[H^{-1}]}=\frac{pd^{-\alpha}}{\eta\mathbb{E}[H^{-1}]}, the effective interference-free SNR after taking into account the power cost of inversion, is larger than the SINR threshold β\beta. The validity of the outage lower bound/density upper bound as well as of the Jensen’s approximation are evaluated in the numerical and simulation results in Section IV.

When the thermal noise can be ignored, these results simplify to the expressions given in the following corollary:

Corollary 2

When η=0\eta=0 the expressions in Theorem 2 simplify to:

qcp(λ)\displaystyle q^{\rm cp}(\lambda) \displaystyle\geq qlcp(λ)=1𝔼[exp{λπd2βδ𝔼[Hδ]H00δ}]\displaystyle q^{\rm cp}_{l}(\lambda)=1-\mathbb{E}\left[\exp\left\{-\lambda\pi d^{2}\beta^{\delta}\mathbb{E}\left[H^{\delta}\right]H_{00}^{-\delta}\right\}\right]
\displaystyle\approx q~lcp(λ)=1exp{λπd2βδ𝔼[Hδ]𝔼[Hδ]},\displaystyle\tilde{q}^{\rm cp}_{l}(\lambda)=1-\exp\left\{-\lambda\pi d^{2}\beta^{\delta}\mathbb{E}\left[H^{\delta}\right]\mathbb{E}\left[H^{-\delta}\right]\right\},
qci(λ)\displaystyle q^{\rm ci}(\lambda) \displaystyle\geq qlci(λ)=1exp{λπd2βδ𝔼[Hδ]𝔼[Hδ]},\displaystyle q^{\rm ci}_{l}(\lambda)=1-\exp\left\{-\lambda\pi d^{2}\beta^{\delta}\mathbb{E}\left[H^{\delta}\right]\mathbb{E}\left[H^{-\delta}\right]\right\},
λcp(ϵ)\displaystyle\lambda^{\rm cp}(\epsilon) \displaystyle\approx λ~cp(ϵ)=log(1ϵ)1πd2βδ1𝔼[Hδ]𝔼[Hδ],\displaystyle\tilde{\lambda}^{\rm cp}(\epsilon)=-\log(1-\epsilon)\frac{1}{\pi d^{2}\beta^{\delta}}\frac{1}{\mathbb{E}\left[H^{\delta}\right]\mathbb{E}\left[H^{-\delta}\right]},
λci(ϵ)\displaystyle\lambda^{\rm ci}(\epsilon) \displaystyle\leq λuci(ϵ)=log(1ϵ)1πd2βδ1𝔼[Hδ]𝔼[Hδ].\displaystyle\lambda^{\rm ci}_{u}(\epsilon)=-\log(1-\epsilon)\frac{1}{\pi d^{2}\beta^{\delta}}\frac{1}{\mathbb{E}\left[H^{\delta}\right]\mathbb{E}\left[H^{-\delta}\right]}. (17)

Note that these expressions match Theorem 3 and Corollary 3 of the SIR-analysis performed in [16].

In the absence of noise the constant power outage probability approximation equals the channel inversion outage probability lower bound: q~lcp(λ)=qlci(λ)\tilde{q}_{l}^{\rm cp}(\lambda)=q_{l}^{\rm ci}(\lambda). As a result, the constant power transmission attempt intensity approximation equals the channel inversion transmission attempt intensity upper bound: λ~cp(ϵ)=λuci(ϵ)\tilde{\lambda}^{\rm cp}(\epsilon)=\lambda_{u}^{\rm ci}(\epsilon). Comparing λ~cp(ϵ)=λuci(ϵ)\tilde{\lambda}^{\rm cp}(\epsilon)=\lambda_{u}^{\rm ci}(\epsilon) in (2) with λupl(ϵ)\lambda_{u}^{\rm pl}(\epsilon) in (6) it is evident that the impact of fading on the transmission capacity is measured by the loss factor, Lcp=LciL^{\rm cp}=L^{\rm ci}, defined as

Lcp=Lci=1𝔼[Hδ]𝔼[Hδ]<1.L^{\rm cp}=L^{\rm ci}=\frac{1}{\mathbb{E}\left[H^{\delta}\right]\mathbb{E}\left[H^{-\delta}\right]}<1. (18)

The inequality is obtained by applying Jensen’s inequality to the convex function 1/x1/x and the random variable HδH^{\delta}. If constant power is used, the 𝔼[Hδ]\mathbb{E}[H^{-\delta}] term is due to fading of the desired signal while the 𝔼[Hδ]\mathbb{E}[H^{\delta}] term is due to fading of the interfering links. Fading of the interfering signal has a positive effect while fading of the desired signal has a negative effect. If channel inversion is performed the 𝔼[Hδ]\mathbb{E}[H^{-\delta}] term is due to each interfering transmitter using power proportional to Hii1H_{ii}^{-1}. When the path loss exponent, α\alpha, is close to 2 then δ=2/α\delta=2/\alpha is close to one, so the term 𝔼[Hδ]\mathbb{E}[H^{-\delta}] is nearly equal to the expectation of the inverse of the fading, which can be extremely large for severe fading distributions such as Rayleigh. As a less severe example, α=3\alpha=3, the loss factor for Rayleigh fading is Lcp=Lci=0.41L^{\mathrm{cp}}=L^{\mathrm{ci}}=0.41.

III-B Transmission capacity under fractional power control

In this section we generalize the results of Theorem 2 by introducing fractional power control (FPC) with parameter s[0,1]s\in[0,1]. Under FPC the transmission power is set to Pi=p𝔼[Hs]HiisP_{i}=\frac{p}{\mathbb{E}[H^{-s}]}H_{ii}^{-s} for each ii. The received power at receiver ii is then PiHiidα=p𝔼[Hs]Hii1sdαP_{i}H_{ii}d^{-\alpha}=\frac{p}{\mathbb{E}[H^{-s}]}H_{ii}^{1-s}d^{-\alpha}, which depends upon ii aside from s=1s=1. The expected transmission power is pp, ensuring a fair comparison with the results in Theorems 1 and 2. Note that constant power corresponds to s=0s=0 and channel inversion corresponds to s=1s=1. The following theorem gives good approximations on the outage probability and maximum allowable transmission intensity under FPC.

Theorem 3

Fractional power control. Consider a network where the SINR at the reference receiver is given by (2) with Pi=p𝔼[Hs]HiisP_{i}=\frac{p}{\mathbb{E}[H^{-s}]}H_{ii}^{-s} for all ii, for some s[0,1]s\in[0,1]. Then the following expressions give good approximations of the outage probability and maximum transmission attempt intensity for λ,ϵ\lambda,~\epsilon small

qfpc(λ)\displaystyle q^{\rm fpc}(\lambda) \displaystyle\geq qlfpc(λ)=1(H00κ(s))×\displaystyle q^{\rm fpc}_{l}(\lambda)=1-\mathbb{P}\left(H_{00}\geq\kappa(s)\right)\times
𝔼[exp{λπd2𝔼[Hsδ]𝔼[Hδ](H001sβ𝔼[Hs]SNR)δ}|H00κ(s)]\displaystyle\hskip 71.13188pt\mathbb{E}\left[\exp\left\{-\lambda\pi d^{2}\mathbb{E}[H^{-s\delta}]\mathbb{E}[H^{\delta}]\left(\frac{H_{00}^{1-s}}{\beta}-\frac{\mathbb{E}[H^{-s}]}{\rm SNR}\right)^{-\delta}\right\}\Big{|}H_{00}\geq\kappa(s)\right]
\displaystyle\approx q~lfpc(λ)=1(H00κ(s))×\displaystyle\tilde{q}^{\rm fpc}_{l}(\lambda)=1-\mathbb{P}\left(H_{00}\geq\kappa(s)\right)\times
exp{λπd2𝔼[Hsδ]𝔼[Hδ]𝔼[(H001sβ𝔼[Hs]SNR)δ|H00κ(s)]}\displaystyle\hskip 71.13188pt\exp\left\{-\lambda\pi d^{2}\mathbb{E}[H^{-s\delta}]\mathbb{E}[H^{\delta}]\mathbb{E}\left[\left(\frac{H_{00}^{1-s}}{\beta}-\frac{\mathbb{E}[H^{-s}]}{\rm SNR}\right)^{-\delta}\Big{|}H_{00}\geq\kappa(s)\right]\right\}
λfpc(ϵ)\displaystyle\lambda^{\rm fpc}(\epsilon) \displaystyle\approx λ~fpc(ϵ)=log(1ϵ(H00κ(s)))1πd21𝔼[Hsδ]𝔼[Hδ]×\displaystyle\tilde{\lambda}^{\rm fpc}(\epsilon)=-\log\left(\frac{1-\epsilon}{\mathbb{P}\left(H_{00}\geq\kappa(s)\right)}\right)\frac{1}{\pi d^{2}}\frac{1}{\mathbb{E}[H^{-s\delta}]\mathbb{E}[H^{\delta}]}\times
(𝔼[(H001sβ𝔼[Hs]SNR)δ|H00κ(s)])1\displaystyle\hskip 184.9429pt\left(\mathbb{E}\left[\left(\frac{H_{00}^{1-s}}{\beta}-\frac{\mathbb{E}[H^{-s}]}{\rm SNR}\right)^{-\delta}\Big{|}H_{00}\geq\kappa(s)\right]\right)^{-1}

where κ(s)=(βSNR𝔼[Hs])11s\kappa(s)=\left(\frac{\beta}{\rm SNR}\mathbb{E}[H^{-s}]\right)^{\frac{1}{1-s}}.

Proof:

Under FPC, the transmit power for each user is constructed as Pi=p𝔼[Hs]HiisP_{i}=\frac{p}{\mathbb{E}[H^{-s}]}H_{ii}^{-s}. Substituting this value into the proof for Theorem 2 immediately gives the expression for qlfpc(λ)q^{\rm fpc}_{l}(\lambda). Again, the transmission attempt intensity approximation is obtained by solving q~l(λ)=ϵ\tilde{q}_{l}(\lambda)=\epsilon for λ\lambda. ∎

As with Theorem 2, the approximation qlfpc(λ)q~lfpc(λ)q^{\rm fpc}_{l}(\lambda)\approx\tilde{q}^{\rm fpc}_{l}(\lambda) is accurate when the exponential term in qlfpc(λ)q^{\rm fpc}_{l}(\lambda) is approximately linear in its argument and thus Jensen’s is tight. In other words, this approximation utilizes the fact that ex{\rm e}^{-x} is nearly linear for small xx. Looking at the expression for qlfpc(λ)q^{\rm fpc}_{l}(\lambda) we see that this reasonable when the relative density λπd2\lambda\pi d^{2} is small. If this is not true then the approximation q~lfpc(λ)\tilde{q}^{\rm fpc}_{l}(\lambda) is not sufficiently accurate, as will be further seen in the numerical results presented in Section IV. The FPC transmission attempt intensity approximation, λ~fpc(ϵ)\tilde{\lambda}^{\rm fpc}(\epsilon), is obtained by solving q~lfpc(λ)=ϵ\tilde{q}^{\rm fpc}_{l}(\lambda)=\epsilon for λ\lambda. The following corollary illustrates the simplification of the above results when the noise may be ignored.

Corollary 3

When η=0\eta=0 the expressions in Theorem 3 simplify to:

qfpc(λ)\displaystyle q^{\rm fpc}(\lambda) \displaystyle\geq qlfpc(λ)=1𝔼[exp{λπd2βδ𝔼[Hδ]𝔼[Hsδ]H00(1s)δ}]\displaystyle q^{\rm fpc}_{l}(\lambda)=1-\mathbb{E}\left[\exp\left\{-\lambda\pi d^{2}\beta^{\delta}\mathbb{E}\left[H^{\delta}\right]\mathbb{E}\left[H^{-s\delta}\right]H_{00}^{-(1-s)\delta}\right\}\right]
\displaystyle\approx q~lfpc(λ)=1exp{λπd2βδ𝔼[Hδ]𝔼[Hsδ]𝔼[H(1s)δ]},\displaystyle\tilde{q}^{\rm fpc}_{l}(\lambda)=1-\exp\left\{-\lambda\pi d^{2}\beta^{\delta}\mathbb{E}\left[H^{\delta}\right]\mathbb{E}\left[H^{-s\delta}\right]\mathbb{E}\left[H^{-(1-s)\delta}\right]\right\},
λfpc(ϵ)\displaystyle\lambda^{\rm fpc}(\epsilon) \displaystyle\approx λ~fpc(ϵ)=log(1ϵ)1πd2βδ1𝔼[Hδ]𝔼[Hsδ]𝔼[H(1s)δ].\displaystyle\tilde{\lambda}^{\rm fpc}(\epsilon)=-\log(1-\epsilon)\frac{1}{\pi d^{2}\beta^{\delta}}\frac{1}{\mathbb{E}\left[H^{\delta}\right]\mathbb{E}\left[H^{-s\delta}\right]\mathbb{E}\left[H^{-(1-s)\delta}\right]}. (19)

The loss factor for FPC, LfpcL^{\rm fpc}, is the reduction in the transmission capacity approximation relative to the pure pathloss case:

Lfpc(s)=1𝔼[Hδ]𝔼[Hsδ]𝔼[H(1s)δ].L^{\rm fpc}(s)=\frac{1}{\mathbb{E}\left[H^{\delta}\right]\mathbb{E}\left[H^{-s\delta}\right]\mathbb{E}\left[H^{-(1-s)\delta}\right]}. (20)

Clearly, the loss factor LfpcL^{\mathrm{fpc}} for FPC depends on the design choice of the exponent ss.

III-C Optimal Fractional Power Control Exponent

Fractional power control represents a balance between the extremes of no power control and channel inversion. The mathematical effect of fractional power control is to replace the 𝔼[Hδ]\mathbb{E}[H^{-\delta}] term with 𝔼[Hsδ]𝔼[H(1s)δ]\mathbb{E}[H^{-s\delta}]\mathbb{E}[H^{-(1-s)\delta}]. This is because the signal fading is softened by the power control exponent s-s so that it results in a leading term of H(1s)H^{-(1-s)} (rather than H1H^{-1}) in the numerator of the SINR expression, and ultimately to the 𝔼[H(1s)δ]\mathbb{E}[H^{-(1-s)\delta}] term. The interference power is also softened by the fractional power control and leads to the 𝔼[Hsδ]\mathbb{E}[H^{-s\delta}] term.

The key question of course lies in determining the optimal power control exponent. Although it does not seem possible to derive an analytical expression for the exponent that minimizes the general expression for qlfpc(λ)q^{\rm fpc}_{l}(\lambda) given in Theorem 3, we can find the exponent that minimizes the outage probability approximation in the case of no noise.

Theorem 4

In the absence of noise (η=0\eta=0), the fractional power control outage probability approximation, q~lfpc(λ)\tilde{q}^{\rm fpc}_{l}(\lambda), is minimized for s=12s=\frac{1}{2}. Hence, the fractional power control transmission attempt intensity approximation, λ~fpc(ϵ)\tilde{\lambda}^{\rm fpc}(\epsilon) is also maximized for s=12s=\frac{1}{2}.

Proof:

Because the outage probability/transmission density approximations depend on the exponent ss only through the quantity 𝔼[Hsδ]𝔼[H(1s)δ]\mathbb{E}\left[H^{-s\delta}\right]\mathbb{E}\left[H^{-(1-s)\delta}\right], it is sufficient to show that 𝔼[Hsδ]𝔼[H(1s)δ]\mathbb{E}\left[H^{-s\delta}\right]\mathbb{E}\left[H^{-(1-s)\delta}\right] is minimized at s=12s=\frac{1}{2}. To do this, we use the following general result, which we prove in the Appendix. For any non-negative random variable XX, the function

h(s)=𝔼[Xs]𝔼[Xs1],h(s)=\mathbb{E}\left[X^{-s}\right]\mathbb{E}\left[X^{s-1}\right], (21)

is convex in ss for ss\in\mathbb{R} with a unique minimum at s=12s=\frac{1}{2}. Applying this result to random variable X=HδX=H^{\delta} gives the desired result. ∎

The theorem shows that transmission density is maximized, or equivalently, outage probability is minimized, by balancing the positive and negative effects of power control, which are reduction of signal fading and increasing interference, respectively. Using an exponent greater than 12\frac{1}{2} over-compensates for signal fading and leads to interference levels that are too high, while using an exponent smaller than 12\frac{1}{2} leads to small interference levels but an under-compensation for signal fading. Note that because the key expression 𝔼[Hsδ]𝔼[H(1s)δ]\mathbb{E}\left[H^{-s\delta}\right]\mathbb{E}\left[H^{-(1-s)\delta}\right] is convex, the loss relative to using s=12s=\frac{1}{2} increases monotonically both as s0s\to 0 and s1s\to 1.

One can certainly envision “fractional” power control schemes that go even further. For example, s>1s>1 corresponds to “super” channel inversion, in which bad channels take resources from good channels even more so than in normal channel inversion. Not surprisingly, this is not a wise policy. Less obviously, s<0s<0 corresponds to what is sometimes called “greedy” optimization, in which good channels are given more resources at the further expense of poor channels. Waterfilling is an example of a greedy optimization procedure. But, since 𝔼[Hsδ]𝔼[H(1s)δ]\mathbb{E}\left[H^{-s\delta}\right]\mathbb{E}\left[H^{-(1-s)\delta}\right] monotonically increases as ss decreases, it is clear that greedy power allocations of any type are worse than even constant transmit power under the SINR-target set up.

The numerical results in the next section show that FPC is very beneficial relative to constant transmit power or channel inversion. However, fading has a deleterious effect relative to no fading even if the optimal exponent is used. To see this, note that x12x^{-\frac{1}{2}} is a convex function and therefore Jensen’s yields 𝔼[X12](𝔼[X])12\mathbb{E}[X^{-\frac{1}{2}}]\geq(\mathbb{E}[X])^{-\frac{1}{2}} for any non-negative random variable XX. Applying this to X=HδX=H^{\delta} we get (𝔼[Hδ2])2(𝔼[Hδ])1\left(\mathbb{E}\left[H^{-\frac{\delta}{2}}\right]\right)^{2}\geq\left(\mathbb{E}[H^{\delta}]\right)^{-1}, which implies

Lfpc(1/2)=1𝔼[Hδ](𝔼[Hδ2])21.\displaystyle L^{\rm fpc}(1/2)=\frac{1}{\mathbb{E}\left[H^{\delta}\right]\left(\mathbb{E}\left[H^{-\frac{\delta}{2}}\right]\right)^{2}}\leq 1.

Therefore, fractional PC cannot fully overcome fading, but it is definitely a better power control policy than constant power transmission or traditional power control (channel inversion).

IV Numerical Results and Discussion

In this section, the implications of fractional power control are illustrated through numerical plots and analytical discussion. The tightness of the bounds will be considered as a function of the system parameters, and the choice of a robust FPC exponent ss will be proposed. As default parameters, the simulations assume

α=3,β=1(0dB),d=10m,SNR=pdαη=100(20dB),λ=0.0001usersm2.\begin{array}[]{cccccc}\alpha=3,&\beta=1~(0~{\rm dB}),&d=10{\rm m},&{\rm SNR}=\frac{pd^{-\alpha}}{\eta}=100~(20~{\rm dB}),&\lambda=0.0001~\frac{\rm users}{{\rm m}^{2}}.\end{array} (22)

Furthermore, Rayleigh fading is assumed for the numerical results.

IV-A Effect of Fading

The benefit of fractional power control can be quickly illustrated in Rayleigh fading, in which case the channel power HH is exponentially distributed and the moment generating function is therefore

𝔼[Ht]=Γ(1+t),\mathbb{E}[H^{t}]=\Gamma(1+t), (23)

where Γ()\Gamma(\cdot) is the standard gamma function. If fractional power control is used, the transmission capacity loss due to fading is

Lfpc=1𝔼[Hδ]𝔼[Hsδ]𝔼[H(1s)δ]=1Γ(1+δ)Γ(1sδ)Γ(1(1s)δ)\displaystyle L^{\rm fpc}=\frac{1}{\mathbb{E}\left[H^{\delta}\right]\mathbb{E}\left[H^{-s\delta}\right]\mathbb{E}\left[H^{(1-s)\delta}\right]}=\frac{1}{\Gamma(1+\delta)\cdot\Gamma(1-s\delta)\cdot\Gamma(1-(1-s)\delta)} (24)

In Fig. 1 this loss factor (LL) is plotted as a function of ss for path loss exponents α={2.1,3,4}\alpha=\{2.1,3,4\}. Notice that for each value of α\alpha the maximum takes place at s=12s=\frac{1}{2}, and that the cost of not using fractional power control is highest for small path loss exponents because Γ(1+x)\Gamma(1+x) goes to infinity quite steeply as x1x\rightarrow-1. This plot implies that in severe fading channels, the gain from FPC can be quite significant.

It should be noted that the expression in (24) is for the case of no thermal noise (η=0\eta=0). In this case the power cost of FPC completely vanishes, because the same power normalization (by 𝔼[Hs]\mathbb{E}[H^{-s}]) is performed by each transmitting node and therefore this normalization cancels in the SIR expression. On the other hand, this power cost does not vanish if the noise is strictly positive and can potentially be quite significant, particularly if SNR\rm SNR is not large. A simple application of Jensen’s shows that the power normalization factor 𝔼[Hs]\mathbb{E}[H^{-s}] is an increasing function of the exponent ss for any distribution on HH. For the particular case of Rayleigh fading this normalization factor is Γ(1s)\Gamma(1-s) which makes it prohibitively expensive to choose ss very close to one; indeed, the choice s=1s=1 requires infinite power and thus is not feasible. On the other hand, note that Γ(.5)\Gamma(.5) is approximately 2.52.5 dB and thus the cost of a moderate exponent is not so large. When the interference-free SNR\rm SNR is reasonably large, this normalization factor is relatively negligible and the effect of FPC is well approximated by (24).

IV-B Tightness of Bounds

There are two principle approximations made in attaining the expressions for outage probability and transmission capacity in Theorem 3. First, the inequality is due to considering only dominant interferers; that is, an interferer whose channel to the desired receiver is strong enough to cause outage even without any other interferers present. This is a lower bound on outage since it ignores non-dominant interferers, but nevertheless has been seen to be quite accurate in our prior work [18, 20, 16]. Second, Jensen’s inequality is used to bound 𝔼[exp(X)]exp(𝔼[X])\mathbb{E}[\exp(X)]\geq\exp(\mathbb{E}[X]) in the opposite direction, so this results in an approximation to the outage probability rather than a lower bound; numerical results confirm that this approximation is in fact not a lower bound in general. Therefore, we consider the three relevant quantities: (1) the actual outage probability qfpc(λ)q^{\rm fpc}(\lambda), which is determined via Monte-Carlo simulation and does not depend on any bounds or approximations, (2) a numerical computation of the outage probability lower bound qlfpc(λ)q^{\rm fpc}_{l}(\lambda), and (3) the approximation to the outage probability q~lfpc(λ)\tilde{q}^{\rm fpc}_{l}(\lambda) reached by applying Jensen’s inequality to qlfpc(λ)q^{\rm fpc}_{l}(\lambda). Note that because of the two opposing bounds (one lower and one upper), we cannot say a priori that method (2) will produce more accurate expressions than method (3).

The tightness of the bounds is explored in Figs. 2 - 5. Consider first Fig. 2 for the default parameters given above. We can see that the lower bound and the Jensen approximation both reasonably approximate the simulation results, and the approximation winds up serving as a lower bound as well. The Jensen’s approximation is very accurate for large values of ss (i.e., closer to channel inversion), and while looser for smaller values of ss, this “error” actually moves the Jensen’s approximation closer to the actual (simulated) outage probability. The Jensen’s approximation approaches the lower bound as s1s\to 1 because the random variable H(1s)δH^{(1-s)\delta} approaches a constant, where Jensen’s inequality trivially holds with equality (see, e.g., (3)). Changing the path loss exponent α\alpha, the SNR, the target SINR β\beta, or the density λ\lambda can have a significant effect on the bounds, as we will see. With the important exception of high density networks, the approximations are seen to be reasonably accurate for reasonable parameter values.

Path loss. In Fig. 3, the bounds are given for α=2.2\alpha=2.2 and α=5\alpha=5, which correspond to much weaker and much stronger attenuation than the (more likely) default case of α=3\alpha=3. For weaker attenuation, we can see that the lower bound holds the right shape but is less accurate, while the Jensen’s approximation becomes very loose when the FPC exponent ss is small. For path loss exponents near 22, the dominant interferer approximation is weakened because the attenuation of non-dominant interferers is less drastic. On the other hand, both the lower bound and Jensen’s approximation are very accurate in strong attenuation environments as seen in the α=5\alpha=5 plot. This is because the dominant interferer approximation is very reasonable in such cases.

SNR. The behavior of the bounds also varies as the background noise level changes, as shown in Fig. 4. When the SNR is 10 dB, the bounds are quite tight. However, the behavior of outage probability as a function of ss is quite different from the default case in Fig. 2: outage probability decreases slowly as ss is increased, and a rather sharp jump is seen as ss approaches one. When the interference-free SNR is only moderately larger than the target SINR (in this case there is a 10 dB difference between SNR\rm SNR and β\beta), a significant portion of outages occur because the signal power is so small that the interference-free received SNR falls below the target β\beta; this probability is captured by the (H00κ(s))\mathbb{P}\left(H_{00}\geq\kappa(s)\right) terms in Theorem 3. On the other hand, if SNR\rm SNR is much larger than the target β\beta, outages are almost always due to a combination of signal fading and large interference power rather than to signal fading alone (i.e., (H00κ(s))\mathbb{P}\left(H_{00}\geq\kappa(s)\right) is insignificant compared to the total outage probability). When outages caused purely by signal fading are significant, the dependence on the exponent ss is significantly reduced. Furthermore, the power cost of FPC becomes much more significant when the gap between SNR\rm SNR and β\beta is reduced; this explains the sharp increase in outage as ss approaches one. When SNR=30\rm SNR=30 dB, the behavior is quite similar to the 20 dB case because at this point the gap between SNR\rm SNR and β\beta is so large that thermal noise can effectively be neglected.

Target SINR. A default SINR of β=1\beta=1 was chosen, which corresponds roughly to a spectral efficiency of 1 bps/Hz with strong coding, and lies between the low and high SINR regimes. Exploring an order of magnitude above and below the default in Fig. 5, we see that for β=0.1\beta=0.1 the bounds are highly accurate, and show that s=12s^{*}=\frac{1}{2} is a good choice. For this choice of parameters there is a 30 dB gap between SNR\rm SNR and β\beta and thus thermal noise is essentially negligible. On the other hand, if β=10\beta=10 the bounds are still reasonable, but the outage behavior is very similar to the earlier case where SNR=10\rm SNR=10 dB and β=0\beta=0 dB because there is again only a 10 dB gap between SNR\rm SNR and β\beta. Despite the qualitative and quantitative differences for low SNR and high target SINR from the default values, it is interesting to note that in both cases s=12s=\frac{1}{2} is still a robust choice for the FPC exponent.

Density. The default value of λ=0.0001\lambda=0.0001 corresponds to a somewhat low density network because the expected distance to the nearest interferer is approximately 5050 m, while the TX-RX distance is d=10d=10 m. In Fig. 6 we explore a density an order of magnitude lower and higher than the default value. When the network is even sparser, the bounds are extremely accurate and we see that s=12s^{*}=\frac{1}{2} is a near-optimal choice. However, the behavior with ss is very different in a dense network where λ=.001\lambda=.001 and the nearest interferer is approximately 1717 m away. In such a network we see that the nearest neighbor bound is quite loose because a substantial fraction of outages are caused by the summation of non-dominant interferers, as intuitively expected for a dense network. Although the bound is loose, it does capture the fact that outage increases with the exponent ss. On the other hand, the Jensen approximation is loose and does not correctly capture the relationship between ss and outage. The approximation is based on the fact that the function exe^{-x} is approximately linear for small xx. The quantity xx is proportional to πλd2\pi\lambda d^{2}, which is large when the network is dense relative to TX-RX distance dd, and thus this approximation is not valid for relatively dense networks.

IV-C Choosing the FPC exponent ss

Determining the optimum choice of FPC exponent ss is a key interest of this paper. As seen in Sect. III-C, s=12s^{*}=\frac{1}{2} is optimal for the Jensen’s approximation and with no noise, both of which are questionable assumptions in many regimes of interest. In Figs. 710, we plot the truly optimal choice of ss^{*} for the default parameters, while varying α\alpha, SNR, β\beta, and λ\lambda, respectively. That is, the value of ss that minimizes the true outage probability is determined for each set of parameters. The FPC exponents sl(Δ)s_{l}(\Delta) and su(Δ)s_{u}(\Delta) are also plotted, which provide Δ\Delta% error below and above the optimum outage probability. For the plots, we let Δ=1\Delta=1 and Δ=10\Delta=10.

The key findings are: (1) In the pathloss (α\alpha) plot, s=12s^{*}=\frac{1}{2} is a very robust choice for all attenuation regimes; (2) For SNR, s=12s^{*}=\frac{1}{2} is only robust at high SNR, and at low SNR constant transmit power is preferable; (3) For target SINR β\beta, s=12s^{*}=\frac{1}{2} is robust at low and moderate SINR targets (i.e. low to moderate data rates), but for high SINR targets constant transmit power is preferred; (4) For density λ\lambda, s=12s^{*}=\frac{1}{2} is robust at low densities, but constant transmit power is preferred at high densities.

The explanation for findings (2) and (3) is due to the dependence of outage behavior on the difference between SNR\rm SNR and β\beta. As seen earlier, thermal noise is essentially negligible when this gap is larger than approximately 20 dB. As a result, it is reasonable that the exponent shown to be optimal for noise-free networks (s=12s=\frac{1}{2}) would be near-optimal for networks with very low levels of thermal noise. On the other hand, outage probability behaves quite differently when SNR\rm SNR is only slightly larger than β\beta. In this case, power is very valuable and it is not worth incurring the normalization cost of FPC and thus very small FPC exponents are optimal. Intuitively, achieving high data rates in moderate SNR or moderate data rates in low SNR are difficult objectives in a decentralized network. The low SNR case is somewhat anomalous, since the SNR is close to the target SINR, so almost no interference can be tolerated. Similarly, to meet a high SINR constraint in a random network of reasonable density, the outage probability must be quite high, so this too may not be particularly meaningful.

To explain (4), recall that the Jensen-based approximation to outage probability is not accurate for dense networks and the plot shows that constant power (s=0s=0) is preferred at high densities.111Based on the figure it may appear that choosing s<0s<0, which means users with good channels transmit with additional power, outperforms constant power transmission. However, numerical results (not shown here) indicate that this provides a benefit only at extremely high densities for which outage probability is unreasonably large. Intuitively, a user with a poor channel in a dense network is extremely unlikely to be able to successfully communicate and global performance is improved by having such a user not even attempt to transmit, as done in the threshold-based policy studied in [16]. Fractional power control softens signal fading at the expense of more harmful interference power, and this turns out to be a good tradeoff in relatively sparse networks. In dense networks, however, there generally are a large number of nearby interferers and as a result the benefit of reducing the effect of signal fading (by increasing exponent ss) is overwhelmed by the cost of more harmful interference power. Note that this is consistent with results on channel inversion (s=1s=1) in [16], where s=0s=0 and s=1s=1 are seen to be essentially equivalent at low densities (as expected by the Jensen approximation) but inversion is inferior at high densities.

V Possible Areas for Future Study

Given the historically very high level of interest in the subject of power control for wireless systems, this new approach for power control opens many new questions. It appears that FPC has potential for many applications due to its inherent simplicity, requirement for only simple pairwise feedback, and possible a priori design of the FPC parameter ss. Some areas that we recommend for future study include the following.

How does FPC perform in cellular systems?. Cellular systems in this case are harder to analyze than ad hoc networks, because the base stations (receivers) are located on a regular grid and thus the tractability of the spatial Poisson model cannot be exploited. On the other hand, FPC may be even more helpful in centralized systems. Note that some numerical results for cellular systems are given in reference [17], but no analysis is provided.

Can FPC be optimized for spectral efficiency?. In this paper we have focused on outage relative to an SINR constraint as being the metric. Other metrics can be considered, for example maximizing the average spectral efficiency, i.e. max𝔼[log2(1+SINR)]\max\mathbb{E}[\log_{2}(1+{\rm SINR})], which could potentially result in optimal exponents s<0s<0, which is conceptually similar to waterfilling.

What is the effect of scheduling on FPC? If scheduling is used, then how should power levels between a transmitter and receiver be set? Will s=12s=\frac{1}{2} still be optimal? Will the gain be increased or reduced? We conjecture that the gain from FPC will be smaller but non-zero for most any sensible scheduling policy, as the effect of interference inversion is softened.

Can FPC be used to improve iterative power control? At each step of the Foschini-Miljanic algorithm (as well as most of its variants), transmitters adjust their power in a manner similar to channel-inversion, i.e., each transmitter fully compensates for the current SINR. While this works well when the target SINR’s are feasible, it does not necessarily work well when it is not possible to satisfy all users’ SINR requirements. In such a setting, it may be preferable to perform partial compensation for the current SINR level during each iteration. For example, if a link with a 10 dB target is currently experiencing an SINR of 0 dB, rather than increasing its transmit power by 10 dB to fully compensate for this gap (as in the Foschini-Miljanic algorithm), an FPC-motivated iterative policy might only boost power by 5 dB (e.g., adjust power in linear units according to the square root of the gap).

VI Conclusions

This paper has applied fractional power control as a general approach to pairwise power control in decentralized (e.g. ad hoc or spectrum sharing) networks. Using two approximations, we have shown that a fractional power control exponent of s=12s^{*}=\frac{1}{2} is optimal in terms of outage probability and transmission capacity, in contrast to constant transmit power (s=0s=0) or channel inversion (s=1s=1) in networks with a relatively low density of transmitters and low noise levels. This implies that there is an optimal balance between compensating for fades in the desired signal and amplifying interference. We saw that a gain on the order of 50%50\% or larger (relative to no power control or channel inversion) might be typical for fractional power control in a typical wireless channel.

We prove that for any non-negative random variable XX, the function

h(s)=𝔼[Xs]𝔼[Xs1],h(s)=\mathbb{E}\left[X^{-s}\right]\mathbb{E}\left[X^{s-1}\right], (25)

is convex in ss for ss\in\mathbb{R} with a unique minimum at s=12s=\frac{1}{2}. In order to show h(s)h(s) is convex, we show hh is log-convex and use the fact that a log-convex function is convex. We define

H(s)=logh(s)=log(𝔼[Xs]𝔼[Xs1]),H(s)=\log h(s)=\log\left(\mathbb{E}\left[X^{-s}\right]\mathbb{E}\left[X^{s-1}\right]\right), (26)

and recall Hölder’s inequality:

𝔼[XY](𝔼[Xp])1p(𝔼[Yq])1q,1p+1q=1.\mathbb{E}[XY]\leq\left(\mathbb{E}[X^{p}]\right)^{\frac{1}{p}}\left(\mathbb{E}[Y^{q}]\right)^{\frac{1}{q}},~~~~\frac{1}{p}+\frac{1}{q}=1. (27)

The function H(s)H(s) is convex if H(λs1+(1λ)s2)λH(s1)+(1λ)H(s2)H(\lambda s_{1}+(1-\lambda)s_{2})\leq\lambda H(s_{1})+(1-\lambda)H(s_{2}) for all s1,s2s_{1},s_{2} and all λ[0,1]\lambda\in[0,1]. Using Hölder’s with p=1λp=\frac{1}{\lambda} and q=11λq=\frac{1}{1-\lambda} we have:

H(λs1+(1λ)s2)\displaystyle H(\lambda s_{1}+(1-\lambda)s_{2}) =\displaystyle= log(𝔼[X(λs1+(1λ)s2)]𝔼[X(λs1+(1λ)s2)1])\displaystyle\log\left(\mathbb{E}\left[X^{-(\lambda s_{1}+(1-\lambda)s_{2})}\right]\mathbb{E}\left[X^{(\lambda s_{1}+(1-\lambda)s_{2})-1}\right]\right) (28)
=\displaystyle= log(𝔼[Xλs1X(1λ)s2]𝔼[Xλ(s11)X(1λ)(s21)])\displaystyle\log\left(\mathbb{E}\left[X^{-\lambda s_{1}}X^{(1-\lambda)s_{2}}\right]\mathbb{E}\left[X^{\lambda(s_{1}-1)}X^{(1-\lambda)(s_{2}-1)}\right]\right)
\displaystyle\leq log(𝔼[Xs1]λ𝔼[Xs2]1λ𝔼[Xs11]λ𝔼[Xs21]1λ)\displaystyle\log\left(\mathbb{E}\left[X^{-s_{1}}\right]^{\lambda}\mathbb{E}\left[X^{s_{2}}\right]^{1-\lambda}\mathbb{E}\left[X^{s_{1}-1}\right]^{\lambda}\mathbb{E}\left[X^{s_{2}-1}\right]^{1-\lambda}\right)
=\displaystyle= λlog(𝔼[Xs1]𝔼[Xs11])+(1λ)log(𝔼[Xs2]𝔼[Xs21])\displaystyle\lambda\log\left(\mathbb{E}\left[X^{-s_{1}}\right]\mathbb{E}\left[X^{s_{1}-1}\right]\right)+(1-\lambda)\log\left(\mathbb{E}\left[X^{s_{2}}\right]\mathbb{E}\left[X^{s_{2}-1}\right]\right)
=\displaystyle= λH(s1)+(1λ)H(s2).\displaystyle\lambda H(s_{1})+(1-\lambda)H(s_{2}).

This implies H(s)H(s) is convex, which further implies convexity of h(s)h(s). The derivative of hh is

h(s)=𝔼[Xs]𝔼[Xs1logX]𝔼[Xs1]𝔼[XslogX],h^{\prime}(s)=\mathbb{E}\left[X^{-s}\right]\mathbb{E}\left[X^{s-1}\log X\right]-\mathbb{E}\left[X^{s-1}\right]\mathbb{E}\left[X^{-s}\log X\right], (29)

and it can easily be seen that s=12s^{*}=\frac{1}{2} is the unique minimizer satisfying h(s)=0h^{\prime}(s)=0.

References

  • [1] N. Jindal, S. Weber, and J. G. Andrews, “Fractional power control for decentralized wireless networks,” in Proc., Allerton Conf. on Comm., Control, and Computing, Monticello, IL, Sept. 2007, available at http://www.ece.umn.edu/users/nihar/Publications.html.
  • [2] A. Goldsmith and P. Varaiya, “Capacity of fading channels with channel side information,” IEEE Trans. on Info. Theory, pp. 1986–92, Nov. 1997.
  • [3] G. J. Foschini and Z. Miljanic, “A simple distributed autonomous power control algorithm and its convergence,” IEEE Trans. on Veh. Technology, vol. 42, no. 8, pp. 641–646, Nov. 1993.
  • [4] R. D. Yates, “A framework for uplink power control in cellular radio systems,” IEEE Journal on Sel. Areas in Communications, vol. 13, no. 7, pp. 1341–47, Sept. 1995.
  • [5] N. Bambos, S. Chen, and G. J. Pottie, “Radio link admission algorithms for wireless networks with power control and active link quality protection,” in Proc., IEEE INFOCOM, Boston, MA, Apr. 1995, pp. 97–104.
  • [6] J. Herdtner and E. Chong, “Analysis of a class of distributed asynchronous power control algorithms for cellular wireless systems,” IEEE Journal on Sel. Areas in Communications, vol. 18, no. 3, Mar. 2000.
  • [7] J. F. Chamberland and V. V. Veeravalli, “Decentralized dynamic power control for cellular CDMA systems,” IEEE Trans. on Wireless Communications, vol. 2, no. 3, pp. 549–59, May 2003.
  • [8] M. Schubert and H. Boche, QoS-Based Resource Allocation and Transceiver Optimization. NOW: Foundations and Trends in Communications and Information Theory, 2005.
  • [9] M. Chiang, P. Hande, T. Lan, and C. W. Tan, Power Control in Cellular Networks. NOW: Foundations and Trends in Networking, To appear 2008.
  • [10] T. ElBatt and A. Ephremides, “Joint scheduling and power control for wireless ad hoc networks,” in Proc., IEEE INFOCOM, June 2002, pp. 976–84.
  • [11] R. Cruz and A. V. Santhanam, “Optimal routing, link scheduling and power control in multihop wireless networks,” in Proc., IEEE INFOCOM, Apr. 2003, pp. 702– 711.
  • [12] M. Haenggi, “The impact of power amplifier characteristics on routing in random wireless networks,” in Proc., IEEE Globecom, San Francisco, CA, Dec. 2003, pp. 513–17.
  • [13] S. Agarwal, S. V. Krishnamurthy, R. H. Katz, and S. K. Dao, “Distributed power control in ad-hoc wireless networks,” in Proc., IEEE PIMRC, Oct. 2001, pp. 59–66.
  • [14] V. Kawadia and P. R. Kumar, “Power control and clustering in ad hoc networks,” in Proc., IEEE INFOCOM, 2003.
  • [15] M. Chiang, “Balancing transport and physical layers in wireless multihop networks: Jointly optimal congestion control and power control,” IEEE Journal on Sel. Areas in Communications, vol. 23, no. 1, pp. 104–16, Jan. 2005.
  • [16] S. Weber, J. G. Andrews, and N. Jindal, “The effect of fading, channel inversion, and threshold scheduling on ad hoc networks,” IEEE Trans. on Info. Theory, vol. 53, no. 11, pp. 4127 – 4149, Nov. 2007.
  • [17] W. Xiao, R. Ratasuk, A. Ghosh, R. Love, Y. Sun, and R. Nory, “Uplink power control, interference coordination and resource allocation for 3GPP E-UTRA,” in Proc., IEEE Veh. Technology Conf., Sept. 2006, pp. 1–5.
  • [18] S. Weber, X. Yang, J. G. Andrews, and G. de Veciana, “Transmission capacity of wireless ad hoc networks with outage constraints,” IEEE Trans. on Info. Theory, vol. 51, no. 12, pp. 4091–4102, Dec. 2005.
  • [19] F. Baccelli, B. Blaszczyszyn, and P. Muhlethaler, “An Aloha protocol for multihop mobile wireless networks,” IEEE Trans. on Info. Theory, pp. 421–36, Feb. 2006.
  • [20] S. Weber, J. G. Andrews, X. Yang, and G. de Veciana, “Transmission capacity of wireless ad hoc networks with successive interference cancellation,” IEEE Trans. on Info. Theory, vol. 53, no. 8, pp. 2799–2814, Aug. 2007.
Refer to caption
Figure 1: The loss factor LL vs. ss for Rayleigh fading. Note that LcpL^{\mathrm{cp}} and LciL^{\mathrm{ci}} are the left edge and right edge of the plot, respectively.
Refer to caption
Figure 2: The outage probability (simulated, lower bound, and Jensen’s approximation) vs. FPC exponent ss for the default parameters.
Refer to caption
Refer to caption
Figure 3: The outage probability (simulated, lower bound, and Jensen’s approximation) vs. FPC exponent ss for α=2.2\alpha=2.2 (left) and α=5\alpha=5 (right).
Refer to caption
Refer to caption
Figure 4: The outage probability (simulated, lower bound, and Jensen’s approximation) vs. FPC exponent ss for SNR=10{\rm SNR}=10 dB (left) and SNR=30{\rm SNR}=30 dB (right).
Refer to caption
Refer to caption
Figure 5: The outage probability (simulated, lower bound, and Jensen’s approximation) vs. FPC exponent ss for β=10\beta=-10 dB (left) and β=10\beta=10 dB (right).
Refer to caption
Refer to caption
Figure 6: The outage probability (simulated, lower bound, and Jensen’s approximation) vs. FPC exponent ss for λ=0.00001\lambda=0.00001 (left) and λ=0.001\lambda=0.001 (right).
Refer to caption
Figure 7: The optimal choice of FPC exponent ss vs. PL exponent α\alpha, with ±1\pm 1% and ±10\pm 10% selections for ss.
Refer to caption
Figure 8: The optimal choice of FPC exponent ss vs. transmitter SNR =ρη=\frac{\rho}{\eta}, with ±1\pm 1% and ±10\pm 10% selections for ss.
Refer to caption
Figure 9: The optimal choice of FPC exponent ss vs. SINR constraint β\beta, with ±1\pm 1% and ±10\pm 10% selections for ss.
Refer to caption
Figure 10: The optimal choice of FPC exponent ss vs. density λ\lambda, with ±1\pm 1% and ±10\pm 10% selections for ss.