Open AccessArticle

Satellite-Assisted Cell-Free Massive MIMO Systems with Multi-Group Multicast

Jiamin Li

^1,†,

Lingling Chen

²,

Pengcheng Zhu

^1,*,

Dongming Wang

¹ and

Xiaohu You

National Mobile Communications Research Laboratory, Southeast University, Nanjing 210096, China

School of Information Science and Engineering, Southeast University, Nanjing 210096, China

Author to whom correspondence should be addressed.

^†

Current address: Wireless Valley, Nanjing 210096, China.

Sensors 2021, 21(18), 6222; https://doi.org/10.3390/s21186222

Submission received: 2 August 2021 / Revised: 13 September 2021 / Accepted: 15 September 2021 / Published: 16 September 2021

(This article belongs to the Special Issue Cell-Free Ultra Massive MIMO in 6G and Beyond Networks)

Download

Browse Figures

Versions Notes

Abstract

In this paper, we use satellite-assisted and multi-group multicast mechanisms to relieve ground traffic pressure and improve data transmission efficiency of cell-free massive MIMO systems. We propose to estimate channel state information (CSI) by common pilot scheme. Given the estimated CSI, we derive the closed-form expressions of achievable rate with maximum ratio transmission (MRT) and zero-forcing (ZF) precoding. The correctness of the closed-form expressions is verified through simulations. The results show that with the help of satellite and multicast, the average system spectrum efficiency (SE) can be significantly improved.

Keywords:

satellite-assisted; cell-free massive MIMO; multi-group multicast

1. Introduction

With the rapid development of mobile communication, more and more people and machines are connected to the network, and the exponential growth of data traffic has brought tremendous pressure to the terrestrial cellular network. This promotes the integration of the satellite network and terrestrial cellular network. Using satellites to assist terrestrial base stations, the data traffic of the ground cellular networks is shunted to provide users with better services. Moreover, satellite network can cope with the situation where limited local resources affect service quality by virtue of its global awareness. Therefore, the integration of satellite network and terrestrial network can effectively relieve ground traffic pressure and promote global information dissemination [1,2].

Recently, the concept of cell-free massive multiple-input multiple-output (MIMO) has attracted widespread attention in the academic community [3,4,5]. In cell-free massive MIMO systems, multiple access points (APs) jointly serve users within the same time-frequency resource, without existing cell. Ideally, all APs are connected to a central processing unit (CPU) to exchange information. The estimated channel state information (CSI) obtained by uplink pilot is used to precode the transmitted data in the downlink and perform data detection in the uplink. Moreover, the concept of cell-free massive MIMO is to combine massive MIMO [6] with distributed MIMO [7]. Therefore, it has the advantages of both massive MIMO and distributed MIMO [8].

Faced with the content-centric form of communication, hot information needs to be transmitted to multiple users, resulting that repeated transmission of information will bring a waste of resources [9]. Point-to-multipoint transmission using multicast can provide services to more users under the same resource conditions, therefore improving efficiency. In multicast scenario, users needing the same information are grouped into the same group, then multiple independent data streams are transmitted to different groups. Each group is called a multicasting group. Through multicast transmission, we can use the CSI at the transmitter to optimize the transmission based on a desired performance metric [10]. Ref. [11] introduced mixed unicast and multicast transmission under centralized massive MIMO systems, and performed multi-objective optimization to improve system performance. In [12], a non-orthogonal transmission framework based on layered-division multiplexing (LDM) is proposed to support both multicast and unicast services in cooperative multi-cell cellular networks. The above articles are all aimed at traditional centralized massive MIMO multicast systems. Considering the scenario of multi-group multicast cell-free massive MIMO, ref. [13] compared and analyzed the spectral efficiency (SE) after channel estimation and it is found that the use of downlink pilot can enhance the user’s ability to detect signals. A conjugate beamforming precoder was used in [14] to analyze the performance of multi-group multicast cell-free massive MIMO systems. They all consider only terrestrial cell-free massive MIMO systems with estimated CSI, not taking satellite-assisted into consideration.

Due to the growth of data traffic and the need to expand signal coverage, there have been many studies on the performance of integrated satellite ground multicast systems in recent years. Ref. [15] studied user access and resource allocation in the integrated satellite and terrestrial multi-cell multicast scenario, and formulated an optimal resource schedule scheme. In [16], a time-division multi-group multicast transmission scheme was proposed in the scenario of integrating satellites and terrestrial multi-cells, dividing the signal transmission into ground transmission phase and satellite transmission phase, which effectively improved system SE. The above research on satellite communication is based on ideal CSI which is not practical and will inevitably affect the performance of the system. Therefore, in this article, we will consider the non-ideal CSI and reasonably analyze the SE of the satellite-assisted terrestrial multicast systems.

In this paper, we consider satellite-assisted cell-free massive MIMO multi-group multicast scenarios. The users who need the same information are divided into the same group. APs and satellite multiplex the same bandwidth to provide multicast services for each group of users on the ground, and they are all equipped with multiple antennas. Assuming users in each multicast group receive the signals sent from the satellite and APs at the same time, then we use the estimated CSI to analyze the SE of the system. The main contributions of this paper are summarized as follows:

We consider a multi-group multicast scheme in satellite-assisted cell-free massive MIMO systems. We use common pilot scheme to estimate satellite and terrestrial channel.
With the estimated CSI, the closed-form expressions of the SE of multicast users with maximum ratio transmission (MRT) and zero-forcing (ZF) precoding schemes in satellite-assisted cell-free massive MIMO are derived. The SE of satellite-assisted cell-free massive MIMO multicast systems is analyzed.
The accuracy of the derived closed-form expressions and the effectiveness of the multi-group multicast transmission scheme in satellite-assisted cell-free massive MIMO systems are verified. Insightful conclusion can be drawn that satellite-assisted can relieve terrestrial pressure and multicast can improve system SE effectively.

The remainder of the paper is organized as follows. In Section 2, we describe the system model including system configuration, satellite channel model, terrestrial channel model and signal transmission model. Section 3 consider the non-ideal CSI, using minimum mean-square error (MMSE) method to estimate satellite channel and terrestrial channel and the closed-form expressions of the MRT and ZF downlink achievable rate are derived. Representative numerical results are given in Section 4 before we conclude the paper in Section 5.

The following notations are used. All boldface letters stand for vectors (lower case) or matrices (upper case).

I_{N}

is the size-N identity matrix. Italic letters (e.g., X or x) denote scalars. The transpose, Hermitian transpose and trace operators are denoted by

{(\cdot)}^{T}

{(\cdot)}^{H}

and

tr (\cdot)

, respectively.

C^{m \times n}

denotes the set of

