Determining Macro Cell Coverage Based on Okumura Hata in Matlab

Note

This is English translated of my undergraduate assignment in Wireless Network course at the Department of Electrical Engineering, Faculty of Engineering, Udayana University. This assignment has never been published anywhere and I, as the author and copyright holder, license this assignment as customized CC-BY-SA where anyone can share, copy, republish, and sell on condition to state my name as the author and notify that the original and open version available here.

Chapter 1 Introduction

1.1 Background

The signal strength received on the MS (Mobile Station) from the BS (Base Station) is generally affected at 3 stages. The first stage is signal generation at the BS (Base Station), the second stage is propagation, and the third stage is the MS (Mobile Station). Until now, there has not been an empirical calculation of signal propagation, which is currently used is prediction based on experiment and research (no theory). One of them is Okumura-Hata where researchers from Japan surround urban, suburban, and rural areas, trying certain distances and measuring the signal obtained by MS (Mobile Station) from BS (Base Station). To measure signal strength using Okumura must provide many tables based on the situation of the area. Even though using Okumura is accurate, it wastes time because we have to check the tables. Hata derived the calculation formula from Okumura but it also took a long time to compute manually. In this paper, calculations will be made from BS (Base Station) to MS (Mobile Station) based on the Hata formula in MATLAB software based on some data. What is calculated is the RX (Received) Power at MS (Mobile Station) in Urban, Suburban, and Rural areas from a distance of 1 km to 20 km.

1.2 Problem

What are the results of the calculation of receiver power in MS (Mobile Station) from transmitter BS (Base Station) in urban, suburban, and rural areas based on the Hata formula using MATLAB software?

1.3 Objective

To make signal propagation calculations from BS (Base Station) to MS (Mobile Station) based on Okumura-Hata in MATLAB.

1.4 Benefit

It is more energy and time efficient to find signal propagation calculations.

1.5 Scope and Limits

  • Calculation at a distance of 1 - 20 Km.
  • Signal at 900 MHz.
  • The tool specifications have been determined.

Chapter 2 Literature Review

2.1 Signal Strength

The unit of signal strength is the watt, which is the unit of power in electricity. To facilitate the calculation of the signal strength from watts it is converted into decibell units (Kolawole, 2002).

Where:

P = Electrical power (watt)2.2 Free Space Loss

Free Space Loss is the signal strength loss in a vacuum, depending on wavelength and distance. This is a loss that must exist (Kolawole, 2002).

Where:

FSL = Free Space Loss

r = the distance between the transmitter and the receiver (meter)

lamda = wavelength (meter)

c = speed of light (meter/second)

2.3 Hata Path Loss

Until now, the empirical formula for calculating the signal strength from BS (Base Station) to MS (Mobile Station) has not been found in the real world environment. What does exist is an approximate prediction based on research. Okumura and Hata surround an area at certain intensities and record the signal strength received by the MS (Mobile Station) from the BS (Base Station) at certain distances. The calculation of the Hata model is as follows (ETSI, 1999):

Limitation:

Frequency (f): 150 - 1000 MHz

Base Station Height (Hb): 30 - 200 m

Mobile Height (Hm): 1 - 10 m

Distance (d): 1 - 20 km

Untuk wilayah metropolitan:

a(Hm)= f ≤ 200 MHz

a(Hm)= f ≥ 400 MHz

For nonmetropolitan areas:

a(Hm)=

Loss for Urban area: 6.55log(Hb)log(d)

Loss for Suburban area

Loss for Rural Quasi-Open area

Loss for Rural Open-Area

2.4 Link Budget

In general, there are 3 steps to calculate the signal strength from the BS (Base Station) to the MS (Mobile Station). The first stage is the signal strength generated by the BS (Base Station) affected by the smear loss and gain that occurs at the BS (Base Station). The second stage is the path loss, it can be calculated with Hata. The third stage is the gain and loss that occurs in MS (Mobile Station). Therefore it can be calculated (ETSI, 1999):

Where:

Pm = Signal strength on MS

Pb = Signal strength on BS

Gb = Gain Atenna BS

Ld = Duplexer Loss BS

Lj = Jumper Loss BS

Lf = Feeder Loss BS

Ltf = Filter Loss BS

Mf = Fade Margin

Ab = Body Atennuation

Av = Vehide Atennuation

Abd = Building Atennuation

Lpj = Path loss

Lm = Feeder Loss MS

Gm = Gain Atenna MS

Lo = Other Loss

2.5 MATLAB

