1. Introduction
These days the world is experiencing global warming, and one of the ways of tackling this issue is to create smart cities and autonomous vehicles. A smart city is a city equipped with different types of active and passive sensors communicating with each other. Cyber and signal security are tremendously important to keep a smart city safe and prevent any cyber and signal attacks. This is even more important for autonomous vehicles and their navigation. Therefore, a smart city should be equipped well with security networks for the localisation of any signal interference device. In this article, the design of such a security network based on the received signal strength difference (RSSD) is presented.
Two known types of signal interference are jamming and spoofing. The former means to transmit a signal into the same or near band as the satellite navigation band to disable navigation and spoofing stands for transmitting a fake signal [
1]. Studies have shown that a simple and cheap spoofer can overtake for e.g., a ship navigation without being detected; see [
2,
3]. Since the power level of the global navigation satellite systems (GNSS) signal is low, therefore, a weak interference signal can jam a receiver [
1], see some real examples in [
4,
5,
6,
7,
8,
9,
10,
11].
Localisation of an interferer means estimating its coordinates in a pre-defined coordinate system. Some points with known coordinates, which are so-called anchor nodes or monitoring stations (MSs) have some sensors to receive information from the interferer. If the sensors can detect an interfering signal and measure its time of arrival (TOA), angle of arrival (AOA) or time difference of arrivals (TDOA) the interfering device can be localised. Numerous studies have been conducted for localisation based on these observables, e.g., see [
12,
13,
14].
Another type of information, which can be received from the interferer is known as pathloss information such as received signal strength (RSS), RSSD, or received signal voltage. Thompson et al. [
14] tested an outdoor localisation of a Wifi with unknown transmission power based on the RSS and concluded that in their 140 m
2 test area the source can be localised with a root mean squares error of about 37 m. Thompson et al. [
15] used automatic gain control [
16] for the detection and estimation of jammer to noise of an interference and concluded that their method results in the capability of detecting jammer-to-noise ratios down to −8 dB. Wang and Inkol [
17] developed an efficient algorithm for solving the nonlinear problem of localisation using the RSSD. Thompson et al. [
18] developed the localisation procedure based on the RSS in the presence of ground reflection. Thompson [
19] investigated the interference device detection and localisation by analysing the dilution of precision (DOP), from the RSS, TDOA and concluded that the TDOA is superior to the RSS. Bekcibasi and Tenruh [
20] presented a method for increasing the accuracy of localisation based on the RSS and named it dynamic distance reference anchor method. Egenbråten [
21] studied radio frequency emitter localisation based on the power difference of arrivals. Robertson et al. [
22] performed a similar study on a Network of GNU Radio Sensors. Hu and Leus [
23] presented a robust localisation method based on the RSSD and the semidefinite programming. Nyström [
24] also studied the localisation of the GNSS interference using the RSSD. Zhou et al. [
25] studied the error propagation for both mobile and target tracking and sensor node location calibration. Localisation using the RSS was also studied by Niu et al. [
26] with focus on practical issues and solutions. Pu and You [
27] proposed a location fingerprinting algorithm based on the general and weighted k-nearest neighbour algorithms to estimate the position of the target node using the RSS with an estimated position error of about 1.8 m. Li et al. [
28] used the orientation compensation method and the RSSD for localisation and mentioned that their approach was effective in enhancing localisation that had multiple devices with three-dimensional orientation diversity. Wu et al. [
29] presented a positioning approach based on linear regression of the AOA and the RSSD and showed that the proposed algorithm can achieve better accuracy than existing Radio Frequency Identification (RFID) positioning approaches. Xu and Dogancay [
30] suggested a method for estimating the optimal geometric configuration of sensors for the three-dimensional localisation of targets using RSS measurements. The proposed evaluation function has been solved using successive optimisation algorithms to obtain optimal positions of the sensors. Xu [
31] investigated the optimal sensor placement strategy for single static target localisation using the hybrid RSS, AOA and TOA measurements on the two-dimensional plane. Lee et al. [
32] proposed a hybrid localisation algorithm to boost the accuracy of range-based localisation; they replaced the ranging part of the rule-based localisation method with a deep regression model that uses data-driven learning with the dual-band RSS. Bo et al. [
33] studied the configuration of sensors for the cooperative localisation of autonomous underwater vehicles using the RSS. Alanezi et al. [
34] proposed a method to accurately determine the position and distance of the wireless sensors linked in a local network. The method utilises the RSS at the target node to identify its location in the localised grid system.
The location of the MSs or in other words, their geometric configuration has a significant role in the quality of localisation. Xu [
31] studied this issue for the hybrid RSS, AOA, and TOA, but Eshagh [
35] applied quadratic optimisation with constraint for obtaining an optimal configuration for the MSs based on the AOA and the TDOA, and concluded the optimal design based on the TDOA is highly dependent on the geometric form of the area, and his efforts for design optimal configuration for the Landvetter airport of Sweden was not successful. However, in another study Eshagh [
36] showed that optimisation can be performed without problem if the area is close to a square and applied to the Arlanda International airport in Sweden.
This study aims to develop a method, provide all mathematical tools, and show how to apply them for the optimal design of a wireless localisation security network based on the RSSDs from a signal interference device. To do so, a configuration of a four-MS localisation network is considered and optimised using the predefined variance-covariance matrix (VCM) or the dilution of precision (DOP) for the coordinates of all nodes of a grid covering a study area. In addition, this paper shows that simply by selecting a suitable place for the MSs, the quality of localisation can significantly be improved without any extra cost. In
Section 2, the principle of the RSSD is presented, and after that, in
Section 3 the problem of localisation using the RSSD is provided.
Section 4 deals with mathematical developments of the optimisation criteria VCM and DOP for a localisation network based on the RSSD.
Section 5 discusses the optimisation model and the required limiting constraints. In
Section 6, the presented design method is applied to the Arlanda international airport in Sweden.
2. Received Signal Strength Difference
RSS depends on how long the signal has been on its way to the receiver, in other words, the distance from the transmitter. This RSS (
) has the following mathematical formula [
17]:
where
is the signal strength at the transmitter,
is a parameter, which depends on the transmitter and receiver antenna’s heights
and
and frequency of the radio signal
f.
is called pathloss exponent and
d is the distance between the transmitter and receiver. In a two-dimensional network, the following well-known relation between the Cartesian coordinates of the transmitter and receiver exists:
where
are the coordinates of the transmitter and
those of the
ith MS with the receiver. By taking logarithm from both sides of Equation (1), the known equation of the RSS is derived e.g., [
17]:
The first term on the right-hand side of Equation (3) can be a constant, because , f, and are constant and if is specified in such a way that all antennas have the same height.
There are two points regarding the localisation process using Equation (3). First, C is unknown, and second, is also unknown and dependent on the environment. In practice, our goal is to estimate the coordinates of an interfering transmitter, which sends no information about its signal power and antenna’s height. In a two-dimensional localisation problem, at least, two equations are required for estimating the x- and y-coordinate of the transmitter if the parameters C and γ are known; otherwise, four equations, or MSs, are needed to simultaneously estimate the coordinates with C and γ.
An alternative way is to perform the localisation process in two steps. The first, known as calibration, is to estimate C and γ with from a priori RSS measurements over the study area. In the second step, the estimated C and γ are used and considered as constants for the RSS mathematical model (3) and the localisation process is conducted.
By writing a differential from of Equation (3), the mathematical model of the RSSD is obtained, which is the difference between the RSSs at the MSs
i and
j (see e.g., [
20] or [
17]):
As Equation (4) shows, the RSSD is a function of coordinates of the MSs
i and
j and the transmitter in addition to
. Here,
can be also estimated with the coordinates of the transmitter if at least three RSSDs are measured from three MSs. Some closed-form solutions for localisation using the RSSD with three and four MSs are available in [
37]. Generally,
ranges between 1 and 6 [
20];
= 2 is used for a free space, and
= 4 for a flat environment. It varies between 4 and 6 for indoor environments and in some cases, such as tunnels, it is less than 2.
3. Localisation Based on Received Signal Strength Differences
By assuming that
is a known constant, six observation equations of Equation (4) type can be constructed for four MSs. These mathematical models need to be linearised and solved iteratively. Their matrix form is represented by the following Gauss-Markov model (cf. [
38]):
where
x is the vector of the coordinate updates to the approximate coordinates of the transmitter, and
the vector of random errors with
, where
stands for the statistical expectation,
the VCM of the observations, and finally
a priori variance of unit weight. The matrix
A contains partial derivatives of the observables, the Formula (4), with respect to the unknown
x- and
y-coordinates of the transmitter. The general formulae for these derivatives are:
For the four MSs M, N, O and P, and the two unknown coordinates of the transmitter,
A has six rows and two columns. According to the formulae of the partial derivatives presented in Equation (6), the structure of
A is
In Equation (5),
L is the vector of differences between actual and approximate observations computed from the coordinates of MSs and the approximate coordinates of the transmitter
, i.j = M, N, O, P are the RSSDs generated from the coordinates of the MSs and the approximate coordinates of the transmitter. , i = M, N, O and P are the computed approximate distances from the transmitter.
The least-squares solution of Equation (5) for the vector of coordinate updates is (see e.g., [
39]):
where ()
T stands for the transposition operator of matrix algebra.
is added to the approximate coordinates of the transmitter, and after that new
A and
L are computed leading to another new
is estimated. This process is iterated until the coordinate updates
do not significantly change the transmitter coordinates; in other words, the solution converges. Note that the choice of
has no effect on
. The matrix
A of the last iteration is used for the estimation of the VCM of the coordinates. By assuming
and
Q =
I, this VCM is:
can also be estimated via
where the denominator 4 is the degree of freedom of the localisation problem since there are six observables and two unknowns in the created system of equations. However, the main goal of this study is to optimally design a localisation network and not localisation of a transmitter, therefore, no observations exist in our system, and
can be simply used in our computations, because it does not affect the geometry of the localisation network.
4. Optimisation Criteria for Localisation Networks
For designing or analysing a localisation network a grid of points covering the study area is designed, where each node of the grid is a representative of the probable location of the interfering transmitter. A local coordinate system is simply defined for this grid, e.g., the node in the lowest left corner of the grid can be considered as the origin, and based on the area, the orientation of the grid can be specified. By such definitions, the coordinates of the nodes are derived according to the grid resolution. The transmitter is located amongst these nodes in the localisation process; see [
35,
36]
Generally, for computing the VCMs of these nodes, the coordinates of the MSs are needed. The initial coordinates for these stations are specified in the defined local coordinate system. Therefore, many A matrices, created between the MSs and nodes, are computed solely from their coordinates, and therefore the VCMs from A matrices; see Equation (10). The total number of VCMs is equal to the number of grid nodes. Each VCM is a tool connecting the configuration, created between the transmitter and the MSs, to the quality of the estimated coordinates for the transmitter. This means that by varying the coordinates of the MSs the quality of the transmitter coordinates changes.
In a two-dimensional network, each VCM has four elements, two variances as diagonal and two equal covariances as off-diagonal. The desired VCMs or the square root of their traces (DOPs) for the nodes, can be regarded as two criteria for the optimisation of the network. In fact, the coordinates of the MSs vary until the estimated VCMs or DOPs are fitted to the desired ones.
In the following subsections, the mathematical models of these criteria and their relationship with the coordinates of the MSs are presented.
4.1. Variance-Covariance Matrix as a Criterion
When the VCM is selected for optimisation, its direct mathematical relation with the coordinates of the MSs is needed. To derive it, let us expand the VCM by the Taylor series for a node and four MSs of M, N, O and P (see [
35,
36,
40,
41,
42,
43]):
where
is the estimated VCM of the node of the initial design,
and
are the coordinate updates to the initial locations of the MSs, and
and
are, respectively, the partial derivatives of the VCM with respect to the coordinates of the
ith MS, with the following expressions (see [
35,
36]):
where the elements of
and
, are partial derivatives of
A or the second-order partial derivatives of the observables. Note that the first-order derivative is taken with respect to the transmitter coordinates and the second-order one with respect to
xi and
yi -coordinates of the
ith MSs. The following general formulae are presented for the second-order partial derivatives of observable Equation (4) with respect to coordinates of
ith and
jth MSs:
From these general formulae, structures of the derivatives of
A with respect to coordinates of M, N, O and P, are:
4.2. DOP as a Criterion
DOP is the square root of the trace of the VCM of a node and its Taylor expansion with respect to the coordinates of the MSs is:
where trace() stands for the trace operator or the sum of the diagonal elements of a square matrix,
and
are respectively partial derivatives of DOP with respect to the coordinates of the
ith MS with the following formula:
The structures of and have been already discussed in the previous section.
5. Optimisation Model and Required Constraints
The expanded VCM and DOP by the Taylor series, Equations (12) and (22), are our mathematical models for optimal estimation of updates to the coordinates of the MSs. On the left-hand side of these equations, the desired VCMs and DOPs are, and their right-hand sides are computed from the coordinates of a node and the initial coordinates of the MSs. The estimated initial VCMs and DOPs do not have a good fit for the desired ones, therefore, by changing the coordinates of the MSs using quadratic optimisation, we try to fit them.
It is important to note that the number of unknowns, or the coordinate updates, depends on the number of the MSs. Since four MSs are used in this study, then, the number of unknowns becomes eight as the network is two-dimensional. Each VCM has four elements meaning that four equations are created for one node. Such a system is under-determined, and more equations are needed to solve the system. By adding the VCM of another node, the number of equations increases by four. Consequently, considering each extra node leads to four additional equations in the system. As mentioned before, a grid of nodes, covering the area, is designed in the earliest stage where the resolution and size of the area specify the number of nodes. In short, there are many nodes over the area and the number of equations is four times larger than the number of nodes. Therefore, there is enough redundancy for optimisation.
The situation is rather similar for optimisation based on DOP with only one difference. For each node, only one DOP is defined, and the number of equations is equal to the number of nodes. The system of equations is not as redundant as that one based on the VCM, but both have enough redundancies.
Let us present the overdetermined system in the following Gauss-Markov form:
where
stands for the MSs’ coordinate updates
is the vector of residuals, and
the vector of differences between the elements of the criterion and initial VCMs or their DOPs,
stands for the coefficient’s matrix containing the partial derivatives of the VCM, or DOP with respect to the MSs’ coordinates. In the case of using the VCM for optimisation, we have
where operator “vec” inserts the columns of the VCM below each other and converts the 2 × 2 matrices to 4 × 1 vectors,
n means the number of the nodes.
When the DOP is used for optimisation, we have:
The system of Equation (24) should be solved for the coordinate updates in a least-squares sense, but Eshagh [
35] has mentioned that applying some constraints in the quadratic optimisation model is necessary, otherwise, instability and improper configurations are obtained. Therefore, the following optimisation model for solving the system (24) is suggested:
The first constraint (31) controls the coordinate updates in such a way that the MSs move towards some specified directions. The mathematical derivations of them are available in [
36].
D and
d have the following structures:
where Az stands for the specified azimuth in which the MSs should move during the optimisation process. M’, N’, O’ and P’ are, respectively, some helpful points that M, N, O and P should move towards them.
The second constraint (32) define the search area around each MS, and the coordinate updates are computed in such a way that the MSs remain inside the search area; for details see [
35,
36]:
and , i = M, N, O and P, are, respectively, the lower and upper bounds of the search area for the x-coordinate of the ith MS, and and , are similar ones for the y-coordinate.
6. Design of Interference Localisation Security Network for the Arlanda International Airport Based on Received Signal Strength Differences
For testing our methodology, we select the Arland international airport of Sweden, in the northern part of Stockholm. The airport is almost square in shape and has three runways. The goal is to optimally design a four-MS interference localisation network on the runways. The VCM and DOP criteria are defined and applied to optimise the geometric configuration of MSs. After that the resolution and precision are presented and discussed.
We consider rectangular search areas around each MS; see
Figure 1. As the figure shows a two-dimensional coordinates system is defined with an origin outside the airport in the south-west part of the area with the geodetic coordinates of φ = 59°37′10′′ and λ = 17°53′50′′. The
y-axis is parallel to the western runway having an azimuth of 10°, and the
x-axis is perpendicular to the
y-axis with an azimuth of 100°.
The updates to the
x-coordinates of M, N and O are estimated in such a way that these points move ±250 m from their initial values according to the width of the runways.
yM,
yN and
yO are limited, respectively, between −2000 to 500 m, +1000 and −2000 m, +1000 and −2000 m.
xP is limited between −1000 and 1000 m, and
yP between −500 and 500 m.
Figure 1 shows the photo of the Arlanda international airport (taken from Google Earth), the local coordinate system, and the initial positions of M, N, O, P, shown by the small red circles and their rectangular search areas. Our goal is to keep the MSs on the runways during the optimal design procedure.
In addition, as observed the y-axis of the system is chosen parallel to the western runway for simplification. However, the choice of the coordinate system is not important as the design can be conducted based on any chosen system and later the whole network is georeferenced. Since the eastern and western runways have an azimuth of about 10° then the y-axis of the system has the same azimuth.
Table 1 shows the coordinates of the MSs before and after optimisation based on the VCM and DOP criteria for
γ = 2, 4 and 6. As mentioned before,
γ can be estimated in the calibration step of localisation. However, since there is no measurement in the design step, then we consider these three values to test the role of
γ in the design.
γ is shown in the left most column of the table, and in the column before, the criteria of VCM and DOP are specified. Three resolutions of 20, 40 and 80 m are considered to check their effects on the design.
As seen in
Table 1,
xM,
xN and
xO do not change by the resolutions, criteria and
γ because M, N and O are in the northward runways of the airport, almost parallel to the
y-axis of the defined coordinate system. The directional constraints (31) keep the movements of the MSs along these runways, therefore, it is normal to see no change in their
x-coordinates but large changes in the
y-coordinates.
From
Table 1, based on both the VCM and DOP criteria, the optimised coordinates do not show any sensitivity to the choice of
γ for the resolutions 40 and 80 m. For the resolution 20 m with the DOP criterion the optimal coordinates are not significantly different for the selected pathloss exponents. For the VCM criterion, the optimal results are the same for
γ = 4 and 6.
Generally, we can conclude that lower resolutions than 40 m and γ have no significant influence on the optimal coordinates. For the resolution 20 m and higher in the free space γ = 2, the results are different.
Table 2 illustrates the statistics of the DOP values before and after optimisation for
γ = 2, 4 and 6 and the resolutions of 10, 20 and 40 m. The chosen value for
γ has a significant influence on the DOP of the network. Convergence is important in any optimisation process and in the optimisation of a security network for localisation, having a convergent solution is not straightforward. In addition, expecting convergence for any selected values for the VCM or DOP criterion is neither logical nor realistic. For example, selecting the minimum DOP of the initial design is not realistic and the optimisation process diverges. The statistics of the initial DOPs are indicators of the right choice of criterion values. In our study, selecting the mean value of the initial DOP and optimising the network based on that in all three values of
γ leads to convergence. According to
Table 2, when
γ = 2, the mean DOP is 4.8, this means that the network is optimised by changing the positions of the MSs in such a way that the DOPs are fitted to 4.8. The mean DOPs are 2.4 and 1.6, respectively, for
γ = 4 and 6. A diagonal matrix with diagonal elements equal to the square of the mean DOP is considered the VCM criterion for all nodes.
The VCM criterion, as shown in
Section 4.1, is the expected VCM for all points, and the DOP criterion is the square root of a trace of this VCM. As
Table 2 shows, the statistics are almost the same for both chosen criteria, the reason is that in both cases the DOP of the network is presented. However, when the DOP criterion is used, there is no control over the covariances between the estimated coordinates. In addition, for each point, four equations are created in the system when the VCM criterion is applied but only one equation when the DOP criterion is utilised. Therefore, the VCM criterion adds a larger number of equations, considers covariances, and has the possibility of selecting a special structure for the VCM. For example, two different values for the variances of the
x- and
y-coordinates can be selected, as well as any correlation between coordinates, if required.
The numerical optimisation process is time-consuming and has low rates of convergence for high-resolution grids. In other words, the process needs a large number of iterations to converge. In this study, when the norm of the coordinate updates becomes less than 1 m, the iteration is stopped. Considering smaller values is also possible but with the costs of many iterations and lower convergence rate.
By assuming that the Arlanda international airport is flat, we select γ = 4.
Figure 2 is the map of the initial DOPs of the network with the locations of the MSs M, N, O and P. The maximum DOP reaches 15, and the large values are along an ellipse form going through the MSs. Most of these large values are outside the airport and in the surrounding forests, except those between stations N and O.
Figure 3 is the map of the DOPs after optimisation of the network based on the VCM criterion for the Arlanda international airport. A significant reduction is observed in the DOP values and as
Table 2 shows, the maximum reaches 4. The stations M moved slightly northward and N southward and both along the specified azimuth of 10°. Point O is located at the northern part of the runway and P with slight displacement along the runway. As observed, the high values are in the form of an ellipse passing through the MSs, this ellipse is smaller than the one created based on the initial design. Some large DOPs are seen in the southeast of the area, but they are less than 4.
Figure 4 is the map of DOPs after optimisation based on the DOP criterion. Again, large values are seen along an ellipse passing through the MSs. The DOP values are significantly reduced so that their maximum does not exceed 4. Station M is in the northern part of the runway, and N is moved to the most southern part of it. P and Q are found closer to the middle of their runways.
A comparison of
Figure 2, which is the map of the DOP of the initial design, and
Figure 3 and
Figure 4, shows that the optimisation process is successful, because a maximum DOP of 15 in the initial design reduces to 4 in both optimal designs.
Figure 2,
Figure 3 and
Figure 4 show ellipse forms for the large DOPs passing through the four MSs, but with different sizes and orientations because the optimal positions of the MSs are not the same.
Figure 2 shows large DOPs, reaching 15, from the station P towards Q and M, and smaller values, about 6, between M and N; see
Figure 2, but after optimisation using the VCM criterion, these values decrease to 4 and even smaller around Q; see
Figure 3. However,
Figure 4 shows values around 4 amongst all MSs and the map is more symmetric compared to the map presented in
Figure 3. Interestingly,
Figure 4 shows small values of DOP over the buildings in the area unlike
Figure 2 and
Figure 3, but this is accidental. If the goal of optimisation is to have small DOPs over buildings, the optimisation criterion can be designed specifically for this purpose.