The Effect of Mutual Coupling on a High Altitude Platform Diversity SystemUsing Compact Antenna Arrays

We analyze the destructive effects of mutual coupling and spatial correlation between the separate antenna elements on a combined diversity system consisting of multiple HAPs (High-Altitude Platforms) employing various compact MIMO (Multiple-Input Multiple-Output) antenna array configurations, in order to enhance the mutual information in HAP communication links. In addition, we assess the influence of the separation angle between HAPs on system performance, and determine the optimal separation angles that maximize the total mutual information of the system for various compact MIMO antennas. Simulation results show that although the mutual information is degraded by mutual coupling and spatial correlation, the proposed HAP diversity system still provides better performance compared to a nondiversity system for all tested scenarios.


Introduction
High-Altitude Platforms (HAPs) are quasistationary aerial platforms operating in the stratosphere.This emerging technology is preserving many of the advantages of both satellite and terrestrial systems [1][2][3][4][5][6] and presently started to attract considerable attention in Europe through the European Union COST 297 Action [1].Using narrow bandwidth repeaters on HAPs for high speed data traffic has several advantages compared to using satellites [1], especially when operating in a local geographical area.In this paper, we propose a multiple HAP diversity system in order to enhance the mutual information in HAP communication links.To further enhance the system performance, independent spatial or polarization channels can be acquired by means of multiple antenna elements at both the transmitter and the receiver, and thus this technique is referred to as Multiple-Input Multiple-Output (MIMO) communications.The MIMO wireless communication is also combined with OFDM (Orthogonal Frequency Division Multiplexing) to form a combined HAP MIMO-OFDM diversity system.
The proposed system consists of virtual MIMO spatial channels (created by multiple HAP diversity) in conjunction with the polarization and antenna pattern diversity (formed via special compact MIMO antenna arrays), as shown in Figure 1.Two compact MIMO antenna array configurations (the MIMO-Cube and MIMO-Tetrahedron) will be investigated, and their performance, in term has of mutual information, is analysed.The enhanced mutual information that is obtained through use of these compact MIMO antenna arrays, and of multiple HAPs, will degrade if the signals of the virtual MIMO spatial channels are correlated.Correlation of the signals can be caused by either low multipath richness in the radio propagation environment or by mutual coupling between the compact antenna array elements.Therefore, the effect of mutual coupling and spatial correlation will be taken into account when performing the simulations for this combined HAP diversity system.We will also analyse the influence of the separation angle between the multiple HAPs on system performance, and determine the optimal separation angles that maximize the total mutual information of the system for the various compact MIMO antennas.
The paper is organized as follows: in Section 2 the proposed combined MIMO-HAP diversity system is introduced.The mathematical analysis of the spatial and polarization channels is presented in Section 3. In Section 4, the various compact MIMO antenna array configurations, their characteristics, and modelling are introduced.The analysis of mutual coupling and spatial correlation effects on these channels is explained in Section 5. Simulation results and comparisons are presented in Section 6. Finally the conclusions are presented in Section 7.

The Proposed Multiple HAP Diversity System
In this paper, we propose an application for high data rate transmissions using a system employing multiple HAPs.This system consists of a number of uncorrelated MIMO channels created in the space/polarization domain using HAP diversity in combination with a special type of compact MIMO antenna arrays [7][8][9] as discussed in Section 4.
Figure 1 shows the location diversity setup for the case of three HAPs separated by the angle θ.
The propagation channel H ch from the transmitter to the receiver can be seen as a simple transformation that contains the distance dependent-free space gain of the signal being transmitted where |r m − r n | is the distance along the path between transmitting antenna element n and receiving antenna element m.There are no atmospheric interferences and the noise in the system is modelled as uncorrelated Gaussian noise.The total MIMO channel can then be assembled as where M rx and M tx are transformations of the transmitter and receiver antenna signals into and out of the mode domain as described in Section 3.

