Abstract—The analysis of electromagnetic environment using
deterministic mathematical models is characterized by the
impossibility of analyzing a large number of interacting network
stations with a priori unknown parameters, and this is characteristic,
for example, of mobile wireless communication networks. One of the
tasks of the tools used in designing, planning and optimization of
mobile wireless network is to carry out simulation of electromagnetic
environment based on mathematical modelling methods, including
computer experiment, and to estimate its effect on radio
communication devices. This paper proposes the development of a
statistical model of electromagnetic environment of a mobile
wireless communication network by describing the parameters and
factors affecting it including the propagation channel and their
statistical models.
Keywords—Electromagnetic Environment, Statistical model,
Wireless communication network.
I. INTRODUCTION
HE electromagnetic environment (EME) is the sum of
electromagnetic fields from various radio communication
devices and natural electromagnetic processes at a given
geographical location in space. The intensive development of
radio communication systems has led to a significant
concentration of radio communication devices which are
sources of electromagnetic radiation, especially in large cities.
The consequence of this is the increase of both intra-system
and intersystem interferences and the complication of
electromagnetic environment. The ability of a radio
communication device to adequately function at a given
location is completely defined by the EME and its
characteristics and specifications. The prospects of radio
communication systems development, including mobile
wireless communication networks, depend to a considerable
degree on correct and rational planning. However, the
technological development of radio communication systems
planning lags behind the developmental rate of these systems,
and this complicates further systems development and
consequently leading to accumulation of planning errors.
C. Temaneh Nyah is with the Electronics and Telecommunications
Engineering Department, University of Namibia P.O.Box 3624 Ongwediva,
Namibia (phone: +264652324111; fax: +264652324071; e-mail:
cntemaneh@ unam.na, clementtemaneh@yahoo.com).
The available tools for frequency terrestrial planning
considerably simplify and increase the effectiveness of mobile
wireless network design, planning and optimization processes.
One of the tasks of the above tools is to use mathematical
models of radio transmitter stations emissions, radio receiver
stations receptivity, antenna feeder units, radio wave
propagation, and different noise and interference mechanisms
to simulate the electromagnetic environment (EME) and to
estimate its effect on network stations. In general, it is
necessary to evaluate the aggregate action of many
independent signals on the network stations characterized by
different operational structures and algorithms, and also by
the presence of a set of random parameters like their random
number, their operating time, the position of mobile stations,
random physical processes in the radio channel. The principal
unavoidable, drawback in the deterministic approach [1 – 4]
to the description of EME is the impossibility of analyzing a
large number of interacting network stations with a priori
unknown parameters, and this is characteristic, for example,
of mobile wireless communication networks.
The statistical approach to the description of EME is given
in [5-7]. It is based on defining the statistical distribution of
the parameters of network stations (coordinate, frequency, the
power of radiations, etc.), the calculation of the statistical
characteristics of electromagnetic environment and statistical
evaluation by analytical methods of the action of
electromagnetic environment on the network stations. The
main disadvantage of the above mentioned works is the
essential simplification in the models of the distribution of the
stations random parameters for the purpose of obtaining their
statistical characteristics by analytical methods, which in
practice leads to incorrect statistical conclusions.
Therefore, developing a statistical model of EME while
adequately accounting for the set of the random parameters of
the network stations in practice is impossible without the use
of special statistical methods, one, of which is the Monte
Carlo method. By describing the parameters and factors
affecting the electromagnetic environment, including the
propagation channel and their statistical models, the well
known Monte-Carlo technique [8] can be used to develop a
statistical model of EME according to the steps shown in
figure 1 and this is the main focus of this paper.
Developing a Statistical Model for
Electromagnetic Environment for Mobile Wireless
Networks
C. Temaneh Nyah
T
World Academy of Science, Engineering and Technology 61 2012
744
Fig. 1 Structure for the formation of statistical model of EME
The remainder of the paper is organized as follows: The
next section considers the traffic model developed using the
Erlang-B traffic model. Section III outlines how the statistical
models used in describing the MS and BS parameters are
developed. In section IV, an algorithm of obtaining the
statistics of the propagation loss by modeling the slow and
fast fading as log-normal and Nakagami distribution
respectively is proposed. In section V, the structure for
developing the statistical model of EME is presented.
II. TRAFFIC MODEL
The Erlang-B model probability ( )AnP , of occupation of
n channels for a given BS load A is given by the expression
( ) ,,
!
!
,
0
MSn
m
m
n
Nn
m
A
n
A
AnP ≤=
∑
=
(1)
Where MSN — Maximum number of MS, providing a
specified blocking probability value for a given BS load, A.
When MSNn = , the Erlang-B model probability ( )AnP ,
equals the blocking probability.
The recurrent expression (1) can be expressed in the form
represented by expression (2) which is used to compute the
required number of MS MSN for a given blocking probability.
( )
( )
( )
( )
=
−
+
−
=
1,0
,1.1
,1.
,
AP
n
AnPA
n
AnPA
AnP
, (2)
The distribution of active channels when the maximum
number of BS channels 16
max
=N
and load of 7,5,3=A
is shown in figure 2.
Fig. 2 Distribution of active channels when 16
max
=N and
7,5,3=A
The distribution of active channels for different values of
max
N for a given value of A is shown in Fig. 3.
Fig. 3 Distribution of active channels when ErlA 5= and
16,8,4
max
=N
The probability of occupying n out of
max
N channels for
different values of A are shown in Fig. 4.
Fig. 4 Distribution of active channels when 16
max
=N and
7,5,3=A
0 2 4 6 8 10 12 14 16
0
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0.16
0.18
Number of channels
Pr
ob
ab
ilit
y
0 2 4 6 8 10 12 14 16
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
Number of channels
Pr
ob
ab
ilit
y
Nmax = 4
Nmax = 8
Nmax = 16
0 2 4 6 8 10 12 14 16
0
0.05
0.1
0.15
0.2
0.25
Number of channels
Pr
ob
ab
ilit
y
A = 3
A = 5
A = 7
Input network and system data
Traffic model
Model of MS and BS parameters
Frequency
distribution
MS
Spatial
distribution
Transmitter
emission
model
Antenna
feeder
device
model
Statistical model of the radio communication channel
Model of EME
World Academy of Science, Engineering and Technology 61 2012
745
The total number of MS MSN is determined by summing
all the maximum allowable number of MS for all the BS or
BS sectors in the network as follows:
∑
=
=
BSN
i
islotsMS NnN
1
, (3)
slotsn - number of time slots in each frequency channel (used
for TDMA network). For FDMA and CDMA networks,
1=slotsn . Note that the value of MSN must not be greater
than the number of grid points in the region to be analysed.
III. MODEL OF MS AND BS PARAMETERS
A. Frequency Distribution
Cellular network based on FDMA / TDMA technology
consists of a collection of BS grouped in clusters. A cluster
represents a group of BS, in which any frequency channel is
used by only one BS.
The transmission frequency for a given signal or
interference transmitter is obtained as a random variable from
a set of possible frequencies { }nfff ...,, 21 for that transmitter
according to a discrete uniform distribution [7].
B. Mobile Stations Spatial Distribution
In any given region of analysis, the spatial coordinates of
the location of the MS (latitude, longitude) are defined as
pairs of independent random numbers ( )ii yx , distributed
over a uniform distribution with probability density functions:
( )
( )
<<
−
><
=
<<
−
><
=
maxmin
minmax
maxmin
maxmin
minmax
maxmin
1
,,0
1
,,0
yyyfor
yy
yyyyfor
yf
xxxfor
xx
xxxxfor
xf
i
ii
i
i
ii
i
, (4)
In a city, however, there is non-uniform distribution of MS
(users) depending on the time of the day. During working
hours, the maximum number of users is located at the city
business center, while in the evening they are redistributed in
the living areas of the city. Therefore, in this case we use the
exponential distribution model of MS [9]. The probability
density function in this case is express as:
( ) ( )
<<
<<
=
−−
M
Rr RrRfore
Rrfor
rf
ασ
σ
0
0 0,
, (5)
where
0σ –– user’s constant density at the center of the region of
analysis;
R –– radius of the circular part defining the center of the
region of analysis;
α –– parameter defining the decrease of the number of users
with distance from center.
This model is most appropriate to describe the distribution
of MS within a city with radial structure of buildings.
As an example, the city of Moscow has the following
parameters [9]; 20 5.168 kmErl=σ , kmR 45.1= , 48.0=α .
C. Transmitter Model
The transmitter emission’s power in practice is not constant
but depends on many factors, amongst which includes the
frequency dependence of the channels being used to its
instability, variations in transmitter supply voltage, the impact
on the transmitter output stage of the radiation of other closely
located transmitters etc. Some of these factors are random,
and this is confirmed by experimental measurements [1].
Therefore, in the description of the transmitter emission’s
power a statistical distribution of power is used. The statistical
model of the BS transmitter’s main emission mainBSP , is based
on the assumption of normal distribution of power
( ) ( )
−
−= 2
2
,
,
, 2
exp
2
1
P
mainBSmainBS
P
mainBS
PP
Pf
σσpi
and is defined as:
pPP mainBSmainBS ∆+= ,, , (6)
where, mainBSP , –– mean value of BS main emission power
indicated in the technical documentation of transmitters; p∆
–– Gaussian variable with zero mean.
The statistical model of the BS side frequency emission at
harmonics
sideBSP , is defined as:
sideBSmainBSsideBS XPP ,,, += , (7)
where, sideBSX , –– attenuation (relative to carrier’s power) of
the BS side frequencies emissions at harmonics.
Usually the emissions at the side frequencies at harmonics of
the fundamental frequency have the highest power level
compared to the emissions of other side frequencies [10].
Therefore, in most cases, the emissions of the other side
frequencies can be neglected.
D. Antenna Model
Considering that there are no perfect matching feeders with
the antenna, and also that the antenna radiation pattern is
affected by various surrounding objects, the BS transmitter
antenna gain BSG is defined as a random variable with
normal distribution express as:
gGG BSBS ∆+= , (8)
were, g∆ –– Gaussian variable with zero mean and standard
deviation ––
aσ , numerically equal to the value the spread,
which is contained in the technical characteristics of antennae;
BSG –– average gain of BS transmitter antenna and can be
expressed by the formula:
HCGG BS ++= 0 , (9)
0G –– Receiver antenna gain for the working frequency
range and given polarization with radiation pattern taken into
account;
C –– correction for the frequency dependence of antenna
pattern (for the frequencies of side emissions and receiver
World Academy of Science, Engineering and Technology 61 2012
746
spurious channel); H – correction for different polarizations
of the transmitting (interference) and receiving antennae.
IV. CHANNEL MODEL
The well known propagation models such as Hata model,
ITU-R P.529-3 and ERC Report 68 with sufficient certainty
predict the median loss for mobile communications systems in
urban areas, but do not give the statistics, due to the effects of
slow and fast fading. We propose a technique of obtaining the
statistics of propagation losses in the channel by numerical
simulation, taking into account the effects of both slow
(lognormal distribution) and fast (Rayleigh, Rice) fading thus
improving the accuracy of the estimated losses.
For the case when the lognormal distribution, the
Nakagami distribution characterizes the slow and fast fading
respectively, the loss statistics is obtained by the generation
of the instantaneous value of propagation loss ijL from the j -
th BS to the location of the i -th MS is as follows:
1. Calculate the median value of the loss ijmL from BS with
index j to the location of MS with index i .
2. Generate the slowly varying local mean value of the loss
sL according to normal distribution with mean ijmL and
given standard deviation
sLσ .
3. Generate the instantaneous value of propagation loss gL
according to gamma distribution with fading parameter m
and scale parameter ( ) ( )( )[ ]2101021 sLmmm ⋅+ΓΓ=Ω .
4. Calculate the instantaneous loss by the formula:
( )gij LL 10log10= (10)
V. CONCLUSION
The proposed procedure for developing a statistical model
for EME can be used in the algorithm for the statistical
estimation of electromagnetic compatibility of a
communication network and also in organisations ranging
from telecommunication authorities and network planners
carrying out services requiring computation in software tools
for Radio Network Planning, Optimisation of Radio Networks
and Spectrum Management.
REFERENCES
[1] Donald R. J. White. Electromagnetic compatibility of radio electronic
equipment and unwanted interference. Moscow: Soviet Radio, 1977.
[2] Knyazev A.D. Elements of the theory and practice for electromagnetic
compatibility of radio electronic equipment. - M.: Radio and
Communications, 1984.
[3] Vladimirov V.I, Doctorov A.L, Elizarov F. V. et al. Electromagnetic
compatibility of radio electronic equipment and systems. Edited by
N.M. Tsarkova. - M.: Radio and Communications, 1985.
[4] Petrovsky, V. I., Sedelnikov U.E. Electromagnetic compatibility of
radio electronic equipment: A manual for higher institutions. M.: Radio
and Communications, 1986.
[5] ITU-R M.1634.
[6] Alter L.S. The probability of intermodulation interference in mobile
radio communication systems receivers. Mobile systems. 2003. 12.
pp.55-58.
[7] Aporovich A.F. The statistical theory of electromagnetic compatibility
of radio electronic equipment, Edited by V.Y. Averyanova. - Mn.:
Science and Technology, 1984.
[8] Nyah C. Temaneh. Monte-Carlo technique estimation of a probability of
intermodulation interference in a cellular wireless communication
network. Proceedings from the IEEE Region 8 International Conference
on Computational Technologies in Electrical and Electronics
Engineering SIBIRCON-2010, Irkutsk – Russia, Volume 1, pp. 329-33.
[9] A. Dementiev. Reliable methods of estimating the spatial distribution of
the load in cellular networks. Telecommunication journal № 2 in 2002.
[10] Petrovsky V.I., Sedelnikov U.E. Electromagnetic compatibility of radio
electronic devices. Moscow: Radio and communication, 1986.
C. Temaneh Nyah (M’2009) was born in Cameroon. He received the
Master of Science (MSc.) in Engineering from the Faculty of Radio
communication, Radio broadcasting and Television at the Ukrainian State
Telecommunications Academy – Odessa in 1995 and a Ph.D. degree in
Telecommunications from the Yaroslavl State University – Russia in 2007.
Presently he is Senior lecturer at the University of Namibia. His research
interest is on Telecommunication’s network performance simulation models,
for network management, planning, optimisation and analysis including
electromagnetic compatibility analysis.
World Academy of Science, Engineering and Technology 61 2012
747