Abstract

Multiuser multiple input multiple output (MU-MIMO) wireless communication system provides substantial downlink throughput in millimeter wave (mmWave) communication by allowing multiple users to communicate at the same frequency and time slots. However, the design of the optimum beam-vector for each user to minimise interference from other users is challenging. In this paper, based on the concept of signal-to-leakage plus noise ratio (SLNR), we analyze the ergodic sum-rate capacity using statistical Eigen-mode (SE) and zero-forcing (ZF) models with Ricean fading channel. In the analysis, the orthogonality of channel vectors between users is assumed to guarantee interference cancelation from other cochannel users. The impact of the number of antenna elements on the achievable sum-rate capacity obtained by dirty paper coding (DPC) method considered as a nonlinear scheme for approximating average system capacity is studied. A power iterative precoding scheme that iteratively finds the most dominant eigenvector (optimum weight vector) for minimising cochannel interference (CCI), that is, maximising the SLNR for all users simultaneously, is designed resulting in enhancement of average system capacity. The average system capacities achieved by the proposed power iterative technique in this study compared with the singular value decomposition (SVD) method are in the ranges of 5–11 bps/Hz and 1–6 bps/Hz, respectively. Therefore, the proposed power iterative method achieves higher performance than the SVD regarding achievable sum-rate capacity.

1. Introduction

Millimeter wave (mmWave) communication which explores shorter propagation distance in frequency band of 20–40 GHz is a key enabler for the fifth generation (5G) mobile communication systems [1]. The mmWave communication also provides significant benefits to a variety of applications such as vehicular communication, wire-able networks, and autonomous robots [2]. In downlink transmission for a MU-MIMO system, a base station (BS) serves multiple users simultaneously in the same frequency and time slot. Thus the throughput can be enhanced by spatial multiplexing. However, the cochannel interference (CCI) becomes a dominant factor in capacity due to nonorthogonal signalling [3]. Hence, minimising CCI as much as possible at the end users is of importance. CCI can be suppressed by using linear precoders and decoders at transmitter and receiver sides [4]. Additionally, the channel state information (CSI) as well as partial information on the transmitter side can be used to improve the system’s performance [5]. For all users served by the BS, the CSI should be known in advance at the BS to support the CCI minimisation which is not available at user’s end. Explicitly, the overhead on the system is reduced in this approach because the channel information feedback is not required [6, 7]. To perfectly cancel CCI at each end user, the restriction on the system configuration is necessary. For instance, the number of antenna elements at the BS should be larger than the total number of antennas at the end users [7]. This assumption is usually valid for single antenna per user which is generally adopted in the MU-MIMO system.

There are several studies that have been undertaken in MU-MIMO downlink systems on how to minimise the CCI problem related to CSI and analyze the achievable sum-rate per user as well as the average system capacity. In [8], an in-depth capacity analysis for nonorthogonal multiple access (NOMA) mmWave massive MIMO systems was provided. A simplified mmWave channel model was also explored by extending the uniform random single-path (UR-SP) model with the angle-of-arrival (AOA). Furthermore, the capacity analysis was divided into high and low signal to noise ratio (SNR) regimes, where the dominant factors for the signal to interference plus noise ratio (SINR) were determined as interference and noise. A multicast beamforming approach was proposed in [9], where users located within the proximity of the BS can receive two different data streams simultaneously, while those away can receive only one data stream. In [7], a model was designed to transmit beamforming vectors to maximise the signal to leakage ratio (SLR), minimising transmission power which may cause interference to other users. The model is not restricted to a certain number of transmission antenna elements and can be extended to more general scenarios. An asymptotic deterministic SLNR optimisation approach for regularized zero-forcing (RZF) considering perfect CSI and antenna correlation was proposed [10]. It was found that when the users are homogeneously distributed and the number of antenna elements is large enough , the SLNR is asymptotically equal to SINR. In [11], an efficient statistical Eigen-mode space division multiple access (SE-SDMA) scheme for downlink ergodic sum-rate analysis based on SLNR was proposed. To maximise the approximate ergodic sum-rate capacity, authors in [12] designed an optimal beamforming vector for each user based on a three-dimensional beamforming algorithm. The authors in [3] focused on imperfect CSI scenario and developed a robust SLNR approach for compensating performance degradation caused by random CSI errors. In [13], a precoding scheme for heterogeneous networks (Het-Nets) was proposed based on SLNR under imperfect CSI scenario. In this case, regularization parameter was used to “weight” the precoding information for other cells (BSs).