MATLAB is a high-level language and an interactive environment for numerical computation, visualization, and programming. Using MATLAB, you can analyze data, develop algorithms, and create models and applications. Languages, tools, and built-in math functions allow you to explore multiple approaches and reach solutions faster than traditional spreadsheets or programming languages, such as C / C ++ or Java. MATLAB can be used for a wide variety of applications, including signal processing and communications, image and video processing, control systems, test and measurement, computational finance, and computational biology. More than one million engineers and scientists in industry and academia use MATLAB, the language of technical computing (Little, 2013).

Chapter 3 Experimental Method

3.1 Place and Time of Experiment

The experiment was carried out at home on May 22, 2013.

3.2 Tools

3.3 Program

Created 3 programs. The first program to calculate the path loss according to Hata for urban, suburban, and rural areas with distances of 1km, 2km, 3km, 4km, 5km, 6km, 7km, 8km, 9km, 10km, 11km, 12km, 13km, 14km, 15km, 16km, 17km, 18km, 19km, 20km. The second program is product loss, which is the sum of the path loss and indoor propagation loss. The third program for calculating the RX signal strength based on TX, gain, and loss on the device. Can also be made in simulink.

Dalam bentuk simulink
Figure 3.1 In simulink form
%Hata Model%Frequency 150 – 1000 (MHz)
f = 900;

%Base station height 30 – 200 (m)
Hb = 40;

%Mobile Height 1 – 10 (m)
Hm = 1.5;

%Distance 1 – 20 (km)
d = 1:20;

%Urban Area Loss small medium city (dB)
asm = (((1.1*log10(f))-0.7)*Hm)-((1.56*log10(f))-0.8);
Lusm = 69.55 + (26.16*log10(f)) – (13.82*log10(Hb)) – asm + ((44.9 – (6.55*log10(Hb)))*log10(d));

%Urban Area Loss large city (dB)
if f <= 200
    al = (8.29*(log10(1.54*Hm).^2))-1.1;
    Lul = 69.55 + (26.16*log10(f)) – (13.82*log10(Hb)) – al + ((44.9 – (6.55*log10(Hb)))*log10(d));
else if f >= 400
    al = (3.2*(log10(11.75*Hm).^2))-4.97;
    Lul = 69.55 + (26.16*log10(f)) – (13.82*log10(Hb)) – al + ((44.9 – (6.55*log10(Hb)))*log10(d));
else
    ‘frequency range undefine’;
end

%Suburban Area Loss small medium city (dB)
Lsusm = Lusm – (2*(log10(f/28).^2))-5.4;

%Suburban Area large city (dB)
Lsul = Lul – (2*(log10(f/28).^2))-5.4;

%Rural Area small medium city (Quasi-Open) (dB)
Lrqosm = Lusm – (4.78*(log10(f).^2))+(18.33*log10(f))-35.94;

%Rural Area large city (Quasi-Open) (dB)
Lrqol = Lul – (4.78*(log10(f).^2))+(18.33*log10(f))-35.94;

%Rural Area small medium city (Open Area) (dB)
Lrosm = Lusm – (4.78*(log10(f).^2))+(18.33*log10(f))-40.94;

%Rural Area large city (Open Area) (dB)
Lrol = Lul – (4.78*(log10(f).^2))+(18.33*log10(f))-40.94;

subplot(2,1,1);
plot(d,Lusm,d,Lsusm,d,Lrqosm,d,Lrosm);
grid on;
subplot(2,1,2);
plot(d,Lul,d,Lsul,d,Lrqol,d,Lrol);
grid on;

Code 3.1 Path Loss Hata Model

%Product Loss%u = urban %s = surban
%r = rural

%Total Path Loss (dB)
Lptu = Lul;
Lpts = Lsul;
Lptr = Lrol;

%Building Attenuation (dB)
Abdu = 15;
Abds = 12;
Abdr = 0;

%Vhide Attenuation (dB)
Avu = 0;
Avs = 0;
Avr = 0;

%Body Attenuation (dB)
Abu = 2;
Abs = 2;
Abr = 2;

%Fade Margin (dB)
Mfu = 5.6;
Mfs = 5.6;
Mfr = 5.6;

%Feeder loss per meter (dB/m)
Lf_m = 0.0646;

%Feeder loss (dB)
Lf = Hb*Lf_m;

%Product Loss (dB)
Lpu = Lptu + Abdu + Avu + Abu + Mfu + Lf;
Lps = Lpts + Abds + Avs + Abs + Mfs + Lf;
Lpr = Lptr + Abdr + Avr + Abr + Mfr + Lf; 

Code 3.2 Product Loss

%MS_BS RX Power%u = urban %s = surban
%r = rural

