MIMO Antenna Array Design with Polynomial Factorization

One of the main advantages of multiple-input multiple-output (MIMO) antenna is that the degrees-of-freedom can be significantly increased by the concept of virtual antenna array, and thus theMIMO antenna array should be carefully designed to fully utilize the virtual antenna array. In this paper, we design the MIMO antenna array with the polynomial factorization method. For a desired virtual antenna array, the polynomial factorization method can optimally design the specifiedMIMO transmitter and receiver.The array performance is examined by analyzing the degrees-of-freedom and statistical output signal-to-interference-plus-noise ratio (SINR) performance. Design examples and simulation results are provided.


Introduction
Multiple-input multiple-output (MIMO) antenna array has received much attention in communication, radar, and navigation systems [1][2][3][4].MIMO idea is not new; its origin in control systems can be traced back to 1970s [5].The early 1990s saw an emergence of MIMO ideas into the field of wireless communication systems [6][7][8].More recently, the ideas of MIMO appear in sensor and radar systems [9][10][11][12].Different from conventional phased-array systems [13], in MIMO antenna systems each antenna transmits a unique waveform, orthogonal to the waveforms transmitted by other antennas, and a matched filter bank is used to extract the orthogonal waveform components.When the orthogonal signals are transmitted from different antennas, the return of each orthogonal signal will carry independent information.The phase difference caused by different transmitting antennas along with the phase differences caused by different receiving antennas can form a new virtual antenna array steering vector.With optimally designed antenna array geometry, we can create a very long array steering vector with a small number of antennas [14,15].More importantly, this provides good flexibility and reconfigurability in array configuration, thus enabling a flexible and reconfigurable system.
Although MIMO antenna arrays are widely used in various applications, radar-related application is assumed in this paper.Even so, the rich literature of research in this area appears in threefold: (a) the comparison between conventional phased-array and MIMO array, and the development of MIMO architecture [16][17][18]; (b) MIMO waveform diversity and transmitter design [19,20]; and (c) the system performance and the receiver design for enhanced performance [21,22].In fact, one of the main advantages of MIMO antenna array is that the degrees-of-freedom can be significantly increased by the concept of virtual antenna array, and thus the MIMO antenna array should be carefully designed to fully utilize the virtual antenna array.This paper focuses on the MIMO sparse antenna array design to obtain a uniform virtual array for the system.
Many array design methods have been proposed for conventional phased arrays [23][24][25][26].Specially, linear programming-based algorithms were proposed by [27,28] for reduced sidelobe thinned array design.Genetic algorithm was proposed by [29] to optimally design thin large phased arrays.Simulated annealing-based technique was proposed by [30] to synthesize the position and weight coefficients of a linear array that minimizes peak sidelobes.A hybrid approach combining particle swarm optimization with combinational techniques was proposed by [31] to synthesize planar thinned arrays.Recently a general polynomial factorization-based design of sparse periodic linear arrays was proposed by [32] to design sparse periodic arrays.In the method, transmit and receive aperture polynomials are selected such that their product results in a polynomial representing the desired combined transmit/receive effective aperture function.
In this paper, we use the polynomial factorization method proposed in [32] to design MIMO antenna array and comparatively analyze the designed array performance.To obtain the desired virtual MIMO antenna array, we use the polynomial factorization-based method to design the transmit and receive antenna arrays.Since there may be multiple solutions for the transmit and receive antenna arrays and different solutions may have different array performance, we use adaptive beamforming algorithm to comparatively analyze the statistical output signal-to-interference-plus-noise ratio (SINR) performance of the designed arrays.Extensive design examples and simulation results are also provided.
The remaining sections are organized as follows: after introducing the polynomial factorization in Section 2, we describe the polynomial factorization-based MIMO antenna array design in Section 3. Next, the output SINR performance of the designed MIMO antenna arrays is analyzed by adaptive beamforming on the receiver in Section 4, and design examples and simulation results are provided in Section 5. Finally, this paper is concluded in Section 6 with a short discussion of future work.

