Next Article in Journal
Integrated Assessment of Coastal Subsidence in Nansha District, Guangzhou City, China: Insights from SBAS-InSAR Monitoring and Risk Evaluation
Next Article in Special Issue
Ship Target Detection in Optical Remote Sensing Images Based on E2YOLOX-VFL
Previous Article in Journal
A Comprehensive Assessment of the Pansharpening of the Nighttime Light Imagery of the Glimmer Imager of the Sustainable Development Science Satellite 1
Previous Article in Special Issue
A Lightweight Man-Overboard Detection and Tracking Model Using Aerial Images for Maritime Search and Rescue
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sensing and Deep CNN-Assisted Semi-Blind Detection for Multi-User Massive MIMO Communications

1
School of Software Engineering, Tongji University, Shanghai 201804, China
2
Engineering Research Center of Key Software Technologies for Smart City Perception and Planning, Ministry of Education, Shanghai 200003, China
3
School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China
4
Beijing Institute of Technology Chongqing Innovation Center, Chongqing 401120, China
5
School of Information and Communication Engineering, Beijing University of Posts and Telecommunications, Beijing 100876, China
6
School of Electronic and Information Engineering, Tongji University, Shanghai 201804, China
*
Author to whom correspondence should be addressed.
Submission received: 14 December 2023 / Revised: 1 January 2024 / Accepted: 5 January 2024 / Published: 8 January 2024

Abstract

:
Attaining precise target detection and channel measurements are critical for guiding beamforming optimization and data demodulation in massive multiple-input multiple-output (MIMO) communication systems with hybrid structures, which requires large pilot overhead as well as substantial computational complexity. With benefits from the powerful detection characteristics of MIMO radar, we aim for designing a novel sensing-assisted semi-blind detection scheme in this paper, where both the inherent low-rankness of signal matrix and the essential knowledge about geometric environments are fully exploited under a designated cooperative manner. Specifically, to efficiently recover the channel factorizations via the formulated low-rank matrix completion problem, a low-complexity iterative algorithm stemming from the alternating steepest descent (ASD) method is adopted to obtain the solutions in case of unknown noise statistics. Moreover, we take one step forward by employing the denoising convolutional neural network (DnCNN) to preprocess the received signals due to its favorable performance of handling Gaussian denoising. The overall paradigm of our proposed scheme consists of three stages, namely (1) target parameter sensing, (2) communication signal denoising and (3) semi-blind detection refinement. Simulation results show that significant estimation gains can be achieved by the proposed scheme with reduced training overhead in a variety of system settings.

Graphical Abstract

1. Introduction

With the vigorous development of the wireless communication industry and the deep integration with intelligent information processing technology, various remote sensing applications, e.g., vehicle-to-everything, smart manufacturing and environmental monitoring, come to the force, which puts forward increasing demand for network capacity and spectrum utilization [1,2,3]. As an important solution to support ultra-high transmission rate and intelligence requirements of the future mobile communications, massive multiple-input multiple-output (MIMO) inherently allows for simultaneous serving of multiple users by deploying large-scale antennas at a base station (BS) and facilitates the generation of directed beams towards target users, thereby expanding connected density and reducing the mutual interference in scenarios involving multiple users [4,5]. To obtain larger antenna gains and alleviate the pressure of spectrum scarcity, exploring novel transceiver design and frame structure operating at higher frequencies, such as the millimeter wave (mmWave) band [6,7], paves the way to guarantee reliable communications while enabling satisfactory user experiences.

1.1. Background

Considering that major obstacles preventing the deployment of massive MIMO are related to huge hardware costs and power consumptions imposed on the excessive number of radio frequency (RF) chains, hybrid analog/digital structures have been utilized to canonically route the RF chains via a well-designed phase-shifter network [8,9,10]. Furthermore, to compensate for prominent path loss encountered in mmWave channels and mitigate severe interference among multiple users, optimization on hybrid beamforming is expected to strike a favorable tradeoff between cost and rate [11,12,13]. One of the critical premises for obtaining the optimal beamforming is to acquire channel state information. Nevertheless, accurate channel estimation is challenging due to increased dimensionality of signal matrix and the application of hybrid architecture that limits the available observations. In this way, developing an efficient receiver scheme to enhance the transmission performance is essential in hybrid massive MIMO systems.
Recently, based on the fact that radar and communication systems share similar underline signal processing and gradually work in approximate frequency bands, there emerges a technological trend to integrate both functionalities into a single system [14,15,16], which offers an exciting opportunity to implement radar sensing by inherent wireless infrastructures. Given the potential ability of rapid and wide-range detection, it would be beneficial to leverage the radar sensing to fully extract geometric information from the surrounding environment [17,18] and thus assist in obtaining more precise parameter estimations. Toward this end, the coexistence of radar and communication components will constitute the primary architecture of this paper for enhanced performance gains.

1.2. Problem Being Addressed

Against the above background, as the dimensionality of antennas and carrier frequencies proliferate, the task of effectively balancing detection performance with computational cost becomes increasingly arduous. Consequently, the critical problem to be addressed is to propose novel solutions for comprehensive extraction of valuable information from the received signals in hybrid massive MIMO communication systems, thereby contributing to more reliable and cost-effective channel estimators.

1.3. Existing Solutions

In the literature, various channel estimators have been conducted for guiding the follow-up decoding and beamforming operations in hybrid massive MIMO systems. Among them, many techniques refer to the non-blind detections, which perform explicit channel estimation basically by transmitting and dealing with pilot symbols [19,20,21]. However, as the number of antennas increases, the pilot overhead could be enlarged accordingly and makes a negative impact on spectral efficiency. To address such issue, several research works began to focus on leveraging the sparse representation exhibited in delay-domain or angle-domain of mmWave channels [22,23,24], based on which a number of compressed-sensing (CS) tools have been adopted to figure out the formulated sparse recovery problem. Considering the grid-based CS estimator may lead to effects of basis mismatch, another gridless method appeared to construct continuous parameter domain by employing the atomic norm minimization [25]. Moreover, in order to avoid substantial computational complexity caused by high-dimensional optimizations and exploit more information from the payload data, our prior work [26] has proposed a semi-blind detection scheme for a hybrid massive MIMO system over frequency-selective fading channels, which utilized the uplink data as virtual pilots to improve the accuracy of channel estimation and data detection. However, the existing research has not sufficiently capitalized on the intrinsic features of massive MIMO systems and realistic channel characteristics, which constitute a significant challenge for achieving better tradeoff between training overhead and detection performance.
In addition to searching for rich scattering assumptions at higher frequencies, the low-rankness of massive MIMO systems has also been widely considered. In [27], according to the fact that both the number of BS antennas and the number of symbols in one coherence block are typically larger than the number of active users, the joint channel and data estimation was modeled as a low-rank matrix completion problem, for which two iterative algorithms were put forward to effectively recover the original channel and data. Moreover, schemes exploiting both the low-rank property and the knowledge of array responses have been proposed in [28,29,30], where the initial channel estimates provided by solutions to an inductive matrix completion would be further refined through multiple stages. Without loss of generality, most of the existing methods need to recover the matrix factorizations following an alternating minimization procedure [31,32,33]. Adhering to the existing methodologies, the intricate matrix multiplication or matrix inversion would induce unacceptable computational costs, especially when generalized to large-scale problems. As an effort to relieve the computational burden, an alternating steepest descent (ASD) algorithm was introduced in [34] to replace the least squares (LS) subproblem solutions with exact line-search updates. In view of the simplicity and moderate accuracy in the presence of non-negligible noise, we are inspired to tailor the ASD algorithm for estimating channels in hybrid massive MIMO systems.
Instead of simply obeying stochastic distributions, more realistic channels consisting of multiple propagation paths usually display a certain geometric structure. Since the process of communication channel estimation is similar to the radar target detection to some extent, research on making full use of sensing parameters to support communication functionality has attracted growing attention [35,36,37]. For example, by combining both functionalities of phased-MIMO radar and hybrid communication in the mmWave band, a novel strategy for joint target detection and channel estimation was proposed in [35]. To reduce the spectral resources as well as hardware cost, the employment of radar sensing for assisting the prediction of motion parameters in vehicular communication networks was investigated in [36]. Moreover, a robust MIMO-radar-aided channel estimation scheme deployed in multi-user (MU) MIMO communication systems has proven to be advantageous in improving the estimation accuracy with fewer training overhead [37]. Given that channel estimation may suffer from severe performance degradation in low signal-to-noise ratio (SNR) regions due to conceivable noise enhancement along with the LS-based estimation process, a supervised deep-learning (DL)-based signal denoiser was also involved in [37] to eliminate the noise on received signals before estimating the channel gains. However, training a neural network with high-dimensional signals and large datasets seems not tolerable in terms of complexity and latency [38]. Therefore, developing a cost-effective signal denoising module to make the state-of-the-art intelligent algorithms practically implementable for massive MIMO channel estimation is worth exploring.

