A Novel Joint Transmitting and Receiving Antenna Selection for Spatial Multiplexing Systems

With the progress of mobile technologies, iPADs, portable computers, and smart phones have been widely used in our life and work. People always hope that they can get easily information and access quickly to the Internet in anytime and anywhere, the momentum will be stronger in the future, it needs higher transmission data rate and better quality of service (QoS) of communication networks. The traditional mobile system installed single antenna at the transmitter-receiver, which is defined as SISO, has been far from meeting the actual needs. For the purpose, we erect multiple transmitting antennas and multiple receiving ones at both transmitter and receiver of mobile systems. Compared with SISO mobile communications, we need not increase signals bandwidth and their transmitting power, the channel capacity and QoS performance of the transmission links can be greatly enhanced by using antennas array at the transmitting side and receiving side [1,2,3]. In the past few years, the theoretical and applied investigation on MIMO technology has attracted the attention of experts, scholars and operators all over the world in the past 20 years. The 3rd Generation Partnership Project (3GPP) working group initially introduced MIMO technology to 3G, and the current 4G, 5G and the future 6G are full of application examples of large-scale MIMO technology [4].

From the perspective of the manufacturing cost of MIMO systems, antenna is usually composed of small or discrete monopole antenna and bipolar antenna, which has low cost and can be widely used. The transmitting and receiving radio frequency (RF) links are complex in structure, each channel including modulation and demodulation, amplifying circuit, digital to analog to digital conversion circuit (D/A/D), filter and other components. Therefore, the manufacture cost of multi-antenna wireless systems are high and the energy expenditure are large. The AS application is undoubtedly the first technology to comprehensively consideration both the high manufacturing complexity and good performance for the mobile MIMO systems [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28]. In terms of average capacity and bit error rate (BER), both theory and experiment have showed that the performance with AS system is much higher than one without AS under the same RF channels.

Although there is a good performance of the AS MIMO systems, but they are still a big conundrum, how to quickly select a set of antenna subsets from multiple antennas array, its computational complexity (CC) is very high. In the past two decades, many researchers have extensively studied how to reduce the CC of AS algorithm at the transmitter or receiver, a lot of relevant research literatures have been published [5,6,7,8,9,10,11,12,13,14]. Antenna selection is a milestone on the application of MIMO technologies, the early studies on the AS mainly focus on the transmitting side or the receiving side. There are a few research literatures on simultaneous joint AS at both transmitter and receiver (JTRAS). Especially under the background of large-scale MIMO technology application becoming development trend in the future, it will be a challenging task to study JTRAS algorithm. Research difficulties on JTRAS mainly have two aspects, first, considering the working band of millimeter wave, and the current RF channels of transmitter and receiver can achieve 4 or 8. When antenna array is large scale, no matter how many antennas are selected, if the optimal exhaustive search AS algorithm (ESAS) [15] is adopted, its CC is too high to have practical value. Second, the wireless mobile channel is rapidly changing, if the AS algorithm cannot quickly select the optimal antenna combination during the transmission of a symbol or a group of symbols, the communication system will be in a chaotic state. In the last decade, due to the above two reasons, the literatures on the research of JTRAS sub-optimal algorithm are relatively scattered. More research is how to apply MIMO technology to communication systems and how to combine it with other communication enhancement technologies. At present, the research direction of AS is to find the sub-optimal scheme as an alternative algorithm, which makes some performance concessions in exchange for the reduction of computational complexity [15,16,17,18,19]. The research on JTRAS sub-optimal algorithm mainly starts from two aspects, one is based on the simplification of matrix characteristics, the other is the parallel computing capability of evolutionary algorithm. Their common characteristic is to reduce the CC. The suboptimal scheme studied in this paper is a decreasing algorithm based on matrix characteristics. Gorokhov A et al. [15] proposed for the first time a joint AS method for transmitting and receiving antennas (JTRAS) of MIMO mobile systems, which is called exhaustive search AS (ESAS) in this paper. The ESAS algorithm requires to compute [CNRLr⋅CNTLt][C_{{N_R}}^{{L_r}} \cdot C_{{N_T}}^{{L_t}}] determinants with size L_r × L_t sub-matrix, because of its exhaustive search for the optimal antenna subset, its CC is very high. In order to filter out an approving antenna combination, it is necessary to further reduce the CC of antenna subset. For the reason, an AS scheme in literature [5] based on the channel maximum Frobenius 2 norm (F2 norm) over the row and column of the channel matrix (CHM) is proposed, here we define it as norm-based selection (NBS). To compare with the ESAS, the NBS has lower CC, but waste of its average capacity is higher. In literature [16], a separable transmitting and receiving selection strategy is proposed. This AS algorithm starts at the receiver side and the full CHM with N_R × N_T. One picks out the L_r antennas which contribute a lot to the capacity from the receiver, and they correspond to the L_r rows of the CHM. The same method is applied to select L_t antennas at the transmitter, corresponding to the L_t columns of the CHM. The CC of this algorithm is max {N_T, N_R}²L², it reduce largely the CC compared to the optimal ESAS algorithm [15], but loss of its capacity performance is larger than that of the ESAS.

