Fuzzy MIMO detector for MC-CDMA systems with carrier frequency offset over multipath fading channels

A multistage fuzzy minimum output energy (MOE) detector is proposed for multiple input multiple output (MIMO) multi-carrier code-division-multiple access (MC-CDMA) uplink systems with carrier frequency offset (CFO) over multipath fading channels. The first stage of the receiver uses a novel MIMO receiver model with offset symbols to achieve a fuzzy CFO-constrained MOE detector that suppresses multiple access interference (MAI) and minimizes cancellation of the desired signal. To suppress noise and to enhance signal reception after the fuzzy MOE detector, a signal subspace projection and minimum mean square error weight combiner is proposed to enhance signal-to-interference-plus-noise ratio and bit error rate performances. Simulation results show that the proposed MIMO detector outperforms conventional MOE detectors and achieves the ideal receiver performance against MAI and CFO effects.


Introduction
Multi-carrier code-division-multiple access (MC-CDMA) systems [1][2][3][4] are currently under consideration for application in WiMAX and WiFi systems in future broadband wireless communication systems. The system has recently attracted attention because it combines high spectral efficiency with a robust channel equalizer against multipath fading channels. Like a directsequence CDMA (DS-CDMA) system, MC-CDMA inherently overcomes the intersymbol interference caused by multipath propagation. However, its performance is limited by multiple access interference (MAI) and carrier frequency offset (CFO) [5]. The CFO can cause loss of subcarrier orthogonality, which results in intercarrier interference (ICI). To minimize the CFO effect, the receiver in most uplink MC-CDMA systems [6] must first accurately estimate and compensate for CFO to maintain the orthogonality between multi-carriers. Additionally, the rotated constellation of the desired signal degrades the performance of MC-CDMA systems [7,8].
Popular frequency offset estimators are widely reported in the literature [9][10][11][12][13]. For example, the maximum likelihood estimator is the optimal estimator under additive white Gaussian noise. Existing frequency offset estimators typically assume little or no interference. Unfortunately, frequency offset cannot be estimated accurately when MAI is strong. Therefore, an MC-CDMA receiver must efficiently suppress MAI before performing frequency offset compensation. More importantly, MAI suppression must be independent of frequency offset.
The efficient minimum output energy (MOE) multiuser detector developed earlier for use in suppressing MAI in MC-CDMA systems [14] is designed to minimize the entire output energy of the receiver by constraining the spreading code sequence. However, the ICI effect of the CFO in MC-CDMA systems produces a mismatched spreading code signature, which substantially deteriorates MOE detector performance due to the signal cancellation effect [15]. One proposed solution is the use of subspace MOE detectors [16,17], which are designed to project the weight vector of the MOE scheme onto the signal subspace of the received observation data. However, the subspace MOE still involves the problem of the desired signal cancellation to degrade the system performance. Besides, a well-known technique for improving the link quality for wireless communication in the multiple input multiple output (MIMO) scheme is space-time block coding (STBC) [18]. For example, the WiMAX system [19] adopts the spatial multiplexing and diversity schemes of the matrices A and B MIMO schemes [20] used for user uplink communication. Therefore, this study proposes the use of MIMO STBC MC-CDMA system with multiuser uplink communication for the CFO and multipath fading channel environments. Additionally, since the conventional space-time block decoding scheme in [18] cannot be used for the MAI, ICI, and multipath interference scenarios, the proposed MIMO STBC receiver with fuzzy constrained MOE detector suppresses interference and effectively detects desired signals.
Moreover, to acquire spatial diversity, this study develops a robust MIMO fuzzy MOE detector to suppress MAI, to prevent signal cancellation caused by the CFO effect, and to acquire the multipath and spatial diversity gains. The receiver design procedure is the following. First, a novel MIMO received data model with offset symbol is proposed to prevent the CFO from causing symbol detection failure in the STBC MC-CDMA system. Second, the offset MIMO receiver uses a simple spectral peak search scheme for roughly estimating CFO. Although the CFO estimate need not be accurate, it should lead to the correct CFO region. Next, a MIMO fuzzy MOE detector is then proposed for identifying the novel eigen-based signature with MIMO fuzzy CFO constraint used to suppress MAI. Finally, to maximize suppression of noise and reception of signals after the MIMO fuzzy MOE detector, use of the signal subspace projection and minimum mean square error (MMSE) techniques are proposed to enhance signal-to-interference-plus-noise ratio (SINR) and bit error rate (BER) performances. Simulation results confirm that the proposed MIMO detector has sufficient robustness to suppress CFO and MAI effects and that it outperforms conventional detectors. The rest of this article is organized as follows: Section 2 presents the data model and the ideal receiver. Section 3 then presents the robust MIMO multistage fuzzy MOE detector. Section 4 analyzes the computational complexity of the proposed MIMO multistage fuzzy MOE detector. Next, Section 5 summarizes the simulation results. Finally, Section 6 concludes the study.

