UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS

Fecha: 30/04/2019

MODELOS DE PROPAGACION
Flores Sánchez Manuel Alejandro
[email protected]
Fonseca Rojano Otto Alejandro
[email protected]
Valencia Barragán Karina Lizeth
[email protected]

Carrera: Ingeniería Electrónica

Mención: Telecomunicaciones
Curso: Décimo Grupo 02

Deber N° 2

en el cálculo de las pérdidas en el trayecto y se ha

RESUMEN: En el presente informe se da a conocer el convertido en la planificación de sistemas móviles en
análisis de los modelos de propagación Okumura-Hata Japón. [1]
y Walfish-Ikegami (COST 231), el cual se encuentra
implementado en una simulación en el software Con el objetivo de hacer que este método fuera más
Matlab. fácil de aplicar, Hata estableció una serie de relaciones
numéricas que describen el método gráfico propuesto
PALABRAS CLAVE: Modelos, perdidas, Matlab, por Okumura. Dichas expresiones de carácter
ganancias. empírico, son conocidas bajo el nombre de modelo de
Okumura-Hata, también llamado modelo de Hata. [1]
1. OBJETIVOS
El modelo de Hata es una formulación empírica de las
curvas de pérdidas proporcionadas por Okumura. Hata
1.1 OBJETIVO GENERAL.
estableció una serie de relaciones numéricas que
describen el método gráfico propuesto por Okumura
 Mediante el software Matlab realizar una
con el objetivo de hacer que éste fuera más fácil de
simulación de los modelos de propagación
aplicar. Dichas expresiones son conocidas bajo el
nombre de modelo de Okumura-Hata.
1.2 OBJETIVOS ESPECIFICOS.
Sus límites son:
 Determinar la potencia de transmisión y recepción
de los modelos de propagación ya mencionados.  Frecuencia: 150 a 1500 MHz.
 Altura antena fija: 30 a 200 m.
 Utilizar cualquier modelo propagación de acuerdo  Altura antena móvil: 1 a 10 m.
a los límites que tienen dichos modelos, así el  Distancia: 1 a 20 km.
software decidirá cuál es el mejor modelo para
realizar el análisis. La ecuación 1 describe las pérdidas de este tipo de
modelo de propagación.
 Conocer los parámetros de cada modelo de
predicción de pérdidas de propagación para
aplicarlo en el programa. 𝐿𝑂−𝐻 = 𝐴 + 𝐵 log(𝑑) + 𝐶

2. MARCO TEÒRICO 𝐴 = 69.55 + 26.16 log(𝐹) − 13.82 log(ℎ𝑏 ) − 𝑎(ℎ𝑚 )

2.1 Modelo Okumura-Hata. 𝐵 = 44.9 − 6.55 log(ℎ𝑏 )

Dónde:
A partir de una extensa campaña de medidas llevada (ℎ𝑚 ) = 𝐴𝑙𝑡𝑢𝑟𝑎 𝑑𝑒𝑙 𝑚𝑜𝑣𝑖𝑙
a cabo en Tokio, en el rango de 100 MHz a 1920 MHz, (ℎ𝑏 ) = 𝐴𝑙𝑡𝑢𝑟𝑎 𝑑𝑒 𝑙𝑎 𝑒𝑠𝑡𝑎𝑐𝑖ó𝑛 𝑏𝑎𝑠𝑒
Okumura publica uno de los modelos más utilizados
𝐹 = 𝐹𝑟𝑒𝑐𝑢𝑒𝑛𝑐𝑖𝑎 𝑒𝑛 (𝑀𝐻𝑧)
para la predicción de la pérdida de propagación en
𝑑 = 𝐷𝑖𝑠𝑡𝑎𝑛𝑐𝑖𝑎 (𝐾𝑚)
áreas urbanas. El principal resultado del trabajo de
𝑎(ℎ𝑚 ) = 𝐹𝑎𝑐𝑡𝑜𝑟 𝑑𝑒 𝑐𝑜𝑟𝑟𝑒𝑙𝑎𝑐𝑖ó𝑛
Okumura fue un conjunto de curvas que proporcionan
el nivel de atenuación media relativa al espacio libre,
El factor de correlación y C depende del lugar o área
en función de la frecuencia, la distancia entre
donde se encuentre. Y de aquí se derivan las
transmisor y receptor, la altura de las antenas de la
siguientes formulas. [2]
estación base y la estación móvil, además de varios
factores de corrección específicos para diferentes
Factor de correlación y C en un área metropolitana.
tipos de trayecto. Este modelo está considerado entre
los más simples y mejores en términos de su precisión
a(hm)
UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

en el tejado adyacente a éste; y la pérdida por

8.29(log(1.54 ∗ ℎ𝑚))2 − 1.1 𝑝𝑎𝑟𝑎 𝑓 ≤ 200 𝑀𝐻𝑧 difracción multipuntilla Lmsd, producida por múltiples
difracciones en los tejados de los edificios situados a
3.2(log(11.75 ∗ ℎ𝑚))2 − 4.97 𝑝𝑎𝑟𝑎 𝑓 ≥ 400 𝑀𝐻𝑧 lo largo del trayecto. El modelo distingue además
casos LOS y NLOS. [3]
C=0

Factor de correlación y C en un área de medianas El modelo se también se utiliza para ambientes urbano
o pequeñas ciudades. denso y se base en diversos parámetros como lo son:

𝑎(ℎ𝑚) = (log(𝑓) − 0.7)ℎ𝑚 − (1.56 log(𝑓) − 0.8)  Densidad de los edificios

 Altura promedio de los edificios
C=0  Altura de antenas menor a los edificios.
 Anchura de las calles (w)
Factor de correlación y C en unas áreas  Separación entre los edificios (b)
suburbanas.  Dirección de las calles con respecto a la
trayectoria de la antena transmisora y el
𝑎(ℎ𝑚) = (log(𝑓) − 0.7)ℎ𝑚 − (1.56 log(𝑓) − 0.8) móvil.
En general, los límites para el modelo son las
siguientes:
𝑓 2
𝐶 = −2 (log ( )) − 5.4
28  f=800–2000MHz
 hb=4–50m (altura de transmisor)
Factor de correlación y C en unas áreas rurales.  hm=1–3m (altura de receptor)
 d=0.02–50km (distancia entre transmisor y
𝑎(ℎ𝑚) = (log(𝑓) − 0.7)ℎ𝑚 − (1.56 log(𝑓) − 0.8) receptor)
 Δhb > 0 m (altura relativa del transmisor respecto
𝐶 = −4.78(log(𝑓))2 + 18.33 log(𝑓) − 40.84 a los edificios).

Las pérdidas se modelan en las ecuaciones

2.2 MODELO WALFISCH-IKEGAMI (COST 231). siguientes:

