Interference-free SDMA for FBMC-OQAM

Filter-bank multi-carrier (FBMC) modulations have recently been considered for the emerging wireless communication systems as a means to improve the utilization of the physical resources and the robustness to channel time variations. FBMC divides the overall frequency channel in a set of subchannels of bandwidth proportionally decreasing with the number of subchannels. If the number of subchannels is high enough, the bandwidth of each subchannel is small enough to assume that it is approximately flat. On the other hand, space-division multiple access (SDMA) is a recognized technique to support multiple access in the downlink of a multi-user system. The user signals are precoded at the base station equipped with multiple antennas to separate the users in the spatial domain. The application of SDMA to FBMC is unfortunately difficult when the channel is too frequency selective (or when the number of subchannels to too small) to assume flat subchannels. In that case, the system suffers from inter-symbol and inter-subchannel interference, besides the multi-user interference inherent to SDMA. State-of-the art solutions simply neglect the inter-symbol/subchannel interference. This article proposes a new SDMA precoder for FBMC capable of mitigating the three sources of interference. It is constructed per subchannel in order to keep an acceptable complexity and has the structure of a filter applied on each subchannel and its neighbors at twice the symbol rate. Numerical results demonstrate that the precoder can get rid of all the interference present in the system and benefit therefore from the diversity and power gains achievable with multiple antenna systems.


