A Novel Linear Antenna Synthesis for Linear Dispersion Codes Based on an Innovative HYBRID Genetic Algorithm

As far as taking-away of the symmetry constraints is concerned, as a scientifically symmetry problem, the global synthesis for antenna arrays that produce the desired radiation pattern is also a highly nonlinear optimization issue in fact. Besides this, the built criteria offer the reasonable power patterns. The consequent synthesis could be implemented by looking for a nominal pattern. When the criteria are already sufficient, it can simply do the whole synthesis process. To utilize multiple antennae, a method to choose a transmit antenna for the linear dispersion codes (LDC-TAS) is implemented in this paper. The authors used the max–min-postsignal to noise ratio (SNR) criteria to select these optimal transmitting antennae while this dependent, linear receiver is applied to the varying and slow channel. The simulated results illustrate that this max–min-post-SNR criterion outperforms the Bell Labs layered space time transmitting antenna selection (BLAST-TAS) applying the same spectral efficiency than space–time block codes (STBC)-TAS in the environment with low SNR. Furthermore, once the M antennae are selected under the selection criteria, a max–min-post-SNR rule, a novel linear antenna synthesis to linear dispersion codes on the basis of an innovative HYBRID (of mixed characters or solutions) genetic algorithm has been presented and evaluated to formulate and address the optimal problem to non-uniformly spaced and linear arrays. The restricted side-lobes level, the main-lobe width, and the shaped beam pattern are contemporarily concerned via maximizing a pretty suitable cost function through the innovational advanced genetic-algorithm-based algorithm. The method proposed in this paper can provide flexibility and a simple insertion of the a priori knowledge under a small computing pressure. At last, a computing simulation is completed well and the results are shown. It should be noticed that some extensions of the presented method could also be easily utilized without an obvious increase in the algorithm complexity.