1.4. Motivations and Objectives

A summary of the existing literature is provided in Table 1. The motivation for this paper stems from the notable limitations in current receiver techniques for hybrid massive MIMO communication systems. As one can see, challenges of large pilot overhead and substantial computational complexity will arise due to increased dimensionality and complicated signal nature in hybrid architectures, which restrict observation capabilities. Additionally, the growing demand for employing radar sensing in extracting environment information motivates us to seek advisable coordination between sensing and communication modules. In light of these challenges, the objective of this paper is to develop an efficient detection scheme tailored for hybrid massive MIMO systems with the aid of radar sensing, aimed at enhancing estimation performance and alleviating computational burdens.

1.5. Our Proposal and Contributions

In this paper, we propose a sensing and deep CNN-assisted semi-blind detection refinement scheme for a hybrid massive MU-MIMO system, where different functional modules, i.e., radar and communication, are integrated into the BS by partitioning multiple antennas. Specifically, both the low-rankness and geometric characteristic of high-frequency channels are fully exploited in a progressive manner. To reduce the large training overhead and substantial computational complexity, a novel transmission frame and an efficient receiver scheme are presented that can be split into three stages, namely (1) target parameter sensing, (2) communication signal denoising and (3) semi-blind detection refinement. Below, we crisply summarize the main contributions in this paper:
  • A novel time-division duplex (TDD) transmission frame capable of coordinating the radar and communication operations is designed, based on which the root multiple signal classification (MUSIC) algorithm is firstly applied to environment sensing and then the extracted target angle information is utilized for refining the subsequent communication detection results.
  • A generic representation for analog combining with phase shifters is considered, and the signal recovery problem is transformed into a low-rank matrix completion. To obtain the matrix factorizations with lower complexity, an iterative algorithm modified from ASD is proposed without any prior knowledge of noise statistics. In addition, different from the conventional pilot-only method, the semi-blind detection scheme is employed with reduced training overhead.
  • A pre-trained denoising convolutional neural network (DnCNN) is adopted to preprocess the received signals before performing the semi-blind detection, which attempts to handle Gaussian noise removal with unknown noise level and shows powerful ability of improved accuracy especially in low SNR regions.
The remainder of this paper is organized as follows. In Section 2, the system model is presented and the conventional pilot-only estimation method is reviewed. By transforming the signal recovery problem into a low-rank matrix completion, our proposed sensing and deep CNN-assisted semi-blind detection scheme is detailed with three stages in Section 3. Simulation results are provided in Section 4. Finally, conclusions are given in Section 5.
Notations: The list of acronyms used in this paper is summarized in Table 2. Unless otherwise specified, bold uppercase letters are used to represent matrices, bold lowercase letters are used to represent vectors and scalars are denoted by normal font. A ( i , : ) and A ( : , j ) denote the i-th row and j-th column of A , respectively. A i , j denotes the ( i , j ) -th entry of A . A T , A H and A * denote the operations of transpose, Hermitian transpose and element-wise conjugate, respectively. diag ( a ) stands for a diagonal matrix with diagonal elements given by the vector a . E and vec ( · ) represent the expectation operation and vectorization operation, respectively.

2. System Model and Background

In this paper, we consider a sensing-assisted hybrid massive MU-MIMO communication system, as shown in Figure 1, where K single-antenna users are served by a BS with two modules: one for the uplink communication and the other for the target sensing. Similarly as introduced in [37], the system operates in TDD mode and the multi-functional transmission frame structure is designed as shown in Figure 2, which consists of two parts. In the former part, the sensing module transmits probing signals and receives the echoes to detect the target directions. In the latter part, the communication module receives the pilots and payload data to complete the channel estimation and recover the uplink signals.
Suppose the uniform linear arrays (ULAs) are employed and are parallel with the users such that the angle-of-arrival (AoA) is identical to the angle-of-departure (AoD). Given the massive MIMO scenario and asymptotical orthogonality of steering vectors, we assume the sensing signals are emitted through line-of-sight (LoS) channels. Meanwhile, for the uplink communication, the transmitted signals from each single-antenna user will be scattered by the others around it when regarding them as point targets [36]. In addition, we consider the block-fading channel model, whose coefficients stay constant during the coherence block T.
Specifically, we first introduce the sensing model for assisting the AoA/AoDs estimation of target users. Then, the uplink communication model is presented, based on which the channel estimation problem for the hybrid massive MU-MIMO system is formulated. Finally, the conventional linear minimum mean squared error (LMMSE) estimator for solving the corresponding problem is reviewed.

2.1. Sensing Model

For the sensing module, we consider the mono-static radar case, where both the transmit and receiver arrays are equipped with N s antenna elements and we assume there is sufficient isolation between them. Let X s C N s × T s denote the probing signal sent by the BS, then the corresponding echoes reflected by target users can be expressed as [35]
Y s = k = 1 K β k b ( θ k ) b H ( θ k ) X s + N s ,
where N s CN ( 0 , σ s 2 I ) represents the additive white Gaussian noise plus interference, β k CN ( 0 , σ β k 2 ) denotes the reflection coefficient of the k-th target user, θ k denotes the k-th user’s azimuth angle relative to the BS and b ( θ k ) C N s × 1 represents the steering vector satisfying
b ( θ k ) = 1 N s [ 1 , e j 2 π λ d sin ( θ k ) , , e j 2 π λ d ( N s 1 ) sin ( θ k ) ] T ,
where λ and d represent the signal wavelength and inter-element spacing, respectively. Generally, we set d = λ / 2 .
For notational brevity, we arrange all steering vectors into a uniform matrix B ( Θ ) = [ b ( θ 1 ) , , b ( θ K ) ] C N s × K with Θ = { θ 1 , , θ K } and set Λ s = diag ( β 1 , , β K ) . Then, we rewrite the sensing signal model in (1) as a matrix form given by
Y s = B ( Θ ) Λ s B H ( Θ ) X s + N s .

2.2. Uplink Communication Model

For the uplink communication module, suppose the BS adopts a hybrid analog–digital architecture with N a antennas and N r   ( N r < N a ) RF chains to combine the incoming signals as depicted in Figure 1. Let x k [ t ] be the t-th transmitted symbol from user k with unit average power. Then, the received signals across N a BS antennas is given by
r c [ t ] = k = 1 K h k x k [ t ] + n c [ t ] , 1 t T c ,
where n c [ t ] CN ( 0 , σ c 2 I ) and h k C N a × 1 denotes the channel between user k and the BS.
Considering the hybrid arrays deployed, r c [ t ] will further pass through the RF phase shifter network. In particular, let W [ t ] C N r × N a denote the analog combining matrix, then the signal model combined in the RF band is consequently formulated as
y c [ t ] = W [ t ] H x c [ t ] + W [ t ] n c [ t ] ,
where y c [ t ] C N r × 1 indicates the received signal after passing through the analog structure, x c [ t ] C K × 1 denotes the transmitted signal vector from K users at time slot t and H [ h 1 , , h K ] C N a × K denotes the channel matrix between K users and the BS.
Next, we elaborate on the high-frequency channel model containing the intrinsic geometric structure. Instead of simply modeling the channels by stochastic distributions as [27], we adopt the scattering model here to describe the geometric environment over which the communication takes place. Suppose the users equipped with single-antenna are randomly distributed within communication distance, then the transmitted signal from each user will be scattered by its surrounding users. In this case, according to the widely used Saleh–Valenzuela channel model as presented in [7], the multi-path channel vector between user k and the BS can be provided as
h k = N a K l = 1 K α l , k a ( θ l , k ) ,
where α l , k and θ l , k denote the complex scattering coefficient and the azimuth angle relative to the BS for the l-th path of user k, respectively. Specifically, the path l = k contributes to the direct path from user k to the BS [36]. In addition, a ( θ l , k ) C N a × 1 denotes the steering vector expressed as
a ( θ l , k ) = 1 N a [ 1 , e j 2 π λ d sin ( θ l , k ) , , e j 2 π λ d ( N a 1 ) sin ( θ l , k ) ] T .
By referring to (5)–(6), we arrive at the input–output relationship for uplink communication signals, which is given by
y c [ t ] = W [ t ] A ( Θ ) G x c [ t ] + v [ t ] ,
where v [ t ] W [ t ] n c [ t ] and v [ t ] CN ( 0 , σ c 2 W [ t ] W [ t ] H ) , A ( Θ ) = [ a ( θ 1 ) , , a ( θ K ) ] C N a × K corresponds to the steering vector matrix and G C K × K merges the channel path gains and scattering coefficients that satisfies
g i , j = N a K α l , k , i = l , j = k , i , j K .
According to the communication signal model in (8), even with sufficient resolution for target sensing parameters, i.e., the geometric pattern of the channel matrix A ( Θ ) is known to the BS, the stochastic channel gains of different paths are still unable to infer. Therefore, a certain amount of pilots is required for completing the channel estimation.

