This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Joint Precoding and Fronthaul Compression for Cell-Free MIMO Downlink With Radio Stripes

Sangwon Jo, Student Member, IEEE, Hoon Lee, Member, IEEE,
and Seok-Hwan Park, Member, IEEE
This work was supported by the National Research Foundation (NRF) of Korea funded by the MOE under Grant 2019R1A6A1A09031717, and MSIT under Grants 2021R1I1A3054575 and RS-2023-00238977, and in part by by Korea Research Institute for defense Technology planning and advancement (KRIT) grant funded by the Korea government (DAPA(Defense Acquisition Program Administration)) (21-106-A00-007, Space-Layer Intelligent Communication Network Laboratory, 2022). S. Jo and S.-H. Park are with the Division of Electronic Engineering, Jeonbuk National University, Jeonju, Korea (email: [email protected], [email protected]). H. Lee is with the Department of Information and Communications Engineering, Pukoyng National University, Busan, Korea ([email protected]). Copyright (c) 2023 IEEE. Personal use of this material is permitted. However, permission to use this material for any other purposes must be obtained from the IEEE by sending a request to [email protected].
Abstract

A sequential fronthaul network, referred to as radio stripes, is a promising fronthaul topology of cell-free MIMO systems. In this setup, a single cable suffices to connect access points (APs) to a central processor (CP). Thus, radio stripes are more effective than conventional star fronthaul topology which requires dedicated cables for each of APs. Most of works on radio stripes focused on the uplink communication or downlink energy transfer. This work tackles the design of the downlink data transmission for the first time. The CP sends compressed information of linearly precoded signals to the APs on fronthaul. Due to the serial transfer on radio stripes, each AP has an access to all the compressed blocks which pass through it. Thus, an advanced compression technique, called Wyner-Ziv (WZ) compression, can be applied in which each AP decompresses all the received blocks to exploit them for the reconstruction of its desired precoded signal as side information. The problem of maximizing the sum-rate is tackled under the standard point-to-point (P2P) and WZ compression strategies. Numerical results validate the performance gains of the proposed scheme.

Index Terms:
Cell-free massive MIMO, radio stripe fronthaul network, finite-capacity fronthaul.

I Introduction

Cell-free massive MIMO network is envisioned to mitigate the impact of inter-cell interference signals of cell boundary users. This can be made possible by moving the channel encoding and decoding functionality from distributed access points (APs) to a central processor (CP) [1]. An essential requirement for cell-free MIMO systems is fronthaul networks which realize coordination among APs and CP through, for instance, optical fiber cables. Many early works have considered a star fronthaul topology where each AP has a dedicated fronthaul connection to the CP [1, 2]. Albeit the success in the communication performance, the star fronthaul topology poses practical challenges as it requires a number of high-capacity long-length cables to realize individual fronthaul interconnection for all APs.

A bus fronthaul topology, also referred to as radio stripes, has been regarded as a promising solution to implement practical cost-effective fronthaul systems [3, 4, 5, 6]. In this setup, distributed APs are sequentially connected to the CP with a single cable. Compared to conventional star fronthaul topology, the radio stripe enables an easier deployment of cell-free MIMO systems, particularly for dense areas such as sports arenas and railway stations.

The viability of the radio stripe architecture has been investigated for uplink cell-free MIMO networks [3, 4]. The design of linear processing at APs was addressed in [3]. The fronthaul links are simply assumed to deliver soft estimates of the uplink transmit signals and additional side information. Compress-and-forward relaying strategies were adopted in [4] to exploit finite-capacity radio stripes. The downlink cell-free MIMO systems with the radio stripes have also been studied [5, 6], but they have been confined to specific application scenarios such as energy transfer [5] and a single receiving user equipment (UE) [6]. These scenarios have no adequate investigations on inter-UE interference management. Thus, the signal processing strategies of these conventional works cannot be straightforwardly extended to generic multi-UE downlink cell-free MIMO networks.

For the first time, this paper considers a generic downlink cell-free MIMO system with the radio stripes where separate APs, connected through the radio stripe fronthaul networks, provide donwlink communication services for multiple UEs. Unlike the scenarios studied in [5, 6], which focused on energy transfer or single-UE setups, the design of fronthaul and wireless links should account for the impact of inter-UE interference signals. To this end, the CP first performs a cooperative precoding across APs and sends the compressed information of precoded signals. Since the compressed blocks are serially transferred on the radio stripes through APs, each AP has a chance of decompressing all the blocks that pass through it. To realize this novel fronthauling policy, we exploit the notion of the Wyner-Ziv (WZ) compression [7, 8], in which each AP decompresses all the received compressed blocks so as to utilize them as side information for decompressing its transmit signal. We tackle the sum-rate maximization problems with the standard point-to-point (P2P) compression [9, Sec. 3.6] and WZ compression techniques [7, 8] under the constraints on the per-AP transmit powers and fronthaul capacity. To efficiently address this nonconvex problem, iterative algorithms based on the weighted minimum mean squared error (WMMSE) approach [10] and a convexification method [11] are derived. The performance gains of the WZ compression scheme compared to the P2P compression are investigated via numerical results.

II System Model

Refer to caption
Figure 1: Illustration of the downlink of cell-free MIMO system with a radio stripe fronthaul network. A single fronthaul cable connects the CP to AP LL, AP L1L-1, \ldots, AP 2, and AP 1.

We consider the downlink of a cell-free MIMO system in which KK single-antenna UEs are served by a CP through LL APs each equipped with NN antennas. Let 𝒦={1,2,,K}\mathcal{K}=\{1,2,\ldots,K\} and ={1,2,,L}\mathcal{L}=\{1,2,\ldots,L\} be sets of UE and AP indices, respectively. As shown in Fig. 1, LL APs and the CP are connected by the radio stripes [3] so that the fronthaul connection goes from the CP to AP LL, AP L1L-1, …, to AP 11. For simplicity, we denote the CP as AP L+1L+1. Let 𝙵i\mathtt{F}_{i} be the fronthaul link from AP ii to AP i1i-1 by 𝙵i\mathtt{F}_{i} for i{2,3,,L,L+1}i\in\{2,3,\ldots,L,L+1\}. Then, 𝙵L+1\mathtt{F}_{L+1} represents the fronthaul link from the CP to AP LL. Every fronthaul link 𝙵i\mathtt{F}_{i} has a finite capacity CFC_{F} in bits/symbol, i.e., the number of bits delivered per symbol period.

Let 𝐡k,i𝒞𝒩(𝟎,ρk,i𝐈)N\mathbf{h}_{k,i}\sim\mathcal{CN}(\mathbf{0},\rho_{k,i}\mathbf{I})\in\mathbb{C}^{N} be the Rayleigh fading channel vector from AP ii to UE kk with the path-loss ρk,i\rho_{k,i}. Then, received signal yky_{k}\in\mathbb{C} of UE kk is given as

yk=𝐡kH𝐱+zk,\displaystyle y_{k}=\mathbf{h}_{k}^{H}\mathbf{x}+z_{k}, (1)

