# Hybrid Multicast Transmission for Public Safety Network in 5G

Authors: Fei Qi; Xuetian Zhu; Xinran Zhang; Songlin Sun
Journal: Mobile Information Systems
 (2016)

---

## Abstract

We investigate the application of wireless multicast technology in public safety network (PSN) in future wireless communication system. The hybrid unicast/multicast transmission system is proposed and analyzed in 3D massive multi-input multioutput (MIMO) channel. The mutual coupling channel model is adopted under the different antenna array configuration scenarios. The proposed hybrid system adopts multicast beamforming in the multicast groups as well as multiuser-MIMO (MU-MIMO) linear precoding in the unicast group to increase system throughput. The null space method based interference cancellation is further performed between each group to eliminate signal leakage generated from each group. Comparisons between two types of antenna array configurations, different channel models, linear precoding as well as multicast beamforming, and user grouping strategies for multicast services are presented and analyzed by simulation.


---

## Body

## 1. Introduction

The public safety network (PSN) refers to wireless communication networks used by first responders (FRs) such as police, firefighters, and emergency medical technicians (EMTs) for public protection and disaster relief [1, 2]. It is an important issue in the design of commercial wireless communication system [3]. Although the network used by FRs is generally private system, the long-term evolution (LTE) system is designed to provide functions in enabling FRs agencies to collaborate with each other as well as provide access service for these agencies when needed. In future fifth-generation (5G) system which features in outstanding physical-layer transmission technologies such as millimeter-wave (mmWave) communication and massive multiple-input multipleoutput (MIMO) transmission [4, 5], the PSN and FR service model has to be reconsidered and reevaluated under the new scenarios. For example, in LTE system and future 5G system, different antenna array configuration schemes in a massive MIMO manner allow the transmitting evolved Node B (eNB) to control both the azimuth and elevation angles of wireless signal propagation and are consequently referred as the three-dimensional (3D) MIMO as well [6]. Its flexibility in network deployment and adequate throughput volume advantage will help FRs to facilitate better with cellular network in PSN.In PSN, the multicast is a promising technology to provide services for FRs [1, 2, 7]. It has been proved to be an efficient and important communication technology and has been standardized in 3rd generation partner project (3GPP) protocols in the form of evolved multimedia broadcast and multicast services (eMBMS). Public Safety Networks (PSNs) have been recognized as the key to a successful response to emergency and disaster situations. In a PSN, first responders often cooperatively accomplish one mission to release as many survivors as possible. Example applications include audio/video conferencing as well as one-to-many data dissemination in critical situations. Multicast is a very efficient way to support group-oriented applications, especially in mobile environments where the spectrum is valuable. This work, correspondingly, addresses that the wireless communication in 5G PSN involves both unicast and multicast traffic and thus formulates a hybrid transmission scenario. In this work we investigate the multicast for PSN in 5G scenario by focusing on the 3D massive MIMO channel model based on the mutual coupling channel model. In the form of physical-layer multicast beamforming, multicast can significantly increase system capacity and spectrum efficiency [7–13]. In terms of theoretical analysis, the general wireless multicast beamforming problem is NP-hard problem. In massive MIMO scenario, this problem can be solved in asymptotic analysis with closed-form solutions [10, 11], which is important for the design for PSN. Existing literatures have focused on the asymptotic analysis for multicast problem under both perfect channel state information (CSI) scenario and imperfect CSI scenario, which are performed under Rayleigh fading channel model.To explore the technical insight of multicast in PSN of 5G scenario, we adopt the 3D MIMO mutual coupling channel model [14, 15] to describe the multicast beamforming problem. To proceed, we propose a hybrid transmission system in which the multicast groups adopt multicast beamforming while the unicast group adopts multiuser-MIMO (MU-MIMO) linear precoding. The null space method based interference cancellation between each group is also performed in the proposed system. The optimization of hybrid unicast/multicast precoding is barely investigated and highly challenging. As indicated in [10], even the multicast beamforming problem by itself is already NP-hard, and its optimal solution cannot be obtained efficiently in general. For this reason, this work focuses on the development of a heuristic approach based on null space interference cancellation. We find our proposed hybrid precoding scheme outperforms other heuristic approaches with better efficiency and complexity tradeoff. The system is simulated and analyzed under linear antenna array configuration and rectangular antenna array configuration. We present the comparison amongst different antenna array configurations, channel models, and precoding schemes.The rest of the paper is organized as follows. Section2 presents the 3D mutual coupling channel model for the linear antenna array configuration and rectangular antenna array configuration. Section 3 presents the proposed system model. The multicast problem under the mutual coupling channel model and interference cancellation method are modeled and described. The hybrid MU-MIMO linear precoding/multicast beamforming transmission scheme is also described in this section. Section 4 presents the illustrated simulation results, and Section 5 concludes the paper.Notation. Boldface uppercase and lowercase letters denote matrices and column vectors, respectively. N+ and C denote the set of positive integer and complex number, respectively. ·∗, ·T, and ·H denote complex conjugate, transpose, and Hermitian transpose, respectively. IN denotes the N×N identity matrix.


## 2. System Architecture

### 2.1. Architecture for PSN