To reduce the CC of the optimal ESAS algorithm and retain its capacity, some JTRAS algorithms were proposed in [17,18,19,20]. Blum R. S. et al. [17] proposed a greedy algorithm to combine AS at the transmitter and receiver. The algorithm starts with an empty selection matrix (SM), in each iterative step, we choose a pair of antennas with maximum capacity increment to add to the rows and columns of the SM until all antennas are selected. Therefore, the algorithm is also called the incremental JTRAS (IJTRAS). The capacity performance based on IJTRAS is very close to the optimal ESAS scheme and exceeds the other suboptimal AS schemes. By channel capacity formula derivation and resorting to the matrix inversion lemma, an efficient JTRAS (EJTRAS) algorithm was proposed [18]. On the basis of IJTRAS algorithm, assuming that L represents the number of selected antennas, the EJTRAS scheme can reduce the CC of the IJTRAS in the factor L and retain its capacity. Based on the thinking of literature [17] and [18], Chen C. E. [20] proposed a new calculable efficient JTRAS algorithm, which can effectively decrease the CC of the algorithm by modifying the AS's judgment criterion, and at the same time, there is no overhead in the statistical capacity. Developing the theory of the matrix partition, a concise joint AS judgment formula (CJAS) was presented [19], which reduces the AS's CC by a factor of 2L, where L is same definition in literature [18], and keeping the average capacity is near to the ESAS method. The authors in literature [21,22] believe that it is an inevitable development of mobile communication to have massive antennas on the base station. It is assumed that the mobile terminal of the cell is still equipped to a single antenna, that is, the point-to-point transmission between the base station and the mobile terminal is a MISO system or a SIMO system. This paper integrates all the mobile terminals in the service cell into a multi-antenna system and proposes an AS system with the largest transmitting and receiving capacity of the base station. AS always takes place at the base station side, and the selection criterion is the capacity optimization of the whole cell, rather than the capacity maximization of a single user. In order to choose quickly antenna subset, the concepts of square maximum volume and rectangular maximum volume are introduced to simplify the capacity calculation of MIMO system, which greatly reduces the computational complexity of the AS algorithm. However, the system is equivalent to MIMO single-ended AS and is not a JTRAS. Two JTRASs based on genetic algorithm (GA) have been presented in MIMO systems [23, 24]. Similar to GA, Particle swarm optimization (PSO) is an important evolutionary algorithm. By changing the position and speed of the particle population, the adaptive function is calculated to search the antenna subset. The average capacity obtained from the simulation experiment shows that PSO algorithm is a sub-optimal selection. Based on the maximum capacity selection criteria, several AS algorithms using binary PSO (BPSO) are presented in literature [25,26,27,28]. The simulating experiments show that these AS algorithms can get high capacity. They are especially suitable for massive MIMO system.