Introduction
The emerging wireless communication systems are carefully designed to simultaneously optimize the delivered user capacity, the average spectral efficiency and the cell coverage.Nowadays most of the emerging communications systems rely on the orthogonal frequencydivision multiplexing (OFDM) modulation technique [1].By adding a redundant cyclic prefix to the symbol sequences, OFDM creates orthogonal and flat subchannels in the frequency domain that can be equalized at a low complexity by single coefficient inversion.The OFDM modulation has also been selected because it is easily combined with multi-input multi-output (MIMO) techniques exploiting the extra spatial dimension besides the time and frequency dimensions to increase the link reliability and capacity [2,3].The MIMO techniques are applied on each subchannel independently in order to keep a low computational complexity.
Filter-bank multi-carrier (FBMC) modulations, and more specifically FBMC-offset quadrature amplitude modulation (OQAM), are seen as an interesting alternative to OFDM for future wireless communication systems [4,5].The time/frequency resolution of the waveforms is increased and can be traded-off resulting in a better utilization of the physical resources and potentially in an improved robustness to time-variant channel characteristics and carrier frequency offsets.Like OFDM, FBMC-OQAM decomposes the communication channel in a set of lower-bandwidth subchannels that can therefore also be compensated at a low complexity with a single-tap equalizer.Contrary to OFDM, FBMC-OQAM does not require the addition of a cyclic prefix and the created subchannels are only approximately flat and orthogonal.When the channel frequency selectivity increases, the FBMC-OQAM system suffers from both inter-subchannel interference and inter-symbol interference on each subchannel, making it necessary to use advanced equalizer structures [6].Furthermore the combination of FBMC-OQAM with MIMO techniques results in an unmanageable interference term appearing between the antenna streams on adjacent subchannels, that makes the design of the system challenging.The design of MIMO FBMC-OQAM systems has triggered a lot of research recently.
The FBMC-OQAM system can first be combined with spatial division multiplexing (SDM) used to increase the link capacity.At the transmitter, different symbol streams are sent independently on the subchannels of the different antennas.Article [7] proposes to compensate for the inter-subchannel interference generated by the multiantenna streams at the receiver with an equalizer working per subchannel at a fractional symbol rate.While [7] focuses on linear and successive interference cancelation equalizers, article [8] extends the results to the adaptive decision feedback equalizer.To better cope with intersubchannel interference in the SDM system, a modified symbol mapping is also proposed in [9,10] and combined with maximum-likelihood (ML) detection.
The FBMC-OQAM is secondly combined with the space-time block codes (and specifically with the Alamouti scheme) and with the space-time trellis codes to improve the communication reliability.Article [11] proposes to apply the Alamouti scheme designed initially for a flat fading channel to each subchannel independently at the transmitter and to handle the inter-subchannel interference generated by the multi-antenna streams at the receiver with an additional equalizer.Article [12] rather combines FBMC-OQAM with code-division multiple access (CDMA) to restore the complex orthogonality among the subchannels and be therefore able to apply and decode the Alamouti scheme independently per-subchannel.Finally [13] applies the space-time trellis codes to FBMC-OQAM and combines it with iterative decoding at the receiver.
In the article [14] and in the report of the European project Phydyas [15], it is proposed to perform transmit beamforming by multiplying each subchannel at the transmitter with a complex coefficient.An equalizer is further foreseen at the receiver to compensate for the inter-symbol and inter-subchannel interferences.
In systems where the base station is equipped with multiple antennas, it is also possible to provide multiple access in the spatial domain.This multi-user technique is referred to as space-division multiple access (SDMA).In line-of-sight environments, SDMA can be conceived as a form of beamsteering where different beams are aimed at different users, thereby separating their signals.In non-line-of-sight environments, wideband beamforming is required.Combined with OFDM, wideband SDMA reduces to narrowband beamforming per-subchannel [16,17].Combined with FBMC-OQAM, wideband SDMA can also be implemented with a singletap beamformer per-subchannel by neglecting the intersymbol and inter-subchannel interferences.However, the level of remaining interference can significantly reduce the performance if an additional equalizer structure is not present.
In this article, we design a wideband SDMA linear precoder that separates the users spatially in the downlink of an FBMC-OQAM system.Compared with the existing literature, the main contribution of this article is to propose a structure that is able to suppress all sources of interference present in the FBMC-OQAM system: inter-symbol, inter-subchannel, and multi-user interference; even when the subchannels are frequency selective (highly selective channels, low number of subchannels).In order to do so, we present a structure where the precoder is a filter applied at twice the symbol rate, on each subchannel and its two adjacent subchannels.The computational complexity of the resulting precoder is linearly increasing with the number of subchannels and is thus still moderate.The SDMA precoder is built based on the knowledge of the channel impulse responses according to the zero-forcing criterion (inversion of the frequency selective channel).
The article is organized as follows.Section 2 introduces the model of the FBMC-OQAM system.The SDMA precoder is built in Section 3 based on a unified system model, detailing the structure of the precoder and showing how it can be computed in practice.Finally Section 4 demonstrates the efficiency of the scheme numerically.In the sequel, we use bold lowercase letters for the vectors and bold uppercase letters for the matrices.The operators (.) T and (.) H are used to denote the transpose and conjugate transpose of a vector or a matrix.Matrix 0 M,N is the size-M × N matrix composed of zeros, matrix I N is the size-N identity matrix.

