PROPAGACIÓN DE ONDA
UTPL
UNIVERSIDAD TÉCNICA PARTICULAR DE LOJA
TALLER 01 (2B)
Integrantes: Jennyfer Quezada, Freddy Jiménez
Fecha: 9 de junio de 2021
Paralelo: “A”
ATENUACIÓN POR VEGETACIÓN Y LLUVIA
Calcular la atenuación excesiva por vegetación cuando una de las terminales se
encuentra en el interior de una zona boscosa, si se está trabajando en polarización
oblicua, en la frecuencia de 1852.2 MHz, cuando el dispositivo está dentro de la zona
boscosa 15 metros. Si el PIRE es 30 dBm. Calcular la potencia de recepción en el
receptor si tiene una antena de 1 dBi. La distancia entre el receptor y la estación base
es de 250 metros.
𝑑 = 15 𝑚𝑒𝑡𝑟𝑜𝑠
𝐷 = 250 𝑚𝑒𝑡𝑟𝑜𝑠
Datos:
• Polarización oblicua
o 𝑓 = 1852.2𝑀𝐻𝑧
o 𝛾 = 0.30 (𝑑𝐵⁄𝑚)
o 𝐴𝑚 = 29(𝑑𝐵)
• 𝑑 (𝑑𝑒𝑛𝑡𝑟𝑜 𝑑𝑒 𝑙𝑎 𝑧𝑜𝑛𝑎) = 15𝑚
• 𝑃𝐼𝑅𝐸 = 30𝑑𝐵𝑚
• 𝐺𝑟 = 1𝑑𝐵𝑖
• 𝐷(𝑟𝑒𝑐𝑒𝑝𝑡𝑜𝑟 − 𝑒𝑠𝑡𝑎𝑐𝑖ó𝑛 𝑏𝑎𝑠𝑒) = 250𝑚
a) Calcular la atenuación excesiva por vegetación cuando una de las terminales
se encuentra en el interior de una zona boscosa
𝑑𝛾
−
𝐴𝑒𝑣 = [1 − 𝑒 𝐴𝑚 ] → Atenuación excesiva debido a la presencia de vegetación
(15𝑚)∗0.30(𝑑𝐵⁄𝑚)
−
𝐴𝑒𝑣 = [1 − 𝑒 29(𝑑𝐵) ]
𝐴𝑒𝑣 = 4.16824
b) Calcular la potencia de recepción en el receptor si tiene una antena de 1 dBi.
𝐷 = 250 𝑚𝑒𝑡𝑟𝑜𝑠
𝑃𝑟 = 𝑃𝐼𝑅𝐸 − 𝑳𝒃 + 𝐺𝑟 + 𝐴𝑒𝑣 → Potencia recibida a la entrada del receptor (dBm)
𝐿𝑏𝑓 = 32.45 + 20𝑙𝑜𝑔 (𝑓𝑀𝐻𝑧 ∗ 𝐷𝑘𝑚 ) → Pérdida básica de espacio libre (dB)
𝑳𝒃 = 𝐿𝑏𝑓 → Pérdida básica de propagación (dB)
Por lo tanto:
𝐿𝑏𝑓 = 32.45 + 20𝑙𝑜𝑔 (1852.2 𝑀𝐻𝑧 ∗ 0.25 𝑘𝑚) = 85.7626𝑑𝐵
𝐿𝑏 = 85.7626𝑑𝐵
𝑃𝑟 = 30𝑑𝐵𝑚 − 85.7626𝑑𝐵 + 1𝑑𝐵𝑖 + 4.16824
𝑃𝑟 = −50.59436 𝑑𝐵𝑚
Calcular la atenuación 𝐴0,01 por precipitación en el 0.01% del tiempo, si se tiene el
siguiente escenario:
𝜃 = 10°
𝑑 = 45 𝑘𝑚
Datos:
• Polarización vertical
o 𝑓 = 18𝐺𝐻𝑧
o 𝑘𝐻 = 0.07078
o 𝛼𝐻 = 1.0818
o 𝑘𝑉 = 0.07708
o 𝛼𝑉 = 1.0025
o 𝜏 = 90°
• 𝑑 = 45 𝑘𝑚
• 𝑅0,01 = 100 𝑚𝑚⁄ℎ
a) Obtención de la intensidad de la lluvia R
𝑅0,01 = 100 𝑚𝑚⁄ℎ
b) Cálculo de atenuación especifica por hidrometeoros
𝛾𝑅 = 𝒌 ∗ 𝑅0,01 𝛼 → Atenuación especifica
𝑘𝐻 +𝑘𝑉 + (𝑘𝐻 −𝑘𝑉 ) 𝑐𝑜𝑠 𝜃 2 ∗ 𝑐𝑜𝑠 2𝜏
𝒌=
2
𝑘𝐻 𝛼𝐻 +𝑘𝑉 𝛼𝑉 + (𝑘𝐻 𝛼𝐻 −𝑘𝑉 𝛼𝑉 ) 𝑐𝑜𝑠 𝜃 2 ∗ 𝑐𝑜𝑠 2𝜏
𝛼=
2∗𝑘
Por lo tanto:
0.07078 + 0.07708 + (0.07078 − 0.07708) 𝑐𝑜𝑠(10)2 ∗ 𝑐𝑜𝑠 2(90)
𝑘=
2
𝑘 = 0.07698
0.07078(1.0818) + 0.07708(1.0025) + (0.07078(1.0818) − 0.07708(1.0025)) 𝑐𝑜𝑠(10)2 ∗ 𝑐𝑜𝑠 2(90)
𝛼=
2(0.07698)
𝛼 = 1.0036
𝛾𝑅 = 𝑘 ∗ 𝑅0,01 𝛼
𝛾𝑅 = 𝟎. 𝟎𝟕𝟔𝟗𝟖 ∗ 100 (1.0036)
𝛾𝑅 = 𝟕. 𝟖𝟐𝟔𝟔 (𝒅𝑩⁄𝒌𝒎)
c) Cálculo distancia efectiva
𝑑𝑒𝑓𝑓 = 𝑟 ∗ 𝑑 → Distancia efectiva
1
𝑟= 0.073∗𝛼
0.477𝑑 0.633 ∗ 𝑅0,01 ∗ 𝑓 0.123 − 10.579(1 − 𝑒𝑥𝑝(−0.024 ∗ 𝑑 ))
1
𝑟= 0.073∗(1.0036)
0.477(45𝑘𝑚)0.633 ∗ 100 𝑚𝑚⁄ℎ ∗ (18𝐺𝐻𝑧)0.123 − 10.579(1 − 𝑒𝑥𝑝(−0.024 ∗ (45𝑘𝑚)))
𝑟 = 0.2755
El valor máximo recomendado de r es 2.5, de manera que, si el denominador de la
ecuación es menor de 0.4, utilizar 𝑟 = 2.5
0.2755 < 0.4
𝑟 = 2.5
Por lo tanto:
𝑑𝑒𝑓𝑓 = 2.5 ∗ (45𝑘𝑚)
𝑑𝑒𝑓𝑓 = 112.5 𝑘𝑚
d) Cálculo de atenuación del trayecto excedida durante el 0.01% del tiempo
𝐴0,01 = 𝛾𝑅 ∗ 𝑑𝑒𝑓𝑓
𝐴0,01 = 7.8266 (𝑑𝐵⁄𝑘𝑚) ∗ 112.5 𝑘𝑚
𝐴0,01 = 880.4925 𝑑𝐵
Implementar un código en MATLAB u OCTAVE, que permita calcular la A0.01, en
cualquier escenario, deberá permitir ingresar parámetros como frecuencia,
intensidad de lluvia, distancia del enlace, tipo de polarización, e inclinación del
trayecto.
Código en MATLAB
%%Universidad Técnica Particular de Loja
%%Electrónica y Telecomunicaciones
%%Propagación de Ondas
%%Integrantes: Freddy Jiménez y Jennyfer Quezada
%%Tema: Cálculo de la atenuación excesiva
%-------------------------------------------
clc;
clear all;
disp("Cálculo de la atenuacíon excesiva")
disp(" ")
%%Ingresar Datos
disp("*** ESCENARIOS ***")
disp(" ")
disp(" Atenuación por presencia de vegetación --> 1")
disp(" Atenuación por gases atmosféricos --> 2")
disp(" Atenuación por hidrometeoros o precipitación --> 3")
disp(" ")
texto = 'Elija el esecenario: ';
escenario = input(texto); %ingresar escenario
%%Calcular la atenuación en escenario elegido
if escenario==1
disp(" ")
disp("*** Atenuación por presencia de vegetación ***")
disp(" ")
texto2 = 'Ingrese la frecuencia (MHz): ';
f= input(texto2); %frecuencia
texto3 = 'Ingrese la distancia del trayecto dentro de la zona boscosa (m): ';
d= input(texto3); %distancia
parametro=[105.9 466.475 949.0 1852.2 2117.5; 0.04 0.12 0.17 0.30 0.34; 9.4
18.0 26.5 29.0 34.1];
A=parametro(1,:);
[j,k]= min(abs(A-f)); %encontrar la frecuencia más cercana
B=parametro(2,:);
Y=B(k); %atenuación específica
C=parametro(3,:);
Am=C(k); %atenuación adicional total máxima
%Atenuación excesiva (Aev)
Aev = Am*(1-exp(-((d*Y)/(Am))));
disp("Variables:")
disp(" y= "+y+" dB/m, Am= "+Am+" dB")
disp(" ")
disp("--La atenuación excesiva es: "+Aev+" dB--")
else if escenario==2
disp(" ")
disp("*** Atenuación por gases atmosféricos ***")
disp(" ")
texto1 = 'Ingrese la atenuación del vapor de agua: ';
Yw= input(texto1); %atenuación debido al vapor de agua
texto2 = 'Ingrese la atenuación del oxígeno: ';
Yo= input(texto2); %atenuación debido al oxígeno
texto3 = 'Ingrese la distancia (km): ';
d= input(texto3); %distancia
%Atenuación excesiva (Lm)
Lm=(Yw+Yo)*d;
disp("--La atenuación excesiva es: "+Lm+" dB--")
else if escenario==3
disp(" ")
disp("*** Atenuación por hidrometeoros o precipitación ***")
disp(" ")
texto1 = 'Ingrese la frecuencia (entre 1 y 20)GHz: ';
f= input(texto1); %frecuencia
texto2 = 'Ingrese el ángulo de elevación del trayecto: ';
angulo= input(texto2); %elevación del trayecto
texto3 = 'Ingrese la intensidad de lluvia (mm/h): ';
R= input(texto3); %intensidad de lluvia
disp(" ")
disp("Polarización horizontal --> 1")
disp("Polarización vertical --> 2")
disp("Polarización circular --> 3")
disp(" ")
texto4 = 'Ingrese el tipo de polarización: ';
P= input(texto4); %polarización
texto5 = 'Ingrese la distancia del enlace (km): ';
d= input(texto5); %distancia
%Valor de tao
if P==1
t=2*pi;
else if P==2
t=pi/2;
else if P==3
t=pi/4;
end
end
end
%Parámetros P.838
frecuencias=[1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 7 8 9 10 11 12 13 14 15 16
17 18 19 20];
factorkH=[0.0000259 0.0000443 0.0000847 0.0001321 0.0001390 0.0001155
0.0001071 0.0001340 0.0002162 0.0003909 0.0007056 0.001915 0.004115 0.007535
0.01217 0.01772 0.02386 0.03041 0.03738 0.04481 0.05282 0.06146 0.07078 0.08084
0.09164];
factoraH=[0.9691 1.0185 1.0664 1.1209 1.2322 1.4189 1.6009 1.6948
1.6969 1.6499 1.5900 1.4810 1.3905 1.3155 1.2571 1.2140 1.1825 1.1586 1.1396 1.1233
1.1086 1.0949 1.0818 1.0691 1.0568];
factorkV=[0.0000308 0.0000574 0.0000998 0.0001464 0.0001942 0.0002346
0.0002461 0.0002347 0.0002428 0.0003115 0.0004878 0.001425 0.003450 0.006691
0.01129 0.01731 0.02455 0.03266 0.04126 0.05008 0.05899 0.06797 0.07708 0.08642
0.09611];
factoraV=[0.8592 0.8957 0.9490 1.0085 1.0688 1.1387 1.2476 1.3987
1.5317 1.5882 1.5728 1.4745 1.3797 1.2895 1.2156 1.1617 1.1216 1.0901 1.0646 1.0440
1.0273 1.0137 1.0025 0.9930 0.9847];
[j,k]= min(abs(frecuencias-f));
kH=factorkH(k);
aH=factoraH(k);
kV=factorkV(k);
aV=factoraV(k);
%Atenuación específica (yR)
K=(kH+kV+(kH-kV)*(cos((angulo*(pi/180))^2))*(cos(2*t)))/2;
a=((kH*aH)+(kV*aV)+((kH*aH)-
(kV*aV))*(cos((angulo*(pi/180))^2))*(cos(2*t)))/(2*K);
yR=(K*R^a);
%Distancia efectica (deff)
den=((0.477*d^0.633)*(R^(0.073*a))*(f^0.123)-(10.579*(1-exp(-
0.024*d))))^(-1);
if den<0.4
r=2.5;
else
r=(1/den);
end
deff=r*d;
%Atenuación excedida
A=yR*deff;
disp("Variables:")
disp(" kH= "+kH+", aH= "+aH+", kV= "+kV+", aV= "+aV)
disp(" k= "+K+", a= "+a)
disp(" k= "+K+", a= "+a+" r= "+r)
disp(" yR= "+yR+" dB/km")
disp(" deff= "+deff+" km")
disp(" ")
disp("--La atenuación excedida durante el 0.01% del trayecto es: "+A+"
dB--")
end
end
end
RESULTADOS
• Escenarios
1. Cálculo de la atenuación por presencia de vegetación.
2. Cálculo de la atenuación excesiva por gases atmosféricos.
3. Cálculo de la atenuación excesiva por hidrometeoros o precipitación con una
polarización vertical.