Public safety networks (PSNs) support communications for special services including police, fire, and ambulance. In 5G PSN, the requirement has been to develop systems that are highly robust and can address the specific communication needs of emergency services. 5G system features in providing a massive amount of miscellaneous devices to connect into the packet data network (PDN), the core network. For example, in case of an emergent disaster the operator has to provide the first responders with reliable information as soon as possible. To achieve high network efficiency, 5G PSN will establish a heterogeneous system by incorporating the software defined network (SDN) and radio access network (RAN). Operators may provide multicast services for PSN to efficiently utilize limited bandwidth in the confined RAN coverage area by delivering the identical emergency information for FRs.For the 5G PSN considering heterogeneity and multicast services, we proposed the PSN architecture as shown in Figure1. In a future PSN scenario, the heterogeneous PDN in the form of SDN is connected to the RAN of PSN as illustrated. The macrocells are coupled with multiple low power nodes (LPNs) in a ultra-dense manner. The LPNs may apply mmWave communication to provide a much higher speed of data service access to their serving FRs. The FRs within the coverage of the cloud RAN may subscribe to a specific multicast group and utilize multicast beamforming transmission from the heterogeneous access point to receive identical emergency information or subscribe to no multicast groups and utilize MU-MIMO linear precoding to receive high-speed data transmission. In this hybrid service architecture, interference cancellation technology is essential to the system performance [16], since the multicast beamforming transmission and MU-MIMO linear precoding transmission occupy the same time-frequency resources within the coverage of cloud RAN. Consequently, it is crucial to perform the null space method based interference cancellation algorithm to ensure that the manipulated channels of each multicast group as well as the unicast group are orthogonal to each other, resulting in no interference from theoretical point of view.Figure 1
Hybrid multicast transmission in PSN.


### 2.2. Multicast in PSN

One of the most challenging characteristics of PSN is the abrupt intense communication traffic. Although 5G features in providing seamless coverage by densified deployment of LPNs, the traffic peak introduced by PSN can easily paralyse the core network as well as the heterogeneous access points. Also, physical damage to the access points in PSN will contribute to the traffic burden for the whole system. As a result, multicast has been considered as a crucial solution for PSN communication since it can fast spread identical message to the coverage area without introducing further communication costs. In practical system, the multicast can be performed in IP-layer by utilizing various multicast protocols. The IP multicast provides a more efficient method for many-to-many communication and has become more and more important for both commercial networks and dedicated systems. IP multicast allows the source to send a single copy of data, using a single address for the entire group of recipients. The most used multicast routing protocol is the Protocol Independent Multicast Sparse Mode (PIM-SM). Routers between the source and recipients use the group address to route the data. The routers forward duplicate data packets only when required; that is, the path to recipients diverges.Although IP-layer multicast has been successfully implemented in practical system, it introduces extra communication complexity for routers and signaling overhead. The wireless communication system, however, has multicast nature in it and is ideal to perform multicast transmission for PSN. In the form of physical-layer wireless multicast, the multicast FRs share the same access medium, that is, wireless spectrum and time; hence it will greatly increase system capacity and lighten the peak traffic problem in PSN. Combined with large-scale antenna array technology in 5G, the multicast is a suitable and effective solution to deliver communication services for FRs. We model the hybrid multicast system for PSN in the following sections.


## 2.1. Architecture for PSN

Public safety networks (PSNs) support communications for special services including police, fire, and ambulance. In 5G PSN, the requirement has been to develop systems that are highly robust and can address the specific communication needs of emergency services. 5G system features in providing a massive amount of miscellaneous devices to connect into the packet data network (PDN), the core network. For example, in case of an emergent disaster the operator has to provide the first responders with reliable information as soon as possible. To achieve high network efficiency, 5G PSN will establish a heterogeneous system by incorporating the software defined network (SDN) and radio access network (RAN). Operators may provide multicast services for PSN to efficiently utilize limited bandwidth in the confined RAN coverage area by delivering the identical emergency information for FRs.For the 5G PSN considering heterogeneity and multicast services, we proposed the PSN architecture as shown in Figure1. In a future PSN scenario, the heterogeneous PDN in the form of SDN is connected to the RAN of PSN as illustrated. The macrocells are coupled with multiple low power nodes (LPNs) in a ultra-dense manner. The LPNs may apply mmWave communication to provide a much higher speed of data service access to their serving FRs. The FRs within the coverage of the cloud RAN may subscribe to a specific multicast group and utilize multicast beamforming transmission from the heterogeneous access point to receive identical emergency information or subscribe to no multicast groups and utilize MU-MIMO linear precoding to receive high-speed data transmission. In this hybrid service architecture, interference cancellation technology is essential to the system performance [16], since the multicast beamforming transmission and MU-MIMO linear precoding transmission occupy the same time-frequency resources within the coverage of cloud RAN. Consequently, it is crucial to perform the null space method based interference cancellation algorithm to ensure that the manipulated channels of each multicast group as well as the unicast group are orthogonal to each other, resulting in no interference from theoretical point of view.Figure 1
Hybrid multicast transmission in PSN.


## 2.2. Multicast in PSN

