SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Green Radio.

Open Access Open Badges Research

Distributed beam forming with phase-only control for green cognitive radio networks

Xiaohua Lian*, Homayoun Nikookar and Leo P Ligthart

Author Affiliations

Microwave Technology and System for Radar (MTS-Radar), Department of Electrical Engineering, Mathematics and Computer Science, Delft University of Technology, Mekelweg 4, 2628CD Delft, The Netherlands

For all author emails, please log on.

EURASIP Journal on Wireless Communications and Networking 2012, 2012:65  doi:10.1186/1687-1499-2012-65

The electronic version of this article is the complete one and can be found online at: http://jwcn.eurasipjournals.com/content/2012/1/65

Received:1 October 2011
Accepted:27 February 2012
Published:27 February 2012

© 2012 Lian et al; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Cognitive radio (CR) is an intelligent radio system and is able to share the spectrum with licensed users (LU). By adopting adaptive beam forming techniques, CR can reuse the spectrum with LU via directing main beams towards CR users while displaying nulls towards LU. In this article, we present a new distributed beam forming (DB) technique, and study the performance of its application to decentralized CR networks. The presented DB method only controls the phase of the transmitted signal of each CR node in the CR network, and therefore, it is called phase-only DB (PODB) method. It can be implemented at each node independently with only prior knowledge of its own location, and directions of distant CR (DCR) users and LU. The average beam pattern and the average gains for an arbitrary number of CR nodes of PODB method are discussed, showing that CR nodes can constructively transmit signals to DCR users with less interferences to LU via employing the PODB method to form and direct main beams towards directions of DCR users while nulls towards directions of LU. PODB proves to be a green DB method by prolonging the lifetime of the CR network due to effective battery power consumptions at CR nodes. The cumulative distributed function of the beam pattern for a large number of CR nodes is derived and analysed to demonstrate that the PODB method increases the possibility that the transmitted power of the whole CR network at the directions of LU is lower than a certain threshold, which guarantees that CR network causes less disturbing effect to LU.

distributed beam forming; collaborative beam forming; cognitive radio networks; phase-only distributed beam forming; null-steering distributed beam forming; green cognitive radio networks

1. Introduction

Cognitive radio (CR) is a new concept radio presented in [1,2], which is a promising approach to solve the intensive usage of the precious natural source-spectrum [3]. As described in [3], CR is capable of sensing the communication environment and adapting to it by adjusting its parameters. By applying a beam forming technique, CR can direct its main beams towards CR users while putting nulls to licensed users (LU) in uplink in order to share the same spectrum with LU without disturbing them [4,5].

Distributed beam forming (DB), which is also referred to as collaborative beam forming, is originally employed as an energy-efficient scheme to solve the long distance transmission by wireless sensors networks (WSN), in order to reduce the amount of required energy and consequently to extend the utilization time of the sensors. The basic idea of DB is that a set of nodes in wireless networks act as a virtual antenna array and then form a beam towards a certain direction to collaboratively transmit a signal. DB is proposed in [6] and it shows that by employing K collaborative nodes, the collaborative beam forming can result in up to K-fold gain in the received power at a distance access point. A cross-layer approach to DB for wireless adhoc networks is discussed in [7] with more complicated models and two time phases of communication steps. The improved beam pattern and connectivity properties are shown in [8], and a reasonable beam forming performance affected by a node synchronization error is discussed in [9]. DB has also been introduced in relay communication systems [10-14]. In [10-13], several DB techniques for relay networks with flat fading channels have been developed, while in [14] frequency selective fading has been considered. DB requires accurate synchronization; in other words, the nodes must start transmitting at the same time, synchronize their carrier frequencies, and control their carrier phases so that their signals can be combined constructively at the destination. Several synchronization techniques for CB can be found in [15,16], and a review is given in [17].

As discussed in [4,5], the CR base station (BS) needs to be equipped with array antenna (AA) so that it can form and direct main beams to its own users while steering nulls towards the directions of LU via adaptive beam forming techniques. In this study, instead of requesting CR BS to be deployed together with antenna arrays, we consider decentralized CR networks, which is formed by distributed CR users in a certain area. Those CR users are therefore regarded as CR nodes. Meanwhile DB is suggested to be employed by the CR network to form beams towards distant CR (DCR) users so that the CR network is able to forward the signals to DCR users cooperatively. Thus, the whole network decreases the energy consumption effectively, and also increases the communication range.

However, the CR network distinguishes itself from WSN in that it is designed to achieve two goals under one restricted condition, i.e., the received power at DCR users should be maximized, while that at LU be minimized. The restriction limits that CR nodes are only aware of its own location, as well as the directions of DCR users and LU. The constraint supports DB for a CR network functional even without cluster nodes. If cluster nodes are selected, each CR node is able to have the information of locations of all CR nodes during the process of information exchanging and sharing. Then, the DB problem can be described as designing adaptive beam formers with an irregular AA structure. In this article, we consider no CR cluster nodes, which is the basic and simplest network application. For DB, the aim of a CR network can directly be translated into generating main beams towards DCR users while displaying nulls towards LU. In addition, the weights calculation can be performed by each node, and for each node it should be independent from the prior knowledge of other nodes due to the restriction introduced above. This article also suggests generating a certain angular range of nulls instead of point nulls towards LU. The reason is due to the study of the probability density function (pdf) of the angle of arrival (AOA) of the scattered wave at the mobile station in wireless communications [18-20], which will be explained later in this article.