El modelo COST 231 es un modelo semi-empírico de 1) Cuando hay línea de vista en las antenas:
predicción de las pérdidas en un trayecto, resultado de
la combinación de los modelos Walfisch-Bertoni e 𝐿𝑜 = 42.6 + 26 log(𝑑) + 20 log(𝑓)
Ikegami.
2) Cuando no hay línea de vista
Es recomendado para macro-células en escenarios
urbanos y suburbanos, con buenos resultados de las 𝐿𝑏 = 𝐿𝑜 + 𝐿𝑟𝑠𝑡 + 𝐿𝑚𝑠𝑑 𝑝𝑎𝑟𝑎 𝐿𝑦𝑡𝑠 + 𝐿𝑚𝑠𝑑 > 0
pérdidas en el trayecto para antenas transmisoras 𝐿𝑏 = 𝐿𝑜 𝑝𝑎𝑟𝑎 𝐿𝑦𝑡𝑠 + 𝐿𝑚𝑠𝑑 ≤ 0
situadas por encima de la altura media de los tejados.
Sin embargo, el error en las predicciones aumenta Dónde:
considerablemente a medida que la altura del
transmisor se acerca a la altura de los tejados, Lo: Pérdidas por el espacio libre
llegando a tener un rendimiento muy pobre para Lrts: Pérdidas por difracción de múltiples esquinas de
transmisores situados por debajo de ese nivel. los techos de los edificios.
Lmsd: Pérdidas debido a una única difracción final
cuando la onda se propaga hacia la calle.
 Amplia el rango de frecuencias en el cual puede
usarse (800 - 2000 MHz). El modelo realiza un 𝐿𝑟𝑡𝑠 = −8.2 −10𝑙𝑜(w)+ 10𝑙𝑜𝑔(𝑓) + 20𝑙𝑜𝑔(∆ℎm) +
cálculo más detallado de la atenuación,
𝐿𝑜𝑟i
basándose en cuatro parámetros adicionales:
Donde Lori es el factor de correlación y depende
 Altura de los edificios
de 𝝋.
 Ancho de las calles
 Separación entre edificios
−10 + 0.3571 𝜑. 0 ≤ 𝜑.35°
 Orientación de la calle respecto a la dirección de
propagación
2,5 + 0,075(𝜑. − 35) 35 ≤ 𝜑. < 55°

4 − 0,114(𝜑. − 55) 55 ≤ 𝜑. < 90°

La pérdida básica de propagación, Lb, es calculada
como la suma de tres componentes: la pérdida de Donde Lmsd
propagación en condiciones de espacio libre, Lo; la
pérdida por difracción del tejado a la calle, Lrts, 𝐿𝑚𝑠𝑑=𝐿𝑏𝑠ℎ + 𝑘𝑎 + 𝑘𝑑 log (𝑑) + 𝑘𝑓 log (𝑓) − 9 log
producida en el interior de la calle en la que se
(d)
encuentra el receptor, como resultado de la difracción
UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

Donde Lbsh  𝝋 : ángulo formado por la dirección y

propagación y el eje de la calle (grados)
𝐿𝑏𝑠ℎ = −18 log (1 + ∆ℎ𝐵) 𝑑𝑜𝑛𝑑𝑒 ℎb > ℎr  hb: altura de la antena estación base
 hm: altura de la antena dispositivo móvil
 ∆hm= hr-hm (m)
𝐿𝑏𝑠ℎ =0 𝑑𝑜𝑛𝑑𝑒 ℎb < ℎr  ∆hb= hb –hr (m)
 L: distancia total entre el primer y el ultimo
Para Ka edificio (m)
 d: distancia entre estación basen y
54 ℎb > ℎ𝑟 dispositivo móvil
 f: frecuencia (MHz)
54 − 0,8 ∆ℎb ≤ ℎ𝑟 𝑦 𝑑 ≥ 0,5

54 – 0.8∆ℎb 𝑑/0.5 ℎb ≤ ℎ𝑟 𝑦 𝑑 ≥ 0,5 2.3 MODELO WALFISCH-BERTONI

Para Kd Este modelo, fue propuesto por Joram Walfisch y Henri

Bertoni, tiene en cuenta las pérdidas producidas por
18 ℎb > ℎ𝑟 las difracciones que se producen en las azoteas de los
edificios. Es un modelo que no considera la existencia
18 − 15 ∆ℎb/hr ℎb ≤ ℎ𝑟 de línea de vista entre el transmisor y el receptor,
utiliza el fenómeno de la difracción para describir las
pérdidas que sufre la señal antes de alcanzar al
Para la constante 𝑘𝑓se hace la distinción de ciudades receptor ubicado a baja altura sobre la calle. Es
grandes y medianas, así en el caso de ciudades de despreciada la contribución de los rayos que penetran
tamaño medio y centro suburbano con densidad en los edificios y de los que sufren múltiples
moderada de vegetación: [3] difracciones.

Para kf Los rangos para este modelo son los siguientes:

𝑘𝑓 = −4 + 0,7 (𝑓/ 925 − 1) f = 300 – 3000 MHz

hb = 4 – 50 m (altura de transmisor)
hm = 1 – 3 m (altura de receptor)
𝑘𝑓 = −4 + 1,5 (𝑓/ 925 − 1) d = 0.2 – 5 km (distancia entre transmisor y receptor)

Fig 3. Escenario típico del modelo Walfish-Bertoni

 d: distancia entre Tx y Rx en KM
 f: frecuencia (MHz)
 H: altura promedio de la antena con respecto
Fig 1. Escenario típico NLOS en vista de perfil a la altura de los edificios.
 A: variable que expresa la influencia de los
edificios en la señala
 hb: altura de la antena estación base
 hm: altura de la antena dispositivo móvil
 b: espacio entre edificios

en donde la fórmula para la variable A esta expresada

por:

𝑏 2
𝐴 = 5 log [( ) + (ℎ𝑏 − ℎ𝑟)2 ]-9 log b +120
2
Fig 2. Escenario típico NLOS en vista superior. 2(ℎ𝑏−ℎ𝑟)
log{tan−1 [ ]}
𝑏
Los parámetros definidos son los siguientes:

 hr: altura media de los edificios (m). Y la fórmula para las perdidas es:
 w: anchura de la calle (m).
 b: separación entre edificios(m)
UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

Lp= 89.55 + 21log(𝑓𝑀ℎ𝑧)+38 log (d km)-18 log H + A planificación de sistemas móviles y sistemas de radio
𝑑2 terrestre en ambientes poblados a nivel mundial.
- 18 log{1 − }
17𝐻 Proporciona el nivel de atenuación media relativa al
espacio libre en función de la frecuencia, la distancia
Tx y Rx la h altura de las antenas de la estación base
2.4 MODELO MULTIPLE DE ERICSSON. y la estación móvil.

Características: 𝐿𝑝 𝑑𝐵 = 𝐿𝐹 + 𝐴𝑚𝑢 − 𝐺(ℎ𝑡𝑒 ) − G(ℎ𝑟𝑒 ) − 𝐺𝐴𝑅𝐸𝐴

Fue desarrollado por ingenieros de Ericsson Donde
basándose en el modelo de Okumura Hata extendido. • 𝐿𝑝 perdidas por trayectoria en decibeles.
es un modelo muy sencillo donde su exactitud queda • 𝐿𝐹 perdidas en el espacio libre
determinada por el correcto ajuste de los parámetros • 𝐺(ℎ𝑡𝑒 ) ganancia de la altura de la antena de
libres en base a mediciones para cada región. transmisión (dB)
• G(ℎ𝑟𝑒 ) ganancia de la altura de la antena en
Los rangos para este modelo son los siguientes: recepción
• 𝐺𝐴𝑅𝐸𝐴 ganancia del entorno
Frecuencia (150-2000 MHz) • 𝐴𝑀𝑈 Atenuación relativa promedio
Distancia (0,2-100 km)
Altura antenas estación base (20-200m) Okumura encontró que 𝐺(ℎ𝑡𝑒 ) tiene una variación de
Altura antenas estación móvil (1-5m) perdidas de 20 dB / década y que G(ℎ𝑟𝑒 ) tiene una
variación de 10 dB /década para alturas menores a 3m.
El modelo puede ser descrito por contribuciones a las
pérdidas.

 Pérdidas adicionales que se presentan