In this paper based on the SLNR concept and Ricean fading channel model with perfect CSI in MU-MIMO systems, the statistical Eigen-mode (SE) and zero-forcing (ZF) models are derived. In addition, the average achievable sum-rate by each user and the overall average system capacity with different number of users are analyzed. Besides, the expression of dirty paper coding (DPC) method as a function of number of users and number of antenna elements is addressed. Moreover, the impact of on the average achievable sum-rate is demonstrated. The Ricean fading channel with the line-of-sight (LOS) and scattering components as channel model are adopted. We find that the ergodic capacity of DPC is approximately free from channel matrix which mainly depends on the number of BS antenna elements and average transmission power.

Furthermore, by exploring the leakage signal criteria, we propose a new design solution for precoding based on the proposed power iteration technique. The technique finds the optimum weight vector which maximises the SLNR, that is, enhances the system capacity by obtaining the dominant eigenvector for minimising the CCI. Finally, we compare the proposed method with the conventional solutions such as SVD. The numerical results demonstrate that significant throughput performance can be achieved with the proposed technique.

A comparative study of the proposed and SVD methods regarding computation complexity and storage is also carried out. The simplicity of the proposed algorithm is observed in approximating only one eigenvalue of a matrix in a sequence which is considerably more efficient as the number of iterations increase. In comparison, the SVD has to calculate all eigenvectors and pick the one with the most significant eigenvalue (the maximum eigenvalue), which requires larger computation capability and storage capacity.

The contributions of this paper are in twofold:

The SE and ZF beam former models in Ricean fading channel are examined based on SLNR instead of the commonly used SINR and orthogonal condition of the beam weight vector in MU-MIMO systems. The linear achievable sum-rate capacity is investigated using these models to validate the proposed method. Moreover, in nonlinear achievable sum-rate capacity, the DPC technique regarding the impact of the number of users and antenna elements on the ergodic capacity is addressed. The results show that the number of antenna elements has a significant contribution to the ergodic capacity, and the DPC sum-rate capacity logarithmically increases with the number of users.

A new optimum weight vector is developed based on the proposed power iteration method, which allows each user to maximise the SLNR and minimise the CCI from other users. The proposed method compared with the conventional SVD method shows that the proposed method can achieve better performance with relatively low cost.

This paper is organized as follows. In Section 2, the MU-MIMO downlink system is introduced. The beamforming and leakage signal approaches are presented in Section 3. In Section 4, the results and discussion of the study are presented. Finally, the conclusion is drawn in Section 5. The significant proofs of the algorithm are addressed in the appendix.

In this paper, the superscript , and denote the conjugate-transpose (Hermitian), transpose, and conjugate, respectively. , , and are the expectation operator, square identity matrix of size , and Gaussian random complex numbers, respectively.

2. System Model

Consider the MU-MIMO system model shown in Figure 1. The BS is equipped with uniform linear transmission antennas that serve users simultaneously. Assume that each user is equipped with a single antenna element [3]. Additionally, the number of the BS antenna elements is assumed to be larger than that of users ; and the power is equally allocated among all users [11, 12].

2.1. Signal Model

From the previous assumptions, the received signal at th user is given byThe expression in (1) can be rewritten in a more simplified form aswhere is the flat Ricean downlink channel vector between the BS and the th user. The weight vector is unit normalised beamforming vector of user and satisfies . Furthermore, the vector is the transmitted data symbol of user with . is unit normalised complex additive white Gaussian noise (AWGN). Moreover, is the transmission power of the th user under the constraint . Since the total transmission power is divided among all the users equally, the average power can be written as for each user [3, 11, 12]. Consequently, the signal magnitude of each user can be written as while the interference of the th user is determined by the sum of interference from all other users. Hence, the interference can be denoted as . Thus the SINR of the th user can be expressed as in (3) [12, 14]:

(Note: for simplicity we drop the “subscript ” from in the rest of the paper).The corresponding average achievable data rate for user can be written asBy substituting (3) into (4), the average achievable sum rate per user is given as in (5) [15, 16].Subsequently, the overall system ergodic sum-rate capacity can be written in bps/Hz as

