Abstract
For massive multipleinput multipleoutput (MIMO) systems, the realtime channel state information (CSI) acquisition is crucial but difficult in fast timevarying scenarios, especially for downlink (DL) channels in frequency division duplex (FDD) systems. This paper proposes an adaptive trainingfeedback scheme to estimate the CSI of DL channels. Specifically, first, base station (BS) determines a subspace containing uplink (UL) channel and an orthogonal normal basis of this subspace by using received signals in UL channel and the principle component analysis (PCA) technique. Second, by using the spatial reciprocity, it is found that the obtained subspace above also contains the DL channel vector. Thus, regarding the basis as pilots, the BS transmits pilots to mobile user (MU), and the user estimates received signals, feeds back them to the BS. Finally, according to information of the feedback, the BS can construct the DL channel vector. In this scheme, the pilots are adaptive to the change of the UL channel. Furthermore, the times of training is the dimension of the subspace, and feedback overhead is the coefficients of linear combination of DL channel vector under the basis only. Thus, cost of training and feedback can be greatly reduced. The simulation results show that the performance of the proposed scheme can approach the optimal scheme with very few training times and feedback overhead at high speed.
1. Introduction
Massive multipleinput multipleoutput (MIMO) technology occupies a pivotal position in wireless communication systems, owing to its provision of excellent spectrum and energy efficiency to systems. It also brings challenges, one of which is the acquisition of accurate channel state information (CSI). The problem is severe especially for frequency division duplex (FDD) systems in mobile scenarios, since the CSI obtained through the traditional pilotbased training method is outdated immediately. The estimated CSI in one time slot is not equal to the CSI of the following time slot due to rapid time variant of the channel.
Widely efforts have been made by scholars to investigate this problem. For example, [1, 2] both performed channel prediction and channel tracking under the autoregressive (AR) model, of which [1] proposed a linear finite impulse response (FIR) Wiener predictor. [2] used Kalman filter and the minimum meansquare error (MMSE) decision feedback equalizer to solve the tracking and equalization problem of the timevarying channel and can achieve sound channel tracking performance, but the complexity is higher than traditional adaptive algorithms. [3, 4] analyzed the overall rate performance of the FIR Wiener predictor in the case of delayed CSI.
Papers [5–7] also used Kalman filters to propose their own schemes. In [5], based on an assumption that the uplink (UL) downlink (DL) conversion problem is a linear estimation problem, an ULDL transform algorithm by using Kalman filtering was proposed. [6] proposed a spatialtemporal basis expansion model based on the characteristics of the largescale antenna arrays. This model can reduce the channel dimension. Furthermore, the reciprocity of the physical angle between the UL and DL channels was used to reduce the complexity of DL channel tracking. However, the timevarying spatial information tracking of UL is built based on the Kalman filter and Taylor series expansion, which leads to estimation accuracy deterioration.
Under the assumption of channel sparsity in massive MIMO systems, [7] proposed a sparse Bayesian learning framework based on the expectation maximization. [8] proposed a structured compressed sampling matching tracking algorithm, which can obtain reliable CSI. In [9], a channel prediction algorithm was proposed based on the firstorder Taylor expansion channel model, and the interval of effective prediction was derived. [10] proposed a dynamic turbo orthogonal approximate message passing algorithm based on a twodimensional Markov model, which can recursively track dynamic channels. However, the computational complexity is high in all these mentioned methods.
In this paper, based on the principal component analysis (PCA) technology and the spatial reciprocity of adjacent frequency bands in the FDD system, we propose an adaptive trainingfeedback scheme to estimate and track CSI. The main steps of our scheme are as follows: the first, the base station (BS) utilizes the PCA technology to construct a covariance matrix from the received signals, and the singular value decomposition (SVD) is performed on this matrix. The second, the left singular vectors corresponding to the first largest singular values are selected to form an orthonormal basis of the subspace of UL channel. According to the spatial reciprocity, this subspace is same as the subspace generated by the DL channel. Thus, an orthogonal basis to linearly represent the DL channel is obtained. The third, elements of this basis are transmitted to the mobile station (MS) as pilots, and the MS can estimate the received signals. These estimations can be regarded as coefficients of linear representation of DL channel under the orthonormal basis. The fourth, the MS finds a precoder in predesigned precodebook to match the coefficients, and sends the index of this precoder to the BS. The last, the BS utilizes the index, the precodebook and the orthonormal basis to construct a precoder and transmits data signals. It is worth mentioning that the pilots designed in this paper can adapt to the channel, and the proposed PCAbased adaptive trainingfeedback scheme has few training times, and the required feedback overhead is also very small.
The paper is organized as follows: in the next section, the 3D timevarying channel model, the spatial reciprocity in the FDD system and the PCA principle are briefly described. In Section 3, a PCAbased adaptive trainingfeedback scheme is proposed, and the specific steps of this scheme are described. The simulation results are shown in Section 4 and Section 5 summarizes the paper.
The following notations are given: , , and denote scalar, vector and matrix, respectively. , denote transpose and conjugate transpose, respectively. denotes the element in the matrix and denotes its rank. and are the absolute value and normed space, respectively. “” denotes the tensor product of matrices.
2. Preliminary
In this section, we will describe the system model, spatial reciprocity in FDD systems, and the principles of PCA.
2.1. System Model
We consider a massive MIMO system with a singleantenna user in FDD transmission mode, where the frequency of the UL channel is , the frequency of the DL channel is , and the frequency interval between them is . Assume that the antennas provided by the BS are arranged by an antenna array of , where and are the number of columns and rows of the BS antenna array and .
The channel model involved in this paper is given by the standard timevarying channel model proposed by the 3rd Generation Partnership Project (3GPP) Organization, see the document [11] for details. Assume that there are clusters per channel and propagation paths per cluster, we only consider the main propagation path of each cluster; that is, . Then the UL channel from the MS antenna to BS antennas at time slot , denoted as , can be expressed as where
is a steering matrix, and are the UL channel gains of different clusters. In this paper, we assume that they are independently and identically distributed (i.i.d.) complex Gaussian random variables with zeromean and unit variance. In Eq. (1) where is the wavelength of UL carrier frequency, and is the speed of MS, and are the travel elevations angle and azimuth angle of MS, respectively. The definition of the elevation angle in this paper is the direction of the signal propagation with respect to the zaxis in the global coordinate systems. The definition of is as follows where and are the elevation arrival angle (ZOA) and the azimuth arrival angle (AOA) in th cluster, respectively.
Correspondingly, the DL channel at time from BS to MS is obtained by where and are the DL channel gains of different clusters. Matrix is a steering matrix. Also, assume that they are i.i.d. complex Gaussian random variables with zeromean and unit variables. In Eq. (5) where is the wavelength of DL carrier frequency, and the transmit spherical unit vector is defined as where and are the elevation departure angle (ZOD) and the azimuth departure angle (AOD) in th cluster, respectively.
In Eq. (2), the steering matrix is composed of steering vectors , , where and are the elevation angles (ZOD, ZOA) and the azimuth angles (AOD, AOA) of the cluster, respectively. , where and are the MS antenna 3D position coordinates and the BS antennas 3D position coordinates. Since MS has only one antenna, where is a complex number. Matrix given in (6) is a steering matrix, and its components are then similar with the above. The details can be found in the following subsection.
2.2. The Spatial Reciprocity in FDD Systems
In FDD systems, since the reflections and propagation paths experienced by signals transmitted between the UL and DL channels are the same, the angles (AOA, ZOA, AOD, and ZOD) and delay of the UL and DL channels are all the same, which has been confirmed in [12]. The phenomenon is also called spatial reciprocity or channel reciprocity. This phenomenon is utilized to estimate the DL channel at different conditions, see [13–15]. Recently, specific details on the spatial reciprocity can be found in [16].
For the timevarying channel in this paper, the initial positions of the BS antennas and the MS antenna are determined, and the subsequent changes in antenna positions are mainly reflected in speed and delay, which means that and do not change with time.
According to the above analysis and the definition of the steering matrix, the difference in steering matrices of the UL and DL channels are only caused by the frequency interval between the UL and DL channels. In most cases, the frequency of the UL and DL channels is very high (such as GHz), and the frequency interval between the UL and DL is very small compared to the carrier frequency, then, according to the calculations given in [16], the difference between the UL and DL steering matrices caused by the frequency interval can basically be ignored. Consequently, we can make the assumption that the following holds
3. PCABased Adaptive TrainingFeedback Scheme
From equations (1), (5), and assumption we can see that the UL channel vector and the DL channel vector lying on the same subspace, the subspace generated by the columns of matrix . Thus, the precoder used in the DL transmissions should be in this subspace, and hence, it is worth getting this subspace.
To determine this subspace at the BS, it is useful to take the advantage of PCA based on the UL transmissions. After that, an orthonormal base of this subspace is transmitted on the DL channel as pilots, and from these received signals, the MS can evaluate the coefficients of combination of on this base. Using the designed codebook, information of these coefficients can be sent back to the BS, and the BS makes the precoder and data transmissions.
Based on the above, we give specific steps of the PCAbased adaptive trainingfeedback scheme, as shown in following algorithm.