One of the most challenging characteristics of PSN is the abrupt intense communication traffic. Although 5G features in providing seamless coverage by densified deployment of LPNs, the traffic peak introduced by PSN can easily paralyse the core network as well as the heterogeneous access points. Also, physical damage to the access points in PSN will contribute to the traffic burden for the whole system. As a result, multicast has been considered as a crucial solution for PSN communication since it can fast spread identical message to the coverage area without introducing further communication costs. In practical system, the multicast can be performed in IP-layer by utilizing various multicast protocols. The IP multicast provides a more efficient method for many-to-many communication and has become more and more important for both commercial networks and dedicated systems. IP multicast allows the source to send a single copy of data, using a single address for the entire group of recipients. The most used multicast routing protocol is the Protocol Independent Multicast Sparse Mode (PIM-SM). Routers between the source and recipients use the group address to route the data. The routers forward duplicate data packets only when required; that is, the path to recipients diverges.Although IP-layer multicast has been successfully implemented in practical system, it introduces extra communication complexity for routers and signaling overhead. The wireless communication system, however, has multicast nature in it and is ideal to perform multicast transmission for PSN. In the form of physical-layer wireless multicast, the multicast FRs share the same access medium, that is, wireless spectrum and time; hence it will greatly increase system capacity and lighten the peak traffic problem in PSN. Combined with large-scale antenna array technology in 5G, the multicast is a suitable and effective solution to deliver communication services for FRs. We model the hybrid multicast system for PSN in the following sections.


## 3. Channel Model for Multicast in PSN

In preparation for multicast transmission modeling, we first present the channel modeling work. Considering downlink massive MIMO multicast system where eNB equipped withN antenna broadcasts grouped messages to K FRs (here and after referred to as FRs) equipped with single antenna within eNB’s coverage, N,K∈N+, N≥K. Assuming the non-line-of-sight (NLOS) scenario, the channel model can be described by [6](1)G=HD,where G=[g1,g2,…,gK] denotes the channel matrix and gk∈CN×1(k=1,2,…,K) denotes the channel vector between eNB and the kth FR. D=diag⁡β1,β2,…,βK denotes the large-scale propagation matrix with βk=κdk-γξk. κ is a constant value determined by antenna characteristic and carrier frequency, dk denotes the distance between eNB and the kth FR, γ denotes the path loss exponent, and ξk denotes the log-normal shadow fading factor.Consider massive MIMO antenna array and 3D MIMO channel model. Since the number of antennas rises to a larger order of magnitude than tradition MIMO antenna array, the mutual impedance of adjacent antennas has to be considered for a more accurate channel model. We adopt the mutual coupling channel model in whichH=[h1,h2,…,hK] denotes the fast fading matrix, and the column vector hk∈CN×1(k=1,2,…,K) is described by the following equation [14, 15]:(2)hk=ZRkvk.Z∈CN×N denotes the constant mutual coupling matrix related to antenna configuration and has the following form:(3)Z=ZA+ZLΓ+ZLI-1,where(4)Γ=ZAZM0⋯0ZMZAZM⋯00ZMZA⋯0⋮⋮⋱⋱⋮00⋯ZMZA.The complex value of ZA, ZL, and ZM denotes the antenna impedance, load impedance, and mutual impedance, respectively. Rk denotes steering matrix and vk is the Gaussian stochastic factor.The form and size ofRk and vk depend on antenna array configuration [6]. For linear antenna array scenario, suppose that the kth FR has Ak different angles of arrivals (AoAs) denoted by θk,l(l=1,2,…,Ak). In that case Rk has Ak column of steering vectors; that is, Rk=[rk,1,rk,2,…,rk,Ak], where(5)rk,l=aθk,l=1Ak1,ej2πd/λsin⁡θk,l,ej2πd/λ2sin⁡θk,l,…,ej2πd/λN-1sin⁡θk,lT.d denotes the distance between adjacent antennas, and λ denotes the wave length. a(θ) is defined as the steering vector function for later convenience.For the rectangular antenna array scenario, thekth FR has Ak different azimuths of arrivals and elevations of arrivals (EoAs) denoted by θk,i and ϕk,i, respectively. In that case the column steering vectors in Rk=[rk,1,rk,2,…,rk,Ak] is described by the following equation:(6)rk,l=1Ak×vec⁡aθk,l⊗aϕk,lT=1Ak×vec⁡1,ej2πd/λsin⁡θk,l,ej2πd/λ2sin⁡θk,l,…,ej2πd/λN-1sin⁡θk,l⊗1,ej2πd/λsin⁡ϕk,l,ej2πd/λ2sin⁡ϕk,l,…,ej2πd/λN-1sin⁡ϕk,l,where the symbol ⊗ represents the Kronecker product of matrices, and the vec⁡(·) function denotes vectorization of matrix.With the established channel model above, the transmission and service model can be accurately described. In the following section we present the proposed hybrid transmission system for multicast services.


## 4. Hybrid Multicast Transmission for PSN

### 4.1. FR Grouping