2.3. The Conventional LMMSE Estimator

The conventional pilot-only estimation methods may cause large training overhead, especially with massive MIMO scenarios. To be specific, let X p [ x c [ 1 ] , , x c [ T p ] ] C K × T p denote the transmitted pilots, where T p indicates the number of pilots. Then, the signal model of (8) can be transformed as [27]
y c [ t ] = x p T [ t ] W [ t ] I K A ( Θ ) g + v [ t ] ,
where g vec ( G ) C K K × 1 and ⊗ denotes the Kronecker product.
In terms of whether the azimuth angles information is known or unknown to the BS, we define the channel counterparts to be estimated as g and h u I K A ( Θ ) g , respectively. By further defining y p [ y c T [ 1 ] , , y c T [ T p ] ] T C N r T p × 1 , v p [ v T [ 1 ] , , v T [ T p ] ] T C N r T p × 1 and W p [ x p T [ 1 ] W [ 1 ] ; ; x p T [ T p ] W [ T p ] ] C N r T p × N a K , the LMMSE estimates of the uplink channels can be given by [39]
h ^ u = ( W p H W p + σ c 2 I ) 1 W p H y p ,
and
g ^ = A I H ( Θ ) W p H W p A I ( Θ ) + σ c 2 I 1 A I H ( Θ ) W p H y p ,
for which A I ( Θ ) = I A ( Θ ) and we suppose W [ t ] W [ t ] H = I according to popular choices of analog combining matrix.
Upon determining the estimated channels, let X d [ x c [ T p + 1 ] , , x c [ T c ] ] C K × T d with T d = T c T p indicating the length of payload data, then the LMMSE-based data detection can be performed by [39]
x ^ d [ t ] = ( W [ t ] H ^ ) H W [ t ] H ^ + σ c 2 I 1 ( W [ t ] H ^ ) H y c [ t ] ,
where H ^ = vec 1 ( h ^ u ) or H ^ = vec 1 ( A I ( Θ ) g ^ ) depends on whether the azimuth angles of target users are already known to the BS.
From the above discussions, we can conclude that for the conventional pilot-only estimator, the number of transmitted pilots and the variance of noise have significant influences on the accuracy of channel estimation. To be specific, for the most common case without knowledge of geometric information [27], since there are involved N a × K unknowns in H ^ , the number of required pilots needs to scale with T p N a K / N r , which may occupy a large portion of transmission frame given the massive receive antennas. For the case where the channel estimation is conducted after resolving the AoA/AoDs [37], although the number of required pilots can be reduced to T p K K / N r , the accuracy of the LMMSE estimator will be greatly affected by non-negligible noise even with super-resolution of angle detections [40]. Therefore, to reduce the large training overhead while maintaining acceptable estimation accuracy, we will propose a sensing and deep denoising-assisted semi-blind detection scheme to improve the performance achieved by the existing methods.

3. Proposed Sensing and Deep CNN-Assisted Semi-Blind Detection Scheme

In this section, we introduce the proposed scheme to recover the signals formulated in Section 2. Based on the transmission frame structure designed in Figure 2, the proposed estimator will successively address the AoAs estimation, low-rank matrix completion and semi-blind detection refinement to exploit both the low-rank property and the geometric knowledge of channels in hybrid massive MIMO systems. In particular, by integrating the low-complexity spectrum estimation algorithm and high-efficiency matrix completion algorithm, MIMO radar’s rapid detection characteristics can be fully utilized to assist in obtaining more accurate channel estimation as well as taking up less training overhead.

3.1. Target Parameter Sensing

Considering the potential superiority of MIMO radar in achieving high-resolution angular estimations of target users, especially in identifying user positions, scattering paths and channel scenarios, the primary task of the sensing module is to detect and acquire the geometric directions of target users.
As shown in Figure 2, the angular estimation phase at the sensing module consists of two parts, during which the BS firstly sends probing signals and then receives the echo signals. Suppose the azimuth angles of K users are randomly distributed as Θ U [ π 2 , π 2 ] . To achieve the optimal estimation performance, we assume the probing signals sent by different antennas are mutually orthogonal such that the energy is evenly dispersed at each angle, i.e.,
R s = 1 T s X s X s H = P s N a I N s ,
where P s denotes the total energy of the BS. Moreover, under the assumption that the additive white Gaussian noise is unrelated to the probing signals and the averaged received power at each antenna being fixed as one, we can obtain the covariance matrix of the received echoes as  
R Y s = 1 T s Y s Y s H = B ( Θ ) R s ˜ B H ( Θ ) + σ s 2 I N a ,
where
R s ˜ = 1 T s Λ s B H ( Θ ) B ( Θ ) Λ s * .
To recover the AoAs of different target users, subspace decomposition-based methods can be invoked to separate the effective signals from background noise. Specifically, the signal subspace contains information from different directions, while the noise subspace comprises unrelated noise information. Then, following the standard subspace decomposition process [41], we present the signal and noise subspaces by taking eigenvalue decomposition of (15) as follows
R Y s = ( U s , U n ) Λ r ( U s , U n ) H ,
where Λ r C N s × N s is a diagonal matrix with eigenvalues in descending order and U s C N s × K and U n C N s × ( N s K ) correspond to eigenvector matrices that constitute the signal and noise subspaces, respectively. It can be seen that the columns of B ( Θ ) and U s span the same subspace. Without loss of generality, given the orthogonal relationship between the signal and noise subspaces, the pseudo-spectrum used for locating the AoAs of K target users can be readily formulated as
P a ( θ ) = 1 | b H ( θ ) U n U n H b ( θ ) | .
Instead of applying the classical approach that finds the angles at peaks through tedious spectrum search, we here employ the root-MUSIC approach, which transforms the problem of spectral peak detection into polynomial root finding, thus eliminating the dependence on a predefined angular grid. As a result, the proposed approach can avoid the grid mismatch phenomenon induced by superficial quantization and identify distinct source signals even when their angular differences are minimal, which allows for super-resolution angular estimation with low complexity. Specifically, according to [30], we rewrite the denominator of (18) as
b H ( θ ) U n U n H b ( θ ) = n = 0 N s 1 m = 0 N s 1 e j ξ m   sin θ C m , n e j ξ n   sin θ = l = N s + 1 N s 1 C l e j ξ l   sin θ ,
where ξ = 2 π d λ , C = U n U n H represents the correlation matrix with C l denoting the sum of elements along its l-th diagonal. Let z e j ξ   sin θ , then we can further simplify (19) to a polynomial
f ( z ) = l = N s + 1 N s 1 C l z l .
Then, the 2 ( N s 1 ) roots of f ( z ) exactly correspond to the poles of the MUSIC spectrum given by
z i = | z i | e j arg ( z i ) , i = 1 , , 2 ( N s 1 ) ,
K roots out of which lying closest to the unit circle are chosen to yield the AoAs as
θ ^ i = sin 1 ( arg ( z i ) ξ ) , i = 1 , , K ,
where arg ( z i ) denotes the phase angle of z i .

3.2. Low-Rank Matrix Completion

