GSM based CommSense system to measure and estimate environmental changes

by   Abhishek Bhatta, et al.

Facilitating the coexistence of radar systems with communication systems has been a major area of research in radar engineering. The current work presents a new way to sense the environment using the channel equalization block of existing communication systems. We have named this system CommSense. In the current paper we demonstrate the feasibility of the system using Global System for Mobile Communications (GSM) signals. The implementation has been done using open-source Software Defined Radio (SDR) environment. In the preliminary results obtained in our work we show that it is possible to distinguish environmental changes using the proposed system. The major advantage of the system is that it is inexpensive as channel estimation is an inherent block in any communication system and hence the added cost to make it work as an environment sensor is minimal. The major challenge, on which we are continuing our work, is how to characterize the features in the environmental changes. This is an acute challenge given the fact that the bandwidth available is narrow and the system is inherently a forward looking radar. However the initial results, as shown in this paper, are encouraging and we intend to use an application specific instrumentation (ASIN) scheme to distinguish the environmental changes.



There are no comments yet.


page 9

page 12

page 13

page 14

page 15

page 16

page 22


Frequency-Hopping MIMO Radar-Based Communications: An Overview

Enabled by the advancement in radio frequency technologies, the converge...

Channel Estimation for Visible Light Communications Using Neural Networks

Visible light communications (VLC) is an emerging field in technology an...

Vehicle Detection and Classification using LTE-CommSense

We demonstrated a vehicle detection and classification method based on L...

MAJoRCom: A Dual-Function Radar Communication System Using Index Modulation

Dual-function radar communication (DFRC) systems implement both sensing ...

Interfering Channel Estimation in Radar-Cellular Coexistence: How Much Information Do We Need?

In this paper, we focus on the coexistence between a MIMO radar and cell...

AI Based Landscape Sensing Using Radio Signals

In many sensing applications, typically radio signals are emitted by a r...

Some aspects of physical prototyping in Pervasive Computing

This document summarises the results of several research campaigns over ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

Classical radar systems have been designed primarily for military operations. Of late many interesting ways of using radio-frequency spectrum for radar purpose have been under research. One such concept is commensal111The word commensal has been borrowed from biology in which this represents co-existence of two species out of which one is benefited and the other remains unaffected. or passive radar, which uses signals of opportunity to detect targets without affecting the functionality of the parent system [13, 14, 12, 28, 29, 6, 26].

One implementation of commensal system is the recently active area of Passive Bistatic Radar (PBR) [9, 16, 23, 7, 11].

In our work we shall describe a kind of commensal radar which we call communication sensor (CommSense) [18, 17] system. It uses the channel estimation processing in communication systems to estimate changes in the environment. This system can, potentially, be used to monitor land terrain, sea state or even natural disasters. The novelties of the proposed system are as follows.
First of all, this system is built upon existing communication systems. The cost of implementing such a system is therefore considerably low. Secondly, unlike PBR systems it does not process the information using correlation, rather it uses a known training sequence to estimate changes in the environment. It therefore eradicates the need for a reference antenna pointing towards the transmitter. We have implemented the prototype system using open source GNURadio software and Nuand BladeRF SDR hardware[8, 20].

It must, however, be noted that there are two major challenges in implementing this system. First of all, the conventional concepts of resolution will not be valid here mainly because of the fact that the bandwidth available to us is very narrow. This makes the measurement system a severely ill-posed inverse problem. Secondly, this is a forward looking and non-coherent radar system which theoretically limits the amount of information it can capture. We hypothesize that these two problems can be solved using the application specific instrumentation framework [19, 24]. This will be our future work, whereas in this paper we demonstrate the feasibility of the concept of sensing the environment from communication signals.

The rest of the paper is organized as follows. Section 2 gives a basic understanding of GSM protocol and an introduction to the design of this system. Section 3 provides details about the real time implementation of the system. Section 4

shows how the data analysis is done, specifically focussing on the Probability Distribution Function (PDF) analysis and clustering. Finally Section

5 concludes the work and shows the scope for future developments.

2 GSM Channel Estimation and CommSense System Design

Any signal transmitted over the channel gets affected by the channel itself. This effect is minimized in wireless communication systems using post-reception processing blocks called channel estimation and equalization [22, 21]. The hypothesis of our work is that the channel information from the received signal can be harnessed and used to statistically monitor the changes in the environment. For example we will be able to differentiate between different weather conditions or produce an alternate model (which, potentially, can give more information) for sea waves and mountainous terrains. We call this system CommSense system. The reader is suggested to refer [18] for system level information of a generic CommSense system.