Introduction
An array of the antennae can offer an efficient way for processing and detecting the signals coming from some different orientations.Contrasted to a single antenna, a specific array of sensors is able to have its own beam pattern modified through the amplitude and phase distribution named the weights of this array.In order to achieve an array beam pattern of the antennae, the signals could be power weighted and accumulated once the pre-process on the antenna outputs is over.The problem on the antenna array pattern synthesis includes how to get the weights to satisfy a series of specifications over the beam pattern.
In this paper, we make the purpose of proposing a specific modular algorithm on the basis of the genetic algorithms (GAs) for various restrictions, just like the array thinned, the linear dimension minimization, the side-lobes peak minimization, BP (beam pattern) shape modeling, etc.Some completely experimental instances show that the method presented can work well.The presented GAs-based processing method is suitable and flexible to the antenna array.
Contrasted to the traditional SISO system (single-input single-output system), the MIMO system (multiple-input multiple-output system) can provide a huge capacity gain [1,2].Assuming that only the receiver has known the channel state information, the MIMO system has two modulation methods: Diversity and multiplexing [3,4].The orthogonal space-time block codes do not achieve the overall channel capacity in the MIMO channels in spite of its advantage of maximizing the diversity.Scientists Hassibi and Hochwald [5] firstly presented the revolutionary linear dispersion codes (LDCs).The researcher Robert W. Heath [6] proposed a family of LDC designs on the basis of the frame theory, which certifies good performance in terms of the diversity advantage.While the MIMO system can increase system capacity and assure communication quality, but an introduction of multiple RF (Radio Frequency) links requires additional requirements according to volume, power consumption, and costs of the hardware.Therefore, antenna selection becomes a cheap and low-complexity application method [5].Space-time coding schemes, involved in STTC (space-time trellis codes) [6,7] and STBC (space-time block codes) [8,9], have also been studied.Reference [10] proposes a BLAST-TAS (Bell Labs layered space time transmitting antenna selection) system for zero-force (ZF) receivers.
The multiple antenna system can improve the reliability and increase the capacity in the mobile communication.However, the RF chains are expensive depending on size, the hardware and the power.Consequently, the antenna selection intends to be the lower cost and lower complexity alternative scheme to catch a lot of the merits of the MIMO systems [11].In order to apply multiple antennae, some articles [11][12][13][14][15] have already come out lately in the research report in which the conception over the antenna selection was proposed for the space-time codes.Other authors [16,17] proposed a complete performance analysis on orthogonal STBCs using the antenna selection.The quasi-orthogonal STBCs were researched in some papers [18,19].
In order to increase the transmission rate and improve the diversity gain of the system, the article used a transmit antenna selection [11,12] scheme based on linear dispersion code (LDC-TAS).This system scheme adopts the minimum posterior signal to noise ratio (SNR) criterion, which can be applied to any configuration of the number of antennae at both transmitters and receivers, so as to provide diversity capability [13,14] under the condition of maximizing mutual information between transmitters and receivers.The computer simulation results illustrate that the minimum posterior SNR criterion outperform the other criterions under the different modulation ways just like 8PSK, 16QAM, and QPSK.Additionally, Contrasted to BLAST-TAS [15], this presented method still has significant diversity gain at the same spectral efficiency as well as low signal-to-noise ratio.
The global combining for antenna arrays that produce the desired radiation pattern [16][17][18] is also a nonlinear optimization issue.A number of mathematical algorithms have already been presented on this solution.Those examples of the analytical technologies including the Chebishev [19] method and Taylor algorithm [20], which are both well-known.Nevertheless, the calculus or analytical ways commonly used cannot optimize the power weights and positions of the array elements at the same time.To treat huge nonlinear searching spaces and enlarge the research to include an efficient placement of elements, stochastic methods are essential.
The problem-oriented genetic algorithms (GAs) have been extended to utilize some different experimental instances [21][22][23].As for the symmetrical array synthesis, the real-coded GA-based process was presented in the article [24] to optimize the array power weights while sensors are λ/2-spaced.Moreover, the binary genetic method [25] was utilized to address the problem of isophorical array thinning, and the stochastic approach was proposed [26] to optimize isophorical arrays with a certain number of the elements.The methods above consider symmetric arrays to reduce the calculating time.However, then this arrays optimizing [27,28] have already been researched via considering one higher number of degrees of freedom just like the elements position as well as weighting.From then on, the optimization algorithm based on a simulation annealing procedure was used in order to optimize both the positions and power weights, and in the article [29,30] to the beam pattern shaping and the global array combinations [31][32][33].
The arrangement of the article is organized as the following parts.After authors offer the specific background of this paper in advance.Firstly, the linear dispersion code is reviewed in Section 2. Secondly, a model of the system over the LDC-TAS is introduced, then some assumptions are also offered inside Section 3, and the full description about the antenna selection methods for LDCs based on ZF receiver is shown in detail in Section 4. Thirdly, in Section 5, this paper focuses on proposing that the modular algorithm, according to the genetic algorithm, can combine the linear real-weighted arrays based on different restrictions, including arrays thinning, the side-lobe peak minimization, the beam pattern (BP) shape modeling, and the minimization of the linear dimension.Then, Section 6 includes the numerical computing simulation results and the complete analysis to those results.Finally, Section 7 presents our conclusions.

Linear Dispersion Code
The linear discrete code was first proposed by [34][35][36], which maximizes the mutual information of transceiver ends.The LDC encoder uses a linear combination of a series of base matrices to construct space-time codewords [37,38].Supposing that {S n } N−1 n=0 is a set of the scalar symbols to be sent selected from a specific complex constellation set C, {M n } N−1 n=0 is a M × T base matrix codeword, and the base matrix set is stored at both ends of the transceiver.In it, M is an actual number of transmitting antennae, and T is the duration time of the codeword.Supposing that the mean and variance of S n satisfy: and the formula ε{•} is the expected value.Then, the power limitation condition should be met by the basis matrices: The corresponding transmitter code corresponding to {s n } N−1 n=0 can be constructed by a linear combination of the base matrix and s n weighting coefficients.