When the sensing module obtains the knowledge of all the AoAs of target users, the angle information will be transferred to the communication module, which is devoted to identifying the coefficients of different scattering paths by using the uplink pilot and payload data. For the sake of saving the training overhead while making full use of available information at hand, we consider to explore the low-rankness of the massive MIMO system and then recover the uplink channels based on low-rank matrix completion.
Before constructing the signal recovery problem, we show a generic model of the analog combining matrix and discuss how the transmitted signals from K users are related to the received signals at the BS. Specifically, given the fully-connected hybrid structure as presented in Figure 1 such that each RF chain is connected to all antennas, various realizations of the phase shifter with constant modulus can be modeled by ω e j ϕ and then make up the combining matrix W [ t ] . To this end, we can regard the generation of analog combining matrix at each time slot as randomly selecting N r rows from a DFT matrix F C N a × N a with elements given by [42]
F m , n = 1 N a e j 2 m n π N a , m , n = 0 , , N a 1 .
According to the signal model formulated in (5), by further defining H ˜ F H C N a × K and z [ t ] H ˜ x c [ t ] C N a × 1 , we can find that the received signal y c [ t ] is exactly composed of N r elements from z [ t ] . Then, let Ω C N r × T c represent the index set indicating each independent selection of W [ t ] during the communication period, and denote Y c Ω C N a × T c as the observation matrix with entries in Ω given by { y c [ t ] } t = 1 T c and otherwise zero, then the received signals can be recast as a matrix form
Y c Ω = P Ω ( Z + N ) ,
where Z [ z [ 1 ] , , z [ T c ] ] C N a × T c , N CN ( 0 , σ c 2 I ) and P Ω ( · ) represents a linear operation that preserves the entries in Ω while filling them not in Ω with zero.
Based on (24), we plausibly view Y c Ω as an incomplete observation of the matrix Z corrupted by the Gaussian noise N , where Z = H X c with X c [ x c [ 1 ] , , x c [ T c ] ] C K × T c . Noting the fact that in massive MIMO scenarios, the number of users is normally much smaller than the number of BS antennas and the length of coherence block, i.e., N a , T c K , the matrix Z displays an inherent low-rank property, i.e., rank [ Z ] K min { N a , T c } . Moreover, since multiplying by unitary matrix does not change the distributions, H ˜ and H are interchangeable for the channel estimation. Consequently, as derived in [27], the channel and data estimation is formulated as the following low-rank matrix completion problem
min U C N a × K , V C K × T c rank [ U V ] , s . t . Y c Ω = P Ω ( U V + N ) ,
where U and V are related to the solutions of channel and data estimates H ^ and X ^ c , respectively.
Different from the LMMSE estimator, which requires channel and noise second order statistics, we here propose to employ a low-complexity algorithm invoked in [34] to solve (25) without any prior information and exclusively based on matrix-wise updates. In this vein, considering the original rank-minimization problem is non-convex and generally NP-hard, an approximate relaxation scheme such as Frobenious norm minimization comes to the scene, which suggests to recover the matrix factorizations from
( U * , V * ) = arg min U C N a × K V C K × T c 1 2 Y c Ω P Ω ( U V ) F 2 .
To solve the above optimization problem, alternative minimization techniques, i.e., successively finding out U and V with the minimized residual, have been widely used due to their simplicity and flexibility. As an effort to circumvent the computational burden, a better choice called ASD algorithm incorporates a simple line-search to update solutions instead of dealing with the least squares subproblems. To be specific, let f V ( U ) and f U ( V ) define the gradient of the function f ( U , V ) = | | Y c Ω P Ω ( U V ) | | F 2 with respect to U and V , respectively, the stepsizes for updating solutions along the steepest gradient descent direction can be explicitly computed as [34]
t U = | | f V ( U ) | | F 2 | | P Ω ( f V ( U ) V ) | | F 2 ,
and
t V = | | f U ( V ) | | F 2 | | P Ω ( U f U ( V ) ) | | F 2 ,
where
f V ( U ) = ( Y c Ω P Ω ( U V ) ) V H ,
f U ( V ) = U H ( Y c Ω P Ω ( U V ) ) .
For initializations, the intuition behind (8) tells that the random realization of the channel mainly depends on the distribution of G assuming perfect estimates of angles. In view of this, given the angle information Θ ^ estimated by Section 3.1, we initialize U 0 = A ( Θ ^ ) G 0 C N a × K and V 0 C K × T c with elements of G 0 C K × K and V 0 drawing from i.i.d. Gaussian random variables with zero mean and unit variance. Then, the key idea of our proposed algorithm is to alternatively update U and V until the termination criteria are reached. Specifically, for the ( t + 1 ) -th iteration, after obtaining the steepest descent stepsizes and gradient descent directions following (27)–(30), we propose to update U t to U t + 1 in a two-step way, i.e.,
U t + 1 = U t t U t f V t ( U t ) , U t + 1 = A ( Θ ^ ) G t + 1 ,
where
G t + 1 = ( A H ( Θ ^ ) A ( Θ ^ ) ) 1 A H ( Θ ^ ) U t + 1 .
As a similar procedure, V t will be updated to V t + 1 by
V t + 1 = V t t V t f U t + 1 ( V t ) ,
where
t V t = | | f U t + 1 ( V t ) | | F 2 | | P Ω ( U t + 1 f U t + 1 ( V t ) ) | | F 2 .
Note that obtaining the gradient and stepsize will involve the product of a residual matrix and a N a × K or K × T c matrix. However, such operation only implements once at the beginning of each iteration and then can be efficiently used to update variables. Moreover, since the above approach replaces the least square solutions with an exact line-search step and thus avoids complex operations of matrix inversions, the computational cost for each iteration is much smaller than that of alternating least squares (ALS) and the conventional singular value decomposition (SVD) methods.

3.3. Refined Semi-Blind Detection

3.3.1. Ambiguity Removal

Although the above approach is able to provide solutions to (26) directly from the received observation Y c Ω , the non-uniqueness of matrix factorization may lead to great estimation bias. More clearly, suppose the initial channel and data estimates obtained from the ASD iterations are denoted by U ^ and V ^ , respectively. If there exists an invertible matrix Σ C K × K satisfying H X c = U ^ Σ Σ 1 V ^ , then solving (26) would not guarantee a unique recovery for H and X . Therefore, a certain amount of pilots are required to resolve the ambiguity.
On this basis, considering the number of elements that need to be determined is K × K , we develop a pilot and data placement as X c [ X p , X d ] , where X p C K × T p with T p K . Then, we have
V ^ = Σ X c = Σ [ X p , X d ] .
Due to the fact that pilot sequences are usually selected to be orthogonal, i.e., X p X p H = T p · I K , we can obtain Σ as
Σ = V ^ ( : , 1 T p ) · X p H T p ,
where V ^ ( : , i j ) refers to the sub-matrix consisting of i-th to j-th column vectors of V ^ . Upon resolving Σ , the data and channel estimation can be accordingly given by [26]
X ^ d = Σ 1 V ^ ( : , T p + 1 T c ) ,
and
H ^ = U ^ V ^ ( : , 1 T p ) X p H / T p .
The process of our proposed sensing-assisted semi-blind detection scheme is detailed in Algorithm 1. By carefully inspecting the workflow of the proposed approach, we can observe that the resulting estimation is not directly obtained from matrix factorization solutions; indispensable pilots are needed to facilitate the overall estimation as shown in Line 14, which explains the origin of the name “semi-blind detection”. The advantage of employing semi-blind detection lies in the fact that while the conventional non-blind detection implements channel estimation based on pilots only as displayed in (11)–(12), semi-blind detection can reduce the number of required pilots from an order of N a to K, naturally leading to the consideration of its application in large-scale MIMO systems.
Since the core idea of semi-blind detection hinges on full exploitation of prior information carried by both known pilots and estimated payload data, the efficiency of semi-blind detection mainly depends on the efficiency of the preceding matrix completion processes as illustrated in Line 4–13. As the scale of antennas or users increases, some computational delays may arise. Fortunately, instead of adopting the conventional LMMSE/ALS, the proposed ASD also offers advantages in terms of cost-effectiveness, which will be validated by experimental results in Section 4.
Algorithm 1 The proposed sensing-assisted semi-blind detection with reduced pilot overhead
1:Input  A ( Θ ^ ) C N a × K , Y c Ω C N a × T c , X p C K × T p
2:Randomly initialize G 0 C K × K and V 0 C K × T c .
3:Set U 0 = A ( Θ ^ ) G 0 and t = 0 .
4:repeat
5:Step (1) Update U t to U t + 1 as follows
6: Compute f V ( U ) = ( Y c Ω P Ω ( U V ) ) V H
 and t U = | | f V ( U ) | | F 2 / | | P Ω ( f V ( U ) V ) | | F 2 .
7: U t + 1 U t t U t f V t ( U t ) .
8: U t + 1 A ( Θ ^ ) ( A H ( Θ ^ ) A ( Θ ^ ) ) 1 A H ( Θ ^ ) U t + 1 .
9:Step (2) Update V t to V t + 1 as follows
10: Compute f U ( V ) = U H ( Y c Ω P Ω ( U V ) )
 and t V = | | f U ( V ) | | F 2 / | | P Ω ( U f U ( V ) ) | | F 2 .