%MS TX Power (dBm)
Pm = 30;

%MS Antenna Gain (dBi)
Gm = 2;

%MS Feeder Loss (dB)
Lm = 0;

%BS TX Power (dBm)
Pb = 47;

%BS Antenna Gain (dBi)
Gb = 20;

%BS Diversity Gain (dB)
Gd = 3.5;

%BS Duplexer Loss (dB)
Ld = 0.8;

%BS Jumper/Connector Loss (dB)
Lj = 0.9;

%BS TX Filter Loss (dB)
Ltf = 2.3;

%Other Loss (dB)
Lo = 0;

%BS RX Power
Sbu = Pm + Gm – Lm + Gb + Gd – Ld – Lj – Lpu – Lo;
Sbs = Pm + Gm – Lm + Gb + Gd – Ld – Lj – Lps – Lo;
Sbr = Pm + Gm – Lm + Gb + Gd – Ld – Lj – Lpr – Lo;

%MS RX Power
Smu = Gm – Lm + Pb + Gb – Ld – Lj – Ltf – Lpu – Lo;
Sms = Gm – Lm + Pb + Gb – Ld – Lj – Ltf – Lps – Lo;
Smr = Gm – Lm + Pb + Gb – Ld – Lj – Ltf – Lpr – Lo;

plot(d,Smu,d,Sms,d,Smr);
legend(‘Urban’, ‘Suburban’, ‘Rural’);
title(‘Received Power’);
xlabel(‘Distance (km)’);
ylabel(‘Power (dBm)’);
grid on; 

Code 3.3 BS MS Power

a = 20; b = 20; ang=0:0.01:2*pi;
x1=d(1)*cos(ang);
y1=d(1)*sin(ang);
x2=d(2)*cos(ang);
y2=d(2)*sin(ang);
x3=d(3)*cos(ang);
y3=d(3)*sin(ang);
x4=d(4)*cos(ang);
y4=d(4)*sin(ang);
x5=d(5)*cos(ang);
y5=d(5)*sin(ang);
x6=d(6)*cos(ang);
y6=d(6)*sin(ang);
x7=d(7)*cos(ang);
y7=d(7)*sin(ang);
x8=d(8)*cos(ang);
y8=d(8)*sin(ang);
x9=d(9)*cos(ang);
y9=d(9)*sin(ang);
x10=d(10)*cos(ang);
y10=d(10)*sin(ang);
x11=d(11)*cos(ang);
y11=d(11)*sin(ang);
x12=d(12)*cos(ang);
y12=d(12)*sin(ang);
x13=d(13)*cos(ang);
y13=d(13)*sin(ang);
x14=d(14)*cos(ang);
y14=d(14)*sin(ang);
x15=d(15)*cos(ang);
y15=d(15)*sin(ang);
x16=d(16)*cos(ang);
y16=d(16)*sin(ang);
x17=d(17)*cos(ang);
y17=d(17)*sin(ang);
x18=d(18)*cos(ang);
y18=d(18)*sin(ang);
x19=d(19)*cos(ang);
y19=d(19)*sin(ang);
x20=d(20)*cos(ang);
y20=d(20)*sin(ang);

%plot(a+x1,b+y1,a+x2,b+y2,a+x3,b+y3,a+x4,b+y4,a+x5,b+y5,a+x6,b+y6,a+x7,b+y7,a+x8,b+y8,a+x9,b+y9,a+x10,b+y10,a+x11,b+y11,a+x12,b+y12,a+x13,b+y13,a+x14,b+y14,a+x15,b+y15,a+x16,b+y16,a+x17,b+y17,a+x18,b+y18,a+x19,b+y19,a+x20,b+y20);
fill(a+x20,b+y20,[0.4 0 0],a+x19,b+y19,[0.6 0 0],a+x18,b+y18,[0.9 0 0],a+x17,b+y17,[1 0.2 0],a+x16,b+y16,[1 0.4 0],a+x15,b+y15,[1 0.7 0],…
a+x14,b+y14,[1 0.9 0],a+x13,b+y13,[0.8 1 0],a+x12,b+y12,[0.6 1 0],a+x11,b+y11,[0 1 0.3],a+x10,b+y10,[0 1 0.5],a+x9,b+y9,[0 1 0.7],…
a+x8,b+y8,[0 1 0.9],a+x7,b+y7,[0 1 1],a+x6,b+y6,[0 0.7 1],a+x5,b+y5,[0 0.6 1],a+x4,b+y4,[0 0.4 1],a+x3,b+y3,[0 0.3 1],…
a+x2,b+y2,[0 0.2 1],a+x1,b+y1,[0 0 1]);
colorbar(‘YTickLabel’,{‘0’,Smu(2),Smu(4),Smu(6),Smu(8),Smu(10),Smu(12),Smu(14),Smu(16),Smu(18),Smu(20)});
grid on;
xlabel(‘X-distance (km)’); %# Add an x axis label
ylabel(‘Y-distance (km)’); %# Add a y axis label
subplot(3,1,1); 