Notation
The following notation is used. Bold uppercase letters, e. g., A, are used to denote matrices, bold lowercase letters, e.g., a, are used to denote column vectors, and non-bold letters in italics, e.g., a or A, are used to denote scalar values. The symbols A H , A T , A * , and tr(A) denote the transpose-conjugate, the transpose, the complex conjugate, and the trace operations of A, respectively. The diag(a) represents a diagonal matrix whose diagonal entries are the elements of the vector a.

Data model for MIMO MC-CDMA systems
Consider an uplink multiuser MIMO MC-CDMA system over multipath channels with frequency offset. Figure 1 is a block diagram of the proposed MIMO MC-CDMA transceiver. The MIMO structure in the diagram is the simplest case of two transmitting antennas and two receiving antennas, which can be extended to multiple receiver antennas with increased spatial diversity gain. That is, for the scenario of subscriber user with two transmit antennas and base station with N R receiver antennas, the proposed scheme can simultaneously be processed by each receiver. Next, the combiner can then collect the equalized signal from different receivers, which enhances detection performance. In the transmitter, assume a K active user transmission in the uplink MIMO MC-CDMA system, each user is assigned a unique spreading code in the frequency domain such that the complex transmitted signal vectors with N × 1 dimension ðs ðmÞ k ðiÞ; s ðmÞ k ði þ 1ÞÞ; m ¼ 1; 2 of the kth user for the ith and (i + 1)th data symbol after the STBC can be expressed by where the super-scripts (1) and (2) denote the first and second transmitter antennas, respectively, and σ σ k is the transmitted amplitude of two transmitter antennas with total power of σ 2 k . The c k is an N × 1 vector with the spreading code of length N, where c 1 is the spreading code of the desired user 1. The Q H denotes the N × N orthogonal inverse fast Fourier transform matrix. The d k (i) and d k (i + 1) are the ith and (i + 1) th data symbols, respectively. Next, to prevent interblock interference, a cyclic prefix (CP) longer than the multipath channel response is then inserted into each transmitted data block. For fixed wireless communication, the channel impulse response (CIR) L k × 1 vector of user k from the n T th transmit antenna to the n R th receive antenna can be modeled as h h ðn R ;n T Þ k ðL k À 1Þ T , n T , n R = 1 or 2, where L k is the channel length. The assumptions are that CIR h ðn R ;n T Þ k ðiÞ; i ¼ 0; . . . ; L k À 1 is an independent, identically distributed (i.i. d.) complex Gaussian random variable with zero-mean and unit-variance and that h ðn R ;n T Þ k is a quasi-static channel that remains constant during each packet and is independent between packets. Next, let Δε k ¼ Δf k T k denote the unknown normalized CFO of user k between the transmitter and receiver, let Δf k denote the CFO of the kth user, and let T k denote the MC-CDMA symbol duration of the kth user [21]. Figure 1b is a MIMO receiver block diagram showing that, after discarding CP and applying fast Fourier transform (FFT), the ith and (i + 1)th post-FFT received signal blocks with N × 1 dimension from the first and second receiver antennas in the presence of CFO can be expressed by where H , n T , n R = 1 or 2, denotes the N × N frequency domain channel matrix from the n T th transmit antennas to the n R th receive antenna. The N × N CFO effect matrix in (2) is modeled by e jϕ k ðqÞ G k ðΔε k Þ, q = i -1, i, with G k ðΔε k Þ ¼ Qdiag e j2πΔε k N g =N ; . . . ; È e j2πΔε k ðN g þNÀ1Þ=N gQ H and ϕ k ¼ j2πΔε k ðN g þ NÞ=N .
The v ðn R Þ ðiÞ and v ðn R Þ ði þ 1Þ, n R = 1, 2, denote the N × 1 complex white Gaussian noise with zero mean and covariance matrix σ 2 n I NÂN . Without loss of generality, the signal model in (2) then assumes that user 1 is the desired user and that the others are MAI. Thus, the received signals of the ith and (i + 1)th symbol blocks in (2) can be rewritten as respectively, for n R = 1, 2, where C 1 is the N × L 1 frequency shift spreading code matrix of user 1 due to CFO and multipath delay time, i.e., C 1 ¼ G 1 ðΔε 1 ÞC 1 with the N × L 1 composite spreading code matrix C 1 ¼ ½c 1;0 ⋯c 1;l ⋯c 1;L 1 À1