2.2. Channel Model

In this subsection, we describe the channel model as in [11, 12]. In the most recent studies on mmWave communication, the Rayleigh fading channel model is employed which faces difficulty in capturing the fading variation in the presence of the LOS component. Moreover, mmWave MU massive antenna system can be used mainly in 5G hotpot scenarios in the LOS case. To overcome the capturing problem in Rayleigh fading channel, the Ricean fading channel is applied in this paper. Considering the Ricean fading channel model, the channel vector has two components: LOS and Rayleigh random distribution components. Thus, the channel vector can be expressed as [11]where is the power ratio between the LOS component and scattering components, also called Ricean -factor of the th user. Moreover, the channel vector is the random component, and the entries are independent and identical distribution complex Gaussian random variables. Vector is the deterministic component of the channel mean vector of the th user. Based on the assumption that the BS is equipped with uniform linear array antennas, the deterministic channel vector for each user can be written as [17, 18]where is the space between two antenna elements, λ is the wavelength of the signal, and is the angle of departure (AoD) for the th user, .

3. Beamforming and Leakage Signals

It is assumed that all users have a perfect instantaneous knowledge of their channel vectors while the BS knows all users channel vector including the channel mean vector. Additionally, assuming users, it is difficult to maximise the average system capacity due to a couple of calculations (capacity of each user) required. Therefore, it is hard to use SINR directly to obtain the optimum beamforming vector . To overcome this problem, the SLNR can be easily implemented to control the leakage signals from a specific user to other users as shown in Figure 1. The SLNR of the th user can be written as in (9) by assuming unit normalised Gaussian noise [4, 16]:Equations (3) and (9) have different denominators to define SINR and SLNR. In (3), the user channel vector is used with other users’ beamforming vectors to calculate SLNR of the th user. While in (9) the th user beamforming vector is used with other users’ channel vectors to calculate SLNR of the th user. Besides, the lower bound (LB) on the average SLNR of the th user can be represented as demonstrated in [11, 16]:In (10), we can have :where is a Hermitian matrix (channel correlation matrix) which can be constructed as in (12), where , and is calculated in (8). Hence, of the th user can be decomposed regarding the deterministic component of channel vector after normalisation by a factor [6]. Therefore, we can use the normalised version in the decomposition of the channel correlation matrix as inwhere is the orthogonal subspace of . Then we can obtain the optimum beamforming vector which maximises the lower bound of the SLNR as follows:On the other hand, to minimise the denominator of (9) we havewhere vector is orthogonal to . Accordingly, from (12) we haveThe maximum value of the lower bound can be achieved if is orthogonal to . By substituting (15) in (10), we have In (16), the lower bound is directly affected by the number of antenna elements , the Ricean factor , and the number of users . Assuming that the number of users is increased as , taking and , the lower bounds of SLNR are expressed approximately as and , respectively. Therefore, it is seen that the lower bound is independent of when .

3.1. Analysis of Linear Achievable Rates by SE and ZF

In this section, we evaluate the Ergodic sum-rate capacity of SE obtained by (13) and (14) and assume that ZF capacity has perfect CSI when calculating the mean gap loss between and . Based on the orthogonal beamforming condition in (13) and (14), the achievable sum-rate capacity obtained by the th user is given byThe achievable sum-rate capacity by ZF is expressed aswhere is the unit normalised beamforming vector, which is selected as the th column of the normalised matrix with to suppress the interference from all other users for , if . Accordingly, the mean gap loss with perfect instantaneous CSI is given byFrom the orthogonal condition which is addressed in (14), we neglect the interference from the other users to the th user regarding the signal component. In this case the sum-rate capacity used in (14) changes to the loose bound sum-rate capacity of the th user expressed as [11]Using (18) and (20), we rewrite (19) as the mean gap loss in (21) as

3.2. Nonlinear Achievable Rate Analysis by DPC

The BS can serve all users simultaneously and achieve maximum system capacity as much as possible as the CSI is known at the BS. In MIMO systems, the system capacity obtained by downlink strategy is called DPC. Practically, it is challenging to implement DPC because the encoding and decoding have high computations which is ineffective for large number of users. The achievable sum-rate capacity of the DPC can be written as [19]When the number of users is large, the DPC achievable sum-rate capacity is approximately expressed as [19, 20]