ConsiderM different multicast groups indexed as 1,2,…,M, and define the group indexed with 0 as the unicast group for nonmulticast FRs. Each user can subscribe to no more than one specific multicast group to receive multicasted information. To characterize that, define the subscription vector m={m1,m2,…,mK}, where the element mk∈{0,1,2,…,M}(k=1,2,…,K) represents that the kth user is subscribed to multicast group indexed with mk (for unicast FRs, mk=0). For the grouping method, define gm={g1,g2,…,gK} for group m(m=0,1,2,…,M), where gk∈{0,1} (k=1,2,…,K) denotes grouping scheme with gk=1 being the kth user subscribed to group m and gk=0 being otherwise. For FRs in group m, let Km={k∣mk=m} and let Km∈N+ denote the set of indices of FRs in group m. Let Km denote the number of FRs in group m with ∑m=0MKm=K.


### 4.2. Interference Cancellation

Define channel grouping matrixGm=gm,1,gm,2,…,gm,Km, where vector gm,j∈CN×1(j=1,2,…,Km) denotes the channel vector for the jth user in group m. Before describing the interference cancellation, the receiving signal model is addressed in (7) firstly for the full view of the whole system.(7)Y=PGJωS+z,z=z0,z1,…,zMT,S=s0,s1,…,sMT,ω=diag⁡Ω0,w1,…,wM,J=J0,J1,…,JM,G=G0,G1,…,GMT,Y=y0,y1,…,yMT,where P is the signal power in the transmitter and Y denotes the receiving signal of the transmitted signal S, which has been processed by the projection matrix J and precoding matrix ω. The white Gaussian noise is denoted as z. For the matrix wm, it can be described as wm=1/λ[g^m,1∗,g^m,2∗,…,g^m,Km∗], in which λ is the normalization factor, which can make sure ωm=1/λ∑j=1Kmg^m,1∗ and Ω0 is the precoding matrix of the unicast users.Considering downlink multicast service transmission, in order to cancel the interferences between all groups, we perform the null space method based interference cancellation algorithm [17] in the system. In order to cancel the interferences between all M multicast groups as well as the unicast group, the projection matrix for the channel matrix of group m must satisfy that it spans a subspace of the null space of the channel vectors of underlying victim FRs in all other M groups. For group m, the joint channel matrix for FRs in all other M groups is defined as follows:(8)Gm¯=G0,G1,…,Gm-1,Gm+1,…,GM.By singular value decomposition (SVD), the projection matrix for group m denoted by Jm∈CN×N is derived by the following equations:(9)Gm¯=Um¯Λm¯Vm¯,Jm=Vm¯0Vm¯0H,Um¯∈CN×N,Vm¯∈CK-Km×K-Km,Λm¯∈CN×K-Km.Um¯, Vm¯, and Λm¯ are SVD factor matrices. Vm¯0∈CK-Km×N is a subspace of column vector space of Vm¯. By multiplying Jm, the new channel matrix for group m denoted as G^m and new channel vector for the jth user in group m denoted as g^m,j are transformed as follows:(10)G^m=JmGm,g^m,j=Jmgm,j.After the manipulation, the signal leakage from each group can be eliminated, guaranteeing the performance of the hybrid multicast/unicast transmission. In the following parts we will adopt the new channel matrices for analysis.


### 4.3. Formulation of FR Multicast

In this part we model the multicast problem based on previous channel modeling work. For simplicity, consider local cell multicast scenario only. Letωm∈CN×1 be the multicast beamforming vector for multicast group m with ωm2=1. Let sm be the stochastic information with unit power multicasted in group m. Then the transmitted signal in group m is given by(11)xm=pmωmsm,where pm is the transmission power for group m. Assume that all groups share the same spectrum, then the received signal for the jth user in group m is given by(12)ym,j=g^m,jHxm+∑n≠mMg^m,jHxn+zm,j=pmg^m,jHωmsm+∑n≠mMpng^m,jHωnsn+zm,j,where zm,j is the additive Gaussian noise with zero mean and variance being σ2. The unperfect precoding scheme causes the interference from intragroup to be the main interference as intergroup interference would be 0 after proper precoding. Then the received SINR for the jth user in group m is(13)SINRm,j=pmg^m,jHωmsm2∑n≠mMpng^m,jHωnsn2+σ2.In massive MIMO scenario, it is assumed that transmission power per antenna scales inversely to the number of antennas [5]; that is, pm=ρmE/N, where E denotes the total transmission power of the base station and ρm∈[0,1] denotes the power ratio of group m. Consider the perfect CSI scenario for simplicity; the multicast problem is formulated as the max–min fairness (MMF) problem [10] as follows:(14a)P:maxωm,ρm,m⁡min0⩽m⩽Mmin∀j∈Km⁡ρmE/Ng^m,jHωmsm2∑n≠mMρnE/Ng^m,jHωnsn2+σ2(14b)s.t.ωm2=1,∀m(14c)∑m=0Mρm⩽1,ρm∈0,1(14d)m=m1,m2,…,mK,mk∈0,1,…,M,k∈1,2,…,K(14e)Km=k∣mk=m.


### 4.4. Asymptotic Optimal Multicast Beamforming