In the wireless mobile environment, to satisfy the different service requirements of various customers, it is necessary to ensure that the mobile communication equipments can transmit data stream at high speed at all times. Therefore, it is necessary to study the fast AS algorithm matching the random and rapid changes of the mobile environment. Especially when the NSA is large, such as more than half of the TNA, the CC of the several AS schemes analyzed above have no obvious advantages [15,16,17,18,19,20]. From what has been discussed above, applying the block theory of matrix, we propose a novel decreasing JTRAS (DJTRAS) algorithm to choose quickly antenna subset. The novel algorithm starts from a full MIMO channel matrix. Based on a simplified channel capacity formula, in every iterative period we delete repeatedly a row and a column of the equivalent decrement CHM with minimal capacity loss. The simulation results show that the new DJTRAS is better than NBS and CJAS. In terms of the CC, the numerical result curves show that it is better than CJAS. The proposed AS scheme achieves a good balance in both capacity performance and computational overhead.

In addition to part 1, the manuscript has five parts. We will explain MIMO channel model and capacity formula, discuss AS criteria in part 2. Part 3 presents a AS scheme with integrating transmitter antennas and receiver antennas. Part 4 analyses the complexity of the proposed selection algorithm. Part 5 presents the capacity performance curves calculated by computer simulation using the proposed AS algorithm, and some discussions and analyses on the results are given. Finally, some conclusions are generalized in Part 6. Several notations are often appeared in the following part of the paper: matrices and vectors are typed in boldface and italic, their difference is that uppercase letter is matrix and lower case letter denotes vector. I_k and det(·) represent a k × k unit matrix and a k × k determinant, respectively. A complex conjugate transpose is defined by the superscript (·)^H.

Fig. 1 plots a block diagram for MIMO system, it equips with transmitter antennas N_T and receiver antennas N_R, here for convenience, Let’s denote it {N_T, N_R}. It is assumed that the antenna spatial interval is large enough to meet the spatial multiplexing requirements of MIMO system. Furthermore, it is assumed that the MIMO system is only equipped with L_t and L_r RF circuits at the sending side and accepting side and satisfies the relation of L_t ≤ N_T and L_r ≤ N_R. Therefore, to meet the antenna and RF circuit matching, it is important that the MIMO mobile system performs AS operation, the selecting criterion of antenna is based on the maximum capacity. The MIMO fading channel is modeled as a slowly varying flat Rayleigh distribution, the MIMO channel matrix is expressed by H. The modulated signals are sent by N_T at the transmitter simultaneously and go through the same distribution of wireless fading channels. It is assumed that the receiver adopts coherent demodulation without frequency deviation and they are completely synchronized. In this way, these received RF signals of antenna N_R at the receiver are demodulated and the baseband signals after matching filtering are given by the following equation (1). (1)y=ρNTHS+Nnoi{\bf{y}} = \sqrt {{\rho \over {{N_T}}}} {\bf{HS}} + {{\bf{N}}_{noi}}

Fig. 1

Where S is the signals sent by each transmitting antenna after matching modulation, which is a vector S=[s1,s2.⋯,sNT]T{\bf{S}} = {[{s_1},{s_2}. \cdots,{s_{{N_T}}}]}^T with N_T dimension. Assume that these modulation signals are uniformly taken from constellation points both MPSK and MQAM and that they satisfy the constraint relation E[S^HS] = P. N_noi = [n₁,n₂.…,n_NR]^T denotes the noise vector corrupted in N_R receiving antennas, each noise element is subject to independent and identically Gaussian distribution, ni~ℕ(0,σn2)n_i \sim {\rm{\mathbb N}}(0,\sigma _n^2) . The channel matrix H for a {N_T, N_R} system is in equation (2). (2)H=(h11h12⋯h1NTh21h22⋯h2NT⋮⋮⋮⋮hNR1hNR2⋯hNRNT){\bf{H}} = \left({\matrix{{{h_{11}}} & {{h_{12}}} & \cdots & {{h_{1{N_T}}}} \cr {{h_{21}}} & {{h_{22}}} & \cdots & {{h_{2{N_T}}}} \cr \vdots & \vdots & \vdots & \vdots \cr {{h_{{N_R}1}}} & {{h_{{N_R}2}}} & \cdots & {{h_{{N_R}{N_T}}}} \cr}} \right)