cuando la propagación es modificada debido
a picos de montaña.
 Para distancias mayores a 10 km aparecerá
pérdidas adicionales debido a los disturbios
causados por la curvatura de la tierra.
 Perdidas por la topología de la zona. 2.6 MODELO LONGLEY RICE.

Formulas: El modelo Longley Rice, también llamado de ITM

(Irregular Terrain Model), es un modelo matemático
a0, a1, a2, a3 - parámetros del modelo desarrollado por Longley y Rice, en 1968, para el
estudio de propagación de señales de TV, modela
a0 a1 a2 a3 obstáculos lejanos y cercanos.
Urbano 36.2 30.2 -12 0.1 Permite realizar cálculos de propagación para enlaces
Suburbano 43.2 68.93 -12 0.1 punto a punto y punto a zona, debido a sus
Rural 45.95 100.6 -12 0.1 características, es utilizado como motor de cálculos de
varios programas de telecomunicaciones como SPLAT
𝑃𝐿 = 𝑎0 + 𝑎1 log(𝑑) y Radiomobile.
+ 𝑎2 log(ℎ𝑏) Se aplica a esquemas de comunicación en el rango de
+ 𝑎3 log(ℎ𝑏) log(𝑑) frecuencias desde VHF hasta EHF, es decir desde lo
− 3.2(log(11.75 ∗ ℎ𝑚))2 + 𝑔(𝑓) 40MHz hasta los 100GHz, sobre diferentes tipos de
terrenos.
Donde:
𝑔(𝑓) = 44.49 log(𝑓) − 4.78(log(𝑓))2 Para el cálculo de las perdidas el modelo usa la teoría
de la difracción, la refracción troposférica y el
scattering del terreno.
2.5 MODELO OKUMURA.
Las perdidas adicionales están basadas en medidas
tomadas en varias situaciones.
El modelo de Okumura es utilizado para predecir la
potencia en un receptor ubicado en un área urbana Los rangos para este modelo son los siguientes:
para comunicaciones móviles. Este modelo es
aplicable para el rango de frecuencias entre 150 a f =200 MHz-20 GHz
1920 MHz es decir comprende la banda de VHF Y d = 1-2000 Km
UHF. Según este modelo, la distancia máxima de Altura de antena=0.5 m – 3 km
separación que puede existir entre el transmisor y el
receptor es de hasta 100 km. 1uede ser usado para L(dB)=Lel+Acr
alturas de la antena de la
estación base en el rango de 30 m a 100 m. Donde Lel se refiere a las pérdidas en espacio libre,
calculadas por medio de la ecuación:
Este modelo está considerado entre los más simples y Lel(dB)=32,45+20logf+20logd
mejores en términos de su precisión en el cálculo de
las perdidas en el trayecto y se ha convertido en la Mientras que Acr representa la atenuación relativa del
espacio libre y puede ser calculada, empleando algún
UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

mecanismo de propagación, de acuerdo a los GTx=18

siguientes rangos de distancia. GRx=12

Dentro de la línea de vista: se emplea la ecuación de • Perdidas Hata

pérdidas de propagación del modelo de dos rayos (un
rayo directo y un rayo reflejado), también conocido LoH=A+B*log(d “km.”)+ C
como modelo de tierra plana:
• Determinar el valor de A
L(dB)=120−20loghT+20loghR+40logd

Donde “d” está dada en km. hT, hR (m) es la altura de A=69.55+26.16*log(f”MHz.”)-13.62*log(hb)-

la antena transmisora y receptora. a(hm)

Fuera de la línea de vista: se recurre al mecanismo • Se procede a calcular a(hm)

de propagación de difracción, calculando un valor
promedio estimado de todas las atenuaciones por a(hm)=(1.1*log(f)-0.7)*2-(1.56*log(f)-0.8)
difracción causadas por el terreno irregular o por
obstáculos tipo “filo de cuchillo” interpuestos en la
a(hm)=(1.1*log(1000)-0.7)*2-
trayectoria directa de la señal entre transmisor y
(1.56*log(1000)-0.8)
receptor.

En el caso de objetos agudos se puede determinar a(hm)= 1.32 dB

teóricamente las perdidas debidas a la difracción.
A=69.55+26.16*log(1000)-13.62*log(85)-
La diferencia de caminos se expresa con la diferencia 1.32
de fase entre las señal emitida y recibida:
A=120.43 dB

• Determinar el valor de B

La diferencia de fase se expresa como: B=44.9 – 6.55log(hb)

B=44.9 – 6.55log(85)

Se expresa en función del termino: B=32.26 dB

• Determinar el valor de C

C=-4.78 log f^2+18.33log f-40.94

con v el parámetro de difracción de Fresnel-Kirchhoff: C=-4.78 log (1000)^2+18.33log (1000)-40.94

C=-28.97

• LoH=120.43+32.26*log(15)+ (-28.97)

• LoH=129.40
3. MATERIALES Y EQUIPO
• Determinar la potencia con la que se debe
 Software Matlab transmitir

• Prx= -26.98
4. DESARROLLO Y PROCEDIMIENTO
Ptx=Prx-Gtx-Grx+LoH
 Como primer punto se analizó los modelos
con un ejercicio práctico matemáticamente. Ptx=-30-0-0+129.4

Modelo Okumura Hata

Ptx=43.08 dB.

Datos: Modelo Cost 231

f=500MHz
d= 15 Km Un sistema inalámbrico NLOS en un área
hb= 35 m. metropolitana trabajando a f=1500MHZ, y antenas con
hm= 1.5m. ganancia de 18dB en transmisión y recepción, para
Prx=0.002 dB poder establecer una comunicación fiable necesita
UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

10μW. Si la altura de la antena para el móvil es de 2m. Ptx(dB) = -50 – 18 – 18 + 209.90037

Si se conoce que el primer edificio se encuentra a
20km del último edificio, los edificios tienen un Ptx = = -50 + 173.90037= 123.90037 dB
promedio de 6 pisos y techo inclinado, la separación
promedio entre edificios es de 20m. Calcular la
potencia de transmisión cuando la antena de la BS Modelo Walfish-Bertoni
está ubicada 3m sobre el tejado de un edificio, si el
receptor se encuentra a una distancia de 25Km. Datos:

• DATOS H=15 m,
f=1500MHz hb=16m,
d=25Km hr=12 m
G=18dB b=2 m
Prx=10uW f=350 MHz
hm=2m d =5 Km
b=20m
w=b/2 = 10m 2 2
l=20Km = 20000m 𝐴 = 5log [( ) + (16 − 12)2 ]-9 log 2 + 20
2
hr= 21 m 2(16−12)
log{tan−1 [ ]}
fhi = 90° 2
hb=24m
hm =2m 𝐴 = 6.15-2.71 + 37.61=41.0

