
Advanced Wireless Networks - 4G Technologies
.pdf
48 |
CHANNEL MODELING FOR 4G |
|
|
|
Table 3.1 Summary of macro cellular measurement environments |
||
|
|
|
|
Class |
|
BS antenna height |
Description of environment |
|
|
|
|
TU, typical urban |
10 and 32 m |
The city of Aarhus, Denmark. Uniform |
|
|
|
|
density of buildings ranging from four to |
|
|
|
six floors. Irregular street layout. |
|
|
|
Measurements were carried out along six |
|
|
|
different routes with an average length of |
|
|
|
2 km. No line-of-sight between MS and BS. |
|
|
|
MS–BS distance varies from 0.2 to 1.1 km |
TU |
|
21 m |
Stockholm city, Sweden (area 1). Heavily |
|
|
|
built-up area with a uniform density of |
|
|
|
buildings, ranging from four to six floors. |
|
|
|
Ground is slightly rolling. No line-of-sight |
|
|
|
between MS and BS. MS–BS distance |
|
|
|
varies from 0.2 to 1.1 km |
BU, bad urban |
21 m |
Stockholm city, Sweden (area 2). Mixture of |
|
|
|
|
open flat areas (river) and densely built up |
|
|
|
zones. Ground is slightly rolling. No |
|
|
|
line-of-sight between MS and BS. MS–BS |
|
|
|
distance varies from 0.9 to 1.6 km |
SU suburban |
12 m |
The city of Gistrup, Denmark. Medium-sized |
|
|
|
|
village with family houses of one to two |
|
|
|
floors and small gardens with trees and |
bushes. Typical Danish residential area. The terrain around the village is rolling with some minor hills. No line-of-sight between MS and BS. MS–BS distance varies from 0.3 to 2.0 km
along a certain route. The local average power azimuth-delay spectrum is given as
P (φ, τ ) = E |
L |
|
|αl |2 δ (φ − φl , τ − τl ) |
(3.2) |
|
|
l=1 |
|
From Equation (3.2), local power azimuth spectrum (PAS) and the local power delay spectrum (PDS) are given as
PA (φ) = |
P (φ, τ ) dτ |
(3.3) |
PD (τ ) = |
P (φ, τ ) dφ |
(3.4) |
The radio channels’ local azimuth spread (AS) σA and the local delay spread (DS)σD are defined as the root second central moments of the corresponding variables. The values of the local AS and DS are likely to vary as the MS moves within a certain environment. Hence, we can characterize σA and σD as being random variables, with the joint pdf f (σA, σD).

MACROCELLULAR ENVIRONMENTS (1.8 GHz) |
49 |
Their individual PDFs are
fA(σA) = |
f (σA, σD) d σD |
(3.5) |
fD (σD) = |
f (σA, σD) d σA |
(3.6) |
The function f (σA, σD) can be interpreted as the global joint PDF of the local AS and DS. If the expectation in Equation (14.2) is computed over the radio channel’s fast fading component, we can furthermore apply the approximation
P (φ, τ ) dφ dτ |
= |
= |
(3.7) |
|
hchannel hloss (d) hs |
where hchannel is the radio channel’s integral path loss, hloss (d) is the deterministic long-term distance-dependent path loss, while hs is the channel’s shadow fading component, which is typically modeled by a log–normal distributed random variable [5.6]. The global PDF of hs is denoted fs (hs). The global degree of shadow fading is described by the root second central moment of the random shadow fading component expressed in decibel, i.e.
σs = Std 10 log10 (hs) |
(3.8) |
where Std {} denotes the standard deviation. Empirical results for cumulative distribution functions (CDF) for σA and σD are given in Figures 3.1 and 3.2, respectively. The log–normal fit for σA results is given as
|
σA = 10εA X+μA |
where X |
is a zero-mean Gaussian distributed random variable with unit |
μA = E |
log10 (σA) is the global logarithmic mean of the local AS, and |
log10 (σA) is the logarithmic standard deviation of the AS.
(3.9)
variance, εA = Std
Cumulative distribution
1.0 |
|
|
|
|
|
|
0.9 |
|
|
|
|
|
|
0.8 |
|
|
|
|
|
|
0.7 |
|
|
|
|
|
|
0.6 |
|
|
|
|
|
|
0.5 |
|
|
|
CDF of the AS |
|
|
0.4 |
|
|
|
Empirical results TU-32 |
|
|
|
|
|
|
|
||
0.3 |
|
|
|
Log--normal distribution |
|
|
|
|
|
Empirical results BM |
|
||
|
|
|
|
|
||
0.2 |
|
|
|
Log--normal distribution |
|
|
|
|
|
|
|
|
|
0.1 |
|
|
|
|
|
|
0.0 |
5 |
10 |
15 |
20 |
25 |
30 |
0 |
Azimuth spread (degree)
Figure 3.1 Examples of empirical CDF of AS obtained in different environments. The CDF of a log–normal distribution is fitted to the empirical results for comparison. (Reproduced by permission of IEEE [4].)