The communication technology infrastructure used for the current implementation of the CommSense is GSM. This is because it is currently the most widely used wireless communication technology in South Africa. We are also working towards investigating its feasibility using LTE infrastructure [25]. GSM provides coverage of almost the entire nation. It is a wireless communication protocol developed by European Telecommunications Standard Institute (ETSI) to provide a standard in the wireless telecommunication industry. CommSense system relies on the broadcast signals transmitted by the base stations for phones to detect and connect to the particular base station, thus not hampering the communication system in any way. GSM transmits a known training sequence in every frame so that the receiver is able to detect and cancel the channel interference. Using the information from this, the channel states can be estimated.

Figure 1 shows the concept diagram of the application. Signal from the base station reaches the mobile station through multiple different paths and is affected by the physical properties of each path. Out of all these returns, the signal with highest Signal to Noise Ratio (SNR) is extracted by the method of preliminary channel equalization and used for the purpose of communication. In our work we shall be using this preliminary channel equalization block. It, however, can be noted that there are further blocks in the processing chain in the receiver which use a range of algorithms to take care of channel noise. For the purpose of our work we shall use channel equalization to mean this preliminary channel equalization processing. In order to implement channel equalization the communication system transmits a known sequence of bits with every frame. These known bits, also known as training sequence, are extracted from the channel convoluted received signal and used to extract The difference between the known and expected signal. This difference between the known and the received signal gives the state of the channel through which the signal has travelled. In the proposed system this estimated information is used to find the changes in the channel state in different environmental conditions.

Figure 1: System Overview.

2.1 System Design

CommSense system acts as a commensal radar, eliminating the need of a new transmitter. As the processing power of smart-phones has significantly improved, it is planned to eventually implement CommSense as a mobile application. For the purposes of proving the hypothesis it has been implemented on the SDR platform with open source software. The hardware used for the implementation is BladeRF manufactured by Nuand [20]. It is a Universal Serial Bus (USB) powered SDR board with an on-board transmitter and receiver port. It can operate over a frequency range of 300 MHz to 3.8 GHz which gives the opportunity to implement the proposed system using GSM technology. The software used for implementing the real time channel estimation system is GNU Radio [8]. It is an open source software tool-kit, distributed under GNU general public license, that provides the opportunity to implement signal processing blocks using open source RF hardware for research purposes.

GSM is a wireless communication protocol that acts as a host to this system [1, 2, 10]. It uses Frequency Division Multiple Access (FDMA) and Time Division Multiple Access (TDMA) in order to accommodate multiple users. The various frequency bands of operation for this technology are 800 MHz, 900 MHz, 1800 MHz, 1900 MHz. GSM works in frequency duplex mode thus having a different frequency for uplink and downlink. GSM900, used here for data collection, uses 880-915 MHz for uplink and 925-960 MHz for downlink communications. These bands are further divided into 200 kHz bands and are separated by a number called Absolute Radio Frequency Channel Number (ARFCN). Each wireless channel is divided into TDMA frames with 8 time slots of 577 s each. Each timeslot is time shared between the mobile and the base station. The logical channels are piggybacked on the physical channel, described in [1].

Each frame in the GSM system transmits a known training sequence for the receiver to detect and reduce the possibility of transmission error. This functionality provided by GSM is used in this project to estimate the channel effects on the signal and analyse it to find the differences in the various types of terrain and environmental conditions. A frame in GSM system is shown in Figure 3. This structure shows the normal burst that has been used for the extraction of the parameters. This burst has 114 bits of information transmitted with tail bits, flags and a guard period to identify the start and stop of the burst and different datasets. The training sequence is 26 bits long out of which the central 16 bits are used for channel equalization. In the current implementation these equalized channel values are used to determine the surrounding environment. A guard period is added at the end of each frame to provide a window of error against distortions that occur due to the rise and fall time of the signal.

Figure 2: Channel Estimation Block Diagram.

2.2 Channel Equalization Filter