m \times n

complex valued matrices.

| x |

is the absolute value of a scalar x,

∥ X ∥

is the spectral norm of a matrix

X

x \sim CN (0, σ^{2})

means that x is a circularly symmetric complex Gaussian random variable with mean zero and variance

σ^{2}

E [\cdot]

denotes the expectation operator.

2. System Model

2.1. System Configuration

We consider a satellite-assisted cell-free massive MIMO system where

N_{B}

APs and a LEO multi-beam satellite provide multicast services to K single-antenna users in a cooperative manner. Assuming that the satellite and all APs use the same frequency spectrum during the entire transmission process, and all APs are distributed within the coverage area of the satellite. Each AP has

M_{B}

antennas and the satellite has

M_{S}

antennas. All users are grouped into G groups according to the requested contents. In multicast group g, the number of users is

K_{g}

. These users are arbitrarily distributed in the system. As shown in Figure 1, there are 34 users in total, which are divided into 3 multicast groups, then,

G = 3, K_{1} = 11, K_{2} = 11, K_{3} = 12

2.2. Channel Model

For the satellite channel, we assume that the satellite link is quasi-static. We do not take the influence of rain fading into consideration. The channel vector between the satellite and the k-th user in the g-th group is [17]

g_{g, k} = L_{g, k}^{\frac{1}{2}} b_{g, k}, g = 1, 2, \dots, G; k = 1, 2, \dots, K_{g},

(1)

where

L_{g, k}

is large-scale attenuation and can be defined by

L_{g, k} = \frac{G_{t} G_{r} B_{g, k}}{P L_{g, k}},

(2)

G_{t}

and

G_{r}

represent satellite and user antenna gain, respectively.

P L_{g, k}

is free space loss in dB [17],

P L_{g, k} = 90.45 + 20 \lg (f_{c}) + 20 \lg (d_{g, k}),

(3)

f_{c}

is satellite carrier frequency,

d_{g, k}

is the direct distance between the satellite and the user on the ground,

A_{g, k}

is beam gain,

B_{g, k} = {|\frac{J_{1} (u_{g, k})}{2 u_{g, k}} + 36 \frac{J_{3} (u_{g, k})}{u_{g, k}^{3}}|}^{2}

, where

u_{g, k} = 2.07123 \frac{sin φ_{g, k}}{sin φ_{3 d B}}

φ_{g, k}

indicates the angle between the k-th user in the g group and the center of the corresponding beam,

φ_{3 d B}

represents the angle corresponding to the 3dB power loss of each spot beam,

J_{1}

and

J_{3}

are the first-order and third-order Bessel functions of the first kind, respectively.

b_{g, k} = {[b_{g, k, 1}, b_{g, k, 2}, \dots, b_{g, k, M_{S}}]}^{T}

is the small-scale attenuation between the satellite and the user, using the Shadowed Rice (SR) model [18],

b_{g, k, i} = A_{g, k, i} e^{j ψ_{g, k, i}} + Z_{g, k, i} e^{j ϕ_{g, k, i}}

, where A and Z are the amplitudes of the scattering and line of sight (LoS) parts respectively, which obey the Rayleigh distribution and Nakagami-m distribution respectively. Assuming that

ψ

is a static random phase, uniformly distributed in

[0, 2 π)

ϕ

is the determined phase of the LoS part, which is independent of A and Z. Therefore, the SR distribution finally can be determined by the three parameters

(b, m, Ω)

. The parameter values corresponding to different shadow levels can be referred to [18].

Ω = E [Z^{2}]

indicates the average power of the LoS part,

2 b = E [Z^{2}]

is the average power of the scattering part.

m = E {[Z^{2}]}^{2} / Var \{Z^{2}\} \geq 0

represents the Nakagami parameter determined according to the severity of the attenuation. Then, the probability density function of

{|b_{g, k, i}|}^{2}

can be expressed as

f_{{|b_{g, k, i}|}^{2}} (x) = {(\frac{2 b m}{2 b m + Ω})}^{m} \frac{x}{b} e^{(- \frac{x^{2}}{2 b})} B (b, m, Ω, x),

(4)

where

x \geq 0

B (b, m, Ω, x) =_{1} F_{1} (m, 1, \frac{Ω x^{2}}{2 b (2 b m + Ω)})

is Kummer function. Let

η_{g, k}

be the statistical character of satellite channel.

It is worth noting that the channel model we consider here is a little different from the actual satellite channel. First, the satellite channel is actually constantly changing, and the influence of rain attenuation is also difficult to avoid, especially in the Ka-band. It becomes difficult for us to obtain real-time channel state information. Additionally, it may be that the currently acquired and estimated channel state information is no longer the current channel state. The channel aging problem also needs to be solved urgently. Secondly, if low-orbit satellites are considered, the Doppler frequency shift caused by the high-speed movement of the satellites actually needs to be taken into account in the channel model. Doppler frequency shift not only causes inter-symbol interference between signals. Doppler frequency shift spreading will also cause rapid time change of the channel, making it difficult for the receiving terminal to obtain instantaneous channel state information. Therefore, the estimation range of the Doppler frequency offset and the rate of change of the Doppler frequency offset during the integration process of the satellite system and the terrestrial 5G system need to be emphatically considered and urgently needed to be resolved.

For the terrestrial channel, we apply uncorrelated Rayleigh fading channel. The channel vector between APs and the k-th user in the g group can be expressed as

h_{g, k} = β_{g, k}^{\frac{1}{2}} u_{g, k}, g = 1, 2, \dots, G; k = 1, 2, \dots, K_{g},

(5)

where

β_{g, k} = diag (β_{1, g, k}, β_{2, g, k}, \dots, β_{N_{B}, g, k}) \otimes I_{M_{B}}

β_{n, g, k} = d_{n, g, k}^{- a}

represents the large-scale attenuation between the n-th AP and the k-th user in group g, a is the path loss factor, generally between 3.0 and 5.0,

d_{n, g, k}

represents the distance between the n-th AP and the k-th user in group g,

u_{g, k}

is uncorrelated small-scale fast fading vector,

u_{g, k} \sim CN (0, I_{N_{B} M_{B}})

, including independent identically distributed cyclic symmetric complex Gaussian (ZMCSCG) random variables with a mean of 0 and a variance of 1.

3. Performance Analysis

3.1. Channel Estimation

Obtaining accurate channel state information is particularly important in wireless communication transmission. In an actual communication scenario, neither the transmitter nor the receiver may be able to fully understand the channel state. At this time, it is necessary to estimate the channel state information. The channel state information on the receiver side can be obtained through pilot transmission. The channel state information on the transmitter side can be obtained through the feedback link or channel reciprocity. According to [19], the estimation accuracy of the channel state information is inconsistent, resulting in inconsistent system capacity. It can be seen that when the estimation error increases, the system capacity decreases. It can be seen that accurate estimation of channel state information is particularly important.