where 𝐡k=[𝐡k,1H𝐡k,LH]HNL×1\mathbf{h}_{k}=[\mathbf{h}_{k,1}^{H}\,\cdots\,\mathbf{h}_{k,L}^{H}]^{H}\in\mathbb{C}^{NL\times 1} denotes the channel vector from all APs to UE kk, 𝐱=[𝐱1H𝐱LH]HNL×1\mathbf{x}=[\mathbf{x}_{1}^{H}\,\cdots\,\mathbf{x}_{L}^{H}]^{H}\in\mathbb{C}^{NL\times 1} stacks the transmit signal vectors all APs with 𝐱iN×1\mathbf{x}_{i}\in\mathbb{C}^{N\times 1} being the signal transmitted by AP ii, and zk𝒞𝒩(0,σz2)z_{k}\sim\mathcal{CN}(0,\sigma_{z}^{2}) represents the additive noise. The transmit signal of each AP ii is subject to a power constraint

𝔼[𝐱i2]Ptx.\displaystyle\mathbb{E}\left[\|\mathbf{x}_{i}\|^{2}\right]\leq P_{\text{tx}}. (2)

We define the transmit signal-to-noise ratio (SNR) of the downlink channel as Ptx/σz2P_{\text{tx}}/\sigma_{z}^{2}

The CP designs the signal processing strategy of all APs based on the estimated channel state information (CSI). Let 𝐡^k,i\hat{\mathbf{h}}_{k,i} be an estimate of 𝐡k,i\mathbf{h}_{k,i}, which can be written as

𝐡k,i=𝐡^k,i+𝐞k,i,\displaystyle\mathbf{h}_{k,i}=\hat{\mathbf{h}}_{k,i}+\mathbf{e}_{k,i}, (3)

with the CSI error vector 𝐞k,i\mathbf{e}_{k,i}. With the linear minimum mean squared error estimator, 𝐞k,i\mathbf{e}_{k,i} becomes uncorrelated with 𝐡^k,i\hat{\mathbf{h}}_{k,i} and distributed as 𝐞k,i𝒞𝒩(𝟎,βk,i𝐈)\mathbf{e}_{k,i}\sim\mathcal{CN}(\mathbf{0},\beta_{k,i}\mathbf{I}) [12]. The error variance βk,i[0,ρk,i]\beta_{k,i}\in[0,\rho_{k,i}] increases as the pilot length gets shorter or the SNR of the uplink training channel decreases. Let us define the nominal channel vector 𝐡^k=[𝐡^k,1H𝐡^k,2H𝐡^k,LH]H\hat{\mathbf{h}}_{k}=[\hat{\mathbf{h}}_{k,1}^{H}\hat{\mathbf{h}}_{k,2}^{H}\cdots\hat{\mathbf{h}}_{k,L}^{H}]^{H} and CSI error vector 𝐞k=[𝐞k,1H𝐞k,2H𝐞k,LH]H\mathbf{e}_{k}=[\mathbf{e}_{k,1}^{H}\mathbf{e}_{k,2}^{H}\cdots\mathbf{e}_{k,L}^{H}]^{H} for UE kk. Then, the error vector 𝐞k\mathbf{e}_{k} is distributed as 𝐞k𝒞𝒩(𝟎,𝐄k)\mathbf{e}_{k}\sim\mathcal{CN}(\mathbf{0},\mathbf{E}_{k}) with 𝐄k=blkdiag({βk,i𝐈}i)\mathbf{E}_{k}=\text{blkdiag}(\{\beta_{k,i}\mathbf{I}\}_{i\in\mathcal{L}}).

The downlink transmission from the CP to the UEs consists of two phases: i) CP-to-APs fronthauling phase; and ii) APs-to-UEs wireless accessing phase. In the fronthauling phase, the CP sends bit streams, which represent quantized versions of the precoded baseband signals, to the APs via the radio stripes. In the wireless accessing phase, the APs simultaneously transmit the downlink signals over wireless channels toward UEs, where the downlink signals are obtained by decompressing the bit streams received through the fronthaul links. We denote the time durations of the fronthauling and wireless accessing transmission phases by nFn_{F} and nAn_{A}, respectively, which are measured in the symbol period. Thus, nAn_{A} becomes the number of channel uses for downlink transmission which corresponds to the blocklength of the channel coding. The ratio between two phases is defined as r=nF/nAr=n_{F}/n_{A}.

III Precoding and Fronthaul Compression With Radio Stripes

This section proposes a sequential fronthauling strategy to compute efficient precoding solutions through the radio stripes. Let sk𝒞𝒩(0,1)s_{k}\sim\mathcal{CN}(0,1) be the data symbol intended to UE kk. To obtain the transmit signal vector 𝐱\mathbf{x} in (1), the CP first generates an intermediate vector 𝐱~NL×1\tilde{\mathbf{x}}\in\mathbb{C}^{NL\times 1} given by

𝐱~=k𝒦𝐯ksk[𝐱~1H𝐱~LH]H,\displaystyle\tilde{\mathbf{x}}=\sum\nolimits_{k\in\mathcal{K}}\mathbf{v}_{k}s_{k}\triangleq[\tilde{\mathbf{x}}_{1}^{H}\ \cdots\ \tilde{\mathbf{x}}_{L}^{H}]^{H}, (4)

where 𝐯kNL×1\mathbf{v}_{k}\in\mathbb{C}^{NL\times 1} is the linear precoding vector for UE kk and the ii-th subvector 𝐱~iN×1\tilde{\mathbf{x}}_{i}\in\mathbb{C}^{N\times 1} of 𝐱~\tilde{\mathbf{x}} indicates the transmit signal at AP ii. For the sequential fronthauling, the CP delivers the signal vector 𝐱~i\tilde{\mathbf{x}}_{i} to AP ii through APs L,L1,,i+1L,L-1,\cdots,i+1. Since the fronthaul links have finite capacity CFC_{F} bits/symbol, each precoded signal 𝐱~i\tilde{\mathbf{x}}_{i} is quantized and compressed. Thus, the signal vector 𝐱i\mathbf{x}_{i} transmitted by AP ii is a quantized version of 𝐱^i\hat{\mathbf{x}}_{i} and is modeled as

𝐱i=𝐱~i+𝐪i.\displaystyle\mathbf{x}_{i}=\tilde{\mathbf{x}}_{i}+\mathbf{q}_{i}. (5)

Under the Gaussian test channel [11], the quantization noise 𝐪i\mathbf{q}_{i} is independent of 𝐱~i\tilde{\mathbf{x}}_{i} and distributed as 𝐪i𝒞𝒩(𝟎,𝛀i)\mathbf{q}_{i}\sim\mathcal{CN}(\mathbf{0},\boldsymbol{\Omega}_{i}). The power constraint (2) is rewritten as

