9 Pages
English

L Pisani F Rapetti C Vittoli| A Kirchho integral approach for decimetric radiowave propagation in urban areas

Gain access to the library to view online
Learn more

Description

L. Pisani, F. Rapetti, C. Vittoli| A Kirchho integral approach for decimetric radiowave propagation in urban areas 1 5

  • eld

  • scattered elds

  • computational aspects

  • kirchho integral

  • only over

  • using maxwell

  • building

  • lic buildings

  • plane wave

  • tor kirchho integral


Subjects

Informations

Published by
Reads 9
Language English

decimetric
L
urban
Pisani
v
F
h
Rap
w
etti
propagation
C
1
Vittoli
for
|
radio
A
a
Kirc
e
hho
in
in
areas
tegral
5
approacpages
A
pro
Kirc
(
hho
v
in
!
tegral
is
approac
to
h
discussed
for
ld
decimetric
(2)
radio
(3)
w
t
a
a
v
com
e
determine
propagation
the
in
e
urban
r
areas
By
L
relation
Pisani
4
F
olume
Rap
as
etti
hho
C
to
Vittoli
ha
Cen
approac
tro
F
di
d
Ricerca
urban
Svilupp
tegral
o
Kirc
e
ector
Studi
are
Sup
a
eriori
x
in
v
Sardegna
r
RS
equation
Via
x
N
r
Sauro
0
10,
p
09123
S
Cagliari
.
Italy
the
.
to
Abstract
using
W
h
e
a
consider
er
a
een
threeimensional
complexit
approac
alidated
h
erimen
based
related
on
scaling
the
as
Kirc
co
hho
are
's
ector
metho
section
d
the
in
relation
order
ose
to
harmonic
predict
When
electromagnetic
olume
w
v
a
inside
v
2
es
(1)
propagation
is
in
elo
urban
ergence
en
S
vironmen
to
ts
erid
In
tit
particular
S
w
r
e
)
are
n
in
dx
terested
x
here
is
in
b
the
unit
ev
the
aluation
h
of
in
the
at
electro
has
magnetic
t
ld
metho
on
iterativ
v
where
ery
corresp
large
rection
threeimensional
In
domains
um
(t
appro
ypically
e
with
tro
linear
the
dimensions
.
of
is
the
n
order
with
of
measuremen
h
re
un
the
dreds
and
of
the
meters
w
generated
application
b
an
y
erage
a
vironmen
high
ted
frequency
The
source
hho
(t
In
ypically
e
of
basis
the
ec
order
in
of
kson
Hz
435).
whic
the
h
E
corresp
dep
onds
t
to
ld
a
the
w
,
a
satiss
v
Helmholtz
elength
e
of
,
ab
+
out
E
30
=
cen
k
timeters
and
Some
w
n
propagation
umerical
y
tests
the
and
I
comparisons
=
with
n
exp
is
erimen
that
tal
is
measuremen
follo
ts
an
ha
E
v
=
e
[
b
0
een
e
done

to
ik
v
r
alidate

this
E
approac
(3)
h
=
1
,
In
t
tro
the
duction
b
Eien
the
t
and
planning
n
of
in
mobile
olume
comm
(3),
unication
kno
systems
Kirc
is
gral
based
ld
on
p
an
inside
accurate
led
determination
implemen
of
Kirc
the
's
co
d
v
an
erage
e
re
cedure
gion
eac
of
iteration
the
onds
an
a
tenna
on
In
surface
urban
addition
en
n
vironmen
b
ts
of
this
ximations
can
v
b
b
e
in
quite
duced
complicated
limit
since
computational
electromagnetic
y
w
This
a
h
v
v
es
through
are
umerical
absorb
parisons
ed
exp
rected
tal
transmitted
ts
and
urther
diacted
sults
b
to
y
con
build
ergence
ings
time
The
of
sub
metho
ject
as
of
ell
our
its
researc
to
h
the
activit
tenna
y
v
is
in
the
en
analysis
ts
and
presen
dev
and
elopmen
2
t
v
of
Kirc
sim
in
ulation
relation
to
this
ols
w
to
review
ectiv
theoretical
ely
of
predict
v
the
tor
propagation
hho
of
tegral
the
ac
electromagnetic
1975,
w
432-
a
Supp
v
that
es
v
and
ld
the
has
an
time
tenna
endence
co
i
v
.
erage
the
in
sources
threeimensional
outside
urban
v
domains
V
whic
the
h
E
ha
the
v
ector
e
w
h
v
uge
equation
dimensions
V
with
(
resp
2
ect
k
to
)
the
(
w
)
a
0
v
where
e
=
length
c
T
c
ec
the
hniques
a
used
e
to
v
solv
cit
e
.
this
applying
problem
div
in
theorem
a
V
reason
dV
ably
I
fast
V
w
d
a
it
y
easy
include
see
ra
when
y
(1)
tracing
v
k
the
egami
wing
et
is
al
iden
,
y
1991,
(
Rossi
)
et
I
al
V
,
E
1992,
n
K


0
urner
ik
et
4
al
r
,
e
1993,
r
Rizk

et
(
al
0
,
r
1994)
)
and
]
t
0
w
In
oimensional
r
T
x
ransmission
0
Line
the
Matrix
oin
LM
x
uthi
inside
et
v
al
V
,
ounded
1996).
y
W
surface
e
V
prop
the
ose
normal
a
0
threeimensional
directed
ap
to
proac
v
h
V
based
Equation
on
whic
the
is
Kirc
wn
hho
the
in
hho
tegral
te
metho
expresses
d
electric
The
E
presence
a
of
oin
m
x
ultiple
a
rectionsonly
L
com
Pisani
ld
F
the
Rap
b
etti
b
C
ld
Vittoli
than
|
is
A
the
Kirc
elength
hho
terer
in
ects
tegral
b
approac
to
h
surface
for
b
decimetric
resp
radio
w
w
there
a
or
v
are
e
enerated
propagation
as
in
go
urban
e
areas
to
2
computational
closed
e
v
OS
olume
iterativ
V
al
in
t
terms
using
of
w
the
obstacle
v
illuminated
alues
shado
of
ob
the
equal
ld
the
and
with
its
region
normal
on
deriv
inciden
ativ
aluated
e
B
on
the
the
metho
b
around
oundary
w
surface
Optic
S
the
V
e
.
allo
By
n
using
other
Maxw
on
ell
pro
equations
More
w
one
e
corresp
can
duced
rewrite
n
equation
from
(3)
the
in
second
terms
rected
of
the
the
length
electric
linear
and
of
magnetic
ximately
lds
and
E
s
and
the
B
is
:
W
E
on
(
to
x
e
)
alls
=
of
i
the
4
in

eien
I
aluate
S
This
V
follo
e
i
ik
an
r
S
r
lds

are
k
o

the
(
Kirc
n
d
0
of

ject
E
w
)
the

y
1
for
+
assuming
i
scattered
k
the
r
a

this
(4)
sible
+
it
k
of
(
that
n
asp
0
in

w
B
cedure
)
the
k
w
(
When
n
pro
0
e

where
E
rec
)
inciden

the
1
h
+
the
i
an
k
the
r
an

h

the
dx
ev
0
in
F
the
or
dissipated
the
Because
Helmholtz
v
equation
small
(1)
to
to
of
b
surface
e
building
satisd
divided
inside
to
V
S
,
shad
the
one
medium
In
con
ed
tained
ld
in
S
V
zero
m
is
ust
assume
b
scattered
e
s
homogeneous
opp
Therefore
inciden
when
addition
w
ose
e
external
need
either
to
a
ev
ature
aluate
ect
the
a
electromagnetic
a
scattering
illumi
from
resnel
ob
for
jects
to
whic
scattered
h
i
ha
cedure
v
summarized
e
s
a
lds
dielectric
B
constan
y
t
are

the
6
the
=
the

s
0
on
,
xi
w
ed
e
e
can
lds
iden
V
tify
using
V
in
with
This
the
es
free
d
space
scattered
b
opaque
et
the
w
the
een
are
the
can
scatterers
o
The
cedure
b
a
oundary
shado
surface
inciden
S
b
V
b
of
t
V
are
is
b
tak
When
en
dealing
equal
n
to
of
S
ximation
[
a
S
in
1
b
where
duces
S
b
is
teractions
the
of
surface
eac
of
teract
the
will
scattering
in
ob
4.1.
jects
pap
and
will
S
st
1
the
is
and
the
as
surface
Optics
t
cedure
init
one
y
with
It
building
can
dure
b
can
e
using
pro
tec
v
h
en
to
ac
on
kson
:
1975,
ld
pages
y
433-434)
is
that
on
the
surface
in

tegral
iteration
o
is
v

er
second
S
of
1
are
v
and
anishes
faces
so
v
the
illuminated
in
iteration
tegral
tribution
in
b
(4)
Kirc
can
not
b
from
e
surface
computed
largely
only
inside
o
building
v
the
er
a
S
e
.
is
It
with
is
ect
useful
the
to
dimensions
sp
the
ecialize
the
(4)
S
to
the
a
can
scattering
e
situation
appro
and
in
to
an
exhibit
region
a
i
formal
a
expression
o
for
ed
the
S
scattering
.
am
the
plitude
w
as
region
an
total
in
on
tegral
surface
of
s
the
nearly
scattered
he
lds
ject
o
opaque
v
e
er
that
S
the
.
ld
Since
S
the
is
scattering
and
ob
osite
jects
the
are
t
supp
In
osed
w
to
supp
b
that
e
building
outside
w
the
are
v
t
olume
with
V
radius
,
curv
equation
large
(4)
resp
holds
to
for
w
the
v
scattered
As
lds
result
(
the
E
nated
s
F
;
co
B
ts
s
rection
),
used
that
ev
is
the
the
ld
total
S
lds
.
(
pro
E
can
;
e
B
as
)
ws
min
the
us
t
the
E
in
and
ciden
i
t
b
lds
the
(
tenna
E
ev
i
on
;
surface
B
of
i
scat
)
s
(i.e.,
scattered
the
E
lds
and
generated
s
b
S
y
appro
ac
mated
tiv
describ
e
ab
sources
v
suc
s
h
scattered
as
throughout
an
space
tennas
are
In
puted
equation
the
(4)
hho
E
tegral
(
d
x
metho
)
giv
dep
a
ends
o
explicitly
estimate
on
the
the
ld
outgoing
an
direction
ob
of
If
k
diaction
.
on
The
building
dep
edges
endence
negligible
on
e
the
simplify
inciden
ab
t
v
direction
pro
is
b
implicit
assuming
in
eometrical
the
s
scattered
w
lds
the
E
t
and
i
B
y
.
that
3
oth
Description
inciden
of
and
the
lds
mo
equal
del
zero
The
eside
ev
obstacle
aluation
w
of
are
the
with
ld
large
E
um
through
er
equation
scatterers
(4),
appro
re
can
quires
w
the
sen
kno
reduction
wledge
the
of
time
the
ecause
lds
re
E
the
s
um
and
er
B
pairn
s
nly
on
pairs
the
scatterers
sur
ee
face
h
S
in
.
This
3.1
ect
A
b
single
discussed
scatterer
detail
In
section
the
Later
absence
this
of
er
kno
e
wledge
refer
ab
the
out
pro
the
as
correct
eneral
lds
cedure
E
to
s
second
and
the
B
eometrical
s
Shado
on
pro
the
3.2
surface
than
S
scatterer
,
dealing
w
more
e
one
m
the
ust
ce
mak
for
e
scatterer
some
b
appro
extended
xi
an
mations
e
First
hnique
w
eac
e
iteration
neglect
onds
the
a
ld
tion
transmitted
a
through
a
the
The
building
t
This
pro
is
b
equiv
all
alen
sources
t
ev
to
uated
assuming
eac
that
building
the
S
build
:
ing
at
is
st
an
the
paque
source
ob
the
ject
tenna
This
starting
assumption
the
is
iteration
reasonable
sources
at
inciden
the
ld
frequencies
the
w
tenna
e
all
are
sur
dealing
whic
with
ha
as
e
the
een
part
during
of
previous
the
The
w
con
a
is
v
aluated
e
y
energy
the
whic
hho
h
tegral
issection
120.0 z=0
building 1 building 2
100.0
receiver
80.0
source
60.0
direct path
path involving one reflection 40.0
axispath involving two reflections

k
20.0
E
0.0
0.0 20.0 40.0 60.0 80.0 100.0 120.0
;
decimetric
onen
for
o
h
e
approac
ulations
tegral
e
Figure
cated
1:
k
Example
n
of
at
sourceeceiver
w
paths
while
involving
is
up
in
to
A
t
i
w
er
o
n
rections
t
These
external
paths
The
a
in
re
et
accounted
free
fo
w
r
ottom
b
along
y
Rap
the
buildings
metho
ergence
d
n
after
(
t
;
w
;
o
3
iterations
refer
b
or
:
the
When
t
the
rep
inciden
computed
t
t
lds
t
E
an
i
a
and
in
B
an
i
w
are
eac
kno
<
wn
sim
at
p
eac
t
h
wn
p
a
oin
the
t
the
of
re
the
,
surface
d
S
follo
n
max
,
E
the
;
scattered
;
lds
E
E
;
s
N
and
of
B
v
s
h
on
o
S
are
n
terial
are
In
ev
is
aluated
part
as
ed
describ
ure
ed
alues
in
of
section
of
3.1.
of
A
m
t
of
the
This
frequencies
b
considered
tric
the
an
diaction
h
lds
metallic
are
ca
small
electric
with
a
resp
e
ect
that
to
eep
the
further
Geometrical
(i.e.
Op
outside
tics
the
lds
in
and
measured
they
ts
can
a
b
1.8
e
axis
neglected
the
if
e
the
C
in
cated
tensit
line
y
w
of
w
inciden
dimensions
t
in
radiation
iteration
is
con
lo
the
w
c
Therefore
through
w
parameters
e
i
ha
E
v
k
e
1
c
)
ho
=
sen
r
to
=
ev
er
aluate
max
the
E
scattered
=
lds
g
on
is
S
b
n
ers
b
orts
y
of
using
at
the
The
eneral
t
pro
t
cedure
the
when
either
the
with
source
dielectric
of
=
the
second
inci
v
den
b
t
h
lds
the
E
is
i
y
and
alls
B
w
i
the
is
t
the
onen
an
electromagnetic
tenna
the
and
lic
the
v
OS
comp
pro
y
cedure
h
elsewhere
than
E
comp
s
x
and
b
B
ely
s
noting
will
irected
b
can
e
ciated
used
irected
in
e
the
rected
st
een
step
alls
of
resonat
the
y
next
x
iteration
related
to
irected
compute
e
the
lea
inciden
system
t
exp
ld
y
on
will
the
hanging
other
iteration
surfaces
b
c
the
:
<
When
more
the
and
electromagnetic
b
lds
of
do
sho
not
4
signian
are
tly
at
c
oin
hange
lo
from
at
one
heigh
iteration
of
to
meters
the
the
F
sho
in
in
3
ure
areas
propagation
Kirc
v
Pisani
w
urban
2:
hho
used
A
test
|
radio
of
Vittoli
metho
etti
Receivers
L
re
next
Figure
one
Geometry
the
to
pro
convergence
cedure
the
ends
d
d
a
:
lo
A
along
t
dashed
the
b
end
t
of
een
the
t
iterativ
o
e
ll
pro
a
cedure
given
w
meters
e
t
kno
n
w
the
the
v
electromagnetic
of
lds
metho
on
is
all
hec
the
ed
in
the
v
wing
olv
er
ed
E
surfaces
=
and
fj
w
n
e
(
can
)
use
n
once
i
again
k
the
j
Kirc
k
hho
1
in
N
tegral
g
to
i
ev
1
aluate
3
the
n
electromagnetic
=
ld
f
throughout
n
the
i
domain
i
V
1
.
3
The
where
ab
r
o
the
v
um
e
er
pro
receiv
cedure
Figure
is
rep
sk
the
etc
alues
hed
er
in
E
ure
eac
1
iteration
in
lines
the
to
case
w
of
diren
t
sim
w
where
o
buildings
buildings
considered
and
metallic
up
built
to
ma
the
of
second
constan
iteration

It
4.
is
the
seen
case
that
con
eac
ergence
h
faster
iteration
ecause
is
eac
equiv
iteration
alen
of
t
inciden
to
energy
include
absorb
new
b
paths
the
connecting
w
the
In
source
4
to
e
the
ort
receiv
v
er
of
T
w
o
comp
analyze
ts
the
the
con
ld
v
in
ergence
case
of
metal
the
buildings
metho
con
d
ergence
w
the
e
onen
consider
E
the
is
conuration
uc
describ
slo
ed
er
in
that
ure
the
2.
onen
A
E
plane
.
w
can
a
e
v
tuitiv
e
explained
of
y
w
that
a
y
v
elec
elength
ld

b
=
asso
6
to
:
x
28
w
meters
v
and
whic
w
is
a
b
v
w
e
the
v
w
ector
as
k
a
=
ing
(0
vit
:
,
89
the
;
irected
0
ld
:
to
42
y
;
w
0
v
:
is
15),
to
whose
v
pro
the
jection
So
on
e
to
ect
z
the
=
comp
0
t
is
k
dra
c
wn
at
in
h
the
nd
picture
rection
has
et
b
een
een
buildings
tak
65
en
y
as
100)
the
than
inciden
them
t
this
ld
conmed
A
y
t
results
eac
the
h
ulation
iteration
wn
the
ure
ld
(b
v
aluesEx
2.5
1 iter.
2 iter.
non metallic buildings 10 iter.
2.0 metallic buildings
0.9
1.5
1.0
-0.1
0.5
0.0 -1.1
1 2 3 4 5 6 7 8 9 10 40.0 60.0 80.0 100.0 120.0
Number of iterations y
Ey
3.0
1 iter.
2 iter.
3 iter.
10 iter.2.0
1.0
0.0
-1.0
-2.0
40.0 60.0 80.0 100.0 120.0
y
(
Convergence
0
of

the
x
metho
r
d
0
White
v
dots
e
refer
equiv
to
0
p
e
erfectly
j
conducting
h
surfaces
e
Black
regions
dots
.
refer
(
to
(4)
dielectric
r
surfaces
+
(
with

a
=
r
4
r
).
0
The
|
quantit
that
y
0
aximum
d
va
o
riation
phase
rep
e
o
8
rted
r
on
far
y
=
axis
F
is
2
deed
0
as
x
er
r
n
x
E
using
=
=
max
0
f
(5)
er
00
n
0
E
2
i
:
;
hho
i
Figure
=
x
1
w
;
=
3

g
Let
,
the
where
S
er
titativ
n
the
E
b
i
at
is
T
er
d
n

E
d
i
(8)
=
the
max
this
fj
s
E
e
n
k
i
)
(

k
2
)
(
E
00
n
x
1
j
i
2
(
large
k
0
)
,
j
rewrite
;
suitable
k
expansion
=
0
1
(
;
00
N
+
r
T
g
0
;
2
i
[
=
(
1
00
;
r
3
:
and
decimetric
N
tegral
r
A
is
C
the
a
numb
on
er
Pro
of
is
receivers
ha
4
ik
Computational
ik
asp
i
ects
x
In
)
this
0
section
b
w
of
e
1
presen
.
t
a
the
estimate
main
ose
computational
on
fea
b
tures

of
lo
the
(6),
mo
see
del
negligible
4.1
0
Scaling
)
of
0
the
=
computer
0
time
)
with
0
the
h
dimension
t
of
ellno
the
relation
domain
ximation
The
ecomes
mo
x
del
4
describ
r
ed
e
in
x
section
x
3
and
requires
=
on
r
eac
0
h
r
sur

face
x
the
x
ev
)
aluation
j
of
0
the
00
ld
2
scattered
1
from
When
the
is
other
compared
surfaces
x
Th
and
us
00
at
w
st
can
sigh
r
t
a
the
T
computational
ylor
time
r
gro
r
ws
r
as

N
x
2
x
with
)
N
0
the
T
n
where
um
=
b
x
er
x
of
j
surface
2
elemen
0
ts
r
radio

urban
x
4:
x
onents
)]
x
2
and
3
y
+
ottom
:
dif
(6)
iterations
for
the
approac
sho
in
in
Kirc
2.
k
spatial
Vittoli
areas
etti
3:
w
(
(9)
Pisani
o
e
propagation
in
4
o
Rap
along
F
dashed
L
is
and
Figure
as
Comp
D
E
4
op
with
E
D
(b
the
at
diameter
ferent
of
in
the
example
domain
wn
This
ure
b
The
eha
co
vior
rdinate
is
the
ev
line
en
rep
w
rted
orse
the
than
axis
the
vided
one
T
sho
negligible
wn
e
b
v
y
e
the
r
metho
e
ds
r
whic
e
h
k
require
(
the
0
discretization
00
of
(7)
the
d
whole
and
space
00
(
e
D
diameters
3
the
),
S
but
and
as
2
w
T
e
get
will
quan
sho
e
w
w
can
imp
b
that
e
error
easily
the
impro
is
v
ounded
ed
y
Consider
8
the
Then
situation
oking
describ
equation
ed
w
in
can
ure
that
5.
is
S
when
2
d
is
+
a
00
region
2
of
r
the
<
surface
16
on
)
whic
d
h
+
w
00
e
<
need
r
to
2
compute
whic
the
is
ld
alen
scattered
to
b
w
y
wn
a
ld
region
In
S
appro
1
equation
of
b
another
E
surface
(
w
)
e
i
can

write
ik
r
0
=
0
r
i
0

x
00
0
+
00
maximum variation
Intensity Intensityo’ r0
x’ o’’ k
Kirchhoff
x’’S1 N^1.5r
N
S2
21.0
11.0
1.0
4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36
N
d
etti
and
Rap
the
F
hho
Figure
the
5:
giv
Scattering
to
geometry
Optics
F
b
(
y
k
are
)
tribution
=
b
I
propagation
S
e
1
w
e
refer
i
t
k
no

discon
x
it
0
ground
f
tions
k
b
(
tenna
n
small
0
In

treatmen
B
oin
s
ground
)
5
+
times
[
from
k
The

ligh
(
only
n
inciden
0
all

but
E
allround
s
t
)
alues
k
case
(
t
n
t
0
to

t
E
ts
s
The
)]
a

of
1
er
+
Rections
i
of
k
areas
r
ground

from
dx
exp
0
the
(10)
computed
In
|
equation
propagation
(10),
rre
the
the
in
e
tegral
w
do
for
es
whic
not
noticeable
dep
oundaries
end
These
on
deal
the
or
co
lds
ordinate
F
x
discon
00
domain
in
ensated
ternal
due
to
the
the
90
region
tin
S
seen
2
w
.
notice
Therefore
conducting
when
diaction
w
t
e
dis
decomp
onding
ose
alls
the
tion
surfaces
attributed
in
is
to
ure
regions
represen
whic
resp
h
e
are
p
small
the
compared
y
to
scattering
their
um
relativ
buildings
e
further
distance
the
w
sim
e
radio
need
ur
to
n
ev
of
aluate
the
the
tegral
in
computational
tegral
of
only
e
once
therefore
for
from
eac
b
h
y
region
Vittoli
Since
h
the
urban
area
v
of
decimetric
the
ond
surface
fo
regions
integral
is
cannot
prop
for
ortional
build
to
W
r
ure
0
brief
ee
of
equation
are
(8)),
GO
with
near
this
w
appro
the
ximation
rected
the
oundaries
computa
w
tional
rections
time
all
for
rections
the
and
ev
v
aluation
tin
of
the
the
rection
in
u
teraction
inside
b
ure
et
partially
w
y
een
uit
t
the
w
rection
o
all
surfaces
form
is
angle
in
o
v
coincide
ersely
b
prop
the
ortional
the
to
con
their
W
relativ
in
e
p
distance
alls
r
ects
0
nonex
.
fact
Moreo
o
v
only
er
con
the
o
probabilit
single
y
the
that
corresp
one
double
build
h
ing
b
A
a
1
source
an
ted
see
GO
another
p
building
and
A
a
2
receiv
(i.e.
ely
there
can
aren
the
buildings
an
b
t
et
surface
w
er
een
either
A
oin
1
diren
and
a
A
n
2
b
)
of
decreases
and
exp
N
onen
4.2
tially
on
with
ground
r
the
0
ulation
.
the
Then
electric
it
in
can
ban
b
a
e
umerical
demonstrated
t
that
the
the
using
compu
Kirc
tational
in
time
is
gro
a
ws
p
as
t
N
view
3
ensiv
2
W
with
assume
the
that
n
rection
um
the
Pisani
can
tegral
e
6:
b
in
Geometrical
hho
This
C
approac
=T
A
fo
Kirc
a
areas
r
in
of
e
N
a
p
radio
buildings
diamonds
for
L
of
t
relative
T
b
Figure
er
Scaling
N
the
of
time
trans
N
paren
4
t
r
buildings
regula
and
distribution
linearly
p
with

the
N
n
Black
um
co
b
sp
er
to
of
measured
opaque
r
ones
Kirchho
In
metho
order
assumption
to
b
c
done
hec
rection
k
the
these
ings
results
alls
w
e
e
to
ha
7
v
a
e
explanation
p
ects
erformed
diaction
some
h
computations
neglected
on
y
conurations
are
of
only
n
the
2
thado
buildings
b
with
of
n
inciden
=
and
2
ld
;
b
:
exist
:
when
:
e
;
with
6
from
illuminated
w
b
F
y
ground
the
the
ld
t
emitted
rected
from
ha
an
e
an
discon
tenna
uities
The
or
buildings
ground
ha
double
v
the
e
tin
b
it
een
is
uniformly
the
distrib
ee
uted
8),
o
is
v
comp
er
b
a
the
square
tin
region
y
(
to
n
w
p
double
er
Since
direction
w
and
and
the
ground
an
a
tenna
degree
has
the
b
w
een
discon
lo
uities
cated
and
at
can
the
e
cen
that
ter
v
of
of
the
t
la
o
y
tributions
out
similar
The
e
results
that
are
the
sho
of
wn
erfectly
in
w
ure
and
6
the
where
of
N
are
=
isten
n
In
2
the
represen
w
ts
surfaces
the
e
n
three
um
tinct
b
tributions
er
w
of
corresp
buildings
to
and
rec
the
on
relativ
w
e
one
computational
onding
time
a
is
rec
deed
Eac
as
con
the
can
ratio
e
T
to
N
diren
=T
image
4
and
b
accoun
et
for
w
y
een
In
the
9,
compu
oin
tational
A
time
B
required
t
for
source
N
a
buildings
er
and
ectiv
that
.
required
source
for
b
4
either
buildings
an
W
or
e
y
can
oin
see
of
that
scattering
the
while
scaling
receiv
of
can
time
e
is
an
similar
p
to
t
N
a
3
2
for
Relative computational timeIncident + reflected
field
Building
Incident
field
Shadow
Incident field


1.0
Incident fieldGround-wall
reflection
0.0
Wall-ground
reflection
explicit floor
implicit floor
-1.0
40.0 60.0 80.0 100.0 120.0
y
A r
B
r
1
r
2
C ground
=
some
oundary
urban
at
with
b
urban
are
mem
explicitly
(
describ
to
ulation
of
rections
decimetric
one
for
in
h
v
approac
decrease
tegral
obtained
in
rep
hho
In
Figure
and
7:
approac
Lighthado
ts
w
5.1
b
here
ounda
b
ries
white
of
section
the
treatmen
incident
umerical
and
obtained
re
and
cted
a
ld
computational
v
at
diaction
when
a
itly
gral
obtained
ground
v
the
of
v
in
erimen
w
results
v
Kirc
W
A
a
|
of
Vittoli
tenna
C
In
Figure
represen
8:
y
Geometrical
with
description
The
of
of
the
to
ground
whic
all
to
dou
an
ble
ure
rection
allo
paths
same
surface
reduc
or
oth
an
and
y
.
p
it
oin
the
t
ground
where
of
w
whic
e
computed
w
rections
an
e
t
radio
to
y
compute
ground
the
Kirchho
ld
dashed
v
the
alue
y
GO
5
3
4,
e
presen
related
v
complexit
6
metho
areas
here
pap
of
b
describ
sim
with
is
Figure
comparing
9:
lation
Geometrical
n
scheme
of
to
measure
illustrate
this
ground
are
rec
t
tions
these
treatment
electric
F
in
or
vironmen
eac
giv
h
results
direct
umerical
path
the
r
ld
from
an
A
out
to
vironmen
B
11,
there
c
is
the
another
real
path
parameters
r
3,
1

and
4.
r
implicit
2
t
whic
ground
h
brings
in
n
v
results
olv
h
es
similar
one
those
rection
with
on
explicit
the
ee
ground
10)
at
it
the
ws
p
the
oin
time
t
sensible
C
tion
The
b
con
the
tribution
time
of
storage
the
ory
rec
Moreo
tion
er
in
a
C
oids
to
generation
the
the
ld
b
v
,
alue
titious
in
terms
etti
h
urban
the
10:
ld
onent
ground
x
are
in
considered
propagation
v
F
w
Pisani
the
the
Rap
ld
their
continuous
L
rep
B
Figure
is
Comp
computed
E
considering
of
the
electromagnetic
lengths
The
r
line
1
resents
,
ld
r
b
2
explic
and
computing
b
rections
y
the
using
inte
F
the
resnel
line
co
resents
eien
ld
ts
b
These
implicit
considerations
purely
ha
)
v
treatment
e
Results
b
sections
een
and
v
w
alidated
ha
b
e
y
ted
imple
results
men
to
ting
con
b
ergence
oth
time
the
y
explicit
the
v
d
aluation
ed
of
The
the
alidation
Kirc
the
hho
h
in
ed
tegral
this
and
er
implicit
done
eometrical
y
Optics
the
treatmen
u
ts
outcome
of
a
rections
um
on
er
the
exp
ground
tal
and
men
comparing
In
the
section
results
e
on
going
a
presen
test
the
conuration
of
The
comparisons
la
Radio
y
co
out
erage
(p
an
osition
en
of
t
buildings
e
and
e
observ
the
ation
of
axis
n
is
sim
sc
of
hematically
propagation
describ
the
ed
emitted
in
y
ure
an
2.
through
The
an
input
en
data
t
of
ure
the
the
sim
blo
ulation
ks
are
t
equal
buildings
to
a
those
cit
giv
and
en
in
Intensitythe
L
and
Pisani
the
F
and
Rap
are
etti
g
C
one
Vittoli
metho
|
where
A
in
Kirc
Optics
hho
obtained
in
f
tegral
the
approac
tenna
h
ld
for
Geometrical
decimetric
Kirc
radio
oth
w
using
a
It
v
er
e
Kirc
propagation
the
in
hhof
urban
and
areas
the
7
large
ielectric
a
constan
v
t
to

ligh
,
ab
heigh
d
t
are
h
the
b
in
,
,
etc
ere
ha
12,
v
threshold
e
v
b
with
een
ts
c
measured
hosen
v
uniformly
-70
o
1
v
ery
er
only
the
under
la
it
y
co
out
due
b
and
ecause
The
of
uni
the
v
lac
buildings
k
cedure
of
that
detailed
prosp
measuremen
een
ts
dees
In
o
particular
threshold
w
Kirc
e
y
ha
w
v
e
e
d
c
w
hosen
threshold

input
=
parame
4
equal
and
erimen
h
y
b
bionics
=
v
10
represen
meters
er
The
seen
frequency
aluated
of
20%
the
metho
ld
of
emitted
f
b
it
y
those
the
d
source
13
n
the
tenna
v
with
p
a
3
nonisotropic
are
radiation
b
pat
v
tern
concerns
is
measured
1.8
b
GHz
calcula
whic
region
h
's
corresp
clearly
onds
on
to
of
a
the
w
to
a
y
v
t
elength
an
of
ork
16.7
the
cen
around
timeters
y
Input
er
p
an
o
Kirc
w
v
er
is
is
problem
250
elopmen
mW
e
The
en
source
gra
is
area
lo
puted
cated
er
at
v
a
b
heigh
tics
t
's
of
dark
4.5
the
meters
p
and
v
the
o
electric
only
ld
's
is
area
computed
p
at
v
a
elo
heigh
represen
t
k
of
for
1.6
ulations
meters
source
The
diameter
cross
those
in
5.1.
the
data
cen
tained
ter
DECT
of
tester
the
In
ure
computed
marks
area
the
tenna
p
for
osition
o
of
-85
the
b
source
the
and
area
the
Geo
gra
ab
y
than
scale
hho
indicates
Moreo
the
w
sim
p
ulated
ered
ld
;
p
lie
o
all
w
alues
er
aluated
decreasing
's
from
ab
white
threshold
to
b
blac
y
k
up
F
The
or
fall
an
in
area
ts
of
2
ab
4
out
of
300
or
meters
to
diameter
of
the
's
whole
area
sim
t
ulation
oin
has
g
tak
alue
en
threshold
24
Geometrical
hours
hho
on
lo
a
the
RISC
on
560
Kirc
IBM
of
Figure
erage
11:
v
Example
Optics
of
other
the
the
p
input
o
y
w
is
er
is
emitted
a
b
of
y
6
an
p
anisotropic
the
antenna
tenance
obtained
mobile
b
net
y
based
the
diction
Kirchho
w
integral
atten
metho
source
d
the
5.2
of
Comparison
this
with
e
exp
prop
erim
e
en
on
tal
in
mea
Results
suremen
sho
ts
metho
Comparisons
ell
with
e
exp
ors
erimen
of
tal
in
measuremen
v
ts
b
in
tak
realis
The
tic
t
conditions
y
are
the
quite
where
diult
com
b
p
ecause
w
of
is
the
o
absence
e
of
for
accurate
oth
data
Op
on
and
the
hho
computational
metho
domain
The
The
gra
dielec
indicates
tric
area
constan
the
t
o
of
er
the
alues
in
ab
v
v
olv
threshold
ed
for
materials
hho
and
metho
the
The
heigh
where
t
computed
of
o
the
er
buildings
alues
are
b
not
w
kno
is
wn
ted
Moreo
blac
v
The
er
data
the
b
map
sim
of
uildings
the
ters
to
frequency
wn
area
is
are
a
to
v
used
ailable
section
with
Exp
a
tal
precision
w
of
ob
one
b
meter
a
therefore
system
the
propagation
direction
Sym
error
SP
for
ure
shorter
the
w
co
alls
erage
can
of
b
an
e
is
of
ted
the
a
order
p
of
w
ab
of
out
dbm
ten
can
degrees
e
whic
that
h
co
doubles
erage
when
ev
w
with
e
metrical
consider
is
the
out
direction
less
of
obtained
the
Kirc
rected
's
ld
d
F
v
urther
t
due
o
to
the
in
oin
terfer
consid
ence
(
phenomena
5
the
7
ld
)
can
outside
v
while
ary
the
considerably
v
within
and
one
ev
meter
with
Finally
hho
,
metho
no
are
information
o
is
e
a
Figure
v
has
ailable
een
ab
b
out
raising
the
threshold
presence
to
of
dbm
ob
measured
jects
alues
other
under
than
threshold
buildings
four
a
oin
park
f
ed
;
car
;
can
;
easily
g
in
three
tercept
them
the
outside
groundected
v
ld
close
In
the
these
oundaries
conditions
Kirc
the
f
computed
co
p
erage
o
The
w
eviden
er
discrepancy
ma
p
y
t
easily
4
dir
whose
b
v
y
is
10
the
to
while
20
oth
dbm
Optics
from
Kirc
the
's
measured
tions
one
cate
and
inside
p
illuminated
oin
Therefore
t
the
to
hand
oin
hho
t
estimate
comparisons
the
cannot
v
b
area
e
impro
used
es
to
Geometrical
v
one
alidate
the
or
hand
in
to
v
inaccuracy
alidate
the
the
data
co
reliabilit
de
of
On
results
the
limited
other
it
hand
advisable
the
use
errors
fairly
are
margin
fairly
safet
lo
.
calized
Conclusions
and
starting
the
oin
co
for
de
planning
can
main
b
of
e
y
still
comm
used
cation
with
w
some
is
conence
on
to
pre
estimate
of
the
radio
co
a
v
e
erage
uation
areas
the
In
n
the
in
follo
la
wing
out
ures
the
the
In
cross
pap
marks
w
the
ha
p
e
osition
osed
of
iterativ
the
pro
an
based
tenna
the
while
hho
the
tegral
n
d
um
ha
b
e
ered
wn
dots
this
indicate
d
the
w
p
suited
ositions
solv
of
this
the
and
p
promising
oin
ects
ts
dev
where
t
measuremen
this
ts
The
ha1 1
4 48 8
2 2
14 149 9
6 6
5 5
13 13
12 15 12 15
16 16
11 117 73 3
17 17
10 10
Pro
urban
Areas
in
L
areas
t
propagation
Y
e
urner
v
Micro
a
Ph
w
in
radio
rans
decimetric
Classical
Figure
terrainased
12:
O
The
Notes
coverage
hho
region
tracing
of
A
the
egami
antenna
strength
fo
in
r
pp
a
New
thresh
k
old
an
p
7,
o
W
w
Approac
er
pp
of
ev
-85
zero
dbm
J
The
o
light
The
gra
for
y
ol
a
euc
rea
prediction
dees
bile
the
,
zone
|
where
12
b
J
oth
John
Geometrical

Optics
and
and
results
Kirchho
e
's
IEEE
metho
v
d
1993.
p
J
redict
Prop
a
ely
value
TML
fo
,
r
L
the
the
ld
scattering
p
metho
o
h
w
.
er
Gardiol
ab
prediction
ove
viron
threshold
C
The

da
ra
rk
a
gra
Radio
y
pp
a
T
rea
and
co
heoret
rresp
mean
onds
urban
to
IEEE
the
tenn
zone
approac
where
Kirc
8
etti
only
F
the
as
value
ol
p
1991.
redicted
Jac
b
dynamics
y
&
Kirchho
ork
's
K
metho
Cic
d
Wiesb
is
cepts
ab
3D
ove
a
threshold
mo
The
v
black
Se
a
un
rea
11,
dees
1002-
the
P
zone
B
where
agen
b
v
oth
in
Geometrical
a
Op
aral
tics
Using
and
d
Kirchho
Computer
's
ol
metho
1996.
d
\
p
reduction
redict
umer
values
of
b
the
elo
in
w
to
threshold
sional
The
to
antenna
Comm
lo
K
cation
agen
is
Ra
indicated
path
b
t
y
cellular
the
ts
cross
IEEE
while
pp
the
Netherlands
numb
.
ered
,
dots
mo
indicate
radio
the
e
lo
urban
cation
,
of
no
receivers
1992.
fo
T
r
ak
exp
hi
er
S
imental
oshida
measurements
ical
mo
of
for
ld
Figure
for
The
mo
region
radio
the
T
fo
An
a
Propagat
old
h
o
tegral
er
hho
-70
A
The
Vittoli
of
C
symb
Rap
and
Pisani
scale
the
del
13:
is
coverage
curren
of
tly
antenna
under
r
revision
thresh
to
p
impro
w
v
of
e
dbm
its
meaning
time
the
requiremen
ols
ts
grey
isani
is
1997).
same
7
ure
Ac
v
kno
39,
wledgemen
299-302,
ts

W
D
e
kson
w
Electro
ould
,
lik
Wiley
e
Sons
to
Y
thank
(1975).
Prof
T
G

Mazzarella
D
of
hon
the
W
Uni
ec
v
on
ersit
and
y
for
of
digital
Cagliari
w
for
v
helpful
propagation
discussions
dels
and
o
commen
erview
ts
J
on
lect
the
Comm
man
,
uscript
ol
W
no
e
pp
ac
12,
kno

wledge
.
A
Luthi
Lo
Chopard
c
W
he
\
for
a
pro
e
vid
agation
ing
Urban
us
cells
with
Massiv
a
P
graphic
lel
in
h
terface
the
for
Metho
the
Lecture
visualization
in
of
Science
the
v
radio
1041,
electric
408,
co

v
Pisani
erage
Dimensionali
and
y
Ing
in
F
n
T
ical
allone
aluation
of
electroagnetic
the
from
Cen
Kirc
tro
's
Studi
tegral
e
d
Lab
a
oratori
dimen
T
approac
elecom
submitted
unicazioni
Computer
SEL
ys
T
un
of

T
Rizk
urin
W
for
F
pro
\
viding
y
us
based
with
loss
exp
in
erimen
w
tal
micro
measuremen
en
ts
men
This
in
w
ceedings
ork
PIMR
has
,
b
210-4,
een
Hague
carried
1994.
out
J
with
Rossi
the
Levy
ancial
\A
supp
y
ort
del
of
decimetric
T
w
elecom
v
Italia
propagation
and
an
the
area
Sardinia
Science
Region
v
Authorities
27,
References
6,

971-9,
F
Ik