The multicast problemP can be solved under massive MIMO scenario by asymptotic solution. Similar to existing literatures [10, 11], we readily have the following theorem.Theorem 1.
The asymptotic optimal beamforming vectorωm∗, optimal power ratio ρm∗, and asymptotic minimum SINR achieved in group m are given by(15a)ωm∗=αm∑j∈Kmg^m,jHβj,∀m,(15b)αm=N∑j∈Kmβj-1-1/2,(16)ρm∗=∑j∈Kmβj-1∑k=1Kβk-1,∀m,SINRm∞=Eσ2∑k=1Kβk-1,∀m.The theorem indicates that whenN→∞, the SINR achieved in each group is irrelevant to the user grouping scheme and only depends on large-scale channel attenuation and total number of FRs in the macrocell. In practical massive MIMO scenario the antenna number is large enough to observe the asymptotic behavior above. Existing literatures also indicate that this phenomenon can be affected by pilot contamination in a time-division duplexing (TDD) system. Unlike previous works that focus on the analysis of problem (14a), (14b), (14c), (14d), and (14e) in different scenarios, we focus on the hybrid transmission scheme in practical massive MIMO system and directly adopt the asymptotic results from Theorem 1 as the multicast beamformer and power ratio for the simulation system.


### 4.5. Hybrid Multicast Transmission

With the above multicast model, we investigate several types of MU-MIMO linear precoding schemes [17, 18] for the unicast group that is widely adopted in MIMO system for comparison. For the K0×1 multiuser transmit vector, the precoding matrix Ω0∈CN×K0 is multiplied to the transmit vector to perform MU-MIMO linear precoding as follows.MRT Precoding. The maximum ratio transmission (MRT) precoding algorithm constructs a matched filter in spatial domain. The MRT precoding matrix for the unicast group is denoted as follows:(17)Ω0,MRT=1β0G^0∗,β0=trB0B0Hp0,Bm=G0∗.BD Precoding. The block diagonalization (BD) precoding algorithm is derived with the same manner as the null space method based interference cancellation method from (9) and (10), where the precoding matrix Ω0,BD guarantees that the precoded channel matrix of a specific FR lies within the null space of the joint channel matrix of other FRs in the same group; hence interferences from all other FRs are eliminated. The mathematical notations are similar to (9) and (10) and are omitted for simplicity.These MU-MIMO linear precoding schemes are performed in the unicast group while the multicast beamforming scheme from Theorem1 is performed in the multicast groups. The proposed hybrid system can provide high-speed transmission for unicast FRs and simultaneously provide multicast services to multicast FRs, increasing the spectrum efficiency of the system greatly.


## 4.1. FR Grouping

ConsiderM different multicast groups indexed as 1,2,…,M, and define the group indexed with 0 as the unicast group for nonmulticast FRs. Each user can subscribe to no more than one specific multicast group to receive multicasted information. To characterize that, define the subscription vector m={m1,m2,…,mK}, where the element mk∈{0,1,2,…,M}(k=1,2,…,K) represents that the kth user is subscribed to multicast group indexed with mk (for unicast FRs, mk=0). For the grouping method, define gm={g1,g2,…,gK} for group m(m=0,1,2,…,M), where gk∈{0,1} (k=1,2,…,K) denotes grouping scheme with gk=1 being the kth user subscribed to group m and gk=0 being otherwise. For FRs in group m, let Km={k∣mk=m} and let Km∈N+ denote the set of indices of FRs in group m. Let Km denote the number of FRs in group m with ∑m=0MKm=K.


## 4.2. Interference Cancellation

Define channel grouping matrixGm=gm,1,gm,2,…,gm,Km, where vector gm,j∈CN×1(j=1,2,…,Km) denotes the channel vector for the jth user in group m. Before describing the interference cancellation, the receiving signal model is addressed in (7) firstly for the full view of the whole system.(7)Y=PGJωS+z,z=z0,z1,…,zMT,S=s0,s1,…,sMT,ω=diag⁡Ω0,w1,…,wM,J=J0,J1,…,JM,G=G0,G1,…,GMT,Y=y0,y1,…,yMT,where P is the signal power in the transmitter and Y denotes the receiving signal of the transmitted signal S, which has been processed by the projection matrix J and precoding matrix ω. The white Gaussian noise is denoted as z. For the matrix wm, it can be described as wm=1/λ[g^m,1∗,g^m,2∗,…,g^m,Km∗], in which λ is the normalization factor, which can make sure ωm=1/λ∑j=1Kmg^m,1∗ and Ω0 is the precoding matrix of the unicast users.Considering downlink multicast service transmission, in order to cancel the interferences between all groups, we perform the null space method based interference cancellation algorithm [17] in the system. In order to cancel the interferences between all M multicast groups as well as the unicast group, the projection matrix for the channel matrix of group m must satisfy that it spans a subspace of the null space of the channel vectors of underlying victim FRs in all other M groups. For group m, the joint channel matrix for FRs in all other M groups is defined as follows:(8)Gm¯=G0,G1,…,Gm-1,Gm+1,…,GM.By singular value decomposition (SVD), the projection matrix for group m denoted by Jm∈CN×N is derived by the following equations:(9)Gm¯=Um¯Λm¯Vm¯,Jm=Vm¯0Vm¯0H,Um¯∈CN×N,Vm¯∈CK-Km×K-Km,Λm¯∈CN×K-Km.Um¯, Vm¯, and Λm¯ are SVD factor matrices. Vm¯0∈CK-Km×N is a subspace of column vector space of Vm¯. By multiplying Jm, the new channel matrix for group m denoted as G^m and new channel vector for the jth user in group m denoted as g^m,j are transformed as follows:(10)G^m=JmGm,g^m,j=Jmgm,j.After the manipulation, the signal leakage from each group can be eliminated, guaranteeing the performance of the hybrid multicast/unicast transmission. In the following parts we will adopt the new channel matrices for analysis.


