SER/BER expression for M-QAM OFDM systems with imperfect channel estimation and I/Q imbalance

We investigate the joint effect of channel estimation and frequency flat transmitter and receiver I/Q imbalance on an Orthogonal Frequency Division Multiplex (OFDM) system. We assume independent fading with identically and independently distributed (i.i.d) mirror carrier channel coefficients. Closed form expressions for symbol error rate and bit error rate of M-QAM modulation are derived. We consider join and separate channel estimation, as well as joint and separate equalization of the mirror carriers. Performance is evaluated by treating the I/Q interference as a non-Gaussian random variable. The results show that a Gaussian approximation of I/Q interference is very good, especially for low order modulations.


Introduction
Orthogonal frequency division multiplexing (OFDM) based transceivers have been selected for several wireless standards such as IEEE 802.11,LTE and WiMax, and are the likely selection for 4th generation systems fulfilling the International Mobile Telephony (IMT)-Advanced requirements [1].Low-cost implementation of transceiver devices is challenging in view of impairments associated with imperfect analog components.Imperfections due to power amplifier (PA) non-linearity, frequency offset, phase noise, timing jitter, quantization noise and I/Q imbalance distort the baseband signal [2].The most challenging IMT-Advanced requirements are for local area scenarios [1].In local area, transmission powers are comparatively small, which alleviates problems with PA non-linearity.In OFDM systems, carrier frequency offset, phase noise and timing jitter cause inter-carrier interference (ICI) which is strongest between neighboring subcarriers, and their relative strength is related to the subcarrier bandwidth.Their effect may thus be reduced by selecting a suitable OFDM parametrization.With increasing computational capabilities, the number of bits in analog-to-digital conversion may be increased, reducing quantization noise.Thus the only impairment left without a straight forward engineering cure is I/Q imbalance.
Different approaches for digital compensation of the I/Q imbalance have been proposed [3][4][5].The goal of compensation is to provide an improved image rejection, which by nature depends on the accuracy of the digital estimation and compensation approaches used.
Knowledge about the quantitative relationship between transceiver parameters such as the image rejection ratio (IRR) and system parameters such as the symbol error probability is essential for the design and dimensioning of communications systems.In practice, channel estimation cannot be done perfectly in a fading environment.
Several solutions have been proposed in the literature to estimate the I/Q imbalance parameters, and to compensate in analog and digital domains [4][5][6][7][8][9][10][11].Blind estimation and compensation techniques [5], interference cancelation based [6], pilot aided and preamble based techniques have been discussed to estimate the channel and the I/Q imbalance parameters [8][9][10].Joint estimation of channel and the I/Q parameters can be found in [11].http://jwcn.eurasipjournals.com/content/2012/1/303 The I/Q imbalance compensation for different fading environments and modulation orders, have been investigated in [3].Performance degradation due to I/Q imbalance under perfect channel knowledge is evaluated in [12].It is a challenging task to implement computationally efficient estimation techniques while minimizing the estimation error with the impact of I/Q imbalance.
Krondorf and Fettwis [13,14] numerically evaluate the bit error rate performance of a coded and uncoded OFDM systems with all impairments including carrier frequency offset (CFO), I/Q imbalance, frequency selectivity and channel estimation errors on direct conversion ZF receiver with Least Square (LS) channel estimation.The inter-carrier interference (ICI) due to CFO and I/Q imbalance on effective transmitted symbol is assumed to be Gaussian distributed.
In this article, we concentrate on digital compensation of the I/Q imbalance.The goal is to understand system performance of three different estimation and equalization schemes which are commonly used in communication systems.The baseline solution is separate channel estimation and separate equalization of mirror subcarriers.This is contrasted with joint channel estimation with separate or joint equalization.We present closed form expressions for both the symbol error rate (SER) and the bit error rate (BER) with Gray coded bit mapping in independently Rayleigh fading channels under the joint effects of least squares channel estimation and frequency flat transmitter and receiver I/Q imbalance.Comparing to the Gaussian BER approximation results of [13], we observe that the Gaussian approximation works well, especially for small constellation sizes.
The remainder of this article is organized as follows.Section 'System model' describes the system model under the effect of I/Q imbalance and imperfect channel estimation.In Section 'Signal models with channel estimation errors' , we evaluate the error displacement vector for three different channel estimation and equalization schemes related to joint/separate treatment of the mirror subcarriers.In Section 'M-QAM symbol error rate' , we calculate the SER based on displacement error vector analysis.It is further simplified using the power distribution pattern of the desired and mirror signals.Then, the previous result is extended towards BER in Section 'Bit error rate for M-QAM symbol transmission' with a closed form solution for square M-QAM modulations.Analytical limits for perfect channel knowledge are addressed in Section 'Analytical limits' .Ultimately, simulation results and conclusions are made in Sections 'Simulation results' and 'Conclusion' , respectively.
Notation.E(•), Q(•), erfc(•), and P(E) symbolize the expectation operator, the Gaussian Q-function, complementary error function and the error probability, respectively.The superscripts * stands for complex conjugate and det(•) and | • | refer the matrix determinant and absolute value in that order.