11: V t + 1 V t t V t f U t + 1 ( V t ) .
12:Step (3) Update t to t + 1 .
13:until Convergence
14:Get Σ V ^ ( : , 1 T p ) · X p H / T p .
15:Output  H ^ = U ^ Σ

3.3.2. Signal Denoising

Given that the estimation accuracy of the semi-blind method partially relies on assumptions about the signal structure, it may suffer from pilot contamination and exhibit performance degradation, especially in the presence of high noise levels. In this way, it may be beneficial to firstly recover the latent clean signal from the noisy observations before performing the semi-blind detection for enhancing the reliability of payload data estimates.
With the recent development of deep learning, approaches based on CNNs have achieved remarkable success in the task of signal denoising, among which DnCNN formulated in [43] provides the state-of-the-art performance. Due to the integration of residual learning and batch normalization in the network architecture design, DnCNN has been validated to be highly effective in removing the Gaussian-type noise, which is exactly suitable for our considered signal models, especially at high noise levels. Specifically, the residual learning makes the inputs to each layer Gaussian-like distributed and less correlated with the image content, therefore avoid oversmoothing at high noise level. Furthermore, batch normalization alleviates the internal covariate shift by incorporating a normalization step and a re-scaling step before the nonlinear layer, which makes the model less sensitive to the noise level and generalizes well to inputs with strong noise. Moreover, DnCNN is shown to be able to extend to general denoising tasks, for example, the blind Gaussian denoising, which means the capacity of signal denoising without being informed of the noise level. Based on the above reasons, we propose to apply the DnCNN as a preprocessing procedure before performing the semi-blind detection, which aims to denoise the received signal Y c C N r × T c and expects better estimation accuracy in noisy environments.
Specifically, we choose to use the pre-trained model provided in [43] for Gaussian denoising, which takes the noisy signals as the input and generates the residual noise as the output. By subtracting the residual from the input, we then obtain the denoised signal Y c . It is worth noting that the model is without fine-tuning on our dataset because our signals exhibit similar distributions with the Gaussian noise, which can be seen by comparing Figure 3d,e, making it difficult for the DnCNN model to distinguish our signal from the noise. In this way, if the model is trained using our dataset, it tends to overestimate the noise level and over-smooth the input signal. On the contrary, the pre-trained model is trained with natural images with distinct distribution from the Gaussian noise, which can be seen by comparing Figure 3d,f, making it easier to implement proper denoising strength and avoid over-smoothing the input signal. Therefore, the pre-trained model is adopted, which not only saves the effort of the time-consuming training process but also enhances the denoising accuracy.
The overall framework of our proposed sensing and deep CNN-assisted semi-blind detection scheme for hybrid massive MU-MIMO systems is summarized in Figure 4, which can be divided into three stages, i.e., target parameter sensing, communication signal denoising and semi-blind detection refinement. In the sensing stage, the BS detects and acquires the azimuth angles of target user equipments (UEs) by sending probing signals and receiving echoes. Then, the uplink communication signals transmitted by UEs are processed in the RF band and denoised by using the pre-trained DnCNN model. Finally, based on the estimated angles and denoised signals, the refined semi-blind detection will be performed at the communication module.

4. Simulation Results

Simulation results are provided in this section to demonstrate the superior performance of the proposed sensing and deep CNN-assisted semi-blind detection scheme. We consider a hybrid massive MU-MIMO system as given in Figure 1, where K = 6 single-antenna users tend to communicate with the BS-integrated radar and communication modules. For the sensing module, the transmit and receive arrays of radar are equipped with N s = 64 , and the length of snapshots are set to be T s N s . For the communication module, the receiver is equipped with N a = 64 antennas and N r = 32 chains. Here, we assume the pilot sequences are mutually orthogonal with the pilot length satisfying T p K . In addition, we choose the payload data independently from the QPSK constellation with unit average power. Moreover, all the scattering and reflection coefficients are supposed to obey the standard complex Gaussian distributions.

4.1. The Proposed Transmission Frame

We firstly compare the channel estimation performance of adopting the multi-functional transmission frame as shown in Figure 2 with some existing schemes, i.e., the conventional pilot-only estimator [39] and the original semi-blind estimator relying solely on the communication module without aid of radar sensing [27]. Specifically, the channel estimation performance is evaluated by the normalized mean squared error (NMSE) defined as NMSE = E [ | | H ^ H | | F 2 ] / E [ | | H | | F 2 ] . For fair comparison, we assume the process of solving the low-rank matrix completion in Algorithm 1 is temporarily replaced by the same ALS iteration as derived in [27], and the noise variance is approximated according to SNR = 10 lg ( P ¯ / σ 2 ) with P ¯ being the average received power. The other configurations for the transmission frame are set as T s = 100 , T p = 2 K and T d = 100 . Unless otherwise stated, each simulation result is obtained through Monte-Carlo with 1000 independent channel realizations. Figure 5 depicts the NMSE performance versus SNR by using different transmission schemes. We can observe that the proposed sensing-assisted scheme shows significant performance gains especially with low SNRs compared to the alternative estimators, which thanks to that more geometric information about the scattering environments can be gleaned from the target sensing stage, while the other two schemes neglect the structural characteristics of the multi-path channel models.
To further validate the superiority of utilizing radar sensing function to yield the target angles rather than extracting the path angles from the roughly estimated channels as conducted in [30], which designed a three-stage estimator (denoted by “TSTE”) and suggested to obtain the knowledge of array responses directly from an initial low-rank channel estimate, Figure 6 illustrates the root mean squared error (RMSE) defined as RMSE = 10 lg E [ k = 1 K | | θ ^ k θ k | | 2 2 / K ] to measure the performance of AoAs resolutions of target users by using different detection schemes. Since the estimation deviation of the subspace decomposition-based method mainly depends on the dimension of received covariance matrix and some preset parameters such as the array spacing, we consider two settings of d / λ . It can be seen that the proposed sensing-assisted detection scheme achieves better resolutions, especially when the inter-element spacing is larger. Moreover, it is noteworthy that the RMSE under our proposed scheme would improve as the SNR increases, while that obtained by the TSTE estimator only fluctuates slightly versus different SNR, which is reasonable because their path angles are learned from the low-rank matrix solutions instead of straight from the received echoes, and the noise components have been removed before constructing the signal subspaces.

4.2. Refined Semi-Blind Detection

Next, we illustrate the performance gain of our proposed Algorithm 1. To show the advantage of performing semi-blind detection at the receiver rather than the LS/LMMSE-based detection [37], suppose the comparison experiment is set as performing the LMMSE channel estimation after obtaining the knowledge of target angles from the sensing module. Considering two settings of pilot length, T p = K and T p = 1.5 K , Figure 7 plots the NMSE performance of channel estimation with different receiver schemes. From Figure 7, we observe that the NMSE performance of Algorithm 1 outperforms the LMMSE method even if the sensing module is involved, which indicates the necessity of exploiting low-rank characteristics in massive MIMO systems to improve the estimation accuracy with relatively lower pilots. In addition, we can note that the NMSE performance of Algorithm 1 improves as the length of data sequences T d increases, while the performance of LMMSE remains almost the same in all settings. The reason is that by incorporating data into the decomposition process of received signals, Algorithm 1 can take full advantage of payload data as “virtual pilots” to glean more useful information over the conventional pilot-only method.
We then characterize the performances of Algorithm 1 in terms of the computational efficiency and the estimation accuracy without any prior knowledge of channel or noise statistics. By fixing T p = K and T d = 100 , Figure 8 depicts the convergence rates by employing different types of iterations, where “R-ALS iteration” refers to the regularized alternating least squares approach proposed by [27] for solving the formulated low-rank matrix completion problem (26). As shown in Figure 8, both the ASD iterations presented in Algorithm 1 and the R-ALS iterations will converge within limited iterations, and their convergence behaviors are similar for different SNR levels. Although the convergence rate of R-ALS iterations seems faster that that of ASD iterations, the calculations of U t + 1 and V t + 1 in Algorithm 1 avoid two matrix inversion operations for each iteration, which would otherwise involve additional complexity with cubic order of the matrix dimension, i.e., O ( N a 3 ) . Therefore, the total computational overhead incurred by Algorithm 1 can be relatively reduced especially with high-dimensional matrices.
To further validate the computational complexity of the proposed ASD approach, we conduct more experiments compared with the conventional LMMSE and R-ALS, focusing on their execution time and memory consumption. Following the principles of controlled variable experimental methodology, and with reference to the results illustrated in Figure 7 and Figure 8, we initially quantify the averaged computational cost for all three approaches, under the conditions of SNR = 2   dB and the number of iterations set to 6. Subsequently, given that the LMMSE approach would induce a single execution at the receiver adhered to Equations (12) and (13) without convergence issues, we only record the time and memory expenditures for the proposed ASD and R-ALS iterations by establishing the termination criterion as ϵ = | | H ^ t H ^ t 1 | | F 2 / | | H ^ t | | F 2 < 0.01 . The outcomes are presented in Table 3, from which we can observe that the proposed ASD approach exhibits lower time and memory consumption compared to its existing alternatives, thereby demonstrating the advantages in reduced complexity and higher computational efficiency.
In order to visualize the influence with unknown prior knowledge on the estimation accuracy, Figure 9 compares the NMSE performance versus different SNRs by adopting Algorithm 1 and the R-ALS-based estimator, which mimics and modifies the process in [27] given the target angles. For the former scenarios without available noise, the maximum iteration number is set to be t max = 20 ; while for the latter, the regularization parameters used to incorporate the prior information are set to be consistent with the noise variance. It can be seen that benefiting from taking noise effect into account, the R-ALS estimator performs better with low SNR values, which is owing to the non-negligible noise that imposes great challenges on accurate estimations in such SNR ranges. Nevertheless, the performance differences between two algorithms become narrow as the SNR increases, and Algorithm 1 would even achieve better performance under the condition of higher SNRs. By combining the results in Figure 8 and Table 3, we can conclude that our proposed Algorithm 1 is a good choice for refining semi-blind detection due to its comparable performance and lower complexity.

