Advanced Noncoherent Detection in Massive MIMO Systems via Digital Beamspace Preprocessing

: Noncoherent detection in massive multiple-input/multiple-output (MIMO) uplink systems provides a low-complexity alternative to its coherent counterpart. Requiring no actual channel knowledge but the per-user induced power at the base station, comparable performance as channel-estimation-based detection can be achieved when the users are located in the near-ﬁeld of the base station. However, noncoherent detection fails in scenarios where users are in the far-ﬁeld due to an insufﬁcient capability to separate the users in terms of their spatially induced power. For this purpose, a dielectric lens or an analog beamforming structure can be employed, which are capable to focus the power of the incident waves to a smaller subset of the antennas at the base station. These so-called analog beamspace techniques have been demonstrated to enable again the noncoherent detection scheme. Analogous to a spatial Fourier transform, beamspace techniques can be also realized in the digital domain offering more ﬂexibility. Its applicability to noncoherent detection is studied in this paper. It is shown numerically that by means of digital beamspace preprocessing, considerable performance gains can be achieved. Applied in dominant line-of-sight channels, a large number of users can be accommodated and the residual performance gap to coherent detection with perfect channel knowledge is minimal.


Introduction
The massive multiple-input/multiple-output (MIMO) technology embodies one of the key elements in future wireless communication systems [1][2][3][4][5][6]. Equipped with a very large number of antennas at the base station (BS), massive MIMO systems offer the potential to meet the ever-increasing demand for high throughput and energy efficiency. To this end, sufficiently accurate channel state information (CSI) is essential, which is gathered via pilot signaling using predefined (orthogonal) pilot symbols. Since the number of channel coefficients to be estimated is huge in massive MIMO systems, CSI acquisition can become a challenging task and may suffer from pilot contamination due the re-use of pilot symbols in neighboring communication cells [7]. Moreover, outdated CSI or estimation inaccuracies have a severe impact on the system performance [8,9]. Additionally, the channel-training overhead reduces the overall time, where data can be transmitted resulting in a loss of spectral efficiency [10].
In order to overcome the need of channel estimation and pilot symbol allocation, noncoherent detection schemes can be applied in the uplink as proposed in [11,12]. Here, actual channel knowledge is not required to recover the individual transmitted data streams but the per-user induced power over the receive antennas of the BS, namely the user-specific power-space profile (PSP) of the users. As long as the users can be separated at the BS in terms of their spatially induced power, noncoherent detection performs comparably to its CSI-based (coherent) counterpart [13,14]. However, the spatial separation can only be obtained in scenarios where the users are located in the Telecom 2020, 1 212 near-field of a physically-large BS [15][16][17][18]. Consequently an adaption is needed to handle far-field users and to keep the overall size of the BS in an acceptable range for a practical implementation. To this end, (analog) beamspace techniques can be employed focusing the energy of the incoming waves to a smaller subset of the BS antenna array depending on the angle of arrival. This process can also be interpreted as a spatial Fourier transform. Exploiting the angular sparsity of the propagation channel, where the uplink transmit signals of the users impinge at the BS from sufficiently different angular directions, the users can be separated enabling again the noncoherent detection scheme. This has been successively demonstrated for a lens-embedded and an analog beamforming-based system [19,20]. So far, lens antenna arrays were investigated solely in the context of coherent massive MIMO transceivers aiming at a reduction of hardware costs in both uplink and downlink [21][22][23][24][25][26]. The same applies to the studies in [27][28][29][30][31], where an alternative representation of the lens is used on the concept of beamspace MIMO.
To achieve the crucial user separation for noncoherent detection in far-field scenarios, the proposed analog beamspace techniques in [19,20] require additional bulky hardware, which can suffer from non-idealities, high complexity, and lacking adaptability to specific channel conditions. Thus, an all digital implementation is of major interest avoiding the stated issues and providing the potential of even improved user separability via adapted signal processing at the same time.
This work studies beamspace techniques in the digital domain and its application to the noncoherent detection scheme. To this end, the relevant requirements are identified and detailed to conceive a highly advanced and competitive noncoherent detection scheme compared to the classical and initial approach proposed in [11,12] by means of digital beamspace preprocessing. Originally, a high-quality phase reference and a perfect synchronization of the local oscillators (LOs) is not required in case of noncoherent detection. Consequently it can benefit from a low-complexity LO distribution network, where independent LOs for each receive branch can be employed for instance [17]. This applies also for the analog implementation of the beamspace techniques as in [19,20]. However, the application of the digital beamspace preprocessing requires a coherent demodulation of the receive branches to preserve the phase progression between adjacent receive channels. This can become a very challenging task, especially in massive MIMO systems, where the number of receive paths is quite large. Therefore, an extension to a sub-array architecture is explored in the context of digital beamspace preprocessing as well offering the opportunity to relax the design requirements of the LO distribution network. In summary, the contribution of this work is to combine the individual (known) concepts of noncoherent detection scheme, sub-array approach, and digital beamspace processing for a joint examination under proper channel conditions. To the authors' best knowledge, this is presented for the first time in the area of massive MIMO transmission. This paper is organized as follows. The system model is introduced in Section 2, which includes a generalized system overview when digital beamspace preprocessing is applied. After a brief recapitulation of the underlying channel model and the noncoherent detection scheme denoting the most important aspects, a detailed description of the digital beamspace preprocessing tailored to noncoherent detection is presented. Numerical results are shown in Section 3, where full-and sub-array architectures and the influence of the channel properties are discussed. Finally, Section 4 concludes the paper.
Notations: Bold-faced capital letters and bold-faced small letters denote matrices and column vectors, respectively. Superscripts (·) T , (·) H , and (·) −1 indicate the transpose, the conjugate transpose, and the inverse of a matrix, respectively. The operators E{·}, ||·|| 2 , and ||·|| F denote the expectation, the Euclidean norm of a vector, and the Frobenius norm of a matrix, respectively. I N is the N × N identity matrix and diag (·) represents a diagonal (block) matrix.

System Overview
In this paper, a multi-user uplink system is considered as depicted in Figure 1a. Here, N u single-antenna users simultaneously transmit to a central BS, where a uniform linear array is employed with a very large number of receive antennas N rx N u and an antenna distance of d a . The users are assumed to be located in the far-field, i.e., the (radial) user distance d u meets the far-field condition with respect to the physical antenna array size D = N rx d a at the BS and the wavelength λ according to d u > 2D 2 /λ [32]. Consequently, plane waves impinge at the receiver from various angular directions φ (measured with respect to the normal of the array axis) depending on the user positioning and characteristics of the propagation channel. In order to recover the individual data streams transmitted by the users, noncoherent detection schemes are applied at the BS. This requires the per-user received power induced at the antenna array of the BS to be spatially separated, which is not the case due to the far-field assumption. To this end, beamspace preprocessing is utilized in the digital domain after coherent downconversion and analog-to-digital conversion. Digitally and simultaneously generated beams pointing into different angular directions within a predefined field of view (FoV) focus the energy of the incoming waves to a smaller subset of the overall receive paths (see Figure 1a).   In the following, a joint digital beamspace preprocessing of all N rx receive branches is denoted as the full-array architecture. This system can exploit its full potential in terms of angular resolution and power-focusing capability. However, to form beams into a desired direction, a fully coherent demodulation of the received branches is required. The complexity of such a high-quality LO distribution network scales with the number of receive branches. In order to relax the design requirements, the full-array architecture can be divided into N sub sub-arrays of equal size (see Figure 1b), where coherency of the receive channels can be provided more easily and where one LO with proper distribution network can be employed instead. Consisting of N rx,sub = N rx /N sub antennas or rather receive branches, each sub-array can perform digital beamspace preprocessing individually, i.e., forming beams in the digital domain within a predefined FoV. Noteworthy, due to the reduced aperture size D = D/N sub , the (radial) user distance could be less to meet the far-field condition in the sub-array case.
Further details regarding digital beamspace preprocessing can be found in Section 2.4. First, the underlying channel model and the noncoherent detection scheme are briefly reviewed.

Channel Model
The massive MIMO propagation channel is modeled using a geometry-based stochastic approach, more precisely using an adapted version of the COST 2100 channel model [33,34]. Here, randomly placed clusters in the topological simulation environment emulate physical scattering objects, where each cluster consists of groups of multi-path components (MPC). Its parameters, such as complex amplitude and angular direction, are computed on the basis of the simulation area. Generally, the clusters can be assigned to the following four categories: (i) local clusters (scatterers surrounding the users and/or the BS), (ii) single clusters (scattering on one object), (iii) twin clusters (scattering on multiple objects), and (iv) common clusters (propagation through the same set of clusters). In addition to multi-path propagation, a line-of-sight (LOS) between the users and the BS can be considered as well. In order to reduce the complexity of the channel and to ease the interpretation of the subsequent numerical results, the propagation via LOS and local clusters around the users are taken into account.
Considering omni-directional single-antenna users located in the far-field and a BS employing a uniform linear array of directional antennas with pattern C rx (φ), the complex channel vector T of the u-th user in case of a flat-fading channel is given by where a LOS u is the complex amplitude of the LOS component and φ LOS u its corresponding direction of arrival (DOA). The p-th MPC originating from N p local clusters is denoted by a MPC u,p with DOA φ MPC u,p . a N rx (φ) represents the steering vector of a uniform linear array with N rx elements of distance d a , i.e., where φ denotes the DOA, k = 2π/λ is the wavenumber, and λ is the wavelength.
Since the local clusters are uniformly distributed in the azimuthal plane around the users, a MPC u,p and its DOA φ MPC u,p becomes random for different channel realizations, where φ MPC u,p can be generally characterized by an angular spread with standard deviation σ s . On the contrary, φ LOS u is deterministic and the power of the LOS component is adjusted by a LOS-to-MPC ratio (LMR), which is defined as

Noncoherent Detection
The received signal r k at the N rx BS antennas at time instant k can be written as where T comprises the transmit symbols of each user u, and n k = [n 1,k , . . . , n N rx ,k ] T gathers the circular-symmetric complex Gaussian noise n m,k with zero mean and variance σ 2 n at the m-th receive branch. The transmit symbols are differentially encoded unit-magnitude phase-shift keying (DPSK) symbols, which are generated as follows: where a k,u are the data symbols drawn from an M-ary PSK constellation M = {e j2πi/M | i = 0, 1, . . . , M−1}. Due to the PSK signaling, the total received power per user is ||h u || 2 2 and the signal-to-noise ratio (SNR) is given as ||h u || 2 2 /σ 2 n . Assuming the channel to be constant during a transmission burst (blockfading channel), noncoherent detection of user u is performed block-wise over a temporal block R of N bl time steps given by and on the basis of the N bl × N bl correlation matrix using a user-specific diagonal weighting matrix There are two detection strategies based on (7). Each user can be detected individually using decision-feedback differential detection (DFDD) [11], where decisions are made successively and previous detected symbols are considered. To cope with multi-user interference, DFDD can be combined with noncoherent decision-feedback equalization (DFDD/nDFE) [12], where the average interference of already detected users is subtracted additionally.
Both presented schemes take advantage of the fact that the per-user induced power known as power-space profile (PSP), exhibits a high concentration at some spatial location at the BS. Thus, the users can be separated for detection and W u in (8) acts as a filter extracting only the relevant entries in R of each user in the spatial domain. The windowing can be numerically optimized (via a global search) to obtain the maximum signal-to-interference-plus-noise ratio (SINR) in (7), when the PSPs in (9) are known [12,35]. However, due to users located in the far-field, P m,u is equally distributed over the receive antennas. Consequently, this demands an alternative solution to redistribute the power, which is presented in the following.

Digital Beamspace Preprocessing
In order to obtain an advantageous power distribution for noncoherent detection, digital beamspace preprocessing is applied, which is commonly used for DOA estimation in the literature [36][37][38][39][40]. To this end, the receiving block R in (6) is modified according tô where is a N rx × N b beamforming (BF) matrix. N b determines the number of simultaneous generated beams within a predefined FoV and has to be equal to the number of receive paths or rather BS antennas N rx . The received signals are superimposed using the respective weight vectors w m for the m-th receive branch and the resultant spatial (total) power spectrum can be expressed aŝ where Telecom 2020, 1 216 denotes the covariance matrix of the initial received signal r k . The PSPP m,u (spatial per-user power) after digital beamspace preprocessing can be defined aŝ and the resulting SINR of the u-th user at the m-th receive branch is given as The beamforming matrix W BF should be designed to maximize the SINR in (15) for the respective users and to fulfill the condition to ensure the noise to be spatially white. In the following, the composition of W BF is discussed for a full-and sub-array architecture.

Full-Array
Considering a full-array architecture, a first choice of the beamforming matrix is to select N rx steering vectors defined in (2) as the respective columns, i.e., and where φ m ∈ [−Φ/2, Φ/2] is the corresponding beam direction of the m-th receive branch chosen from a given FoV Φ. This method is also known as conventional or Bartlett beamforming (BBF) [36]. It can also be interpreted as an N rx -point discrete Fourier transform (DFT), where the weights w B,m represent a bank of narrowband filters each tuned to the spatial frequency which is related to the respective angles φ m or rather the spatial location m. A uniform sampling in the frequency domain given a FoV Φ leads to a spatial frequency spacing Considering the entire visible range Φ = π and the general antenna distance d a = λ/2, the N rx spatial beams become orthogonal with spacing ∆θ = 1/N rx . This corresponds to an angular spacing of ∆φ = sin −1 (λ/D), which is a measure of the angular resolution of a uniform linear array with aperture size D = N rx d a also known as Rayleigh resolution limit [39]. Strictly speaking, this does not hold for reception at endfire directions, where the resolution drops due to a reduced array aperture size. Furthermore, the FoV is limited by the antenna pattern in practice. Thus, the orthogonality of the beams vanishes when uniformly chosen from a smaller angular range Φ < π. Alternatively, the antenna distance d a could be optimized (increased) for a specific FoV such that ∆θ = 1/N rx holds in (20) in order to restore the orthogonality.
The resultant spatial power distribution of a plane-wave phase front associated with direction φ can be expressed as [39] w H B,m a N rx (φ) which is approximately a squared sinc-function sampled at the chosen beam directions φ m of the m-th receive branch. Its peak is shifted towards the spatial location, where the direction of the incoming wave matches one of the generated beams, i.e., where φ m = φ. This is exemplary depicted in Figure 2, where two plane waves separated by two adjacent beams impinge at an array of N rx = 12 elements considering an FoV of Φ = π/2. Apparently, the signals are not entirely resolvable, since the selected beams are not orthogonal in case of d a = λ/2. This in turn results into (strong) spatial interference. Hence, when the angular separation of the impinging signals is small, the spatial interference rejection capability using Bartlett beamforming is limited. The optimal solution in terms of spatial filtering and maximizing the SINR is obtained via Capon's method (or a minimum variance distortionless response filter (MVDR)) [36,40,41], where the signal of interest is passed undistorted while the power from all other spatial frequencies and directions is minimized. Based on the knowledge of the covariance matrix in (13), this data-dependent approach leads to a higher resolution and thus power focusing capability than the conventional beamforming. However, mismatches in the array response due to multi-path propagation and imprecise knowledge of the covariance matrix, which has to be estimated in practice, lead to a severe performance degradation [40]. The same applies to subspace-based methods such as multiple signal classification (MUSIC) [36,42] or rotational invariance techniques (ESPRIT) [43], which offer high resolution capabilities.
Alternatively, a combination of Bartlett and eigenspace-based processing [44], denoted as Bartlett-eigenspace beamforming (BEBF) in the following, can provide a more robust approach against aforementioned issues. The corresponding overall weighting matrix can be expressed as where denotes the eigenspace-based beamforming matrix. In an initial step, conventional beamforming is performed for spatial prefiltering according tō Due to the beamspace transformation, the array response vector changes tō The covariance matrix of (24) can be computed via its sample estimate and an eigenvalue decomposition is applied as follows:Qr where the diagonal matrices Λ s and Λ n contain the N s largest (signal-subspace) and the N rx − N s smallest (noise-subspace) eigenvalues, respectively. V s and V n denote the corresponding eigenvectors. The respective weight vectors in (23) for eigenspace-based post-processing are composed according to Comparing the resultant spatial power distribution to conventional beamforming as illustrated in Figure 2, the additional eigenspace-based post-processing restores the orthogonality of the selected beams and achieves a higher resolution as expected. This provides more flexibility since the beams do not need to be selected from a regular grid anymore. Most importantly, the energy of the incoming wave is focused almost entirely to one receive branch reducing the spatial interference. Noteworthy, this approach is more computational complex than conventional beamforming, which can be implemented via a fast Fourier transform (FFT) later on. Moreover, the knowledge of the number of sources N s is required, i.e., the number of users N u , and this method can suffer from instabilities when performing the inverse of the estimated covariance matrix. Diagonal loading can counteract this issue [40].

Sub-Array
All above mentioned beamforming schemes for a full-array architecture can be applied analogously to the sub-array architecture. Instead of N rx beams according to the number of receive antennas/branches, the number of beams is limited to the per-sub-array elements N rx,sub = N rx /N sub and, when applied, the eigenspace-based post-processing is performed separately at each sub-array. The beam directions can be chosen from a different set and/or angular sector since each sub-array can process the incoming signals individually. However, the condition in (16) might be not fulfilled in this case. Therefore, the same set of beams is used per sub-array selected from a predefined FoV Φ. The corresponding weight matrix of the l-th sub-array is W BF sub ,l = [w 1 , w 2 , . . . , w N rx,sub ] (28) and the total beamforming matrix is constructed via

Numerical Results
Numerical simulations are conducted for an N u = 3 user scenario and a BS employing a uniform linear array of N rx = 128 directional patch antennas [15] with inter-element distance d a = λ/2. The FoV of the BS, where the incoming signals can be processed from, is assumed to be Φ = π/2. The users equipped with single omni-directional antennas are located under a given angle φ u in the far-field of the BS as illustrated in Figure 3. In the following, users 1 and 3 are positioned at a fixed angle of A propagation channel is considered, where both LOS and multi-path propagation are taken into account and its respective power ratio is set according to (3). The multi-path components of each user, which are in total N p = 3, impinge at the BS having an angular spread with standard deviation σ s from the nominal angle φ u . The respective channel matrices from the cluster-based channel model are normalized in such a way that ||H|| 2 F = N u . This power control ensures that the total induced sum power at the BS is equal to the number of users and the difference in channel attenuation between the users remain. If the total receive power of the users is equal, the total SNR per user reads 1/σ 2 n , which is used for setting the noise power σ 2 n at each antenna branch. The performance measure is the symbol error rate (SER), which is evaluated for noncoherent detection using DFDD/nDFE with (and without) digital beamspace preprocessing. To this end, transmit symbols from a quaternary DPSK constellation and a transmission burst of length N bl = 201 are considered. The SER is averaged over N ch = 50,000 channel realizations acquired from the adapted COST 2100 channel model.
The system configuration and the parameterization of both propagation channel and noncoherent detection are summarized in Table 1.

Full-Array vs. Sub-Array Architecture
First, the full-and sub-array architectures are investigated. To this end, a dominant LOS channel is considered, where the power of the LOS component is significantly larger than the sum power of the multi-path components originating from local clusters around the users. The LMR is set to 10 dB and the spatial spread of the multi-path components is σ s = 2 • . The PSPs of the users are assumed to be perfectly known and the beamforming weights are chosen optimally, i.e., the beams are evenly distributed in the given FoV Φ = π/2, whereby one of the generated beams is perfectly aligned to a respective user. For comparison, the SER performance of coherent detection via Bell Laboratories layered space-time (BLAST) [45,46] with perfect channel knowledge is included in the following considerations.
Applying the digital beamspace preprocessing, the induced power of the users at the BS, which has been uniformly distributed over the receive antennas due the to far-field location, is spatially redistributed and focused onto a smaller subset of the N rx receive branches.
First, regarding the full-array architecture, the resulting PSPs show a very narrow peak as depicted in Figure 4a, where user 2 is located at boresight (φ 2 = 0 • ). In case of conventional beamforming, the PSPs follow a squared sinc function, which is minimally broadened due to the low-power multi-path components. A combination with eigenspace-based beamforming yields even narrower power peaks, where most of the user power is concentrated at one receive branch. A difference in the total power of the users can be observed, which is due to a loss in received power for non-boresight users since directional patch antennas are employed at the BS.
The SER performance as function of the SNR is illustrated in Figure 4b when user 2 is again at the boresight location. While noncoherent detection fails without beamspace preprocessing due to uniform PSPs, the spatial redistribution of the user power in the digital domain via beamforming allows almost identical performance compared to the coherent detection with perfect channel knowledge. Here, no performance gains are observed when eigenspace-based processing is applied in addition. The worse performance of the users 1 and 3 is again related to the power loss using directional antennas.
Analyzing the SER as function of the angular position φ 2 of user 2 (see Figure 4c), the performance of the noncoherent detection scheme deteriorates when approaching one of the other users due to the overlapping of the respective PSPs. Here, additional eigenspace-based beamforming exhibits a significant improvement in terms of user separability since the digital redistributed user power is almost located at one receive branch. Noteworthy, coherent detection using BLAST exhibits a minor degradation in performance when the users are exactly at the same angular position.  Next, the sub-array architecture is considered allowing a low-complexity LO distribution network since each sub-array can employ its own LO. To this end, the antenna array of the BS is splitted into N sub = 8 sub-arrays each comprising N rx,sub = 16 antennas. Consequently, 16 different spatial beams can be generated, which are assumed to be the same for each sub-array. This leads to a periodic pattern of the PSPs as depicted in Figure 5a. The SER as function of the SNR remains almost unchanged when compared to the full-array architecture using Bartlett beamforming (see Figure 5b). A saturation is observed at high SNR, where the residual inter-user interference due to non-orthogonal beams begins to dominate, deteriorating the performance. A slight performance degradation is observed for Bartlett-eigenspace beamforming. Analyzing the SER vs. the angle φ 2 , the user separability is diminished due to a loss in angular resolution (see Figure 5c). Thus, closely spaced users cannot be separated anymore as compared to the full-array case since the individual PSPs of the users already start to overlap at larger angular distance. Here, eigenspace-based processing reveals again a superior performance than just conventional beamforming.  Results of different sub-array configurations are illustrated in Figure 6 considering various angular positions of user 2. The smaller the number of antennas per sub-array N rx,sub , the smaller the angular resolution and thus the user separability. At small N rx,sub , the conventional beamforming suffers from performance variations, which originates from strong inter-user interference caused by non-orthogonal beams. On the contrary, eigenspace post-processing offers a better user separation at all times.

Influence of Propagation Channel
In the following, the influence of the propagation channel is analyzed. To this purpose, the dominance of the LOS is varied by means of the LMR and the spatial spread σ s of the multi-path components is modified. The analysis is restricted to the full-array architecture (observations apply to the sub-array in the same manner). The stronger the LOS component compared to the multi-path components, the narrower are the individual PSPs enhancing the angular user separation for both conventional beamforming and the combination with eigenspace-based processing due to lower inter-user interference (see Figure 7). At the same time, higher channel correlations deteriorate the overall performance. Increasing the spatial spread σ s , the multi-path components become more relevant and cause a broadening of the PSPs in case of conventional beamforming. Consequently, a higher inter-user interference occur and the user separability is worse in these scenarios (see Figure 8). On the contrary, additional eigenspace-based processing is insensitive to higher spatial spreads regarding the user separation and exhibits a performance improvement due to lower channel correlations. In summary, noncoherent detection with beamspace preprocessing performs best when channel correlations are low, i.e., multi-path propagation is prevailing and the multi-path components impinge at a large angular spread at the BS. However, in order to achieve best performance in terms of user separability, a dominant LOS channel needs to be present, where the standard deviation of the multi-path components from the nominal LOS angle is preferably small at least for conventional beamforming. This allows to accommodate a large number of users. Ideally, one generated beam or rather receive branch can serve one user. Thus, the upper bound is N rx in the full-array and N rx,sub in the sub-array case. However, this limit might be reached in pure LOS channels or/and for eigenspace-based beamforming, where the orthogonality of the beams is given as compared to conventional beamforming.

Conclusions
In this paper, the application of digital beamspace preprocessing to a massive MIMO system is assessed, where noncoherent detection is performed in the uplink. First, the fundamentals and the requirements are discussed when both concepts are treated jointly and when a full-array or a sub-array architecture is considered. Then, the symbol error rate performance is analyzed for a three-user scenario on the basis of a cluster-based channel model capturing its principal behavior under physically accurate conditions. It is demonstrated that digital beamspace preprocessing transforms the unsuitable uniform PSPs, when the users are located in the far-field of the BS, into PSPs with a narrow peak located at the receive branch associated with the corresponding angle of arrival of the incoming wave. Hence, the users can be spatially separated when the signals impinge at the BS from sufficiently different angular directions enabling again the noncoherent detection scheme. An improved performance in terms of angular user separability is achieved when the conventional processing is extended by means of eigenspace-based beamforming. At the costs of higher computational complexity, this procedure reduces the inter-user interference in case of closely spaced users. Regarding a sub-array configuration, the PSPs exhibit a periodic pattern and the angular user separation is less due the reduced aperture size. Furthermore, a small performance degradation is observed.
The characteristics of the propagation channel can have a large impact on the system performance. Best results in terms of user separability are obtained in dominant LOS channels, which allows an accommodation of a large number of users. On the contrary, better overall SER performance is achieved in multi-path scenarios since the channel correlations are lower. However, a higher level of inter-user interference can occur for larger spatial spreads of the multi-path components reducing the total number of users, which can be served optimally.
Generally, the knowledge of the PSPs for noncoherent detection are unknown in practice. However, the estimation is straightforward when the angular position of the users is known beforehand. Otherwise, the common known concept of DOA estimation can be employed extracting the respective angular directions via a peak search in the angular spectrum. This works particularly in dominant LOS channels. The extracted knowledge of the angular position can also be used for the downlink transmission afterwards.
Comparing the results to coherent detection with perfect channel knowledge, the performance gap is minimal. However, the estimation overhead and potential channel estimation errors are not taken into account, which severely degrades its performance. When considered, the presented combination of the noncoherent detection scheme and digital beamspace preprocessing outperforms the CSI-based detection.