k𝒦𝐃iH𝐯k2+tr(𝛀i)Ptx,\displaystyle\sum\nolimits_{k\in\mathcal{K}}\left\|\mathbf{D}_{i}^{H}\mathbf{v}_{k}\right\|^{2}+\text{tr}\left(\boldsymbol{\Omega}_{i}\right)\leq P_{\text{tx}}, (6)

where 𝐃iNL×N\mathbf{D}_{i}\in\mathbb{C}^{NL\times N} is filled with zeros except for the rows from N(i1)+1N(i-1)+1 to NiNi being an identity matrix of size NN.

𝙵L+1(CP  AP L)\!\!\!\!\!\!\begin{array}[]{c}\mathtt{F}_{L+1}\\ \text{\scriptsize(CP $\to$ AP $L$)}\end{array}\!\!\!\!\!\! 𝙵L(AP L  AP L1)\!\!\!\!\!\!\begin{array}[]{c}\mathtt{F}_{L}\\ \text{\scriptsize(AP $L$ $\to$ AP $L-1$)}\end{array}\!\!\!\!\!\! \cdots 𝙵1(AP 2  AP 1)\!\!\!\!\!\!\begin{array}[]{c}\mathtt{F}_{1}\\ \text{\scriptsize(AP $2$ $\to$ AP $1$)}\end{array}\!\!\!\!\!\!
Slot 1 𝙱1\mathtt{B}_{1} \cdots
Slot 2 𝙱2\mathtt{B}_{2} 𝙱1\mathtt{B}_{1} \ddots \vdots
\vdots \vdots \vdots \ddots
Slot LL 𝙱L\mathtt{B}_{L} 𝙱L1\mathtt{B}_{L-1} \cdots 𝙱1\mathtt{B}_{1}

Table I: Compressed blocks 𝙱1,,𝙱L\mathtt{B}_{1},\ldots,\mathtt{B}_{L} transferred on the fronthaul links 𝙵1,,𝙵L\mathtt{F}_{1},\ldots,\mathtt{F}_{L} at time slots t=1,2,,Lt=1,2,\ldots,L

As illustrated in Table I, the time-division multiple access protocol is adopted for the fronthauling phase. We define 𝙱i\mathtt{B}_{i} as the compressed block of the quantized signal 𝐱i\mathbf{x}_{i}, which needs to be delivered to AP ii. The time duration of the fronthauling phase is divided into LL time slots, where each time slot spans nF/Ln_{F}/L symbol periods. In time slot 1, the fronthaul link 𝙵L+1\mathtt{F}_{L+1} connecting the CP to AP LL carries the block 𝙱1\mathtt{B}_{1} which is intended for AP 11. Subsequently, in time slot 2, 𝙱1\mathtt{B}_{1} is transferred from AP LL to AP L1L-1 on the fronthaul 𝙵L\mathtt{F}_{L}. At the same time, AP L1L-1 receives 𝙱2\mathtt{B}_{2} from the CP on 𝙵L+1\mathtt{F}_{L+1}. In general, in time slot tt (t=1,2,,Lt=1,2,\ldots,L), the block 𝙱i\mathtt{B}_{i} with iti\leq t moves on the fronthaul link 𝙵L+it+1\mathtt{F}_{L+i-t+1} from AP L+it+1L+i-t+1 to AP L+itL+i-t.

Once the fronthaul transmission over the LL slots is completed, each AP ii has an access to its coded block 𝙱i\mathtt{B}_{i}. Thus, it can decompress 𝙱i\mathtt{B}_{i} to obtain the transmit signal vector 𝐱i\mathbf{x}_{i}, which is a quantized version of 𝐱~i\tilde{\mathbf{x}}_{i} as in (5), and transmits it over the downlink wireless channel. Thus, the received signal yky_{k} in (1) can be expressed as

yk\displaystyle y_{k} =𝐡^kH𝐯ksk+l𝒦{k}𝐡^kH𝐯lsl+l𝒦𝐞kH𝐯lsl+𝐡kH𝐪+zk,\displaystyle=\hat{\mathbf{h}}_{k}^{H}\mathbf{v}_{k}s_{k}+\sum\nolimits_{l\in\mathcal{K}\setminus\{k\}}\hat{\mathbf{h}}_{k}^{H}\mathbf{v}_{l}s_{l}+\sum\nolimits_{l\in\mathcal{K}}\mathbf{e}_{k}^{H}\mathbf{v}_{l}s_{l}+\mathbf{h}_{k}^{H}\mathbf{q}+z_{k}, (7)

where 𝐪=[𝐪1H𝐪LH]H\mathbf{q}=[\mathbf{q}_{1}^{H}\cdots\mathbf{q}_{L}^{H}]^{H} stands for the quantization noise in the fronthauling phase distributed as 𝐪𝒞𝒩(𝟎,𝛀¯)\mathbf{q}\sim\mathcal{CN}(\mathbf{0},\bar{\boldsymbol{\Omega}}) with 𝛀¯=blkdiag({𝛀i}i)\bar{\boldsymbol{\Omega}}=\text{blkdiag}(\{\boldsymbol{\Omega}_{i}\}_{i\in\mathcal{L}}). In (7), the second and third terms stand for the inter-UE interference and the noise induced by the CSI estimation error 𝐞k\mathbf{e}_{k}, respectively. The achievable data rate of UE kk can be written as [12]

Rk=log2(1+|𝐡^kH𝐯k|2𝙸𝙽k(𝐯,𝛀)),\displaystyle R_{k}=\log_{2}\left(1+\frac{|\hat{\mathbf{h}}_{k}^{H}\mathbf{v}_{k}|^{2}}{\mathtt{IN}_{k}\left(\mathbf{v},\boldsymbol{\Omega}\right)}\right), (8)

where the interference-plus-noise power, denoted by 𝙸𝙽k(𝐯,𝛀)\mathtt{IN}_{k}(\mathbf{v},\boldsymbol{\Omega}), is defined as

𝙸𝙽k(𝐯,𝛀)\displaystyle\mathtt{IN}_{k}\left(\mathbf{v},\boldsymbol{\Omega}\right) =l𝒦{k}|𝐡^kH𝐯l|2+l𝒦𝐯lH𝐄k𝐯l\displaystyle=\sum\nolimits_{l\in\mathcal{K}\setminus\{k\}}\left|\hat{\mathbf{h}}_{k}^{H}\mathbf{v}_{l}\right|^{2}+\sum\nolimits_{l\in\mathcal{K}}\mathbf{v}_{l}^{H}\mathbf{E}_{k}\mathbf{v}_{l}
+tr((𝐡^k𝐡^kH+𝐄k)𝛀¯)+σz2\displaystyle+\text{tr}\left(\left(\hat{\mathbf{h}}_{k}\hat{\mathbf{h}}_{k}^{H}+\mathbf{E}_{k}\right)\bar{\boldsymbol{\Omega}}\right)+\sigma_{z}^{2} (9)

with 𝛀={𝛀k}i\boldsymbol{\Omega}=\{\boldsymbol{\Omega}_{k}\}_{i\in\mathcal{L}}.