Lp= 89.55 + 21log(350)+38 log (5)-18 log 15 + 41.05

Lo = 32,4 + 20 log(d) + 20 log f 52
- 18 log{1 − }
17𝑥15
Lo = 32,4 + 20 log(25) + 20 log 1500
Lp= 89.55 + 53.42+ 26.56- 21.16 + 41.05 – (-0.806) =
Lo = 123.8806 190.226

Modelo múltiple de Ericsson

Lrts = -8,2 – 10 log w + 10 log f + 20 log (Δhm) + Lori
hb=30m
Lori= 4 – 0.114 (fhi - 55) = 0,01 hr=3m
d=5km
Lrts = -8,2 – 10 log(10) + 10 log (1500) + 20 log (19) + f=1900MHz
0,01 hr=15m
a0=43,20
Lrts = 39,14598 a1=68,93
a2=-12.0
Lmsd= Lbsh + ka + kd log(d) + kf log (f) -9 log (b) a3=0.1

𝑃𝐿 = 𝑎0 + 𝑎1 log(5)
Lbsh = -18 log (1 + Δhb) = -18 log (1+ 3) = -10,8370
+ 𝑎2 log(30)
hb > hr 24m > 21 m
+ 𝑎3 log(30) log(5)
− 3.2(log(11.75 ∗ 15))2 + 𝑔(𝑓)
Ka = 54
𝑃𝐿 = 150.68 𝑑𝑏
Kd = 18
Modelo Longley Rice
Kf = -4 + 1,5(f/925-1) = -4 + 1,5 (1500/925 – 1 ) = -
3,06756

Lmsd= -10,8370 + 54 + 18 log(25) + (-3.06756) log

(1500) -9 log (20)

Lmsd = 46,87369

Lrts + Lmsd = 39,14598 + 46,87369 = 86,019779 >

0

Lb = Lo + Lrts + Lmsd = 123.8806 + 39,14598 +

46,87369 = 209.90037 (dB)

Prx = 10uW = -50 dB

Ptx(dB) = Prx(dB) – Gtx(dB) – Grx(dB) + Lb(dB)

UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

if (f>=150 && f<=2000)&&(hb>=30 && hb<=200) &&(hm>=1 &&

hm<=10) &&(d>=1 && d<=20)
zf=1;
end
if (f>=150 && f<=1920)&&(hb>=30 && hb<=100) &&(hm>=1 &&
hm<=100) &&(d>=10.02 && d<=100)
zg=1;
end
if (za==1)&&(zb==1)&&(zc==1)&&(ze==1)&&(zf==1)&&(zg==1)
disp('El modelo que se puede utilizar es de:');
disp('1.Okura-Hata')
disp('2.Wlafish-Ikegami')
disp('3.Dos Rayos')
disp('4.Difraccion Knife- Edge')
disp('5.Wlafish-Bertoni')
disp('6.Ericsson')
disp('7.Okumura')
dato=input('');

elseif (za==1)&&(zb==1)&&(zc==1)&&(ze==1)&&(zg==1)
disp('El modelo que se puede utilizar es de:');
disp('1.Okura-Hata')
disp('2.Wlafish-Ikegami')
disp('3.Dos Rayos')
disp('4.Difraccion Knife- Edge')
disp('5.Wlafish-Bertoni')
disp('7.Okumura')
dato=input('');
Modelo de okumura
elseif (za==1)&&(zb==1)&&(zc==1)&&(zf==1)&&(zg==1)
disp('El modelo que se puede utilizar es de:');
Datos: disp('1.Okura-Hata')
disp('2.Wlafish-Ikegami')
disp('3.Dos Rayos')
d = 50 Km disp('4.Difraccion Knife- Edge')
Htx = 100 m disp('6.Ericsson')
disp('7.Okumura')
Hrx = 10 m dato=input('');
Entorno suburbano
elseif (za==1)&&(zb==1)&&(ze==1)&&(zf==1)&&(zg==1)
PIRE= 1 Kw = 30 db disp('El modelo que se puede utilizar es de:');
f= 900 Mhz disp('1.Okura-Hata')
disp('2.Wlafish-Ikegami')
Grx= 1 w = 0 db disp('4.Difraccion Knife- Edge')
disp('5.Wlafish-Bertoni')
disp('6.Ericsson')
L50(db) = FLS + Amu(f,d) - G (Ht) - G (Hr) – Garea disp('7.Okumura')
FLS = 32.5 +20 log f + 20 log d dato=input('');

FLS = 32.5 +20 log 900 + 20 log 50 elseif (za==1)&&(zc==1)&&(ze==1)&&(zf==1)&&(zg==1)

FLS = 32.5 + 59.08 + 33.98 disp('El modelo que se puede utilizar es de:');
disp('1.Okura-Hata')
FLS = 125.56 db disp('3.Dos Rayos')
G (Ht) = 20 log (Ht/200) = 20 log (100/200) = -6 db disp('4.Difraccion Knife- Edge')
disp('5.Wlafish-Bertoni')
G (Hr) = 20 log (Ht/3) = 20 log (10/3) = 10.46 db disp('6.Ericsson')
L50(db) = FLS + Amu(f,d) - G (Ht) - G (Hr) – Garea disp('7.Okumura')
dato=input('');
L50(db) = 125.5 + 43 –(-6) - 10.46 – 9
elseif (zb==1)&&(zc==1)&&(ze==1)&&(zf==1)&&(zg==1)
L50(db) = 155.04 db disp('El modelo que se puede utilizar es de:');
Prx = PIRE – L50 + Grx
disp('2.Wlafish-Ikegami')
Prx = 30 – 155.04 + 0 disp('3.Dos Rayos')
Prx = -125.04 dbw disp('4.Difraccion Knife- Edge')
disp('5.Wlafish-Bertoni')
disp('6.Ericsson')
disp('7.Okumura')
dato=input('');
Como segundo punto se procedió a realizar el
elseif (za==1)&&(zb==1)&&(zc==1)&&(ze==1)&&(zf==1)
programa en el software Matlab. disp('El modelo que se puede utilizar es de:');
disp('1.Okura-Hata')
disp('2.Wlafish-Ikegami')
disp('Modelos de propagación') disp('3.Dos Rayos')
f=input('Ingrese la frecuencia que se utilizara en Mhz: '); disp('4.Difraccion Knife- Edge')
hb=input ('ingrese la altura de la base en m:'); disp('5.Wlafish-Bertoni')
hm=input ('ingrese la altura del movil en m:'); disp('6.Ericsson')
d=input('Ingrese la distacia en km:'); dato=input('');
za=0;
zb=0;
zc=0; elseif (za==1)&&(zb==1)&&(zc==1)&&(zg==1)
zd=1;%cuchillas disp('El modelo que se puede utilizar es de:');
ze=0; disp('1.Okura-Hata')
zf=0; disp('2.Wlafish-Ikegami')
zg=0; disp('3.Dos Rayos')
disp('4.Difraccion Knife- Edge')
if (f>=150 && f<=1500)&&(hb>=30 && hb<=200) &&(hm>=1 && disp('7.Okumura')
hm<=10) &&(d>=1 && d<=20) dato=input('');
za=1;
end
if (f>=800 && f<=2000)&&(hb>=4 && hb<=50) &&(hm>=1 && elseif (za==1)&&(zb==1)&&(zf==1)&&(zg==1)
hm<=3) &&(d>=0.02 && d<=50) disp('El modelo que se puede utilizar es de:');
zb=1; disp('1.Okura-Hata')
end disp('2.Wlafish-Ikegami')
disp('4.Difraccion Knife- Edge')
if (f>=150 && f<=3000)&&(hb>=50 && hb<=200) &&(hm>=1 && disp('6.Ericsson')
hm<=100) &&(d>=5 && d<=100) disp('7.Okumura')
zc=1; dato=input('');
end
if (f>=300 && f<=3000)&&(hb>=4 && hb<=50) &&(hm>=1 && elseif (za==1)&&(ze==1)&&(zf==1)&&(zg==1)
hm<=3) &&(d>=0.2 && d<=5) disp('El modelo que se puede utilizar es de:');
ze=1; disp('1.Okura-Hata')
end disp('4.Difraccion Knife- Edge')
disp('5.Wlafish-Bertoni')
UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