In this paper, we use pilot-assisted channel estimation. Due to the channel reciprocity, the CSI obtained by uplink channel estimation can be used for the precoding of downlink data transmission. In multicast transmission, each multicast group has many users. Take a TV channel for multicast transmission as an example, when hundreds of users watch the same channel, the resources in the coherent range will be exhausted. Therefore, for multicast users, it is more appropriate to adopt the new common pilot allocation scheme mentioned in [20]. Users in each multicast group share the same pilot, so a total of G pilots are required instead of

\sum_{g = 1}^{G} K_{g}

pilots. The length of pilot sequence is

τ_{p}

For the satellite channel, the pilot sequence transmitted by the g-th user group is

ψ_{g}

ψ_{g} \in C^{τ_{p} \times 1}

, and the pilot signal received by the satellite is expressed as

r_{p} = \sqrt{τ_{p} ρ_{p, S}} \sum_{g = 1}^{G} \sum_{k = 1}^{K_{g}} g_{g, k} ψ_{g}^{H} + N_{p, S} .

(6)

The MMSE estimate of the satellite channel is [21]

{\hat{g}}_{g, k} = \frac{\sqrt{τ_{p} ρ_{p, S}} η_{g, k}}{τ_{p} ρ_{p, S} \sum_{k^{'} = 1}^{K_{g}} η_{g, k^{'}} + 1} (\sqrt{τ_{p} ρ_{p, S}} \sum_{k = 1}^{K_{g}} g_{g, k} + n_{S, g}) = γ_{g, k}^{1 / 2} {\hat{b}}_{g, k},

(7)

where

n_{S, g} \sim CN (0, I_{M_{S}})

is standardized additive noise,

ρ_{p, S}

is the power of uplink pilot,

γ_{g, k} = \frac{τ_{p} ρ_{p, S} η_{g, k}^{2}}{τ_{p} ρ_{p, S} \sum_{k^{'} = 1}^{K_{g}} η_{g, k^{'}} + 1}

is the equivalent large-scale fading and

{\hat{b}}_{g, k} \sim CN (0, I_{M_{S}})

Let

g_{g} = \sum_{t = 1}^{K_{g}} \sqrt{τ_{p} ρ_{p, S}} g_{g, t}

represent the linear superposition of the user channels used in group g, then the estimated value can be expressed as

{\hat{g}}_{g} = \frac{\sum_{t = 1}^{K_{g}} τ_{p} ρ_{p, S} η_{g, t}}{τ_{p} ρ_{p, S} \sum_{k^{'} = 1}^{K_{g}} η_{g, k^{'}} + 1} (\sqrt{τ_{p} ρ_{p, S}} \sum_{k = 1}^{K_{g}} g_{g, k} + n_{S, g}) = ε_{g}^{1 / 2} {\hat{b}}_{g},

(8)

where

ε_{g} = \frac{{(\sum_{t = 1}^{K_{g}} τ_{p} ρ_{p, S} η_{g, t})}^{2}}{τ_{p} ρ_{p, S} \sum_{k^{'} = 1}^{K_{g}} η_{g, k^{'}} + 1}

is the equivalent large-scale fading and

{\hat{b}}_{g} \sim CN (0, I_{M_{S}})

Using a

M_{S} \times G

matrix

\hat{G} = [{\hat{g}}_{1}, {\hat{g}}_{2}, \dots, {\hat{g}}_{G}]

to represent the channel estimation vector between the satellite and users in all G multicast groups, as a result, there is only one scalar coefficient between

{\hat{g}}_{g, k}

and

{\hat{g}}_{g}

\begin{matrix} {\hat{g}}_{g, k} & = c_{g, k} {\hat{g}}_{g} \end{matrix}

(9)

\begin{matrix} c_{g, k} & = \frac{\sqrt{τ_{p} ρ_{p, S}} η_{g, k}}{\sum_{t = 1}^{K_{g}} τ_{p} ρ_{p, S} η_{g, t}} . \end{matrix}

(10)

Similarly, the MMSE estimate of the terrestrial channel is

{\hat{h}}_{g, k} = χ_{g, k}^{1 / 2} {\hat{u}}_{g, k},

(11)

where

χ_{g, k} = diag (χ_{1, g, k}, χ_{2, g, k}, \dots, χ_{N_{B}, g, k}) \otimes I_{M_{B}}

is the equivalent large-scale fading,

χ_{n, g, k} = \frac{τ_{p} ρ_{p, B} β_{n, g, k}}{τ_{p} ρ_{p, B} \sum_{k^{'} = 1}^{K_{g}} β_{n, g, k^{'}} + 1}

and

{\hat{u}}_{g, k} \sim CN (0, I_{N_{B} M_{B}})

The estimated value of

h_{g}

can be expressed as

{\hat{h}}_{g} = α_{g}^{1 / 2} {\hat{u}}_{g},

(12)

where

α_{g} = diag (α_{1, g}, α_{2, g}, \dots, α_{N_{B}, g}) \otimes I_{M_{B}}

is the equivalent large-scale fading,

α_{n, g} = \frac{{(\sum_{t = 1}^{K_{g}} τ_{p} ρ_{p, B} β_{n, g, t})}^{2}}{\sum_{k^{'} = 1}^{K_{g}} τ_{p} ρ_{p, B} β_{n, g, k^{'}} + 1}

and

{\hat{u}}_{g} \sim CN (0, I_{N_{B} M_{B}})

. Using a

N_{B} M_{B} \times G

matrix

\hat{H} = [{\hat{h}}_{1}, {\hat{h}}_{2}, \dots, {\hat{h}}_{G}]

to represent the channel estimation vector between APs and users in all G multicast groups, as a result, there is a diagonal matrix difference between

{\hat{h}}_{g, k}

and

{\hat{h}}_{g}

{\hat{h}}_{g, k} = κ_{g, k} {\hat{h}}_{g},

(13)

where

κ_{g, k} = diag (κ_{1, g, k}, κ_{2, g, k}, \dots, κ_{N_{B}, g, k}) \otimes I_{M_{B}}

and

κ_{n, g, k} = \frac{\sqrt{τ_{p} ρ_{p, B}} β_{n, g, k}}{\sum_{t = 1}^{K_{g}} τ_{p} ρ_{p, B} β_{n, g, t}}

3.2. Downlink Spectral Efficiency Analysis