IV Fronthaul Compression Strategies and Problem Definition

In this section, we describe two fronthaul compression strategies: i) Wyner-Ziv (WZ) compression [7]; and ii) point-to-point (P2P) compression [9, Sec. 3.6]. It is then followed by the description of the joint optimization task for fronthauling and accessing phases.

IV-A WZ Compression

Due to the nature of the sequential fronthauling phase, AP ii receives the preceding blocks 𝙱1,𝙱2,,𝙱i1\mathtt{B}_{1},\mathtt{B}_{2},\ldots,\mathtt{B}_{i-1} before receiving its desired block 𝙱i\mathtt{B}_{i}. Consequently, AP ii has a chance of recovering the signals 𝐱1,𝐱2,,𝐱i1\mathbf{x}_{1},\mathbf{x}_{2},\ldots,\mathbf{x}_{i-1}. As depicted in Fig. 2(a), these preceding signals can be exploited as side information for the decompression of 𝙱i\mathtt{B}_{i} The WZ compression, also called distributed compression, can be applied to the design of compressor and decompressor for this scenario which enables the compressor at the CP to use a finer quantizer [7, 8]. The number of bits, i.e., the minimal size of block 𝙱i\mathtt{B}_{i}, needed for successful decompression with the WZ compression scheme is given as nAgWZ,i(𝐯,𝛀)n_{A}\cdot g_{\text{WZ},i}(\mathbf{v},\boldsymbol{\Omega}) bits, where gWZ,i(𝐯,𝛀)g_{\text{WZ},i}(\mathbf{v},\boldsymbol{\Omega}) measures the compression rate in the number of bits per sample. The function gWZ,i(𝐯,𝛀)g_{\text{WZ},i}(\mathbf{v},\boldsymbol{\Omega}) is defined as

gWZ,i(𝐯,𝛀)\displaystyle g_{\text{WZ},i}\left(\mathbf{v},\boldsymbol{\Omega}\right) =I(𝐱~i;𝐱i|𝐱1,𝐱2,,𝐱i1)\displaystyle=I\left(\tilde{\mathbf{x}}_{i};\mathbf{x}_{i}\big{|}\mathbf{x}_{1},\mathbf{x}_{2},\ldots,\mathbf{x}_{i-1}\right)
=(a)log2det(𝐃¯iH𝐕Σ𝐃¯i+𝛀¯i)log2det(𝐃¯i1H𝐕Σ𝐃¯i1+𝛀¯i1)\displaystyle\overset{(\text{a)}}{=}\log_{2}\det\left(\bar{\mathbf{D}}_{i}^{H}\mathbf{V}_{\Sigma}\bar{\mathbf{D}}_{i}+\bar{\boldsymbol{\Omega}}_{i}\right)-\log_{2}\det\!\left(\bar{\mathbf{D}}_{i-1}^{H}\mathbf{V}_{\Sigma}\bar{\mathbf{D}}_{i-1}+\bar{\boldsymbol{\Omega}}_{i-1}\right)
log2det(𝛀i),\displaystyle-\log_{2}\!\det\!\left(\boldsymbol{\Omega}_{i}\right), (10)

where 𝐕Σl𝒦𝐯l𝐯lH\mathbf{V}_{\Sigma}\triangleq\sum_{l\in\mathcal{K}}\mathbf{v}_{l}\mathbf{v}_{l}^{H}. Here, we have defined the matrices 𝐃¯i=[𝐃1𝐃2𝐃i]NL×Ni\bar{\mathbf{D}}_{i}=[\mathbf{D}_{1}\,\mathbf{D}_{2}\,\cdots\,\mathbf{D}_{i}]\in\mathbb{C}^{NL\times Ni} and 𝛀¯i=blkdiag(𝛀1,𝛀2,,𝛀i)Ni×Ni\bar{\boldsymbol{\Omega}}_{i}=\text{blkdiag}(\boldsymbol{\Omega}_{1},\boldsymbol{\Omega}_{2},\ldots,\boldsymbol{\Omega}_{i})\in\mathbb{C}^{Ni\times Ni}. The equality (a) can be shown by using the chain rule of mutual information.

The block 𝙱i\mathtt{B}_{i} of length nAgWZ,i(𝐯,𝛀)n_{A}\cdot g_{\text{WZ},i}(\mathbf{v},\boldsymbol{\Omega}) bits is delivered via each fronthaul link of capacity CFC_{F} bit/symbol over nF/Ln_{F}/L symbol periods. Thus, for successful decompression of 𝙱i\mathtt{B}_{i} with the WZ compression, the beamforming vector 𝐯={𝐯k}k𝒦\mathbf{v}=\{\mathbf{v}_{k}\}_{k\in\mathcal{K}} and the quantization noise covariance matrix 𝛀\boldsymbol{\Omega} should satisfy the following fronthaul capacity constraint:

gWZ,i(𝐯,𝛀)rCFL.\displaystyle g_{\text{WZ},i}(\mathbf{v},\boldsymbol{\Omega})\leq\frac{rC_{F}}{L}. (11)

IV-B P2P Compression

As shown in Fig. 2(b), in the P2P compression scheme, AP ii only decodes its desired block 𝙱i\mathtt{B}_{i} while simply bypassing others 𝙱j\mathtt{B}_{j}, j<i\forall j<i. Therefore, the P2P compression strategy, which does not exploit the side information about the preceding APs, can be viewed as a special case of the WZ approach. Likewise (11), in the P2P compression scheme, the fronthaul capacity constraint can be expressed as

gP2P,i(𝐯,𝛀)rCFL.\displaystyle g_{\text{P2P},i}\left(\mathbf{v},\boldsymbol{\Omega}\right)\leq\frac{rC_{F}}{L}. (12)

where the compression rate function for the P2P compression scheme gP2P,i(𝐯,𝛀)g_{\text{P2P},i}(\mathbf{v},\boldsymbol{\Omega}) is given by

gP2P,i(𝐯,𝛀)\displaystyle g_{\text{P2P},i}\left(\mathbf{v},\boldsymbol{\Omega}\right) =I(𝐱~i;𝐱i)\displaystyle=I\left(\tilde{\mathbf{x}}_{i};\mathbf{x}_{i}\right)
=log2det(𝐃iH𝐕Σ𝐃i+𝛀i)log2det(𝛀i).\displaystyle=\log_{2}\det\left(\mathbf{D}_{i}^{H}\mathbf{V}_{\Sigma}\mathbf{D}_{i}+\boldsymbol{\Omega}_{i}\right)-\log_{2}\det\left(\boldsymbol{\Omega}_{i}\right). (13)
Refer to caption
Refer to caption
Figure 2: Comparison of the decompression operations at the APs with WZ and P2P compression schemes

IV-C Problem Description