FBMC-OQAM system
The article investigates the downlink of a SDMA system: a multi-antenna base station transmits independent signals to multiple users.The number of antennas N T at the base station is assumed equal to or larger than the number of users N R .To generate the signal at the antenna n T , the base station implements a synthesis filterbank composed of M subchannels.Figure 1 details the model of the SDMA FBMC-OQAM system, focusing on antenna n T at the base station (n T = 1, . . ., N T ), on receiver n R (n R = 1, . . ., N R ), and on subchannel k and its two adjacent subchannels k − 1 and k + 1 (k = −M/2, . . ., M/2 − 1).
In the conventional FBMC-OQAM system, the sequence of QAM information symbols s (n T ,k) [ l] transmitted on each subchannel k at the symbol rate R s = 1/T s (T s is the QAM symbol duration) is separated into its real http://asp.eurasipjournals.com/content/2013/1/46 To implement the OQAM modulation, the sequence k) [ m] is delayed by one sample (half a symbol), multiplied by j, and summed with the sequence s R (n T ,k) [ m].In this article, we design a SDMA precoder that generates the sequences s R (n T ,k) [ m] and s I (n T ,k) [ m] directly at its output.Because, the sequences s R (n T ,k) [ m] and s I (n T ,k) [ m] are defined at twice the symbol rate, the SDMA prefilter is built on the frequency band [ −1/T s , 1/T s ].
After upsampling by M/2 to reach the sampling rate 1/MT s , the symbols are convolved with the pulse shaping filter u [ n].Without loss of generality, we assume that the pulse shaping filter u[ n] is a halfroot Nyquist filter of cut-off frequency equal to the symbol rate.In practice, the pulse shaping filter is often causal and of finite length.The model focuses on the subchannel k and its two adjacent subchannels k − 1 and k + 1 as they are the only ones to impact the subchannel k.A factor ±j multiplies the signal on the odd subchannels.The M subchannels are added together after signal shifting on the frequency grid.Two adjacent subchannels are spaced by a frequency shift equal to the input symbol rate.
The sequence transmitted at antenna n T is therefore given by: ( Each user receives the sum of the signals transmitted from the N T base station antennas.The signal from antenna n T is received at the receiver n R after propagating through the frequency selective channel w corrupts the received sequence y n R [ n] at the receiver n R .At each receiver, the analysis filterbank implements the counterpart of each operation performed at the transmitter.After frequency shifting to process each subchannel around the origin, the received signal is convolved with the filters http://asp.eurasipjournals.com/content/2013/1/46matched to the pulse shaping filters.The half symbol delay applied on the imaginary branch is compensated so that both the real and imaginary signals are received synchronously when downsampling by M to reach back the symbol rate 1/T s .Only the real part on the real branch and the imaginary part on the imaginary branch are kept at the receiver. In the case of a single-antenna/single-user system (N T = N R = 1), no interference is left at the output of the receiver when the channels are flat or slightly frequency selective thanks to the OQAM modulation and to the perfect reconstruction property of the pulse shaping filter.Therefore, a single slightly frequency selective channel can be compensated at a low complexity by simple channel coefficient inversion (like in OFDM).In the case of a multi-antenna system communicating over highly frequency selective channels, the system unfortunately suffers from three sources of interference: (i) inter-symbol interference on each subchannel; (ii) interference between the adjacent subchannels; (iii) inter-antenna interference.The goal of the SDMA precoder is to compensate for the multi-antenna frequency selective channel by prefiltering the information symbols so that the three sources of interference disappear.

Equivalent system model
The real output sequences r R k) [ l] at the receiver n R on subchannel k can be expressed as a function of the real input sequences s R (n T ,i) [ m] and s I (n T ,i) [ m] coming from all transmit antennas (n T = 1, . . ., N T ) on the subchannel k and its adjacent subchannels (i = k−1, k, k+ 1).Since we are developing the precoder at twice the symbol rate, we need to define the polyphase components of the input sequences: for ρ = 0, 1 so that the overall model is eventually working at the symbol rate: where ⊗ is the convolution operator.The parameter ρ is relative to the two polyphase components.The composite channel impulse responses result from the convolution of the transmitter filter, channel response and receiver filter, downsampled by a factor M, as defined in: are the synthesis and analysis filters.They include the pulse shaping filter, the frequency shifts and the potential j factor: and: (15) http://asp.eurasipjournals.com/content/2013/1/46 The noise sequences k) [ l] are obtained by filtering the input noise sequence w n R [ n] with the analysis filters, downsampling by M, and keeping the real part of the result on the real branch or the imaginary part of the result on the imaginary branch.
It is noteworthy to observe by inspecting Figure 1 that: it is impossible to build a unified complex model equivalent to ( 6), (7).Therefore, we have to stick to the real convolution model ( 6), (7) to correctly represent the FBMC-OQAM system (and the precoder that will be derived in the next sections).
The equivalent system model is illustrated in the right part of Figure 2. In order to simplify the diagram, the real/imaginary parts of the input sequence polyphase components at all transmit antennas have been arranged in a sequence of length-4N T vectors (i = k − 1, k, k + 1): the real/imaginary parts of the noise and received sequences at all receive antennas have been arranged in sequences of length-2N R vectors: and the real/imaginary parts of the channel impulse responses corresponding to all transmit/receive anten- with: for X, Y = R, I. Expressions ( 6), ( 7) are equivalently written: where the definition of the convolution ⊗ of a vector sequence x[ l] with a matrix sequence G[ l] is traditionally defined as:

SDMA precoder
The SDMA precoder could be jointly designed for all subchannels taking their mutual impact into account.However this option would be of prohibitive computational complexity and it is in fact possible to build a simple zero-forcing precoder independently for each subchannel.It is based on the property that each subchannel only interferes with its two adjacent subchannels (it has no frequency overlap with the others).Then, each precoder for each subchannel is designed in such a way that the interference it causes on the other subchannels is canceled out, as well as the inter-antenna interference and the intersymbol interference inside its own subchannel.In order to have enough degrees of freedom, the precoder is built at twice the symbol rate.This way, the system is ensured to be free of interference.This section derives the SDMA precoder for subchannel k (k = −M/2, . . ., M/2 − 1).It is performed in two steps: we first detail the structure of the SDMA precoder for subchannel k based on the equivalent system model and, we secondly show how it can be computed in practice.

Structure of the SDMA precoder
The SDMA precoder for the subchannel k is illustrated in the left part of Figure 2. It precodes the real input symbol vector sequence: with the real matrix filters to the subchannels k − 1, k and k + 1, respectively, of the FBMC-OQAM system.The SDMA precoder is designed to fulfill the following conditions (zero-forcing criterion): • The input symbols on subchannel k do not interfere with the subchannels k − 2 and k + 2: • The input symbols on subchannel k do not interfere with the subchannels k − 1 and k + 1: • There is no inter-symbol/inter-antenna interference on subchannel k : where δ[ l] is the Dirac pulse.
It should be noted that the usual condition required to design the SDMA precoder for a flat fading channel as proposed in [14,15] is G k,k ⊗ P k,k = I 2N R that involves only the second term of (29) independent of the time index.The SDMA precoder for the FBMC-OQAM system has to fulfill the five conditions (25)-(29) together mostly because of the inter-subchannel interference inherent to the system.

Expression of the SDMA precoder
The overall system relative to subchannel k, including the SDMA precoder and the channel as illustrated in Figure 2, can be mathematically expressed in a unified matrix model: where the received and noise sequences on the different subchannels have been stacked in length-10N R received and noise vectors: and the size-12N T × 2N R precoder matrix is: Model (30) focuses on the contribution of the symbols transmitted on subchannel k to the received vectors on subchannels k − 2 to k + 2, and therefore neglects the contribution from the other subchannels.Since the SDMA precoder is designed to prevent interference from subchannel k on the subchannels k ± 1, k ± 2, model ( 30) is however sufficient.A closed-form expression of the zero-forcing (ZF) precoder, of infinite length in the general case, is conveniently given in the z-domain.Let G k [ z] denotes the z-transformed channel matrix.The ZF precoder matrix is expressed as [18]: where Matrix is composed of the 2N R central columns of the identity matrix I 10N R so that the part of the precoder relative to the information symbols on the subchannel k is selected.Among all precoders satisfying the zero forcing criterion (25)-(29), the precoder (35) taking the form of the channel matrix pseudo-inverse maximizes the sum capacity if a global power constraint is assumed at the transmitter [19].We normalize the precoder when its performance is assessed to satisfy the power constraint.The zero-forcing precoder applies a filter matched to the channel impulse responses preceded by the inversion of the overall channel responses to mitigate the interference.The matched filter can be seen as a wideband beamformer and the ZF precoder may hence bring the corresponding power and diversity gains (when there is only one receiver, it is equivalent to the maximum-ratio combining scheme that is well-known to benefit from the two types of gain).
Different methods exist to practically compute the ZF filter for a MIMO channel [18].In this article, we have followed the time-domain method first described in [20] that relies on the construction of a matrix model of the convo-lution MIMO channel.It is detailed in the Appendix 1 for the particular case considered in this article.The method always produces a finite impulse response solution of accuracy increasing with the matrix model dimension (the parameter N in the Appendix 1).The computation of the SDMA precoder requires the inversion of the channel correlation matrix which is always well-conditioned when the number of antennas at the transmitter N T is larger than or equal to the number of users N R .It is also interesting to note that the SDMA precoder can be developed at the symbol rate (instead of at twice the symbol rate) when N T ≥ 2N R .However the number of users has to be much smaller than the number of antennas at the base station in that case which drastically reduces the system capacity.