From (23), we can observe that the sum-rate capacity increases linearly with the number of antenna elements. However, the approximate achievable sum-rate capacity is a nonlinear function of . Therefore, the DPC is regarded as a nonlinear model.

3.3. Optimum Weight Vector Formulation

In this subsection, we find the solution for optimum weight vector , in which the maximum SLNR can be achieved to guarantee that all users have the ability to access the limited resources. For equal transmission of power to all users in a given case, we need to design to maximise the SLNR for each user at the same time minimise interference. Hence, the optimum weight vector can be written as follows:Equation (24) is subjected toConsequently, the SLNR is written asEquation (27) is further written aswhere is congregate interfering channel transfer matrix (CICTM) of the th user, which is an extended channel matrix that excludes only for the th user. Then is written aswhere and . The optimum weight code vector can be derived from (28) as follows [21]:where is square identity matrix.

The optimum weight code vector in (30) is called SVD solution, which is determined from the matrix in (29). Sometimes, when calculating the maximum eigenvalues for the given matrix, the dominant eigenvector does not correspond to the optimum code weight vector. Hence, the eigenvalues are put in ascending or descending order; then an optimum vector corresponding to the most significant eigenvalue is easily selected but at the expense of extra overhead to the system [4, 7]. To enhance the system performance and minimise the cost, a new approximation method expected to improve the average user achievable sum-rate capacity is proposed.

3.4. Proposed Optimum Weight Vector Design

A new method is developed based on power iteration method to find the optimum weight vector. The method requires only one weight vector during implementation. Equation (30) is modified and adopted as the main matrix in the proposed power approximation method as follows:

3.4.1. Power Approximation Method

The power approximation method generates a sequence of vectors , where is a nonzero initial selected vector. By normalising these sequence vectors under conditions stated in Section 3.4.3, the vector converges to the dominant eigenvector corresponding to the most significant eigenvalue. The normalisation is used to ensure that the most significant component of the given iteration is equal to one [22]. The steps in the proposed power iteration method are given in Algorithm 1.

Step 1: Initial Inputs. number of users ;
number of iteration ;
initial non-zero vector ;
Tolerance .
Step 2: For  : //to construct the CICTM for each user;
Calculate and as defined in Section 2.2;
Calculate for as in (7);
Construct the CICTM for each user as in (29);
End
Step 3: For : //to compute for each user.
Calculate matrix for th user as in (31)
Step 4: While  ;
set: , and ; (normalisation)
“the value of is the scaled version of
if , then ;
else, set ;
end
Step 5: calculate the dominant eigenvector ;
the dominant eigenvalue ;
End;
Step 6: SLNR maximise weight vector ;
End;
Step 7: Output: the optimum weight vector for each user.

This algorithm is summarised in a sequence of iterations as follows:where is the component (magnitude) of vector . Before proceeding to the next iteration, it is necessary to scale down (normalised) the sequence vectors in each approximation so as to keep the largest component at unity.

3.4.2. Convergence of the Algorithm

The matrix is square matrix with eigenvalues , which are in descending order as . With the initial vector chosen, the sequences and are recursively generated by the relations and , where and . These sequences converge to the dominant eigenvector and eigenvalue , respectively, as and . The proof of convergence is referred to in Appendix A.

3.4.3. The Speed of Convergence

Referring to (A.6) in Appendix A, we observe that the coefficient of in the sequence goes to zero, which is proportional to ; the convergence speed of sequence to is governed by . Therefore, the rate of convergence and the convergence of the constants to are linear. For any linear convergent sequence , we can use the Aitken technique which is used to make the linearly convergent sequences fast. The new convergence of the sequence can be written asThe convergence property of the algorithm is given by the following.

Theorem 1. Assume that there is one and only one eigenvalue of , and is semisimple; then either the initial vector which has no component in the invariant subspace associated with or the sequence of vectors generated by the algorithm converges to the eigenvector associated with and and converges to . The proof is shown in Appendix B.

Definition 2. Eigenvalue of matrix is semisimple with the degree of if it has a geometric multiplicity one and algebraic multiplicity . The vector is normalised by a specific scalar to make the most significant component of the vector unity. The initial vector is decomposed as follows:where , are the spectral projectors associated with the eigenvalues .