To optimize the fronthauling and accessing phases jointly, we tackle the problem of optimizing the precoding vector 𝐯\mathbf{v} and the fronthaul quantization policy 𝛀\boldsymbol{\Omega}, with the goal of maximizing the sum-rate k𝒦Rk\sum_{k\in\mathcal{K}}R_{k}. The sum-rate maximization problem is stated as

max.𝐯,𝛀,𝐑\displaystyle\underset{\mathbf{v},\boldsymbol{\Omega},\mathbf{R}}{\mathrm{max.}}\,\,\, k𝒦Rk\displaystyle\sum\nolimits_{k\in\mathcal{K}}R_{k}\, (14a)
s.t.\displaystyle\mathrm{s.t.}\,\,\,\,\,\, Rklog2(1+|𝐡^kH𝐯k|2𝙸𝙽k(𝐯,𝛀)),k𝒦,\displaystyle R_{k}\leq\log_{2}\left(1+\frac{|\hat{\mathbf{h}}_{k}^{H}\mathbf{v}_{k}|^{2}}{\mathtt{IN}_{k}\left(\mathbf{v},\boldsymbol{\Omega}\right)}\right),\,k\in\mathcal{K}, (14b)
gX,i(𝐯,𝛀)rCFL,i,\displaystyle g_{X,i}\left(\mathbf{v},\boldsymbol{\Omega}\right)\leq\frac{rC_{F}}{L},\,i\in\mathcal{L}, (14c)
k𝒦𝐃iH𝐯k2+tr(𝛀i)Ptx,i,\displaystyle\sum\nolimits_{k\in\mathcal{K}}\left\|\mathbf{D}_{i}^{H}\mathbf{v}_{k}\right\|^{2}+\text{tr}\left(\boldsymbol{\Omega}_{i}\right)\leq P_{\text{tx}},\,i\in\mathcal{L}, (14d)

where 𝐑={𝐑k}k𝒦\mathbf{R}=\{\mathbf{R}_{k}\}_{k\in\mathcal{K}} denotes the achievable rate vector and left-hand-side gX,i(𝐯,𝛀)g_{X,i}\left(\mathbf{v},\boldsymbol{\Omega}\right) of the fronthaul capacity constraint in (14b) can be chosen as X{WZ,P2P}X\in\{\text{WZ},\text{P2P}\}, i.e., it is set to either (11) or (12). Due to the constraints in (14b) and (14c), the above problem becomes nonconvex in general, and thus it is difficult to obtain the globally optimal solution. It is not difficult to see that gWZ,i(𝐯,𝛀)gP2P,i(𝐯,𝛀)g_{\text{WZ},i}(\mathbf{v},\boldsymbol{\Omega})\leq g_{\text{P2P},i}(\mathbf{v},\boldsymbol{\Omega}), and thus (12) is stricter than (11). Consequently, the feasible set of problem (14) for the WZ scheme always contains that of the P2P compression, implying that the WZ compression is superior to the P2P scheme.

V Optimization Algorithms

This section presents optimization algorithms to solve problem (14) for the WZ and P2P compression schemes.

V-A WZ Compression

We first consider the WZ compression where constraint (14c) boils down to (11). The problem (14) is nonconvex due to the constraints (14b) and (14c). We handle the nonconvex access link rate constraint (14b) with the WMMSE approach (see, e.g., [10]). The constraint (14b) is satisfied if there exist uku_{k}\in\mathbb{C} and wk>0w_{k}>0 that satisfy

Rklog2wk1ln2wkek(𝐯,𝛀,uk)+1ln2,k𝒦,\displaystyle R_{k}\leq\log_{2}w_{k}-\frac{1}{\ln 2}w_{k}\,e_{k}(\mathbf{v},\boldsymbol{\Omega},u_{k})+\frac{1}{\ln 2},\,k\in\mathcal{K}, (15)

where uku_{k} and wkw_{k} indicate scalar receiver filter and weight for UE kk, respectively, and the mean squared error (MSE) function ek(𝐯,𝛀,uk)e_{k}(\mathbf{v},\boldsymbol{\Omega},u_{k}) is defined as

ek(𝐯,𝛀,uk)\displaystyle e_{k}(\mathbf{v},\boldsymbol{\Omega},u_{k}) =𝔼[|skukHyk|2]\displaystyle=\mathbb{E}\left[\big{|}s_{k}-u_{k}^{H}y_{k}\big{|}^{2}\right]
=|1ukH𝐡^kH𝐯k|2+|uk|2𝙸𝙽k(𝐯,𝛀).\displaystyle=\big{|}1-u_{k}^{H}\hat{\mathbf{h}}_{k}^{H}\mathbf{v}_{k}\big{|}^{2}+|u_{k}|^{2}\mathtt{IN}_{k}\left(\mathbf{v},\boldsymbol{\Omega}\right). (16)

The condition (15) becomes equivalent to the original rate constraint (14b) when the receive filter uku_{k} and the weight wkw_{k} are respectively set to

uk=𝐡^kH𝐯k|𝐡^kH𝐯k|2+𝙸𝙽k(𝐯,𝛀) and wk=1ek(𝐯,𝛀,uk).\displaystyle u_{k}=\frac{\hat{\mathbf{h}}_{k}^{H}\mathbf{v}_{k}}{\big{|}\hat{\mathbf{h}}_{k}^{H}\mathbf{v}_{k}\big{|}^{2}+\mathtt{IN}_{k}(\mathbf{v},\boldsymbol{\Omega})}\text{ and }w_{k}=\frac{1}{e_{k}(\mathbf{v},\boldsymbol{\Omega},u_{k})}. (17)

The nonconvex constraint in (14c) with X=WZX=\text{WZ} imposes that (11) is satisfied for all ii\in\mathcal{L}. By substituting the constraint (11) with i=ji=j into that with i=j+1i=j+1 for j=1,2,,L1j=1,2,\ldots,L-1 sequentially, it becomes

log2det(𝐃¯iH𝐕Σ𝐃¯i+𝛀¯i)j=1ilog2det(𝛀j)irCFL,i.\displaystyle\log_{2}\det\left(\bar{\mathbf{D}}_{i}^{H}\mathbf{V}_{\Sigma}\bar{\mathbf{D}}_{i}+\bar{\boldsymbol{\Omega}}_{i}\right)-\sum\nolimits_{j=1}^{i}\log_{2}\det\left(\boldsymbol{\Omega}_{j}\right)\leq\frac{irC_{F}}{L},\,i\in\mathcal{L}. (18)

It has been reported in [11] that (18) holds if there exist a group of positive definite matrices 𝚯iNi×Ni\boldsymbol{\Theta}_{i}\in\mathbb{C}^{Ni\times Ni}, ii\in\mathcal{L}, that satisfy