Performance and implementation complexity analysis
The performance and the implementation complexity of the proposed SDMA precoder is evaluated numerically.
The FBMC system is composed of a varying number of subchannels (M = 32, 64, 128), each modulated with offset quaternary phase shift keying (OQPSK) symbols shaped with a halfroot Nyquist filter.A 1 rolloff factor is considered and there is hence significant spectral overlap between the adjacent subchannels.The length of the halfroot Nyquist filters is fixed to 8 symbols.The sample rate at the output of the subchannels is equal to 64 Msamples/s, so that the symbol rate at the input of each subchannel is equal to 2, 1, and 0.5 Msymbols/s for M = 32, 64, and 128, respectively.The base station is equipped with 4 antennas and the number of receivers is varying from 1 to 4. A frequency selective channel of exponentially decreasing power delay profile is assumed.The channel delay spread is equal to 0.0525 μs (3.36 samples).The results are averaged over 1000 channel realizations that are normalized in power on the average.

Implementation complexity
The computational complexity required to implement the FBMC-OQAM system with SDMA is assessed and compared to the one required to implement the SDMA-OFDM system [16,17].A distinction must be made between the operations required to compute the precoder and those required to precode the sequence of symbols.
Regarding the complexity required to compute the precoder: • Since the precoder computation involves mainly the inversion of the size 12N R N real channel correlation matrix for each subchannel (see the Appendix 1 for the details) and is therefore O(N 3 R N 3 ), the parameter N is a direct measure of the computational complexity required for the computation of the precoder.
• In comparison, the SDMA-OFDM system requires the inversion of a size N R complex channel correlation matrix for each carrier and is therefore Regarding the precoding complexity: • The SDMA precoder on each subchannel can be seen as a size 12N T × 2N R matrix filter of length equal to N + L − 1, where L is the fixed length of the channel impulse responses (see the Appendix 1 for the details).Therefore, 24N T N R (N +L −1) multiply-andaccumulate (MAC) real operations are performed for each set of N R symbols on each subchannel.• In comparison, the SDMA-OFDM system requires the multiplication with a size N T × N R matrix on each subchannel so that N T N R MAC complex operations are performed for each set of N R symbols on each subchannel.
Therefore a special care needs to be taken in order to select the dimension of the matrix model used to compute the precoder (the parameter N).On the one hand, N should be as large as possible to ensure the accuracy of the precoder computation.On the other hand, N should be as small as possible to limit the computational complexity.Figure 3 illustrates the mean square error (MSE) as a function of the parameter N for a varying number of subchannels.The MSE is computed by assessing the error remaining between the overall channel impulse response (precoder and channel impulse response) and the dirac pulse a .The number of receivers has been fixed to 3. Figure 3 first confirms that the accuracy of the precoder computation improves when N is increased as the remaining MSE becomes negligible.The minimum value N necessary to achieve the desired MSE can therefore be inferred from the figure.By comparing the solid curves in Figure 3, it is secondly observed that the MSE is only slightly reduced when the number of subchannels increases as each subchannel becomes approximately frequency flat (no improvement is observed for M larger than 64).Therefore the computational complexity is approximately independent of the number of subchannels.The dashed curve in Figure 3 illustrates finally the MSE loss obtained when the roll-off (ro) factor of the halfroot Nyquist filters is equal to 0.5 (instead of 1).In that case, the dimension of the matrix model necessary to compute the SDMA precoder (and therefore the computational complexity) must significantly be increased to ensure the same accuracy.This is explained by the fact that the filters are less localized in the time domain.
Even if the complexity of the proposed SDMA FBMC-OQAM system is much larger than the complexity of the corresponding SDMA-OFDM system, it increases linearly with the number of subchannels and is therefore acceptable compared to the complexity that would be necessary when the precoder is jointly designed on all subchannels.
In the following performance results, the SDMA precoder is computed by building a dimension 20 time-domain matrix model to ensure the accuracy of the precoder computation.The length of the resulting precoder filters is approximately equal to 30.
Another important system design parameter is the peak-to-average-power ratio (PAPR) that measures the implementation feasibility of the power amplifiers at the transmitter.Figure 4 illustrates the cumulative distribution of the PAPR for the FBMC-OQAM system with and without SDMA precoding.The comparison is first made by considering only 1 transmit antenna and 1 user, and then by extending the system to 4 transmit antennas and 3 users.The PAPR is assessed by comparing the peak power to the average power when bursts of 10 symbols are transmitted.It is observed in both cases that the SDMA precoder tends to decrease slightly the PAPR by introducing correlation between the signals on the adjacent channels.