Where h_i,j (i=1,2,...N_R; j=1,2,...N_T) represents the channel coefficient between the jth transmitting antenna and the ith receiving antenna, which is a random process, its amplitude follows Rayleigh distribution, where its average power takes 0.5, and the phase follows uniformly distributed at (0,2π). Assuming that each transmitting antenna transmits the same power each time, we can write the capacity formula of {N_T, N_R} system by equation (3a,3b) [6]. (3a)C=log2det(INR+ρNTHHH) NR≥NTC = \mathop {\log}\nolimits_2 \det ({{\bf{I}}_{{N_R}}} + {\rho \over {{N_T}}}{\bf{H}}{{\bf{H}}^H})\quad \quad {N_R} \ge {N_T}(3b)C=log2det(INT+ρNTHHH) NR≤NTC = \mathop {\log}\nolimits_2 \det ({{\bf{I}}_{{N_T}}} + {\rho \over {{N_T}}}{{\bf{H}}^H}{\bf{H}})\quad \quad {N_R} \le {N_T}

Where the variable ρ denotes the received Signal-to-Noise Ratio (SNR), assuming that we can estimate the channel state information (CSI) at the receiver, the transmitter obtains the selected antenna subset and partial CSIs in real time through a special feedback channel. Similarly, we assume that the actual RF channels of {N_T, N_R} system is L_t and L_r at the sending and receiving ends, respectively. In this way, we can modify equation (3) into the system capacity calculation equation (4) after antenna selection. (4)C(Hsel)=log2[det(ILr+ρLtHselHselH)]C({{\bf{H}}_{sel}}) = \mathop {\log}\nolimits_2 [\det ({{\bf{I}}_{{L_r}}} + {\rho \over {{L_t}}}{{\bf{H}}_{sel}}{\bf{H}}_{sel}^H)]

For a selected subset by particular AS algorithm, its channel matrix defined as H_sel. Equation (4) is used to calculate the capacity and compare the capacity results of all possible selected subsets. Finally, the antenna subset combination with the largest capacity is selected. For example, we see that the search number of subset combination is CNRLr×CNTLtC_{{N_R}}^{{L_r}} \times C_{{N_T}}^{{L_t}} for the ESAS scheme. Equation (5) gives a judgment formula of target adaptation function for any AS methods based on the decision of maximum capacity. (5)Hsel:argmax︸Hsel⊂H{log2[det(ILr+ρLtHselHselH)]}{{\bf{H}}_{sel}}:\underbrace {\arg \max}_{{{\bf{H}}_{sel}} \subset {\bf{H}}}\{\mathop {\log}\nolimits_2 [\det ({{\bf{I}}_{{L_r}}} + {\rho \over {{L_t}}}{{\bf{H}}_{sel}}{\bf{H}}_{sel}^H)]\}

MIMO system channel is modeled as a flat and slow Rayleigh fading, it will remain constant when every frame or block of signals is sent. But it will change from one frame cycle to the next, so the AS operation must be completed in a frame period. Therefore, the key of AS research is to find an efficient and fast AS algorithm, they comprehensively consider the speed, performance and complexity of the AS algorithm, which will be explored in this paper.

