I Introduction
A key challenge for 5G and beyond5G deployments will be, without doubt, a seamless integration of Internet of Things (IoT) services, as these are expected to span a very broad range of application scenarios and technical requirements. The envisioned high IoT device diversity and density make the problem even more challenging. It is also expected that large antenna arrays, massive MIMO, and the centralized radio access network (CRAN) architecture will be an integral part of future deployments and can be exploited to enable efficient network access to IoT devices.
While future networks may need to accommodate humongous numbers of IoT devices, the communication requirements of most IoT devices are expected to be intermittent and sporadic. As a result, random access (RA) protocols are ideally suited for serving IoT devices. Most existing RA works consider cellular RA, i.e., access to a singlesite [1, 2]. RA protocols for CRAN architectures, however, are also gaining attention [3, 4].
Recent advances in cellular RA include coded and slotted versions of ALOHA, and compressed sensing (CS) techniques aiming to exploit the inherent sparsity in user activity (see e.g., [1]). Some of the benefits offered by massive MIMO have already been exploited in cellular RA, and include spatial collision resolution (see e.g., [2]).
Several CRAN RA methods have emerged recently. [3] leverages CS for sparse user detection, and clustering of nearby base stations (BS) to reduce computational complexity. [4] also exploits CS techniques for user detection, but, unlike [3]
, it performs user detection and channel estimation jointly by modeling a twodimensional sparsity, i.e., sparsity in both device activity and devicesite connectivity. While the twodimensional sparsity model is wellsuited to CRAN architectures with IoT devices, the computational complexity in
[4] scales up with the number of antennas at both BSs and users.Contentionbased cellular RA is preferred in LTE as a means for accommodating IoT devices with sparse and intermittent activity [5]. Uplink (UL) pilot collisions, detected by a BS, are resolved using additional message exchanges between BS and users [5]. The inherent delays in resolving collisions reduce the appeal of these schemes, especially in the presence of networks with small cells and high IoT densities.
In this work, we consider CRAN RA protocols that enable high detection rates of active devices in the network and study the problem of learning the largescale channel gains of the links between detected devices and CRAN sites. Such link information can be leveraged for scheduling, loadbalancing, and interference suppression to greatly improve performance in dense massive MIMO deployments [6], [7].
We focus on a usercentric architecture of the type introduced in [8] that leverages dense deployments of largearray remote radio heads (RRHs) and virtual sectorization. This architecture is ideal for IoT RA, as it allows onthefly lowlatency detection of active devices across the CRAN coverage area. It also greatly outperforms conventional cellular RA in terms of density of simultaneously detected devices [8, 9].
In the context of these RA protocols, we present techniques for learning the strengths of links between detected devices and RRH sectors in the network. By jointly processing the individual RRH sector devicedetection reports via a framework borrowed from recommendation systems [10], the proposed methods are able to effectively classify the strengths of links between detected devices and RRH sectors not directly available in the detection reports. Our proposed schemes only exploit prior knowledge of the distances between RRH sites, but no knowledge of activedevice locations.
The paper is organized as follows. Sec. II describes the CRAN RA schemes of interest. The problem of learning the strengths of as many links as possible is considered in detail in Sec. III. Sec. IV presents a brief simulationbased evaluation of the link classification schemes proposed in Sec. III. Finally Sec. V provides concluding remarks and directions for future investigations.
Ii CRAN Random Access
In this section we present the RA schemes of interest in this work. We consider a CRAN of large antennaarray BSs spanning a wide geographical area and serving a very large number of IoT devices, each with a single antenna. Similar to [11, 8], spatial processing is applied at each BS to create virtual sectors, for a total of sectors in the network. The set of all BS sectors is indexed via the index , where is the BS index, and .
We consider a slotted RA scheme, according to which blocks of timefrequency resources, referred to as RA blocks, are reserved for RA [5]. For convenience, we consider a timeline where RA blocks are periodically reserved. A single frame is formed by one RA block and the following block of resources dedicated to serving the active devices detected during the RA block in the frame. The RA block within the th frame is referred to as the th RA block. See Fig. 1.
We model the device accessrequest arrival process as a Poisson process with a rate of device requests between consecutive RA blocks (over the geographical area spanned by the network). Any device with an “arriving” access request waits for the next RA block to access the network and remains “active” until it is detected. It is served in the serving block following the RAblock it is detected. We also let denote the active users in the th frame.
During the th RA block, active devices in broadcast uniquely distinguishable uplink pilots and a subset of them, , is detected by the network. In the CRAN RA scheme we consider, an active user is considered to be detected by the network as long as at least one sector detects its unique ID during RA.
Much like any RA protocol, not all active devices in an RA block are detected. The subset of active devices, that are not detected are “queued”, i.e., continue to be active in later RAblocks. Thus, the active devices in frame , either arrived in frame , or in an earlier frame but have not yet been detected.
An important property of a good RA protocol is keeping stable queue sizes. By letting denote the device detection rate, it is important that the RA scheme is stable with respect to , that is, . This implies that the active user rate per frame, does not diverge.
For a stable RA protocol, the ratio is inversely proportional to the expected delay (in RA frames) that a device waits before it is detected. A desired average delay of can be achieved if, (and ) is greater than . Stable RA protocols that yield short expected detection delays while keeping RAblock overheads low are highly desirable.^{1}^{1}1In this work, the RAblock overheads are proportional to the number of randomaccess slots, or “pilot dimensions”, reserved in an RAblock.
At this stage it is worth contrasting the slotted RA scheme in consideration against LTEPRACH. LTE employs a form of slotted RA that can be modeled at an abstract level via the schematic in Fig. 1. However, each device in LTE is associated with a single BS, and therefore, PRACH is designed for RA between devices and a single BS. To enable such cellular RA, resource reuse is employed within the RA block in LTE, so that the RAblock pilot dimensions used within a cell are not available in nearby cells. That loss in overhead efficiency in LTE is well known. In practice, it is alleviated by using partially overlapped pilot resources in nearby cells, which inherently results in intercell RAblock interference. An additional difference arises from the fact that, in LTE, collisions are resolved via a sequence of message exchanges between a BS and its colliding devices. In contrast, the CRAN RA schemes we consider rely on instantaneous user identification and collision detection, thereby resulting in lower detection delays and lower protocol complexity.
The proposed CRAN slotted RA enables very flexible operation in the serving block as it enables instantaneous association of the detected devices to one or multiple BS sectors in the network. Efficient network operation, whether cellular, CoMP (coordinated multipoint), or cellfree, can be enabled if channel conditions between detected devices and nearby sectors are available to the network. As an example, in massive MIMO deployments, network performance can be greatly improved via optimized scheduling, loadbalancing, and interference suppression, if the largescale channel gains (link gains) between devices and BS sectors are available [7, 6]. In Sec. IIA, we describe the specific CRAN RA protocol we consider for detecting devices. It is based on an adaptation of coded pilot designs of [9] for RA.^{2}^{2}2The collision detection capability of these codes for RA is also advocated in [2]. Assuming the RA protocol of Sec. IIA, in Sec. III we consider the problem of effectively classifying the strengths of links between detected devices and BS sectors in an effort to enable efficient CRAN operation during the serving blocks in each frame.
Iia RA protocol based on coded pilot designs of [9]
In this section, we describe the protocol according to which active devices engage in RA along with the mechanisms used at the network side to detect active devices. We focus on a fixed but arbitrary RA frame, , and, for convenience, omit the dependence of all parameters and variables on the RA frame index, . We assume that the link gain between any sector and any device , , stays constant for the duration of the frame. We classify the link between sector and device , as strong when the link strength exceeds a predetermined (operator chosen) , and as weak otherwise.
In the protocol we consider, each RA block comprises RA pilot slots during which each active device in
can access the medium. In any pilot slot, each active device accesses the medium with probability
. In particular, the pilot activity of device in pilot slot is captured by a Bernoulli(. The ’s are statistically independent in and . The RA overhead parameter and the access probability are chosen by the network operator, and are assumed fixed and known to all the devices in the network.The RA protocol we consider leverages the UL pilot designs considered in [9] together with largeantenna arrays at each BS. Each device is a priori statically assigned a unique pilot codeword that the device uses in any slot that it sends a pilot. By appropriately processing its pilotslot observations through its large array, each BS is able to obtain the IDs of any active device that is the only one with a strong link to one of its sectors among all the active devices broadcasting a pilot in a given pilot slot [9]. In particular, the combination of the codes in [9] and the BS processing enable a BS sector to gather the following information during a fixed but arbitrary pilot slot:

when multiple active devices with strong links to the BS sector transmit pilots in the slot, the BS sector detects a collision;

when no active devices with strong links to the BS sector transmit pilots in the slot, the BS sector detects silence;

when a single active device with a strong link to the BS sector transmits a pilot in the slot, the sector detects the user ID and obtains an accurate estimate of the link gain.
Each BS sector provides the information it has gathered through its pilot slots to the central unit. By merging this information the central unit is able to obtain , the set of detected active devices across the whole network. It is also able to obtain additional valuable information, including: (i) : the set of sectors that detected silence on pilot dimension ; (ii) : the set of sectors that detected the th active device (assuming all active devices in are reindexed from 1 to ); (iii) : the set of sectors that on pilot dimension detected a device different from ; (iv) link gains between BS sectors and the active devices they detected.
Iii Classifying the active device links
In this section we focus on classifying the strengths of as many links between active devices and BS sectors as possible based on the information provided in the RA reports of all the BS sectors. In particular, focusing on a fixed but arbitrary RA frame, the goal is to classify at the central unit the links between the BS sectors and the detected devices into two classes: strong and weak. Each BS sector’s RA report contains entries, one per pilot dimension. In accordance to the RA protocol in Sec. II, for each , each sector reports one of the following: (i) a pilot collision; (ii) pilot silence; (iii) device detection, along with the identity of the detected device and the associated link gain.
In Sec. IIIA we formulate the strong link detection problem of interest. In Sec. IIIB we present a baseline scheme that classifies link strengths by processing the RA report information per detected device (i.e., individually). Specifically, the central unit processes the sector RA reports separately per device, and, in the process, infers all the information it can regarding links between the given device and BSsectors. In Sec. IIIC, we develop a class of methods motivated by what are referred to as recommendation systems. These online methods jointly classify the links of all detected devices across the network, exploiting in the process the inherent spatial correlation in the link strengths of nearby detected devices without making use of device location information.
Iiia Problem Formulation
We focus on the problem of classifying the strengths of all links in the network (i.e., between all detected devices and all BS sectors) as strong and weak, during a fixed but arbitrary RA frame and omit the dependence of all parameters and variables on the RA frame index, . The link gains between the BS sectors and the set of detected devices in the frame are compactly represented via a matrix . By letting denote the ID of the th detected device in the frame, the th entry of contains the link gain between the th BS sector and the th detected device, i.e., .
Links are classified as weak or strong based on the comparison against a predefined classification threshold . To this end, we refer to the matrix , with th entry
(1) 
as the hypothesis matrix. The task of interest amounts to forming a classifier at the central control unit using the BS sector RA reports, so that is as close to as possible.
The classification performance is measured in terms of probability of detection and probability of false alarm, viz.,
(2a)  
(2b) 
Our performance evaluations rely on estimates of the quantities in (2) obtained over RA frame realizations with large :
(3a)  
(3b) 
and where we made the dependence of and on the RA frame index, , explicit.
IiiB Baseline Scheme
In this section, we present a baseline scheme, which obtains an estimate by processing the RA reports separately per detected device. The baseline classifier exploits the fact that some entries of are directly available or can be indirectly inferred from the BSsector RA reports. Let denote the subset of entries of the matrix (i.e., subset of links) where the value is known to the central controller. Also, let denote the complement of , i.e., the set comprising of all entries for which is not known to the controller.
Naturally, the baseline scheme sets for all . Furthermore, the baseline scheme applies a randomized decision to determine the value of for each . In particular, for each the baseline scheme flips a biased coin whose probability of heads equals . The scheme sets if the coinflip outcome is heads, and sets , otherwise. In summary,
(4) 
and where denotes the biased coinflip outcome associated with entry .
Due to the randomization in the baseline scheme, its Receiver Operating Characteristic (ROC) curve is piecewise linear. With , all links in are classified as strong, yielding the point in the ROC with the maximum false alarm and detection rates . For the scheme with , all links in are classified as weak, yielding the point in the ROC with the minimum false alarm and detection rates . Varying the value of in the range yields all the points on the line segment connecting and .
We next focus on determining the set , or equivalently the set of links that can be obtained from the BSsector RA reports. It is convenient to express as , with () denoting the sets of locations of the 1’s (0’s) in that can be obtained from these reports.
Obtaining the known strong links from the RA reports: Some entries of are directly available from the RA reports. Focusing, in particular, on the th detected device, for an arbitrary but fixed index , we note that is available to the controller for all sectors that detected device . By letting
with denoting the BS sectors that have detected the th device, all the 1 entries in with coordinates in are directly available at (known to) the central controller.
Inferring zeros: By preprocessing the RA report per device, , the baseline scheme is also able to infer additional entries of , and, in particular, zero entries (i.e, weak links).
First, consider any pilot dimension during which the th detected device has been detected by a sector in the network: Any sectors that detected silence or a different device (from ) in the same pilot dimension must have a weak link to the th device. This implies that the following set of entries (all along the th column) of are known to be zero:
(5) 
where is the set of pilot slots on which the th (detected) device has been detected by at least one sector in the network.
Secondly, the knowledge of the BS geographical locations can be exploited to determine some additional weak links across the network. For two BSs which are sufficiently far apart, no device can have strong links simultaneously to both BSs. Let denote the distance beyond which two BSs are sufficiently far from one another so that no device can simultaneously have strong links to both BSs.^{3}^{3}3The distance can be experimentally obtained. Note that with this approach a tiny fraction of the “known” entries of are misclassified, that is, , while . Hence, although, in principle, we should have , in our experiments is slightly less than 1. Using knowledge of the distances between any pair of BSs in the network, the following set of entries (all along the th column) of are also known to be zero:
(6) 
where denotes the BS sectors that have detected the th device, and where denotes the set of all BS sectors whose distance from BS sector exceeds .
In summary, the set of all known zero entries of is given by , where .
IiiC Matrix Completion Method
In this section we consider an online method for improving upon the stronglink prediction performance provided by the baseline scheme. The online method makes use of the intuitive spatialconsistency assumption that there is an underlying lowdimensional subspace describing the links between BS sectors and devices in the network. Consequently, it uses the set of strong and weak links over the set of (BSsector, device) pairs (that are available to the central controller via the RA reports) as data to determine the “best” lowdimensional model and uses this model to predict and classify links in the complement set of (BSsector, device) pairs.
The online methods we consider fall within the class of matrix completion problems, which have recently gained attention in a broad range of applications. One important application is found in recommendation systems, and includes the “Netflix Challenge”, according to which, movie recommendations are made to customers, based on their previous ratings and otheruser ratings [10]. At its core, this problem amounts to estimating the missing entries of a rating matrix where the matrix dimensions are the number of movies and the number of customers. While some entries of the matrix are already known, as some customers have rated some of the movies, rating estimates are required for some moviecustomer pairs. Since the Netflix Challenge, various solutions have been devised, which take into account various practical aspects of the problem [12]
. Recommendation systems research is very active, with many stateoftheart works exploiting deep learning algorithms. As a proof of concept, in this paper we exploit a vanilla solution based on matrix factorization (collaborative filtering with latent factor models)
[10]. However, more advanced techniques can be also used and can, in principle, yield additional performance benefits.In adapting the matrix completion problem to our setting, we rely on the following analogy: BS sectors are analogous to movies, devices are analogous to customers, and link gains between devices and sectors are analogous to movie ratings.
The problem of interest can be readily mapped into an intermediate matrix completion problem by noting that: (i) the central unit has available for ; (ii) the fact that for . As a result, we can formulate the following matrix completion problem:
(7)  
for some appropriately chosen to optimize performance.^{4}^{4}4Practical aspects of the optimization problem, regarding rank selection, regularization etc. are discussed in detail in Appendix A.
The solution to (7) can be used to classify links not in . Indeed, by letting , a hypothesis test can be employed for all , of the form
(8) 
whereby is chosen so that the falsealarm probability does not exceed a preassigned value .
Varying between and can control the detection and false alarm rates, similar to the way the value of controls these rates in the randomized baseline decision scheme. A block diagram for the proposed onlinemethod is shown in Fig. 2.
Iv Numerical Results
In this section, we compare the performance of the matrix completion based sectordevice link classification against that of the baseline scheme in terms of detection and false alarm rates as defined in Sec. III.
We consider a network layout involving
BSs and large numbers of devices, both uniformly distributed over a square geographical area of size
. Each BS has sectors. The four sector orientation is chosen randomly and independently per BS. Fig. 3 shows such a sample layout: BSs are shown with squares and devices are shown with dots (only devices are shown for representational clarity).The link strengths between BSsectors and devices are determined via a generalization of the onereflection pathloss model in [8] (described in detail Appendix B. Fig. 4 provides a plot of the strengths of all the links between the BSsectors and the devices in Fig. 3 as a function of link distance. The horizontal line in the figure, set at dB, reflects the a priori networkdefined threshold, separating strong links from weak ones. As the figure reveals, links at similar distances can yield a broad range of link strengths. At the same time, the probability that a given link is strong decreases as a function of the distance.
Our comparisons are based on the use of the RA protocol in Sec. IIA, and are tuned such that an arrival rate of devices can be accommodated with an expected delay requirement of , where . The RAblock within each frame is assumed to have an overhead of pilot dimensions, where denotes the relative overheads of the RA scheme. Given and , the target operating point is and . With this configuration 99% of the active devices are detected on average during the RA block of a single frame. Implicit in this operatingpoint specification, however, is the choice of the parameter, which, given a particular overhead, guarantees , with delay requirements .
In principle the value of the parameter can be determined and tuned via simulations. In particular, as discussed in Sec. II, in order to stabilize the queues and, in addition, meet the delay requirements, the scheme should guarantee the following:
(9) 
Hence, given and , we can estimate and as a function of via simulations, and subsequently use the value that satisfies the constraints in (9).
As an alternative to extensive simulations, we develop a simple method for predicting the value of yielding (9) and investigate its match in the parameter range of interest based on simulations. The method is based on modeling the device detection process at each sector. This model assumes that, at any given slot out of the slots in the RA training block, each active user not yet detected has a probability equal to some value to be detected in that slot (independent of all other users). Based on this rudimentary model, we can derive the following ruleofthumb formula between and :
(10) 
For the network in Fig. 3, Fig. 5 shows a simulationbased evaluation of the efficacy of the formula in (10) using and , as a function of and . As the figure reveals, the value predicted by the model matches simulations especially well for the values of that maximize and the values of that can reach . As a result, for any given operating point of interest given by a triplet of (, , ) values, we set .
Fig. 6 shows the probability of detection vs. false alarm of strong links for various values. For each value, the dotted line depicts the performance of the baseline scheme, while the associated dashed curve corresponds to the matrixcompletion online method. The red solid curve in the figure shows the performance of the matrixcompletion method at , operating over a window of the two last frames (current and most recent past frame) in predicting the links of currentframe active devices. As the figure reveals, the online methods yield sizable advantages in predicting strong links. Furthermore, the performance improves when the online method is applied over a window of the two most recent frames.
V Conclusion
We have studied the efficacy of matrixcompletion algorithms in classifying the strengths of links between detected devices and RRH sectors, based on CRAN RA observations collected across a widearea network. Our analysis and simulations reveal that such purely datadriven schemes can significantly improve link classification. More important, link classification becomes more accurate when the proposed methods are applied over windows of consecutive RA frames, allowing the algorithm to work on larger numbers of devices.
Our work indicates the existence of an efficient representation for the collection of linkstrength device vectors in terms of a lowdimensional “feature” subspace. One can foresee that, with lots of IoT devices, these feature vectors can be accurately learned and subsequently used to make link predictions per device. Such techniques are left for future work.
a Notes on the practical issues of matrix completion problem
For simplicity, we let in (7) and add regularization terms:
(11)  
where denotes the Frobenius norm. In (11), the first two terms indicate the sum square of estimation error, while the third term represents the regularization factor to reduce the overfitting effect [13]. Although (11) is not jointly convex with and , it is a convex quadratic function of one variable if we fix others. Therefore, a simple block descent algorithm can be considered to achieve a local optimal solution:

Following the first order optimality condition, find the update equations of and , respectively.

With a random initialization, iteratively evaluate variables by update equations until the estimation error converges or the maximum number of iterations is reached.
To achieve a suboptimal solution of the matrix factorization problem and estimate unknown links, we fix the parameter settings in Tab. I empirically,
Step size of gradient descent  

Maximum no. of iterations  
where is the parameter to lower bound the normalized sum square of estimation error over the training data (known link entries):
(12) 
Therefore, either (12) is satisfied or the number of iterations reaches its maximum, the algorithm stops and outputs the result of matrix factorization.
B Channel model
We consider a radio access network (RAN) with RRHs and UEs. Each RRH is equipped with an Melement uniform linear array (ULA), while each UE is equipped with a single antenna. With the assumption of the block fading channel [14], the channel response vector between the th RRH and the th UE, i.e., , can be represented as
(13) 
which is a superposition of MPCs. In (13), is the carrier frequency, indicates the number of MPCs between and , and is the large scale loss for the th MPC of link to.
reflects the small scale fading following complex Gaussian distribution, and
is the propagation time of MPC. represents the direction of departure (DOD) from the perspective of RRH, and is the steering vector with respect to :(14) 
where and indicate the antenna spacing and wavelength, respectively. Note that we only consider the twodimensional channel model in the azimuth domain in (13), but it can be easily generalized to incorporate the elevation domain.
Within the stationarity region of the second order channel statistics, longterm channel parameters, including DODs of MPCs and large scale loss remain approximately the same for a relatively long period of time, which can be equivalent to hundreds of coherence blocks [15]. On the other hand, small scale variables, e.g., , remain the same within the coherence time and coherence bandwidth, but vary independently across different blocks. Therefore, based on the assumption of uncorrelated scattering, similar to [16], we can achieve the closeform expression of channel covariance as a function of longterm large scale parameters:
(15) 
We consider the DFT based virtual sectorization within each RRH, where each sector is formed by spatially filtering the ULA signal through a contiguous set of sectorspecific DFT beams. Define as the normalized DFT matrix (each column has unit norm), we acquire the average sector gain between the th sector of and as
(16) 
where indicates the sector gain matrix, is the number of sectors per RRH, and denotes set of DFT beam indices for the th sector of . Splitting all DFT beam tones equally among sectors of an RRH, we have .
To generate synthetic channel profiles, we implement a geometric stochastic channel model (GSCM) [8], which captures the propagation characteristics of direct LOS path and nondirect reflected path. Since MPCs interacting with multiple objects will usually suffer severe propagation loss, especially at high frequencies [17], we only consider the nondirect path with singlebounce scattering. Blockage effects [18] are also incorporated by adding circleshaped blockers so that the existence of direct path is determined geometrically.
For the attenuation of the direct path, we utilize a smooth transition function [19]:
(17) 
where denotes the separation distance between two link ends, and and indicate the breakpoint distance and pathloss exponent, respectively. Regarding to the nondirect paths, we combine the separate loss from RRHtoscatterer and scatterertoUE:
(18) 
which is commonly used for modeling reflected path attenuation [14]. In (18), is the attenuation coefficient of NLOS path, and and denote the distance of UE toscatterer and scatteretoRRH, respectively. More discussions on the validity of (18) can be found in [8].
In conclusion, we briefly summarize the procedures of channel modeling below:

In the operating area, randomly drop RRHs, UEs, scatterers, and blockers satisfying separate Poisson point process (PPP).

For the direct path of a particular link, if it is blocked by any blocker, there is no LOS path for that link.

Remove weak NLOS paths, whose attenuation is below than a threshold .
Since scatterers of a link are usually not too far away from either the RRH or the UE, the last step will eliminate nonrealistic paths, and maintain the set of effective scatterer for every link, respectively. Note that this GSCM model provides an implementation of the “common scatterer” concept used, e.g., in COST 2100 [20].
References
 [1] V. Boljanovic, D. Vukobratovic, P. Popovski, and C. Stefanovic, “User activity detection in massive random access: Compressed sensing vs. coded slotted ALOHA,” arXiv preprint arXiv:1706.09918, 2017.
 [2] E. D. Carvalho, E. Bjornson, J. H. Sorensen, P. Popovski, and E. G. Larsson, “Random access protocols for massive MIMO,” IEEE Commun. Mag., vol. 55, no. 5, pp. 216–222, May 2017.
 [3] X. Xu, X. Rao, and V. K. N. Lau, “Active user detection and channel estimation in uplink CRAN systems,” in IEEE Intern. Conf. on Commun. (ICC), Jun. 2015, pp. 2727–2732.
 [4] Q. He, T. Q. S. Quek, Z. Chen, and S. Li, “Compressive channel estimation and multiuser detection in CRAN,” arXiv preprint arXiv:1702.0638, 2017.
 [5] M. Hasan, E. Hossain, and D. Niyato, “Random access for machinetomachine communication in LTEadvanced networks: issues and approaches,” IEEE Commun. Mag., vol. 51, no. 6, pp. 86–93, Jun. 2013.
 [6] D. Bethanabhotla, O. Y. Bursalioglu, H. Papadopoulos, and G. Caire, “Optimal usercell association for massive MIMO wireless networks,” IEEE Trans. on Wireless Commun., vol. 15, no. 3, pp. 1835–1850, Mar. 2016.
 [7] Q. Ye, O. Y. Bursalioglu, H. Papadopoulos, C. Caramanis, and J. G. Andrews, “User association and interference management in massive MIMO HetNets,” IEEE Trans. on Commun., vol. 64, no. 5, pp. 2049–2065, May 2016.
 [8] O. Y. Bursalioglu, C. Wang, H. Papadopoulos, and G. Caire, “A novel alternative to cloud RAN for throughput densification: Coded pilots and fast userpacket scheduling at remote radio heads,” in 50th Asilomar Conf. on Signals, Systems and Computers, Nov. 2016.
 [9] O. Y. Bursalioglu, C. Wang, H. Papadopoulos, and C. Caire, “RRH based massive MIMO with “on the fly” pilot contamination control,” in 2016 IEEE Intern. Conf. on Commun. (ICC), May 2016.
 [10] Y. Koren, R. Bell, and C. Volinsky, “Matrix factorization techniques for recommender systems,” Computer, vol. 42, no. 8, pp. 30–37, Aug. 2009.
 [11] Z. Li, N. Rupasinghe, O. Y. Bursalioglu, C. Wang, H. Papadopoulos, and G. Caire, “Directional training and fast sectorbased processing schemes for mmwave channels,” in IEEE Intern. Conf. on Commun. (ICC). IEEE, May 2017, pp. 1–7.
 [12] B. Smith and G. Linden, “Two decades of recommender systems at amazon.com,” IEEE Internet Computing, vol. 21, no. 3, pp. 12–18, May 2017.

[13]
A. Ng,
Recommender Systems, Machine Learning
. Coursera.  [14] D. Tse and P. Viswanath, Fundamentals of Wireless Communication. New York, NY: Cambridge University Press, 2005.
 [15] A. Adhikary, J. Nam, J.Y. Ahn, and G. Caire, “Joint spatial division and multiplexing?the largescale array regime,” IEEE transactions on information theory, vol. 59, no. 10, pp. 6441–6463, 2013.
 [16] A. Adhikary, E. Al Safadi, M. K. Samimi, R. Wang, G. Caire, T. S. Rappaport, and A. F. Molisch, “Joint spatial division and multiplexing for mmwave channels,” IEEE J. Sel. Areas Commun., vol. PP, no. 99, 2014.
 [17] A. Molisch, Wireless Communications. WileyIEEE Press, 2005.
 [18] T. Bai, R. Vaze, and R. W. Heath, “Analysis of blockage effects on urban cellular networks,” IEEE Transactions on Wireless Communications, vol. 13, no. 9, pp. 5070–5083, 2014.
 [19] A. B. Constantine et al., “Antenna theory: analysis and design,” MICROSTRIP ANTENNAS, third edition, John wiley & sons, 2005.
 [20] L. Liu, C. Oestges, J. Poutanen, K. Haneda, P. Vainikainen, F. Quitin, F. Tufvesson, and P. De Doncker, “The cost 2100 MIMO channel model,” IEEE Wireless Communications, vol. 19, no. 6, pp. 92–99, 2012.
Comments
There are no comments yet.