disp('6.Ericsson') disp('El modelo que se puede utilizar es de:');

disp('7.Okumura') disp('1.Okura-Hata')
dato=input(''); disp('2.Wlafish-Ikegami')

disp('4.Difraccion Knife- Edge')

elseif (zc==1)&&(ze==1)&&(zf==1)&&(zg==1) disp('5.Wlafish-Bertoni')
disp('El modelo que se puede utilizar es de:'); disp('6.Ericsson')

disp('3.Dos Rayos') dato=input('');

disp('4.Difraccion Knife- Edge')
disp('5.Wlafish-Bertoni')
disp('6.Ericsson') elseif (za==1)&&(zb==1)&&(zc==1)&&(zf==1)
disp('7.Okumura') disp('El modelo que se puede utilizar es de:');
dato=input(''); disp('1.Okura-Hata')
disp('2.Wlafish-Ikegami')
elseif (zb==1)&&(ze==1)&&(zf==1)&&(zg==1) disp('3.Dos Rayos')
disp('El modelo que se puede utilizar es de:'); disp('4.Difraccion Knife- Edge')

disp('2.Wlafish-Ikegami') disp('6.Ericsson')

disp('4.Difraccion Knife- Edge') dato=input('');

disp('5.Wlafish-Bertoni')
disp('6.Ericsson') elseif (za==1)&&(zb==1)&&(zc==1)&&(ze==1)
disp('7.Okumura') disp('El modelo que se puede utilizar es de:');
dato=input(''); disp('1.Okura-Hata')
disp('2.Wlafish-Ikegami')
disp('3.Dos Rayos')
elseif (zb==1)&&(zc==1)&&(zf==1)&&(zg==1) disp('4.Difraccion Knife- Edge')
disp('El modelo que se puede utilizar es de:'); disp('5.Wlafish-Bertoni')

disp('2.Wlafish-Ikegami') dato=input('');
disp('3.Dos Rayos')
disp('4.Difraccion Knife- Edge') elseif (za==1)&&(zb==1)&&(zg==1)
disp('El modelo que se puede utilizar es de:');
disp('6.Ericsson') disp('1.Okura-Hata')
disp('7.Okumura') disp('2.Wlafish-Ikegami')
dato=input('');
disp('4.Difraccion Knife- Edge')
elseif (zb==1)&&(zc==1)&&(ze==1)&&(zg==1)
disp('El modelo que se puede utilizar es de:'); disp('7.Okumura')
dato=input('');
disp('2.Wlafish-Ikegami')
disp('3.Dos Rayos') elseif (za==1)&&(zf==1)&&(zg==1)
disp('4.Difraccion Knife- Edge') disp('El modelo que se puede utilizar es de:');
disp('5.Wlafish-Bertoni') disp('1.Okura-Hata')

disp('7.Okumura') disp('3.Dos Rayos')

dato=input('');
disp('6.Ericsson')
elseif (zb==1)&&(zc==1)&&(ze==1)&&(zf==1) disp('7.Okumura')
disp('El modelo que se puede utilizar es de:'); dato=input('');
elseif (ze==1)&&(zf==1)&&(zg==1)
disp('2.Wlafish-Ikegami') disp('El modelo que se puede utilizar es de:');
disp('3.Dos Rayos')
disp('4.Difraccion Knife- Edge') disp('4.Difraccion Knife- Edge')
disp('5.Wlafish-Bertoni') disp('5.Wlafish-Bertoni')
disp('6.Ericsson') disp('6.Ericsson')
disp('7.Okumura')
dato=input(''); dato=input('');

elseif (za==1)&&(zc==1)&&(zf==1)&&(zg==1) elseif (zb==1)&&(zf==1)&&(zg==1)

disp('El modelo que se puede utilizar es de:'); disp('El modelo que se puede utilizar es de:');
disp('1.Okura-Hata')
disp('2.Wlafish-Ikegami')
disp('3.Dos Rayos')
disp('4.Difraccion Knife- Edge') disp('4.Difraccion Knife- Edge')

disp('6.Ericsson') disp('6.Ericsson')
disp('7.Okumura') disp('7.Okumura')
dato=input(''); dato=input('');

elseif (za==1)&&(zc==1)&&(ze==1)&&(zg==1) elseif (zb==1)&&(zc==1)&&(zg==1)

disp('El modelo que se puede utilizar es de:'); disp('El modelo que se puede utilizar es de:');
disp('1.Okura-Hata')
disp('2.Wlafish-Ikegami')
disp('3.Dos Rayos') disp('3.Dos Rayos')
disp('4.Difraccion Knife- Edge') disp('4.Difraccion Knife- Edge')
disp('5.Wlafish-Bertoni')
disp('7.Okumura')
disp('7.Okumura') dato=input('');
dato=input('');
elseif (zb==1)&&(zc==1)&&(ze==1)
elseif (za==1)&&(zc==1)&&(ze==1)&&(zf==1) disp('El modelo que se puede utilizar es de:');
disp('El modelo que se puede utilizar es de:');
disp('1.Okumura-Hata') disp('2.Wlafish-Ikegami')
disp('3.Dos Rayos')
disp('3.Dos Rayos') disp('4.Difraccion Knife- Edge')
disp('4.Difraccion Knife- Edge') disp('5.Wlafish-Bertoni')
disp('5.Wlafish-Bertoni')
disp('6.Ericsson') dato=input('');

dato=input(''); elseif (za==1)&&(zc==1)&&(zg==1)

disp('El modelo que se puede utilizar es de:');
elseif (za==1)&&(zb==1)&&(ze==1)&&(zg==1) disp('1.Okura-Hata')
disp('El modelo que se puede utilizar es de:');
disp('1.Okura-Hata') disp('3.Dos Rayos')
disp('2.Wlafish-Ikegami') disp('4.Difraccion Knife- Edge')

disp('4.Difraccion Knife- Edge') disp('7.Okumura')