Code 3.4 Coverage 1/2

fill(a+x20,b+y20,[0.4 0 0],a+x19,b+y19,[0.6 0 0],a+x18,b+y18,[0.9 0 0],a+x17,b+y17,[1 0.2 0],a+x16,b+y16,[1 0.4 0],a+x15,b+y15,[1 0.7 0],…a+x14,b+y14,[1 0.9 0],a+x13,b+y13,[0.8 1 0],a+x12,b+y12,[0.6 1 0],a+x11,b+y11,[0 1 0.3],a+x10,b+y10,[0 1 0.5],a+x9,b+y9,[0 1 0.7],… a+x8,b+y8,[0 1 0.9],a+x7,b+y7,[0 1 1],a+x6,b+y6,[0 0.7 1],a+x5,b+y5,[0 0.6 1],a+x4,b+y4,[0 0.4 1],a+x3,b+y3,[0 0.3 1],…
a+x2,b+y2,[0 0.2 1],a+x1,b+y1,[0 0 1]);
colorbar(‘YTickLabel’,{‘0’,Smu(4),Smu(8),Smu(15),Smu(17),Smu(20)});
grid on;
xlabel(‘X-distance (km)’); %# Add an x axis label
ylabel(‘Y-distance (km)’); %# Add a y axis label
subplot(3,1,2);

fill(a+x20,b+y20,[0.4 0 0],a+x19,b+y19,[0.6 0 0],a+x18,b+y18,[0.9 0 0],a+x17,b+y17,[1 0.2 0],a+x16,b+y16,[1 0.4 0],a+x15,b+y15,[1 0.7 0],…
a+x14,b+y14,[1 0.9 0],a+x13,b+y13,[0.8 1 0],a+x12,b+y12,[0.6 1 0],a+x11,b+y11,[0 1 0.3],a+x10,b+y10,[0 1 0.5],a+x9,b+y9,[0 1 0.7],…
a+x8,b+y8,[0 1 0.9],a+x7,b+y7,[0 1 1],a+x6,b+y6,[0 0.7 1],a+x5,b+y5,[0 0.6 1],a+x4,b+y4,[0 0.4 1],a+x3,b+y3,[0 0.3 1],…
a+x2,b+y2,[0 0.2 1],a+x1,b+y1,[0 0 1]);
colorbar(‘YTickLabel’,{‘0’,Sms(4),Sms(8),Sms(15),Sms(17),Sms(20)});
grid on;
xlabel(‘X-distance (km)’); %# Add an x axis label
ylabel(‘Y-distance (km)’); %# Add a y axis label
subplot(3,1,3);

fill(a+x20,b+y20,[0.4 0 0],a+x19,b+y19,[0.6 0 0],a+x18,b+y18,[0.9 0 0],a+x17,b+y17,[1 0.2 0],a+x16,b+y16,[1 0.4 0],a+x15,b+y15,[1 0.7 0],…
a+x14,b+y14,[1 0.9 0],a+x13,b+y13,[0.8 1 0],a+x12,b+y12,[0.6 1 0],a+x11,b+y11,[0 1 0.3],a+x10,b+y10,[0 1 0.5],a+x9,b+y9,[0 1 0.7],…
a+x8,b+y8,[0 1 0.9],a+x7,b+y7,[0 1 1],a+x6,b+y6,[0 0.7 1],a+x5,b+y5,[0 0.6 1],a+x4,b+y4,[0 0.4 1],a+x3,b+y3,[0 0.3 1],…
a+x2,b+y2,[0 0.2 1],a+x1,b+y1,[0 0 1]);
colorbar(‘YTickLabel’,{‘0’,Smr(4),Smr(8),Smr(15),Smr(17),Smr(20)});
grid on;
xlabel(‘X-distance (km)’); %# Add an x axis label
ylabel(‘Y-distance (km)’); %# Add a y axis label

Code 3.5 Coverage 2/2

3.4 Input