For the proposed algorithm, we give further explanation as follows:
Step 1. Design a dimensional precoding codebook containing codewords. There are a lot of methods to design the codebook in the existing literatures, for example, see [17].
Step 2. The MS sends a signal to BS, and the received signal at the th time slot of BS is where is the signaltonoise ratio (SNR), is the UL channel at time slot , is a Gaussian standard noise.
Step 3. The BS combines the received signals at different times with a time interval in the period as , and takes its row mean .
Step 4. Calculate the covariance matrix of the received signals and normalize as
Step 5. Implement SVD on the covariance matrix as follows: where , , satisfying with , , , and .
Step 6. Take the left singular vectors corresponding to the first largest singular values of as pilots and send them to the MS through the DL channel. The received signals are where is the DL channel at time slot , , denote the noises. The Minimum Mean Square Error (MMSE) estimation is implemented on Eq. (13) to estimate values (). Denote the estimation of as . Put
Step 7. Find the index of the codeword corresponding to in the precoding codebook , which is given by where is the th dimensional codeword in the precoding codebook. Then, MS sends the index to the BS.
Step 8. The BS constructs the following precoder according to the index obtained by the feedback: where , .
Step 9. With the precoder constructed above, the BS first sends a pilot (namely 1) to the MS. The signal received by the MS can be expressed as Furthermore, the MS can estimate .
Step 10. The subsequent data signal is transmitted using the precoder and the estimated value of .
4. Simulations
In this section, we will demonstrate the performance of the proposed PCAbased adaptive trainingfeedback scheme. The channel established in the simulation of this paper is established according to Section 7 of [11]. The RootMeanSquare (RMS) delay spread is 10°, the delay distribution proportionality factor is 3, and the per cluster shadowing std is 3 dB. The scaling factors for azimuth angles generation and the scaling factors for elevation angles generation are 1.273 and 1.1764, respectively. The RMS azimuth spread of arrival angles (ASA) and departure angles (ASD) are 9° and 10°, respectively. The RMS elevation spread of arrival angles (ZSA) and departure angles (ZSD) are both 10°. The cluster ASA is 22°, the cluster ASD is 10°, and the cluster ZSA is 7°. In addition, we select the values of the remaining parameters as shown in Table 1.
Based on the above channel model, for the case of 8 trainings (), we first designed an 8dimensional precoding codebook containing 16 codewords, (see [17]), defined as where , and . .
The calculation formula of the throughput is as follows:
where denotes the precoder. In the following, when , it is called as “ideal feedback”, and when the BS uses the coefficients given in (14) to construct , it is called as “full feedback”.
Simulation 1. According to the above precoding codebook, when , in the case of 32 BS antennas and the speed of the MS being 120 km/h, as shown in Figure 1, we simulated the throughput in the three feedback situations: ideal feedback, full feedback, and limited feedback based on the above precoding codebook. Notice that the ideal feedback involved in simulations refers to full feeding back the true DL channel at the current time slot to the BS. It can be seen from the figure that at a speed of 120 km/h, the throughput of the proposed scheme under full feedback is close to the throughput under ideal feedback.
Simulation 2. In order to further demonstrate the effect of this scheme at high speed, taking 120 km/h and 8 trainings as examples, we compare the full feedback throughput of this scheme and the ideal feedback throughput under different BS antenna numbers (), as shown in Figure 2. This figure illustrates that even with a large number of BS antennas (), the throughput of the proposed scheme under full feedback is also close to that of ideal feedback.
5. Conclusions
In this paper, by utilizing the spatial reciprocity between the UL and DL channels in the FDD system, we propose a PCAbased adaptive trainingfeedback scheme for timevarying massive MIMO systems. The proposed scheme combines PCA technology and SVD method to extract the principal singular vectors of the covariance matrix of UL received signals. According to the spatial reciprocity of the FDD system, it is shown that singular vectors can constitute a subspace of the DL covariance matrix. Therefore, the pilots designed are based on these singular vectors and could be adaptive, which greatly reduces the number of training times and feedback bits. Simulation results show that the throughput of the proposed scheme is close to that of ideal feedback even under high speed, with low training and feedback overhead.
Data Availability
Data is available upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.