The Space-Polarization Channel
The polarization and radiation pattern of the electromagnetic field can be expressed as a multipole expansion [10] of the spherical wave field emanating from a virtual sphere encompassing the antenna array that is being analyzed.This expansion consists of weighted orthogonal vector base functions on the surface of the virtual sphere that allow for a solution to Maxwell's equations in spherical coordinates which can be written as where η is the intrinsic impedance of vacuum and underscoring denotes a field vector.The base functions X lm (ϕ, θ) represent mode (l, m) of the electromagnetic field in spherical coordinates when the far-field of the antenna is projected onto the virtual sphere.The functions g l and f l in (3) are spherical Hankel functions of order l representing an outgoing (transmitted) or incoming (received) wave.The scalar weights a E and a M are the corresponding coefficients which will give the gain of each orthogonal function (mode) for a particular electromagnetic far-field pattern, as shown by Using (4) we can calculate which modes are active on any arbitrary antenna enveloped by the virtual sphere by knowing the electric and magnetic field patterns of the antenna.These modes are theoretically orthogonal to each other and therefore represent independent ports of the antenna.The transmitting antenna-channel transform matrix M tx is then assumed to be the linear transformation of the input signals x into the mode domain a tx according to a tx = M tx x, and for the receiving channel-antenna, we have a similar transformation matrix from the mode domain a rx into the output signals y of the system following y = M rx a rx , where a tx and a rx are vectors containing the mode gains for a specific antenna type.

The Compact MIMO Antenna Array Configurations
Each transmit and receive antenna arrangement of the proposed system consists of a special compact MIMO antenna array configuration, which possesses a structure of different designs and complexities as shown in Figure 2.
The first compact MIMO antenna array is known as the MIMO-Cube [10], which consists of twelve electric dipoles positioned on the twelve edges of a cube.It can be seen from Figure 3 that this antenna array theoretically has twelve independent antenna ports.
The second compact antenna that we have investigated is known as the MIMO-Tetrahedron, which has six electric dipoles placed at the six edges of a tetrahedron (see Figure 2), and so theoretically has six independent ports as can be seen in Figure 4.
International Journal of Antennas and Propagation   In Figure 5 we show an example of the activation of the first 30 modes calculated by (4) for the MIMO-Tetrahedron antenna array.
In order to separate the modes in Figure 5 into independent MIMO channels, we employ a singular value decomposition (SVD) operation where U and V are unitary eigenvector matrices of the MIMO channel modes and Σ comprises the subchannel gains of the independent ports of the antenna array.The orthogonalization of the channel into the independent subchannels is accomplished by multiplying the signal to be transmitted x with the matrix V on the transmitter side of the channel, and multiply the received signal y on the receiver side of the channel with the matrix U H according to where y is a vector containing the received and decoded signal from the independent MIMO subchannels.Inserting (5) into (6) then orthogonalizes the channels into where w = U H v is the received noise.system mutual information as will be shown in the results section.
Through the combination MIMO-OFDM (Orthogonal Frequency Division Multiplexing), we also achieve orthogonality over the frequency domain by using the Discrete Fourier Transform (DFT) with a Cyclic Prefix (CP).Assuming a cyclic prefix CP > τ d • f s , where f s = 1/T s and T s is the symbol time, we can write the received signals as a circular convolution over the DFT-frame where K is the size of the DFT, m = {1, . . ., M} is receiving antennas, n = {1, . . ., N} is transmitting antennas, and v m (t) is AWGN.If we have a MIMO antenna system with N transmitting antennas and M receiving antennas, we can then write the signal for subchannel k in the frequency domain between any pair of transmitting and receiving antennas as where k denotes each separate subcarrier.If we write (9) in vector notation, we get where R(k) is an M × 1 vector and S(k) is an N × 1 vector of the received and transmitted signals, respectively.H(k) is the M × N frequency response matrix of the channel between N transmitters and M receivers.The noise in the system W(k) is an M ×1 vector assumed to be additive white Gaussian noise.The correlation matrix of the noise vector W(k) is thereby where σ 2 n is the variance of the noise and I M is the M × M identity matrix.Since we are using singular value decomposition (SVD) technique, we can now write the channel frequency response matrix H(k) as where Σ(k) is an M × N matrix containing singular values that are larger than zero where r is the rank of the matrix H(k), and the M ×M matrix U(k) and the N × N matrix V(k) are unitary and contain the corresponding eigenvectors as matrix column vectors.To obtain a diagonalized system, we then define where Since the MIMO-OFDM channels in (12) are uncorrelated and the correlation of the noise I M , then we can write the mutual information of the system [11] as where σ 2 xm (k) is the variance of the separate uncorrelated input signals in vector X(k) and r is the rank of the system.The mutual information in ( 14) is constrained by the total radiated power from the transmitting antennas, defined as To maximize the total sum of mutual information in all the subchannels, we use the so-called "water-filling" technique in which we allocate more power to the subchannels  with high eigenvalues.The optimal "water-filling" solution [12] is then given by where k = {0, . . ., K − 1}, m = {1, . . ., r}, σ 2 n is the variance of the noise, and γ is a predefined threshold level of the signal-to-noise ratio in the system and it is dependent on the total transmitted power P rad in (15).The average signal-tonoise ratio SNR avg in the simulations is calculated as where N active is the number of active channels used for which the variance of the input signals σ 2 xm (k) > 0.

Mutual Coupling and Spatial Correlation
The interaction between the elements of the compact antenna array will cause mutual interference in the radiation patterns of the array elements, and this effect is known as "mutual coupling" [13][14][15][16][17].The mutual coupling between the elements of the compact antenna array can be viewed as a multiport network on the transmitter and receiver sides, respectively, in which the relationships between the terminal voltages and currents of the different ports can be define as [13], where Z tx and Z rx are the impedance matrices of the antenna array.Z tx and Z rx have been calculated numerically by using the finite element method (FEM) for the antenna arrays shown in Figure 2.An infinite approximation of the array is assumed by not taking the edge effects of the finite array into account.The terminal voltages are then associated with the source voltages through the source impedances Z S , as described in [13],

S
, where u (TX)

S
is the source voltage at the transmitter and u (RX) S is the source voltage at the receiver generated by the impinging electromagnetic wave.The impedance coupling matrices are denoted as Z TX and Z RX , and are defined by where Z tx and Z rx are the impedance matrices of the transmitting and receiving compact antenna arrays, Z S and Z L are diagonal matrices containing the source and load impedances which have been chosen as the complex conjugate of the self-impedance of the antenna given by the diagonal [Z tx ] ii and [Z rx ] ii of the impedance matrices Z tx and Z rx , and ) are normalizing factors [13].In (8) it was shown that the complex envelope of the received signal at the receiving compact array after matched filtering can be expressed in vector notation as where s is the N × 1 vector containing the signals from the N transmitting antenna elements, v is the M × 1 zero-mean AWGN vector, and H is the normalized M × N channel matrix that, assuming separable transmit and receive spatial correlation matrices, is modeled using the Kronecker model [18,19], where H 0 is an i.i.d.channel response without spatial correlation and mutual coupling.R rx and R tx are the separable spatial correlation matrices on the receiving and transmitting sides, respectively.
International Journal of Antennas and Propagation  The spatial correlation matrices R tx and R rx are calculated between the antenna element positions and polarization states according to [20] ρ p,q = ϕ θ ζ κ α q (Θ)α * p (Θ) sin(θ)p(Θ)dκ dζ dθ dϕ ϕ θ ζ κ α q (Θ) where element ρ p,q of matrix R tx or R rx is the spatial correlation between antenna elements p and q. ϕ and θ are the spherical coordinates expressing the spatial domain and ζ and κ are the polarization angle and phase difference, respectively, that account for the effects in the polarization domain [20].In (23), α(Θ) is the steering vector and p(Θ) is the joint probability distribution function of the parameter vector Θ = [θ ϕ ζ κ] T , where T denotes a transpose operator.It is assumed that all the parameters are independent and that p(ζ) = u{0, π} and p(κ) = u{−π, π} are uniformly distributed.