Table 3.1 Tools used
NO Tools Spesification Software
1 PC Intel® Pentium® core-i5 processor 4 GB RAM 2 GB VGA Operating System Windows 7 Ultimate 2 MATLAB 2010
Table 3.2 Input 1
Frequency (f) 900 MHz
Base Station Height (Hb) 40 m
Mobile Station Height (Hm) 1.5 m
Distance (d) 1:20 km
BS Power (dBm) 47
BS Antenna Gain (dBi) 20
MS Antenna Gain (dBi) 2
BS Duplexer Loss (dB) 0.8
BS Jumper/Connector Loss (dB) 0.9
BS Filter Loss (dB) 2.3
Table 3.3 Input 2
Description Urban Suburban Rural
Fade Margin (dB) 5.6 5.6 5.6
Feeder Loss / Meter (dB/m) 0.0646 0.0646 0.0646
Body Atennuation (dB) 2 2 2
Building Atennuation (dB) 15 12 0

3.5 Output

Grafik jarak terhadap MS power
Figure 3.2 Graph of distance to MS power
MS power pada jarak 1:20 km
Figure 3.3 MS power at a distance of 1:20 km

Chapter 4 Discussion

4.1 Signal Strength in MS (Mobile Station)

It can be seen in Figure 3.1 that the signal strength decreases exponentially with the distance between the MS (Mobile Station) and the BS (Base Station). According to Hata and Free Space Loss, the signal loss will get bigger as the distance value increases (d). A clearer graphic can be seen in Figure 4.1.

Grafik jarak terhadap MS (Mobile Station) diperbesar
Figure 4.1 Enlarged distance graph to MS (Mobile Station)

From the graph above, it can be seen that the Hata formula shows that urban MS (Mobile Station) has the lowest signal strength value compared to suburban and rural areas. Meanwhile, rural areas have the highest signal strength. If using the Hata formula at a distance of 20 km, the Mobile must have a minimum sensitivity of -130 dBm in urban areas, -117 dBm in suburban areas, and -86 dBm in rural areas. Using MATLAB can also illustrate the value of the MS signal strength in the form of a map at a certain radius as follows with the BS in the middle:

MS power pada jarak 1:20 km diperbesar
Figure 4.2 MS power at a distance of 1:20 km is enlarged

The graph and image above are the values:

6.55log(Hb))log(d)

With input variations d = 1 km to d = 20 km. The formula above is Loss according to the Hata Model. The formula is based on the ETSI 1999 standard.

4.2 Increase MS Power

4.2.1 Minimizes Hata Path Loss

At first glance the Hata formula does not appear whether the frequency or high parameters of the BTS should be increased or decreased. The distance and height variables MS are natural so they cannot be adjusted. If you use MATLAB, you can see the effect of frequency and height of BTS on loss (with d = 20 km and Hm = 1.5 m).

Grafik Path Loss terhadap Base Station Height
Figure 4.3 Graph of Path Loss to Base Station Height

From the graph above, it can be seen that the Path Loss will decrease when the BS (Base Station) gets higher. Therefore increasing the BTS height is one solution to reduce loss. The graph above shows the values:

6.55log(Hb))log(d)

With a variation of input Hb = 30 m to d = 200 m. The formula above is Loss according to the Hata Model. The formula is based on the ETSI 1999 standard.

Grafik Path Loss terhadap Frekuensi
Figure 4.4 Graph of Path Loss to Frequency

From the chart above, it can be seen that the Path Loss will increase the higher the frequency. The graph above shows the values:

6.55log(Hb))log(d)

With input variations f = 150 MHz to f = 1000 MHz. The formula above is Loss according to the Hata Model. The formula is based on the ETSI 1999 standard.

4.2.2 Increasing Tool Specifications

Based on the link budget formula from ETSI 1999:

A higher MS (Mobile Station) signal strength value can be achieved by increasing antenna gain, increasing BS (Base Station) transmit power, and minimizing device loss.

Chapter 5 Closing

5.1 Conclusion

Using MATLAB is much faster than manual counting and can graph and illustrate a BS (Base Station) radius illustration with MS (Mobile Station). From the MATLAB results at a distance of 1 km, the signal strength of MS is -85 dBm in urban areas, -72 dBm in suburban areas, and -8-41 dBm in rural areas. At the farthest distance of 20 km, the signal strength of MS is -130 dBm in urban areas, -117 dBm in suburban areas, and -86 dBm in rural areas.

Bibliography

  • ETSI. 1999. ETSI TR 101 362 V8.3.0 (2000-04). Sophia Antipolis Cedex : France.
  • Kolawole, M. 2002. Satellite Communication Engineering. Marced Deker, inc : New York
  • Little, J. 2013. http://www.mathworks.com/products/matlab. Access 27 Mei 2013.

Mirror

Comments