log2det(𝚯i)+1ln2tr(𝚯i1(𝐃¯iH𝐕Σ𝐃¯i+𝛀¯i))Niln2\displaystyle\log_{2}\det\left(\boldsymbol{\Theta}_{i}\right)+\frac{1}{\ln 2}\text{tr}\left(\boldsymbol{\Theta}_{i}^{-1}\left(\bar{\mathbf{D}}_{i}^{H}\mathbf{V}_{\Sigma}\bar{\mathbf{D}}_{i}+\bar{\boldsymbol{\Omega}}_{i}\right)\right)-\frac{Ni}{\ln 2}
j=1ilog2det(𝛀j)irCFL,i.\displaystyle\,\,\,\,\,\,-\sum\nolimits_{j=1}^{i}\log_{2}\det\left(\boldsymbol{\Omega}_{j}\right)\leq\frac{irC_{F}}{L},\,i\in\mathcal{L}. (19)

We can make the constraint (19) equivalent to the original one (11) by setting 𝚯i\boldsymbol{\Theta}_{i} to

𝚯i=𝐃¯iH𝐕Σ𝐃¯i+𝛀¯i.\displaystyle\boldsymbol{\Theta}_{i}=\bar{\mathbf{D}}_{i}^{H}\mathbf{V}_{\Sigma}\bar{\mathbf{D}}_{i}+\bar{\boldsymbol{\Omega}}_{i}. (20)

Thus, for the WZ compression scheme, the sum-rate maximization task (14) can be equivalently transformed into

max.𝐯,𝛀,𝐑,𝐮,𝐰,𝚯\displaystyle\underset{\mathbf{v},\boldsymbol{\Omega},\mathbf{R},\mathbf{u},\mathbf{w},\boldsymbol{\Theta}}{\mathrm{max.}}\,\,\, k𝒦Rk\displaystyle\sum\nolimits_{k\in\mathcal{K}}R_{k}\, (21)
s.t.\displaystyle\mathrm{s.t.}\,\,\,\,\,\, Rkln2lnwkwkek(𝐯,𝛀,uk)+1,k𝒦,\displaystyle R_{k}\ln 2\leq\ln w_{k}-w_{k}\,e_{k}(\mathbf{v},\boldsymbol{\Omega},u_{k})+1,\,k\in\mathcal{K},
lndet(𝚯i)+tr(𝚯i1(𝐃¯iH𝐕Σ𝐃¯i+𝛀¯i))Ni\displaystyle\ln\det\left(\boldsymbol{\Theta}_{i}\right)+\text{tr}\left(\boldsymbol{\Theta}_{i}^{-1}\left(\bar{\mathbf{D}}_{i}^{H}\mathbf{V}_{\Sigma}\bar{\mathbf{D}}_{i}+\bar{\boldsymbol{\Omega}}_{i}\right)\right)-Ni
j=1ilndet(𝛀j)irCFLln2,i,\displaystyle\,\,\,\,\,\,-\sum\nolimits_{j=1}^{i}\ln\det\left(\boldsymbol{\Omega}_{j}\right)\leq\frac{irC_{F}}{L}\ln 2,\,i\in\mathcal{L},
k𝒦𝐃iH𝐯k2+tr(𝛀i)Ptx,i,\displaystyle\sum\nolimits_{k\in\mathcal{K}}\left\|\mathbf{D}_{i}^{H}\mathbf{v}_{k}\right\|^{2}+\text{tr}\left(\boldsymbol{\Omega}_{i}\right)\leq P_{\text{tx}},\,i\in\mathcal{L},

where we have defined 𝐮={uk}k𝒦\mathbf{u}=\{u_{k}\}_{k\in\mathcal{K}}, 𝐰={wk}k𝒦\mathbf{w}=\{w_{k}\}_{k\in\mathcal{K}}, and 𝚯={𝚯i}i\boldsymbol{\Theta}=\{\boldsymbol{\Theta}_{i}\}_{i\in\mathcal{L}}.

We summarize an optimization process for solving (21) in Algorithm 1. Note that if the auxiliary variables 𝐮\mathbf{u}, 𝐰\mathbf{w}, and 𝚯\boldsymbol{\Theta} are fixed in the problem (21), the problem becomes a convex problem whose solution can be efficiently found with standard convex solvers such as, e.g., the CVX software [13]. Also, the optimal 𝐮\mathbf{u}, 𝐰\mathbf{w}, and 𝚯\boldsymbol{\Theta} for fixed primary variables 𝐯\mathbf{v} and 𝛀\boldsymbol{\Omega} are given in closed-forms as (17) and (20). Therefore, we can monotonically increase the objective sum-rate k𝒦Rk\sum_{k\in\mathcal{K}}R_{k} by alternately optimizing one of the variable sets {𝐯,𝛀}\{\mathbf{v},\boldsymbol{\Omega}\} and {𝐮,𝐰,𝚯}\{\mathbf{u},\mathbf{w},\boldsymbol{\Theta}\} for fixed other. Such an alternating optimization process is repeated until the convergence.

Algorithm 1 Alternating optimization algorithm for problem (14) under the WZ compression

1 Initialize: Set t1t\leftarrow 1, and initialize {𝐯,𝛀}\{\mathbf{v},\boldsymbol{\Omega}\} to an arbitrary point that satisfies (14b)-(14d);

2 repeat

3    tt+1t\leftarrow t+1;

4    update {𝐮,𝐰,𝚯}\{\mathbf{u},\mathbf{w},\boldsymbol{\Theta}\} according to (17) and (20) with {𝐯,𝚯}\{\mathbf{v},\boldsymbol{\Theta}\} being fixed;

5   update {𝐯,𝛀}\{\mathbf{v},\boldsymbol{\Omega}\} as a solution of the problem (21) with {𝐮,𝐰,𝚯}\{\mathbf{u},\mathbf{w},\boldsymbol{\Theta}\} being fixed;

6 until convergence

V-B P2P Compression

Next, we tackle (14) for the P2P compression. The access link rate constraint in (14b) can be readily addressed by the WMMSE formulation (15)-(17). Also, similar to the WZ compression case, it can be shown that (12) is fulfilled if there exists a positive definite matrix 𝚺iN×N\boldsymbol{\Sigma}_{i}\in\mathbb{C}^{N\times N} such that

log2det(𝚺i)+1ln2tr(𝚺i1(𝛀i+𝐃iH𝐕Σ𝐃i))Nln2\displaystyle\log_{2}\det\left(\boldsymbol{\Sigma}_{i}\right)+\frac{1}{\ln 2}\text{tr}\left(\boldsymbol{\Sigma}_{i}^{-1}\left(\boldsymbol{\Omega}_{i}+\mathbf{D}_{i}^{H}\mathbf{V}_{\Sigma}\mathbf{D}_{i}\right)\right)-\frac{N}{\ln 2}
log2det(𝛀i)rCFL,i.\displaystyle\,\,\,\,\,\,-\log_{2}\det\left(\boldsymbol{\Omega}_{i}\right)\leq\frac{rC_{F}}{L},\,i\in\mathcal{L}. (22)

Constraint (22) boils down to (14c) with the auxiliary matrix 𝚺i\boldsymbol{\Sigma}_{i} set to

