On the Reduction of Transmission Complexity in MIMO-WCDMA Frequency-Selective Fading Orientations via Eigenvalue Analysis

In this paper, a novel transmission strategy for Mutliple Input Multiple Output Wideband Code Division Multiple Access (MIMO-WCDMA) orientations operating in frequency-selective fading environments is investigated, in terms of overall algorithmic complexity reduction. To this end, Principal Component Analysis (PCA) is employed on the received data matrix, in order to define the significant terms that are taken into account during transmission matrix formulation. According to the presented results, feedback information of only the primary eigenvector of the corresponding covariance matrix of the received data matrix is required, in order to maintain the mean Bit Error Rate (BER) at acceptable levels. In particular, a complexity reduction of up to 10% can be achieved, when comparing BER values derived by the selection of all components of the received covariance matrix during transmission matrix formulation, and the corresponding BER when selecting half of the components. This reduction is maintained to 10%, when considering a realistic four-element antenna design; however, in this case mean BER inaccuracy is further reduced to 1%.


Introduction
The rapid growth in the demand for high data rate services in wireless networks over the last decades has boosted scientific research towards the implementation and evaluation of efficient transmission strategies, as well as the reduction of overall transceiver complexity.To this end, in 1998 the Wideband Code Division Multiple Access (WCDMA) physical layer protocol was adopted for the third generation (3G) wireless networks [1].The key idea was to allow mobile users to transmit simultaneously, occupying the whole transmission bandwidth, by assigning different spreading codes per user during transmission phase.Proper code design should ideally result in minimum auto-correlation and cross-correlation properties for any pair of the considered sequences.Therefore, both Multiple Access Interference (MAI) as well as Intersymbol Interference (ISI) can be significantly mitigated in multiuser scenarios [2].In addition, Bandwidth on Demand (BoD) can also be supported, either by assigning multiple codes per user or by adjusting the corresponding spreading factor.
Although 3G networks can provide bit rates per user that reach up to 2 Mbps, it soon became apparent that even higher transmission rates should be supported, as in the era of broadband wireless networks various bandwidth demanding services (e.g., video streaming) are quite popular among mobile subscribers.In this context, a practical solution with no additional spectrum requirements is the deployment of multiple antennas at both ends of the wireless transceiver, also known as Multiple Input Multiple Output (MIMO) [3].The increased degrees of freedom offered by such topologies, can be used either to improve mean Bit Error Rate (BER) or overall network throughput [4].In the first case, also known as diversity combining (DC) [4], link outage is significantly minimized, since there are multiple independent replicas of the transmitted signal, especially in rich scattering environments.In the second case, also known as spatial multiplexing (SM), independent data streams are transmitted from different antennas.Therefore, bit rate per user can be increased, at the cost however of increased transceiver complexity, since the diversity order of the system is reduced by a factor equal to the number of transmit antennas used for independent stream transmission.
In realistic wireless networks however, the existence of both MAI and ISI can significantly degrade the performance of MIMO transceivers [5].Therefore, advanced signal processing techniques should be employed at the receiver, in order to maximize DC or SM gain.To this end, various studies in literature examine appropriate transceiver techniques for MIMO orientations, as well as practical limitations from the deployment of MIMO architecture in next generation broadband wireless networks.In this context, in [5] the outage probability and capacity of MIMO-CDMA cellular systems is derived.According to the presented results, intercell interference can significantly degrade the performance of the MIMO orientation in SM scenarios.In [6], the performance of various transmission techniques in MIMO-WCDMA networks is investigated, with and without channel knowledge at the transmitter.As results indicate, adjacent sector cooperation can significantly improve the performance of the cellular orientation, at the cost however of increased signal burden for channel estimation and feedback to the transmitter.In [7], the performance of various precoding matrices is analyzed, when receiver operations are based on Minimum Mean Square Error (MMSE) detection.As results indicate, maximum likelihood detector can improve overall performance when the spatial correlation among the two input streams is taken into account.In [8], the performance of various beamforming algorithms for CDMA systems was evaluated, in terms of algorithmic complexity.In [9], the performance of MIMO-WCDMA downlink systems is investigated with antenna selection at the transmitter in correlated Nakagami fading channels.In [10], the multiuser case is examined in MIMO-WCDMA/HSDPA systems, where an effective user scheduling technique is proposed.In [11], mean BER is improved by applying initial-phase optimized Gold codes in a quasi-static fading channel.In [12], the performance of various beamforming algorithms for CDMA systems is investigated, in terms of achievable algorithmic complexity reduction.
In this work, the performance of a novel approach for MIMO-WCDMA networks is investigated, that is based on eigenvalue decomposition of the received signal matrix.This technique is also known as Principal Component Analysis (PCA) [13].PCA has been used in various scientific applications (e.g., [14,15]), where the main goal is to reduce signal dimension and consequently algorithmic complexity for signal processing and analysis.Therefore, in this study our main purpose is to reduce signal burden by sending to the transmitter only the primary eigenvector of the corresponding covariance matrix.As in [6], all multipath components contribute to signal estimation at the receiver.This work extends the study initially presented by the authors in the conference paper in [16], as the proposed approached is also evaluated considering a realistic four-element transmit antenna.
In the next section, signal transmission and reception in MIMO-WCDMA networks is described, in terms of equivalent signal matrices.In addition, a theoretical analysis for mean BER estimation is provided as well.In Section 3, the proposed approach based on PCA is analyzed, while results are presented in Section 4. Proposals for future work along with concluding remarks are provided in Section 5.

