Spectrum Assignment Algorithm for Cognitive Machine-to-Machine Networks

A novel aggregation-based spectrum assignment algorithm for Cognitive Machine-To-Machine (CM2M) networks is proposed. The introduced algorithm takes practical constraints including interference to the Licensed Users (LUs), co-channel interference (CCI) among CM2M devices, and Maximum Aggregation Span (MAS) into consideration. Simulation results show clearly that the proposed algorithm outperforms State-Of-The-Art (SOTA) algorithms in terms of spectrum utilisation and network capacity. Furthermore, the convergence analysis of the proposed algorithm verifies its high convergence rate.


Introduction
Today, there are around 4 billion M2M devices in the world, while in 2022, the number is expected to reach 50 billion [1].According to Cisco systems, currently a single M2M device can generate as much traffic as 3 basic-feature phones; in addition, emerging applications and services of M2M networks are expected to increase average traffic per device from 70 MB per month in 2014 to 366 MB per month in 2018 [2].Because of the growth rate of the number of devices and high demand of data traffic, future M2M networks will face many challenges, especially with the so-called spectrum scarcity problem.
Cognitive Radio (CR) is introduced as a promising solution to tackle spectrum scarcity problem in M2M networks.CR has become one of the most intensively studied paradigms in wireless communications.In CR, unlicensed users exploit CR technology to opportunistically access licensed spectrum as long as interference to LUs is kept at an acceptable level [3].A number of M2M applications (such as smart grid, healthcare, and car parking) can benefit from the combination of CR and M2M communications [1].CM2M networks can improve spectrum utilisation and energy efficiency in M2M networks [4].The CM2M device can interact with the radio environment by either performing spectrum sensing or accessing spectrum databases or both of them to detect spectrum opportunities [4].After sensing, CM2M device utilises the discovered unused spectrum according to the device requirements.
Furthermore, TV bands (VHF/UHF), which have highly favourable propagation characteristics, are traditionally reserved to broadcasters.But after the transition from the analogue broadcast television system to the digital one, a huge number of TV channels (also known as TV White Spaces (TVWS)) are freed up and unused.In September 2010, the Federal Communications Commission (FCC) released significant rule to enable unlicensed broadband wireless devices to use TVWS.Unfortunately, due to spectrum fragmentation and as a result of an inefficient command and control spectrum management approach, a continuous wide segment of TVWS is rare in many countries including the United Kingdom.
As CM2M network can sense and be aware of its radio environment, the aggregation of narrow spectrum opportunities becomes possible.Spectrum aggregation provides wider bandwidth and higher throughput for the CM2M devices.CM2M devices can access discontinuous portions of the TVWS simultaneously by means of Discontinuous Orthogonal Frequency Division Multiplexing (DOFDM) [5,6].
DOFDM is a multicarrier modulation technique and is a variant of OFDM used to aggregate discontinuous segments of spectrum.The main difference between OFDM and DOFDM is ON/OFF subcarrier information block [7].A multiple segments of spectrum can be occupied by other CM2M devices or LUs.As a result, these subcarriers are offlimits to the CM2M devices [6].Thus, to avoid interfering with these other transmissions, the subcarrier within their vicinity is turned off and unusable for CM2M devices, as shown in Figure 1.Moreover, available (usable) subcarriers are located in the unoccupied segments of spectrum, which are determined by spectrum broker.
Spectrum aggregation is one of the most important LTEadvanced technologies from physical layer perspective and standardised in LTE Release 10 [8].However, in spite of standardisation of spectrum aggregation, little effort has been made to optimise spectrum aggregation by exploiting CR technology in M2M networks.There is limited literature available on spectrum assignment among CM2M devices having spectrum aggregation capabilities.
In [9], an Aggregation-Aware Spectrum Assignment Algorithm (AASAA) is proposed to aggregate discrete spectrum fragments in a greedy manner.The algorithm in [9] utilises the first available aggregation range from the low frequency side and assumes that all users have the same bandwidth requirement.
Huang et al. [10] proposed a prediction based spectrum aggregation scheme to increase the capacity and decrease the reallocation overhead.The proposed scheme is referred to as Maximum Satisfaction Algorithm (MSA) for spectrum assignment.The main idea is to assign spectrum for the user with larger bandwidth requirement first, leaving better spectrum bands for remaining users, while taking into consideration different bandwidth requirements of users and channel state statistics.However, MSA does not enhance spectrum utilisation by reusing spectrum within unlicensed network; that is, CCI is neglected in MSA.
Recently, genetic algorithm (GA) is used for spectrum allocation [11].Ye et al. [11] introduced a GA based spectrum assignment in CR networks; but spectrum aggregation capability of users is not considered.
For CM2M networks, existing spectrum assignment and aggregation solutions are not applicable directly as practical issues such as Maximum Aggregation Span (MAS) must be taken into account.Furthermore, in aggregation-based spectrum assignment a major challenge is to manage CCI among CM2M devices which is not taken into account in the existing literature.The major contributions of this study are twofold.
(1) To prevent multiple CM2M devices from colliding in the overlapping portions of the spectrum, a centralised approach is applied.Furthermore, an integer optimisation problem to maximise cell throughput is formulated, considering CCI and MAS in an aggregation-aware CM2M network.(2) As the spectrum assignment problem is inherently seen as an NP-hard optimisation problem, evolutionary approaches can be applied to solve this challenging problem.In this article, GA is used to solve the aggregation-aware spectrum assignment because of its simplicity, robustness, and fast convergence of the algorithm [12].
This article is organised as follows.In Section 2, the spectrum assignment and aggregation models are presented.The proposed algorithm is explained in Section 3. Simulation results are discussed in Section 4, followed by conclusions in Section 5.

System Model
2.1.Spectrum Assignment Model.We assume a CM2M network consisting of  CM2M devices defined as Φ = { 1 ,  2 , . . .,   } competing for  nonoverlapping orthogonal channels Γ = { 1 ,  2 , . . .,   } in uplink.All spectrum assignment and access procedures are controlled by a central entity called spectrum broker.We assume that distributed sensing mechanism and measurement conducted by each device is forwarded to the spectrum broker [13].A spectrum occupancy map that is constructed at the spectrum broker and CCI among CM2M devices is determined.Furthermore, the spectrum broker can lease single or multiple channels for   ∈ Φ in a limited geographical region for a certain amount of time.Finally, a base station can transmit data to   in the assigned channels.Figure 2 depicts system model used in this article.
We define the channel availability matrix L = { , |  , ∈ {0, 1}} × as an  ×  binary matrix representing channel availability, where  , = 1 if and only if   is available to   and  , = 0 otherwise.Each   is associated with a set of available channels at its location defined as Γ  ⊂ Γ; that is, Due to the different interference range of each LU (which depends on LU's transmit power and the physical distance) at the location of each CM2M device, Γ  of different CM2M devices may be different [14].According to the sharing agreement, any   ∈ Γ can be reused by a group of CM2M devices in the vicinity defined by Φ  such that Φ  ⊂  In a dynamic environment, channels availability and interference constraint matrix both vary continually; in this study, we assume that spectrum availability is static or varies slowly in each scheduling time slot; that is, all matrices remain constant during the scheduling period.In our proposed solution, a subset of CM2M devices is scheduled during each time slot and the available spectrum is allocated among them without causing interference to LUs.

Spectrum Aggregation Model.
In the traditional spectrum assignment, each channel is composed of a continuous spectrum fragment; thus it is not feasible for users to utilise small spectrum fragments which are smaller than the users bandwidth demand.For instance, assume a CM2M network where every machine requires 4 MHz channel bandwidth, and the available spectrum consists of two spectrum fragments of 4 MHz and four spectrum fragments of 2 MHz (Figure 3).For continuous spectrum allocation, the 2 MHz spectrum fragments cannot be utilised by any machine.Therefore, a continuous spectrum assignment mode can only support two devices for communication (2 × 4 MHz).However, spectrum aggregation-enabled device can exploit fragmented segments of the spectrum by using specialised air interface techniques, such as DOFDM.In Figure 3, if a number of small spectrum fragments are aggregated into a wider channel, then 16 MHz of unused spectrum is available to support four CM2M devices (4 × 4 MHz).
Due to the limited aggregation capabilities of the RF front-end, only channels that reside within a range of MAS can be aggregated.With this constraint, some spectrum fragments may not be aggregated because their span is larger than MAS.Our proposed algorithm takes MAS into consideration.For the sake of simplicity, we make following assumptions.
(1) All CM2M devices have the same aggregation capability (i.e., MAS for all devices is the same).( 2) Guard band between adjacent channels is neglected.
(3) Bandwidth requirement of each device and bandwidth of each channel are an integer multiple of subchannel bandwidth Δ, which is the smallest unit of bandwidth (in fact, the smaller fragments would demand excessive filtering to limit adjacent channel interference); that is, where N is the set of natural numbers,   is the number of requested subchannels by   ,   is the number of subchannels in   , and BW  is the bandwidth of   .
The total available spectrum (i.e.,  channels) is subdivided into multiple number of subchannels.If the available spectrum band consists of C subchannels (i.e., total available bandwidth is C ⋅ Δ), then where 1 ≤  ≤ , where   has   subchannels and γ, represents the th subchannel of   .Each γ, can be represented in an interval defined as [F  , , F  , ], where F  , and F  , are the lowest and highest frequency of γ, : (3) Based on this new subchannel indexing, matrices L and C can be rewritten as Aggregating spectrum

Available spectrum
Unavailable spectrum where c represents index of each subchannel within the available spectrum.
The subchannel assignment matrix A = { ,c |  ,c ∈ {0, 1}} ×C is an ×C binary matrix representing subchannels assigned to CM2M devices for aggregation such that  ,c = 1 if and only if subchannel c is available to   and 0 otherwise.We define the reward vector B = {  = Δ ⋅ ∑ C c  ,c } ×1 to represent total bandwidth that is allocated to each CM2M device during scheduling time period for a given subchannel assignment.

Problem Formulation
3.1.Optimisation Problem.One of the key objectives of the deployment of CM2M network is to enhance the spectrum utilisation.To consider this crucial goal, we define network utilisation to maximise the total bandwidth that is assigned to CM2M devices and referred to as Maximising Sum of Reward (MSR): To maximise MSR the spectrum aggregation problem can be defined as a constrained optimisation problem as follows: subject to Expression (8) (9) guarantees that the range of allocated spectrum is equal to or less than MAS.A must satisfy the interference constraints ( 10) and (11); expressions (10) and (11) guarantee that there is no harmful interference to LUs and other CM2M devices, respectively.

Spectrum Aggregation Algorithm Based on Genetic
Algorithm.Traditionally the spectrum assignment problem has been classified as an NP-hard problem [12].Herein, GA is employed to solve the aggregation-based spectrum assignment problem in order to obtain faster convergence.GA is a stochastic search method that mimics the process of natural evolution.In addition, it is easy to encode solutions of spectrum assignment problem to chromosomes in GA and compare the fitness value of each solution.The specific operations of the proposed algorithm, referred to as MSR Algorithm (MSRA), can be described through the following steps.
(1) Encoding.In MSRA, a chromosome represents a possible conflict-free subchannel assignment.In order to decrease search space (by reducing redundancy in the data) and obtain faster solutions, similar approach as described in [12] is adopted in this article.We apply a mapping process between A and the chromosomes, based on the characteristics of L * and C * .Only those elements of A are encoded whose corresponding elements in L * take the value of 1; that is,  ,c = 0, where (, c) satisfies  * ,c = 0.As a result of this mapping, the chromosome length is equal to the number of nonzero elements of L * and the search space is greatly reduced.Based on a given L * , length of the chromosome can be calculated as ∑  =1 ∑ C =1  * , .(2) Initialisation.During initialisation process, the initial population is randomly generated based on a binary coding mechanism as applied in [12].The size of the population depends on |Φ| and |Γ|; for larger |Φ| and |Γ|, population size should be increased, where | ⋅ | indicates cardinality of a set.
(3) Selection.The fitness value of each individual of the current population according to MSRA criteria defined in ( 6) is computed.According to the individuals fitness value, excellent individuals are selected and remain in the next generation.The chromosome with largest fitness value replaces the one with a small fitness value by the selection process.
(4) Genetic Operators.To maintain high fitness values of all chromosomes in a successive population, the crossover and mutation operators are applied.Two randomly selected chromosomes are chosen in each iteration as the parents and the crossover of the parent chromosomes is carried out at probability of crossover rate.In addition to selection and crossover operations, mutation at certain mutation rate is performed to maintain genetic diversity.
(5) Termination.The stop criteria of GA are checked in each iteration.If they can not be satisfied, step (3) and step (4) are repeated.The number of maximum iterations and the difference of fitness value are used as the criteria to determine the termination of GA.
The population of chromosomes generated after initialisation, selection, crossover, and mutation may not satisfy the given constraints defined in ( 8)- (11).To find feasible chromosomes that satisfy all constraints, a constraint-free process is applied that has the following steps (in order).
(1) Bandwidth Requirements.The vector B as given in Section 2.2 is calculated.  should be equal to either   or zero; otherwise all genomes related to   are changed to zero.
(2) MAS.To satisfy the hardware limitations of the transceiver, expression (9) should be satisfied; otherwise all genomes related to   are changed to zero.
(3) No Interference to LUs.Expression (10) guarantees that CM2M devices transmissions do not interfere LUs transmissions, ensuring that CM2M network does not harm LUs performance.If expression (10) is not satisfied, all genomes related to   are changed to zero.
(4) CCI.Expression (11) guarantees that there is no harmful interference to other CM2M devices.If expression (11) is not satisfied, one of two conflicted devices is chosen at random, and then all genomes of the selected device are changed to zero.
To achieve higher spectrum utilisation and faster convergence, after each generation, MSRA assigns all unassigned spectra to remaining CM2M devices randomly, whenever possible.At the same time, MSRA guarantees that all the constraints defined in ( 8)-( 11) are satisfied at all time.

Simulation Results
In this section, a set of system-level performance results are presented in order to compare and show the efficiency of MSRA over MSA [10], AASAA [9], and RCAA.The simulation results demonstrate high potential of the proposed Due to the random nature of the channel bandwidth and the devices bandwidth demand, Monte Carlo simulations are performed and each simulation scenario is repeated 100,000 times.The default parameters used in the simulations are listed in Table 1, where (1, 20) represents the discrete uniform random integer numbers between 1 and 20.Each of the channels is modeled as flat Rayleigh channel with path loss model of PL = 128.1 + 37.6 log 10  ( is in km) and penetration loss of 20 dB.The mean and standard deviation of log-normal fading are zero and 8 dB, respectively.In our simulation model, the CM2M devices located randomly without restrictions within a rectangular area of 2 km × 1 km.All channels are randomly selected between 54 MHz and 806 MHz television frequencies (channels 2-69).Typically, the number of M2M devices is very high in each cell, but in this study, because of high computational complexity of SOTA solutions, smaller number of M2M devices is considered for comparison purposes.
To investigate the simulation results effectively, the following terms are defined and used in our analysis.
(1) Spectrum Utilisation.It is referred to as U which is defined as the ratio of the sum of rewarded bandwidth to the sum of all available bandwidths; that is, (2) Network Load.It is referred to as L which is defined as the ratio of the sum of all CM2M devices bandwidth requirements to the sum of all available bandwidths; that is, (3) Number of Rejected Devices.Rejected devices are those machines that are not assigned any spectrum in a certain scheduling time slot.
4.1.Scenario I: Without CCI.In this scenario, the performance of MSRA is compared with the SOTA algorithms including MSA [10], AASAA [9], and RCAA when among CM2M devices is not considered.Therefore, we assume that CM2M devices transmissions do not overlap with the transmission of other CM2M devices using the same channel.
For  = 30, L increases by increasing the number of CM2M devices from 5 to 60. Figure 4 shows that when the number of CM2M devices increases, the spectrum utilisation also increases in all three methods, but MSRA utilises all available whitespaces in various network loading conditions more efficiently than MSA, AASAA, and RCAA.This can be explained by the fact that, in case of higher L, network can allocate better segments of spectrum to users because of higher multiuser diversity.In addition, because of using stochastic search method, MSRA achieves near to optimum solution in comparison to other SOTA solutions which are based on approximate algorithms.For MSRA, when L is higher than 3, CM2M network becomes saturated due to the lack of available spectrum.However, for the rest of the methods, there are still unassigned spectrum slices.

Scenario II: With CCI.
In this scenario, CCI exists among CM2M devices and we compare our algorithm, MSRA, with AASAA and RCAA.As MSA inherently does not consider CCI, for that reason, we do not include MSA for comparison.Figure 5 shows the spectrum utilisation, according to different network loads by increasing the number of CM2M devices from 5 to 55 when there are only seven available channels (i.e.,  = 7).As shown in Figure 5, MSRA outperforms AASAA and RCAA for different network loads.Similar to Scenario I, MSRA utilises TVWS even better than previous scenario, because some CM2M devices in network may reuse spectrum that is used by other devices in CM2M network.
Figure 6 represents the number of rejected CM2M devices when the network load increases.The number of rejected CM2M devices increases with the network load; MSRA has fewer numbers of rejected CM2M devices (or more satisfied devices) than AASAA and RCAA of different network loads.MSRA optimises spectrum utilisation by admitting devices with better channel quality to the network and allocates the spectrum resources effectively.Furthermore, MSRA does not assign any spectrum resources to the devices that has least contribution to overall network throughput.Figure 6 implies that MSRA increases the capacity of network (which is very vital for M2M networks because of a very large number of devices).Our approach may starve some of devices which are located far from the base station; in our future work, we will optimise network performance based on proportional fairness objective function to guarantee the fairness among devices.

Convergence of MSRA.
Because of the nature of genetic programming, it is arguably impossible to make formal guarantees about the number of fitness evaluations needed for an algorithm to find an optimal solution.However, herein, computer experiments are performed to show the impact of  the number of generations on the performance of MSRA.
The system parameters used in the section for simulation are listed in Table 2.For the purpose of convergence studies, we assume  = 200 and  = 10.
Figure 7 shows the best fitness value (MSRA) for a population in a different number of generations.As shown in Figure 7, the performance of algorithm is enhanced, when the number of generations increases; however this is at the cost of increased processing time.After roughly 34 generations, the fitness value saturates at optimal value which shows the effectiveness of using GA for spectrum assignment using spectrum aggregation.
Moreover, Figure 8 illustrates distribution of processing time for MSRA to find an optimal solution.As shown in Figure 8, at 85% of time, MSRA finds an optimum solution in less than scheduling time slot (1 ms) and 15% takes more than scheduling time slot.Additionally, MSRA can be optimised to use fewer processor resources, so that it can execute more rapidly.
Furthermore, Lobo et al. [15] provided a theoretical and empirical analysis of the time complexity of traditional simple GAs.According to [15], GA has time complexities of O(∑  =1 ∑ C =1  * , ) which is dependent on length of each chromosome.The linear time complexity for GA occurs because the population sizing grows with the square root of chromosome length.The time complexity presented herein is for the worst-case scenario when the population size is assumed to be fixed and maximum of rest of generations.

Figure 4 :
Figure 4: The impact of varying network load conditions on spectrum utilisation (scenario I: without CCI).

Figure 5 :
Figure 5: The impact of varying network load conditions on spectrum utilisation (scenario II: with CCI).

Figure 6 :
Figure 6: The impact of varying network load conditions on the number of rejected CM2M devices (scenario II: with CCI).

Figure 7 :FrequencyFigure 8 :
Figure 7: The impact of the number of generations on MSRA results.
assures that rewarded bandwidth   to each accepted   must be equal to   's bandwidth demand   ; if CM2M network cannot satisfy   's bandwidth request,   is rejected and   = 0.If F  , (1 ≤  ≤   and 1 ≤  ≤ ) is the lowest frequency of an initial aggregated subchannel and F  , (1 ≤  ≤   and 1 ≤ t ≤ ) is the highest frequency of a terminative subchannel,