## 4.3. Formulation of FR Multicast

In this part we model the multicast problem based on previous channel modeling work. For simplicity, consider local cell multicast scenario only. Letωm∈CN×1 be the multicast beamforming vector for multicast group m with ωm2=1. Let sm be the stochastic information with unit power multicasted in group m. Then the transmitted signal in group m is given by(11)xm=pmωmsm,where pm is the transmission power for group m. Assume that all groups share the same spectrum, then the received signal for the jth user in group m is given by(12)ym,j=g^m,jHxm+∑n≠mMg^m,jHxn+zm,j=pmg^m,jHωmsm+∑n≠mMpng^m,jHωnsn+zm,j,where zm,j is the additive Gaussian noise with zero mean and variance being σ2. The unperfect precoding scheme causes the interference from intragroup to be the main interference as intergroup interference would be 0 after proper precoding. Then the received SINR for the jth user in group m is(13)SINRm,j=pmg^m,jHωmsm2∑n≠mMpng^m,jHωnsn2+σ2.In massive MIMO scenario, it is assumed that transmission power per antenna scales inversely to the number of antennas [5]; that is, pm=ρmE/N, where E denotes the total transmission power of the base station and ρm∈[0,1] denotes the power ratio of group m. Consider the perfect CSI scenario for simplicity; the multicast problem is formulated as the max–min fairness (MMF) problem [10] as follows:(14a)P:maxωm,ρm,m⁡min0⩽m⩽Mmin∀j∈Km⁡ρmE/Ng^m,jHωmsm2∑n≠mMρnE/Ng^m,jHωnsn2+σ2(14b)s.t.ωm2=1,∀m(14c)∑m=0Mρm⩽1,ρm∈0,1(14d)m=m1,m2,…,mK,mk∈0,1,…,M,k∈1,2,…,K(14e)Km=k∣mk=m.


## 4.4. Asymptotic Optimal Multicast Beamforming

The multicast problemP can be solved under massive MIMO scenario by asymptotic solution. Similar to existing literatures [10, 11], we readily have the following theorem.Theorem 1.
The asymptotic optimal beamforming vectorωm∗, optimal power ratio ρm∗, and asymptotic minimum SINR achieved in group m are given by(15a)ωm∗=αm∑j∈Kmg^m,jHβj,∀m,(15b)αm=N∑j∈Kmβj-1-1/2,(16)ρm∗=∑j∈Kmβj-1∑k=1Kβk-1,∀m,SINRm∞=Eσ2∑k=1Kβk-1,∀m.The theorem indicates that whenN→∞, the SINR achieved in each group is irrelevant to the user grouping scheme and only depends on large-scale channel attenuation and total number of FRs in the macrocell. In practical massive MIMO scenario the antenna number is large enough to observe the asymptotic behavior above. Existing literatures also indicate that this phenomenon can be affected by pilot contamination in a time-division duplexing (TDD) system. Unlike previous works that focus on the analysis of problem (14a), (14b), (14c), (14d), and (14e) in different scenarios, we focus on the hybrid transmission scheme in practical massive MIMO system and directly adopt the asymptotic results from Theorem 1 as the multicast beamformer and power ratio for the simulation system.


## 4.5. Hybrid Multicast Transmission

With the above multicast model, we investigate several types of MU-MIMO linear precoding schemes [17, 18] for the unicast group that is widely adopted in MIMO system for comparison. For the K0×1 multiuser transmit vector, the precoding matrix Ω0∈CN×K0 is multiplied to the transmit vector to perform MU-MIMO linear precoding as follows.MRT Precoding. The maximum ratio transmission (MRT) precoding algorithm constructs a matched filter in spatial domain. The MRT precoding matrix for the unicast group is denoted as follows:(17)Ω0,MRT=1β0G^0∗,β0=trB0B0Hp0,Bm=G0∗.BD Precoding. The block diagonalization (BD) precoding algorithm is derived with the same manner as the null space method based interference cancellation method from (9) and (10), where the precoding matrix Ω0,BD guarantees that the precoded channel matrix of a specific FR lies within the null space of the joint channel matrix of other FRs in the same group; hence interferences from all other FRs are eliminated. The mathematical notations are similar to (9) and (10) and are omitted for simplicity.These MU-MIMO linear precoding schemes are performed in the unicast group while the multicast beamforming scheme from Theorem1 is performed in the multicast groups. The proposed hybrid system can provide high-speed transmission for unicast FRs and simultaneously provide multicast services to multicast FRs, increasing the spectrum efficiency of the system greatly.


## 5. Simulation Results