Recently, DB for a CR network has also been discussed in [21-23]. In [21], the effect of noise in phase synchronization on the resultant beam pattern has been analysed. A DB scheme based on zero-forcing beam forming has been presented in [22], where the idea that a CR network can be regarded as a relay network to forward the signal to the destination CR node is adopted. It is worth noticing that in [23], the authors have also addressed the above problem of directing main beam to a required direction while putting nulls to unwanted directions. This fits well the idea of introducing DB to CR network. A novel null-steering DB method has been presented in [23] in a way that each node can calculate its own weight by itself without prior knowledge of others.

Despite that the DB weights derived in [23] achieve our goal of CR network successfully, further improvements of DB in application to green CR networks are appreciated. In this article, we present a phase-only DB (PODB) method by only controlling the phase of the transmitted signal of each node. With the proposed PODB method, the lifetime of the whole CR network is prolonged, which will be proved in this article.

The PODB method is able to direct main beams towards DCR users while nulls are formed towards directions of LU. It can be implemented independently by each node based on very limited knowledge, including directions of DCR users and LU, and its own location. The weight of the PODB method can be regarded as nothing more than a phase perturbation of each node, which has the same amplitude as that of others. By adding certain virtual directions, the proposed PODB method can also generate broadened nulls in a defined range around directions of LU. The average beam pattern is also derived to show what the beam pattern with an arbitrary number of collaborative nodes will converge to. The average gains at both directions of LU and DCR users are calculated approximately. What is more, the cumulative distributed function (CDF) of the beam pattern is derived; it demonstrates the possibility that the transmitted power of the whole CR network at every direction can be lower than a certain threshold. The lifetime of a wireless network with both null-steering DB which has been presented in [23] and the PODB method is discussed and compared, in order to show the "green" aspect of the latter method.

The rest of the article is organized as follows. First in Section 2, the model of the CR network will be presented, as well as several necessary assumptions. Then the PODB method is proposed with and without broadened nulls. The average beam pattern of the PODB method, and its gains at certain directions are discussed in Section 3, as well as the derived CDF of the beam pattern generated by the PODB method. Several simulation results are shown in Section 4. The greenness of the PODB method is illustrated in Section 5 by comparing the lifetime of the wireless networks adopting null-steering DB and the PODB method. Finally, Section 6 concludes the article.

Notations: ||a||, ||a||2 denote the infinitive norm and Euclidean norm of vector a. (·)T, (·)* and (·)H denote the transpose, the conjugate, and the Hermitian transpose, respectively. (A)m, n represents the element of the mth row and nth column of matrix A, and (a)k represents the kth element of vector a. E[·] stands for the statistical expectation, and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M1">View MathML</a> denotes converge with probability one. Jn(·) is the nth-order Bessel function of the first kind, and I0(·) is the zeroth-order-modified Bessel function. Re[·] and Im[·] represent the real and imaginary parts of a variable, respectively.

2. Problem formulations and the proposed PODB method

2.1. Geometrical structure of CR networks

The geometrical structure of the CR network and also distant receiver terminals including LU and DCR users are illustrated in Figure 1. As shown in Figure 1, K CR nodes are uniformly distributed on a disc centred at O with radius R. Let us denote their polar coordinates of the kth CR node as (rk, Ψk). L DCR users are considered as access points, which locate in the same plane at <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M2">View MathML</a>. Meanwhile, M LU are also coexisting with DCR users. Their locations are (Ai, ϕi), i = 1,2,...,M. The CR nodes in the CR network are requested to form a virtual antenna array and collaboratively transmit a common message S(t).

thumbnailFigure 1. CR network with DCR users and LU.

2.2. Necessary assumptions

Without generality, we also adopt the following assumptions:

(1) The number of CR nodes are larger than that of DCR users together with LU, i.e. K > L+M. This is basically required to solve the later matrix equations. In addition, many adaptive beam forming techniques also request that the number of AAs is larger than that of constrained directions where main beams and nulls are dedicatedly directed towards.

(2) All DCR users and LU are located in the far-field of CR network, such that

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M3">View MathML</a>

(3) The bandwidth of S(t) is narrow, so that it guarantees that S(t) is almost constant during 2R/c second where c is the speed of light. It has been discussed in [24] that OFDM scheme is a proper and recommended candidate for CR due to its flexible adaptation of the spectrum. Since OFDM signal can be regarded as a combination of narrow-band modulated signals, the proposed method in this article can also be implemented in a wide band CR OFDM system.