Figure 2 shows the block diagram of the implementation. The received GSM signal can be mathematically represented as equation (1). Here, the received signal is the convolution of the transmitted signal and Channel Impulse Response222CIR is a black-box linear model of the channel effects. (CIR) in the presence of Additive white Gaussian noise (AWGN) . The transmitter sends a known training sequence in each frame as shown in Figure 3, which is divided into reference length of and guard period of bits [22]. This equation can be represented in the matrix form and is shown in equation (2).

y (2)

In equation (2) M represents a Toeplitz like matrix structure of the given training sequence as shown in equation (4). This matrix is made of the m array that is shown in equation (3), which is actually the oversampled information from each received frame after filtering out the Gaussian noise.

Figure 3: GSM Normal Channel Frame Structure.

In equation (4), the matrix M is a matrix where is the reference length of training sequence and is the length of the CIR. equation (5) shows the CIR array, here each value of the array represents the channel parameters of one reflected signal. Thus the length of the array determines how many multi-path signal are being analysed. The Least Squares algorithm finds the CIR by minimising the squared error quantity which in the presence of white Gaussian noise gives equation (6).


Here, denotes Hermitian transpose matrix and denotes matrix inverse. In equation (6) we can observe that the received signal matrix y is multiplied to a matrix also known as the pseudo inverse matrix of M.

Figure 4 shows the plot for simulated CIR in the presence of AWGN. This result is for a single GSM frame being transmitted over a known channel of length 5 taps. The estimated value is very similar to the original CIR value which shows that the algorithm is working properly. Then this is implemented in real-time as shown in Figure 5.

Figure 4: Simulated Least Square Channel Estimation in presence of AWGN.

3 Real Time Implementation

The system is implemented in GNU Radio, which does the channel estimation and plots the CIR in real time as shown in Figure 6. The hardware used is BladeRF, which works at a frequency band of 300 MHz to 3.8 GHz with a maximum physical bandwidth of 28 MHz. The antenna used for this specific implementation is an off the shelf quad band GSM antenna. The accuracy of the system is proved by decoding the received bit-stream and analysing it in Wireshark (a packet analyser tool), which gives the basic information as the location of the base station, carrier service provider and so on.

3.1 Extracting the Channel Values

The steps taken to implement the channel equalizer in real-time are mentioned here. Correlation between the received signal and the known training sequence is used to extract the differences between the two signals and this information is saved in a file for further analysis.

The SDR hardware, bladeRF, receives the GSM signal downconverts it according to its Phase Locked Loop (PLL) clock frequency and gives out the In-phase () and the Quadrature () components of the signal. The receiver first needs to get synchronized with the base station, the Synchronization Channel (SCH) burst is used for this purpose. The SCH burst consists of a 64 bit long training sequence and transmits the same sequence in every timeslot allocated to it, thus it is easily distinguishable and the receiver can get synchronized. Once the receiver synchronizes with the GSM base station, the frequency offset is calculated using the Frequency Correction channel (FCCH).

Figure 5: Estimated Normal channel absolute values in terms of the and component, in this figure is represented by .

The received burst type is identified in two steps first with the help of the carrier index of the burst and then with the burst number. Carrier index zero (C0) is a special case where only Broadcast Channel (BCCH) data is transmitted. BCCH carries a repeating pattern of system information such as identity, configuration and available features for the base station. In case of other carrier indices the burst type is determined by the burst number. The most commonly transmitted burst is the normal burst thus in this implementation we are focussing on extracting the channel information from the normal burst.

In order to receive a burst completely the receiver needs to wait for two consecutive guard periods. Although there will be an overlap between the guard periods of two consecutive bursts, this is necessary to ensure full burst is received. Once a normal burst is identified the start and the stop position of the burst is determined eliminating the guard period as it consists of 8.25 bits on each side and is distinguishable. The central 26 bits of the remaining bits are extracted as it constitutes of the received training sequence. The training sequence extracted here is a complex data so the already available sequence also needs to be converted into a complex number. The two complex training sequences are correlated with each other and the correlation information is saved in a buffer. A pointer is placed at the beginning of the buffer and based on predefined CIR length the pointer is moved from the beginning to the length of the CIR. The values corresponding from the beginning to the CIR length is saved in an array and then moved to a file where it is saved for further analysis. In case of communication systems the maximum absolute value of the CIR is chosen discarding the rest of the values and that information is used to equalize the channel and get better communication. For the purpose of this implementation the entire channel information is necessary thus all of it is saved.