disp('5.Wlafish-Bertoni') dato=input('');
elseif (za==1)&&(zc==1)&&(ze==1)
disp('7.Okumura') disp('El modelo que se puede utilizar es de:');
dato=input('');
disp('2.Wlafish-Ikegami')
elseif (za==1)&&(zb==1)&&(ze==1)&&(zf==1) disp('3.Dos Rayos')
UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

disp('4.Difraccion Knife- Edge') disp('El modelo que se puede utilizar es de:');

disp('5.Wlafish-Bertoni')
disp('2.Wlafish-Ikegami')
dato=input('');

elseif (za==1)&&(zb==1)&&(ze==1) disp('6.Ericsson')

disp('El modelo que se puede utilizar es de:');
disp('1.Okura-Hata') dato=input('');
disp('2.Wlafish-Ikegami')
elseif (zb==1)&&(zg==1)
disp('4.Difraccion Knife- Edge') disp('El modelo que se puede utilizar es de:');
disp('5.Wlafish-Bertoni')
disp('2.Wlafish-Ikegami')
dato=input('');
disp('4.Difraccion Knife- Edge')
elseif (za==1)&&(zb==1)&&(zc==1)
disp('El modelo que se puede utilizar es de:'); disp('7.Okumura')
disp('1.Okura-Hata') dato=input('');
disp('2.Wlafish-Ikegami')
disp('3.Dos Rayos')
disp('4.Difraccion Knife- Edge') elseif (zc==1)&&(zf==1)
disp('El modelo que se puede utilizar es de:');
dato=input('');
disp('3.Dos Rayos')
elseif (zf==1)&&(zg==1) disp('4.Difraccion Knife- Edge')
disp('El modelo que se puede utilizar es de:');
disp('6.Ericsson')
disp('6.Ericsson')
disp('7.Okumura') dato=input('');
dato=input('');
elseif (zc==1)&&(zg==1)
elseif (ze==1)&&(zf==1) disp('El modelo que se puede utilizar es de:');

disp('4.Difraccion Knife- Edge') disp('3.Dos Rayos')

disp('5.Wlafish-Bertoni') disp('4.Difraccion Knife- Edge')
disp('6.Ericsson')
disp('7.Okumura')
dato=input(''); dato=input('');

elseif (zc==1)&&(ze==1) elseif (ze==1)&&(zg==1)

disp('El modelo que se puede utilizar es de:'); disp('El modelo que se puede utilizar es de:');
disp('3.Dos Rayos')
disp('4.Difraccion Knife- Edge')
disp('5.Wlafish-Bertoni') disp('4.Difraccion Knife- Edge')
dato=input(''); disp('5.Wlafish-Bertoni')

elseif (zb==1)&&(zc==1) disp('7.Okumura')

disp('El modelo que se puede utilizar es de:'); dato=input('');
disp('2.Wlafish-Ikegami')
disp('3.Dos Rayos')
disp('4.Difraccion Knife- Edge')
elseif (za==1)&&(zb==1)
dato=input(''); disp('El modelo que se puede utilizar es de:');
disp('1.Okura-Hata')
disp('2.Wlafish-Ikegami')
disp('4.Difraccion Knife- Edge')
elseif (za==1)&&(zc==1)
disp('El modelo que se puede utilizar es de:'); dato=input('');
disp('1.Okura-Hata')
elseif (za==1)
disp('3.Dos Rayos') disp('El modelo que se puede utilizar es de:');
disp('4.Difraccion Knife- Edge') disp('1.Okura-Hata')
disp('4.Difraccion Knife- Edge')
dato=input(''); dato=input('');

elseif (za==1)&&(ze==1)
disp('El modelo que se puede utilizar es de:'); elseif (zb==1)
disp('1.Okura-Hata') disp('El modelo que se puede utilizar es de:');
disp('2.Wlafish-Ikegami')
disp('4.Difraccion Knife- Edge') disp('4.Difraccion Knife- Edge')
disp('5.Wlafish-Bertoni') dato=input('');

dato=input(''); elseif (zc==1)

disp('El modelo que se puede utilizar es de:');
elseif (za==1)&&(zf==1)
disp('El modelo que se puede utilizar es de:'); disp('3.Dos Rayos')
disp('1.Okura-Hata')
disp('4.Difraccion Knife- Edge')
disp('4.Difraccion Knife- Edge') dato=input('');

disp('6.Ericsson')
elseif (ze==1)
dato=input(''); disp('El modelo que se puede utilizar es de:');

elseif (za==1)&&(zg==1) disp('4.Difraccion Knife- Edge')

disp('El modelo que se puede utilizar es de:'); disp('5.Wlafish-Bertoni')
disp('1.Okura-Hata') dato=input('');
elseif (zf==1)
disp('4.Difraccion Knife- Edge') disp('El modelo que se puede utilizar es de:');
disp('4.Difraccion Knife- Edge')
disp('7.Okumura') disp('6.Ericsson')
dato=input('');
dato=input('');
elseif (zg==1)
elseif (zb==1)&&(ze==1) disp('El modelo que se puede utilizar es de:');
disp('El modelo que se puede utilizar es de:'); disp('4.Difraccion Knife- Edge')
disp('7.Okumura')
disp('2.Wlafish-Ikegami') dato=input('');
elseif (zd==1)
disp('4.Difraccion Knife- Edge') disp('El modelo que se puede utilizar es de:');
disp('5.Wlafish-Bertoni') disp('4.Difraccion Knife- Edge')
dato=input('');
dato=input(''); else
disp('No se adapta a ningun modelo')
elseif (zb==1)&&(zf==1) end
UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

end