The spatial effect of signal scattering and reflection at LU are taken into account. As a result, the angular spread at LU is suggested to be eliminated by displaying spread nulls around each direction of LU rather than point nulls. However, other channel effects, such as multipath fading and shadowing, are ignored in this article.

For simplicity, the proposed PODB is explained and presented when L = 1. In case L > 1, we show in [25] how to generate more than one main beam towards DCR users. We simplify <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M4">View MathML</a> by (A0, ϕ0) when there is only one DCR user.

2.3. System model of CR network

Let xk(t) denote the transmitted signal from kth node,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M5">View MathML</a>


where wk is the weight adopted by the kth node. The received signal at an arbitrary point (A, ϕ) in the far-field due to the kth node transmission is [23],

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M6">View MathML</a>


where dk is the distance between the kth node and the access point (A, ϕ), and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M7">View MathML</a> is the signal path loss with γ donating the path loss exponent. In [23], βk and dk are approximated by

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M8">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M9">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M10">View MathML</a>. The two approximations above are due to assumption 2. It also ensures <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M11">View MathML</a>, and then βk β.

Substituting Equations (3) and (4) into Equation (2) it follows that [23],

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M12">View MathML</a>


If we adopt the initial phase of each node as <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M13">View MathML</a>[6], the received signal at (A, ϕ) is

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M14">View MathML</a>


We define

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M15">View MathML</a>


The power received at (A, ϕ) is denoted by p(ϕ). If S(t) has normalized power, p(ϕ) is shown as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M16">View MathML</a>


The goal of our proposed PODB method is to find w, where w = (w1, w2, ..., wK)T, which satisfies

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M17">View MathML</a>


Based on Equation (7), it is easy to verify that <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M18">View MathML</a>. Therefore, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M19">View MathML</a> is always satisfying.

The constraint in (9) shows that each wk is no more than a phase adjustment. In other words, we only control the phase modification of each node.

A different optimization problem is discussed in [23], which has been proposed based on null-steering beamformers. It has the constraint of the total power of w, i.e. wHw PT, instead of |wi| = 1, i = 1,2,...,K.

2.4. The proposed PODB method

We propose the solution of Equation (9) in this section. First we assume W = 0, where W is the width of the spread nulls towards each LU, which requires only a point null at each LU direction, and then we give further results of how to generate spread nulls by modifying the angles.

Equation (9) also appeared in the application of adaptive array techniques to solve the problems of array pattern nulling, which was studied in [26]. The proposed solution in [26] linearized the nonlinear problem by considering small position perturbation of each antenna element, and then derived analytic results. It was also proved in [26] that good nulling performance could be achieved, even though the solutions were based on approximations. Unfortunately, the results shown in [26] could be only adopted by uniformly linear arrays (ULA).

In the case of CR network, the structure and the motivation of the beam forming problem are very different from that in [26]. As mentioned before, our goal is to find a set of weights (w1, w2, ..., wK), which can be implemented by each node under the constraint that only the information of (rk, Ψk), the directions of DCR (ϕ0) and LU (ϕi, i = 1,2,..., M) are available.

We assume that wk has the form of

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M20">View MathML</a>


Then, we perform a Taylor expansion in Equation (7) and retain the first two terms. We obtain

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M21">View MathML</a>


FC(ϕ) represents the second item of (11) and is a cancellation pattern which can be used to achieve M nulls.

Let we define the following variables

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M22">View MathML</a>

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M23">View MathML</a>

where μ ∈ ℜK × 1, x= Re[u], and y= Im[u].

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M24">View MathML</a>

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M25">View MathML</a>

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M26">View MathML</a>

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M27">View MathML</a>

where e is a K × 1 vector with one as its entities.

Thus, Equation (11) can be simplified as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M28">View MathML</a>


According to Equation (9), when W = 0, in order to steer nulls at direction of ϕm, = 1,2,..., M, we require

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M29">View MathML</a>


Due to assumption 1, Equation (13) has a solution and can be solved as the following. Using the results in [6],

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M30">View MathML</a>


Equation (13) can be further written as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M31">View MathML</a>


If we define

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M32">View MathML</a>

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M33">View MathML</a>

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M34">View MathML</a>

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M35">View MathML</a>

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M36">View MathML</a>

where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M37">View MathML</a>.

The solution of Equation (15) is

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M38">View MathML</a>


where (0) is a vector with length M and 0 as its entries.

To make the calculation of μ be able to be implemented by each node in a distribution fashion, we need to calculate Γ differently, which has a good approximation of original counterparts, and the entries of it depends only on the information that each node is entitled to have (rkk), ϕ0, ϕm, m = 1,2,...,M.

Theorem 1: Give a matrix U, which is defined as

U ≜ (u1, u2, ..., uM) = X+jY, we have the following results:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M39">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M40">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M41">View MathML</a>.

Proof: See Appendix 1

Theorem 2: The entries of Γ have the following results:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M42">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M43">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M44">View MathML</a>