From the formula  , the power can be written as , where is the identity matrix and is the diagonal matrix. Consequently, we haveReferring to the Definition 2, because is a semisimple eigenvalue. Thus, we obtainThe spectral radius of each operator and means that the th power will converge to zero. Theorem 1 is true when . However, if , converge to normalise so that the most significant component is one. Meanwhile, the scalar converges to the eigenvalue , which is an immediate consequence in the formThus, the sequence of vectors is proven to converge to the optimum weight vector which maximises the SLNR in (28). The optimum weight vector is equal to the vector :The result in (38), which is the main objective of this paper, maximises the th user SLNR when substituted in (28) to find .

3.5. Cost Analysis

In this subsection, the cost function is analyzed regarding storage space required for the proposed power iteration technique and SVD method. Referring to (30), the value results in a new square matrix with dimension , where denotes the number of antenna elements. Therefore, finding for each user by using the solution of SVD in (30), we need to find at least two matrices, left eigenvectors matrix, right eigenvectors matrix, and diagonal eigenvalues from the original matrix. On the other hand, the required storage memory for the th user is equal to the size of three matrices in addition to the original square matrix. Hence, the total storage becomes (full SVD), where is the number of users. Then, the determination of the optimum weight vector can be from the right eigenvector or left eigenvector based on the corresponding maximum eigenvalue. Furthermore, we can reduce the size by taking only one side eigenvectors matrix (left or right) to save redundant storage size. Thus the reduced new storage capacity can be written as (reduced SVD).

Similarly, only the original square matrix with dimension and a vector with size is required in the proposed power approximation method. In this case, there is no need for storing the previous iteration result as it can be overwritten up to the last iteration or reaches the determined tolerance. Thus, the required storage for the power iteration technique is which is much less than full and reduced SVDs.

Figure 2 shows the storage capacity required in the proposed power iteration technique and SVD method. From Figure 2, it is observed that the number of antenna elements has a significant effect on the cost function. However, the proposed method has much less storage space compared to SVD hence recommended for a massive MIMO with large antennas.

4. Simulation Results and Discussion

In this section, the numerical results are presented with the following assumptions: the equal power allocation strategy and noise effect for all users are the same [23]. Additionally, the AoD for each user is in a horizontal direction with the uniform distribution in the range . Furthermore, for each user the Ricean -factor has uniform distribution in the range of , where  dB and  dB. Moreover, the number of users is increased up to . For SE and ZF analysis , 32, 64, and 128 are selected as the numbers of transmission antenna elements. The ergodic sum-rate capacities of SE and ZF are obtained in Figures 3 and 4. In Figure 3. It is observed that the average system capacity of SE is, respectively, directly and inversely proportional to the number of antenna elements and number of users. Moreover, when the number of users is small, the impact of on the SE capacity is more significant and conversely when the numbers are large, the impact is less significant. The SE achieves average capacities approximately in the range of 1.6–6.2 bps/Hz. At maximum number of users with , the corresponding SE sum-rate capacity is less than 1.6 bps/Hz. This is because each user faces CCI interference which is considered as sum of leakage signals from all other cochannel users. Referring to (17), it is validated that as the number of users increases, the overall contributed interference to the th user also increases. Therefore, the increase in the denominator in (17) causes a reduction in the sum-rate capacity per user resulting in an average lower system capacity.

In Figure 4, it is shown that the ergodic sum-rate capacities obtained by ZF are directly proportional to the number of antenna elements. By referring to (18) the interference from the other users is perfectly suppressed by the orthogonal beamforming for user channel and weight vectors. These vectors are orthogonal to each other based on the condition in which .

Moreover, the number of antenna elements has much contribution on the average system capacity with ZF than SE as comparatively shown in Figures 4 and 3. The maximum achievable capacities are 105 bps/Hz and 6.2 bps/Hz for ZF and SE, respectively, because of the interference effects on SE.

Figure 5 shows the upper bound of the mean-rate capacity gap loss ) between ZF and SE. Since the number of users tends to be large, the difference in the mean rate capacity is most likely the ZF capacity, because the capacity obtained by ZF is close to 15 times that obtained by SE, which precisely affects the mean-gap loss. Nevertheless, when the number of users increases, the difference tends to be very big. It is explicit that when the number of users is large under any number of antenna elements, the achievable sum-rate capacity obtained by SE almost tends to be zero. However, the ergodic sum-rate capacity achieved by ZF positively increases with increase in the number of antenna elements and tends to a fixed level with large number of users. Meanwhile, the achievable sum-rate capacity is not profoundly affected by the number of users.