Symbol Block
Spreading IFFT

Insert Cyclic Prefix
Data of User 1 (a)  and c 1;l ¼ diagðc 1 Þq l with the q l being the lth N × 1 column vector of Q matrix. Although the original spreading code c 1,0 is used in the frequency domain in the MC-CDMA system, note that the original spreading code in the CFO and multipath environment is distorted by the CFO and different time delay of multipath effect in the frequency domain. Next, the channel response gain vectors with N × , n T , n R = 1 or 2. The N × 1 MAI vectors are denoted by i ðn R Þ ðiÞ and i ðn R Þ ði þ 1Þ , n R = 1, 2. The conventional STBC system can then be used to construct the orthogonal signal structure and to acquire the diversity gain of the transmitter antenna over the single user in the perfect synchronization scenario. However, in (3) with MAI and CFO effects, the conventional STBC receiver design [18] cannot be used for MAI and CFO scenarios when constructed directly by y n R ðiÞ and y n R ði þ 1Þ, n R ¼ 1; 2. Therefore, the conventional receiver design detects the signal in (3) by first compensating for the frequency offset effect and then using linear weight combiner to suppress the MAI interference. That is, after an ideal frequency offset compensation, the post-FFT received signal vector in (3) can be rewritten as for n R = 1, 2. The post-FFT received signal in (4) for the ith and (i + 1)th symbols can also be cascaded asỹðiÞ with 4 N × 1 vector size: whereh 1 andh 2 are the 4 N × 1 composite channel vectors, i.e., iðiÞ andṽðiÞ are the 4 N × 1 composite MAI and noise vectors, respectively. To recover d 1 (i) and d 1 (i + 1) from yðiÞ, a linear receiver can be used to calculate where w 1 and w 2 are the 4 N × 1 weight vectors. The ideal maximum SINR (MSINR) weight vector can be chosen by [22], i.e., where the 4 N × 4 N covariance matrices are Thus, the desired signals d 1 (i) and d 1 (i + 1) can be detected by x 1 (i) and x 1 (i + 1), respectively. It is noted that the ideal MSINR receiver in (6) is implemented by artificially knowing the desired signal components d 1 (i) and d 1 (i = 1) in the data, and using the true composite channel vectorh 1 andh 2 to obtain (7). That is, the ideal MSINR receiver is designed to perform the optimum bound of the proposed MIMO MC-CDMA system, which cannot be realized since the desired signal and the true composite channel are the ideal known assumptions. Therefore, in this article, we will propose a robust fuzzy MIMO MC-CDMA receiver with MAI suppression over CFO and multipath fading channel environment.
Robust MIMO MC-CDMA receiver Figure 1b shows a block diagram of the multistage fuzzy MIMO MC-CDMA receiver. In the first stage of the receiver, a novel received data structure is proposed to overcome the failure of the STBC MC-CDMA system due to CFO. Next, a robust fuzzy CFO-constrained MIMO MOE method with coarse CFO estimation and signal subspace projection is designed to suppress MAI. Finally, the post-MOE output data are constructively combined by MMSE weight vector to extract the symbols transmitted by the desired user.