Performance analysis
Figure 5 compares the bit error rate (BER) performance of the SDMA FBMC-OQAM system proposed in this article (solid curve) to the one of the SDMA-OFDM system investigated in [16,17] (dashed curves).The BER performance is assessed as a function of the signal-to-noise power ratio (SNR) defined as the average transmitted power per antenna divided by the noise power in the system bandwidth (in practice, we normalize the precoder to fix the average transmitted power per antenna and vary the noise power).The number of receivers has been fixed to 3. The trade-off existing between the BER and the spectral efficiency of the OFDM system for a varying cyclic prefix (CP) is investigated.When the CP is null, the spectral efficiency of the OFDM system is maximum but there is significant inter-carrier interference (ICI), explaining why the curves tend asymptotically to a high BER value.When the CP is large (16 samples), there is no ICI anymore and the BER performance is similar to the one of the FBMC-OQAM system, but the spectral efficiency is reduced by 1/5 compared to the FBMC-OQAM system.A small power loss (1/5 of the power) is also observed compared to the FBMC-OQAM system because of the CP removal at the receiver, leading to a shift of the BER curve to the right.It is noteworthy that the ICI significantly limits the BER performance even when the CP is slightly larger than the channel delay spread (4 samples).Compared to the SDMA-OFDM system, the proposed SDMA FBMC-OQAM system is capable of delivering the maximum spectral efficiency without any BER penalty.
Figure 6 illustrates the BER for a varying number of subchannels.The number of receivers is also fixed to 3. The figure compares the performance of the proposed SDMA precoder (solid curves) to the one of the singletap SDMA precoder designed for the OFDM system, but used for the FBMC-OQAM system (dashed curves).Because the channel cannot be assumed to be flat on each subchannel, the single-tap SDMA precoder significantly suffers from the inter-symbol, inter-subchannel and inter- antenna interference.The proposed precoder on the other hand can successfully get rid of all the interference.It delivers therefore a performance approximately independent on the number of subchannels, on the contrary to the single-tap SDMA precoder of performance increasingly limited by the remaining interference when the number of subchannels is reduced.Figure 7 illustrates the BER for a varying number of receivers.The number of subchannels has been fixed to 64.The trade-off existing between the diversity and spatial multiplexing gains is observed by comparing the solid curves: when the number of receivers decreases in the system, the diversity gain progressively replaces the spatial multiplexing gain (the slope of the curves is steeper but the rate is decreased).Furthermore, the system progressively benefits from an improved power gain (the curves are shifted to the left).The dashed curve illustrates the performance obtained when the SDMA precoder is designed at the symbol rate (instead of at twice the symbol rate) for a number of receivers equal to 2 (so that N T ≥ 2N R ).Compared to the system employing a precoder designed at twice the symbol rate for N R = 2, the performance is degraded because the system does not benefit anymore from the additional degrees of freedom brought by the precoder over-sampling factor (OSF).

Conclusion
The objective of the article has been to investigate the combination of the filter-bank multi-carrier (FBMC) modulation, proposed recently for the wireless communications systems because of the improved time/frequency resolution of their waveforms (compared to orthogonal frequency-division multiplexing (OFDM)), with space division-multiple access (SDMA), used in the downlink to separate the users in the spatial domain.The SDMA precoder is built to cancel out the inter-symbol interference on each subchannel, the inter-subchannel interference and the inter-antenna interference.It has the structure of a set of filters used to precode the user streams of information symbols on each subchannel and its adjacent subchannels and on all antennas.Its complexity is therefore much lower than the overall precoder applied on all subchannels simultaneously.The efficiency of the designed scheme has been demonstrated numerically both in terms of multi-access ability and diversity/power gains.http://asp.eurasipjournals.com/content/2013/1/46equalizer of the left-right mirrored MIMO channel and to deduce the ZF precoder of the MIMO channel from it.The left-right mirrored FBMC-OQAM system is obtained by inverting the directions of the arrows in the right part of Figure 2. Neglecting the noise terms, the terms r k [ n] become the inputs and the terms s k [ n] become the outputs of the system.A quick analysis shows that the equations that the filters need to fulfill in order to perform ZF equalization for this left-right mirror systems are identical to the Equations ( 25)-(29) of the ZF precoder of the initial system, which guarantees that the correct solutions to our initial problem are found using this work-around.The matrix model is obtained from the continuous time model (23), after inversion of the inputs/outputs, by assuming that finite length-N input sequences are transmitted on finite length-L subchannel impulse responses.In that case, the convolution with the channel impulse responses can be implemented with the matrix multiplication: where the input vector of length 10N R N is: with: (1,i), (1,j),0 (1,i),(N R ,j),0 G X,Y (1,i),(1,j),1 that improves when the design parameter N increases because the effects on the borders of the finite sequence of symbols become negligible.The SDMA precoder for subchannel k is equal to the equalizer for subchannel k in the left-right mirrored system.The computation of the zero-forcing filter requires the inversion of the correlation matrix (G k ) H • G k which is well-conditioned when the number of rows of G k is larger than its number of columns (12N T (N + L − 1) ≥ 10N R N).This is always the case when N T ≥ N R , independently of the channel length L and the number of symbols N. It can also be observed that the zero-forcing filter can be designed at the symbol rate (instead of at twice the symbol rate) when 6N T (N + L − 1) ≥ 10N R N which is always the case when N T ≥ 2N R .

Figure 1
Figure 1 System model.FBMC-OQAM system model for subchannel k (base station antenna n T , receiver n R ).
and imaginary parts, and upsampled by two to obtain the real sequences: s R (n T ,k) [ m] and s I (n T ,k) [ m]: s R (n T ,k) [ m] := s (n T ,k) [ l] ,

Figure 3
Figure 3 Precoder computation accuracy.Mean square error (MSE) as a function of the dimension N of the matrix model used to compute the precoder for a varying number of subchannels.Solid lines: pulse shaping filter of roll-off 1 for M varying; dashed line: pulse shaping filter of roll-off 0.5 for M = 64.Parameters: N T = 4, N R = 3.

Figure 6 Figure 7
Figure 6 Performance comparison with the state-of-the-art SDMA precoder.Bit error rate (BER) of the SDMA FBMC-OQAM system for a varying number of subchannels M. Solid lines: new SDMA precoder; dashed lines: state-of-the-art single-tap SDMA precoder.Parameters: N R = 4, N T = 3.