Signal Model and Receiver Operations
In all examined MIMO orientations, it is assumed that downlink transmission takes place and each transmitter/receiver is equipped with M t /M r antennas, respectively.The transmitted signal for the kth user (1≤ k ≤ K) can be expressed as a product of individual terms [1]: where p t,k is the total power allocated to the kth user (1 is the spreading code and t k the dim × dim transmission matrix, where dim = min(M t ,M r ) represents the dimension of the MIMO orientation (i.e., the number of operating transmit or receive antennas for a square channel matrix).It is assumed that the transmitted signal undergoes multiple reflections from L dominant clusters [17].Therefore, received signal can be expressed as: where τ l and H k,l,T k are the delay and channel coefficients, respectively, of the lth multipath component and n k is the corresponding Gaussian M r × 1 noise.Each element of H k,l,T k is assumed to be normally distributed with zero mean value and unity variance.The received signal is processed by a 2-D RAKE receiver, where the time of arrival of each multipath component along with the corresponding channel coefficient are estimated.Hence, the Maximal Ratio Combing (MRC) multiplying vectors per RAKE finger [18] can be formulated as: where the conjugate transpose of a matrix A is denoted by A H .The output dim × 1 signal of the lth RAKE finger of the kth user after correlation with the desired spreading sequence is given by: where ρ k,k,|l−l | and ρ k,k,|l−l | are the partial cross correlations [19], b k,0 is the symbol matrix at the current transmission period and b k,−1 is the corresponding transmission matrix at the previous symbol period.The partial cross correlations are given by: where x is the conjugate of x.In Equations ( 2)-( 4), the term T k indicates the selected set of transmit antennas per active user, given defined by: Electronics 2018, 7, 239 4 of 17 In Equation (7), TR is the set of all possible combinations of dim transmit antennas out of M t and λ m (A) is the maximum eigenvalue of matrix A. Considering for example a 4 × 2 MIMO orientation (i.e., M t = 4 and M r = 2), then the set TR will be equal to {(1,2), (1,3), (1,4), (2,3), (2,4), (3,4)}.
Under the assumption of independent symbol matrices among the various users of the MIMO-WCDMA orientation, the signal power of the desired user, ISI and MAI power as well as total noise power will be given by Equations ( 8)- (11), respectively [16]: where N o is the thermal noise level and A F the Frobenius norm of A. All output expressions from Equations ( 8)- (11) are dim × dim matrices.Therefore, all transmission matrices should be appropriately selected so as to provide diagonal expressions for each P s,k matrix.

Theoretical Background for the Calculation of Mean BER in MIMO-WCDMATtopologies
In this section, the Gaussian approach is described for the calculation of mean BER in a MIMO-WCDMA network.In this context, an arbitrary user (e.g., the kth) is considered as the desired one, and MAI from the rest K-1 users as well as ISI are expressed in terms of the desired user signal.With respect to Equation (10), then [19]: The term {A,x} indicates the xth line of matrixA.In addition: The Signal to Interference plus Noise Ratio (SINR) for the kth user will be given by [19]: Electronics 2018, 7, 239 where PG is the system's processing gain [1].In Equations ( 13) and ( 14), p k,d is the transmission power per active user and mode, subject to: The desired user signal can be expressed as: Hence, for BPSK modulation, BER will be given by [18]: where ) is the probability density function (pdf) of the desired user signal of Equation ( 16), while the Q function is given by: The pdf of the desired user signal can be evaluated by taking samples from a sufficient number of independent channel realizations.Afterwards, γ k,d is calculated for each transmission mode.For a 4 × 2 MIMO network with six multipath components, the pdf functions of the two transmission modes follow the Gamma distribution (Γ(x) is the Gamma function): where x = γ k,d .Note that the expressions in Equations ( 19) and ( 20) were derived after 10 7 channel realizations with the help of Matlab©.

Reducing Transmission Complexity in MIMO-WCDMA Networks Using PCA
As it was shown in [6], the selection and processing of all multipath components in frequency-selective MIMO-WCDMA networks can improve mean BER.In this context, an iterative transmission approach was proposed that calculates transmission matrices in Equation ( 8).This algorithm, which is described in Algorithm 1 for the kth user, successively calculates the RAKE multiplying vectors for each multipath component, the total signal power from Equation (8) as well as the transmit weight vector of the desired user, until convergence is achieved.Note that I dim is the dim × dim identity matrix, tr(X) is the trace of matrix X and D is a diagonal matrix indicating the power ratio per transmission mode.However, it is evident from Figure 1, where the pdf of the number of repeats of the Maximization of Signal-to-Noise Ratio (MSNR) algorithm of Algorithm 1 is presented (i.e., mean value of parameter n that indicates total repetitions for convergence), that the mean value can be up to 285 repeats.Note that the curve of Figure 1 is depicted for quanta of L 2 , since the L 2 multip−cation steps to calculate matrix A in Algorithm 1 are also included in the total number of repeats of the MSNR algorithm.

Algorithm 1 Maximization of Signal-to-Noise Ratio (MSNR) Algorithm
Step Step 2: since the L 2 multiplication steps to calculate matrix A in Algorithm 1 are also included in the total number of repeats of the MSNR algorithm.

Algorithm 1 Maximization of Signal-to-Noise Ratio (MSNR) Algorithm
Step Step 2: In our proposed approach, the correlation matrix of the received data is first formulated, with dimensions L × L and sorted afterwards.The (i,j) entry of this matrix indicates the degree of correlation among the signals of two different multipaths.Therefore, instead of considering all L 2 terms in matrix A formulation, C terms are selected that correspond to the first C dominant terms of the correlation matrix.In this approach, it is assumed that all received data for an arbitrary user after RAKE processing can be stacked as an L × S matrix (denoted as Z throughout the rest of the manuscript), where S indicates the samples per observation [13].The received data matrix per transmission mode will be given by: In Equation ( 21), zk,d,l,s is the s th received symbol from the k th user in the l th RAKE finger and d th transmission mode (S denotes the total number of transmitted symbols).Hence, the correlation matrix will be given by:  In our proposed approach, the correlation matrix of the received data is first formulated, with dimensions L × L and sorted afterwards.The (i,j) entry of this matrix indicates the degree of correlation among the signals of two different multipaths.Therefore, instead of considering all L 2 terms in matrix A formulation, C terms are selected that correspond to the first C dominant terms of the correlation matrix.In this approach, it is assumed that all received data for an arbitrary user after RAKE processing can be stacked as an L × S matrix (denoted as Z throughout the rest of the manuscript), where S indicates the samples per observation [13].The received data matrix per transmission mode will be given by: In Equation ( 21), z k,d,l,s is the sth received symbol from the kth user in the lth RAKE finger and dth transmission mode (S denotes the total number of transmitted symbols).Hence, the correlation matrix will be given by: With the help of PCA, a different representation of Z k,d can be formulated (i.e., Y k,d = PZ k,d ), so that the new correlation matrix, namely will be a diagonal one.In this case, the final set of data can be decomposed to individual orthogonal representations, also known as principal components.The importance of each principal direction is now measurable.Namely, the variances associated to each direction (i.e., rows of P) quantify how "principal" each direction is by rank-ordering each basis vector according to the corresponding variances.The reconstruction of matrix C Z k,d is performed according to the following equation: where w (1 ≤ w ≤ L) is the number of dimensions that are considered at matrix reconstruction, while V and Σ correspond to the eigenvector and diagonal eigenvalue matrix, respectively, of C Z k,d .Notation a:b indicates all elements from a to b with unity step.Hence, matrix C Z k,d can now be sorted, where the indexes of each sorted value are stored in the two-dimensional matrix B. Therefore, the modified matrix A can be alternately written as: In this case, the C (1 ≤ C ≤ L 2 ) most important terms can be taken into account during matrix A formulation.In the results that will be presented in the following section, C is chosen to be an even number, since This technique, which will be referred as MSNR-PCA algorithm throughout the rest of the paper, is described in Algorithm 2. Note that where for example the kth pair of rows corresponds to the (i,j) indexes of the kth maximum value of C Z k .
Algorithm 2 Maximization of SNR Algorithm based on Principal Component Analysis (PCA) Step 1:

PCA with Theoretical Antenna Radiation Patterns
In order to evaluate the performance of the proposed approach, Monte Carlo (MC) simulations were performed in various MIMO topologies.In the first set of results (i.e., Figures 2-4), the goal was to calculate mean BER and compare output values with the ones derived from the theoretical model as described in Section 2.2.For this reason, three different MIMO orientations (i.e., 4 × 2, 4 × 3 and 8 × 4) were considered, with two and six resolvable multipath components.The equivalent SNR was set to 10 dB, while a uniform power delay profile regarding the individual powers of the multipath components was taken into account.In figure legend, notation (M t ,M r ,L,PG,SNR,Sim or Th) has been carried out.Assuming independent data streams per different transmit antenna, then user throughput is set to R × min(M t ,M r ) Kbps, where R is the basic transmission rate of 30 Kbps that corresponds to the ratio of the system's bandwidth (3.84 MHz for the WCDMA, [1]) to the processing gain.Moreover, it is assumed that each mobile terminal is equipped with a number of RAKE fingers that equals the number of multipath components.Finally, the number of active users may vary from 1 to 40 with a step of 5.As it is evident from Figures 2-4, there is a good agreement among the Gaussian approximation for the error probability and consequently the mean BER and BER results derived from MC simulations.In addition, the selection of all multipath components (i.e., 6) leads to improved BER compared to the case of two multipath components.As expected, this gain reduces for high values of overall throughput, since to the increased number of active users, MAI can be significantly increased.
to 10 dB, while a uniform power delay profile regarding the individual powers of the multipath components was taken into account.In figure legend, notation (Mt,Mr,L,PG,SNR,Sim or Th) has been carried out.Assuming independent data streams per different transmit antenna, then user throughput is set to R × min(Mt,Mr) Kbps, where R is the basic transmission rate of 30 Kbps that corresponds to the ratio of the system's bandwidth (3.84 MHz for the WCDMA, [1]) to the processing gain.Moreover, it is assumed that each mobile terminal is equipped with a number of RAKE fingers that equals the number of multipath components.Finally, the number of active users may vary from 1 to 40 with a step of 5.As it is evident from Figures 2-4, there is a good agreement among the Gaussian approximation for the error probability and consequently the mean BER and BER results derived from MC simulations.In addition, the selection of all multipath components (i.e., 6) leads to improved BER compared to the case of two multipath components.As expected, this gain reduces for high values of overall throughput, since to the increased number of active users, MAI can be significantly increased.In Table 1, the mean number of repeats of the proposed MSNR-PCA algorithm that are required for convergence is presented, versus the number of power terms that are taken into account i.e., parameters C in Equation ( 24)).Moreover, in Figures 5 and 6 the mean BER is presented for a 4 × 2 MIMO orientation, five active users and SNR equal to 5 dB and 10 dB, respectively, versus the number of power terms.Two simulation cases have been considered: In the first one, correlation matrix per user is formulated considering only the primary eigenvector and eigenvalue from PCA.In the second case, all eigenvalues and eigenvectors are taken into account.In both cases, PG can be either 32 (high data rate services) or 128 (low data rate services).As derived from Figure 5, for high number of power terms there are practically no BER deviations per group of curves.However, even for reduced number of power terms, BER deviations are kept to a minimum.In particular, for 36 power terms, the mean BER is 7.9 × 10 −3 .For one eigenvector in PCA and 18 power terms, the mean BER is 8.1 × 10 −3 .Hence, in this case there is an error of almost 2.6%.However, it follows from Table 1 that there is also a complexity reduction gain of almost 10%, as the mean number of repeats is 108 for 36 power terms and reduces to 98 repeats for 18 power terms.Note that in the corresponding cases for PG equal to 32, there are practically no deviations among the two calculated BER values (1 and 6 eigenvalues and eigenvectors, respectively).Finally, in Figure 6, where SNR is set to 10 dB, convergence is achieved for higher number of power terms compared to the previous case of 5 dB SNR.In particular, for PG equal to 32, there are practically no BER deviations among 24 and 36 power terms.However, for PG equal to 128, convergence is achieved for 30 power terms.for PG equal to 32, there are practically no deviations among the two calculated BER values (1 and 6 eigenvalues and eigenvectors, respectively).Finally, in Figure 6, where SNR is set to 10 dB, convergence is achieved for higher number of power terms compared to the previous case of 5 dB SNR.In particular, for PG equal to 32, there are practically no BER deviations among 24 and 36 power terms.However, for PG equal to 128, convergence is achieved for 30 power terms.

PCA with Antenna Radiation Patterns from Electromagnetic Simulations
In the final set of simulations, a realistic four-element transmit antenna has been considered (Figure 7), as well as two omni antennas per receiver.In this case, the channel between an arbitrary transmit antenna (denoted as q) and an arbitrary receive antenna (denoted as u) for the lth multipath component is given according to 3GPP specifications [17], as: where j is the imaginary unit, P l is the power of the lth path, σ SF is the lognormal shadow fading, M is the number of sub-paths per path, θ l,m,AoD and θ l,m,AoA are the angles of departure (AoD) and arrival (AoA), respectively, for the mth subpath of the lth path, G BS (θ l,m,AoD ) is the transmitter antenna gain for each array element and G MS (θ l,m,AoA ) is the receiver antenna gain for each array element for the AoD and AoA, respectively.Moreover, Φ l,m is the phase of the mth subpath of the lth path, uniformly distributed over [0, 2π].Finally, k w is the wave number 2π/λ, where λ is the carrier wavelength in meters, d q is the distance in meters of the transmit antenna element q from the reference (q = 1) antenna and d u and is the distance in meters of the receive antenna element u from the reference (u = 1) antenna.
Electronics 2018, 7, x FOR PEER REVIEW 12 of 18 antenna and du and is the distance in meters of the receive antenna element u from the reference (u = 1) antenna.The transmit antenna array consists of two elements, each one composed of two similar co-planar elements positioned in close proximity to each other and rotated one to the other, by an angle of 90°, in order to obtain polarization diversity.The elements are fractal shaped, via a modified Hilbert curve configuration and were designed by the technique proposed in [20], the target being the dual frequency operation at ISM (2.4-2.489GHz) and 5-6 GHz bands.The elements are semi-printed on the dielectric substrate RO3003 (εr = 3 and thickness h = 0.812 mm).The ground plane, having the shape shown in Figure 7, is on the back side of the substrate, extended only at a part of the whole antenna strip.This configuration, as proposed in [20], is not enough to ensure low values of all the correlation coefficients among the array's elements when two dual elements are arranged to compose a four-element array, because the four elements have (two by two) the same polarization.The problem was confronted embedding an EBG lattice as shown in Figure 7. EBG lattices are commonly used in antenna structures as they can enhance their performance, increasing The transmit antenna array consists of two elements, each one composed of two similar co-planar elements positioned in close proximity to each other and rotated one to the other, by an angle of 90 • , in order to obtain polarization diversity.The elements are fractal shaped, via a modified Hilbert curve configuration and were designed by the technique proposed in [20], the target being the dual frequency operation at ISM (2.4-2.489GHz) and 5-6 GHz bands.The elements are semi-printed on the dielectric substrate RO3003 (ε r = 3 and thickness h = 0.812 mm).The ground plane, having the shape shown in Figure 7, is on the back side of the substrate, extended only at a part of the whole antenna strip.This configuration, as proposed in [20], is not enough to ensure low values of all the correlation coefficients among the array's elements when two dual elements are arranged to compose a four-element array, because the four elements have (two by two) the same polarization.The problem was confronted embedding an EBG lattice as shown in Figure 7. EBG lattices are commonly used in antenna structures as they can enhance their performance, increasing their frequency bandwidth, contributing to the decrease of their size and affecting their gain either increasing it or keeping its value constant within wide frequency ranges [21][22][23].Besides these ways of influence, to the antennas' operation, the incorporation of EBGs, of mushroom or of other type in antenna arrays, can reduce significantly the isolation among their element reducing their electromagnetic mutual coupling.This result is very useful at MIMO antennas systems as it contributes to the reduction of correlation coefficients among their elements, thus enhancing the MIMO procedure ( [24][25][26][27][28]).In the present work, the performance of the prototype of Figure 7 is presented via Figures 8-11. Figure 8 depicts the scattering coefficients S ii at the input of the antenna elements, confirming that they operate at the ISM band, which is under interest.Figure 9 illustrates the mutual signal scattering coefficients (S ij ) among the elements and, for comparison, the respective coefficients of the same array without EBG are presented.It is observed that significant reduction, more than 5 dB, of the level of those of the S ij coefficients which correspond to pairs of elements that have the EBG between them, as the S 23 , S 24 and S 14 can be achieved.Moreover, indicative results of the 3D and polar gain radiation patterns for one of the elements are shown in Figure 10.The patterns of all four elements are almost similar and the deviations from their maximum values do not exceed 3 dB inside the space angular section of utilization.The suitable orientation of the elements and the reduction of their mutual coupling lead to the reduction of the correlation of the respective coefficients, ρ ij .It is validated by calculating them via the theory of [29] and Figure 11 depicts the results of these calculations for the cases of the array with EBG and without EBG.The comparison between the two arrays, shows a significant reduction of the ρ ij at the array with EBG especially when the elevation angle of the incoming signals are high, namely when they come from the space area at the broadside of the antenna.Among all four elements, of the EBG array, ρ ij values are less than 0.1, being a little higher in some cases.
Electronics 2018, 7, x FOR PEER REVIEW 13 of 18 area at the broadside of the antenna.Among all four elements, of the EBG array, ρij values are less than 0.1, being a little higher in some cases.Frequency In Figure 12, results are provided for two cases, considering a 4 × 2 MIMO orientation, as well as one and six eigenvector and eigenvalues for covariance matrix formulation, respectively.They concern to two antennas.The first of them, Case A, represents an ideal system with an antenna having omni directional radiation pattern and channels modeled according to the Gaussian distribution which corresponds to zero correlation.The second one, Case B, is the realistic antenna array, of Figure 7, designed properly in order the correlation to be minimized.In this case, the channels were modeled according to 3GPP specifications as described in Section 4 and Equation (25).In this equation, for the evaluation of the channel matrix elements, the antenna's radiation patterns are taken into account through the antennas' gain values which are calculated via simulation.The very low space correlation is impressed to these patterns and is inserted substantially to the channel estimation.As it can be observed, in Case B, for 18 power terms the mean BER is 3.142 × 10 −4 , while for 36 power terms the corresponding number is 3.109 × 10 −4 .Therefore, although we have again a complexity reduction equal to 10% as in Case A, the corresponding mean BER inaccuracy is further reduced to 1%.It is also shown that the results of BER for the antennas of both cases are very low and converge to each other proving that a real antenna when properly designed can achieve the records of a theoretically omni directional antenna which operates under ideal conditions.In Figure 12, results are provided for two cases, considering a 4 × 2 MIMO orientation, as well as one and six eigenvector and eigenvalues for covariance matrix formulation, respectively.They concern to two antennas.The first of them, Case A, represents an ideal system with an antenna having omni directional radiation pattern and channels modeled according to the Gaussian distribution which corresponds to zero correlation.The second one, Case B, is the realistic antenna array, of Figure 7, designed properly in order the correlation to be minimized.In this case, the channels were modeled according to 3GPP specifications as described in Section 4 and Equation (25).In this equation, for the evaluation of the channel matrix elements, the antenna's radiation patterns are taken into account through the antennas' gain values which are calculated via simulation.The very low space correlation is impressed to these patterns and is inserted substantially to the channel estimation.As it can be observed, in Case B, for 18 power terms the mean BER is 3.142 × 10 −4 , while for 36 power terms the corresponding number is 3.109 × 10 −4 .Therefore, although we have again a complexity reduction equal to 10% as in Case A, the corresponding mean BER inaccuracy is further reduced to 1%.It is also shown that the results of BER for the antennas of both cases are very low and converge to each other proving that a real antenna when properly designed can achieve the records of a theoretically omni directional antenna which operates under ideal conditions.

Conclusions
The performance of a proposed transmission strategy for MIMO-WCDMA networks based on eigenvalue analysis was evaluated in this work.The goal was to select the appropriate transmission matrices per active user in order to maximize system throughput in a multiuser/multipath environment.As it was shown in the results section, the selection and processing of all multipath components via RAKE receivers can improve system performance in terms of mean BER.Moreover, the adoption of the proposed approach can reduce the overall complexity of transmit weight formulation.Results were also validated for the case of a realistic four-element antenna.In this case, correlation coefficients values among all four elements were less than 0.1.Hence, BER deterioration was only 10% compared to the ideal case of uncorrelated channels among all transmit and receive antennas.
Motivated by the above results, we aim that future work will include the extension of the presented results in cellular networks, where fundamental trade-offs among spatial multiplexing and diversity transmission mode will be investigated.

Conflicts of Interest:
The authors declare no conflicts of interest.

Conclusions
The performance of a proposed transmission strategy for MIMO-WCDMA networks based on eigenvalue analysis was evaluated in this work.The goal was to select the appropriate transmission matrices per active user in order to maximize system throughput in a multiuser/multipath environment.As it was shown in the results section, the selection and processing of all multipath components via RAKE receivers can improve system performance in terms of mean BER.Moreover, the adoption of the proposed approach can reduce the overall complexity of transmit weight formulation.Results were also validated for the case of a realistic four-element antenna.In this case, correlation coefficients values among all four elements were less than 0.1.Hence, BER deterioration was only 10% compared to the ideal case of uncorrelated channels among all transmit and receive antennas.
Motivated by the above results, we aim that future work will include the extension of the presented results in cellular networks, where fundamental trade-offs among spatial multiplexing and diversity transmission mode will be investigated.

Figure 1 .
Figure 1.Probability density function (PDF) curve of the number of repeats of the Maximization of Signal to Noise Ratio (MSNR) algorithm (Algorithm 1).
of repeats of the MSNR algorithm

Figure 1 .
Figure 1.Probability density function (PDF) curve of the number of repeats of the Maximization of Signal to Noise Ratio (MSNR) algorithm (Algorithm 1).

6 :
If tr P k,n+1 − tr P k,n ≤ ε•tr P k,n algorithm terminates, else go to Step 4. Then go to Step 2.

Figure 8 .
Figure 8. Scattering coefficients for the signals, at the input of the array elements, versus frequency.

Figure 8 .
Figure 8. Scattering coefficients for the signals, at the input of the array elements, versus frequency.

Figure 8 .
Figure 8. Scattering coefficients for the signals, at the input of the array elements, versus frequency.

Figure 9 .
Figure 9. Mutual scattering coefficients of the signals, among the inputs of the array elements, versus frequency (With EBG in (a) and without EBG in (b)).

Frequency 24 Figure 9 .Figure 9 .Figure 10 .
Figure 9. Mutual scattering coefficients of the signals, among the inputs of the array elements, versus frequency (With EBG in (a) and without EBG in (b)).

24 Figure 11 .
Figure 11.Correlation coefficients of the array with and without EBG, versus the mean direction of arrival of the waves incoming from the two main planes of the array (f = 2.45 GHz).

Figure 12 .
Figure 12.BER comparisons for a MIMO-WCDMA orientation for varying number of eigenvector/ eigenvalues for signal reconstruction (SNR = 10 dB, PG = 32, Cases A and B).

Table 1 .
Mean number of repeats of the MSNR-PCA algorithm.

Table 1 .
Mean number of repeats of the MSNR-PCA algorithm.