Proof: See Appendix 2

Based on Theorem 2, we can replace Γ by E[Γ], and then calculate (E[Γ])-1 by

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M45">View MathML</a>


If we define

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M46">View MathML</a>

Then by combining (16) and (22), we have

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M47">View MathML</a>


The weight of the kth CR node is represented by

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M48">View MathML</a>


So far we have derived the weight of the PODB method, which shows that the computation of wk for each node only depends on c, Γ2, and (ym)k. The first two can be regarded as constant vectors and a constant matrix, as shown in their definitions and Theorem 2. While in order to calculate (ym)k, each node only requires the prior knowledge of its own location (rk, Ψk), and directions of all LU ϕm, m = 1,2,...,M. It is worth noticing that the information of the direction of DCR user ϕ0 is also needed for deciding the initial phase of each node φk in Equation (6).

2.5. PODB method with spread nulls

The principle of the necessity of generating spread nulls has been explained in [25]. In wireless communication systems, the received signal and power spectra at the mobile station at wireless communication depend on the pdf of the AOA of the scattered wave. Clarke considered a uniform AOA pdf over [-π, π) [18]. However, it has been argued and experimentally demonstrated that the scattering encountered in many environments results into a non-uniform pdf of AOA [19]. Reference [20] suggests the two-parameter Von Mises pdf as a flexible and generalized model for the pdf of AOA, which includes non-isotropic scattering cases, and also the isotropic one as a special case. The pdf is given as [20]

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M49">View MathML</a>


where θp accounts for the main direction of the AOA scatter components. Parameter l ≥ 0 controls the width of the AOA of scatter components. Figure 2 shows the pΘ(θ) with different l for θp = 0. Based on Von Mises model, point nulls which can be generated by CR networks are not sufficient enough for energy depression at directions of LU, because the power around the nulling direction may also leak into LU due to spatial scattering. As a result, instead of generating point nulls, CR network is required to display spread nulls around each direction of the LU.

thumbnailFigure 2. Von Mises pdf for the AOA of scatter components at mobile station.

Null broadening (NB) technique was originally developed as a robust array beam forming technique. It was also regarded as a beam pattern synthesis method. Mailloux [27] presented a NB method by simple modification of the covariance matrix of the received signal. While in [28], a similar NB algorithm was also presented by applying a transform to the same covariance matrix. Both methods were capable of providing notch at locations of interferences signals. In [29], two NB algorithms as covariance matrix taper methods are defined, and considered as effective robust adaptive beam forming techniques, imparting robustness into adaptive pattern by judicious choice of null placement and width. As for wireless communication, NB was employed in a cellular communication system, particularly an example was given in space division multiple access system for downlink beam forming [30].

The above proposed NB methods are all based on ULA. They all require the prior knowledge of the distance between every two AA elements. In our case, this requested information can be translated into the distance between every two CR nodes in CR networks. Consequently, this does not align with our constraint that each node is only aware of its own location.

However, one simple way of generating spread nulls is to add N virtual LU sources ϕi, i = 1,2,..., N around each direction of existing LU with the angle range of W/2, e.g. |ϕi-ϕm| ≤ W/2, m = 1,2,...,M. The number of those virtual directions N can be chosen with a compromising between computation load and the depth of spread nulls. If N is large enough, all the angles around directions of the LU are continuously all regarded as nulling points, thus the spread nulls in the beam pattern are maintaining the same depth with that of the nulling point ϕm. However, with large N, the size of matrix Γ2 in Theorem 2 increases to (MN) × (MN), and consequently each node must have larger computation burden.

3. Properties of the beam pattern generated by the PODB method

In this section, we discuss the properties of the beam pattern that generated by the presented PODB method. The analysed properties include the average beam pattern for an arbitrary number of CR nodes K, the average gains at the directions of DCR user ϕ0, and at those of all the LU ϕm, m = 1,2,...,M. Meanwhile the CDF of the beam pattern with larger K is also derived.

3.1. Average beam pattern

The result of average beam pattern is given by Theorem 3.

Theorem 3: Using the presented PODB method, the average beam pattern, which is defined as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M50">View MathML</a>

when ϕϕ0, can be approximated by

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M51">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M52">View MathML</a>

and the approximation error <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M53">View MathML</a> is bounded by <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M54">View MathML</a>.

Proof: See Appendix 3

3.2. Average gain at direction of DCR

The average gain at ϕ0 is approximated by 1, i.e.

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M55">View MathML</a>


Proof: See Appendix 4

3.3. Average gains at directions of LU

The average gain at ϕm, m = 1,2,...,M is

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M56">View MathML</a>


Proof: See Appendix 5

For CR network, it is very important to limit the transmission power at directions of LU. From Equation (27) and (28), we can see that towards LU directions, the average transmission power are K - 1 times lower than that towards the direction of DCR user.

When there is no LU available, i.e. M = 0. Equation (23) gives the result of μk = 0. Consequently,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M57">View MathML</a>