50 |
CHANNEL MODELING FOR 4G |
|
|
|
|
|
|
|||
|
|
1.0 |
|
|
|
|
|
|
|
|
|
|
0.9 |
|
|
|
|
|
|
|
|
|
distribution |
0.8 |
|
|
|
|
|
|
|
|
|
0.7 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.6 |
|
|
|
CDF of the DS |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Cumulative |
0.5 |
|
|
|
|
|
|
|
|
|
0.3 |
|
|
|
|
|
|
|
|
|
|
|
0.4 |
|
|
|
|
Empirical results TU-32 |
|
|
|
|
|
|
|
|
|
|
Log--normal distribution |
|
||
|
|
|
|
|
|
|
Empirical results BM |
|
|
|
|
|
0.2 |
|
|
|
|
Log--normal distribution |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
0.1 |
|
|
|
|
|
|
|
|
|
|
0.0 |
|
|
|
|
|
|
|
|
|
|
0 |
0.5 |
1.0 |
1.5 |
2.0 |
2.5 |
3.0 |
3.5 |
4.0 |
Delay spread (μs)
Figure 3.2 Examples of empirical CDFs of the DS in different environments. The CDF of a log–normal distribution is fitted to the empirical results for comparison. (Reproduced by permission of IEEE [4].)
Table 3.2 Summary of the first and second central moments of the AS, DS, and shadow fading in the different environments (reproduced by permission of IEEE [4])
Class |
σs |
E {σA} |
μA |
εA |
E {σD} |
μD |
εD |
TU-32 |
7.3 dB |
0 |
0.74 |
0.47 |
0.8 μs |
−6.20 |
0.31 |
80 |
|||||||
TU-21 |
8.5 dB |
80 |
0.77 |
0.37 |
0.9 μs |
−6.13 |
0.28 |
TU-21 |
7.9 dB |
130 |
0.95 |
0.44 |
1.2 μs |
−6.08 |
0.35 |
BU |
10.0 dB |
70 |
0.54 |
0.60 |
l.7 μs |
−5.99 |
0.46 |
SU |
6.1 dB |
8 |
0.84 |
0.31 |
0.5 μs |
−6.40 |
0.22 |
Similarly, |
|
|
|
|
|
|
σD = 10εDY +μD |
|
(3.10) |
where Y |
is a |
zero-mean Gaussian distributed random |
variable with unit |
variance, |
μD = E |
log10 (σD) is the global logarithmic mean |
of the local DS, |
and εD = |
|
Std log10 (σD) |
is the logarithmic standard deviation of the DS. A summary of the re- |
sults for these parameters is given in Table 3.2. For characterization of shedowing fading see References [1–30].
3.2 URBAN SPATIAL RADIO CHANNELS IN MACRO/MICROCELL
ENVIRONMENT (2.154 GHz)
The discussion in this section is based on the experimental results collected with a wideband channel sounder using a planar antenna array [31]. The signal center frequency was 2154 MHz and the measurement bandwidth was 100 MHz. A periodic PN-sequence, 255 chips