Simulation Results
In this section, we investigate the mutual information improvement resulting from the use of multiple HAPs system employing compact MIMO antenna arrays and compare the results to that obtained from a system without diversity.In addition, we will also show the effects of mutual coupling and spatial-polarization correlation (calculated by the Finite Element Method) on the mutual information of the system.The mutual information results shown in Figure 7 are obtained for a system of HAPs operating at an altitude of 22 km and with a separation angle of 20 degrees, and plotted against the average signal-to-noise ratio of the system.The elevation angle is held at 10 degrees throughout the simulations.
It is evident from Figure 7 that the HAP diversity system provides superior performance as compared to the single HAP case, and the MIMO-Cube antenna provides more mutual information than the MIMO-Tetrahedron antenna which is the result of the higher number of acquired HAP platforms.Simulation results also show that although the mutual information is degraded by correlation and mutual coupling, we still achieve significant gain compared to the single HAP case.From Figure 7 it is also clear that for example, for an average SNR of 20 dB we get a 110% more mutual information from the ideal MIMO-Cube antenna than from the ideal MIMO-Tetrahedron antenna, and a 228% more mutual information compared with the nondiversity single HAP system.If we take into account the correlation and the mutual coupling, the MIMO-Cube antenna still has a 108% more mutual information than the MIMO-Tetrahedron antenna and a 195% more mutual information than the ideal nondiversity single HAP system.
In Figure 8, we show the effect of separation angles between HAPs on the mutual information of the two compact MIMO antenna arrays.The separation angle θ (see Figure 2) is varied from 0 to 35 degrees in steps of 5 degrees (due to computational limitations) and plotted for an average SNR of 20 dB.It is clear from this figure that the separation angle between HAPs has a great impact on the mutual information of the system.It is also evident that the optimal separation angle that maximizes the total mutual information of the system at this SNR is found to be approximately 10 degrees for the MIMO-Tetrahedron antenna and approximately 15 degrees for the MIMO-Cube antenna.The degradation of mutual information performance at wider separation angles is due to the increase in distance between the HAPs and the ground stations.Figure 8 also suggests that the vector element antenna array is superior to the MIMO-Tetrahedron array at wide separation angles even though they access the same number of independent subchannels.

Conclusions
In this paper, we have investigated the degradation in mutual information, due to spatial correlation and mutual coupling, when using different compact MIMO antenna arrays in a multiple HAP system employing multiple spatialpolarization modes.Simulation results show that the multiple HAP diversity systems utilizing compact MIMO antenna arrays provide significant mutual information enhancement, even though it can also be noticed that small degradation in mutual information is resulted due to the effects of spatial correlation and mutual coupling between the separate antenna elements of the MIMO-Cube and MIMO-Tetrahedron arrays.It was also shown that the performance of the compact MIMO-Cube antenna is superior to the MIMO-Tetrahedron antenna since it has twice the number of independent channels which will result in more system mutual information.Finally, we have also shown the effects of the separation angle between HAPs on system performance, and determined that the optimal separation angle that maximizes the total mutual information of the system is within a narrow sector of approximately 20 degrees for the compact MIMO antenna arrays investigated in this paper.

Figure 1 :
Figure1: The MIMO-HAP diversity system with three HAPs and the channel paths from the transmitter (Tx) antenna array to the receiver (Rx) antenna array.

Figure 2 :Figure 3 :
Figure 2: The structure of the two compact MIMO antenna array configurations: (a) MIMO-cube and (b) MIMO-tetrahedron.

Figure 4 :
Figure 4: The theoretical antenna pattern when all six channels of the MIMO-Tetrahedron antenna are activated.The three antenna beams at elevation angle 30 degrees have vertical linear polarization and the three beams at elevation −45 degrees have horizontal linear polarization.

Figure 5 :
Figure 5: The activated modes of the MIMO-Tetrahedron antenna array.Each bar-graph represents the modes activated on each of the six antenna array elements.
noise ratio SNR (dB) Ideal MIMO-cube MIMO-cube with mutual coupling and 6 HAPs Single HAP system Vect.elem.ant.or ideal MIMO-tetrahedron. 3 HAPs MIMO-tetrahedron with mutual coupling.3 HAPs Mutual coupling and spatial correlation.3 HAPs MIMO-cube with mutual coupling and spatial correlation, 6 HAPs

Figure 7 :
Figure 7: The mutual information versus average SNR for different configurations of the MIMO-Cube and MIMO-Tetrahedron antenna arrays in comparison to the ideal vector element antenna array and a single HAP system.The elevation angle is 10 degrees and the separation angle between HAPs is 20 degrees.

Figure 8 :
Figure 8: The total system mutual information versus separation angle between HAPs for an average SNR of 20 dB for the two compact MIMO antenna array configurations.