In order to check if the receiver is working as intended the entire receive structure is made referring to the gr-gsm libraries available for GNU Radio and the bursts are printed out using a message printer. Figure 6 shows the plot of the entire implementation including the GSM frequency spectrum, received bursts and the real time estimated channel values.

Figure 6: GNU Radio implementation of Channel Equaliaztion with BladeRF.

3.2 Real Time Channel Estimation and Data Capturing

Figure 2 shows the entire flow of the implementation, where the signal source is the transmitting base station and the processing is done on over the air signals. The receiver filter block removes most of the noise from the signal and shows the GSM frequency spectrum as shown in Figure 6. The data is used to estimate the channel and store the values in a file for further analysis. The darker blocks in Figure 2 are the normal GSM processing blocks which are not affected by this system as it runs parallel to the communication system.

One of the major challenges faced during the implementation is finding the start-point of the burst. This is done with the help of the guard period. The receiver waits for two consecutive guard periods, once received it considers the start point of the burst. In the case of an error in receiving the guard periods, the system drops the current set of data and waits for the next set and then repeats until it gets two consecutive guard periods, thereby eliminating the chances of any ambiguous data.

4 Analysis of Captured Data

To prove the hypothesis the real-time channel information is captured in different locations and classified as Statistical Characterization, Distinguishing Events and Special case. Table


lists details of all the datasets collected and analysed in this work. The estimated channel information is analysed to show the changes in the channel state in different environmental condition. Three different types of analysis are presented in this section, two in statistical domain and one in data/time domain. In the first approach we estimate the PDF of the channel values captured at different scenarios. This shows that the estimated PDF’s vary as the environmental conditions change which means that a “hypothesis test” kind algorithm can be used to distinguish different environment types from the estimated channel values. In rest of the paper we shall use simply PDF to mean estimated PDF. In the second approach we perform Chi-square test to assess the goodness of fit between the observed data and the theoretically expected dataset. In the third approach we check for the clustering of the estimated channel values for different capture scenarios. For easier visualisation we used principal component analysis (PCA) to reduce the dimension of the data. This shows that the data are highly clustered in PC domain.

Captured Data Type Location Comments Statistical Characterization Sea Beach To check the effects of the sea waves on the signal, multiple locations may also be tested Hill next to a road Check the effects on the signal in the presence of a hill and see repeatability (safe location still needs to be decided) Highland Captures taken at the road next to a steep terrain Heavy Rain, slightly humid climate, hot day Captures taken at the same location at different instances of the day showing the differences due to different weather conditions Parking space Captures are taken on different times of the day. One set showing when the place is full and the other showing empty parking lot Distinguishing events Train Station At a train station with and without a train in proximity J stairs At the stairs to the entrance of an approximately 10 meter tall building. Bus With and without the presence of a bus in the proximity, two different sets of captures were taken at different locations Special case Car Captures taken with and without the presence of a car in proximity of the receiver Corner Reflector Captures are taken when a corner reflector is placed at a distance of 2 meter from the antenna in various different configurations ’V’ is vertical, ’H’ is horizontal, ’VH’ refers to a dihedral corner reflector
Table 1: The different sets of captured data.

4.1 PDF Analysis

In this set of analysis we plot the empirical PDF of the data collected from different scenarios as described in Table 1

. In addition to the empirical PDF we also plot the PDF of the data for four different theoretical PDF models, i.e. Rayleigh, Normal, Log-normal and Gamma distributions. A brief description of these distributions and maximum likelihood estimation (MLE) expression of their parameters are given in Appendix I.

Figure LABEL:fig:hist_rep1 and LABEL:fig:hist_rep show the distribution plot of two different datasets in similar conditions. Both the captures are taken near the entrance of a building about 10 meters tall. The difference between the two captures is the wind speed. During the capture of Figure LABEL:fig:hist_rep1 the wind speed was slower compared to the capture of Figure LABEL:fig:hist_rep. The location picture of the J stairs capture set is provided in Appendix II.

The plots display index 12 of two different datasets whose indices go up to a total of 40 which is the number of multipath signals extracted from the received signal. Here we observe the similarities in the distribution pattern of the captures. Although the empirical PDF are not completely matched, the general pattern is similar.

This proves that the data captured from a single location will generally follow a similar pattern in a particular environmental condition (wind, humidity, temperature, etc.), thus aiding in providing enough information to map the surface and the environment at the time of the capture.