In satellite-assisted cell-free massive MIMO multi-group multicast systems, APs and the satellite provide multicast services for G group users. We assume that users in each multicast group receive the signals sent from the satellite and APs at the same time [16]. For downlink signal transmission, the multicast signal received by the k-th user in group g is expressed as

y_{g, k} = h_{g, k}^{H} w_{g} s_{g} + g_{g, k}^{H} v_{g} s_{g} + h_{g, k}^{H} \sum_{i = 1, i \neq g}^{G} w_{i} s_{i} + g_{g, k}^{H} \sum_{i = 1, i \neq g}^{G} v_{i} s_{i} + n_{g, k},

(14)

where

s_{g} \in C

indicates the data symbol sent to users in group g, and

E [{|s_{g}|}^{2}] = 1

n_{g, k} \sim CN (0, σ_{d l}^{2})

is additive noise,

w_{g} \in C^{N_{B} M_{B} \times 1}

and

v_{g} \in C^{M_{S} \times 1}

represent the precoding vectors sent by APs and satellite to group g, respectively.

It is worth noting that we do not take the interference between satellite systems and terrestrial systems into consideration. However, during the integration process of satellite systems and terrestrial systems, the inter-system interference caused by frequency sharing is also worth considering.

Assuming that users only use statistical channel state information for signal decoding, according to [22], the received signal can be rewritten as

