Downlink transmission and channel estimation for cell-free massive MIMO-OFDM with DSDs

In a cell-free massive MIMO system, multiple users arrive at multiple access points at separate times, while in an OFDM system, different delays can be equivalent to delay spread differences (DSDs). Since DSDs are not all the same, in the downlink transmission process, it is necessary to consider its impact on transmission techniques, such as channel estimation and downlink precoding. In this paper, aiming at the performance loss caused by DSDs in cell-free massive MIMO-OFDM system, we propose a multi-RB precoding optimization algorithm that maximizes the downlink sum rate. We derive the sum rate maximization problem into an iterative second-order cone programming form to achieve convex approximation. Then, considering the impact of DSDs on the accuracy of cell-free massive MIMO-OFDM channel estimation, we propose a downlink channel estimation method, which jointly uses channel state information reference signal and demodulation reference signal. The simulation results show that the proposed multi-RB optimal precoding can effectively improve the downlink sum rate, and the proposed downlink channel estimation can obtain accurate multi-RB frequency-domain channel parameters.

RAUs have separate propagation delays. These variations will cause the receiver to synchronize with the transmitted signal that arrives first, and DSDs will be generated on the asynchronous transmitted signal that arrives later. Especially in the OFDM system, accurate symbol time synchronization is an important prerequisite for removing the cyclic prefix (CP) and completing correct demodulation. DSDs will cause the received signal to produce a phase rotation that varies with the subcarrier index during the demodulation process, which will have a negative impact on system performance [5][6][7][8][9][10][11].
The major contributions of this paper are summarized as follows.
1. A multi-RB precoding optimization algorithm to maximize the cell-free massive MIMO system sum rate under the total RAU power constraint is proposed.

A DMRS channel estimation is proposed.
The novelty of this paper is to solve the problem that the performance loss of DSDs to broadband precoding increases as the number of subcarriers increases since per subcarrier precoding is not supported in the current 5G NR protocol [12].
The following notations are used. All boldface letters represent vectors (lower case) or matrices (upper case). The P × P identity matrix is denoted I P . The notations C N and C M×N refer to complex N-dimensional vectors and M × N matrices, respectively. The operators diag(·) , (·) T and (·) H denote diagonalization, transpose and conjugate transpose, respectively. A random vector x ∼ CN (m, ) is a complex Gaussian distribution with a mean vector m and covariance matrix .

Impact of the DSDs
In an OFDM system, obtaining accurate sampling of the transmitted signal within the symbol period is the prerequisite for the receiver to perform FFT. For this reason, it is necessary to perform symbol timing synchronization after removing CP. However, in the cell-free system as shown in Fig. 1, because the distances from RAUs to UEs are different, the receiver performs symbol timing synchronization with the signal that arrives first, and the signal that arrives subsequently will produce τ samples DSDs during the sampling process, as shown in Fig. 2. After FFT, τ samples DSDs result in a phase rotation of e −j2π nτ/N c on the nth subcarrier, where N c is the number of FFT points. It can be seen that phase rotation increases with the increase of subcarrier index, and the impact of phase rotation will be superimposed in the coherence bandwidth, causing additional frequency selective fading of the channel.

System model
We consider a cell-free massive MIMO-OFDM system, including M multi-antenna RAUs and K single-antenna UEs, where each RAU has L antennas. The locations of RAUs and UEs are randomly distributed within the coverage area. The total number of subcarriers in the system is N, UEs are grouped into N RB resource blocks (RBs), and each RB contains RB = N /N RB subcarriers. The number of UEs in each RB is much smaller than the total number of RAUs. Considering the downlink transmission, we model the received signal of UE k on the nth subcarrier as where represents the frequency-domain channel vector on the nth subcarrier from all RAUs to downlink UE k, s k,n ∼ CN (0, 1) is the data symbol sent by all RAUs to downlink UE k, z k,n ∼ CN 0, σ 2 k,n is additive white Gaussian noise, represents the precoding vector on the nth subcarrier for the data stream of downlink UE k, where w m,k ∈ C L×1 . Assume that the precoding granularity is ϕ , the same precoding vector is used on the adjacent N ϕ = ϕ RB subcarriers, and there are a total of R=N /N ϕ different precoding vectors in the system bandwidth. That is, for any n = (r − 1)N ϕ + 1, ..., rN ϕ , where r = 1, ..., R , we have w k,n = w k . Therefore, the received signal of downlink UE k within the system bandwidth can be modeled as 3 Proposed approach

Proposed precoding design
In this paper, we are interested in the problem of the sum rate maximization under total RAU power constraints. On the nth subcarrier, the SINR of downlink UE k is Then, the design problem is given by Since problem (7) is NP-hard, the globally optimal design mainly acts as a theoretical benchmark rather than a practical solution. Herein, motivated by [13,14], we develop a low-complexity algorithm that satisfies the necessary optimal conditions of (7). According to the monotonicity of the logarithmic function, (7) is equivalent to which can be equivalently rewritten as Note that the objective function in (9) admits an SOC representation. Since (9b) and (9d) are already convex forms, we mainly focus on dealing with the constraints in (9c).
First, (9c) can be rewritten as where k,n = h k,n h H k,n . We can see that (10) is also non-convex, since the right side of (10) has the form of quadratic-over-linear, it can be replaced by its first-order expansions [15]. Thus, we define where ≥ 0 and f ≥ α . We obtain the first-order Taylor expansion of (11) about a certain point w (a) , f (a) as From the above analysis, we can transform the constraint of (9c) into a convex form Finally, the original problem (8) can be reformulated as a convex approximate problem (14) which can be solved in Algorithm 1.