This is consistent with the average beam pattern of the DB method presented in [6], which has no specific nulling points.

3.4. CDF of the beam pattern generated by PODB method

The CDF of the beam pattern generated by PODB method with large number of K is given by Theorem 4.

Theorem 4: When K is large enough, the CDF of the beam pattern at different ϕcan be approximated by

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M58">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M59">View MathML</a> for γ ≥ 0 is the CDF of the Rayleigh distribution, sgn(x) = |x|/x and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M60">View MathML</a>.

Proof: See Appendix 6.

We can write when ϕ= ϕm, m = 1,2,..., M

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M61">View MathML</a>


Proof: See Appendix 7.

From the above equation, we can conclude that when ϕ= ϕm, m = 1,2,..., M, the probability of the total transmitted power of CR network lower than P0 can be expressed by a constant, which is decided by P0, K and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M62">View MathML</a>. Based on Equation (21) and the definition of c, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M62">View MathML</a> can be only determined by ϕ0, ϕm, m = 1,2,...,M and R/λ. Usually P0 is defined by LU due to its spectrum interference tolerance threshold. Consequently, if we want to increase the probability of the case that CR network transmitting less power than P0 towards directions of the LU, we need to increase K. In other words, more CR nodes need to be added to the same CR network (maintaining the same network structure R/λ) to collaboratively perform PODB.

4. Simulations and results

The average beam pattern of DB methods are shown in Figure 3. They are DB method without considering nulling directions, which has been presented in [6], the proposed PODB method, and null-steering DB method in [23]. In the simulation, we choose R/λ = 2, K = 100 and ϕ0 = 0°. Meanwhile, two LU are considered, which locate at the directions of ϕ1 = 30°, ϕ2 = -15°. Figure 3 shows that PODB and null-steering DB methods successfully depress transmission power towards the directions of LU. Their performances are similar to each other. The DB method without nulling has a symmetrical average beam pattern without particular nulling points. As mentioned in Equations (27) and (28), the average gain at ϕ0 = 0° approximates to 1, while those at ϕ1 = 30°, ϕ2 = -15° are close to 10 log(1/K) = -20 dB.

thumbnailFigure 3. Average Beampatterns of three DB methods.

Figure 4 shows the average beam pattern with only one LU, which has the direction of ϕ = 30°. The other simulation conditions remain the same with Figure 3. In Figure 4, we also consider the average beam pattern of PODB method with broadened nulls, which assumes two virtual LU directions around ϕ = 30°, which are located at ϕ1 = 20°, ϕ2 = 40°. It can be seen from Figure 4 that PODB method with NB generates spread nulls around ϕ = 30° with the width of W = 20° and depth about -20 dB.

thumbnailFigure 4. Average Beampatterns of PODB methods.

Figure 5 shows the CDF of the beam pattern at different angle ϕ, which has been discussed in Theorem 4. As required in Theorem 4, the number of nodes in CR network needs to be large, so that the sum of random variables can be regarded as Gaussian distribution random variables. Therefore, we consider two cases with K = 100, and K = 500. Meanwhile, only one LU is considered, i.e. ϕ = 5°. The DCR user is at the direction of ϕ0 = 0°. P0 in Theorem 4 is chosen as P0 = 0.01. It can be seen from Figure 5 that, in the angle range of 5° ≤ ϕ ≤ 7°, considering the probability that the transmitted power of CR network lower than P0, PODB method has much higher chance than the DB method without nulling. Meanwhile, as discussed in Section 3.4, Prob{P(ϕm)|P(ϕm) ≤ P0} is a constant which depends on ϕ1, ϕ0 and K. When more nodes participate in CR network to perform DB (K is increasing), Prob{P(ϕm)|P(ϕm) ≤ P0} is also increasing, which is also demonstrated in Figure 5.

thumbnailFigure 5. Probability of with different K.

5. Greenness of the PODB method

The battery consumptions of the network adopting the proposed PODB method and null-steering DB method presented in [23] are illustrated in Figures 6 and 7. Figure 6 shows that to perform PODB, each node of CR network consumes the same amount of battery power at the same time. The power consumption of the whole network is "balanced", because of the constraint in Equation (9), i.e. |wk| = 1. While null-steering DB method utilizes the battery power of each node "unevenly", as shown in Figure 7. When the battery of one node is dying, the others nodes are still power supplied.

thumbnailFigure 6. Power consumption of CR network with PODB method.

thumbnailFigure 7. Power consumption of CR network with null-steering DB method.

For null-steering method, we introduce a new concept of "DB cycle". It describes the process from the moment that all nodes in the network start participating DB till the moment that one node in the network first runs out of its battery. Then, the rest of the nodes of the network will start computing the new weights of DB method with one less number of the nodes, and thus another round of DB cycle begins.

The number of CR nodes K in a DB cycle is requested to satisfy k K0, 1 ≤ K0 K, where K0 is the least number of CR nodes required by CR network to perform DB, due to the restriction of the transmission power towards directions of LU

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M63">View MathML</a>