Figure 7: Empirical PDF comparison of index taken from a 40 index dataset near a building with different wind speeds. In the legend, rayleigh, norm, lognorm and gamma

represents theoretical PDF assuming Rayleigh distribution, Gaussian distribution, Lognormal distribution and Gamma distribution respectively.

Next we analyse the empirical PDF and theoretical PDF of estimated channel values from different scenarios. Some of the example empirical PDFs are shown in Figures 8,9,10,11. Each curve shows the empirical PDF of a particular case. Keeping the -axis constant shows the differences in the given datasets clearly. The datasets presented here include the case of car in proximity, near a big building, bus in proximity, full parking space, train in proximity and the special case with and without the corner reflector.

Dataset Shape Mean Variance Skew Kurtosis Log-normal Gamma jameson_stairs jameson_stairs2 rts_t roof_bcr_loc2 roof_vh_metal_cr north_stop_jammie_loc1 car rhodes_mem parking_full
Table 2: Moments of captured datasets.
Figure 8: Empirical PDF of the captured data in the presence of a bus and a car showing the differenecs in the distributions due to the presence of particular targets.
Figure 9: Empirical PDF of captures at a full parking lot and near a train in a local train station showing the differenecs in the distributions.
Figure 10: Empirical PDF of captures infront of a meter tall building and near a hill. Showing the difference in the distributions between natural objects and man made objects.
Figure 11: Empirical PDF of captures with and without the presence of a dihedral corner reflector.

In Figure 8 and 9 we observe the change in variance with change in the target conditions. In the case of a bus, the distribution has a lower variance which increases as the target changes to a car and even more in case of the train. Looking at Figure LABEL:fig:hist_car and LABEL:fig:hist_parking a similar pattern can be observed with minimum changes in the distribution. This is because both the captures represent cars in a flat space, the only difference is in Figure LABEL:fig:hist_car there is only one car in proximity and in Figure LABEL:fig:hist_parking there are many cars which creates the observable differences in the plots. Figure 10 gives the distribution comparison in case of two different situations, near a hill and near a tall building. The variance in case of the hill is a lot less compared to the variance in case of the building thus differentiating between them. Figure 11 is a special case where the captures are taken with a corner reflector placed 2 meter from the receiver antenna and the differences in the plots are clearly visible. Table 2 provides the various moments extracted directly from the captured data. The shape parameters are extracted by fitting the empirical data to the expected distribution, thus there is a difference between the log-normal and gamma distribution shape parameters. These correlates with the information visible in the plots and gives us a better understanding of the changes in parameters due to the change in environmental conditions.

From the plots and table provided here it is observable that there are differences in the distributions of the channel information due to the changes in the physical environment. This enables us with the opportunity to utilise these differences and model a system that can differentiate between them. Although the differences are quiet clear from the empirical PDF’s, it is always beneficial to further analyse the received data and gain better understanding.

4.2 Chi-square Test

The second set of analysis performed on the captured dataset is the Chi-square test [30, 15]. This is a test which takes two inputs, the observed data and the expected data and gives out a score based of the two following the equation (7

). The Degree of Freedom (DF) in this case is 100 since the length of the array is 101 and DF is defined as one less than the length of the array.


Here is the observed value and is the expected value. In order to perform this analysis PDF of the received data (observed data) and the distribution fitted data (expected data) are passed through equation (7) and the output of the test is recorded in Table 3. The p-value is defined as the probability of obtaining a specific value equal to or higher than what actually is observed, maintaining the assumption that the model is correct. The points where the p-value is null can be rejected because this means the probability of getting a value similar or greater than that is null. All other cases in this test can be accepted and Table 3

shows that Log-normal distribution has the most consistent result. In the next subsection we reduce the dimensionality of the data and gain better visualization by rotating the axes.

Chi-Square Value Capture Set Names Rayleigh Gaussian Log-normal Gamma Chi-Square p-value Chi-Square p-value Chi-Square p-value Chi-Square p-value jameson_stairs jameson_stairs2 rts_t roof_bcr_loc2 roof_vh_metal_cr north_stop_jammie_loc1 car rhodes_mem parking_full
Table 3: Chi-Square test values for the received data matched with the fitted data for each distribution.

4.3 Principle Component Analysis (PCA)