General Polynomial Factorization Method
The polynomial factorization method is based on factorizing a polynomial   () of degree  in positive powers of  and with unity coefficients [33]   () =  ∑ =0   . ( Let the positive integer of the number of coefficients  + 1 be expressible as a product of  + 1 irreducible positive integers   + 1, 0 ≤  ≤ : Considering the product   0 ()  1 (  0 ), we have Therefore,   () can be expressed as [32,33] with where The factorization of a specified  + 1 into a product of irreducible integers can be carried out using Euclid's algorithm [34].When  + 1 can be expressed as a power of two, that is,  + 1 = 3 ] ,   () can be factorized as ) .
Note that the solutions to the cases when  + 1 cannot be expressed as a power of two can be found in [32].The decomposition in ( 6) has been used in designing a sparse antenna array with a uniform aperture function [35] and a linearly tapered aperture function [36].
As an example, we suppose  + 1 = 10; using Euclid's algorithm we then have Then, according to the above polynomial factorization method, we can get or That is to say, there may be multiple possible factorizations for a given polynomial.

Polynomial Factorization-Based MIMO Array Design
Consider a MIMO radar system with a transmitter equipped with  colocated antennas and a receiver equipped with  colocated antennas.Suppose that both the transmit and receive arrays are close to each other in space (possibly the same array) so that they see targets at the same directions.The signal received by each antenna is the weighted summation of all the transmitted waveform where   () is the received signal at the th antenna,   () is the transmitted waveform at the th antenna, and  , is the channel coefficient with the th antenna as input and the th antenna as output.Suppose x() = [ 1 (),  2 (), . . .,   ()] and () = [ 1 (),  2 (), . . .,   ()], and the above equation can be rewritten in a matrix where  is the target direction,   is the complex-valued reflection coefficient of the focal point   , (⋅)  is the transpose operator, and a() and b() are the actual transmit and receive steering vectors associated with the direction .
When the transmitted waveforms are orthogonal, where () * denotes a conjugate operator.The radar return due to the th transmitted waveform can be separated by matched-filtering the received signal with   (),  = 1, 2, . . ., ; that is, where  0 is the radar pulse width.Therefore, there are a total of × independent signals, which can be written as a × 1 virtual data vector where ⊗ is the Kronecker product; that is, the virtual antenna array is Suppose that the th transmitting antenna is located at  , = (/2)  and the th receiving antenna is located at  , = (/2)V  , where  is the wavelength.Considering a farfield point target, the transmitter and receiver steering vectors can be represented, respectively, by Note that the first antenna of a() and b(), respectively, is the reference antenna.According to (15), we can get The array response in the th receiving antenna for the th transmitted waveform can be expressed as It can be noticed that the MIMO antenna array response is the same as the target response received by a receiving array with  antennas located at { , +  , } ,  = 1, 2, . . ., ,  = 1, 2, . . ., . (19) This -element array is just the virtual MIMO array.An utmost of -element virtual array can be obtained by using only  +  physical antenna elements.It is as if we have a receiving array of  elements.
According to (18), the far-field radiation pattern at an angle  away from the broadside of the virtual MIMO array with uniform weighting is expressed as where   () = ∑  =1     sin  and   () = ∑  =1  V   sin  denote, respectively, the far-field radiation pattern for the transmitter and receiver.For a more general format, by substituting  =  2 sin / with  being the interelement spacing, the virtual array, transmit array, and receive array radiation patterns can be expressed, respectively, as where   (),   (), and   () denote the element weighting functions for the virtual array, transmit array, and receive array, respectively.It follows also that [32]   () =   ()   () .
Therefore, the far-field radiation pattern of the virtual array is equal to the product of the transmitter radiation pattern and receiver radiation pattern.This means that we can use different transmit and receive arrays with a combined virtual MIMO array that is a close equivalent to that of a single full array with no missing elements.To reach this aim, like the general polynomial factorization-based design of sparse periodic linear array discussed in [32], we can design the transmit and receive arrays based on factorizing the desired MIMO virtual array polynomial   () as a product of two lower order polynomials with some zero-valued coefficients.
In this paper, we aim to optimally design the transmit and receive arrays for a given virtual MIMO antenna array.The design steps are illustrated in Figure 1.

Array Performance Analysis by Adaptive Beamforming
The output of the matched filter of the MIMO array system can be written as where n denotes the noise plus interference.The output SINR of the beamformer directly impacts the MIMO array performance.The minimum variance distortionless response (MVDR) [37] is a widely employed beamformer choice.The MDVR beamformer for a look angle  is given by where (⋅)  is the conjugate transpose and R is the covariance matrix of the noise and interference signal.The beamformer output can then be expressed as Correspondingly, the output SINR can be expressed as  Since the output SINR is random due to the random covariance matrix R, we can examine the SINR distribution under the standard assumption that n has a complex Gaussian distribution with zero mean and covariance matrix R. In this case, the probability density function (PDF) of SINR can be represented by [38,39] where  is the number of the training signals for obtaining the covariance matrix,  dof is the number of the system degrees-of-freedom, and  max is the maximum output SINR which is achieved when the covariance matrix is known exactly Figure 2 shows the PDF   () as a function of / max .The detailed derivation of ( 27) can be found in the appendix of [38].The corresponding cumulative distribution function is where Γ(; , ) is the incomplete Beta function Target power (dB) The expectation and variance of the output SINR can be derived, respectively, as It can be noticed from (31) that the expectation of the output SINR depends mainly on the ratio  dof /, not the number of the training data, .Equations ( 31) and ( 32) can be used to evaluate the designed MIMO array performance.Figure 3 gives the comparative standard deviation (STD) of the output SINR performance, where  = 20 is assumed.It can be noticed that, for a given number of the training signals, more degrees-offreedom mean a better STD of the SINR performance. Define The corresponding  is To ensure that the mean output SINR of the adaptive beamformer is within 3 dB from its maximum achievable value, the training signals should be satisfactory with [38]  ≥ 2 dof − 3.
If we want to achieve  = 0.95, the  should be

Design Examples and Performance Analysis
We consider two design examples with given virtual MIMO antenna arrays.
According to (2), we factorize  + 1 = 10 into two prime integers 2 and 5, that is, 10 = 2 × 5.Then, according to (4), the virtual array radiation pattern   () can then be factorized as We can use the first factor for the transmit array and the second factor for the receive array or the first factor for the receive array and the second factor for the transmit array; that is, or Suppose that the number of the training signals is  = 25; Figure 4 shows the possible two transmit-receive arrays.It can be noticed that 7 antennas are required for the MIMO antenna arrays and the number of the system degrees-offreedom is 10. Figure 5 gives the statistical mean and STD of the output SINR performance.Note that, although the two arrays have the same statistical SINR performance (because they have the same number of degrees-of-freedom), they should be seen as two different MIMO antenna arrays in the design because they have different transmit array gain and receive array gain (MIMO antenna array has no transmit array gain whereas phased array has both transmit array gain, and receive array gain).According to (2), we can factorize  + 1 = 12 into the following three ways.Factorization 1.We can factorize 12 into three prime integers 2, 2, and 3, that is, 12 = 2 × 2 × 3.Then, according to (4), the virtual array radiation pattern   () can then be factorized as Factorization 2. We can factorize 12 into three prime integers 2, 3, and 2, that is, 12 = 2 × 3 × 2.Then, according to (4), the virtual array radiation pattern   () can then be factorized as Factorization 3. We can factorize 12 into three prime integers 3, 2, and 2, that is, 12 = 3 × 2 × 2.Then, according to (4), the virtual array radiation pattern   () can then be factorized as Then there are nine possible transmit-receive array designs: Suppose that the number of the training signals is  = 64; Figure 6 shows the possible nine transmitter-receiver arrays.It can be noticed that 8 antennas are required for the transmitreceive pairs (a), (b), (d), (f), (h), and (i), and 7 antennas are required for the transmit-receive pairs (c), (e), and (g).All of the nine transmitter-receiver arrays have 10 degrees-offreedom. Figure 7 gives the statistical mean and STD of the output SINR performance.Since they have the same number of degrees-of-freedom, the nine transmit-receive arrays have the same statistical SINR performance.

Conclusion
In MIMO antenna array systems, the transmitter and receiver should be carefully designed to fully utilize the virtual MIMO antenna array.In this paper, we designed the MIMO antenna array with the polynomial factorization method and comparatively analyzed the designed array performance.For a desired virtual antenna array, the polynomial factorization method can optimally design the specified transmitter and receiver antenna arrays.Design examples are also provided, which are verified by simulation results.Note that two MIMO antenna arrays may have the same statistical SINR performance; they often should be seen as two different MIMO antenna arrays in the design because they have different transmit array gain and receive array gain.In this paper, uniform virtual MIMO antenna array elements are assumed.Target power (dB) We are aware that nonuniform virtual MIMO antenna array elements may be desired in actual MIMO antenna arrays.This topic will be further investigated in our subsequent work.

Figure 1 :
Figure 1: Illustration of the MIMO antenna array design steps.

Figure 2 :
Figure 2: The PDF   () as a function of / max .

Figure 3 :
Figure 3: Comparative STD of the output SINR performance.

Figure 4 :Figure 5 :
Figure 4: Illustration of the designed MIMO antenna array positions.The -axis stands for the MIMO antenna position. = 1 denotes an antenna and  = 0 means no antenna.

Figure 6 :
Figure 6: Illustration of the nine possible transmit-receive array positions for the MIMO antenna array.The  axis stands for the MIMO antenna position. = 1 denotes an antenna and  = 0 means no antenna.
Mean of the output SINR STD of the output SINR Statistical output SINR performance (dB)