S(s
As for the MIMO system of M transmitting antennae M r receiving antennae [39].After matched filtering and symbol level sampling, the received signal matrix in the receiving antenna array is: In it, the received antenna signal matrix in Formula (1) is Y, H s is the MIMO channel coefficient matrix [40,41], and the complex AWGN (additive white Gaussian noise) M r × T matrix is V, in which all elements are not dependent and identically distributed complex Gaussian random variables in a complex domain with a distribution is CN(0, N 0 ).E s is the total power transmitted on the M transmitting antenna.The SNR (signal to noise ratio) of those receiving antennae is written by SNR := E s /N 0 .At this point, the spectrum efficiency of the system is defined as: Symmetry 2019, 11, 1176 4 of 17 In above, |C| indicates the potential of the constellation set C, which is the number of all possible constellations [42,43].

System Model
The block diagram of the transmit antenna selection system based on LDC (LDC-TAS) [44,45] is shown in Figure 1, it is easily found that there is a feedback circle link between the receiver with the transmitter in the system.It is supposed that being M t transmitting antennae and M r receiving antennae in it, but just M transmit antennae are selected by the transmitter to transmit.H being a M r × M t channel coefficient array whose elements are absolutely independent and its distribution is CN(0, 1).At the same time, supposing the receiver has known an ideal CSI (channel status information), and given the channel state H and {M n } N−1 n=0 , then the M r × M t channel array Hs can be selected by one function below: : Symmetry 2019, 11, x FOR PEER REVIEW 4 of 19 In above, |ℂ| indicates the potential of the constellation set ℂ, which is the number of all possible constellations [42,43].

System Model
The block diagram of the transmit antenna selection system based on LDC (LDC-TAS) [44,45] is shown in Figure 1, it is easily found that there is a feedback circle link between the receiver with the transmitter in the system.It is supposed that being Mt transmitting antennae and Mr receiving antennae in it, but just M transmit antennae are selected by the transmitter to transmit.H being a

Maximum Minimum Posterior SNR Criterion for the LDC-TAS
As mentioned above, s H is a more effective channel coefficient array chosen from that function  .

( )
For convenience of analysis, defining a linear transferring array: ≔ ( ), ( ), … , ( ) So the channel coefficient matrix can be expressed as In the equation, (•) is the vectorization operator, and ⊗ is the Kronecker product, so Y can be re-written as below:

Maximum Minimum Posterior SNR Criterion for the LDC-TAS
As mentioned above, H s is a more effective channel coefficient array chosen from that function .
For convenience of analysis, defining a linear transferring array: So the channel coefficient matrix can be expressed as In the equation, vec(•) is the vectorization operator, and ⊗ is the Kronecker product, so Y can be re-written as below: In it, , and H s := H s χ being the M r T × N matrix of the equivalent channel coefficient.
From the expression above, any decoding technology, which is suitable for BLAST can be applied for the LDC system decoding.Due to the complexity of the real practical system, only the linear receiver is taken into account in this article.When the N × M r T matrix equalizer G is multiplied to y, an estimated s can be achieved as below: The original information can be recovered from Ŝ.According to the zero-force (ZF) receiver: However, for MMSE receivers, G can be written as: Maximum minimum posterior SNR criterion: Corresponding in ZF receiver, to N multiplexed sub streams, a posteriori SNR is: For each element in the collection H, conjugate transposition of matrices is () H , we should calculate the minimum SNR in the corresponding N multiplexed sub streams.As for each element of the set Ĥ, a minimum SNR of the corresponding N streams is calculated, so an antenna subset using the most maximum minimum SNR can be chosen.
The bit error rate (BER) and frame error rate (FER) profiles of these two systems are displayed in Figure 2. It is obviously seen that the LDC-TAS performance outperformed STBC-TAS, while the SNR was lower than 14 dB.While SNR was at a high level, because STBC could achieve full diversity gain, its performance was indeed better than LDC-TAS.
From the expression above, any decoding technology, which is suitable for BLAST can be applied for the LDC system decoding.Due to the complexity of the real practical system, only the linear receiver is taken into account in this article.When the r N M T × matrix equalizer G is multiplied to y , an estimated s can be achieved as below: The original information can be recovered from  .According to the zero-force (ZF) receiver: However, for MMSE receivers, G can be written as: Maximum minimum posterior SNR criterion: Corresponding in ZF receiver, to N multiplexed sub streams, a posteriori SNR is: For each element in the collection H , conjugate transposition of matrices is H () , we should calculate the minimum SNR in the corresponding N multiplexed sub streams.As for each element of the set  , a minimum SNR of the corresponding N streams is calculated, so an antenna subset using the most maximum minimum SNR can be chosen.
The bit error rate (BER) and frame error rate (FER) profiles of these two systems are displayed in Figure 2. It is obviously seen that the LDC-TAS performance outperformed STBC-TAS, while the SNR was lower than 14 dB.While SNR was at a high level, because STBC could achieve full diversity gain, its performance was indeed better than LDC-TAS.these two systems was exactly R = 3 bit/s/Hz.The FER and BER profiles of these two systems are displayed in Figure 3. LDC-TAS shows a significant diversity gain contrasted to BLAST-TAS under the equal spectral efficiency.It is worth mentioning that LDC-TAS approximately had a 6 dB performance gain of the system when BER reached 10 −2 , moreover, with the increase of SNR, the performance gain was higher.BLAST-TAS used four transmitters and two receivers, QPSK modulation, two antennae for transmission; LDC-TAS used 16QAM modulation, T = 3.It means that their spectral efficiency of these two systems was exactly R = 3 bit/s/Hz.The FER and BER profiles of these two systems are displayed in Figure 3. LDC-TAS shows a significant diversity gain contrasted to BLAST-TAS under the equal spectral efficiency.It is worth mentioning that LDC-TAS approximately had a 6 dB performance gain of the system when BER reached 10 −2 , moreover, with the increase of SNR, the performance gain was higher.

Antenna Synthesis Via an Innovative HYBRID Genetic Algorithm.
When the antenna subset H , which is the

Antenna Synthesis Via an Innovative HYBRID Genetic Algorithm
When the antenna subset H s , which is the M r × M t channel coefficient matrix [46], was selected via the minimum posterior SNR criterion mentioned in part 4 above, we could think about a linear array displayed in the Figure 4, where the M non-uniform spaced elements were all located alongside an exactly straight line L in the positions x k , k = 0, 1, 2, 3, . . ., M − 1. BLAST-TAS used four transmitters and two receivers, QPSK modulation, two antennae for transmission; LDC-TAS used 16QAM modulation, T = 3.It means that their spectral efficiency of these two systems was exactly R = 3 bit/s/Hz.The FER and BER profiles of these two systems are displayed in Figure 3. LDC-TAS shows a significant diversity gain contrasted to BLAST-TAS under the equal spectral efficiency.It is worth mentioning that LDC-TAS approximately had a 6 dB performance gain of the system when BER reached 10 −2 , moreover, with the increase of SNR, the performance gain was higher.

Antenna Synthesis Via an Innovative HYBRID Genetic Algorithm.
When the antenna subset H , which is the The parameter M can be calculated by the Equations ( 6) and ( 14), which are derived in part 4, and the beam pattern (BP) function of the matrix, p(u), being achieved as below: In Equation ( 15), w k being the weighting coefficients of the kth element, λ means the wavelength of the propagation channel, u = sin θ − sin θ 0 , and θ 0 and θ are respectively the steering angles of these arrays and the incident angle of the impinging plane wave.It is very essential to combine an array configuring to produce a BP that meets some requirements (e.g., side lobes level (SLL)) or rebuilds a special shape (p ref (u)).
First of all, we defined the fitness function f that symbolizes a measurement of the difference among desired and combined BP [47,48], as follows: where ζ = [M; x 0 , . . ., x k , . . ., x M−1 ; w 0 , . . ., w k , . . ., w M−1 ] t being the uncertain quantity arrays and where: In the above equations, µ start is a threshold, which allows removing the main-lobe, which could be derived from the computation of the SLL (side-lobes level) [49][50][51].Moreover, D is an aperture of the array, Q is a normalized constant parameter, S is a region of these values for {p dB (u)/Q} > p re f dB (u), but p re f dB (u) is a desired BP shape.Finally, based on the optimization method, k 1 , k 2 , k 3 , and k 4 were all selected as the normalized parameters.
A fitness function f that has been assumed in Expression ( 2) is extremely non-linear with a lot of the local peak values, where the deterministic processes can be captured.Accordingly, in order to address the maximization issue at hand, we applied a stochastic algorithm to remove these local maxima and effectively explore highly nonlinear search spaces.The GAs already displayed their correctness and pretty effectiveness for optimizing those antenna arrays and seem to be the most possible, feasible, and reasonable option.
GAs are optimization algorithms that are basically built on Darwin's Theory of Evolution.The method assumes an evolution process for a population of these individuals (i.e., a set of experimental ways for the solved problem dealt with) in time that are improving some individual characteristics (i.e., the fitness with several constraints assessed using fitness function f ).
Standard GAs (SGAs) are very different from those other methods with the optimization aims for the following reasons [16,17]: These SGAs cooperate with the parameters' coding over computation, not these parameters themselves; • These SGAs being a multiple-agent hunting process (i.e., a multiple sampling of those hunting spaces); • These SGAs do not have to make use of those derivatives indeed; • These SGAs do not apply the deterministic transition schemes, but use some random ones instead.According to the study [18], these SGAs need to be made for every application to obtain the best results.Consequently, an IGA (innovative genetic algorithm) can be suggested to the classes of the problems for antenna combining.Figure 5 shows the specific flow chart of the IGA presented in our paper, and the main characters of the method including:

•
The application for a hybrid coding way;

•
The non-correlation among the chromosomes' genes (i.e., the power weights parameters and the genes, which representing the position of those array elements should be optimized synchronously);

•
The conception over a prior-knowledge-enhanced operator;

•
The design of the scheme for a mind of the adaptive evolution;

•
The definition of a hybridization design using the local hunting way.As for an antenna combining of the linear (λ/2) mean spaced matrix (to bypass the effects of those grating lobes), to optimize the different sorts of the parameters is the key, including: Power weight of the active elements {wk; k = 0, 1, 2, …, M − 1}, and number of the active elements, M. In this part, a hybrid coding based on the GA is adopted to address this problem properly.The chromosome is supposed as the structure below:  = ;  0 ,  1 , …  −1 ,   , …   ,  −1 ;  0 ,  1 , …  −1 ,   , …   ,  −1 .(16) In Equation (17), bk is a Boolean value (i.e., true or false) giving a status (off or on) of the kth array elements, k being an integer number of (λ/2) intervals between the left array limit (xk = k λ/2) and the kth element, wk is the related excitation coefficient, and the parameter N is the number of intervals (λ/2 in length) in which the array length has been dispersed.
As for the adopted representation, in order to achieve an admissible method and improve the In the SGAs, the individual is coded via a chromosome that is represented by a two-dimensional array, so the optimization problems about the pseudo-Boolean (see from examples [6][7][8]) are properly addressed.A coding process is required once the dispersion factors are concerned.Every coefficient can be symbolized via a series of the Q bits.In it, Q = log 2 (L), and L being the number of these values that the discrete parameters are able to suppose [11].However, due to the time consuming coding/decoding processes [12] and the quantization noise, to use the binary coding is disadvantageous and unpractical once the real unknowns are taken in account.Consequently, the real-valued expression is suggested to overtake those drawbacks above [13].As for an antenna combining of the linear (λ/2) mean spaced matrix (to bypass the effects of those grating lobes), to optimize the different sorts of the parameters is the key, including: Power weight of the active elements {w k ; k = 0, 1, 2, . . ., M − 1}, and number of the active elements, M. In this part, a hybrid coding based on the GA is adopted to address this problem properly.The chromosome is supposed as the structure below: In Equation ( 17), b k is a Boolean value (i.e., true or false) giving a status (off or on) of the kth array elements, k being an integer number of (λ/2) intervals between the left array limit (x k = k λ/2) and the kth element, w k is the related excitation coefficient, and the parameter N is the number of intervals (λ/2 in length) in which the array length has been dispersed.
As for the adopted representation, in order to achieve an admissible method and improve the convergence procedure, the suitable genetic operators should be defined.In the presented application, the mutation keeps a means for investigating the novel regions of addressing space and crossover constitutes in a way to mix the best genes of the present population.In the application some innovative choices have been implemented.

The Selection
The roulette wheel selection [14] using fitness scaling is taken into account.According to the scaling, a rule below is used.
In Expression (18), f is a scaled fitness function, f avg being the average fitness value of the present population, f worst shows the lowest fitness values, and i being a generating index.The parameter variables m and n are specified heuristically to prevent premature convergence and speed up the hunting once this population intends to be convergence [15].

The Crossover
The crossover is aimed at interchanging some information about the gene between those chromosomes.It is doubtlessly problem-oriented, which uses the crossover to increase the specific offspring production.The convergence rate of a maximization procedure is greatly increased by the effective contrive.The different strategy is thought about for the Boolean and real part of the chromosome due to the hybrid chromosome symbolization.Let us consider two selected parents, ξ (a) and ξ (b), and a randomly selected crossover port ks: After a crossover computation, the offspring results are going to be: N−1 ; k s +1 , . . ., w Symmetry 2019, 11, 1176 10 of 17 According to the real parts of these chromosomes, a real-crossover is utilized on the basis of the modified version (for those variable length chromosomes) of the algorithm preliminary proposed in the paper [16], so that: In Equation ( 21), r ∈ [0, 1] is a random number such that the resulting genes belong to the acceptance region defined on the basis of the a priori knowledge: In Equation ( 22) above, w min k and are static parameters constants, and their values are selected to avoid the mutual coupling effects arising in those dense arrays.
At that same moment, the Boolean positions intend to obey a semi-probabilistic distribution.The status of these sensors being respectively changed and preserved with probability r for one child and 1 − r for the others.
The crossover being executed with a probability pc and the reproduction with the probability (1 − p c ).

The Mutation
The mutation is executed using a probability pm on a chromosome of the population.A mutation happens using the probability p bm .To let the variations into the chromosome, these values of the string position are changed.The mutation is executed under the different schemes on the basis of the type of the gene to change.Once those genes that are randomly chosen are two-dimensionally-valued, b k , and the standard binary mutation is performed [17] by using the different probabilities for the birth of the array element or the death.
In addition, a randomly specific mutation already has been constructed for the real-gene as well. [ In Expression (25), η is a random parameter variable value such that the obtained solution is physically reasonable and acceptable.

Elitism
The elitism should be utilized to avoid losing the fittest individual from one generation to another [17].At every generation, the best of chromosome achieved until now is reproduced in the novel population.

Result Analyses
To assess the advantage and effectiveness of the proposed algorithm, several different experiments instances were monitored and explored.In the following part, to assess the flexibility and effectiveness of the presented approach, the specific numerical results are shown and compared with some reference ways (obtainable from the literature [20][21][22]).In the first example, to minimize the maximum side lobes level (Φ slp ) via the changing element weight was solved.As for the problem, in [22] the combining process was utilized to the linear matrix of 25 isotropic elements D = 50 λ in length.The power weight parameters were set to change inside the range [0.2, 2.0] and the parameter u start was fixed to 0.04 to compare the performances of the results achieved via normal GA with the IGA optimization.
As far as this application is concerned, we supposed that the length of this array was dispersed in N = 100(λ/2)-spaced steps via pushing the specific number of the elements M, equals to 25.Such that, the architectures of those series of the chromosome results as below: Figure 6 shows the comparison between the element, the BP, and the position layout of such arrays.The most optimal result is the array with a side lobe peak value located at Φ slp = −14.77dB.However, to the authors' knowledge, the best threshold for the side-lobes level obtained in article [26] was just −14.45 dB approximately, and u ml = 0.191 was the width of the half-beam width.
However, then as far as the main lobe width is concerned, one of the most optimal results was the specific array with a beam pattern (BP), which was characterized by a side lobe peak value being Φ slp = −14.67dB and the main-lobe width value u ml = 0.0204.

Result Analyses
To assess the advantage and effectiveness of the proposed algorithm, several different experiments instances were monitored and explored.In the following part, to assess the flexibility and effectiveness of the presented approach, the specific numerical results are shown and compared with some reference ways (obtainable from the literature [20][21][22]).In the first example, to minimize the maximum side lobes level (Φslp) via the changing element weight was solved.As for the problem, in [22] the combining process was utilized to the linear matrix of 25 isotropic elements D = 50 λ in length.The power weight parameters were set to change inside the range [0.2, 2.0] and the parameter ustart was fixed to 0.04 to compare the performances of the results achieved via normal GA with the IGA optimization.
As far as this application is concerned, we supposed that the length of this array was dispersed in N = 100(λ/2)-spaced steps via pushing the specific number of the elements M, equals to 25.Such that, the architectures of those series of the chromosome results as below: ̅ = { , … ,  , … ,  },  = 1  = 0, … ,  − 1,  = 25. ( Figure 6 shows the comparison between the element, the BP, and the position layout of such arrays.The most optimal result is the array with a side lobe peak value located at Φslp = −14.77dB.However, to the authors' knowledge, the best threshold for the side-lobes level obtained in article [26] was just −14.45 dB approximately, and uml = 0.191 was the width of the half-beam width. However, then as far as the main lobe width is concerned, one of the most optimal results was the specific array with a beam pattern (BP), which was characterized by a side lobe peak value being Φslp = −14.67dB and the main-lobe width value uml = 0.0204.In order to meet the accuracy and completeness, Figure 7 illustrated the specific collection of those best results achieved behind the running of the presented method to highlight an obtainable trade-off among the main-lobe width and the side-lobe peak value.As a reminder, it has to be pointed out that, although the side-lobe decrease of approximately 0.2 dB could be lower than the errors' number in the real phased array, these obtained results are pretty good due to the closeness to the theoretically optimal value.For the convenience to the comparison, this figure also shows the outcome with the procedure based on the SA (Simulated Annealing).
In this testing part, the optimizing process of the quantity of these sensors and of the specific length of the array using a fixed beam pattern shape was thought about.The expected pattern was the same as in [11,16] and illustrated in Table 1 below.The algorithm proposed in the paper intended to work, firstly, as every node writes the local interference power value into the HELLO message, and then broadcasts it to its neighbors.In order to meet the accuracy and completeness, Figure 7 illustrated the specific collection of those best results achieved behind the running of the presented method to highlight an obtainable trade-off among the main-lobe width and the side-lobe peak value.As a reminder, it has to be pointed out that, although the side-lobe decrease of approximately 0.2 dB could be lower than the errors' number in the real phased array, these obtained results are pretty good due to the closeness to the theoretically optimal value.For the convenience to the comparison, this figure also shows the outcome with the procedure based on the SA (Simulated Annealing).At last, the isophoric matrices [19] were also thought about as well.In this instance, the power weights of these array elements were statically set (wk; k = 0, 1, 2, 3, ..., N − 2, N − 1) and an optimization procedure took these array element positions into account.Figure 8 draws the beam pattern and the specific location position layout in contact with the best optimal ways, according to these peak values of the minimum side-lobe level (Φslp = −13.06dB as well as uml = 0.0170; IGA (a)) and according to the optimal trade-off among the side-lobe level value and those main-lobe widths (Φslp = −12.32dB as well as uml = 0.0126; IGA (b)), achieved via the ways of those IGA-based methods.These specific characters of the array synthesized using the SA-based procedure are given (Φslp = −12.07dB as well as uml = 0.0133; SA).In this testing part, the optimizing process of the quantity of these sensors and of the specific length of the array using a fixed beam pattern shape was thought about.The expected pattern was the same as in [11,16] and illustrated in Table 1 below.
The algorithm proposed in the paper intended to work, firstly, as every node writes the local interference power value into the HELLO message, and then broadcasts it to its neighbors.At last, the isophoric matrices [19] were also thought about as well.In this instance, the power weights of these array elements were statically set (w k ; k = 0, 1, 2, 3, . . ., N − 2, N − 1) and an optimization procedure took these array element positions into account.Figure 8 draws the beam pattern and the specific location position layout in contact with the best optimal ways, according to these peak values of the minimum side-lobe level (Φ slp = −13.06dB as well as u ml = 0.0170; IGA (a)) and according to the optimal trade-off among the side-lobe level value and those main-lobe widths (Φ slp = −12.32dB as well as u ml = 0.0126; IGA (b)), achieved via the ways of those IGA-based methods.These specific characters of the array synthesized using the SA-based procedure are given (Φ slp = −12.07dB as well as u ml = 0.0133; SA).Due to the random character of the algorithm proposed in this article, the statistical factors that are relevant to the simulations collection were assessed.Table 2 offers all the values over the best, the worst, the average resulting values, and the standard deviation values according to the number of these elements, the minimum peak value of the side-lobe and these values of the main-lobe width.It is beneficial to see that the filling percentage is lower than the value obtained from the paper [7,8]   Due to the random character of the algorithm proposed in this article, the statistical factors that are relevant to the simulations collection were assessed.Table 2 offers all the values over the best, the worst, the average resulting values, and the standard deviation values according to the number of these elements, the minimum peak value of the side-lobe and these values of the main-lobe width.It is beneficial to see that the filling percentage is lower than the value obtained from the paper [7,8] and the side-lobe peak value at all times.In summary, it could be highlighted that some specific taking-away of the symmetry restrictions significantly improves its performances, indeed (these symmetrical arrays were concerned inside the literatures [7,8]).

Complexity Analyses
To evaluate this algorithm, a specific complexity for this method was given in detail in the following part, which is concerned only with the multiplying times of the complexity complex, and the classical matrix operational way.As far as a same set of the parameters of a mobile channel is concerned, it is enough to only analyze one combination of the channel parameters such as H in detail.
In the criterion of the max-min-post-SNR, without considering the operation time of H s , H s = H s χ is a matrix multiplication of order N r T × NT and NT × Q, the multiplying times is N r NQT 2 , the multiplying times in H H j H j −1 is Q 2 N r T, and the inversion times of the Q × Q matrix by the Gauss-Jordan method are about 0.5Q 3 + 0.5Q 2 , so the total complexity is N r NQT 2 + TQ 2 N r + 0.5Q 3 + 0.5Q 2 .

Conclusions
The coming 5G systems need a reliable and safe mobile communication link with a specific large capacity using much easier and cheaper methods.Since the LDCS meet the information-theoretic optimality property and the selection for the transmitting antennae offers the additional order for the diversity, so the LDC-TAS system is introduced in the paper firstly.Contrasted to a multiple antenna-fading channel, the LDC-TAS system is a simple and cheap way to achieve the diversity advantage, and to maintain the high rate for transmission.Authors educed the specific criteria for the antenna selection when the ZF receiver was utilized as well as investigated it by the computer numerical simulations.Finally, the computing simulation results show that this specific max-min-post-SNR criterion, which is specially produced for the ZF receiver, outperforms some other criteria [19,20].Moreover, LDC-TAS has the more significant diversity advantage than BLAST-TAS under the same spectral efficiency and is better than STBC-TAS inside the low SNR environment as well.

Figure 1 .
Figure 1.Transmit antenna selection for linear dispersion codes.

Figure 2 .
Figure 2. Performance of linear dispersion code transmitting antenna selection (LDC-TAS) vs. space-time block code (STBC)-TAS.BLAST-TAS used four transmitters and two receivers, QPSK modulation, two antennae for transmission; LDC-TAS used 16QAM modulation, T = 3.It means that their spectral efficiency of

Figure 4 .
Figure 4.The geometry of the linear array.

Figure 5 .
Figure 5. Flowchart of the innovative genetic algorithm.

Figure 6 .
Figure 6.Optimization of the element power weights-(a) beam pattern (BP) with a side-lobe peak value of −14.77 dB; and (b) positions and power weights of those array elements.

Figure 7 .
Figure 7. Optimization of the element power weights-peak side-lobe value level as a specific function of the main-lobe width ((dB) units; simulated innovative genetic algorithm (IGA); O: Simulated SA).

Figure 7 .
Figure 7. Optimization of the element power weights-peak side-lobe value level as a specific function of the main-lobe width ((dB) units; simulated innovative genetic algorithm (IGA); O: Simulated SA).

Figure 8 .
Figure 8. Optimization of the element location positions and power weights-optimal combinations for isophoric arrays.(a) Beam patterns; and (b) array placement.

Figure 8 .
Figure 8. Optimization of the element location positions and power weights-optimal combinations for isophoric arrays.(a) Beam patterns; and (b) array placement.

Table 2 .
Statistical behavior the quantity of those active elements, side lobe peak value as well as the width of main-lobe behind several dozens of the procedure implement.
Note: In this table, the word "Std.Dev" means "the standard deviation values".