𝚺i=𝛀i+𝐃iH𝐕Σ𝐃i.\displaystyle\boldsymbol{\Sigma}_{i}=\boldsymbol{\Omega}_{i}+\mathbf{D}_{i}^{H}\mathbf{V}_{\Sigma}\mathbf{D}_{i}. (23)

As a consequence, an equivalent formulation of (14) for the P2P compression scheme is given as

max.𝐯,𝛀,𝐑,𝐮,𝐰,𝚺\displaystyle\underset{\mathbf{v},\boldsymbol{\Omega},\mathbf{R},\mathbf{u},\mathbf{w},\boldsymbol{\Sigma}}{\mathrm{max.}}\,\,\, k𝒦Rk\displaystyle\sum\nolimits_{k\in\mathcal{K}}R_{k}\, (24)
s.t.\displaystyle\mathrm{s.t.}\,\,\,\,\,\, Rkln2lnwkwkek(𝐯,𝛀,uk)+1,k𝒦,\displaystyle R_{k}\ln 2\leq\ln w_{k}-w_{k}\,e_{k}(\mathbf{v},\boldsymbol{\Omega},u_{k})+1,\,k\in\mathcal{K},
lndet(𝚺i)+tr(𝚺i1(𝛀i+𝐃iH𝐕Σ𝐃i))N\displaystyle\ln\det\left(\boldsymbol{\Sigma}_{i}\right)+\text{tr}\left(\boldsymbol{\Sigma}_{i}^{-1}\left(\boldsymbol{\Omega}_{i}+\mathbf{D}_{i}^{H}\mathbf{V}_{\Sigma}\mathbf{D}_{i}\right)\right)-N
lndet(𝛀i)rCFLln2,i,\displaystyle\,\,\,\,\,\,\,\,\,\,\,\,\,-\ln\det\left(\boldsymbol{\Omega}_{i}\right)\leq\frac{rC_{F}}{L}\ln 2,\,i\in\mathcal{L},
k𝒦𝐃iH𝐯k2+tr(𝛀i)Ptx,i,\displaystyle\sum\nolimits_{k\in\mathcal{K}}\left\|\mathbf{D}_{i}^{H}\mathbf{v}_{k}\right\|^{2}+\text{tr}\left(\boldsymbol{\Omega}_{i}\right)\leq P_{\text{tx}},\,i\in\mathcal{L},

where 𝚺={𝚺i}i\mathbf{\Sigma}=\{\mathbf{\Sigma}_{i}\}_{i\in\mathcal{L}}. A stationary point of problem (24) can be efficiently found through an alternating optimization procedure between two blocks of variables {𝐯,𝛀}\{\mathbf{v},\boldsymbol{\Omega}\} and {𝐮,𝐰,𝚺}\{\mathbf{u},\mathbf{w},\boldsymbol{\Sigma}\} as detailed in Algorithm 2.

Algorithm 2 Alternating optimization algorithm for problem (24) under the P2P compression

1 Initialize: Set t1t\leftarrow 1, and initialize {𝐯,𝛀}\{\mathbf{v},\boldsymbol{\Omega}\} to an arbitrary point that satisfies (14b)-(14d);

2 repeat

3    tt+1t\leftarrow t+1;

4    update {𝐮,𝐰,𝚺}\{\mathbf{u},\mathbf{w},\boldsymbol{\Sigma}\} according to (17) and (23);

5   update {𝐯,𝛀}\{\mathbf{v},\boldsymbol{\Omega}\} as a solution of the problem (24) with {𝐮,𝐰,𝚺}\{\mathbf{u},\mathbf{w},\boldsymbol{\Sigma}\} being fixed;

6 until convergence

VI Numerical Results

This section presents numerical results assessing the proposed algorithms. The UEs are randomly located within a circular area of radius 200 m, while the APs are positioned evenly along the boundary circle, maintaining equal angular intervals. For a given distance dk,id_{k,i} between UE kk and AP ii, the path-loss ρk,i\rho_{k,i} of wireless access link between UE kk and AP ii is given as ρk,i=ρ0(dk,i/d0)α\rho_{k,i}=\rho_{0}(d_{k,i}/d_{0})^{-\alpha} where dk,id_{k,i} accounts for the distance of the associated link, ρ0=0\rho_{0}=0 dB is the reference path-loss at d0=30d_{0}=30 m, and α=3\alpha=3 indicates the path-loss exponent. Unless stated otherwise, the CSI error power βk,i\beta_{k,i} is assumed to be 0.1ρk,i0.1\rho_{k,i} which corresponds to the case where the error variance equals 10 %\% of the actual channel power. We also assume that the fronthaul and access transmission phases span an equal time duration, i.e., nA=nFn_{A}=n_{F} and r=1r=1.

Refer to caption
Figure 3: Average sum-rate versus the SNR Ptx/σz2P_{\text{tx}}/\sigma_{z}^{2} (K=10K=10 UEs, L{3,5}L\in\{3,5\} APs each with N=2N=2 antennas, and CF=1C_{F}=1 bit/symbol)

In Fig. 3, we plot the average sum-rate with respect to the transmit SNR Ptx/σz2P_{\text{tx}}/\sigma_{z}^{2} with K=10K=10 UEs, L{3,5}L\in\{3,5\} APs each with N=2N=2 antennas, and CF=1C_{F}=1 bit/symbol. The sum-rate performance of both compression schemes increases with the SNR while saturating to finite levels due the finite-capacity fronthaul links. Also, the WZ compression scheme, in which the APs utilize the side information for decompression, provides notable gains compared to the P2P compression scheme. It should be emphasized that increasing the number LL of APs has conflicting impacts on the sum-rate. This is because, although deploying more APs brings a beamforming gain, the effective compression rate available for each sample will decrease (see the right-hand sides of (12) and (11)), since the fronthaul phase is divided into LL time slots, where each slot is used for delivering a single fronthaul block. For this reason, the sum-rate of the P2P compression remains almost same when the number of APs increases from L=3L=3 to 55. Instead, the performance of the WZ scheme is significantly improved with more APs, because the degraded quantization fidelity is overcome by utilizing the side information for decompression. We validate this pattern once again in Fig. 4 by plotting the average sum-rate with respect to the number of APs LL for a cell-free MIMO system with K=10K=10 UEs, N=2N=2 AP antennas, CF=2C_{F}=2 bit/symbol, and Ptx/σz2=15P_{\text{tx}}/\sigma_{z}^{2}=15 dB. Similar to Fig. 3, we observe that as the number of APs LL increases, only the sum-rate of the WZ scheme shows improvement.

Refer to caption
Figure 4: Average sum-rate versus the number of APs LL (K=10K=10 UEs, N=2N=2 AP antennas, CF=2C_{F}=2 bit/symbol, and Ptx/σz2=15P_{\text{tx}}/\sigma_{z}^{2}=15 dB)
Refer to caption
Figure 5: Average sum-rate versus the fronthaul capacity CFC_{F} (K{5,15}K\in\{5,15\} UEs, L=4L=4 APs each with N=2N=2 antennas, and Ptx/σz2=10P_{\text{tx}}/\sigma_{z}^{2}=10 dB)