fprintf('Perdidas de propagacion
es espacio libre (Lrts): %.2f (dB) \n' , Lrts)
if (hb>hr)
Lbsh=-18*log10(1+Var_hb);
switch dato elseif(hb<=hr)
case 1 Lbsh=0;
disp('Modelo de OKUMURA-HATA') end
disp('1 METROPOLITANO') fprintf('Lbsh es igual a: %.2f
disp('2 SMALL/MED') (dB) \n' , Lbsh)
disp('3 SUBURBAN')
disp('4 RURAL') if(hb>hr)
hata=input('Ingrese una opcion: '); ka=54;
switch hata elseif (hb<=hr &&d >=0.5)
case 1 ka=54-0.8*Var_hb;
if(f<=200) elseif (hb<=hr && d<0.5)
a=8.29*(log10(1.54*hm))^2- ka=(54-(0.8*Var_hb)*(d/5));
1.1 end
c=0
elseif (f>=400) if (hb>hr)
a=3.2*(log10(11.75*hm))^2- kd=18;
4.97 elseif (hb<=hr)
c=0 kd=18-15*((Var_hb)/hr);
end
case 2 end
a=(1.1*log10(f)-0.7)*hm-
(1.56*log10(f)-0.8) disp('1 ciudad o entorno
c=0 suburbano')
case 3 disp('2 centro metropolitano')
a=(1.1*log10(f)-0.7)*hm- numero=input('Ingrese una
(1.56*log10(f)-0.8) opcion: ');
c=-2*(log10(f)/28)^2-5.4 switch numero
case 4 case 1
a=(1.1*log10(f)-0.7)*hm- kf=-4+0.7*((f/925)-1);
(1.56*log10(f)-0.8) case 2
c=- kf=-4+1.5*((f/925)-1);
4.78*(log10(f))^2+18.33*log10(f)-40.94 end
fprintf('ka es igual a: %.2f \n'
end , ka)
A=69.55+26.16*log10(f)-13.82*log10(hb)-a fprintf('kd es igual a: %.2f \n'
B=44.9-6.55*log10(hb) , kd)
Lo=A+B*log10(d)+c fprintf('kf es igual a: %.2f \n'
, kf)
PTX=PRX-GRX-GTX+Lo;
fprintf('La Potencia de Transmision
(PTX): %.2f(dB) \n' , PTX) Lmsd=Lbsh+ka+kd*log10(d)+kf*log10(f)-9*log10(b);
case 2 fprintf('El valor de Lmsd es:
disp('Modelo de WALFISH-IKEGAMI') %.2f \n' , Lmsd)
linea=input('El modelo tiene line de Lo=32.4+20*log10(d)+20*log10(f);
vista: ', 's'); fprintf('Las perdidas de
linea1=lower(linea); propagacion en condiciones de espacio libre (Lo): %.2f (dB)
if (linea1=='si') \n' , Lo)
Lo=42.6+26*log10(d)+20*log10(f);
fprintf('Las perdidas de propagacion if((Lrts+Lmsd)>0)
son: %.2f (dB)\n' , Lo); Lb=Lo+Lrts+Lmsd;
elseif (linea1=='no') elseif (Lrts+Lmsd<=0)
%%lori Lb=Lo
resp=input('Sabe el angulo: ', 's'); end
linea2=lower(resp);
if (linea2=='si') fprintf('Las perdidas basicas de
angulo=input('Ingrese el angulo propagacion para el escenario NLOS (Lb): %.2f (dB) \n' , Lb)
en grados: ');
if (angulo >=0 && angulo <=35) PTX=PRX-GRX-GTX+Lb;
lori=10+0.35*angulo; fprintf('La Potencia de
fprintf('Lori es igual a: Transmision (PTX): %.2f \n' , PTX)
%.2f grados \n' , lori); elseif (edificio2=='no')
elseif (angulo>35 && angulo b=input('Ingrese la separacion
<=55) promedio entre edificios (m): ');
lori=2.5+0.075*(angulo-35); p=input('Ingrese el numero de
fprintf('Lori es igual a: pisos (m): ');
%.2f \n' , lori);
elseif (angulo>55 && angulo altura_techo=input('El Techo es
<=90) inclinado: ','s');
lori=40-0.114*(angulo-55);
fprintf('Lori es igual a: altura_techo1=lower(altura_techo);
%.2f \n' , lori); if (altura_techo1=='si')
else m=3;
disp('Ingresar angulo de 0 a
90 grados') elseif (altura_techo1=='no')
end m=0;
end
elseif (linea2=='no') hr=3*p+m;
angulo=90; fprintf('La altura del edificio
lori=40-0.114*(angulo-55); es igual a: %.2f (m) \n' , hr);
fprintf('Lori es igual a: %.2f w=b/2;
\n' , lori); Var_hm=(hr-hm);
end Var_hb=(hb-hr);
%%hr Lrts=-8.2-
edificio=input('Sabe la altura del 10*log10(w)+10*log10(f)+20*log10(Var_hm)+lori;
edificio: ', 's'); if(Lrts<0)
edificio2=lower(edificio); Lrts=0;
end
if (edificio2=='si')
hr=input('Ingrese la altura del fprintf('Perdidas de propagacion
edificio : '); en espacio libre (Lrts): %.2f (dB) \n' , Lrts)
b=input('Ingrese la separacion
promedio entre edificios (m): ') if (hb>hr)
w=b/2; Lbsh=-18*log10(1+Var_hb);
Var_hm=(hr-hm); elseif(hb<=hr)
Var_hb=(hb-hr); Lbsh=0;
Lrts=-8.2- end
10*log10(w)+10*log10(f)+20*log10(Var_hm)+lori; fprintf('Lbsh es igual a: %.2f
if(Lrts<0) (dB) \n' , Lbsh)
Lrts=0;
UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

if(hb>hr)
ka=54; PTX=PRX-GRX-GTX+Lb;
elseif (hb<=hr &&d >=0.5) fprintf('La Potencia de Transmision
ka=54-0.8*Var_hb; (PTX): %.2f \n' , PTX)
elseif (hb<=hr && d<0.5)
ka=(54-(0.8*Var_hb)*(d/5)); case 6
end disp('Modelo de punto de ruptura
múltiple de Ericsson')
if (hb>hr) a0=36.2;
kd=18; a1=30.2;
elseif (hb<=hr) a2=-12;
kd=18-15*((Var_hb)/hr); a3=0.1;
hb=input('Ingrese la Altura del edificio
end (m): ');
while(hb<0)
disp('1 ciudad o entorno disp('Ingrese una altura del
suburbano') edificio mayor a cero ')
disp('2 centro metropolitano') hb=input('Ingrese la Altura del
numero=input('Ingrese una edificio (m): ')
opcion: '); end
switch numero hm=input('Ingrese la Altura del movil
case 1 (m): ');
kf=-4+0.7*((f/925)-1); while(hm<0)
case 2 disp('Ingrese la Altura del receptor
kf=-4+1.5*((f/925)-1); mayor a cero ');
end hm=input('Ingrese la Altura del
fprintf('ka es igual a: %.2f \n' receptor (m): ');
, ka) end
fprintf('kd es igual a: %.2f \n' d=input('Ingrese la distancia entre Tx y
, kd) RX (Km): ');
fprintf('kf es igual a: %.2f \n' while(d<0)
, kf) disp('Ingrese la distancia entre Tx
y RX mayor a cero ');
d=input('Ingrese la distancia entre
Lmsd=Lbsh+ka+kd*log10(d)+kf*log10(f)-9*log10(b); Tx y RX (Km) ');
fprintf('El valor de Lmsd es: end
%.2f \n' , Lmsd) g=44.49*log10(f)-4.78*(log10(f))^2
Lo=32.4+20*log10(d)+20*log10(f);
fprintf('Las perdidas de PL=a0+a1*log10(d)+a2*log10(hb)+a3*log10(hb)*log10(d)-
propagacion en condiciones de espacio libre (Lo): %.2f (dB) 3.2*(log10(11.75))^2+g
\n' , Lo)
fprintf('Las perdida basica de
if((Lrts+Lmsd)>0) propagacion son: %.2f DB \n' , PL)
Lb=Lo+Lrts+Lmsd;
elseif (Lrts+Lmsd<=0) PTX=PRX-GRX-GTX+PL;
Lb=Lo fprintf('La Potencia de Transmision
end (PTX): %.2f \n' , PTX)
case 7
fprintf('Las perdidas basicas de disp('Modelo de OKUMURA')
propagacion para el escenario NLOS (Lb): %.2f (dB) \n' , Lb)
hte=input('Ingrese la Altura del
PTX=PRX-GRX-GTX+Lb; transmisor (m): ');
fprintf('La Potencia de while(hte<0)
Transmision (PTX): %.2f \n' , PTX) disp('Ingrese la Altura del
end transmisor mayor a cero ')
end hte=input('Ingrese la Altura del
case 5 transmisor (m): ')
disp('Modelo de WALFISH-BERTONI') end
hb=input('Ingrese la Altura del edificio hre=input('Ingrese la Altura del
(m): '); receptor (m): ');
while(hb<0) while(hre<0)
disp('Ingrese una altura del disp('Ingrese la Altura del receptor
edificio mayor a cero ') mayor a cero ')
hb=input('Ingrese la Altura del hre=input('Ingrese la Altura del
edificio (m): ') recptor (m): ')
end end
hr=input('Ingrese la Altura del receptor
(m): '); d=input('Ingrese la distancia entre Tx y
while(hr<0) RX (Km): ');
disp('Ingrese la Altura del receptor while(d<0)
mayor a cero '); disp('Ingrese la distancia entre Tx
hr=input('Ingrese la Altura del y RX mayor a cero ');
receptor (m): '); d=input('Ingrese la distancia entre
end Tx y RX (Km) ');
d=input('Ingrese la distancia entre Tx y end
RX (Km): '); Amu=input('Ingrese la atenuacion
while(d<0) relativa promedio (Amu): ');
disp('Ingrese la distancia entre Tx Garea=input('Ganancia debido al tipo de
y RX mayor a cero '); ambiente(Garae):');
d=input('Ingrese la distancia entre Lo=32,4+20*log10(d)+20*log10(f);
Tx y RX (Km) '); if(30<hte && hte<100)
end Ghte=20*log10(hte/200);
b=input('Ingrese la separacion promedio end
entre edificios (m): '); if(hte<3)
while(b<0) Ghte=10*log10(hte/3);
disp('Ingrese la separacion promedio end
entre edificios mayor a cero '); if(3<hte && hte<10)
b=input('Ingrese la separacion Ghte=20*log10(hte/3);
promedio entre edificios (m): '); end
end
H=input('Altura promedio de la antena if(30<hre && hre<100)
con respecto de la altura de los edificios (m):'); Ghre=20*log10(hre/200);
while(H<0) end
disp('Ingrese Altura promedio de la if(hre<3)
antena con respecto de la altura de los edificios mayor a Ghre=10*log10(hre/3);
cero '); end
H=input('Altura promedio de la if(3<hre && hre<10)
antena con respecto de la altura de los edificios (m) '); Ghre=20*log10(hre/3);
end end
A=5*log10((b/2)^2+(hb-hr)^2)- Lp=Lo+Amu-Ghre-Ghte-Garea;
9*log10(b)+20*log(atan(2*(hb-hr))/b);
Lp=89.55+21*log10(f)+38*log(d)- fprintf('Las perdida basica de
18*log10(H)+A-18*log10(1-d^2/17*H); propagacion son: %.2f DB \n' , Lp)
fprintf('Las perdida de trayectoria son: PTX=PRX-GRX-GTX+Lp;
%.2f DB \n' , Lp)
UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

fprintf('La Potencia de Transmision

(PTX): %.2f \n' , PTX)
En la figura 4 nos muestra los tipos de áreas que se
puede escoger, en este caso se escogió un área rural.
case 3
disp('Modelo de 2 RAYOS')
PTX=input('Ingrese la Potencia de Transmision en
(Watts): ');
Gananciatx=10^(GTX/10);
Gananciarx=10^(GRX/10);

PRX=(PTX*Gananciatx*Gananciarx*((hr)^2)*((ht)^2))/(d)^4
fprintf('La Potencia de Recepcion (RTX): %.2f Watts
\n' , PRX)
Lp=40*log10(d)-(10+10+20*log10(ht)+20*log10(hr));
fprintf('Las perdida de trayectoria son: %.2f DB \n'
, Lp)

case 4
disp('Modelo HUYGENS-FRESNEL')
lambda=3e8/f;
V=h*(sqrt((2*(d1+d2))/lambda*(d1*d2)));
if(-0.8< V && V<0)
Lv=-20*log10(0.5-0.62*V);
end

if(0< V && V<1)

Lv=-20*log10(0.5*exp(-0.95*V));
end
if(1< V && V<2.4)
Lv=-20*log10(0.4-(0.1184-(0.38-0.1*V)^2)^1/2);
end
if(V >2.4)
Lv=-20*log10(0.225/V);

end
fprintf('V es igual a: %.2f \n' , V)
fprintf('LV es igual a: %.2f \n' , Lv)
Fig 5. Resultado para el modelo Okumura
end
En la figura 5 se obtiene le resultado de la perdida y de
la potencia de transmisión.
5. ANALISIS Y RESULTADOS
Para el análisis de los resultados se procedió a correr
el programa.

Fig 6. Datos para el modelo COST231

En la figura 6, se ingresan datos con la intención de

que el software indique que se analizara por el modelo
COST231.

Fig 3. Datos para modelos

La figura 3 muestra cómo están puestos los datos y

de acuerdo a esos datos se escogerá el tipo de
modelo.

Fig 7. Selección del modelo COST231

En la figura 7, la simulación dice que se tiene que

trabajar con el modelo COST231

Fig 4. El software escoge el modelo Okumura

UNIVERSIDAD POLITÉCNICA SALESIANA REDES INALÁMBRICAS
Fecha: 30/04/2019

 Se observó claramente que los resultados

realizados a mano en clase, coinciden con los
resultados de la simulación, claro que por el
programa los decimales varían un poco, pero
los resultados son los esperados.

 Se determinó que el mejor modelo para

calcular las perdidas y la potencia es el del
modelo Walfish-Ikegami (COST 231), puesto
que este dispone de más parámetros y
condiciones de la misma manera que se tiene
en el entorno de la vida real.
Fig 8. Selección del ambiente para el modelo
COST231  Al momento de realizar el programa se tuvo
mucho cuidado al rato de realizar cada una
En la figura 8, se escoge el ambiente para el modelo de las validaciones para ello primero se tiene
COST231, en este caso se escogió que haya línea de que conocer cada uno de los parámetros y
vista. rangos que tiene cada modelo aprendidos en
Y a continuación el programa dará el resultado de las clase, también para tener valores coherentes
pérdidas de propagación. en cada variable almacenada y evitar errores
por caracteres no validados.

8. REFERENCIAS
[1] Y. Okumura, E. Ohmori, T. Kawano and K. Fukuda,
"Field strength and its variability in the VHF and UHF
land mobile radio service", Review of the Electrical
Communications Laboratories, vol. 16, no. 9/10, pp.
825-73, Sep 1968.

[2] M. Hata, "Empirical formula for propagation loss in

land mobile radio services," in IEEE Transactions on
Vehicular Technology, vol. 29, no. 3, pp. 317-325, Aug
1980. doi: 10.1109/T-VT.1980.23859

Fig 9. Modelo COST231sin línea de vista [3] J. D. Parsons, "The Mobile Radio Propagation
Channel", New York: John Wiley & Sons, 1992.
En la figura 9, se trabaja sin línea de vista obteniendo
los valores vistos observados en la figura.
Y se obtiene el resultado en este caso la perdida y la
potencia de transmisión.

Fig 10. Rangos en la cual tienen 3 modelos en comun

Con los resultados obtenidos tanto en el programa

como los realizados a mano se obtuvo un 0.423 error
aplicados en los diferentes modelos, siendo un valor
muy aceptable, y comprobando que el programa
realizado en Matlab está bien hecho.

6. CONCLUSIONES