We denote the normalized the weights of the ith DB cycle of null-steering DB method as <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M64">View MathML</a> and it satisfies that <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M65">View MathML</a>. Assume <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M66">View MathML</a> is the length of the time period of the ith DB cycle of the null-steering DB method

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M67">View MathML</a>


where p(i) is the battery power consumption of the CR node, which finishes its battery first in the ith DB cycle. Then, the total lifetime of CR network by adopting null-steering method can be written as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M68">View MathML</a>


The total energy consumption is

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M69">View MathML</a>


The total amount of energy of the whole network is fixed, due to the fact that all nodes will consume all of their batteries power in the end. Let us denote the total energy of the CR network as E.

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M70">View MathML</a>


It is easy to verify that

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M71">View MathML</a>


Only when K0 = 1,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M72">View MathML</a>


However, the concept of DB cycle is not necessary for PODB, because all the nodes will be finishing their batteries at the same time. We define the normalized weights of the PODB method as wpo, and it satisfies <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M73">View MathML</a>. Then, the lasting time of the network with the PODB method using the same amount of energy E is

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M74">View MathML</a>


Using the results of Equations (34), (35), (37) and (39), we can write that

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M75">View MathML</a>


Thus, it has been proved that the life time of the network adopting PODB method is longer than that adopting null-steering method. Only when K0 = 1, the network has the same lifetime by employing both two DB methods.

6. Conclusions

A new DB method, which is called PODB method, has been proposed in this article in the application of a green CR network. It is implemented in a way that CR nodes transmit signals constructively towards directions of DCR users and destructively to the LU. Considering the average beam pattern of the PODB method, it achieves the highest gain at the direction of DCR users and the lowest to the LU. The PODB method can be computed at each node in a so called distributed way, i.e., each node is only aware of its own location, and the directions of LU and DCR users. The PODB is different from the null-steering-based DB method meaning that the weights of the former can be regarded as a phase perturbation of each node. Consequently, it utilizes energy of the whole network evenly and effectively, a characteristic which renders itself to green concept.

In a CR network, it is extremely important to calculate the probability of the case that the power transmitted by the CR network towards directions of LU is lower than a certain threshold. Therefore, the CDF of the beam pattern of the PODB method has been derived in this article assuming a large number of CR nodes. Comparing with the DB method without particular nulling directions, the PODB method achieves much higher chances of transmitting less power towards directions of LU. To increase the chances that the disturbances to LU, which are caused by the CR network are less than the interference power threshold that the LU can tolerate, the presented CDF suggests that more CR nodes should participate in the CR network to perform DB.

Appendix 1

Proof of Theorem 1

It has been discussed in [23] that

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M76">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M77">View MathML</a>, and the pdf of zk has been presented in [6],

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M78','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M78">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M79">View MathML</a>


The above exponential item can be further simplified by the following

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M80">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M81">View MathML</a>

We define

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M82','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M82">View MathML</a>

Equation (43) can be written as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M83">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M84">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M85">View MathML</a> has the same distribution with zk, because the pdf of zk in (42) has nothing to do with ϕm, ϕn. Therefore, using the same method, (45) can be further written as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M86">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M87">View MathML</a>


Especially, when m = n,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M88','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M88">View MathML</a>


Equations (41), (47) and (49) show the results of Theorem 1.

Appendix 2

Proof of Theorem 2

Based on the definition of Γ,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M89','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M89">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M90','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M90">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M91','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M91">View MathML</a>


Using the results in Theorem 1,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M92','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M92">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M93','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M93">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M94','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M94">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M95','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M95">View MathML</a>


Equations (53), (54), (55) and (56) demonstrate the results of Theorem 2.

Appendix 3

Proof of Theorem 3

Based on Equations (7) and (10)

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M96">View MathML</a>


As mentioned before, μk is a small perturbation of the phase of each node. We here adopt the approximation of ex by ex ≈ 1+x , and it holds the following result [31],

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M97','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M97">View MathML</a>


Then we have

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M98','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M98">View MathML</a>


Equation (57) can be approximated by

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M99','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M99">View MathML</a>


We define a new random variable

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M100','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M100">View MathML</a>


Based on the pdf of zk, we have the following results [6]

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M101','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M101">View MathML</a>


Using the result in Theorem 1, it is easy to calculate the expectations of the following variables,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M102','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M102">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M103','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M103">View MathML</a>


Consequently from Equations (63), (64) and (23), the expectation of ξk can be calculated by

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M104','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M104">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M105','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M105">View MathML</a>


From Equation (65), it is easy to verify that

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M106','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M106">View MathML</a>


Meanwhile, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M107','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M107">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M108','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M108">View MathML</a> are identical independently distributed variables. ξk and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M108','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M108">View MathML</a> are independently distributed random variables. Thus, with the result in Equations (62) and (65), considering the general case that ϕϕ0, we can deduce the second item in Equation (60) further by

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M109','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M109">View MathML</a>


Substituting Equation (68) into Equation (60), we have final results

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M110','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M110">View MathML</a>


The error of the approximation in (61) can be bounded using the result of (58) by

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M111','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M111">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M112','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M112">View MathML</a>


Based on Equation (23) and the results of Equation (62), we can deduce that

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M113','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M113">View MathML</a>


In Theorem 2, it holds that

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M114','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M114">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M115','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M115">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M116','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M116">View MathML</a>


By substituting Equations (75), (72) into (70), we can write

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M117','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M117">View MathML</a>


The approximation error in the second item of Equation (60) can be written

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M118','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M118">View MathML</a>


Then the estimation error ε between <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M119','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M119">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M120','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M120">View MathML</a> can be written based on the result of Equation (77)

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M121','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M121">View MathML</a>


Appendix 4

Proof of Equation (27)

When ϕ = ϕ0 , the average beam pattern at ϕ0 can be written as:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M122','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M122">View MathML</a>


Using the same approximation method in Equation (59)

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M123','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M123">View MathML</a>


Based on the result in Equation (72), we can further write

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M124','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M124">View MathML</a>


Substituting Equation (81) into (80), we have the result of Equation (27)

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M125','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M125">View MathML</a>


Appendix 5

Proof of Equation (28)

Theorem 3 and the definition of Γ2 in Theorem 2 illustrate that, when ϕ= ϕm, m = 1,2,...,M,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M126','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M126">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M127','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M127">View MathML</a>


where vm denotes a M × 1 vector with ones as the entities of the mth row.

Meanwhile based on the definition of (q1, q2, ..., qM) we can write when ϕ= ϕm, m = 1,2,...,M,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M128','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M128">View MathML</a>


By substituting (85) into the average beam pattern in Theorem 3,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M129','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M129">View MathML</a>


Which shows the result of Equation (28).

Appendix 6

Proof of Theorem 4

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M130','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M130">View MathML</a>


We define a new random variable

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M131','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M131">View MathML</a>


Definition of ξk is shown in Appendix 3 Equation (61).

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M132','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M132">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M133','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M133">View MathML</a> are identical independently distributed variables. Thus, if K is large enough, according to the Central Limit Theorem, χ has a Gaussian distribution, i.e. χ ~ N(ρ, σ2), where

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M134','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M134">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M135','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M135">View MathML</a>


Using the result discussed in Equations (62) and (65) in Appendix 3

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M136','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M136">View MathML</a>


From Equation (75) in Appendix 3,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M137','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M137">View MathML</a>


Thus, by adopting the results shown in (91) and (92),

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M138','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M138">View MathML</a>


From Equation (93), if we define a new random variable γ by γ ≜ |χ-ρ|, the distribution of γ is Rayleigh distribution, i.e. <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M139','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M139">View MathML</a>, where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M140','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M140">View MathML</a>. We then are able to calculate the probability Prob{P(ϕ)|P(ϕ) ≤ P0}. We next discuss the outcome for the different cases.

If <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M141','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M141">View MathML</a>,


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M142','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M142">View MathML</a>


If <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M143','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M143">View MathML</a>,


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M144','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M144">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M145','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M145">View MathML</a>


By combining (94) and (95), it is easy to verify Theorem 4.

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M146','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M146">View MathML</a>


Appendix 7

Proof of Equation (31)

Using the result of Equation (59) in appendix E, we can write

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M147','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M147">View MathML</a>


As a result,

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M148','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/65/mathml/M148">View MathML</a>


Competing interests

The authors declare that they have no competing interests.


  1. J Mitola, Cognitive radio: making software radios more personal. IEEE Person Commun 6(4), 13–18 (1999). Publisher Full Text OpenURL

  2. J Mitola, Dissertation (Royal Inst. Technol. (KTH), 2000)

  3. S Haykin, Cognitive radio: brain-empowered wireless communications. IEEE J Sel Area Commun 23(2), 201–210 (2005)

  4. X Lian, H Nikookar, J Zhou, Adaptive robust beamformers for cognitive radio. Proc IEEE Wireless Technology (EuWiT) (Amsterdam, The Netherlands), pp. 103–106

  5. X Lian, H Nikookar, LP Ligthart, Adaptive OFDM beamformer with constrained weights for cognitive radio. Proc IEEE 69th Vehicular Technology Conference (VTC Spring) (Bacelona, Spain, 2009), pp. 1–5 (doi:10, 2009), . 1109/VETECS.2009.5073362 OpenURL

  6. H Ochiai, P Mitran, HV Poor, V Tarokh, Collaborative beamforming for distributed wireless ad hoc sensor networks. IEEE Trans Signal Process 53(11), 4110–4124 (2005)

  7. L Dong, AP Petropulu, HV Poor, A cross-layer approach to collaborative beamforming for wireless ad hoc networks. IEEE Trans Signal Process 56(7), 2981–2993 (2008)

  8. K Zarifi, S Affes, A Ghrayeb, Distributed beamforming for wireless sensor networks with random node location. Proc Int Conf Acoustics, Speech, and Signal Process (ICASSP'09) (Taipei, Taiwan, 2009), pp. 2261–2264 (doi:10, 2009), . 1109/ICASSP.2009.4960070 OpenURL

  9. R Mudumbai, G Barriac, U Madhow, On the feasibility of distributed beamforming in wireless sensor networks. IEEE Trans Wirel Commun 6(5), 1754–1763 (2007)

  10. Y Jing, H Jafarkhani, Network beamforming using relays with perfect channel information. Proc Int Conf Acoustics, Speech, and Signal Process (ICASSP'07) (Hawaii, USA, 2007), pp. 473–476 (doi:10, 2007), . 1109/ICASSP.2007.366575 OpenURL

  11. G Zheng, KK Wong, A Paulraj, B Ottersten, Collaborativerelay beamforming with perfect CSI: optimum and distributed implementations. IEEE Signal Process Lett 16(4), 257–260 (2009)

  12. V Havary-Nassab, S Shahbazpanahi, A Grami, ZQ Luo, Distributed beamforming for relay networks based on second-order statistics of the channel state information. IEEE Trans Signal Process 56(9), 4306–4316 (2008)

  13. S Fazeli-Dehkordy, S Shahbazpanahi, S Gazor, Multiple peer-to-peer communications using a network of relays. IEEE Trans Signal Process 57(8), 3053–3062 (2009)

  14. H Chen, AB Gershman, S Shahbazpanahi, Filter-and-forward distributed beamforming in relay networks with frequency selective fading. IEEE Trans Signal Process 58(3), 1251–1262 (2010)

  15. DR Brown III., HV Poor, Time-slotted round-trip carrier synchronization for distributed beamforming. IEEE Trans Signal Process 56(11), 5630–5643 (2008)

  16. JA Bucklew, WA Sethares, Convergence of a class of decentralized beamforming algorithms. IEEE Trans Signal Process 56(6), 2280–2288 (2008)

  17. R Mudumbai, DR Brown, U Madhow, HV Poor, Distributed transmit beamforming: challenges and recent progress. IEEE Commun Mag 47(2), 102–110 (2009)

  18. GL Stuber, Principles of Mobile Communication (Kluwer, Boston, MA, 1996)

  19. JS Sadowsky, V Kafedziski, On the correlation and scattering functions of the WSSUS channel for mobile communication. IEEE Trans Veh Technol 47(1), 270–282 (2008)

  20. A Abdi, JA Barger, M Kaveh, A parametric model for the distribution of the angle of arrival and the associated correlation function and power spectrum at the mobile station. IEEE Trans Veh Technol 51(3), 425–434 (2002). Publisher Full Text OpenURL

  21. Y Deng, AG Buur, D Pearce, D Grace, Distributed beamforming for cognitive radio networks. Proc Communications and Networking in China (CHINACOM'08) (Hangzhou, China, 2008), pp. 1206–1210 (doi:10, 2008), . 1109/CHINACOM.2008.4685243 OpenURL

  22. J Liu, W Chen, Z Cao, Y Zhang, A distributed beamforming approach for enhanced opportunistic spectrum access in cognitive radios. Proc Global Telecommunications Conference (Globecom) (Hawaii, USA, 2009), pp. 1–6 (doi:10, 2009), . 1109/GLOCOM.2009.5425809 OpenURL

  23. K Zarif, S Affes, A Ghrayeb, Collaborative null-steering beamforming for uniformly distributed wireless sensors networks. IEEE Trans Signal Process 58(3), 1889–1903 (2010)

  24. I Budiarjo, MK Lakshmanan, H Nikookar, Cognitive radio dynamic access technique. Wirel Person Commun 45(3), 293–324 (2008). Publisher Full Text OpenURL

  25. X Lian, H Nikookar, LP Ligthart, Efficient radio transmission with adaptive and distributed beamforming for intelligent WiMAX. J Wirel Person Commun 53(3), 405–431 (2011)

  26. H Steyskal, Simple method for pattern nulling by phase perturbation. IEEE Trans Antenna Propagat 31(1), 163–166 (1983). Publisher Full Text OpenURL

  27. RJ Mailloux, Covariance matrix augmentation to produce adaptive array pattern troughs. Electron Lett 31(10), 771–772 (1998)

  28. M Zatman, Production of adaptive array troughs by dispersion synthesis. Electron Lett 31(25), 2141–2142 (1998)

  29. JR Guerci, Theory and application of covariance matrix tapers for robust adaptive beamforming. IEEE Trans Signal Process 47(4), 977–985 (1999). Publisher Full Text OpenURL

  30. J Riba, J Goldberg, G Vazquez, Robust beamforming for interference rejection in mobile communications. IEEE Trans Signal Process 45(1), 271–275 (1997). Publisher Full Text OpenURL

  31. P Billingsley, Probability and Measure (John Wiley & Sons, New York, NY, 1995)