Novel MIMO receiver model with offset symbols
Again, consider the scenario of the uplink MIMO MC-CDMA receiver data model in (3) in a multipath fading channel with frequency offset. The complex signals (3) to provide an alternative receiver data model: and The reconstructed 4 N × 1 data vector _ yðiÞ over the ith and (i + 1)th symbols can then be cascaded by is the effective MIMO CFO shift spreading code matrix with 4 N × 4 L 1 dimension and the 4 L 1 × 2 MIMO channels are The _ iðiÞ and _ nðiÞ are the MAI and the noise vectors with 4 N × 1 dimension, respectively. In (10), the 4 N × 4 N data correlation matrix R _ y _ y can be derived by where N MAI-plus-noise correlation matrix. Notably, the two reconstructed symbols (b 1 (i), b 1 ðiÞ) in (11) involve the orthogonal cross-correlation and the same covariance property, i.e., After substituting (12) into (11), the desired signal correlation matrix R s with 4 N × 4 N dimension can be rewritten as Since H 1 and H 2 are CFO-independent, R s can only be affected by the frequency shift matrixC 1 ðΔεÞ. Therefore, this feature can be used for coarse estimation of CFO as described below.

Coarse CFO estimator
Section 2 showed why the optimum receiver must accurately estimate CFO and compensate for the CFO effect to prevent MAI. Here, however, the CFO estimate need not be accurate as required in the proposed fuzzy MOE receiver, but it should indicate the correct CFO region. A fuzzy frequency offset-constrained MIMO MOE receiver is then proposed for effective suppression of strong MAI. Restated, the frequency offset in the conventional system in (4) requires precise estimation and compensation by the receiver so that the optimum receiver can suppress the MAI. The proposed receiver, however, uses the coarse CFO searcher with signal subspace constraint in (15) to estimate the correct CFO region with residual CFO in the same manner as a global searcher for CFO parameters. After the coarse CFO estimation, the fine residual CFO parameter can be considered for use in the MOE detector of Section 3.3 to span the effective signal subspace, to acquire the desired signal reception, and to suppress MAI interference. The results of performance tests are discussed in Section 5. The proposed coarse CFO estimation procedure is performed in the following steps. First, since spreading code C 1 is given, the CFO shift spreading code matrixC 1 ðΔεÞ with 4 N × 4 L 1 dimension can be constructed by the different frequency offset, i.e., where G 1 (Δε) is an N × N matrix, C 1 is an N × L 1 matrix, and Δε is uniformly increased to search for the coarse CFO. Next, based on the reconstructed data model in (13) and the CFO shift spreading code matrix in (14), a simple CFO estimator is determined as the solution to the following spectral search problem: where the 4 N × 2 K matrix U s contains 2 K dominant eigenvectors of R _ y _ y in (11), which span the signal subspace. Note that the signal subspace U s , which is obtained from the signal covariance matrix R s in (11) and (13), can be used to provide the desired signal subspace and to limit the interference subspace to converge the desired coarse CFO estimate for the spectral search in (15). That is, in (15), for the estimated frequency offset Δε 1 approaching the CFO of desired user,C 1 ðΔεÞ produces a spectral lobe at Δε 1 and suppresses MAI outside Δε 1 . Therefore, if the spectral peak error in the spectrum satisfies the condition S i ðΔε 1 Þ À S iÀ1 ðΔε 1 Þ j j = S i ðΔε 1 Þ < 0:1, the Δε 1 is varied with the ith searching iteration by using Δε 1 as the coarsely estimated CFO of the desired user.