In this subsection we investigate the clustering of the estimated channel values in different scenarios. In order to better visualize different sets of the captured data we used PCA. With this analysis we can reduce the dimensionality of the data and view the datasets from an angle that provides maximum information. Here we have observed different cluster formation of the datasets due to the change in locations and environmental conditions. There are many different ways to reduce the dimensionality of the data and calculate the principle components, for this implementation we have used Singular Value Decomposition (SVD)

[3, 4, 27]. The major goals of PCA are to reduce the dimensionality of the data, extract the important information and analyse the structure of the data. The calculations for PCA are shown in [5].

In order to derive the principle components from the data we first need to generate the SVD equivalent of the dataset A. The input matrix A has sets of data explained by variables, represented by . A has a rank with , then the SVD of A will be given by equation (8).

A (8)
F (9)

Here U is

matrix of left singular vectors,

V is matrix of right singular vectors and is a diagonal matrix of singular values The components for PCA are obtained from the data A using equation (8). The principal component matrix F of dimension is given by equation (9). To get the coefficients of the linear combinations which are used to compute the factor scores the matrix V is used. The matrix can also be interpreted as the projection matrix because gives the projection values of the observations on the principle components as shown in equation (10).


Geometrically the components can also be represented by rotating the original axes and the matrix A can be interpreted as a product of the factor scores given by equation (11).[5] Here I

is the identity matrix.

Figure 12: PCA of different climatic conditions at a single location. The conditions include high humidity but no rain, heavy rain and Hot day with very low humidity.
Figure 13: PCA of different datasets showing the different clusters.

All the plots after PCA are done by stacking the normalized datasets together and comparing the

principle components. The dataset dimensions have been reduced from 40 components to 2 major components and the relationship between the components are plotted to show the best view. The clusters are clearly visible for each data. It can be noted here that taking more components will get us more information and hence will make scene classification even easier.

Figure LABEL:fig:PCALoc1 shows the the clusters of the data captured at sea beach, near the shore, beach rock, at the shore but at a height of around 10 meters from the sea level. J stairs, stairs at the entrance of a building about 10 meters tall. Train station, gives the captures taken at a local train station in two different situations first, when there is a train about 6 meters from the receiver and secondly in the absence of the train. Highland, refers to a location where the captures are taken on the road next to a small hilly terrain. Rhodes memorial, is situated mid way up on to a mountain in Cape Town.

Figure LABEL:fig:PCALoc2 shows a particular set of captures taken with a corner reflector. The different clusters in this figure shows the different set of captures taken at various different configurations. CR is an acronym for corner reflector the different allocations are in this plot. Figure LABEL:fig:PCALoc3 shows the clusters for specific set of captures taken at a bus stop on the campus of University of Cape Town (UCT) at two different locations. The university shuttle service is known as ”jammie”, thus the name given. The different scenarios for this capture are, in the presence of a jammie about 4 meters away from the receiver antenna in location 1, 2 and in the same locations when there is no jammie.

Figure LABEL:fig:PCALoc4 shows the captures taken at a flat parking space when there are no cars and in the presence of cars. This set also shows the plots for a single car dataset, where the captures are taken with and without a car in the same location on two different days (one at the morning 8 am and the other in the evening 6 pm) with different climatic conditions. Figure 12 shows the difference in the clusters at the same location when there is heavy rain, medium rain and on a hot day.

5 Conclusion

In this work we present a novel way of sensing the environment using communication signals, which we call CommSense. The hypothesis under test was that we will be able to sense the environment by analysing the channel equalization values of any communication system. To prove our hypothesis we chose GSM communication. We have implemented a very basic correlation based techniques to estimate the channel values from the received signal using the known training sequence. In the limited drive to prove the hypothesis we collected data from different scenarios. On analysis it was shown that the probability distribution of the collected data differs as the environment changes. Hence a classic Hypothesis-test type algorithm can distinguish the scenes. Secondly, the data for different scenarios are highly clustered. This shows that a naive nearest neighbour algorithm can distinguish the environmental conditions easily. With these we prove our hypothesis. In the future work we intend to gather more data and characterize the environment for robust classification. We also intend to work on increasing the number of receiver nodes as another way towards getting robust environment sensing using CommSense.

Appendix I

In this appendix we shall briefly describe four of the most used PDF models.

5.1 Normal Distribution (Gaussian Distribution)