We evaluate the performance of the proposed hybrid system using MATLAB simulation. We focus on the channel model, antenna array configuration, and precoding/beamforming comparison as our major concerns in the simulation system. The linear antenna array and rectangular antenna array are chosen for analysis. Simulations of MRT/BD precoding as well as the multicast beamforming are run with and without the mutual coupling channel model for comparison. For all the multicast service groups, the illustrated performance metrics are calculated as the average metrics. The beamformer is calculated by Theorem1, while the two MU-MIMO linear precoding schemes are performed in the unicast group. The major parameters of the simulation system are listed in Table 1.Table 1
Simulation parameters.
Parameters Values Transmit antenna number N = 6 4 FR number K=16/32 Antenna interval for linear antenna array 0.5 wavelengths Antenna interval for rectangular antenna array 0.5 wavelengths for horizontal and vertical direction ZA 50Ω ZL 50Ω ZM 50Ω SNR −10 dB to 25 dB Frequency 38 GHz Tx antenna height 36 m Rx antenna height 1.5 m Tx Rx antenna gain 25 / 13 . 3 Scenario NLOS


### 5.1. Analysis for the Mutual Coupling Effect

Figures2 and 3 illustrate the spectrum efficiency of linear antenna array and rectangular antenna array, respectively. The performance of MRT/BD precoding (here and after denoted as MRT and BD in the legends) and multicast beamforming (here and after denoted as multi in the legends) are evaluated with and without mutual coupling (here and after denoted as MC in the legends) based channel model. Since the mutual coupling effect introduces extra correlation between adjacent antennas, the channel becomes less independent in a spatial prospect hence the performance of all precoding and beamforming schemes are expected to be degraded. This is clearly observed in Figures 2 and 3 as the spectrum efficiency degrades for all precoding and beamforming schemes for both antenna array configurations. Also, since the mutual coupling effect of rectangular antenna array is considered to be more prominent than the linear antenna array, the performance degradation of rectangular antenna array is worse compared to linear antenna array as observed.Figure 2
Spectrum efficiency (linear antenna array; number of UEs: 16; +MC: mutual coupling; precoding schemes: BD/MRT/multicast).Figure 3
Spectrum efficiency (rectangular antenna array; number of UEs: 16; +MC: mutual coupling; precoding schemes: BD/MRT/multicast).


### 5.2. Comparison between Linear Precoding and Multicast Beamforming

To compare the linear precoding with multicast beamforming, we first choose the BD precoding for unicast FRs. Figures4 and 5 illustrate the spectrum efficiency comparison between BD precoding and multicast beamforming for linear antenna array and rectangular antenna array, respectively. The total spectrum efficiency for the entire system is illustrated and becomes legend as “all group” for observation convenience. For linear antenna array, the performance of multicast beamforming is worse than BD precoding. For rectangular antenna array, however, despite the fact that the performance of BD precoding and multicast beamforming both degrade, the multicast beamforming suffers much less degradation than BD precoding. Similar phenomena are observed in Figures 6 and 7 which illustrate the comparison between MRT precoding and multicast beamforming. This is because the multicasted information for each beamforming user is the same. With a higher correlation between antennas in rectangular antenna array, the signal leakage of multicast beamforming between each user is more of in-phase component and contributes less to the interference for the intended multicast user than that of the unicast group. Hence the multicast beamforming has better performance than linear precoding in terms of channel degradation introduced by concentrated antenna array configuration.Figure 4
Comparison between BD precoding and multicast beamforming (linear antenna array; number of unicast UEs: 8; number of multicast group UEs: 8; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding Scheme: BD; multicast precoding scheme: multicast precoding).Figure 5
Comparison between BD precoding and multicast beamforming (rectangular antenna array; number of unicast UEs: 8; number of multicast group UEs: 8; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: BD; multicast precoding scheme: multicast precoding).Figure 6
Comparison between MRT precoding and multicast beamforming (linear antenna array; number of unicast UEs: 8; number of multicast group UEs: 8; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: MRT; multicast precoding scheme: multicast precoding).Figure 7
Comparison between MRT precoding and multicast beamforming (rectangular antenna array; number of unicast UEs: 8; number of multicast group UEs: 8; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: MRT; multicast precoding scheme: multicast precoding).


### 5.3. Analysis for the Grouping Strategy

Figures8, 9, and 10 illustrate the trend of total spectrum efficiency with an increasing multicast user number while total user number is fixed to be 32. Specifically, a scenario with a much higher correlation between adjacent antennas for rectangular antenna array is presented in Figure 10. With an increasing number of multicast FRs, the total spectrum efficiency decreases since the multicast beamforming achieves worse performances than linear precoding in low correlation scenario. Meanwhile in Figure 10 we observe that, in high correlation scenario, the multicast beamforming outperforms linear precoding greatly and hence contributes to the growing trend of total spectrum efficiency when assigning more FRs to multicast groups, which conforms with the analysis in previous subsection. This is significant for the design of antenna array in different multicast service model to achieve best performance in practical wireless communication system.Figure 8
Grouping strategy analysis (linear antenna array; number of unicast UEs: 31 to 1; number of multicast group UEs: 1 to 31; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: BD; multicast precoding scheme: multicast precoding).Figure 9
Grouping strategy analysis (rectangular antenna array, low antenna relevance; number of unicast UEs: 31 to 1; number of multicast group UEs: 1 to 31; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: BD; multicast precoding scheme: multicast precoding).Figure 10
Grouping strategy analysis (rectangular antenna array, high antenna relevance; number of unicast UEs: 31 to 1; number of multicast group UEs: 1 to 31; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: BD; multicast precoding scheme: multicast precoding).