4.3. Signal Denoising by DnCNN

Finally, we confirm the effectiveness of signal denoising by DnCNN before performing semi-blind detection in low SNR regions. To evaluate the denoising performance of adopting different DnCNN models for our received signal dataset, we select the pre-trained DnCNN models in [43] for Gaussian denoising with specific noise levels (which are referred to as DnCNN-15 and DnCNN-25) and for blind Gaussian denoising (which is referred to as DnCNN-B) as comparisons. During the denoising phase, we successively obtain the output Y d e n o i s e as the denoised data of the input Y n o i s y = Y c . Suppose the clean signal that refers to Y n o i s e l e s s is related to the combined results of H X c , then the following performance metrics include NMSE and NMSE-Y defined as NMSE-Y = E [ | | Y d e n o i s e Y n o i s e l e s s | | F 2 ] / E [ | | Y n o i s e l e s s | | F 2 ] , which are used to measure the accuracy of channel estimation and signal denoising, respectively.
The averaged NMSE and NMSE-Y results over 1000 samples by adopting different signal denoising models are shown in Table 4, where the best results for each SNR value are highlighted in bold. As one can see, all the selected DnCNN models outperform the competing method without signal denoising when SNR is below 2 dB due to the consideration of noise influence in low SNR regions. In contrast, as the noise level decreases, both DnCNN-25 and DnCNN-B produce better results on NMSE-Y while fail to reach higher estimation accuracy, which is intuitively reasonable because the models with strong denoising power are likely to distort the original signal structures. It is pleasant to find that the model DnCNN-15, which is trained with lower noise level, can not only yield the best denoising performance on most of the discussed SNRs but also generate better channel estimates, which indicates the feasibility of deploying the signal denoising before performing the semi-blind detection.

5. Conclusions

In this paper, we have proposed a sensing and deep CNN-assisted semi-blind detection refinement scheme for uplink hybrid massive MU-MIMO systems. The main idea is to extract the geometric structure corresponding to the communication channels via radar sensing, based on which both the transmitted pilot and payload data are contributed to yield the refined estimates. We have treated the signal recovery problem as acquiring factorization solutions to a low-rank matrix completion, for which an iterative algorithm modified from ASD is presented without any prior knowledge of noise statistics. Specifically, a pre-trained DnCNN has been adopted for signal denoising to alleviate the noise enhancement along with the LS-based estimation process, especially in low SNR regions. The overall approach has been carried out in three progressive stages, i.e., target parameter sensing, communication signal denoising and semi-blind detection refinement. Consequently, simulation results have been provided to validate the superior performance of our proposed scheme. It is seen that the utilization of co-designed multi-functional transmission frame and the refined semi-blind receiver could achieve better estimation accuracy with reduced training overhead, which facilitates the availability of sensing-assisted communications.
Despite the advancements introduced by the proposed architecture, focusing on leveraging radar sensing capabilities to enhance communication gains, there still remain several limitations requiring further investigation. Instead of relying on the coexistence and information sharing between different modules, future work could develop more powerful signal processing techniques to integrate radar and communication modules into a unified hardware platform and decongest the RF environment by using the same signal for both functionalities. In such a case, devising sophisticated signal models that harmoniously incorporate the distinctions between LoS and multi-path channels might unveil new dimensions of system performance.
Another crucial aspect that needs to be considered in future work involves the time-varying property of propagation paths under user mobility scenarios, including time-delay parameters and Doppler shift effects. This aspect necessitates the BS to pre-compensate for these effects before demodulating the received signals for most robust detection performance. Toward this end, we plan to explore the time-varying channel tracking algorithms and compensation methods to mitigate the adverse effects of temporal fluctuations and phase deviations on estimation errors, which should be of great importance to cater to real-time applications.

Author Contributions

Conceptualization, F.H. and J.Z.; methodology, F.H. and L.Z.; software, H.Z. and J.W.; validation, J.Z. and J.W.; formal analysis, H.Z. and L.Z.; investigation, F.H.; resources, L.Z. and H.Z.; data curation, J.Z. and J.W.; writing—original draft preparation, F.H.; writing—review and editing, L.Z. and H.Z.; visualization, J.W.; supervision, J.Z.; project administration, F.H.; funding acquisition, F.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the National Natural Science Foundation of China under Grants 62201388 and 62201389, in part by the Natural Science Foundation of Shanghai under Grant 22ZR1463400 and in part by the Fundamental Research Funds for the Central Universities under Grant 2022-5-YB-01.

Data Availability Statement