|
|
|
RADIO CHANNELS IN MACRO/MICROCELL ENVIRONMENT (2.154 GHz) |
51 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
TX11 |
|
|
|
|
|
|
|
|
100 m |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
W |
|
|
|
|
|
E |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Building height |
|
|
S |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
Theatre tower |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
> 30 m |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
RX2 |
|
|
|
|||
|
26-30 m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
21-25 m |
|
|
|
|
|
Station tower |
|
|
|
|
|
||||||
|
10-20 m |
|
|
|
|
|
|
|
|
|
RX3 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Transmitter
positions |
|
Railway station |
|
RX1 |
|
|
|
|
|
|
|
|
|
RX2 |
|
|
|
|
|
|
|
|
|
|
|
RX1 |
|
|
|
|
|
|
|
RX3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Cathedral |
|
||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
TX17 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
TX28 |
Hotel Torni
Figure 3.3 The measurement area with all three RX sites; TX-positions of the sample plots are marked. (Reproduced by permission of IEEE [31].)
long, was used. The chip rate was 30 MHz, the sampling rate 120 MHz, giving a oversampling factor 4. The correlation technique was used for the determination of the impulse response. Hence, the delay range was 255/30 MHz = 8.5 μs, with a resolution of 1/30 MHz = 33 ns. The transmit antenna at the MS was a vertically polarized omnidirectional discone antenna. The vertical 3 dB beamwidth was 87◦ and the transmit power 40 dBm. Approximately 80 different transmitter positions were investigated.
The receiving BS was located at one of three different sites below, at and above the rooftop level (RX1–RX3, see Figure 3.3). A 16-element physical array with dual polarized λ/2-spaced patch antennas was combined with a synthetic aperture technique to build a virtual two-dimensional (2-D) antenna structure. The patches were linearly polarized at 0◦ (horizontal direction) and 90◦ (vertical direction). With these 16 × 62 elements the direction of arrival (DOA) of incoming waves both in azimuth (horizontal angle) and elevation (vertical angle) could be resolved using the super-resolution Unitary ESPRIT algorithm [32–34]. Note that the number of antenna elements limits the number of identifiable waves, but not the angular resolution of the method. Together with a delay resolution of 33 ns, the radio channel can be characterized in all three dimensions separately for the two polarizations. Array signal processing, including estimation of the DOAs and a comparison of ESPRIT with other algorithms, can be found in Glisic [1].
One prerequisite for the applicability of the synthetic aperture technique is that the radio channel is static during the whole data collection period. To avoid problems, the whole procedure was done at night with minimum traffic conditions.
3.2.1 Description of environment
A typical urban environment is shown in Figure 3.3 [31] with three receiver locations (RX1–RX3) marked by triangles pointing in the broadside direction of the array. Figure 3.3 also shows all the corresponding TX positions. The location RX 1 (height hRX = 10 m)
52 CHANNEL MODELING FOR 4G
was a typical microcell site below the rooftop height of the surrounding buildings, and measurements were performed with 20 different TX positions. RX 2 (height hRX−27 m) was at the rooftop level, and 32 TX positions were investigated. RX 3 (height hRX = 21 m) was a typical macrocell BS position above rooftop heights, and 27 TX positions were measured.
3.2.2 Results
The measurement results show that it is possible to identify many single (particular, different) multipath components, impinging at the receiver from different directions. However, these components are not randomly distributed in the spatial and temporal domain; they naturally group into clusters. These clusters can be associated with objects in the environment due to the high angular and temporal resolution of evaluation. (Sometimes even individual waves, within a cluster, can be associated with scattering objects.) The identification of such clusters is facilitated by inspection of the maps of the environment. A cluster is defined as a group of waves whose delay, azimuth and elevation at the receiver are very similar, while being notably different from other waves in at least one dimension. Additionally all waves inside a cluster must stem from the same propagation mechanism. The definition of clusters always involves a certain amount of arbitrariness. Even for mathematically ‘exact’ definitions, arbitrary parameters (e.g. thresholds or number of components) must be defined. Clustering by human inspection, supported by maps of the environment, seems to give the best results. The received power is calculated within each cluster (cluster power) by means of unitary ESPRIT and a following beam-forming algorithm. The results are plotted in the azimuth-elevation-, azimuth-delay- and elevation-delay-planes
According to the obvious propagation mechanism, each cluster is assigned to one of three different classes:
Class 1, street-guided propagation – waves arrive at the receiver from the street level after traveling through street canyons.
Class 2, direct propagation over the rooftop – the waves arrive at the BS from the rooftop level by diffraction at the edges of roofs, either directly or after reflection from buildings surrounding the MS. The azimuth mostly points to the direction of the transmitter with some spread in azimuth and delay.
Class 3, reflection from high-rise objects over the rooftop – the elevation angles are near the horizon, pointing at or above the rooftop. The waves undergo a reflection at an object rising above the average building height before reaching the BS. The azimuth shows the direction of the reflecting building; the delay is typically larger than for class 1 or class 2. The sum of the powers of all clusters belonging to the same class is called class power. In some cases the propagation history is a mixture of different classes, e.g. street guidance followed by diffraction at rooftops. Such clusters are allocated to the class of the final path to the BS.
For the evaluation of delays we define the vector P containing the powers of the clusters and vector τ of corresponding mean delays. A particular cluster i has mean delay τi , and power Pi . The relation between the delays τ and the powers P is modeled as exponential:
Pn P(τn ) = ae−τn /b |
(3.11) |