## 5.1. Analysis for the Mutual Coupling Effect

Figures2 and 3 illustrate the spectrum efficiency of linear antenna array and rectangular antenna array, respectively. The performance of MRT/BD precoding (here and after denoted as MRT and BD in the legends) and multicast beamforming (here and after denoted as multi in the legends) are evaluated with and without mutual coupling (here and after denoted as MC in the legends) based channel model. Since the mutual coupling effect introduces extra correlation between adjacent antennas, the channel becomes less independent in a spatial prospect hence the performance of all precoding and beamforming schemes are expected to be degraded. This is clearly observed in Figures 2 and 3 as the spectrum efficiency degrades for all precoding and beamforming schemes for both antenna array configurations. Also, since the mutual coupling effect of rectangular antenna array is considered to be more prominent than the linear antenna array, the performance degradation of rectangular antenna array is worse compared to linear antenna array as observed.Figure 2
Spectrum efficiency (linear antenna array; number of UEs: 16; +MC: mutual coupling; precoding schemes: BD/MRT/multicast).Figure 3
Spectrum efficiency (rectangular antenna array; number of UEs: 16; +MC: mutual coupling; precoding schemes: BD/MRT/multicast).


## 5.2. Comparison between Linear Precoding and Multicast Beamforming

To compare the linear precoding with multicast beamforming, we first choose the BD precoding for unicast FRs. Figures4 and 5 illustrate the spectrum efficiency comparison between BD precoding and multicast beamforming for linear antenna array and rectangular antenna array, respectively. The total spectrum efficiency for the entire system is illustrated and becomes legend as “all group” for observation convenience. For linear antenna array, the performance of multicast beamforming is worse than BD precoding. For rectangular antenna array, however, despite the fact that the performance of BD precoding and multicast beamforming both degrade, the multicast beamforming suffers much less degradation than BD precoding. Similar phenomena are observed in Figures 6 and 7 which illustrate the comparison between MRT precoding and multicast beamforming. This is because the multicasted information for each beamforming user is the same. With a higher correlation between antennas in rectangular antenna array, the signal leakage of multicast beamforming between each user is more of in-phase component and contributes less to the interference for the intended multicast user than that of the unicast group. Hence the multicast beamforming has better performance than linear precoding in terms of channel degradation introduced by concentrated antenna array configuration.Figure 4
Comparison between BD precoding and multicast beamforming (linear antenna array; number of unicast UEs: 8; number of multicast group UEs: 8; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding Scheme: BD; multicast precoding scheme: multicast precoding).Figure 5
Comparison between BD precoding and multicast beamforming (rectangular antenna array; number of unicast UEs: 8; number of multicast group UEs: 8; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: BD; multicast precoding scheme: multicast precoding).Figure 6
Comparison between MRT precoding and multicast beamforming (linear antenna array; number of unicast UEs: 8; number of multicast group UEs: 8; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: MRT; multicast precoding scheme: multicast precoding).Figure 7
Comparison between MRT precoding and multicast beamforming (rectangular antenna array; number of unicast UEs: 8; number of multicast group UEs: 8; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: MRT; multicast precoding scheme: multicast precoding).


## 5.3. Analysis for the Grouping Strategy

Figures8, 9, and 10 illustrate the trend of total spectrum efficiency with an increasing multicast user number while total user number is fixed to be 32. Specifically, a scenario with a much higher correlation between adjacent antennas for rectangular antenna array is presented in Figure 10. With an increasing number of multicast FRs, the total spectrum efficiency decreases since the multicast beamforming achieves worse performances than linear precoding in low correlation scenario. Meanwhile in Figure 10 we observe that, in high correlation scenario, the multicast beamforming outperforms linear precoding greatly and hence contributes to the growing trend of total spectrum efficiency when assigning more FRs to multicast groups, which conforms with the analysis in previous subsection. This is significant for the design of antenna array in different multicast service model to achieve best performance in practical wireless communication system.Figure 8
Grouping strategy analysis (linear antenna array; number of unicast UEs: 31 to 1; number of multicast group UEs: 1 to 31; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: BD; multicast precoding scheme: multicast precoding).Figure 9
Grouping strategy analysis (rectangular antenna array, low antenna relevance; number of unicast UEs: 31 to 1; number of multicast group UEs: 1 to 31; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: BD; multicast precoding scheme: multicast precoding).Figure 10
Grouping strategy analysis (rectangular antenna array, high antenna relevance; number of unicast UEs: 31 to 1; number of multicast group UEs: 1 to 31; intergroup interference cancellation: null space; +MC: mutual coupling; unicast precoding scheme: BD; multicast precoding scheme: multicast precoding).


## 6. Conclusion

This work investigates the wireless multicast transmission for PSN in 5G. The proposed hybrid system utilizes multicast beamforming transmission as well as the MU-MIMO linear precoding transmission. Null space method based interference cancellation is performed in the system. The simulation results indicate that under the mutual coupling channel model, the performances of different linear precoding schemes and multicast beamforming scheme degrade, while linear precoding suffers more degradation than multicast beamforming in different scenarios.



---

*Source: 4306356-2016-12-18.xml*