System model
We consider the effect of transmitter and receiver I/Qimbalance on channel estimation in broadband multicarrier communication with direct conversion transceivers.

I/Q imbalance model
In the down-conversion of an ideal direct conversion receiver, the incoming signals in the I (in-phase) and Q (quadrature) branches are perfectly match with gains and orthogonal in phase.However, real receivers introduce gain mismatch and non-orthogonality in phase of I and Q branches due to imperfect electronic implementation.These imperfections at the demodulator can be seen as a receiver I/Q imbalance in baseband domain [3].Similarly, the phase and gain mismatch of the modulator introduce the transmitter I/Q imbalance.
The transmitter I/Q imbalance can be modeled as where s(t) is the output at the modulator after I/Q corruption and x(t) is the ideal baseband equivalent transmitted signal under perfect matching.In a same fashion demodulator output after receiver I/Q imbalance is y(t) = K 1 r(t) + K 2 r * (t) where r(t) is the input signal at the receiver front end which is corrupted by the channel h(t) ∼ CN (0, σ 2 h ) and noise n(t) ∼ CN (0, N 0 ).Both transmitter and receiver I/Q imbalance coefficients )/2 are the complex coefficients which can be determined by the total effective amplitude and phase imbalances of the transceivers, g and φ, respectively.Notice that the above impairment model assumes frequency-flat I/Q imbalance within the signal band which is a rather typical assumption in this context.
The above I/Q imbalance model corresponds to the following transformations of the effective baseband equivalent signals at the receiver output In the equivalent frequency domain signal, there is ICI between mirror carrier pairs which are located symmetrically around the center subcarrier [10] (see Figure 1); where 1 are the I/Q modulated frequency domain channel coefficients and h + and h − are the channel coefficients on mirror carrier pair.The received signal on subcarrier is y + and the corresponding transmitted symbols on mirror carriers are x + and x − .Similarly noise samples can be defined.The mirror-frequency attenuation is described by the image rejection ratio (IRR) of the receiver RF front-end, , which is typically in the order of 25-40 dB [15] in practical receivers.Below, we shall frequently use the inverse IRR, the image leakage ratio, in the analysis.

Channel estimation error model
In general, the estimation error for a channel c is given by e = ĉ − c, where ĉ indicates the channel estimate.Only for minimum mean square error-channel estimation, the error term e and the estimated channel ĉ are uncorrelated; E{e * ĉ} = 0 [16].A versatile channel estimation model in [17] describes the correlated channel estimation error for an arbitrary linear channel estimators.The model decouples the correlation between ĉ and e by introducing a new uncorrelated random variable (RV) z, where a = ρ σ c σc is the biasing factor and σ 2 c and σ 2 c are the variances of the channel and its estimate.The normalized correlation coefficient is ρ = E{ĉ * c}/σ c σc .The uncorrelated estimation error z can be obtained from a complex Gaussian process CN (0, σ 2 z ), where . LS channel estimation is a subset of the above error model and error e in general is the zero mean with a certain variance due to effective noise and it is independent of c over a slow fading channels [16].Using (3), we can decouple the LS estimation error to uncorrelated z and a biasing factor.Then , where σ 2 e is the error variance.Hence, the biasing factor becomes a =

Signal models with channel estimation errors
This article considers different combinations of I/Qimbalance aware channel estimation and equalization.Both estimation and equalization can be performed separately or jointly for the mirror subcarrier pair.Three combinations make sense: (i) Baseline solution with separate channel estimation and separate equalization.(ii) Joint channel estimation but separate equalization.(iii) Joint channel estimation and joint equalization.These can be easily described under the channel estimation error model (3), and a similar error displacement vector model.
Frequency selective channels with independently fading mirror carrier pairs are considered.This is a realistic assumption in broadband communications where the channel coherence bandwidth is typically smaller than the frequency separation between mirror carrier pair.We consider LS estimation throughout.
The pair of I/Q modulated channel coefficients h ++ and h +− at each subcarrier index can be estimated either separately or jointly, using a known sequence of pilots.The I/Q corrupted signal models with channel estimation errors can be constructed based on (2) for the three cases considered.

Separate channel estimation with separate equalization
A known BPSK pilot sequence is inserted in to the OFDM symbol to estimate the I/Q modulated channel the kth subcarrier, assuming the channel is time invariant under several OFDM symbols.The receiver estimates the channel without awareness of the I/Q imbalance, using a LS estimator.Using (2), the estimated channel becomes For known pilots x p+ and x p− on the mirror carrier pair, the estimation error becomes a complex Gaussian RV with zero mean and the variance with same pilot transmission on mirror carrier pair.The noise power spectral density of the pilot transmission is given by N p and σ 2 h is the variance of the physical channel h + (not of the channel to be estimated, h ++ ).
After single tap equalization, the estimate of the data symbol on the subcarrier of interest is The error displacement vector, characterizing the symbol estimate error becomes where the estimation error z ++ is independent from ĥ++ by definition and n + , n − are also independent processes.However, there is a correlation exist between h +− and ĥ++ .Hence, can be modeled as a fraction of two correlated complex Gaussian RVs when conditioned on a given set of transmitted symbols on mirror subcarrier pair.

Joint channel estimation and separate equalization
Here, we consider a receiver where there is a simple pilot structure to mitigate the effect of the I/Q imbalance on the channel estimation, but where I/Q compensation is not considered in equalization.This case was analyzed in [18].We formulate the analysis in terms of estimating the full h until explicit equalization is performed.To remove the effect of I/Q imbalance from channel estimation, pilot transmission from two consecutive OFDM symbols is considered and the pilot symbols transmitted on these subcarriers are where x j , j = 1, 2 are the pilot data symbols in two consecutive OFDM symbols.These are selected so that X p becomes non-singular for each mirror subcarrier pair, i.e., det A simple selection maintaining the orthogonality between subcarriers is based on a Walsh-Hadamard matrix [10]; x j+ = 1, x 1− = −1, and x 2− = 1.The estimated channel coefficients ĥ are functions of these pilot data symbols.The LS estimate of the channel is ĥ = h + X −1 p n.Here the elements of the I/Q corrupted noise vector n are ni = K 1 n i+ +K 2 n * i− .In terms of the channel estimation error, the received signal is Single-tap zero-forcing equalization against a ĥ++ , assuming channel statistics are known to the receiver, gives the error displacement vector where a = . As in the previous case, can be modeled as a ratio of two correlated complex Gaussian RVs.

Joint channel estimation and joint equalization
In addition to the wanted and interfering signal components on subcarrier k, h ++ , and h +− , we now also treat the corresponding quantities on the mirror subcarrier −k.For this, we define These channel coefficients can be estimated as above.All channel coefficients experience channel estimation error with variance The biasing factor for both h ++ and , and for h +− and . The joint signal model for the mirror pair is After I/Q-compensating two-tap zero forcing equalization, the symbol estimates become The estimated data vector x, transmitted symbol vector x, estimated channel matrix Ĥ, uncorrelated estimation error matrix Z and I/Q corrupted noise vector ñ are From ( 11) and ( 12), the estimated symbol on subcarrier can be written as where For the biasing factors, a ≥ b holds for all possible signalto-noise ratios (SNR) values, and ĥ+− ĥ++ ĥ * −− is very small.
Hence, the displacement vector can be approximated by neglecting the latter part of ( 14) to Now, z +− is independent from ĥ++ and then the displacement vector can be modeled as a ratio of two independent complex Gaussian RVs.

Probability density function of the displacement vector
The displacement vectors of first two schemes ( 6) and ( 9) can be modeled as quotients of two correlated complex Gaussian random variables and the last scheme (15) can be treated as a quotient of two independent complex Gaussian random variables.For all cases, let us define the numerator and the denominator as ) and σ 2 D = a 2 σ 2 h ++ , conditioned on a symbol pair x + and x − .The energies of symbols on the subcarrier and its mirror are given by P ± = |x ± | 2 .The numerator and the denominator can be represented in a polar form with corresponding magnitudes by r N and r D , and angles by θ N and θ D .The joint distribution can be found in [19]; where is the inverse covariance matrix of Hence, ϕ 11 /ϕ 22 = σ 2 N /σ 2 D and correlation coefficients for first two schemes discussed above are , respectively.For the last scheme it becomes zero.Let us define the magnitude r = r N /r D and the phase θ = θ D − θ N of the displacement vector.These follow the joint distribution 2 . ( For the derivation of SER, it is convenient to have a PDF of in Cartesian coordinates instead of polar form.The joint PDF f (r, θ) is given by f (x, y) = 1 r f (r, θ) where x = r cos θ and y = r sin θ; When the numerator and the denominator are independent, the joint distribution can be simplified to

M-QAM symbol error rate
In this section, we derive the SER expression in an OFDM system with I/Q imbalance, imperfect channel estimation and independent fading.

Symbol error rate for M-QAM symbol
We focus on the reception of regular M-QAM signals which are most commonly used in multi-carrier systems.
The constellation diagram for 16-QAM with basic constellation boundaries is shown in Figure 2. The amplitude levels of adjacent symbols in both I and Q branches are separated by the distance d.Assuming equi-probable symbol transmission of a rectangular M-QAM consisting of direct product of When the displacement vector moves the received signal to the decision region of another constellation point, a symbol error occurs.To simplify the analysis, we determine two types of error events with corresponding probabilities.In events A 1 , there is an error in the I-or Q-branch, but not in both.For a given constellation point, up to four kinds of A 1 events may occur, due to a displacement in one of the four directions ±I, ±Q.In events A 2 , there is an error in both branches, and similarly up to four kinds of A 2 events may occur for a given constellation point.The constellation points may be split to three classes-corner, edge and middle points.We denote these with subscripts C, E, and M, respectively.The multiplicities of the error events http://jwcn.eurasipjournals.com/content/2012/1/303 of different types are different for the different classes, as shown in Figure 2. The energy of the desired symbol P + = P m,n can be written as where Mi = M i /2, i = 1, 2, m and n are integer indices which represent the desired symbol.Similarly, the energy of the image symbol P − =P p,q can be written by the indices p and q with the same range as m and n, respectively.Given the joint PDF in the Cartesian plane, the error probabilities P( ∈ A 1 ) and P( ∈ A 2 ) can be evaluated for a given pair of desired and mirror symbol characterized by the indices m, n and p, q.
However, for the receiver I/Q imbalance closed form solutions are exist, for generic covariance matrix the integrals above are not solvable in closed forms (when transmitter I/Q imbalance is very small then ρ → 0 and hence the distribution of the displacement vector follows (20)).To progress with the analysis a first order approximation can be applied.Taking the first two terms in the Taylor series expansion around where where and λ m,n,p,q = c 2 /d 2 .In these expressions, The error probability can be evaluated by considering the three classes of constellation points separately which is depicted in Figure 2. The conditional probabilities can be represented using disjoint events A 1 and A 2 .Then http://jwcn.eurasipjournals.com/content/2012/1/303we can derive the error probabilities for the three set of constellation points, corner, edge and middle.
and i m,n,p,q is a function of both desired and image symbols x + = (1 + 2m) + j(1 + 2n) d and x − = (1 + 2p) + j(1 + 2q) d.The multiplicities are The conditioning is removed by averaging over all possible desired and mirror symbol points.First removing the conditioning over the mirror symbol, we define i m,n = M1 −1 p=0 M2 −1 q=0 i m,n,p,q .As the channel coefficients on mirror carriers are independent, the average SER for the given OFDM system equals the average SER per mirror carrier pair.The contribution to the symbol error probabilities due to corner, edge and middle points now become where M is the size of the constellation.The average SER becomes These expressions lead to evaluate the average SER for all three cases, with the corresponding intermediate variables derived in Section 'Signal models with channel estimation errors' .

Closed form expression for SER
Here, we refine the analysis for square QAM constellations with √ M even.The probability density function f (x, y) depends on the transmitted symbol energies on mirror subcarriers.The symbol energy distribution among the four quadrants is identical and symmetric around the diagonal points, then it is sufficient to analyze one quadrant.Thus, (24) can be further simplified to where the pre-factors μ 0 , μ 1 , and μ 2 can be found in Table 1.SNR is denoted by γ = E s /N 0 .The sums i m,n over the mirror symbols can be decomposed in to two by considering the energy of the mirror signal of the first quadrant.
where M = √ M/2.The have proven Proposition 1.The symbol-error probability of an OFDM system with independently fading subcarriers, least squares channel estimation and I/Q imbalance is Example: 4-QAM: From (32), we can derive the SER for 4-QAM easily.There is only one constellation point inside the first quadrant, which is a corner point. where Example: 16-QAM: As shown in Figure 2, in the first quadrant there is one corner point m = 1 with energy P 1 = 2, two edge points m = 2, 3 with energy P 2 = 10, Table 1 μ  No IQ-compensation

Bit error rate for M-QAM symbol transmission
Next, we calculate the bit error rate of an I/Q-distorted system, assuming square QAM symbols, with even √ M, and Gray-coded bits.For Gray-coded M-QAM, the decision boundaries for horizontal signals are independent of the vertical signal levels, and vice versa [20].Thus for perfect channel information it is sufficient to calculate the BER for horizontal √ M-PAM.The probability distribution of the displacement vector in ( 19) is not circularly symmetric over given set of symbol pair x + and x − thus the average error probability along the horizontal and vertical directions are not exactly the same.Therefore, both horizontal and vertical PAM error probabilities need to be evaluated.The individual error probabilities along both directions depend on the marginal densities.Next, the error probability on horizontal vertical PAM can be considered.
Consider a M-QAM symbol x + = x r+ + jx i+ , with distance d between adjacent symbols.Let ¯ = ( ) be the real part of the displacement error vector corresponding to the error probabilities along the real axis.For horizontal PAM, the BER for a given imaginary part of a symbol x i+ and the decision boundary B k can be defined as the probability that ¯ and x r+ fall on the different sides of the decision boundary.Thus the conditional BER for the horizontal part for a particular bit position b k of M-QAM for a given imaginary part of the desired symbol x i+ and the image signal x − is The conditioning on the imaginary part x i+ and the mirror symbol x − is removed by averaging over all possibilities Hence, the total error probability of horizontal PAM becomes Here ∀x ± represents the all QAM constellation points which transmit on either subcarrier or its mirror.Similarly, the bit error rate of vertical PAM can be derived.For square QAM these BERs are the same.The inner sum of the equations above represents the sum of individual bit error probabilities.For a given √ M-PAM signal with √ M even, there are an equal number of points, M = √ M/2, on the positive and negative side of the constellation, and the allocated number of bits per constellation point is L = log 2 √ M. The decision boundaries can be parameterized as B k ±r = ±2 L−k (2r − 1)d, k > 1, with a boundary index r and a bit index k, c.f., Figure 3.
The decision boundaries on right side of B 1 1 have positive boundary index, B k r , and boundaries on the left side have negative indices, B k −r .The number of boundaries per side corresponding to the k th bit position is R = 2 k−2 , k > 1, where r = 1, . . ., R. The sum over the BERs of the √ M bits becomes The coefficients s m ±r , depending on the bit position index m and the boundary index r are Note that for k = 2, the coefficients s m ±r = 1, ∀m.For the above square QAM system with unequal error probabilities on horizontal and vertical PAM, the average BER thus simplifies to [20] The evaluation of the individual probabilities P( ¯ m > κd) is not as simple as in the Gaussian case, because each desired and mirror symbol combination, contributes their own PDF which is centered around the corresponding bit position.Evaluating these, we have Proposition 2. The bit-error probability of an OFDM system with independently fading subcarriers, least squares channel estimation and I/Q imbalance is where

The coefficients s m
±r can be found in (37).The I/Qimbalance affects the terms λ m,n,p,q of Equation (24).Similarly, for the vertical PAM individual error probabilities given by P ¯ n > κ i d = g λ m,n,p,q , κ i which can be derived from its own marginal distribution f (y).
Proof.The marginal distributions of is evaluated first from (19), where The corresponding κ i values for vertical PAM can be obtained by replacing m with n.Then, the average BER can be evaluated from (38)

Analytical limits
From the derived results, known results for error rates in AWGN and Rayleigh channels with and without I/Q imbalance in the absence of channel estimation errors can be derived for square M-QAM constellation.

AWGN channel without I/Q imbalance
As a reference case, in AWGN with perfect receivers 24) can be reduced to where Q(•) is a Gaussian Q-function, and γ is the SNR defined as before.Hence we recover the well known expression [21] P

AWGN channel with receiver I/Q imbalance
In an AWGN channel with receiver I/Q imbalance, the displacement vector can be represented as a complex Gaussian random variable with mean μ = K 2 K 1 x * − and variance σ 2 = (1 + K)N 0 .Hence the PDF of becomes a complex Gaussian distribution.
where μ x = {μ } and μ y = {μ }.The SER is connected to the cumulative distribution functions (CDFs) F X (d) and F Y (d) by where individual CDFs can be obtained by substituting corresponding means μ x and μ y to The argument inside the Q functions can be written This result coincides with [22].

Rayleigh faded AWGN channel with receiver I/Q imbalance
Similarly, for an I/Q corrupted Rayleigh faded system (24) becomes With I/Q imbalance at the receiver, Equation (29) can be reduced to , reproducing the result of [23].

Simulation results
In this section, we compare the above analytical SER/BER expressions to simulation results for an OFDM system in frequency selective fading, with transmitter and receiver I/Q imbalance and channel estimation errors.The OFDM system discussed here is assumed properly synchronize in time and frequency, and the cyclic prefix is sufficient to remove the multipath effect completely.For the simulation, we consider different M-QAM modulations (M = 16, 64, 256).The effect of I/Q imbalance is investigated for different IRR values ranging from 25-55 dB.For the worst case, we select g = 0.9 and φ = 2 o for both transmitter and receiver.The channel coefficients on mirror carrier pairs can be modeled by an independent complex Gaussian process CN (0, 1).The noise is assumed to be white complex Gaussian CN (0, N 0 ) as well.Pilots on mirror carrier pairs are selected from BPSK constellation as in [8,10] for channel estimation.When receiving each information symbol transmitted, a least square channel estimate is used which has been estimated from a pilot symbol transmitted with exactly the same channel.When joint channel estimation is performed, the pilots are chosen from a Walsh-Hadamard code with BPSK symbols, as in (7).
Figure 4 shows the theoretical and simulation results for the three cases of channel estimation and equalization discussed in Section 'Signal models with channel estimation errors' .When there is no IQ-compensation at the receiver (case a), SER performance is poor.The channel estimation error in (4) contains both the I/Q corrupted mirror carrier interference and the channel noise.There is an irreducible error floor visible, due to the mirror carrier interference which dominates at high SNR.When there is IQ-compensation in the channel estimation (case b), the performance is slightly better.However, the effect of I/Q imbalance in the channel noise term, not compensated in equalization, leads to an error floor as in as the previous case.In case c, we compensate for I/Q imbalance both in channel estimation and equalization phases.Hence the performance is better than in the other two cases, and the error floor is removed.The approximation works perfectly and matches the simulation result well.In the two first cases, the closed form analytical result matches simulations perfectly.
Figure 5 shows the behavior of the biasing factors a and b for the above three channel estimation and equalization techniques for a worst I/Q condition (IRR = 25 dB) and a close-to-ideal (IRR = 55 dB).The biasing factor on subcarrier channel estimation can be approximated to a ≈ [ 1+ N p 2σ 2 h ] −1 (for the separate channel estimation and equal- ] −1 ) for the above IRR range.Hence, the effect of IRR is insignificant and the AWGN term dominates the performance.Therefore, at high SNR values a converges to unity Similarly, the biasing factor on mirror carrier channel estimation can be written as For the given channel statistics, the behavior of b for a I/Q corrupted receivers dominates by the SNR.Hence, it converges to unity when SNR increases.But for a ideal receiver with K = 0 it converges to zero irrespective of SNR and the transmitter I/Q imbalance (|G 1 | 2 ≈ 1).
The error rate performance simulations (Figures 6, 7 equalization which is discussed in Section 'Signal models with channel estimation errors' . Figure 6 depicts the SER performance of the ZF equalization on the perfect channel knowledge and the effect of imperfect channel estimation.The effect of channel estimation is simulated with Gaussian approximated mirror carrier interference as well.From Figure 4, it is clear that the biasing factor a does not have any significant impact on system performance.For the relatively small K, a → a ideal =[ 1+ an ideal receiver which again emphasis the above facts.At higher SNR, the estimation error term z ++ [σ 2 z ++ = (1 − a)σ 2 h + ] converges to zero, while the biasing factor goes to unity.Therefore, the imperfect channel estimation follows the equalization with perfect channel knowledge at higher SNR values.Performance of the Gaussian approximated system shows slightly worst error floors with higher order modulations.However, it is reasonable to assume the mirror carrier interference as Gaussian distributed specially at low order modulation schemes.
The effect of transmitter I/Q imbalance is mild compared to the receiver I/Q imbalance and rest of the simulations consider only the receiver I/Q imbalance.There is an irreducible error floor because of the I/Q induced mirror carrier interference which dominates the effective interference at high SNR region.The convergence limit changes according to the IRR value appears at the receiver.Hence, asymptotic SINR approximately equivalents to the IRR and will affects to the SER convergence as depicted in the same plot (projection on to the ideal receiver implies that the Gaussian approximation of I/Q interference).Eventually, the receiver close to ideal receiver perform better.
Then Figure 8 shows the BER performance for 16/64/256-QAM symbols.It validates the expressions derived in (39) with the simulation result for M-QAM symbols with Gray coded bit mapping.Figure 9 interprets the BER performance for 16/64-QAM with IRR = 25, 35.For comparison ideal receiver (without the I/Q imbalance K → 0) also plotted in the same figure.For higher order modulations with higher image rejection follow the ideal receiver for a lower SNR region.But for lower modulation schemes, performance deviates more from the ideal case.But in general, system converges to a particular BER value at high SNR as a result of uncompensated mirror carrier interference.

Conclusion
Closed form expressions for SER/BER performance for the M-QAM modulated OFDM system in frequency selective Rayleigh fading channel under the effect of transmitter and receiver I/Q imbalance and imperfect channel estimation has been derived.Commonly used three channel estimation and equalization schemes are considered and performance is evaluated.If a Gaussian approximation is applied, good results for lower modulations orders are found.For higher order modulation, a Gaussian approximation slightly overestimates error floors.Channel coefficients on mirror subcarriers can be assumed uncorrelated for mirror carrier pairs which are separated more than the coherence bandwidth.In broadband communication this would hold for most of the mirror carrier pairs.The SER/BER expressions obtained here can thus be considered as upper bounds for SER/BER in finite delay spread channels, becoming tight in the limit when the bandwidth becomes large compared to the coherence bandwidth.

Figure 1
Figure 1 Block diagram representation of an OFDM system with transmitter and receiver I/Q imbalances.The subcarrier and mirror carrier channels are h + and h − and the transmitter and receiver I/Q imbalance coefficients are G 1 , G 2 and K 1 , K 2 , respectively.

σ 2 c σ 2 c +σ 2 e= ρ 2 .
For LS estimation, the biasing factor is directly given by the correlation coefficient.If error e and channel c are correlated to each other then the biasing factor becomes a = E{ĉ * c}/ σ 2 c and same channel estimation error model can be still valid.
* − and D = a ĥ++ .The multiplicative factors δ + and δ − , and the biasing factors a and b are different for the different cases described before.Both RVs, N and D have zero mean.Their variances are σ

Figure 2
Figure 2 Distinguished areas for the occurrence of symbol errors in 16-QAM constellation.The distance between constellation points is d.A 1 and A 2 are the corresponding disjoint error probability areas.

Figure 3
Figure 3 Gray coded 8-PAM constellation with decision boundaries B k ±r , where k and r represent the corresponding bit index and boundary index.

, 8 ,Figure 4 Figure 5 Figure 6
Figure 4 SER performance for the three estimation and equalization schemes.(a) Separate channel estimation with separate equalization.(b) Joint channel estimation with separate equalization.(c) Joint channel estimation with joint equalization.Simulation and theoretical results for 16-QAM with gain g = 0.9 and phase φ = 2 o for both transmitter and receiver.

N p σ 2 h] − 1 Figure 7 Figure 8
Figure 7 SER of 16-QAM modulation of ZF equalization with perfect channel estimation and with imperfect channel estimation in Rayleigh fading environment for different receiver IRR values (IRR = 25, 35, 55, ∞ dB).Asymptotic SINR values are marked in red.