MIMO CHANNELS IN MICROAND PICOCELL ENVIRONMENT (1.71/2.05 GHz) |
53 |
Table 3.3 The model parameters a and b for both received polarizations (VP and HP) averaged over all available clusters. The transmitter was VP
|
VP-VP |
VP-HP |
|
|
|
a |
−3.9 dB |
−3.6 dB |
b |
8.9 dB/μs |
11.8 dB/μs |
|
Table 3.4 Average delay of the strongest cluster of RX1, RX2 and RX3 |
|
|
|
|
RX |
Average delay VP-VP μs |
Average delay VP-HP μs |
|
|
|
1 |
0.068 |
0.071 |
2 |
0.38 |
0.28 |
3 |
0.11 |
0.048 |
|
|
|
Experimental data are fit into the model (3.11) using the least square (LS) estimation. The logarithmic estimation error v is defined as
v = 10 log P − 10 log s(θ) |
(3.12) |
||
and its standard deviation σv as |
|
||
σv = |
|
|
(3.13) |
var{v} |
|||
The summary of the results for parameters a and b is shown in Tables 3.3 and 3.4. |
|
||
Pi = ae−τi /b |
|
In equation (14.21) σv was defined as the standard deviation of the logarithmic estimation error in dB. This estimation error was found to be log–normally distributed and, up to a delay of about 1 μs, σv is independent of the delay τ . The value of σv is 9.0 and 10.0 dB (co-and cross-polarization), respectively, averaged over the first microsecond. The average powers for different classes of clusters are shown in Table 3.5. Additional data on the topic can be found in References [31–41].
3.3 MIMO CHANNELS IN MICROAND PICOCELL ENVIRONMENT
(1.71/2.05 GHz)
The model presented in this section is based upon data collected in both picocell and microcell environments [43]. The stochastic model has also been used to investigate the capacity of MIMO radio channels, considering two different power allocation strategies, water filling and uniform and two different antenna topologies, 4 × 4 and 2 × 4. It will be demonstrated that the space diversity used at both ends of the MIMO radio link is an efficient technique in picocell environments, achieving capacities within 14 and 16 b/s/Hz in 80 % of the cases for a 4 × 4 antenna configuration implementing water filling at a signal-to-noise ratio (SNR) of 20 dB.