y_{g, k} = E [h_{g, k}^{H} w_{g} + g_{g, k}^{H} v_{g}] + n_{g, k}^{'},

(15)

where

n_{g, k}^{'} = [h_{g, k}^{H} w_{g} + g_{g, k}^{H} v_{g} - E [h_{g, k}^{H} w_{g} + g_{g, k}^{H} v_{g}]] s_{g} + h_{g, k}^{H} \sum_{i = 1, i \neq g}^{G} w_{i} s_{i} + g_{g, k}^{H} \sum_{i = 1, i \neq g}^{G} v_{i} s_{i} + n_{g, k} .

(16)

E [h_{g, k}^{H} w_{g} + g_{g, k}^{H} v_{g}] s_{g}

is the only expected signal of the k-th user in group g,

n_{g . k}^{'}

is interference, then the achievable SE of the k-th user in group g can be expressed as

R_{g, k} = {log}_{2} (1 + \frac{{|E [h_{g, k}^{H} w_{g} + g_{g, k}^{H} v_{g}]|}^{2}}{\sum_{i = 1}^{G} E [{|h_{g, k}^{H} w_{i} + g_{g, k}^{H} v_{i}|}^{2}] - {|E [h_{g, k}^{H} w_{g} + g_{g, k}^{H} v_{g}]|}^{2} + σ_{d l}^{2}}) .

(17)

In this paper, we consider MRT and ZF precoding. For MRT precoding, the precoding vector of the terrestrial base station for the g-th multicast group is

w_{g}^{mrt} = p_{B, g}^{dl} {\hat{h}}_{g},

(18)

where

p_{B, g}^{dl} = 1 / \sqrt{E [∥ {\hat{h}}_{g} ∥^{2}]} = 1 / \sqrt{M_{B} \sum_{m = 1}^{N_{B}} α_{m, g}}

is power normalization parameter, and

E [∥ w_{g}^{mrt} ∥^{2}] = 1

. The precoding vector of the satellite for the g-th multicast group is

v_{g}^{mrt} = p_{S, g}^{dl} {\hat{g}}_{g},

(19)

where

p_{S, g}^{dl} = 1 / \sqrt{E [∥ {\hat{g}}_{g} ∥^{2}]} = 1 / \sqrt{M_{S} ε_{g}}

is power normalization parameter, and

E [∥ v_{g}^{mrt} ∥^{2}] = 1

Theorem 1.

When the base station and satellite terminals use MRT precoding, the closed-form expression of the downlink achievable rate can be given by

R_{g, k}^{mrt} = {log}_{2} (1 + S I N R_{g, k}^{m r t})

S I N R_{g, k}^{m r t} = \frac{{(p_{B, g}^{dl} M_{B} \sum_{m = 1}^{N_{B}} κ_{m, g, k} α_{m, g} + p_{S, g}^{dl} M_{S} ε_{g} c_{g, k})}^{2}}{\sum_{i = 1}^{G} [\sum_{m = 1}^{N_{B}} (β_{m, g, k} - χ_{m, g, k}) + η_{g, k} - γ_{g, k}] - 2 p_{B, g}^{dl} p_{S, g}^{dl} M_{B} M_{S} ε_{g} c_{g, k} \sum_{m = 1}^{N_{B}} κ_{m, g, k} α_{m, g} + σ_{d l}^{2}} .

(20)

Proof of Theorem 1.

The proof is given in Appendix A. □

For ZF precoding, the precoding vector of the terrestrial base station for the g-th multicast group is

w_{g}^{ZF} = q_{B, g}^{dl} a_{B, g},

(21)

where

a_{B, g}

is the g-th column of

\hat{H} {({\hat{H}}^{H} \hat{H})}^{- 1}

q_{B, g}^{dl} = 1 / \sqrt{E [∥ a_{B, g} ∥^{2}]} = 1 / \sqrt{(M_{B} - G) \sum_{m = 1}^{N_{B}} α_{m, g}}

is power normalization parameter, and

E [∥ w_{g}^{ZF} ∥^{2}] = 1

. The precoding vector of the satellite for the g-th multicast group is

v_{g}^{ZF} = q_{S, g}^{dl} a_{S, g},

(22)

where

a_{S, g}

is the g-th column of

\hat{G} {({\hat{G}}^{H} \hat{G})}^{- 1}

q_{S, g}^{dl} = 1 / \sqrt{E [∥ a_{S, g} ∥^{2}]} = 1 / \sqrt{(M_{S} - G) ε_{g}}

is power normalization parameter, and

E [∥ v_{g}^{ZF} ∥^{2}] = 1

Theorem 2.

When the base station and satellite terminals use ZF precoding, the closed expression of the downlink achievable rate can be given by

R_{g, k}^{ZF} = {log}_{2} (1 + S I N R_{g, k}^{Z F})

S I N R_{g, k}^{Z F} = \frac{{(q_{B, g}^{dl} \sum_{m = 1}^{N_{B}} κ_{m, g, k} + q_{S, g}^{dl} c_{g, k})}^{2}}{\sum_{i = 1}^{G} \sum_{m = 1}^{N_{B}} (β_{m, g, k} - χ_{m, g, k}) + \sum_{i = 1}^{G} (η_{g, k} - γ_{g, k}) - 2 q_{B, g}^{dl} q_{S, g}^{dl} c_{g, k} \sum_{m = 1}^{N_{B}} κ_{m, g, k} + σ_{d l}^{2}} .

(23)

Proof of Theorem 2.

The proof is given in Appendix B. □

Remark 1.

In Theorems 1 and 2, we derive the closed-form expressions of achievable rate with MRT and ZF precoding. In contrast to the derivation process under only terrestrial multicast scenario, satellite-assisted multicast introduces an additional item

g_{g, k}^{H} v_{g}

. It increases the difficulty of derivation. By using large-scale random matrix theory and introducing middle term between terrestrial item and satellite item, we successfully derive the closed-form expressions.

Remark 2.

The closed-form expressions of achievable rate with MRT and ZF precoding presented in Theorems 1 and 2 provide convenience for parameter analyses. Compared with simulation expression, there might be a little different. However, it is relatively tractable and accurate. In particular, with closed-form expressions, there is no need to characterize the distributions of signal and interference powers to obtain an analytical expression, so that we can simplify the analysis process, but achieve relatively accurate results simultaneously.

4. Numerical Results

In this section, the theoretical analysis presented in Section 3 is verified through a set of Monte Carlo simulations, to analyze the performance of the satellite-assisted cell-free massive MIMO multi-group multicast system.

4.1. Simulation Parameter Setting

In the simulation, we consider a circular cell-free massive MIMO system, where all APs and users are randomly distributed and the satellite cover the entire system, so that the satellite and terrestrial base stations provide users with multicast services cooperatively. All the simulation parameters are shown in Table 1.

4.2. Simulation Result Analysis

SE is shown in Figure 2 against signal-noise ratio (SNR) with MRT and ZF precoding under ideal and non-ideal CSI. In addition, the number of AP antennas and satellite antennas is set as 60 and 30, respectively. It can be seen that the performance of ZF precoding is better than MRT precoding. Simultaneously, the SE of the system under ideal CSI is significantly higher than that under non-ideal channels. This also reflects the importance of obtaining accurate channel state information.

SE is shown in Figure 3 against SNR with different number of antennas at APs and satellite. It can be seen that the theoretical results agree well with the simulated results, which verifies the accuracy of the closed-form expressions derived in Theorems 1 and 2. As the number of satellite antennas increases, the SE increases, and the performance of ZF precoding is better than MRT precoding.

Please note that the coexistence of satellite and terrestrial networks may introduce complicated interference problems. When the terrestrial network and the satellite network transmit signals to users at the same time, the repeated use of the same frequency causes inter-system interference which may degrade the system performance severely. How to decrease the inter-system interference will be considered to be a point for further research.

Figure 4 illustrates the average SE against terrestrial downlink transmit power under two scenarios of satellite-assisted multicast and only terrestrial unicast. We set SNR as 10 dB. Meanwhile, we assume that the terrestrial downlink transmitted power of ranging from 30 dBm to 40 dBm and the satellite downlink transmitted power as 30 dBm. Simultaneously, the number of AP antennas and satellite antennas is set as 60 and 30, respectively. It can be seen in Figure 4 that the SE increases with the increasing terrestrial downlink transmitted power. In addition, the system we are studying is satellite-assisted cell-free massive MIMO systems, including satellite networks and terrestrial networks. The satellite network serves as a supplement to the relatively mature terrestrial network. Therefore, when the terrestrial network traffic is overloaded, part of the traffic is offloaded to the satellite network, therefore relieving the communication pressure on the ground. It is reflected in Figure 4. With the help of satellite, when the user’s total rate is constant, part of the rate of ground base station is offloaded to the satellite, resulting that the rate in only terrestrial multicast scene is reduced.

Figure 5 illustrates the influence of different group size on the system SE. The total number of users is set to 60. We change the number of groups, and thus the number of users in each group has also changed. The horizontal coordinate in Figure 5 is the number of users in each group,

\{10, 12, 15, 20, 30\}

correspondingly. Multicast provides a one-time data transmission service for multiple users who need the same resources. It not only avoids repeated transmission of data to improve transmission efficiency, but also ensures that it does not affect other users needing different information to receive data effectively. This can be reflected in Figure 5 that the system SE increases with group size, indicating that multicast scheme helps improve the system performance.

5. Conclusions

In this paper, we analyzed the downlink spectral efficiency of a satellite-assisted cell-free massive MIMO multi-group multicast system, trying to merge the satellite system with the terrestrial cell-free system. First, we use the SR model and the Rayleigh model to establish the satellite channel and the ground channel, and use the common pilot scheme to estimate the satellite and ground CSI. Then, given the estimated CSI, we derived the closed-form expressions of the downlink achievable rate with MRT and ZF beamforming. Finally, the correctness of the derived closed-form expressions was verified, and the simulation results revealed that with the help of multicast, higher system SE can be obtained. Moreover, the integration of satellite and terrestrial network effectively relieved the pressure of terrestrial communication. In the future, we will focus on the study of how to decrease the inter-system interference caused by the coexistence of satellite and terrestrial networks.

Author Contributions

Conceptualization, J.L., L.C., D.W. and P.Z.; writing—original draft preparation, L.C.; writing—review and editing, J.L.; supervision, J.L. and L.C.; project administration, X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China (NSFC) under Grant 61971127, 61871465, 61871122, by the National Key Research and Development Program under Grant 2020YFB1806600, by Natural Science Foundation of Jiangsu Province under Grant BK20180011, and by the Fundamental Research Funds for the Central Universities under grant 2242021k30044.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Hao Wu from Zhongxing Telecommunication Equipment Corporation, Shenzhen 518055, China, also provide conceptual guidance.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of Theorem 1

To obtain the closed-form expression of the downlink achievable rate with MRT beamforming, the following four items need to be calculated:

E [h_{g, k}^{H} w_{g}^{mrt}]

E [g_{g, k}^{H} v_{g}^{mrt}]

\sum_{i = 1}^{G} E [| h_{g, k}^{H} w_{i}^{mrt} |^{2}]

and

\sum_{i = 1}^{G} E [| g_{g, k}^{H} v_{i}^{mrt} |^{2}]

For the term

E [h_{g, k}^{H} w_{g}^{mrt}]

, we have 2.0

\begin{matrix} E [h_{g, k}^{H} w_{g}^{mrt}] & = E [p_{B, g}^{dl} h_{g, k}^{H} {\hat{h}}_{g}] \\ \overset{(a)}{=} E [p_{B, g}^{dl} {\hat{h}}_{g, k}^{H} {\hat{h}}_{g}] \\ \overset{(b)}{=} E [p_{B, g}^{dl} {\hat{h}}_{g}^{H} κ_{g, k}^{H} {\hat{h}}_{g}] \\ \overset{(c)}{=} p_{B, g}^{dl} M_{B} \sum_{m = 1}^{N_{B}} κ_{m, g, k} α_{m, g}, \end{matrix}

(A1)

where

(a)

is because

{\tilde{h}}_{g, k}

and

{\hat{h}}_{g, k}

are independent,

(b)

results from (13) and

(c)

are obtained by the large-scale random matrix conclusion Theorem 1 in [23].

Similarly, for the term

E [g_{g, k}^{H} v_{g}^{mrt}]

, we have

\begin{matrix} E [g_{g, k}^{H} v_{g}^{mrt}] & \overset{(a)}{=} E [p_{S, g}^{dl} c_{g, k} {\hat{g}}_{g}^{H} {\hat{g}}_{g}] \\ \overset{(b)}{=} p_{S, g}^{dl} M_{S} ε_{g} c_{g, k}, \end{matrix}

(A2)

For the term

\sum_{i = 1}^{G} E [{|h_{g, k}^{H} w_{i}^{mrt}|}^{2}]

, we have

\begin{matrix} \sum_{i = 1}^{G} E [{|h_{g, k}^{H} w_{i}^{mrt}|}^{2}] & = \sum_{i = 1}^{G} E [| p_{B, i}^{dl} h_{g, k}^{H} {\hat{h}}_{i} |^{2}] \\ \overset{(a)}{=} \sum_{i = 1}^{G} {(p_{B, i}^{dl})}^{2} E [| {\hat{h}}_{g}^{H} κ_{g, k}^{H} {\hat{h}}_{i} |^{2} + | {\tilde{h}}_{g, k}^{H} {\hat{h}}_{i} |^{2}] \\ \overset{(b)}{=} {(p_{B, g}^{dl})}^{2} E [| {\hat{h}}_{g}^{H} κ_{g, k}^{H} {\hat{h}}_{g} |^{2}] + \sum_{i = 1}^{G} {(p_{B, i}^{dl})}^{2} E [| {\tilde{h}}_{g, k}^{H} {\hat{h}}_{i} |^{2}] \\ \overset{(c)}{=} {(p_{B, g}^{dl} M_{B} \sum_{m = 1}^{N_{B}} κ_{m, g, k} α_{m, g})}^{2} + \sum_{i = 1}^{G} \sum_{m = 1}^{N_{B}} (β_{m, g, k} - χ_{m, g, k}), \end{matrix}

(A3)

where

(a)

results from the independence of

{\tilde{h}}_{g, k}

and

{\hat{h}}_{g, k}

, and the nature of variance.

(b)

is obtained because that

{\hat{h}}_{g}

and

{\hat{h}}_{i}

are uncorrelated random vectors without considering pilot pollution.

(c)

is obtained because of the large-scale random matrix conclusion.

Considering the term

\sum_{i = 1}^{G} E [| g_{g, k}^{H} v_{i}^{mrt} |^{2}]

, we have

\begin{matrix} \sum_{i = 1}^{G} E [{|g_{g, k}^{H} v_{i}^{mrt}|}^{2}] & = \sum_{i = 1}^{G} E [{|p_{S, i}^{dl} g_{g, k}^{H} {\hat{g}}_{i}|}^{2}] \\ \overset{(a)}{=} \sum_{i = 1}^{G} {(p_{S, i}^{dl})}^{2} E [| c_{g, k} {\hat{g}}_{g}^{H} {\hat{g}}_{i} |^{2} + | {\tilde{g}}_{g, k}^{H} {\hat{g}}_{i} |^{2}] \\ \overset{(b)}{=} {(p_{S, g}^{dl})}^{2} E [| c_{g, k} {\hat{g}}_{g}^{H} {\hat{g}}_{g} |^{2}] + \sum_{i = 1}^{G} {(p_{S, i}^{dl})}^{2} E [| {\tilde{g}}_{g, k}^{H} {\hat{g}}_{i} |^{2}] \\ \overset{(c)}{=} {(p_{S, g}^{dl} M_{S} ε_{g} c_{g, k})}^{2} + \sum_{i = 1}^{G} (η_{g, k} - γ_{g, k}), \end{matrix}

(A4)

where

(a)

results from the independence of

{\tilde{g}}_{g, k}

and

{\hat{g}}_{g, k}

, and the nature of variance.

(b)

is obtained because that

{\hat{g}}_{g}

and

{\hat{g}}_{i}

are uncorrelated random vectors without considering pilot pollution.

(c)

is obtained by the large-scale random matrix conclusion.

Substituting the above results into (17) yields the closed-form expression of the downlink achievable rate with MRT beamforming.

Appendix B. Proof of ZF Closed-Form Expression

To obtain the closed-form expression of the downlink achievable rate with ZF beamforming, the following four items need to be calculated:

E [h_{g, k}^{H} w_{g}^{ZF}]

E [g_{g, k}^{H} v_{g}^{ZF}]

\sum_{i = 1}^{G} E [| h_{g, k}^{H} w_{i}^{ZF} |^{2}]

\sum_{i = 1}^{G} E [| g_{g, k}^{H} v_{i}^{ZF} |^{2}]

For the term

E [h_{g, k}^{H} w_{g}^{ZF}]

, we have

\begin{matrix} E [h_{g, k}^{H} w_{g}^{ZF}] & = E [q_{B, g}^{dl} h_{g, k}^{H} \hat{H} {({\hat{H}}^{H} \hat{H})}^{- 1} e_{g}] \\ \overset{(a)}{=} E [q_{B, g}^{dl} {\hat{h}}_{g}^{H} κ_{g, k}^{H} \hat{H} {({\hat{H}}^{H} \hat{H})}^{- 1} e_{g}] \\ \overset{(b)}{=} q_{B, g}^{dl} \sum_{m = 1}^{N_{B}} κ_{m, g, k}, \end{matrix}

(A5)

where

(a)

results from the independence of

{\tilde{h}}_{g, k}

and

{\hat{h}}_{g, k}

and (13),

(b)

is obtained because

E [{\hat{h}}_{g}^{H} \hat{H} {({\hat{H}}^{H} \hat{H})}^{- 1} e_{g}] = 1

Similarly, for the term

E [g_{g, k}^{H} v_{g}^{ZF}]

, we have

\begin{matrix} E [g_{g, k}^{H} v_{g}^{ZF}] & = E [q_{S, g}^{dl} g_{g, k}^{H} \hat{G} {({\hat{G}}^{H} \hat{G})}^{- 1} e_{g}] \\ \overset{(a)}{=} E [q_{S, g}^{dl} c_{g, k} {\hat{g}}_{g}^{H} \hat{G} {({\hat{G}}^{H} \hat{G})}^{- 1} e_{g}] \\ \overset{(b)}{=} q_{S, g}^{dl} c_{g, k}, \end{matrix}

(A6)

For the term

\sum_{i = 1}^{G} E [| h_{g, k}^{H} w_{i}^{ZF} |^{2}]

, we have

\begin{matrix} \sum_{i = 1}^{G} E [{|h_{g, k}^{H} w_{i}^{ZF}|}^{2}] & = \sum_{i = 1}^{G} E [| q_{B, i}^{dl} h_{g, k}^{H} \hat{H} {({\hat{H}}^{H} \hat{H})}^{- 1} e_{i} |^{2}] \\ \overset{(a)}{=} \sum_{i = 1}^{G} {(q_{B, i}^{dl})}^{2} E [| {\hat{h}}_{g}^{H} κ_{g, k}^{H} \hat{H} {({\hat{H}}^{H} \hat{H})}^{- 1} e_{i} |^{2} + | {\tilde{h}}_{g, k}^{H} \hat{H} {({\hat{H}}^{H} \hat{H})}^{- 1} e_{i} |^{2}] \\ \overset{(b)}{=} {(q_{B, g}^{dl})}^{2} E [| {\hat{h}}_{g}^{H} κ_{g, k}^{H} \hat{H} {({\hat{H}}^{H} \hat{H})}^{- 1} e_{g} |^{2}] + \sum_{i = 1}^{G} {(q_{B, i}^{dl})}^{2} E [| {\tilde{h}}_{g, k}^{H} \hat{H} {({\hat{H}}^{H} \hat{H})}^{- 1} e_{i} |^{2}] \\ \overset{(c)}{=} {(q_{B, g}^{dl})}^{2} (\sum_{m = 1}^{N_{B}} κ_{m, g, k}) + \sum_{i = 1}^{G} \sum_{m = 1}^{N_{B}} (β_{m, g, k} - χ_{m, g, k}), \end{matrix}

(A7)

where

(a)

is because

{\tilde{h}}_{g, k}

and

{\hat{h}}_{g, k}

are independent,

(b)

is obtained because that

{\hat{h}}_{g}

and

{\hat{h}}_{i}

are uncorrelated,

(c)

results from the nature of expectation that

E [X_{1} X_{2}] = E [X_{1}] E [X_{2}]

when

X_{1}

and

X_{2}

are independent.

Considering the term

\sum_{i = 1}^{G} E [| g_{g, k}^{H} v_{i}^{ZF} |^{2}]

, we have

\begin{matrix} \sum_{i = 1}^{G} E [{|g_{g, k}^{H} v_{i}^{ZF}|}^{2}] & = \sum_{i = 1}^{G} E [| q_{S, i}^{dl} g_{g, k}^{H} \hat{G} {({\hat{G}}^{H} \hat{G})}^{- 1} e_{i} |^{2}] \\ \overset{(a)}{=} \sum_{i = 1}^{G} {(q_{S, i}^{dl})}^{2} E [| c_{g, k} {\hat{g}}_{g}^{H} \hat{G} {({\hat{G}}^{H} \hat{G})}^{- 1} e_{i} |^{2} + | {\tilde{g}}_{g, k}^{H} \hat{G} {({\hat{G}}^{H} \hat{G})}^{- 1} e_{i} |^{2}] \\ \overset{(b)}{=} {(q_{S, g}^{dl})}^{2} E [| c_{g, k} {\hat{g}}_{g}^{H} \hat{G} {({\hat{G}}^{H} \hat{G})}^{- 1} e_{g} |^{2}] + \sum_{i = 1}^{G} {(q_{S, i}^{dl})}^{2} E [| {\tilde{g}}_{g, k}^{H} \hat{G} {({\hat{G}}^{H} \hat{G})}^{- 1} e_{i} |^{2}] \\ \overset{(c)}{=} {(q_{S, g}^{dl} c_{g, k})}^{2} + \sum_{i = 1}^{G} (η_{g, k} - γ_{g, k}), \end{matrix}

(A8)

where

(a)

is because

{\tilde{g}}_{g, k}

and

{\hat{g}}_{g, k}

are independent,

(b)

is obtained because that

{\hat{g}}_{g}

and

{\hat{g}}_{i}

are uncorrelated,

(c)

is obtained because of the nature of expectation that

E [X_{1} X_{2}] = E [X_{1}] E [X_{2}]

for independent

X_{1}

and

X_{2}

Substituting the above results into (17) yields the closed-form expression of the downlink achievable rate with ZF beamforming.

References

Kuang, L.; Feng, Z.; Qian, Y.; Giambene, G. Interagrated terrestrial-satellite networks: Part two. China Commun. 2018, 15, IV–VI. [Google Scholar]
Wang, P.; Zhang, J.; Zhang, X.; Yan, Z.; Evans, B.G.; Wang, W. Convergence of satellite and terrestrial networks: A comprehensive survey. IEEE Access 2020, 8, 5550–5588. [Google Scholar] [CrossRef]
Zhang, Y.; Yang, L.; Zhu, H. On the spectral efficiency of cell-free large-scale MIMO non-orthogonal multiple access systems. Digit. Signal Process. 2021, 111, 1–11. [Google Scholar] [CrossRef]
Jeong, S.; Farhang, A.; Flanagan, M. Collaborative vs non-collaborative CFO estimation for distributed large-scale MIMO systems. In Proceedings of the 2020 IEEE 92nd Vehicular Technology Conference (VTC2020-Fall), Victoria, BC, Canada, 4–7 October 2020; pp. 1–6. [Google Scholar]
Feng, Y.; Wang, M.; Wang, D.; You, X. Low complexity iterative detection for a large-scale distributed MIMO prototyping system. In Proceedings of the 2019 IEEE International Conference on Communications (ICC), Shanghai, China, 20–24 May 2019; pp. 1–6. [Google Scholar]
Ngo, H.Q.; Larsson, E.G.; Marzetta, T.L. Energy and spectral efficiency of very large multiuser MIMO systems. IEEE Trans. Commun. 2013, 61, 1436–1449. [Google Scholar]
Karakayali, M.K.; Foschini, G.J.; Valenzuela, R.A. Network coordination for spectrally efficient communications in cellular systems. IEEE Wirel. Commun. 2006, 13, 56–61. [Google Scholar] [CrossRef]
Nguyen, L.D.; Duong, T.Q.; Ngo, H.Q.; Tourki, K. Energy efficiency in cell-free massive MIMO with zero-forcing precoding design. IEEE Commun. Lett. 2017, 21, 1871–1874. [Google Scholar] [CrossRef] [Green Version]
Sinky, H.; Khalfi, B.; Hamdaoui, B.; Rayes, A. Responsive content-centric delivery in large urban communication networks: A linknyc use-case. IEEE Trans. Wirel. Commun. 2018, 17, 1688–1699. [Google Scholar] [CrossRef]
Karipidis, E.; Sidiropoulos, N.D.; Luo, Z.Q. Quality of service and max-min fair transmit beamforming to multiple cochannel multicast groups. IEEE Trans. Signal Process. 2008, 56, 1268–1279. [Google Scholar] [CrossRef] [Green Version]
Sadeghi, M.; Bjornson, E.; Larsson, E.G.; Yuen, C.; Marzetta, T. Joint unicast and multi-group multicast transmission in massive MIMO systems. IEEE Trans. Wirel. Commun. 2018, 17, 6375–6388. [Google Scholar] [CrossRef] [Green Version]
Chen, E.; Tao, M.; Liu, Y.F. Joint base station clustering and beamforming for non-orthogonal multicast and unicast transmission with backhaul constraints. IEEE Trans. Wirel. Commun. 2017, 17, 6265–6279. [Google Scholar] [CrossRef]
Doan, T.X.; Ngo, H.Q.; Duong, T.Q.; Tourki, K. On the performance of multigroup multicast cell-free massive MIMO. IEEE Commun. Lett. 2017, 21, 2642–2645. [Google Scholar] [CrossRef]
Zhang, Y.; Cao, H.; Yang, L. Max-min power optimization in multigroup multicast cell-free massive MIMO. In Proceedings of the 2019 IEEE Wireless Communications and Networking Conference (WCNC), Marrakesh, Morocco, 15–18 April 2019; pp. 1–6. [Google Scholar]
Zhang, Y.; Zhu, X.; Jiang, C.; Yin, L. Joint user access and resource association in multicast terrestrial-satellite cooperation network. In Proceedings of the 2018 IEEE Globecom Workshops (GC WKSHPS), Abu Dhabi, United Arab Emirates, 9–13 December 2018; pp. 9–13. [Google Scholar]
Zhu, X.; Jiang, C.; Yin, L.; Kuang, L.; Ge, N.; Lu, J. Cooperative multigroup multicast transmission in integrated terrestrial-satellite networks. IEEE J. Sel. Areas Commun. 2018, 36, 981–992. [Google Scholar] [CrossRef]
Gao, Z.; Liu, A.; Liang, X. The performance analysis of downlink NOMA in LEO satellite communication system. IEEE Access 2020, 8, 93723–93732. [Google Scholar] [CrossRef]
Abdi, A.; Lau, W.; Alouini, M.; Kaveh, M. A new simple model for land mobile satellite channels: First- and second-order statistics. IEEE Trans. Wirel. Commun. 2003, 2, 519–528. [Google Scholar] [CrossRef] [Green Version]
Benkhelifa, F.; Tall, A.; Rezki, Z.; Alouini, M.-S. On the low SNR capacity of MIMO fading channels with imperfect channel state information. IEEE Trans. Commun. 2014, 2, 1921–1930. [Google Scholar] [CrossRef]
Yang, H.; Marzetta, T.L.; Ashikhmin, A. Multicast performance of large-scale antenna systems. In Proceedings of the 2013 IEEE 14th Workshop on Signal Processing Advances in Wireless Communications (SPAWC), Darmstadt, Germany, 16–19 June 2013; pp. 604–608. [Google Scholar]
Li, J.; Wang, D.; Zhu, P.; Wang, J.; You, X. Downlink spectral efficiency of distributed massive MIMO systems with linear beamforming under pilot contamination. IEEE Trans. Veh. Technol. 2018, 67, 1130–1145. [Google Scholar] [CrossRef] [Green Version]
Interdonato, G.; Karlsson, M.; Bjornson, E.; Larsson, E.G. Downlink spectral efficiency of cell-free massive MIMO with full-pilot zero-forcing. In Proceedings of the 2018 IEEE Global Conference on Signal and Information Processing (GlobalSIP), Anaheim, CA, USA, 26–29 November 2018; pp. 1003–1007. [Google Scholar]
Wang, D.; Wang, M.; Zhu, P.; Li, J.; Wang, J.; You, X. Performance of Network-Assisted Full-Duplex for Cell-Free Massive MIMO. IEEE Trans. Commun. 2020, 68, 1464–1478. [Google Scholar] [CrossRef] [Green Version]

Figure 1. System Configuration.

Figure 2. SE against SNR with MRT and ZF beamforming under ideal and non-ideal CSI.

Figure 3. SE against SNR with MRT and ZF beamforming.

Figure 4. SE against downlink transmitted power under satellite-assisted multicast and only terrestrial multicast.

Figure 5. System spectral efficiency of different group size.

Table 1. Simulation parameters.

Parameter Name	Value
Orbit altitude $d_{0}$ 1350 km
Carrier frequency $f_{c}$	20 GHz
The number of AP $N_{B}$	6
The number of multicast group G	3
The number of users per group $K_{g}$	5
Uplink pilot signal transmit power $ρ_{p}$	30 dBm
Cell radius r	1 km
Path loss index a	3.7
Satellite transmitting antenna gain $G_{t}$	24.3 dBi
User receiving antenna gain $G_{r}$	10 dBi
3 dB angle $φ_{3 d B}$	$0 . 4^{\circ}$

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Li, J.; Chen, L.; Zhu, P.; Wang, D.; You, X. Satellite-Assisted Cell-Free Massive MIMO Systems with Multi-Group Multicast. Sensors 2021, 21, 6222. https://doi.org/10.3390/s21186222

AMA Style

Li J, Chen L, Zhu P, Wang D, You X. Satellite-Assisted Cell-Free Massive MIMO Systems with Multi-Group Multicast. Sensors. 2021; 21(18):6222. https://doi.org/10.3390/s21186222

Chicago/Turabian Style

Li, Jiamin, Lingling Chen, Pengcheng Zhu, Dongming Wang, and Xiaohu You. 2021. "Satellite-Assisted Cell-Free Massive MIMO Systems with Multi-Group Multicast" Sensors 21, no. 18: 6222. https://doi.org/10.3390/s21186222

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

Article Menu

Satellite-Assisted Cell-Free Massive MIMO Systems with Multi-Group Multicast

Abstract

1. Introduction

2. System Model

2.1. System Configuration

2.2. Channel Model

3. Performance Analysis

3.1. Channel Estimation

3.2. Downlink Spectral Efficiency Analysis

4. Numerical Results

4.1. Simulation Parameter Setting

4.2. Simulation Result Analysis

5. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Acknowledgments

Conflicts of Interest

Appendix A. Proof of Theorem 1

Appendix B. Proof of ZF Closed-Form Expression

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI