Encuentro de Investigacin en Ingeniera Elctrica

Zacatecas, Zac, Abril 5 7, 2006

Cmputo de Radioenlaces Terrestres en Java

Andrs Calvillo Tllez, Carelia Gaxiola
Centro de Investigacin y Desarrollo de Tecnologa Digital, IPN,
Tijuana BC CP 22510.
TEL: +(664)6231344, ext. 82837, correo-e: [email protected].

Resumen Se presenta el modelado de radioenlaces

terrestres en java, los clculos se realizan sobre un
enlace de microondas que opera en la frecuencia 23
GHz. Se realiza primero el clculo de las prdidas de
propagacin en el espacio libre, posteriormente se
analiza la atenuacin por desvanecimiento, seguido del
clculo de la primera zona de Fresnel, para garantizar
que el nivel de la seal se encuentra dentro del rango
adecuado. Se propone el tipo de antena, cables y
conectores, en base a esto se realizan los clculos.
Abstract We are presenting a terrestrial radiolink
model in Java, it was designed to operate at the 23 GHz
band. We are showing the computations to obtain the
losses in free space. To guarantee that the signal level
stays in an adecuate range, we analized the attenuation
due to scattering, and the 60% clearance of the Fresnel
zone. The analisys was made considering a specific kind
of anntena cables and connectors
Palabras clave Lnea de vista, Microondas, Nivel de
recepcin de la seal, Potencia efectivamente radiada,
Potencia isotrpica radiada efectivamente, Prdidas por
propagacin, Radioenlaces.

I. INTRODUCCIN

n radioenlace terrestre establece la comunicacin

bidireccional entre dos transreceptores ubicados en
dos puntos fijos con coordenadas geogrficas
especficas mediante la emisin de radiofrecuencias de
microondas [1]. Las estaciones de microondas constan de
un par de antenas con lnea de vista que se conectan a un
radio transmisor cuya frecuencia de radiacin se encuentra

en el rango de 1 GHz a 58 GHz. Las frecuencias utilizadas

comnmente estn entre los 18 y 26 GHz. La nica
limitante de estos enlaces es la curvatura de la tierra, aunque
con el uso de repetidores puede extender su cobertura a
miles de kilmetros [2].
El trmino lnea de vista (LOS Line of Sight) no significa
una lnea de transmisin en lnea recta, ya que la estructura
de la atmsfera terrestre es tal que su ndice de refraccin
vara con la altura, lo que causa un efecto de curvatura en la
trayectoria de transmisin. El trmino LOS se refiere a que
la onda es captada por el receptor directamente desde el
transmisor. Este tipo de onda solo se recibe dentro del
horizonte visual, como se ilustra en la Figura 1. En un
enlace debe existir lnea de vista, esto es que no existan
obstrucciones y se debe evitar el efecto del factor de
curvatura de la tierra (k), as como obstrucciones que
permitan una claridad de al menos el 60% del libramiento
de la primera zona de Fresnel. Esta es la regin que encierra
el primer elipsoide que rodea la trayectoria entre el
transmisor y el receptor, la cual contiene la mayor cantidad
de potencia destinada al receptor [3]. El anlisis de un
radioenlace terrestre relaciona algebraicamente la potencia
del transmisor, con las ganancias, prdidas, y efectos de
propagacin de lnea de vista como lo son el porcentaje de
desvanecimiento por trayectorias mltiples y la atenuacin
por lluvia. El cmputo de radioenlaces terrestres utiliza
estos mtodos de cuantificacin de esos parmetros para
que en base a una evaluacin se determine que variables
debern ser modificadas para optimizar el diseo.

32

Encuentro de Investigacin en IE, 5 7 de Abril, 2006

EIRP

Radio de Fresnel

Cobertura

FSL

GTX

GRX

La siguiente expresin nos permite calcular la distancia

entre los dos puntos,

Claridad en la
Linea de Vista

d = P 2 + ( Latitud) 2

(3)

Donde:

PTX
Transmisor

PRX
Receptor

Alcance

Figura 1. Radioenlace Terrestre.

Latitud = diferencia entre la latitud de el punto A y la

latitud del punto B en km.
d = P 2 + ( Latitud) 2 = 9.5914 km

II. ESTIMACIN DE RADIOENLACES

Para conocer la distancia de separacin, primeramente se
deben de localizar los puntos de radio transmisin dentro de
un mapa topogrfico de la regin como se ilustra en la
figura 2. Posteriormente se utilizaran las formulas
respectivas para calcular la distancia y los ngulos de
elevacin [4].

Para conocer el ngulo de apuntalamiento podemos

utilizar cualquiera de las dos expresiones que se muestran a
continuacin
(4)
= cos1 (Pd )
Latitud
(5)
= arc tan Longitud
Sustituimos los valores obtenidos de P y
ecuacin 4
= cos1(Pd ) = 11.5191 = 1131'8.77"

d en la

Para calcular los azimut sustituimos el valor obtenido del

ngulo de apuntalamiento en la siguiente expresin
Az A = 90
(6)
Az B = 360

AzA = 90 1131'8.77"= 7828'51.23"

AzB = 360 1131'8.77"= 34828'51.2"
La figura 3 muestra un diagrama indicando los azimut y
ngulo de apuntalamiento,

Figura 2. Ubicacin topogrfica del enlace.

A. Clculo del Azimut y distancia [5-15]
Las coordenadas geogrficas para la estimacin del enlace
son las siguientes,
Punto A: 32 27 1.62 N 116 58 6.5 E
Punto B: 32 28 13 N 116 51 11.4 E

Azimut B
Latitud

Angulo de apuntalamiento

Azimut A
Longitud

Sabiendo que 1 equivale a 111 Km, tenemos que

Punto A: (3,134.71, 11,299.15)
Punto B: (3,136.63, 11,288.02)

Figura 3. Azimut y ngulo de apuntalamiento.

Utilizando la ecuacin 1 calculamos la latitud promedio,

LatA + LatB
(1)
Lm = Latitud Promedio =
2

donde
LatA = latitud del punto A en km.
LatB = latitud del punto B en km.

FSL =

3,134.71+ 3,136.63
Lm =
= 3,135.67 = 3227' 37.31"
2

P = Longitud (km)cos Lm

(2)

Donde:
Longitud = diferencia entre la longitud del punto A y la
longitud del punto B en km.
P = Long (km) cos(3227 '37.31")
= Long (km) (0.843762937) = 9.3983 km

B. Clculo de las prdidas en el espacio libre (FSL)

La perdida por espacio libre la podemos obtener al
relacionar la distancia de separacin entre los dos puntos y
la frecuencia de transmisin, esta relacin se muestra en la
siguiente expresin,
4 1012
+ 20log dkm + 20log f GHz
3 10 8

dkm = es la distancia del enlace en km.

fGHz = es la frecuencia de operacin en GHz.

FSL = 92.44 + 20log(9.59) + 20log(23) = 139.31 dB

(7)

33

Encuentro de Investigacin en IE, 5 7 de Abril, 2006

C. Clculo de la potencia isotrpica radiada

efectivamente (PIRE)
Para el clculo de la potencia isotrpica radiada
efectivamente se emplea la siguiente expresin

PIRE = PT LL + GT

F. Clculo del ngulo de apertura del lbulo y ngulo

de media potencia del lbulo
Conociendo la frecuencia de operacin del enlace podemos
calcular la longitud de onda por medio de la siguiente
expresin,

(8)

Donde:
PT = Potencia de transmisin en dB
LL = Perdida por lnea de transmisin en dB, la cual se
calcula de acuerdo a las siguientes suposiciones,

(12)

Donde:

Prdidas de conectores: para los conectores se propone

un a perdida de 0.2 dB cada una, se necesitan cuatro
conectores (2 Tx, 2 Rx) = 0.8 dB
Prdidas por cable: 20 metros de 1/2 superflex, 50 ,
0.20 dB/m = 4 dB
GT = Ganancia de transmisin en dB

c = velocidad de la luz
f =frecuencia de operacin
3 10 8
=
= 13.04 103 m
9
23 10
El ngulo de media potencia del lbulo se calcula
utilizando la siguiente expresin,

PIRE = 17 4.8 + 34 = 46.2 dB

BW =

D. Clculo del nivel de recepcin de seal (RSL)

El clculo del nivel de recepcin se efecta por medio de la
siguiente expresin,

RSL = PIRE FSL + GR LL

c
f

(9)

Donde:

PIRE = Potencia isotrpica radiada efectivamente

FSL = Prdidas en el espacio libre
GR = Ganancia en el receptor
LL = Perdida por lnea de receptor en dB, la cual se calcula
igual que las pedidas por lnea de transmisin.

70
D

(13)

Donde:
= longitud de la onda
D = dimetro del plato de la antena de transmisin

D = 305mm = 305 103 m

70(13.04 103 m)
BW =
= 259' 34.03"
305 103

El ngulo de la apertura del lbulo se obtiene

sustituyendo el valor del ngulo de media potencia en la
siguiente ecuacin
= 2BW
(14)

RSL = 46.2 139.31+ 34 4.8 = 63.91 dB

= 2BW = 2(259' 34.03") = 559'8.07"

E. Clculo del margen de desvanecimiento (FM)

El clculo del margen de desvanecimiento se realiza
mediante la siguiente expresin,

G. Clculo de la altura mnima de las antenas

Para el clculo de la primera zona de Fresnel como se
ilustra en la figurase 4 aplica la siguiente ecuacin,

FM = 30log dkm + 10log(6abf Hz ) 10log[(1 )70] (10)

Donde:
dkm = distancia del enlace en km.
a = Factor topogrfico
b = Factor climtico
f Hz = Frecuencia de operacin en Hz
= Confiabilidad del enlace
Para el enlace los valores de los parmetros son los
siguientes
Sustituyendo los valores en la ecuacin 10
FM = 30 log(9.59) + 10 log(6(1)( 12 )(23 109 )) 10 log[(1 0.9998)70]
= 156.38 dB

F1 = 17.32

d1d2
f GHz dkm

(15)

Donde:

dkm = distancia del enlace en km = 9.59

fGHz = frecuencia de operacin del enlace en GHz = 23
F1 = 17.32

(1.65)(7.94)
= 4.22m
(23)(9.59)

Para obtener la altura mnima de las antenas utilizamos

la siguiente expresin

hTx = hRx = hO + h + 0.6F1

(16)

Donde:
h0 = altura del obstculo mas alto = 200 m
F1 = primera zona de Fresnel = 4.22 m

h=

d1d2
1.5k

(17)

34

Encuentro de Investigacin en IE, 5 7 de Abril, 2006

de la cual

k = factor de curvatura de la tierra = 4 3 (d < 80 mi)

d1 = distancia del obstculo mas alto al punto A = 1.65
d2 = distancia del obstculo mas alto al punto B = 7.94
(1.65)(7.94)
h=
= 6.5505m
1.5(0.25)
hTx = hRx = 200 + 6.55 + 0.6(4.22) = 209.082m
La altura mnima para la antena de transmisin y
recepcin es igual. La figura 4 ilustra la altura de cada
antena, tanto de transmisin como recepcin as como la
altura del radio efectivo de la tierra y el obstculo que tiene
que librarse para lograr la claridad del 60% de la primera
zona de Fresnel.
1a ZONA DE FRESNEL

F1 = 17 .32

d1 d 2
( metros )
fD

LINEA DE VISTA

hTX
0.6F1
hRX

Permisionario CITEDI-IPN
Frecuencia 24 GHz
Horario 24 Horas por da
Circuito Cerro las Abejas Tx) a La Gloria (Rx)
Longitud del enlace 9.59 Km
Potencia:
Nominal:+17 dBm
P.I.R.E.: +49 dBm
Antena:
Tipo: Parablica
Polarizacin: Lnear
Angulo de abertura:
Ancho del haz, Plano H 2 59 34.03
Ancho del haz, Plano E 2 59 34.03
Angulo de elevacin 258 28 51.22
Ganancia 34 dBi min.
Azimut de mxima radiacin 348 28 51.2
Altura S.N.M 320 m
Altura S.T. 10 m
No. Canales 30 de 64 Kbps
Latitud 32 28 13
Longitud 116 51 11.4

RADIO EFECTIVO DE LA TIERRA

NIVEL DEL MAR

d1

ae = ka

ho
D

k=
d2

Figura 4. Primera zona de Fresnel y curvatura terrestre.

III. RESULTADOS
Todas estas deducciones matemticas se integraron al
programa en Java para realizar el cmputo de radioenlaces
terrestres obteniendo lo siguiente:
Permisionario CITEDI-IPN
Frecuencia 24 GHz
Horario 24 Horas por da
Circuito La Gloria (Tx) a
Cerro las Abejas (Rx)
Longitud del enlace 9.59 Km
Potencia:
Nominal: +17 dBm
P.I.R.E.: +49 dBm
Antena:
Tipo: Parablica
Polarizacin: Lnear
Angulo de abertura:
Ancho del haz, Plano H 2 59 34.03
Ancho del haz, Plano E 2 59 34.03
Angulo de elevacin 11 31 8.77
Ganancia 34 dBi min.
Azimut de mxima radiacin 78 28 51.23
Altura S.N.M 280 m
Altura S.T. 10 m
No. Canales 30 de 64 Kbps
Latitud 32 27 1.62
Longitud 116 58 6.5

MODULACION FSK de 4 niveles

ANCHO DE BANDA 3.5 MHz
UMBRAL DEL RECEPTOR -87dBm
Para 10-6 BER
PRDIDAS EN ESPACIO LIBRE 139.31 dB
PRDIDAS POR DESVANECIMIENTO 156.38 dB
POTENCIA DE TRANSMISION 17.0 dBm
PRDIDAS EN LOS CONECTORES 0.8 dB
PRDIDAS EN LA LNEA DE TRANSMISION 4.0 dB
PRDIDAS TOTALES 317.49 dB
NIVEL DE SEAL RECIBIDA -61.11 dB
En la figura 5 se presenta el perfil topogrfico del
radioenlace.

La regin de libramiento de la primera zona de Fresnel

se encuentra a 4.22 m.
La altura mnima de las antenas para librar el 60% de la
primera zona de Fresnel es de 209.082 m

Figura 5. Perfil topogrfico del enlace

35

Encuentro de Investigacin en IE, 5 7 de Abril, 2006

IV. CONCLUSIN
Las antenas disponibles en el mercado permiten efectuar
radioenlaces como el propuesto en este trabajo, la antena
utilizada aqu se selecciono por su frecuencia de operacin,
potencia, ganancia, tamao, modulacin y capacidad de tasa
de bits. Estas caractersticas son fundamentales en el
desempeo del enlace.
El equipo de recepcin, del cual en este anlisis solo se
contempla la antena, deber ser suficientemente sensible
para percibir el nivel de seal que llega al punto de
recepcin, de manera que el proceso de demodulacin se
pueda llevar a cabo para la recuperacin de la informacin
transmitida.
En la trayectoria del enlace no hay obstrucciones
considerables a la lnea de vista, de manera que la altura
mnima de las antenas requerida para librar el 60% de la
primera zona de Fresnel no representa un problema si se
desea implementar el enlace.
Partiendo de los clculos de prdidas de propagacin,
nivel de intensidad de la seal de recepcin, margen de
desvanecimiento y el criterio de la primera zona de Fresnel
se establece que el radioenlace punto a punto de frecuencia
de operacin de 23 GHz tiene el nivel de seal adecuado
para garantizar la demodulacin en el equipo de recepcin.
El programa implementado en Java permite simular el
radioenlace. Se selecciono este lenguaje para la
implementacin del programa debido a la portabilidad del
mismo.

REFERENCIAS
[1] Gibson J. (1996), The Communications Handbook
Press / IEEE Press
[2] Wayne T. (2003), Sistemas de comunicaciones
electrnicas. Cuarta edicin. Pearson Educacin.
[3] Townsend A. (1996) Digital Line-of-Sight Radio Links:
a handbook. Prentice Hall. UK.
[4] Mapa de curvas de nivel. INEGI.
[5] ARRL (1990), UHF/Microwave Experimenter's
Manual (American Radio Relay League).
[6] Hall, M.P.M., Barclay, L.W. and Hewitt, M.T. (1996),
Propagation of Radiowaves (Institution of Electrical
Engineers,).
[7] Parsons, J.D (1992), The Mobile Radio Propagation
Channel, Wiley & Sons.
[8] Doble, J., (1996) Introduction to Radio Propagation for
Fixed and Mobile Communications Artech House.
[9] Bertoni, H.L., Honcharenko, W., Maciel, L.R. and Xia,
H.H., (1994), UHF Propagation Prediction for Wireless
Personal Communications, Proceedings of the IEEE,
Vol. 82, No. 9, September, pp. 1333-1359.
[10] Freeman, R.L., Radio System Design for
Telecommunications (Wiley & Sons, 1987).
[11] Andersen, J.B., Rappaport, T.S. and Yoshida, S.,
(1995), Propagation Measurements and Models for

Wireless Communications Channels, IEEE

Communications Magazine, January, pp. 42-49.
[12] Lee, W.C.Y., (1993), Mobile Communications Design
Fundamentals, Second Edition (Wiley & Sons, 1993).

APNDICE
COMPUTO EN JAVA DE RADIOENLACES ENTRE
ESTACIONES TERRESTRES
import java.io.*;
import java.lang.Math.*;
public class enlaceT {
public static void main (String arg[]) throws
IOException {
// CLCULO DEL AZIMUT
float latA;
float latB;
float lonA;
float lonB;
double anguloElevacion;
double azimutA;
double azimutB;
double P;
double Lm;
double d;
String cadena;
float flotante;
//
InputStreamReader lector = new
InputStreamReader(System.in);
BufferedReader entrada = new BufferedReader(lector);
System.out.println("\33cCLCULO DE AZIMUT");
System.out.println();
System.out.println("Latitud Punto A: ");
cadena = entrada.readLine();
latA = new Float(cadena).floatValue();
System.out.println("Longitud Punto A: ");
cadena = entrada.readLine();
lonA = new Float(cadena).floatValue();
System.out.println("Latitud Punto B: ");
cadena = entrada.readLine();
latB = new Float(cadena).floatValue();
System.out.println("Longitud Punto B: ");
cadena = entrada.readLine();
lonB = new Float(cadena).floatValue();
System.out.println();
System.out.println("RESULTADOS ");
Lm = (latA+latB)/2;
System.out.println("Lm: "+Lm);
double angle = Lm * 2.0 * Math.PI/360.0;
P = (lonA-lonB) * 96.6 * Math.cos(angle);
System.out.println("P: "+P);
d = Math.sqrt(Math.pow(P,2)+Math.pow((latAlatB)*96.6,2));
System.out.println("d: "+d);
anguloElevacion = Math.acos(P/d);
angle = (anguloElevacion) /( 2.0 * Math.PI)*360.0;
System.out.println("Angulo de elevacion: "+angle);
azimutA=90-angle;
azimutB=360-angle;
System.out.println("Azimut A: "+azimutA);
System.out.println("Azimut B: "+azimutB);
System.out.println();
System.out.println();
// FSL
float frec;
double FSL;
//
System.out.println("CLCULO DE LAS PRDIDAS EN EL
ESPACIO LIBRE (FSL)");
System.out.println("Frecuencia de operacion (GHz): ");
cadena = entrada.readLine();
frec = new Float(cadena).floatValue();
FSL = 92.44 + 20 * Math.log(frec)/Math.log(10) + 20 *
Math.log(d)/Math.log(10);
System.out.println();
System.out.println("RESULTADOS ");
System.out.println("FSL: "+FSL);
System.out.println();
System.out.println();
// PIRE
float Pt;
float Gt;
float Llt;

36

Encuentro de Investigacin en IE, 5 7 de Abril, 2006

double PIRE;
//
System.out.println("CLCULO DE POTENCIA ISOTROPICA
RADIADA EFECTIVAMENTE (PIRE)");
System.out.println("Potencia de transmision (dB): ");
cadena = entrada.readLine();
Pt = new Float(cadena).floatValue();
System.out.println("Ganancia de transmision (dB): ");
cadena = entrada.readLine();
Gt = new Float(cadena).floatValue();
System.out.println("Prdidas en la lnea de transmision
(dB): ");
cadena = entrada.readLine();
Llt = new Float(cadena).floatValue();
PIRE = Pt-Llt+Gt;
System.out.println();
System.out.println("RESULTADOS ");
System.out.println("PIRE: "+PIRE);
System.out.println();
System.out.println();
// RSL
float Gr;
float Llr;
double RSL;
//
System.out.println("CLCULO DEL NIVEL DE RECEPCION DE LA
SENAL (RSL)");
System.out.println("Ganancia de recepcion (dB): ");
cadena = entrada.readLine();
Gr = new Float(cadena).floatValue();
System.out.println("Prdidas en la lnea de recepcion
(dB): ");
cadena = entrada.readLine();
Llr = new Float(cadena).floatValue();
RSL = PIRE-FSL+Gr-Llr;
System.out.println();
System.out.println("RESULTADOS ");
System.out.println("RSL: "+RSL);
System.out.println();
System.out.println();
// FM
float a;
float b;
float w;
double FM;
//
System.out.println("CLCULO DEL MARGEN DE
DESVANECIMIENTO (FM)");
System.out.println("Factor topografico (a): ");
cadena = entrada.readLine();
a = new Float(cadena).floatValue();
System.out.println("Factor climatico (b): ");
cadena = entrada.readLine();
b = new Float(cadena).floatValue();
System.out.println("Disponibilidad del enlace (w): ");
cadena = entrada.readLine();
w = new Float(cadena).floatValue();
FM = 30 * Math.log(d)/Math.log(10) + 10 *
Math.log(6*a*b*frec*1000000000)/Math.log(10);
FM = FM - 10 * Math.log((1-w)*70)/Math.log(10);
System.out.println();
System.out.println("RESULTADOS ");
System.out.println("FM: "+FM);
System.out.println();
System.out.println();
// MARGEN
float S;
double M;
//
System.out.println("CLCULO DEL MARGEN (S)");
System.out.println("Sensibilidad (dB): ");
cadena = entrada.readLine();
S = new Float(cadena).floatValue();
M = RSL - S;
System.out.println();
System.out.println("RESULTADOS ");
System.out.println("MARGEN: "+M);
System.out.println();
System.out.println();
// APERTURA DEL LOBULO Y ANGULO DE MEDIA POTENCIA DEL
LOBULO
float D;
double aBW;
double BW;
//
System.out.println("CLCULO DEL ANGULO DE APERTURA DEL
LOBULO Y ANGULO DE MEDIA POTENCIA DEL LOBULO");
System.out.println("Diametro del plato de la antena (m):
");
cadena = entrada.readLine();

D = new Float(cadena).floatValue();
double longOnda = 300000000/(frec*1000000000);
BW = longOnda*70/D;
aBW = 2*BW;
System.out.println();
System.out.println("RESULTADOS ");
System.out.println("Angulo de apertura del lobulo:
"+BW);
System.out.println("Angulo de media potencia del lobulo:
"+aBW);
System.out.println();
System.out.println();
// PRIMERA ZONA DE FRESNEL Y ALTURA MINIMA DE LAS
ANTENAS
float h0;
float d1;
float d2;
double F1;
double hAnt;
double h;
//
System.out.println("CLCULO DE LA PRIMERA ZONA DEFRESNEL
Y ALTURA MINIMA DE LAS ANTENAS");
System.out.println("Altura del obstaculo mas alto (m):
");
cadena = entrada.readLine();
h0 = new Float(cadena).floatValue();
System.out.println("Distancia del punto A al obstaculo
mas alto (km): ");
cadena = entrada.readLine();
d1 = new Float(cadena).floatValue();
System.out.println("Distancia del punto B al obstaculo
mas alto (km): ");
cadena = entrada.readLine();
d2 = new Float(cadena).floatValue();
F1 = 17.32 * Math.sqrt(d1*d2/(frec*d));
h = d1 * d2 / ( 1.5 * 4/3 );
hAnt = h0+h+0.6*F1;
System.out.println();
System.out.println("RESULTADOS ");
System.out.println("Primera zona de Fresnel (m): "+F1);
System.out.println("Altura minima de las antenas (m):
"+hAnt);
System.out.println();
System.out.println();
}
}
/*
Latitud Punto A:
32.45045
Longitud Punto A:
116.9684722
Latitud Punto B:
32.47027778
Longitud Punto B:
116.8531667
*/