Fuzzy CFO-constrained MIMO MOE detector
Based on the previous CFO estimator, the CFO can roughly be estimated by the simple spectral searcher in (15). Next, after the reconstructed data is subjected to the CFO effectC 1 ðΔεÞ as described in (10), a fuzzy CFO-constrained algorithm with coarse CFO initialization is proposed to suppress MAI and to collect multipath energy for the linear MOE detector. The optimization problem for the MIMO MOE can be expressed by where i i is the ith column of the identity matrix I and D is the constraint size. Matrix E 1 , which is the MIMO fuzzy CFO-constrained signature matrix with coarse CFO Δε 1 initiation, is then chosen from the dominant eigenvectors of the distorted signature matrix with a dense set of CFOε n ¼ Δε 1 þ ε n , n = 1, 2,. . ., N t , i.e., where {λ i } and {e i } are eigenvalues and 4 N × 1 eigenvectors of the 4 N × 4N R ε matrix. The shift matrixC 1 ðε n Þ ¼ diag G 1 ðε n ÞC 1 ; e jφ 1 G 1 ðε n ÞC 1 ; G 1 ðε n ÞC 1 ; e jφ 1 G 1 ðε n ÞC 1 n o involves the phase shiftφ 1 ¼ j2πε n ðN g þ NÞ=N . Since the D is the dominant mode with large eigenvalues in (17) and N t is the number of frequencies in the dense set of frequenciesε n , the N t shift matricesC 1 ðε n Þ span an effective rank of D for the fuzzy CFO-constrained matrix E 1 with 4 N × D dimension. The 4 N × 1 weight vector w i of the MIMO fuzzy CFO-constrained MOE detector in (16) can be determined by where W is the 4 N × D weight matrix of fuzzy MOE detector. Notably, the MIMO fuzzy CFO constraint enables the proposed MIMO MOE detector to avoid the problem of sensitivity to signature mismatch arising in the conventional MOE [7], which the signature mismatch effect can significantly degrade system performance. Moreover, in (17), the dominant eigenvectors can span the desired signal subspace containing the mismatch signatures induced by CFO effect. Therefore, the new constrained signature E 1 in (16) can consist of the dominant eigenvectors to collect the desired signal energy and to suppress other user interference.

Signal subspace projection of MOE weight
Fuzzy MOE processing gives the proposed receiver sufficient robustness to suppress strong MAI. However, the fuzzy CFO constraint increases noise in the post-MOE received data. Therefore, the constrained MOE weight matrix in (18) is projected onto the signal subspace, which can minimize the effects of noise and improve system performance. The 4 N × D projected weight matrix W s can be expressed by where U s is the signal subspace of R _ y _ y as in (15).

MMSE weight combiner
After the above multistage processing procedures eliminate MAI, only the desired signal is retained in the projected MOE detector output data z 1,l , l = 1, 2,. . .,D.
These output data z 1;l collected according to MMSE criteria are therefore used to extract the desired signal vector x(i) with 2 × 1 dimension. Let V be the D × 2 weight matrix for the combining procedure, i.e., where the projected MIMO MOE detector output data vector z(i) is a D × 1 vector The MMSE weight matrix can be determined by the following criteria where b 1 (i) and b 1 ðiÞ in (9) are the training symbol sequences. This leads to the well-known Wiener solution

Complexity analysis
This section analyzes and compares the computational complexity of the proposed fuzzy MOE method with conventional receivers such as MOE and subspace MOE detectors. For simplicity, only the number of complex multiplications is considered in the comparisons. Specifically, the proposed multistage detector scheme first calculates the computational complexity of coarse CFO estimator, fuzzy MOE detector, signal subspace projection, and MMSE weight combiner step-by-step. Summing the computational complexity of the multistage schemes then obtains the total computational complexity. Moreover, for the coarse CFO estimator, the number of complex multiplications needed to calculate the CFO shift spreading code matrixC 1 ðΔεÞ in (14) is analyzed as follows The number of complex multiplications of the spectral search algorithm in (15) is then computed as Second, for the fuzzy MOE detector, the number of complex multiplications of the weight matrix W of the MIMO fuzzy MOE detector in (18) is computed as where R _ y _ y and E 1 are the 4 N × 4 N and 4 N × D matrices, respectively. Next, the number of complex multiplications for signal subspace projection in (19) is computed as Finally, for MMSE weight combiner, the number of complex multiplications of the linear combiner in (20) and (21) is Then, the number of complex multiplications for the Wiener solution in (23) is The total number of complex multiplications performed by the proposed multistage MOE detector in the above four stages can be summarized as According to (30), the computational complexity of the proposed multistage detector depends on the spreading code length N, the number of multipath channels L 1 , the number of active users K, and the size D of the dominant eigenvalues of the fuzzy constrained matrix E 1 . The equation also shows that the computational load is approximately O(N 3 ). Note that the computational load of the proposed multistage detector is determined mainly by the weight matrix inversion calculated by the fuzzy MOE detector in (18).
Various existing methods are then used to evaluate the computational complexity of the conventional receivers, e.g., MOE and subspace MOE detectors are considered for comparison. The conventional MOE detector with R À1 _ y _ y as reported in [14] is considered first. The total number of multiplication operations is about O(N 3 ) due to the major computations needed for inversion of R_ y _ y . Next, the algorithm for the conventional subspace MOE detector in [17] resembles that for the projected weight matrix in (19). The major computations still require calculation of MOE weight by R À1 _ y _ y . Therefore, the number of complex multiplication is about O(N 3 ). In summary, the complexity analysis shows that the total number of complex multiplication operations required by the proposed multistage detector, the conventional MOE, and subspace MOE detectors is about O(N 3 ). That is, the MOE calculation results in a similar computational complexity in the three MOE calculation methods.
Although the proposed multistage method requires additional computational operations, mainly because of the coarse CFO estimator and MMSE combiner, the computational load is still smaller than that of the MOE detector O(N 3 ). As shown above, the proposed multistage detector has higher computational complexity compared to the conventional MOE and subspace MOE detectors. However, the simulation results confirm that proposed method outperforms the conventional MOE detector and confirms its robust performance over serious CFO environment.

Computer simulations
Simulations are performed to confirm the performance of the proposed MIMO MC-CDMA system over uplink multipath channels. For all users, the QPSK-modulated signals are spread by orthogonal gold codes of length 32 (N = 32). The channel profile assumes L k = 4 independent equal power Rayleigh fading paths with the time delays chosen from {0,3T s }, which is smaller than the guard interval T G =8T s where T s is the sampling period. In this simulation, the selected D is D = 2 L 1 , i.e., two degrees of freedom are used to retain the desired user signal and to suppress MAI. The quasi-static multipath fading channel in the simulations is assumed to be constant during each packet and independent between packets. The fading gains are the i.i.d. complex Gaussian random variables with zero mean and unit variance. The output SINR, which is used as a performance index, is defined as the ratio of the signal power in x(i) to the MAI-puls-noise power in x(i) at the MMSE combiner output. The input SNR is also defined as SNR ¼ 10 log 10 ðσ 2 1 =σ 2 n Þ , and the near-far ratio (NFR) is defined as NFR ¼ 10 log 10 ðσ 2 k =σ 2 1 Þ, k 6 ¼ 1. The analysis is also simplified by assuming all MAI powers are identical. Executing 10,000 Monte-Carlo trials with different fading gains and data/noise sequences in each trial obtains one output SINR and BER value. The normalized CFOs of all users are also assumed to be in the range À0:5 < Δε k < 0:5. The performance comparison included the results obtained by conventional MOE, subspace MOE, and optimal MSINR detectors. The optimal MSINR receiver is implemented by using perfect CFO compensation and by artificially removing the desired user signal component in the data, and then using the ideal composite channel vectorsh 1 andh 2 to obtain (7). Finally, the following "standard" parameters are used throughout the section unless otherwise mentioned: SNR = 10 dB, NFR = 0 dB, K = 10, Δε 1 = 0.25.
The first set of simulations evaluated the output SINR performance as a function of input SNR in the proposed multistage fuzzy MIMO MC-CDMA system. The simulation results in Figure 2 indicate that the performance of the proposed MIMO receiver with three processing stages and MMSE combiner approaches that of the ideal MSINR receiver. The proposed MOE detector outperforms conventional MOE and subspace MOE detectors, which are subject to the effects of CFO and MAI. Moreover, in a high-input SNR region, the first stage of the receiver, which has a coarse CFO estimator and MMSE combiner, also has a lower output SINR compared to the conventional MOE. Since the MAI power increases with input SNR in the NFR = 0 dB scenario, it performs poorly in CFO estimation. Thus, the MMSE combiner cannot effectively suppress MAI and induce output SINR degradation. The second stage with the coarse CFO estimator, fuzzy CFO-constrained MOE detector, and MMSE combiner performs robustly with only about 5 dB degradation due to the fuzzy CFO with large noise sequence. By using the full-stage processor and MMSE combiner described in Section 3, the performance of the third stage of the receiver approaches that of the ideal receiver as shown in Figure 2 since problematic noise is limited by the signal subspace projection. These results confirm that the proposed multistage MIMO MOE receiver indeed successfully suppresses MAI and effectively overcomes the CFO effect.
The second set of simulations compared the proposed MOE, conventional MOE, subspace MOE, and ideal MSINR receivers in terms of output SINR performance with increased numbers of users K. Figure 3 shows that, again, the MIMO multistage MOE detector successively approaches the ideal MIMO MSINR receivers and outperforms the conventional MOE and subspace MOE detectors. Moreover, the output SINRs of the conventional MOE receiver and first stage MOE receiver increase as the number of users increases. Because the power of MAI with different CFOs increases as the number of users increases, the MOE effectively suppress the MAI of overall CFO region. However, the output SINR performances of the conventional MOE and first stage MOE receivers are severely degraded in comparison with the proposed multistage and optimal receivers. Moreover, as the number of users increases, the proposed full-stage detector does not severely degrade the output SINR performance due to lack of degree of freedom effect, which is overcome by the larger code length. The third set of simulations evaluated the tolerance of the proposed receiver to different frequency offsets. Figure 4 plots the resulting output SINR curves as a function of normalized CFO (Δε 1 ). Within the entire range À0:5 < Δε 1 < 0:5, the output SINR of the proposed receiver is fairly constant. However, the conventional MOE and subspace MOE detectors exhibit severely degraded output SINR because of their extreme sensitivity to inaccuracies in frequency synchronization, which cause the desired signal cancellation. Therefore, Figure 4 indicates that the proposed multistage MOE detector with coarse CFO estimate, fuzzy CFO constraint, signal subspace projection, and MMSE combiner achieves the ideal performance and performs robustly against CFO by effectively collecting the desired signal and by canceling MAI interference. Notably, the performance at CFO = 0 case is better than other CFO cases. It is due to  no ICI interferences, which will be easily to suppress the MAI and get more degrees-of-freedom to enhance the output SINR performance.
The fourth set of simulations tested the robustness of the proposed MIMO MOE receiver against channel estimation error caused by time-varying fast fading channels. In this case, ĥ ðn R ;n T Þ 1 ¼ h ðn R ;n T Þ 1 þ σ Δ Δh is varied in (10) and (24), where Δh is a random vector with the entries being i.i.d. complex Gaussian random variables. Note that the entries of h 1 are i.i.d. complex Gaussian random variables with variance σ 2 h1 . Figure 5 shows the output SINR versus σ Δ =σ h1 . The simulation results confirm that the proposed MIMO MOE receiver performs reliably with relative error limited to 20% in channel estimation. Figure 5 also confirms that the proposed receiver outperforms other receivers in the channel estimation error scenario.
The fifth set of simulations investigated the near-far resistance of the proposed MIMO MOE detector under different NFR values. The output SINR results  in Figure 6 show that the proposed robust MIMO MOE receivers provide excellent near-far resistance by multistage MAI cancellation. However, the subspace MOE detector outperforms the coarse method when the input SNR exceeds 14 dB because the strong MAI is suppressed by the signal subspace projection technique.
The sixth set of simulations investigated output SINR performance under varying code lengths. Figure 7 shows that, again, the proposed MIMO MOE receivers successively approach the ideal MIMO MSINR receiver performance for different code lengths. Moreover, the output SINR performance of the detector substantially improves with code length, which confirms that the proposed receiver is sufficiently robust to overcome CFO and MAI effects.
The final set of simulations investigated BER performance under varying input SNR. Figure 8 shows that the performance of the proposed MIMO MOE receiver again approaches that of the optimal MSINR receiver. In terms of BER, the proposed receiver performs much better compared to conventional MOE detectors, which confirms its superior robustness in suppressing CFO and MAI effects.

Conclusions
The novel robust MIMO MC-CDMA transceiver proposed in this study is suitable for uplink systems because of its high-quality performance over multipath fading channels and CFO effects. In the multiuser receiver, a MIMO fuzzy CFO-constrained MOE detector incorporating the rough CFO estimator suppresses MAI and overcomes the degradation problem of the conventional MOE detector with CFO effect. Next, by using signal subspace projection, the proposed detector further reduces noise with its fuzzy MOE filter. Finally, an MMSE weight combiner is proposed to capture the signal multipath diversity gain. Simulation results confirm that the proposed MIMO receiver performs comparably to the ideal receiver, outperforms conventional MOE detectors, and performs reliably over strong MAI and large frequency offsets.