Fig. 5 shows the impact of the fronthaul capacity CFC_{F} on the average sum-rate for the downlink of a cell-free MIMO system with K{5,15}K\in\{5,15\} UEs, L=4L=4 APs each with N=2N=2 antennas, and Ptx/σz2=10P_{\text{tx}}/\sigma_{z}^{2}=10 dB. The sum-rates of both the P2P and WZ schemes grow with the fronthaul capacity CFC_{F}, since with a larger CFC_{F}, the fronthaul quantization causes a smaller distortion. The performance gain of the WZ scheme is most pronounced in the intermediate regime of CFC_{F}, which means that the advantage of the advanced WZ compression is minor when CFC_{F} is sufficiently large or small.

Refer to caption
Figure 6: Average sum-rate versus the normalized CSI error variance βk,i/ρk,i\beta_{k,i}/\rho_{k,i} (K=10K=10 UEs, L=4L=4 APs each with N=2N=2 antennas, CF{2,4}C_{F}\in\{2,4\} bit/symbol and Ptx/σz2=20P_{\text{tx}}/\sigma_{z}^{2}=20 dB)

Fig. 6 assesses the significance of the proposed robust design, which incorporates the CSI error in the optimization process, by comparing the average sum-rate of the proposed WZ scheme with that of a non-robust design. The non-robust design optimizes the precoding vector 𝐯\mathbf{v} and the fronthaul compression strategy 𝛀\boldsymbol{\Omega} without considering CSI error, but its average sum-rate is evaluated using the true CSI error variance. We plot the average sum-rate performance of the proposed robust and the baseline non-robust schemes while gradually increasing the normalized CSI error variance, defined as βk,i/ρk,i[0,1]\beta_{k,i}/\rho_{k,i}\in[0,1], for a cell-free MIMO system with K=10K=10 UEs, L=4L=4 APs each with N=2N=2 antennas, CF{2,4}C_{F}\in\{2,4\} bit/symbol and Ptx/σz2=20P_{\text{tx}}/\sigma_{z}^{2}=20 dB. With the increased CSI errors, the performance of the non-robust scheme experiences a significant degradation. Nevertheless, by incorporating CSI imperfections into the design of precoding and fronthaul compression strategies, it is feasible to mitigate the performance degradation. We also observe that the performance gap between the robust and non-robust schemes becomes more evident as CFC_{F} increases.

VII Conclusion

We have studied the downlink of a cell-free MIMO system with radio stripes fronthaul network. Noting the serial transfer of fronthaul-compressed block through APs, an application of WZ compression was discussed whereby each AP decompresses all the received blocks, that pass through it, and utilize them for decompressing its desired precoded signal. For both the standard P2P and advanced WZ compression, we tackled the sum-rate maximization problems under per-AP power and fronthaul capacity constraints. Via numerical results, the advantages of the WZ compression compared to the P2P compression were validated. As future work, we consider an extension to the segmented fronthaul network in which multiple, not a single, radio stripes are connected in parallel [14]. Furthermore, an intriguing and demanding subject would involve the development of a deep learning-based sequential linear processing of APs that ensures scalability with respect to the numbers of UEs and APs [15].

References

  • [1] H. Q. Ngo, A. Ashikhmin, H. Yang, E. G. Larsson, and T. L. Marzetta, “Cell-free massive MIMO versus small cells,” IEEE Trans. Wireless Commun., vol. 16, no. 3, pp. 1834–1850, Mar. 2017.
  • [2] E. Björnson and L. Sanguinetti, “Making cell-free massive MIMO competitive with MMSE processing and centralized implementation,” IEEE Trans. Wireless Commun., vol. 19, no. 1, pp. 77–90, Jan. 2020.
  • [3] Z. H. Shaik, E. Björnson and E. G. Larsson, “MMSE-optimal sequential processing for cell-free massive MIMO with radio stripes,” IEEE Trans. Commun., vol. 69, no. 11, pp. 7775–7789, Nov. 2021.
  • [4] I. Chiotis and A. L. Moustakas, “On the uplink performance of finite-capacity radio stripes,” in Proc. 2022 IEEE International Mediterranean Conference on Communications and Networking (MeditCom), pp. 118-123, Athens, Greece, Sept. 2022.
  • [5] O. L. A. Lopez, and et al., “Massive MIMO with radio stripes for indoor wireless energy transfer,” IEEE Trans. Wireless Commun., vol. 21, no. 9, pp. 7088–7104, Sept. 2022.
  • [6] P. Zhang and F. M. J. Willems, “On the downlink capacity of cell-free massive MIMO with constrained fronthaul capacity,” Entropy, vol. 22, no. 4, pp. 1–31, Arp. 2020.
  • [7] Z. Xiong, A. Liveris, and S. Cheng, “Distributed source coding for sensor networks,” IEEE Signal Process. Mag., vol. 21, no. 5, pp. 80–94, Sept. 2004.
  • [8] S.-H. Park, O. Simeone, O. Sahin and S. Shamai (Shitz), “Fronthaul compression for cloud radio access networks: Signal processing advances inspired by network information theory,” IEEE Signal Process. Mag., vol. 31, no. 6, pp. 69–79, Nov. 2014.
  • [9] A. E. Gamal and Y.-H. Kim, Network Information Theory, Cambridge, U.K.: Cambridge University Press, pp. 1–685, 2011.
  • [10] S. S. Christensen, R. Agarwal, E. De Carvalho and J. M. Cioffi, “Weighted sum-rate maximization using weighted MMSE for MIMO-BC beamforming design,” IEEE Trans. Wireless Commun., vol. 7, no. 12, pp. 4792–4799, Dec. 2008.
  • [11] Y. Zhou and W. Yu, “Fronthaul compression and transmit beamforming optimization for multi-antenna uplink C-RAN,” IEEE Trans. Signal Process., vol. 64, no. 16, pp. 4138–4151, Aug. 2016.
  • [12] J. Choi, N. Lee, S.-N. Hong, and G. Caire, “Joint user selection, power allocation, and precoding design with imperfect CSIT for multi-cell MU-MIMO downlink systems,” IEEE Trans. Wireless Commun., vol. 19, no. 1, pp. 162–176, Jan. 2020.
  • [13] M. Grant and S.Boyd, “CVX: MATLAB software for disciplined convex programming,” Second Edition, Third Printing, pp. 1–786, Ver. 2.2, Jan. 2020.
  • [14] A. L. P. Fernandes, and et al., “Cell-free massive MIMO with segmented fronthaul: Reliability and protection aspects,” IEEE Wireless Commun. Lett., vol. 11, no. 8, pp. 1580–1584, Aug. 2022.
  • [15] H. Lee, J. Kim and S.-H. Park, “Learning optimal fronthauling and decentralized edge computation in fog radio access networks,” IEEE Trans. Wireless Commun., vol. 20, no. 9, pp. 5599–5612, Sep. 2021.