As shown in Figures 35, the simulation results depict that the ZF seems ideal due to a perfect cancellation of CCI. However, in practice SE is dominant and realisable because it is challenging to ensure that the interference is perfectly cancelled by the orthogonal condition between the channel vector and optimum weight vector. Thus, carefully focusing the beams alignment between the BS and the users a high system capacity can be achieved.

In Figure 6, the antenna elements , 32, 64, and 128 are considered with number of users up to 32 to satisfy and show the effects of . The results show that the ergodic sum-rate capacity achieved by DPC is nonlinearly increased with the number of users. Moreover, if the number of users is fixed to 25, at the number of antenna elements 16, 32, 64, and 128, respectively, the corresponding ergodic sum-rate capacities are 80, 105, 175, and 250 bps/Hz, respectively. As the number of users is increased, the system capacity also increases. Thus, comparing the DPC with ZF and SE achieves much higher capacity with large number of antenna elements and users; hence it can provide multiuser diversity gain.

The results in Figures 79 are used to comparatively evaluate the proposed power iteration method in (38) and the SVD method in (30). In Figure 7, it is observed that the SLNR in the range −10.5 dB to 37 dB achieves the system capacities in the range 5 bps/Hz to 11 bps/Hz for the proposed method. On the other hand, the SVD method has the SLNR in the range −15 dB to 31 dB that achieves system capacities in the range 1 bps/Hz to 6 bps/Hz. It is clear that the proposed method has the average throughput nearly two times that of SVD. The main reason is due to the efficiency of the weight vector obtained in the proposed method that is capable of minimising the CCI much more than in the SVD. It is expected that the system capacity can be enhanced in case the BS has more antenna elements but at a higher simulation time. Additionally, the SVD method has a drawback in finding the exact eigenvector corresponding to the dominant eigenvalue in which the eigenvalues are not always in an orderly way (ascending or descending), which limits the SVD performance. Figures 8 and 9 are the Cumulative Distribution Function (CDF) and Probability Density Function (PDF) of SLNR of the proposed and SVD methods, respectively. It is observed that the proposed method achieves a much better SLNRs than the SVD method.

5. Conclusion

In this paper, based on MU-MIMO with massive antennas, the ergodic capacity for linear sum-rate analysis including SE and ZF beamforming based on SLNR technique is investigated. In addition, the nonlinear sum-rate analysis using DPC in Ricean fading channels based on SLNR was undertaken. A new method is proposed to find an optimum beam weight vector by exploring the power iteration method using eigenvector approximation. The number of antenna elements is selected up to 128 dramatically increasing the system capacity. By way of simulation, the most significant dominant eigenvector to maximise SLNR as well as minimise the CCI is obtained. By comparing with the SVD method, the proposed method achieves higher performance regarding mean achievable sum-rate capacity per user. The throughput of the proposed method is in the range of 5 bps/Hz–11 bps/Hz, while that of SVD in the range of 1 bps/Hz–6 bps/Hz. Therefore, the proposed method can provide significant system capacity enhancement. In future work, the proposed power iteration technique is recommended for a 5G MU massive antenna system.

Appendix

A. Proof of the Algorithm

As we know, matrix has eigenvalues and corresponding eigenvectors which are linearly independent and normalised and form a basis for dimensional space. Thus, let the initial vector be written as the linear combination form asAssume vector is chosen in such a way that , and the elements of are scaled so that . As the vectors represent the eigenvectors of the matrix , the multiplication is followed by normalisation as follows:Moreover, also we haveAfter iterations we arrive atFurthermore, we haveSince we assumed that for , then we haveHence it follows thatWe need both vectors and to be normalised and their most significant component is 1.

B. Proof of the Convergence

For the sequence of vectors which converges to the dominant eigenvector:By substituting with , it yieldsBy merging both into (A.8), we getTherefore, the sequences of the constants converge to the dominant eigenvalue as

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is supported by Key program of Beijing Municipal Natural Science Foundation with no. 17L20052 and also supported by National Key Laboratory of Electromagnetic Environment, China Research Institute of Radiowave Propagation under Grant no. 201600012, and the National Nature Science Foundation of China (NSFC) under Grant no. 61771194.