54 |
CHANNEL MODELING FOR 4G |
|
||
|
|
Table 3.5 Averaged class powers of RX1, RX2 and RX3 |
||
|
|
|
|
|
|
|
|
Horizontal power, percentage |
Vertical power, percentage |
|
RX |
Class |
of total power |
of total power |
|
|
|
|
|
1 |
1 |
96.5 % |
95.7 % |
|
|
|
2 |
2.4 % |
3.8 % |
|
|
3 |
1.1 % |
0.4 % |
2 |
1 |
93.5 % |
97.2 % |
|
|
|
2 |
4.0 % |
2.7 % |
|
|
3 |
2.5 % |
0.1 % |
3 |
1 |
46.7 % |
78.0 % |
|
|
|
2 |
37.2 % |
12.8 % |
|
|
3 |
16.0 % |
9.2 % |
|
|
|
|
|
The basic parameters of the measurements set-up are shown in Figure 3.4. The following notaion is used in the figure: dMS−DS stands for distance between MS and BS; hBSfor the height of BS above ground floor, and AS for azimuth spread [43].
The vector of received signals at BS can be represented as y(t) = [y1(t), y2(t), . . . , yM (t)]T, where ym (t) is the signal at the mth antenna port and [·]T denotes transposition. Similarly, the signals at the MS are s (t) = [s1(t), s2(t), . . . sN (t)]T. The NB MIMO radio channel H X M×N , which describes the connection between the MS and the BS, can be
|
MIMO measurements |
|
|
|
1st measurement set-up |
2nd measurement set-up |
|
|
Rotative mMotion |
|
Linear mMotion |
Tx |
0..5λ |
Tx |
0.5λ |
(MS) |
Fc=1.71 GHz |
(MS) |
|
|
|
Fc=2.05 GHz |
|
|
|
|
|
|
|
Switch |
|
20λ |
Synthesized |
approach |
|
|
11.8λ |
||
|
approach |
|
|
Rx |
Rx |
|
|
|
|
||
(BS) |
(BS) |
|
|
|
0.5λ |
|
1.5λ |
|
|
|
|
|
15 positions |
|
93 positions |
|
(microcell) |
(microcell and picocell) |
Figure 3.4 Functional sketch of the MIMO model. (Reproduced by permission of IEEE [43].)
MIMO CHANNELS IN MICROAND PICOCELL ENVIRONMENT (1.71/2.05 GHz) |
55 |
expressed as [1]:
|
α11 |
α12 |
· · · |
α1N |
|
|
H |
α21 |
α22 |
· · · |
α2N |
(3.14) |
|
. |
. |
. . |
|
. |
||
|
= . |
. |
. |
. |
|
|
|
. |
. |
|
. |
|
|
|
αM1 |
αM2 |
· · · |
αM N |
|
where αmn is the complex transmission coefficient from antenna at the MS to antenna at the BS. For simplicity, it is assumed that αmn is complex Gaussian distributed with identical average power. However, this latest assumption can be easily relaxed. Thus, the relation between the vectors y (t) and s (t) can be expressed as
y (t) = H (t) s (t)
In the sequel we will use the following correlations:
|
|
BS |
= αm1n , αm2n |
|
||||
|
|
ρm1m2 |
|
|||||
|
|
MS |
= αmn1 , αmn2 |
|
||||
|
|
ρn1n2 |
|
|||||
|
|
ρ11BS |
ρ12BS |
· · · |
ρ1BSM |
|
||
RBS |
|
ρ21BS |
ρ22BS |
· · · |
ρ2BSM |
|
||
|
. |
. |
. . |
|
|
. |
|
|
|
= . |
. |
. |
|
. |
|
||
|
|
. |
. |
|
|
. |
|
|
|
|
BS |
BS |
· · · |
BS |
M×M |
||
|
|
ρM1 |
ρM2 |
ρM M |
||||
|
|
ρ11MS |
ρ12MS |
· · · |
ρ1MSN |
|
||
RMS |
= |
ρ21MS |
ρ22MS |
· · · |
ρ2MSN |
|
||
. |
. |
. . |
|
. |
|
|||
|
. |
. |
. |
. |
|
|||
|
|
. |
. |
|
|
. |
|
|
|
|
ρNMS1 |
ρNMS2 |
· · · |
ρNMSN |
N ×N |
(3.15)
(3.16)
(3.17)
(3.18)
(3.19)
The correlation coefficient between two arbitrary transmission coefficients connecting two different sets of antennas is expressed as
n1m1 |
|
, αm2n2 |
(3.20a) |
ρn2m2 = αm1n1 |
|||
which is equivalent to |
|
|
|
n1m1 |
M S |
B S |
(3.20b) |
ρn2m2 |
= ρn1n2 |
ρm1m2 |
provided that Equations (14.24a) and (14.24b) are independent of n and m, respectively. In other words, this means that the spatial correlation matrix of the MIMO radio channel is the Kronecker product of the spatial correlation matrix at the MS and the BS and is given by
RMIMO = RMS RBS |
(3.21) |
where represents the Kronecker product. This has also been confirmed in yu et al. [44].

56 CHANNEL MODELING FOR 4G
3.3.1 Measurement set-ups
The TX is at the MS and the stationary RX is located at the BS. The two set-ups from Figure 3.4 provide measurement results with different correlation properties of the MIMO channel for small antenna spacings of the order of 0.5λ or 1.5λ . The BS consists of four parallel RX channels. The sounding signal is a MSK-modulated linear shift register sequence of a length of 127 chips, clocked at a chip rate of 4.096 Mcs. At the RX, the channel sounding is performed within a window of 14.6 μs, with a sampling resolution of 122 ns (half-chip period) to obtain an estimate of the complex IR. The narowband (NB) information is subsequently extracted by averaging the complex delayed signal components. A more thorough description of the stand-alone testbed (i.e. RX and TX) is documented in References [23, 45]. The description of the measurement environments is summarized in Table 3.6.
A total of 107 paths were investigated within these seven environments. The first measurement set-up was used to investigate 15 paths in a microcell environment, i.e. environment A in Table 3.6. The MS was positioned in different locations inside a building while the BS was mounted on a crane and elevated above roof-top level (i.e. 9 m) to provide direct line-of-sight to the building. The antenna was located 300 m away from the building. The second set-up was used to investigate 92 paths for both microcell and picocell environments,
Table 3.6 Summary and description of the different measured environments (reproduced by permission of IEEE [43])
Cell type |
Environment |
MS locations |
Measurement set-up |
Description |
|
|
|
|
|
Microcell |
A |
15 |
1st |
The indoor environment |
|
|
|
|
consists of small |
|
|
|
|
offices with windows |
|
|
|
|
metallically shielded |
|
|
|
|
–300 m between MS |
|
|
|
|
and BS |
|
B |
13 |
2nd |
The indoor environment |
|
|
|
|
consists of small |
|
|
|
|
offices –31–36 m |
|
|
|
|
between MS and BS |
Picocell |
C |
21 |
2nd |
The indoor environment |
|
|
|
|
is the same as in A |
|
D |
12 |
2nd |
Reception hall – large |
|
|
|
|
open area |
|
E |
18 |
2nd |
Modern open office with |
|
|
|
|
windows metallicaly |
|
|
|
|
shielded |
|
F |
16 |
2nd |
The indoor environment |
|
|
|
|
is the same as in B |
|
G |
12 |
2nd |
Airport – very large |
|
|
|
|
indoor open area |
|
|
|
|
|

MIMO CHANNELS IN MICROAND PICOCELL ENVIRONMENT (1.71/2.05 GHz) |
57 |
i.e. environment B and C–G, respectively, as shown in Table 3.6. The distance between the BS and the MS was 31–36 m for microcell B, with the BS located outside.
3.3.2 The eigenanalysis method
The eigenvalue decomposition (EVD) of the instantaneous correlation matrix R = HHH(not to be confused with RMIMO), where [·]H represents Hermitian transposition, can serve as a benchmark of the validation process. The channel matrix H may offer K parallel subchannels with different mean gains, with K = Rank (R) ≤ min (M, N ) where the functions Rank (·) and min (·) return the rank of the matrix and the minimum value of the arguments, respectively [27]. The kth eigenvalue can be interpreted as the power gain of the kth subchannel [27]. In the following, λk represents the eigenvalues.
3.3.3 Definition of the power allocation schemes
In the situation where the channel is known at both TX and RX and is used to compute the optimum weight, the power gain in the kth subchannel is given by the kth eigenvalue, i.e. the SNR for the kth subchannel equals
Pk |
(3.22) |
γk = λk σN2 |
where Pk is the power assigned to the kth subchannel, λk is the kth eigenvalue and σN2 is the noise power. For simplicity, it is assumed that σN2 = 1. According to Shannon, the maximum capacity normalized with respect to the bandwidth (given in terms of b/s/Hzspectral efficiency) of parallel subchannels equals [46]
K |
|
log2 (1 + γk ) |
|
|
|
|
|||
C = |
|
|
|
|
(3.23) |
||||
k=1 |
|
|
|
|
|
|
|
|
|
K |
|
|
|
|
|
|
Pk |
|
|
= |
|
log2 1 + λk |
|
(3.24) |
|||||
|
|
σN2 |
|
||||||
k=1 |
|
|
|
|
|
|
|
|
|
where the mean SNR is defined as |
|
|
|
|
|
|
|
|
|
SNR = |
E [PRX] |
= |
E |
[PTX] |
(3.25) |
||||
|
σN2 |
|
|
σN2 |
Given the set of eigenvalues {λk }, the power Pk allocated to each subchannel k was determined to maximize the capacity using Gallager’s water filling theorem [27] such that each subchannel was filled up to a common level D, i.e.
1 |
+ P1 = · · · = |
1 |
+ Pk = · · · D |
(3.26) |
λ1 |
λK |
with a constraint on the total TX power such that
K
Pk = PTX |
(3.27) |
k=1