The most general form of distribution used in statistical analysis is the gaussian distribution. The probability density function follows the equation (



= location parameter (mean)

= scale parameter (standard deviation)

Standard Normal Distribution where is given as equation (13)


The maximum likelihood estimation for the gaussian distribution is derived by minimizing the log-likelihood function of the equation (12) given by equation (14),(15)


5.2 Rayleigh Distribution

One of the most common distributions used in wireless communication systems to model the channel is Rayleigh distribution. The probability density function is given by the equation (16).


= scale parameter (mode)

The maximum likelihood estimation for the parameters of rayleigh distribution is given by equation (17)


5.3 Gamma Distribution

Gamma distribution is also a very versatile distribution which can be manipulated in many ways to give other distributions such as k-distribution which is in the scope of the future work. The general formula for the probability density function of gamma distribution is given by equation (18)


= shape parameter
= location parameter
= scale parameter
is the gamma function given by equation (19)


The case where is called the standard gamma distribution and is given by equation (20)


The maximum likelihood estimates for the two parameter gamma distribution are calculated by solving the following equations (21) and (22) simultaneously


with denoting the digamma function given by equation (23), which is the mathematical derivative of the gamma function. These functions are solved by using python’s scikit learn package.

Distribution Shape Mean Variance Skew Kurtosis Gaussian N/A Rayleigh N/A Log-normal Gamma
Table 4: Moments of different distributions shown in Figure 14.
Figure 14: All distributions explained above (y-axis is in log scale).

5.4 Lognormal Distribution


is a random variable distributed lognormally then

is normally distributed where in natural log. It is very helpful in distinguishing the dataset in this work because it gives the best fit so far. The probability density function is given by equation (24)


= shape parameter (standard deviation)
= location parameter
= scale parameter (median)

The standard lognormal distribution is given when and is denoted as equation (25)


The lognormal distribution is commonly characterised with its mean given as , using this parameter the density function changes to equation (26)


Maximum likelihood estimation of this distribution is given by equation (27),(28),(29)


Table 4 and Figure 14 shows the moments of the above mentioned distribution.

Appendix II

The capture named “J Stairs” or “jameson_stairs” or “jameson_stairs2” is taken at a location in UCT as shown in Figure 15.

Figure 15: Location for capture named “J Stairs” or “jameson_stairs” or “jameson_stairs2”


  • [1] Digital cellular telecommunications system (phase 2+); multiplexing and multiple access on the radio path. ETSI TS 145 002 V12.3.0, (2015-01). (3GPP TS 45.002 version 12.3.0 Release 12).
  • [2] Digital cellular telecommunications system (phase 2+); radio transmission and reception. ETSI TS 145 005 V12.4.0, (2015-01). (3GPP TS 45.005 version 12.4.0 Release 12).
  • [3] Hervé Abdi.

    The eigen-decomposition: Eigenvalues and eigenvectors.

    Encyclopedia of measurement and statistics, pages 304–308, 2007.
  • [4] Hervé Abdi. Singular value decomposition (svd) and generalized singular value decomposition. Encyclopedia of measurement and statistics. Thousand Oaks (CA): Sage, pages 907–912, 2007.
  • [5] Hervé Abdi and Lynne J Williams. Principal component analysis. Wiley Interdisciplinary Reviews: Computational Statistics, 2(4):433–459, 2010.
  • [6] Christian R Berger, Bruno Demissie, Jörg Heckenbach, Peter Willett, and Shengli Zhou. Signal processing for passive radar using ofdm waveforms. Selected Topics in Signal Processing, IEEE Journal of, 4(1):226–238, 2010.
  • [7] C Coleman and Heath Yardley. Passive bistatic radar based on target illuminations by digital audio broadcasting. Radar, Sonar & Navigation, IET, 2(5):366–375, 2008.
  • [8] GNU Radio Website. Gnu radio website., accessed October 2015.
  • [9] Hugh Griffiths. Passive bistatic radar. RTO Educational Notes-Lecture Series RTO-EN-SET-133 Waveform Diversity for Advanced Radar Systems, Brno, 2009.
  • [10] Friedhelm Hillebrand. GSM and UMTS: the creation of global mobile communication. John Wiley & Sons, Inc., 2002.
  • [11] J Homer, K Kubik, B Mojarrabi, ID Longstaff, E Donskoi, and M Cherniakov. Passive bistatic radar sensing with leos based transmitters. In Geoscience and Remote Sensing Symposium, 2002. IGARSS’02. 2002 IEEE International, volume 1, pages 438–440. IEEE, 2002.
  • [12] Michael Inggs, Cunsheng Tong, Roaldje Nadjiasngar, Graham Lange, Anadi Mishra, and Francois Maasdorp. Planning and design phases of a commensal radar system in the fm broadcast band. Aerospace and Electronic Systems Magazine, IEEE, 29(7):50–63, 2014.
  • [13] Michael Inggs, Andrew van der Byl, and Cunsheng Tong. Commensal radar: Range-doppler processing using a recursive dft. In Radar (Radar), 2013 International Conference on, pages 292–297. IEEE, 2013.
  • [14] Michael R Inggs, CA Tong, Akhilesh Kumar Mishra, and FDV Maasdorp. Modelling and simulation in commensal radar system design. In Radar Systems (Radar 2012), IET International Conference on, pages 1–5. IET, 2012.
  • [15] MG Kendall and A Stuart. The advanced theory of statistics vol. 2, inference and relationship. charles griffin and co., ltd, 1961.
  • [16] P Krysik, K Kulpa, M Baczyk, Ł Maślikowski, and P Samczynski. Ground moving vehicles velocity monitoring using a gsm based passive bistatic radar. In Radar (Radar), 2011 IEEE CIE International Conference on, volume 1, pages 781–784. IEEE, 2011.
  • [17] Amit Kumar Mishra. Monitoring changes in an environment by means of communication devices, April 20. PCT Application No: PCT/IB2016/052235.
  • [18] Amit Kumar Mishra. Commsense: Radar system using existing communication infrastructure to sense the environment. In Radioelektronika (RADIOELEKTRONIKA), 2015 25th International Conference, pages 276–279. IEEE, 2015.
  • [19] Amit Kumar Mishra and Santu Sardar. Application specific instrumentation and its feasibility for uwb sensor based breast cancer diagnosis. In Power, Control and Embedded Systems (ICPCES), 2010 International Conference on, pages 1–4. IEEE, 2010.
  • [20] Nuand Website. Nuand website., accessed October 2015.
  • [21] Liang Pu, Jian Liu, Yuan Fang, Wei Li, and Zhisen Wang. Channel estimation in mobile wireless communication. In Communications and Mobile Computing (CMC), 2010 International Conference on, volume 2, pages 77–80. IEEE, 2010.
  • [22] Markku Pukkila. Channel estimation modeling. Nokia Research Center, 2000.
  • [23] P Samczynski, K Kulpa, M Malanowski, P Krysik, and L Maślikowski. A concept of gsm-based passive radar for vehicle traffic monitoring. In Microwaves, Radar and Remote Sensing Symposium (MRRS), 2011, pages 271–274. IEEE, 2011.
  • [24] Santu Sardar and Akhilesh Kumar Mishra. ASIN-based UWB radar for sludge monitoring. Access, IEEE, 2:290–300, 2014.
  • [25] Santu Sardar and Amit Kumar Mishra. LTE-CommSense: LTE communication infrastructure based sensing for environment monitoring. In 6th International Conference on Computers and Devices for Communication (CODEC-15). IEEE, 2015.
  • [26] Hongbo Sun, Danny KP Tan, and Yilong Lu. Aircraft target measurements using a gsm-based passive radar. In Radar Conference, 2008. RADAR’08. IEEE, pages 1–6. IEEE, 2008.
  • [27] Yoshio Takane. Relationships among various kinds of eigenvalue and singular value decompositions. In New developments in psychometrics, pages 45–56. Springer, 2003.
  • [28] Danny KP Tan, Hongbo Sun, Yilong Lu, M Lesturgie, and HL Chan. Passive radar using global system for mobile communication signal: theory, implementation and measurements. In Radar, Sonar and Navigation, IEE Proceedings-, volume 152, pages 116–123. IET, 2005.
  • [29] Danny KP Tan, Hongbo Sun, Yilong Lu, and Weixian Liu. Feasibility analysis of gsm signal for passive radar. In Radar Conference, 2003. Proceedings of the 2003 IEEE, pages 425–430. IEEE, 2003.
  • [30] Ali Vardasbi, Mahmoud Salmasizadeh, and Javad Mohajeri. Multiple-chi-square tests and their application on distinguishing attacks. In Information Security and Cryptology (ISCISC), 2011 8th International ISC Conference on, pages 55–60. IEEE, 2011.