The data are not publicly available due to privacy.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Tong, W.; Li, G.Y. Nine challenges in artificial intelligence and wireless communications for 6G. IEEE Wirel. Commun. 2022, 29, 140–145. [Google Scholar] [CrossRef]
  2. Li, X.; Liu, F.; Zhou, Z.; Zhu, G.; Wang, S.; Huang, K.; Gong, Y. Integrated sensing, communication, and computation over-the-Air: MIMO beamforming design. IEEE Trans. Wirel. Commun. 2023, 22, 5383–5398. [Google Scholar] [CrossRef]
  3. Koubaa, A.; Ammar, A.; Abdelkader, M.; Alhabashi, Y.; Ghouti, L. AERO: AI-enabled remote sensing observation with onboard edge computing in UAVs. Remote Sens. 2023, 15, 1873. [Google Scholar] [CrossRef]
  4. Zhu, Y.; Guo, H.; Lau, V.K.N. Bayesian channel estimation in multi-user massive MIMO with extremely large antenna array. IEEE Trans. Signal Process. 2021, 69, 5463–5478. [Google Scholar] [CrossRef]
  5. Guo, J.; Wen, C.-K.; Jin, S.; Li, G.Y. Overview of deep learning-based CSI feedback in massive MIMO systems. IEEE Trans. Commun. 2022, 70, 8017–8045. [Google Scholar] [CrossRef]
  6. Zhang, R.; Cheng, L.; Wang, S.; Lou, Y.; Wu, W.; Ng, D.W.K. Tensor decomposition-based channel estimation for hybrid mmWave massive MIMO in high-mobility scenarios. IEEE Trans. Commun. 2022, 70, 6325–6340. [Google Scholar] [CrossRef]
  7. Wei, X.; Hu, C.; Dai, L. Deep learning for beamspace channel estimation in millimeter-wave massive MIMO systems. IEEE Trans. Commun. 2021, 69, 182–193. [Google Scholar] [CrossRef]
  8. Nguyen, N.T.; Ma, M.; Lavi, O.; Shlezinger, N.; Eldar, Y.C.; Swindlehurst, A.L.; Juntti, M. Deep unfolding hybrid beamforming designs for THz massive MIMO systems. IEEE Trans. Signal Process. 2023, 71, 3788–3804. [Google Scholar] [CrossRef]
  9. Mirzaei, J.; ShahbazPanahi, S.; Sohrabi, F.; Adve, R. Hybrid analog and digital beamforming design for channel estimation in correlated massive MIMO systems. IEEE Trans. Signal Process. 2021, 69, 5784–5800. [Google Scholar] [CrossRef]
  10. Wei, X.; Dai, L. Channel estimation for extremely large-scale massive MIMO: Far-field, near-field, or hybrid-field? IEEE Commun. Lett. 2022, 26, 177–181. [Google Scholar] [CrossRef]
  11. Chafaa, I.; Negrel, R.; Belmega, E.V.; Debbah, M. Self-supervised deep learning for mmWave beam steering exploiting sub-6 GHz channels. IEEE Trans. Wirel. Commun. 2022, 21, 8803–8816. [Google Scholar] [CrossRef]
  12. Wu, K.; Zhang, J.A.; Huang, X.; Guo, Y.J.; Hanzo, L. Simultaneous beam and user selection for the beamspace mmWave/THz massive MIMO downlink. IEEE Trans. Commun. 2023, 71, 1785–1797. [Google Scholar] [CrossRef]
  13. Cheng, Z.; Liao, B. QoS-aware hybrid beamforming and DOA estimation in multi-carrier dual-function radar-communication systems. IEEE J. Sel. Areas Commun. 2022, 40, 1890–1905. [Google Scholar] [CrossRef]
  14. Liu, F.; Cui, Y.; Masouros, C.; Xu, J.; Han, T.X.; Eldar, Y.C.; Buzzi, S. Integrated sensing and communications: Toward dual-functional wireless networks for 6G and beyond. IEEE J. Sel. Areas Commun. 2022, 40, 1728–1767. [Google Scholar] [CrossRef]
  15. Cui, Y.; Liu, F.; Jing, X.; Mu, J. Integrating sensing and communications for ubiquitous IoT: Applications, trends, and challenges. IEEE Netw. 2021, 35, 158–167. [Google Scholar] [CrossRef]
  16. Wang, J.; Varshney, N.; Gentile, C.; Blandino, S.; Chuang, J.; Golmie, N. Integrated Sensing and Communication: Enabling Techniques, Applications, Tools and Data Sets, Standardization, and Future Directions. IEEE Internet Things J. 2022, 9, 23416–23440. [Google Scholar] [CrossRef]
  17. Graff, A.; Chen, Y.; González-Prelcic, N.; Shimizu, T. Deep learning-based link configuration for radar-aided multiuser mmWave vehicle-to-infrastructure communication. IEEE Trans. Veh. Technol. 2023, 72, 7454–7468. [Google Scholar] [CrossRef]
  18. Wang, X.; Zhai, W.; Zhang, X.; Wang, X.; Amin, M.G. Enhanced automotive sensing assisted by joint communication and cognitive sparse MIMO radar. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 4782–4799. [Google Scholar] [CrossRef]
  19. Hassan, K.; Masarra, M.; Zwingelstein, M.; Dayoub, I. Channel estimation techniques for millimeter-wave communication systems: Achievements and challenges. IEEE Open J. Commun. Soc. 2020, 1, 1336–1363. [Google Scholar] [CrossRef]
  20. Zhang, R.; Yang, L.; Tang, M.; Tan, W.; Zhao, J. Channel estimation for mmWave massive MIMO systems with mixed-ADC architecture. IEEE Open J. Commun. Soc. 2023, 4, 606–613. [Google Scholar] [CrossRef]
  21. Wang, Y.; Zhang, Y.; Tian, Z.; Leus, G.; Zhang, G. Super-resolution channel estimation for arbitrary arrays in hybrid millimeter-wave massive MIMO systems. IEEE J. Sel. Top. Signal Process. 2019, 13, 947–960. [Google Scholar] [CrossRef]
  22. Fan, D.; Gao, F.; Liu, Y.; Deng, Y.; Wang, G.; Zhong, Z.; Nallanathan, A. Angle domain channel estimation in hybrid millimeter wave massive MIMO systems. IEEE Trans. Wirel. Commun. 2018, 17, 8165–8179. [Google Scholar] [CrossRef]
  23. Zhou, L.; Dai, J.; Xu, W.; Chang, C. Sparse channel estimation for intelligent reflecting surface assisted massive MIMO systems. IEEE Trans. Green Commun. Netw. 2022, 6, 208–220. [Google Scholar] [CrossRef]
  24. Wang, Y.; Qi, C.; Li, P.; Lu, Z.; Lu, P. Channel estimation for wideband mmWave MIMO OFDM system exploiting block sparsity. IEEE Commun. Lett. 2022, 26, 897–901. [Google Scholar] [CrossRef]
  25. Chu, H.; Zheng, L.; Wang, X. Semi-blind millimeter-wave channel estimation using atomic norm minimization. IEEE Commun. Lett. 2018, 22, 2535–2538. [Google Scholar] [CrossRef]
  26. Han, F.; Wang, X.; Deng, H. A very-low pilot scheme for mmWave hybrid massive MIMO-OFDM systems. IEEE Wirel. Commun. Lett. 2021, 10, 2061–2064. [Google Scholar] [CrossRef]
  27. Liang, S.; Wang, X.; Ping, L. Semi-blind detection in hybrid massive MIMO systems via low-rank matrix completion. IEEE Trans. Wirel. Commun. 2019, 18, 5242–5254. [Google Scholar] [CrossRef]
  28. Li, X.; Fang, J.; Li, H.; Wang, P. Millimeter wave channel estimation via exploiting joint sparse and low-rank structures. IEEE Trans. Wirel. Commun. 2018, 17, 1123–1133. [Google Scholar] [CrossRef]
  29. Vlachos, E.; Alexandropoulos, G.C.; Thompson, J. Massive MIMO channel estimation for millimeter wave systems via matrix completion. IEEE Signal Process. Lett. 2018, 25, 1675–1679. [Google Scholar] [CrossRef]
  30. Masood, K.F.; Hu, R.; Tong, J.; Xi, J.; Guo, Q.; Yu, Y. A low-complexity three-stage estimator for low-rank mmWave channels. IEEE Trans. Veh. Technol. 2021, 70, 5920–5931. [Google Scholar] [CrossRef]
  31. Nguyen, L.T.; Kim, J.; Shim, B. Low-rank matrix completion: A contemporary survey. IEEE Access 2019, 7, 94215–94237. [Google Scholar] [CrossRef]
  32. Giampouras, P.V.; Rontogiannis, A.A.; Koutroumbas, K.D. Alternating iteratively reweighted least squares minimization for low-rank matrix factorization. IEEE Trans. Signal Process. 2019, 67, 490–503. [Google Scholar] [CrossRef]
  33. Chi, Y.; Lu, Y.M.; Chen, Y. Nonconvex optimization meets low-rank matrix factorization: An overview. IEEE Trans. Signal Process. 2019, 67, 5239–5269. [Google Scholar] [CrossRef]
  34. Tanner, J.; Wei, K. Low rank matrix completion by alternating steepest descent methods. Appl. Comput. Harmon. Anal. 2016, 40, 417–429. [Google Scholar] [CrossRef]
  35. Liu, F.; Masouros, C.; Petropulu, A.P.; Griffiths, H.; Hanzo, L. Joint radar and communication design: Applications, state-of-the-art, and the road ahead. IEEE Trans. Commun. 2020, 68, 3834–3862. [Google Scholar] [CrossRef]
  36. Yuan, W.; Wei, Z.; Li, S.; Yuan, J.; Ng, D.W.K. Integrated sensing and communication-assisted orthogonal time frequency space transmission for vehicular networks. IEEE J. Sel. Top. Signal Process. 2021, 15, 1515–1528. [Google Scholar] [CrossRef]
  37. Huang, S.; Zhang, M.; Gao, Y.; Feng, Z. MIMO radar aided mmWave time-varying channel estimation in MU-MIMO V2X communications. IEEE Trans. Wirel. Commun. 2021, 20, 7581–7594. [Google Scholar] [CrossRef]
  38. He, H.; Wen, C.-K.; Jin, S.; Li, G.Y. Deep learning-based channel estimation for beamspace mmWave massive MIMO systems. IEEE Wirel. Commun. Lett. 2018, 7, 852–855. [Google Scholar] [CrossRef]
  39. Kay, S.M. Fundamentals of Statistical Signal Processing, Volume III: Practical Algorithm Development; Prentice Hall: Upper Saddle River, NJ, USA, 2013. [Google Scholar]
  40. Gizzini, A.K.; Chafii, M.; Nimr, A.; Fettweis, G. Enhancing least square channel estimation using deep learning. In Proceedings of the 2020 IEEE 91st Vehicular Technology Conference: VTC2020-Spring, Antwerp, Belgium, 25–28 May 2020; pp. 1–5. [Google Scholar]
  41. Schmidt, R. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef]
  42. Liu, A.; Lau, V. Phase only RF precoding for massive MIMO systems with limited RF chains. IEEE Trans. Signal Process. 2014, 62, 4505–4515. [Google Scholar] [CrossRef]
  43. Zhang, K.; Zuo, W.; Chen, Y.; Meng, D.; Zhang, L. Beyond a gaussian denoiser: Residual learning of deep CNN for image denoising. IEEE Trans. Image Process. 2017, 26, 3142–3155. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Sensing-assisted hybrid massive MU-MIMO communication system. The blue solid lines denote the probing signal sent by the sensing module. The red solid lines denote the corresponding echoes reflected by target users. The green dashed lines denote the uplink communication signals.