Downlink channel estimation
In this section, considering the impact of DSDs on the accuracy of cell-free massive MIMO-OFDM system channel estimation, we propose a high-precision multi-RB downlink channel estimation scheme, which provides the necessary channel information for multi-RB precoding optimization. In the pilot design of the 5G NR [16], CSI-RS is mainly used for channel sounding to obtain the path loss and DSD between UEs and RAUs. DMRS is mainly used for the demodulation of uplink and downlink data. The least squares (LS) algorithm is usually used to estimate the initial channel response of the pilot position, and then, interpolation filtering is performed in the frequency domain to obtain the channel response of the data position [17].
Regarding the CSI-RS of cell-free system, in uplink, sounding reference signal (SRS) is used to collect DSDs and large-scale fading; in downlink, tracking reference signal (TRS), a special CSI-RS, is used to track the phase deviation on the received signal.
For SRS, per UE tracking can be achieved in the 5G NR protocol. Therefore, the uplink channel estimation can be implemented in CPU. For downlink channel estimation, the signal overhead for per AP tracking is too large, and the implementation complexity is too high for UEs. Therefore, to estimate the statistical properties of downlink channels, we choose to use a composite channel tracking method, which is discussed in the following specific steps.
The specific steps of the multi-RB channel estimation we proposed are as follows: First, all RAUs transmit orthogonal CSI-RS, UE k estimates large-scale fading, β pilot, that is, the CSI-RS pilot adopts a single-port design. At this time, the receiver only estimates the composite signal from all RAUs to the UE. Specifically, multiple RAUs are regarded as multipath signals, and the power and delay of each path are estimated on UE k based on the distinguishable time-domain multipath signals. The single-port design can effectively reduce the pilot overhead and design complexity, but it will produce additional channel time/frequency selectivity. In the case of sufficient estimation accuracy, these two methods are theoretically equivalent. Second, considering that there are P equally spaced DMRS in the coherent bandwidth, after LS, the frequency-domain signal collected by UE k is For convenience of representation, the subscript k is omitted in the following notation. On subcarrier p, The precoding vector of the channel between UE k and RAU m is w m . Assuming that the frequency-domain channel remains unchanged on these P resource elements, the frequency-domain signal collected by UE k can be modeled as wherein Subsequently, using minimum mean square error (MMSE) channel estimation, the cross-correlation matrix of y n P and f is and the autocorrelation matrix of y N P is Finally, the MMSE estimate of f is [18] (15) y N P = y N 1 , y N 2 , ..., y N P T ∈ C P×1

Numerical results
In this section, some numerical examples are evaluated to show the performance of the proposed channel estimation and precoding optimization under various system settings.
We consider a cell-free massive MIMO-OFDM system with detailed simulation parameters listed in Table 1. Figures 3 and 4 compare the achievable sum rate of the system when using traditional conjugate precoding and the proposed optimized precoding, where the precoding granularity is 2RB, 4RB, 6RB and 8RB, respectively. In these two precoding schemes, the sum rate increases as the SNR increases and decreases with increasing granularity of the precoding. In addition, as expected, our proposed optimized precoding always performs better than conjugate precoding in terms of sum rate. When the precoding granularity is  2RB and 4RB, the sum rate is increased by approximately 20 percent, and when the precoding granularity is 6RB and 8RB, the sum rate is increased by approximately 10 percent. This result is because as the number of subcarriers increases, the impact of DSDs on system performance becomes increasingly serious. Since the minimum granularity of precoding in the existing 5G NR protocol is 2RB and is usually statically configured as 4RB, the proposed algorithm can reflect the robustness to DSDs in wideband precoding under the existing 5G NR protocol. Figure 5 shows the MSE vs. SNR performance comparison of the proposed downlink channel estimation when the channel estimation granularity is 2RB, 4RB, 8RB and 12RB. It can be observed that as the SNR increases, the MSE of the proposed channel estimation decreases, and the accuracy of the channel estimation increases. In addition, the accuracy of the proposed channel estimation improves with increasing channel estimation granularity. This result is due to the use of DMRS interpolation intervals of the same width. The larger the bandwidth of the estimated channel is, the greater the number of DMRS participating in the channel estimation. Since the CSI-RS pilot can be designed in a single-port mode, the pilot overhead and design complexity are effectively reduced, thereby realizing low-complexity accurate estimation of the wideband channel with DSDs. Figure 6 shows the performance comparison when the UE estimates the path loss and time delay based on the single-port CSI-RS sent by all RAUs, and different numbers of RAUs send DMRS. Assuming that the number of RAUs sending DMRS is 2, 4, 6 and 8 and the granularity of frequency-domain joint channel estimation is 2RB, 12RB and 24RB, it can be seen that the mismatch of CSI-RS and DMRS has little effect on the accuracy of channel estimation. Especially when the signal-to-noise ratio is low and the channel estimation granularity is small, the loss caused by this mismatch can be ignored.

Conclusion
In this paper, we have studied the DSDs problem in a cell-free massive MIMO-OFDM system due to its distributed characteristics, which causes the received signal to produce phase rotation related to the subcarrier index during demodulation. Since the existing 5G NR protocol does not support per subcarrier precoding, to address the impact of DSDs, we designed a multi-RB precoding optimization problem and used convex approximation to solve the problem. Then, to cope with the influence of the DSDs on the accuracy of channel estimation, we propose a downlink channel estimation method, which uses CSI-RS and DMRS jointly to realize multi-RB channel estimation. The simulation results show that, compared with conjugate precoding, the proposed optimal precoding can effectively improve system performance. In addition, the proposed downlink channel estimation can obtain accurate channel information.