Based on the above analysis of the judgment formula for the novel JTRAS algorithm, the CC of the new AS scheme mainly depends on the operation two terms of equation (13), which are gxyBn−1gxyH{{\bf{g}}_{xy}}{\bf{B}}_n^{- 1}{\bf{g}}_{xy}^H and hxHDn−1hx{\bf{h}}_x^H{\bf{D}}_n^{- 1}{{\bf{h}}_x} . These two terms all involve multiplication of complex matrices and complex vectors and the inverse of a matrix. Therefore, to retard the CC of the new proposed AS scheme, we must try to simplify the two terms, so as to reduce the times of multiplication. For arbitrarily complex square matrix P and complex vector u, we assume that u is a (1 × t) row vector, and P is a (t × t) matrix. We define a quadratic form function (QFF) f (u, p) as following [30]: (14)f(u,P)=uPuH=∑i=1tpiiuiuj*+∑i=1,i≠jt∑j=1tpijuiuj*f({\bf{u}},{\bf{P}}) = {\bf{uP}}{{\bf{u}}^H} = \sum\limits_{i = 1}^t {p_{ii}}{u_i}u_j^* + \sum\limits_{i = 1,i \ne j}^t \sum\limits_{j = 1}^t {p_{ij}}{u_i}u_j^*

For a matrix and a vector, the superscript “H” represents their complex conjugate transpose, the superscript “*” denotes complex conjugate operation. Suppose that the symmetric element of the complex matrix P satisfies pij=pji*{p_{ij}} = p_{ji}^* , then pijuiuj*=(pji*ujui*)*{p_{ij}}{u_i}u_j^* = {(p_{ji}^*{u_j}u_i^*)}^* , i, j ∈ 1, …, t. Applying this relation, the second term on the right of equation (14) can be seen that when i > j, pijuiuj*{p_{ij}}{u_i}u_j^* is known, so it is not necessary to calculate the size of (pji*ujui*)*(p_{ji}^*{u_j}u_i^*)^* , so as to reduce the number of multiplications of equation (14). The computational complexity (CC) of the QFF uPu^H goes down to t²/2 complex multiplications.

As both Bn−1{\bf{B}}_n^{- 1} and Dn−1{\bf{D}}_n^{- 1} are conjugate symmetric matrices, the total CC of updating gxyBn−1gxyH{{\bf{g}}_{xy}}{\bf{B}}_n^{- 1}{\bf{g}}_{xy}^H and hxHDn−1hx{\bf{h}}_x^H{\bf{D}}_n^{- 1}{{\bf{h}}_x} can be reduced by half applying the quadratic form function. If Bn−1{\bf{B}}_n^{- 1} and Dn−1{\bf{D}}_n^{- 1} are given, the CC of the two QFF, hxHDn−1hx{\bf{h}}_x^H{\bf{D}}_n^{- 1}{{\bf{h}}_x} and gxyBn−1gxyH{{\bf{g}}_{xy}}{\bf{B}}_n^{- 1}{\bf{g}}_{xy}^H , is as follows: (15)υ=(NR−n)22+(NT−n−1)22\upsilon = {{{{({N_R} - n)}^2}} \over 2} + {{{{({N_T} - n - 1)}^2}} \over 2}

We do not consider the inverse operation of Bn−1{\bf{B}}_n^{- 1} and Dn−1{\bf{D}}_n^{- 1} , in the iterative process of the novel JTRAS algorithm, the algorithm will perform (N_T − n) times of different hxHhx{\bf{h}}_x^H{{\bf{h}}_x} calculation and (N_R − n) times of different gxygxyH{{\bf{g}}_{xy}}{\bf{g}}_{xy}^H calculation. Based on the calculation of the previous QFF, the calculation formula of the total CC of the new AS algorithm is as follows: (16)Vhg=∑n=0NT−Lt−1[(NT−n)(NR−n)22]+∑n=0NR−Lr−1[(NR−n)(NT−n−1)22]{V_{hg}} = \sum\limits_{n = 0}^{{N_T} - {L_t} - 1} \left[ {{{({N_T} - n)({N_R} - n{)^2}} \over 2}} \right] + \sum\limits_{n = 0}^{{N_R} - {L_r} - 1} \left[ {{{({N_R} - n)({N_T} - n - {{1)}^2}} \over 2}} \right]

Besides the CC of QFF, we should consider the inverse operation of Bn−1{\bf{B}}_n^{- 1} and Dn−1{\bf{D}}_n^{- 1} at each iterative step. Apply computational matrix theory [31], we can obtain that the CC of a matrix inversion is O(n³), where n is the invertible matrix’s rank. After considering the inverse calculation of Bn−1{\bf{B}}_n^{- 1} and Dn−1{\bf{D}}_n^{- 1} , the total CC of the new AS algorithm is: (17)VBD=∑n=0NR−Lr−1(NR−n)3+∑n=0NT−Lt−1(NT−n−1)3{V_{BD}} = \sum\limits_{n = 0}^{{N_R} - {L_r} - 1} {({N_R} - n)^3} + \sum\limits_{n = 0}^{{N_T} - {L_t} - 1} {({N_T} - n - 1)^3}

Therefore, the total CC for the proposed decreasing JTRAS scheme is: (18)νC=νhg+υBD{\nu _C} = {\nu _{{\rm{h}}g}} + {\upsilon _{BD}}

The equation (18) can be further written as: (19)νC=νhg+υBD=NTNR3+NRNT36−NTLt3+NRLr36+NR4+NT44−Lr4+Lt44\matrix{{{\nu _C}} \hfill & {= {\nu _{{\rm{h}}g}} + {\upsilon _{BD}}} \hfill \cr {} \hfill & {= {{{N_T}N_R^3 + {N_R}N_T^3} \over 6} - {{{N_T}L_t^3 + {N_R}L_r^3} \over 6} + {{N_R^4 + N_T^4} \over 4} - {{L_r^4 + L_t^4} \over 4}} \hfill \cr}

The derivation process of the last equation in the equation (19) can be found on the Appendix.

Fig. 2 and Fig. 3 show the curves of the complexity of AS algorithms as a function of the NSA L for 8 × 8 and 12 × 12 MIMO systems, respectively. It is convenient to comparison of the CC, we also have drawn curves of the complexity of both CJAS [19] and the optimal AS (OAS) algorithm. When the NSA L takes one, the complexity of the proposed decreasing JTRAS (DJTRAS) scheme is larger than OAS scheme. It is because that the novel DJTRAS algorithm needs to compute (N_T − 1) × (N_R − 1) times of the equation (12). While the OAS algorithm only need to do N_T × N_R times of the vector norm. With the increase of the NSA L, we can observe that the CC of our DJTRA scheme is far less than OAS method. The complexity burden of the proposed DJTRAS is larger than that of the CJAS when the NSA L is less than half of the TNA. While the NSA L more than half of the TNA, the complexity of the proposed algorithm is less than CJAS. For example, if N_T = N_R = N = 8, L_t = L_r = L = 6, the CJAS needs iterate 6 steps to obtain the set number of selected antennas, and the proposed DJTRAS scheme only needs do two iterative periods to obtain it. In this situation the proposed DJTRAS scheme only needs operate two times of matrix multiplication, while CJAS has to perform six times of matrix multiplication. Therefore, under the conditions of L_r = L_t = L and L > max(N_T, N_R)/2, the CC of the novel DJTRAS will be less than that of CJAS. It can be seen that our DJTRAS scheme has obvious advantages when the NSA L is larger.

Fig. 2

Fig. 3

If the transmitting end is configured with N_T antennas, we choose L_t, and the receiving end is configured with N_R antennas and choose L_r on a MIMO system, which defined as {N_T, L_t;N_R, L_r}. It is assumed that the wireless channels between the transmitting and receiving antennas of the MIMO system are mutually independent, and their amplitudes of the channel gain coefficient obey the Rayleigh fading, and their angles is uniformly distributed from 0 to 2π. In this part, we use the method of Monte Carlo simulation experiment to calculate the average capacity and its probability based on decreasing JTRAS (DJTRAS). In order to reflect the advantages of the new AS scheme, the capacity of four common AS algorithms are simulated in the capacity curve, such as rand, ESAS (here defined as the optimal selection algorithm, OS), NBS and CTAS, etc. [19]. According to the capacity comparison of these five different antenna selection algorithms, the paper gives two groups of curves of the {N_T, L_t;N_R, L_r} system under various antennas and RF channels. One group is the curve of the average channel capacity changing with the receiving Signal-to-Noise Ratio (SNR). The other is to use the cumulative distribution function (CDF) of capacity to calculate the statistical capacity. Regardless of which set of curves, these graphs reveal the capacity of the various AS algorithms and how their performance compares. All the simulating curves in this section are obtained from 1000 Monte Carlo simulation experiments as statistical data.

When doing the Monte Carlo experimental simulation, we assume that the SNR range of the receiver is 8dB to 20dB, and an experimental result is simulated every 1dB. From Fig. 4 to Fig. 7, the average capacity curves of {5,2;5,2}, {5,3;5,3}, {7,2;7,2} and {7,3;7,3} AS systems are drawn. In order to contrast the statistical capacity of the novel DJTRAS method, simulation results of other four AS schemes are presented in these figures. The simulation results of the four graphs all show that IEAS algorithm (marked AS “optimal” in the graphs) has the best average capacity, while random AS algorithm (marked AS “rand” in the graphs) has the worst performance. It is important that the statistical capacity of the new decreasing JTRAS is very close to that of the IEAS scheme. Similarly, compared with CJAS algorithm, the new DJTRAS has the advantage of low CC and slightly better average capacity. In addition, when the total antenna number (TAN) both N_T and N_R take 5 and 7, it can be viewed from Fig. 4, Fig. 5 to Fig. 6, Fig. 7 that the larger the number of selected antennas (NSA) L is, the better the median capacity will be. SNR=20dB, it can be observed that the median capacity of the proposed DJTRAS scheme is approximately 12bits/s/Hz and 18bits/s/Hz in figure 4 and figure 5, and approximately 12.5bits/s/Hz and 18.5bits/s/Hz in Fig. 6 and Fig. 7. It can be seen that when L increases from 2 to 3, the capacity will increase by 6bits/s/Hz. When we increase the TAN N_T and N_R, the NSA L is constant, as shown in Fig. 4 and Fig 6, L=2, N_T = N_R = 5 and 7, SNR=20dB, the increment of average capacity is less than 1bit/s/Hz. As shown in Fig. 5 and Fig. 7, L=3, N_T and N_R were 5 and 7 respectively, SNR=20dB, and the average capacity was also increased by less than 1bit/s/Hz. After the above discussion, we come to the conclusion that when the NSA L increases, the improvement degree of the average channel capacity of the AS system is much better than that of the TNA.

Fig. 4

Fig. 5

Fig. 6

Fig. 7

When the receiving SNR remains unchanged, the comparison curves of capacity CDF of various AS algorithms are given in Fig. 8 to Fig. 11. For {5,3;5,3} system, Fig. 8 and Fig. 9 draw the simulation curves of capacity CDF of various AS algorithms with SNR=8dB and SNR=20dB, respectively. In the case of {7,3;7,3} system, Fig. 10 and Fig. 11 also are respectively simulated to calculate the capacity CDF curves with SNR=8dB and SNR=20dB. From the four capacity CDF curves, we can calculate the median channel capacity of various AS algorithms by CDF=0.5. From this set of graphs, we can also get the same conclusions as Fig. 4 to Fig. 7. In short, the larger SNR is, the more the median capacity of AS system is. When L is constant, the more TNA there are, the greater the capacity will be. When the TNA is unchanged, the larger the NSA L is, the greater the median capacity is.

Fig. 8

Fig. 9

Fig. 10

Fig. 11

AS scheme is an effective method to decrease complex structure of the MIMO systems. In order to find a good compromise both the capacity and the CC for AS scheme, a decreasing joint transmitting and receiving AS (JTRAS) technique project is designed in this manuscript, this novel AS strategy is based on achieving maximum capacity for the MIMO systems. The curves of computer simulations on capacity’s indicator have revealed that the new JTRAS scheme can obtain sub-optimal performance. However, its complexity is far less than that of ESAS algorithm. The new JTRAS scheme is suited to large selected number. In addition, the spatial multiplexing of multi-antenna system can also provide a certain diversity gain, so as to it has great advantages in bit error rate (BER).