Figure 1. Sensing-assisted hybrid massive MU-MIMO communication system. The blue solid lines denote the probing signal sent by the sensing module. The red solid lines denote the corresponding echoes reflected by target users. The green dashed lines denote the uplink communication signals.
Remotesensing 16 00247 g001
Figure 2. The structure of multi-functional transmission frame. T s denotes the duration of transmitted waveform by radar module. T p and T d denote the duration of uplink pilots and payload data, respectively.
Figure 2. The structure of multi-functional transmission frame. T s denotes the duration of transmitted waveform by radar module. T p and T d denote the duration of uplink pilots and payload data, respectively.
Remotesensing 16 00247 g002
Figure 3. Comparisons among sample signals of (a) Gaussian noise, (b) our signal and (c) natural image, with corresponding histograms in (d), (e), (f), respectively. (e) exhibits similar distribution with (d), while (f) is distinct from (d).
Figure 3. Comparisons among sample signals of (a) Gaussian noise, (b) our signal and (c) natural image, with corresponding histograms in (d), (e), (f), respectively. (e) exhibits similar distribution with (d), while (f) is distinct from (d).
Remotesensing 16 00247 g003
Figure 4. The overall framework of the proposed sensing and deep CNN-assisted semi-blind detection.
Figure 4. The overall framework of the proposed sensing and deep CNN-assisted semi-blind detection.
Remotesensing 16 00247 g004
Figure 5. The NMSE performance versus SNR by using different transmission schemes.
Figure 5. The NMSE performance versus SNR by using different transmission schemes.
Remotesensing 16 00247 g005
Figure 6. The detection performance of AoAs resolutions under different schemes.
Figure 6. The detection performance of AoAs resolutions under different schemes.
Remotesensing 16 00247 g006
Figure 7. The NMSE performance versus length of data sequences under different receiver schemes with SNR = 2   dB .
Figure 7. The NMSE performance versus length of data sequences under different receiver schemes with SNR = 2   dB .
Remotesensing 16 00247 g007
Figure 8. Convergence rates of adopting different iterative algorithms, where NMSE = E [ | | Z ^ Z | | F 2 ] / E [ | | Z | | F 2 ] with Z ^ = H ^ X ^ c .
Figure 8. Convergence rates of adopting different iterative algorithms, where NMSE = E [ | | Z ^ Z | | F 2 ] / E [ | | Z | | F 2 ] with Z ^ = H ^ X ^ c .
Remotesensing 16 00247 g008
Figure 9. The NMSE performance versus SNR by adopting different iterative algorithms.
Figure 9. The NMSE performance versus SNR by adopting different iterative algorithms.
Remotesensing 16 00247 g009
Table 1. Existing works comparison on receiver techniques.
Table 1. Existing works comparison on receiver techniques.
Existing WorksSystem ModelsReceiver TechniquesLimitations or Features
[19,20,21]mmWave MIMO communicationsnon-blind detection based on pilots onlyPilot overhead can be enlarged greatly with massive antennas.
[22,23,24]mmWave MIMO communicationsCS leveraging the channel sparsityGird-based estimator may lead to effects of basis mismatch.
[25]mmWave MIMO communicationsatomic norm minimizationHigh-dimensional optimizations may cause substantial computational complexity.
[26]mmWave MIMO communicationssemi-blind detection aided with payload dataNeglect of geometric structure of channel distributions.
[27]massive MIMO communicationsLS leveraging the channel low-ranknessUnderexploitation of realistic channel characteristic with multiple propagation paths.
[28,29,30]mmWave MIMO communicationsmultiple stages exploiting both the low-rankness and sparsityHigh computational complexity when generalized to large-scale problems.
[35]joint radar and communicationsmultiple signal classification (MUSIC) for angle estimationJoint signal processing strategy can simultaneously detect targets while estimating the communication channel.
[36]joint radar and communicationsmaximum likelihood (ML) estimatorSensing parameter estimation can promote the dynamic topology construction of surrounding environments.
[37]joint radar and communicationsLS detection vs. DL-based denoiserEfficient estimation with fewer training resources by eliminating noise before recovering channels.
Table 2. List of acronyms.
Table 2. List of acronyms.
AbbreviationsExpansion
ALSAlternating Least Squares
AoAAngle-of-Arrival
AoDAngle-of-Departure
BSBase Station
CNNConvolutional Neural Network
CSCompressed Sensing
DnCNNDenoising Convolutional Neural Network
DLDeep Learning
LMMSELinear Minimum Mean Squared Error
LoSLine-of-Sight
LSLeast Squares
mmWaveMillimeter Wave
MIMOMultiple-Input Multiple-Output
MU-MIMOMulti-User MIMO
MUSICMultiple Signal Classification
NMSENormalized Mean Squared Error
RFRadio Frequency
RMSERoot Mean Squared Error
SNRSignal-to-Noise Ratio
SVDSingular Value Decomposition
TDDTime-Division Duplex
UEUser Equipment
ULAUniform Linear Array
Table 3. Comparisons on computational cost by using different approaches.
Table 3. Comparisons on computational cost by using different approaches.
ApproachesAveraged Computational CostTotal Cost until Convergence
performancememorytimememorytime
proposed ASD0.216 Kb3.11 ms10.22 Kb22.80 ms
R-ALS8.644 Kb43.22 ms46.61 Kb86.96 ms
LMMSE0.636 Kb3.50 ms//
Table 4. Performance comparison by adopting different signal denoising models.
Table 4. Performance comparison by adopting different signal denoising models.
SNR Value−5 dB−4 dB−3 dB
Selected ModelNMSE-YNMSENMSE-YNMSENMSE-YNMSE
Without DnCNN3.17271.65612.54301.29671.99620.9956
DnCNN-151.52231.03621.13520.81590.83200.6679
DnCNN-250.85660.88750.86180.89510.88880.9199
DnCNN-B0.94240.96350.92650.94880.88880.9042
SNR Value−2 dB−1 dB0 dB
Selected ModelNMSE-YNMSENMSE-YNMSENMSE-YNMSE
Without DnCNN1.60250.87331.24250.74770.99510.6051
DnCNN-150.63350.62420.49980.59470.43120.5720
DnCNN-250.90700.93630.93140.95520.93560.9610
DnCNN-B0.93500.86310.82180.88750.77000.8672
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Han, F.; Zeng, J.; Zheng, L.; Zhang, H.; Wang, J. Sensing and Deep CNN-Assisted Semi-Blind Detection for Multi-User Massive MIMO Communications. Remote Sens. 2024, 16, 247. https://0-doi-org.brum.beds.ac.uk/10.3390/rs16020247

AMA Style

Han F, Zeng J, Zheng L, Zhang H, Wang J. Sensing and Deep CNN-Assisted Semi-Blind Detection for Multi-User Massive MIMO Communications. Remote Sensing. 2024; 16(2):247. https://0-doi-org.brum.beds.ac.uk/10.3390/rs16020247

Chicago/Turabian Style

Han, Fengxia, Jin Zeng, Le Zheng, Hongming Zhang, and Jianhui Wang. 2024. "Sensing and Deep CNN-Assisted Semi-Blind Detection for Multi-User Massive MIMO Communications" Remote Sensing 16, no. 2: 247. https://0-doi-org.brum.beds.ac.uk/10.3390/rs16020247

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop