Design of Robust Pulses to Insufficient Synchronization for OFDM/OQAM Systems in Doubly Dispersive Channels

This paper presents a pulse shaping method robust to insufficient synchronization in orthogonal frequency division multiplexing with offset quadrature amplitude modulation (OFDM/OQAM) systems over doubly dispersive (DD) channels.The proposed pulse is designed as a linear combination of several well localizedHermite functions.The coefficients optimization problem ismodeled as a nonconvex constrained fractional programming problem based on the signal-to-interference ratio (SIR) maximization criterion. An efficient iterative algorithm is applied to simplify the problem to a series of quadratically constrained quadratic program (QCQP) problemswhich can be solved by semidefinite relaxation (SDR)method. Simulation results show that the proposed pulse is superior to traditional pulses with respect to SIR performance overDD channels in the presence of carrier frequency offset (CFO) and timing offset (TO).


Introduction
Multicarrier (MC) transmission is employed in various communication systems due to its robustness against time dispersion caused by multipath propagation.The systems are generally lattice structures formed by time-frequency shifts of a prototype pulse.Orthogonal frequency division multiplexing (OFDM) is a conventional MC technology for some advantages of efficient implementation through FFT technique and a simple equalization in frequency domain.However, in OFDM systems, the cyclic prefix (CP) used to prevent intersymbol interference (ISI) will cost a loss of spectral efficiency and increase the power consumption.Moreover, the rectangular prototype pulse which exhibits a sinc-shaped frequency response will result in serious intercarrier interference (ICI) in frequency dispersion channels.
In order to counteract the drawbacks of OFDM systems, some pulse-shaping OFDM (PS-OFDM) systems have been proposed [1][2][3][4].These systems are more robust to frequency dispersion caused by Doppler effect and carrier frequency offset (CFO) by employing prototype pulses with lower sidelobe and faster spectral decay.However, the underlying requirement in PS-OFDM systems is that the prototype pulses must be orthogonal or biorthogonal to each other.From frame theory, a well localized orthogonal pulses set cannot be achieved for unit density structures due to Balian-Low theorem.In the literature, many lattice structures employ half lattice density which results in a loss of bandwidth efficiency.
Recently, Filter Bank Multicarrier (FBMC) technology has been intensively studied as an alternative to OFDM in order to improve the robustness against various dispersive channels and to increase the bandwidth efficiency [5].Among all the FBMC systems, the most popular one is the OFDM/OQAM which has already been introduced in the TIA's Digital Radio Technical Standards and considered in WRAN (IEEE802.22)[6].Its principle is to introduce a half symbol duration time offset between the real and imaginary parts of a QAM constellation and transmit them separately on each subcarrier.With this offset QAM approach, the orthogonality constraint only holds in real field.Therefore, a well localized pulse can be introduced in OFDM/OQAM systems with unit lattice density.

Mathematical Problems in Engineering
The time and frequency dispersion of the DD channel is susceptible to cause ISI and ICI in multicarrier transmission.In addition, OFDM/OQAM systems suffer from high sensitivity to some desynchronization like CFO and TO [7,8].CFO is mainly caused by the different local oscillators of the transmitter and receiver.TO may occur when the actual sampling instant at the receiver is different from the optimal one.Normally, CFO and TO estimation is carried out prior to the demodulation.However, the estimation is usually not perfect.Therefore, design of robust pulses to insufficient synchronization over DD channels is an attractive research topic.In [9], the WSSUS pulse shape optimization is translated to a convex constraint quasi-convex maximization problem.In [10], a pulse shape robust to Doppler effect is constructed with optimized combination of Hermite functions.In [11], a novel cost function for optimization of the filter bank prototype pulse is proposed to achieve robust performance over DD channels.In [12], an OFDM/QAM pulse that is optimal with respect to interference due to CFO is introduced.
In this paper, a robust prototype pulse modeled as the linear combination of normalized Hermite functions is proposed for OFDM/OQAM systems.The proposed pulse is designed to be robust to channel dispersion and insufficient synchronization.Firstly, the exact SIR expression of the demodulated symbol versus CFO and TO is derived in the assumption of wide sense stationary uncorrelated scattering (WSSUS) channel.Then, the coefficients of Hermite functions are optimized based on SIR maximization criterion.Finally, the optimization problem is simplified to a series of QCQP problems by using an effective iterative algorithm.By taking into account the semidefinite relaxation, the QCQP problems can be translated to semidefinite programming (SDP) problems which are easy to be solved.Theoretical analysis and simulation results show that the proposed pulse outperforms traditional pulses especially in the case of seriously dispersive channels with large CFO and TO.
The rest of this paper is organized as follows.The system model of OFDM/OQAM is described in the next section.In Section 3, design of the robust pulse is formulated as a SIR maximization problem.In Section 4, the proposed optimization algorithm is presented and performance of the optimal pulse is tested by numerical simulations in Section 5. Finally, some conclusions are drawn in Section 6.

System Model
Instead of using complex symbols in OFDM/QAM systems, the real-valued symbols modulated by offset QAM are transmitted on each subcarrier in OFDM/OQAM scheme.With the synthesis basis functions obtained by the time-frequency translated version of the prototype pulse, the transmitted baseband signal can be written in the following analytic form [6]: where   = (/2)( + ) is the phase term related to the staggering rule of OQAM,  0 is the intercarrier spacing,  0 is the symbol duration time ( 0  0 = 1/2),  is the subcarrier number which is supposed to be even, and  , denotes the real-valued symbol located on the th subcarrier at time instant  0 .The symbols are obtained by introducing half a symbol spacing between the real and imaginary parts of the QAM complex-valued symbols.They are assumed to be independent and identically distributed with zero mean and average energy  2  .() denotes the real-valued even prototype pulse in this paper.
Assuming a distortion-free channel, the following orthogonality condition is required for complete reconstruction of the data symbols: where R{⋅} is the real part operator and (⋅, ⋅) denotes the 2D Kronecker delta function.And   (, ), the ambiguity function of (), is defined by It is observed that   (, ) is real-valued as long as () is an even function [6].Hence, the orthogonality can be ensured by finding a real and even prototype pulse of which the ambiguity function satisfies The baseband doubly dispersive channel can be modeled as a random time-varying system  with impulse response ℎ(, ) [13].Therefore, considering the noiseless case, the received signal can be expressed as where  max and   are the maximum multipath delay spread and Doppler frequency, respectively.We denote  =  max ×   as the channel spread factor.In general, practical wireless channels are usually underspread and satisfy  ≪ 1.The delay-Doppler spread function (, V) is the Fourier transform of ℎ(, ) with respect to .In the WSSUS assumption, the channel is characterized by the second-order statistics of (, V), that is, the scattering function   (, V): where [⋅] denotes the expectation operator.Moreover, it is assumed that [(, V)] = 0.In this paper, we set ∫  ∫ V   (, V)V = 1, which means that the channel has no overall path loss.
Baseband equivalent of an OFDM/OQAM system with  subcarriers in the presence of CFO and TO.
At the receiver, the demodulated signal    0 , 0 for timefrequency position ( 0 ,  0 ) is performed by a projection of the received signal () onto the set of functions {  0 , 0 ()} with the modified inner product.Consider

Pulse Design and SIR Criterion
3.1.Effects of TO, CFO, and DD Channels.Considering the insufficient synchronization, we introduce a timing offset Δ and a carrier frequency offset Δ.Then, the baseband equivalent of an OFDM/OQAM system with  subcarriers in the presence of CFO and TO is depicted in Figure 1.In general, the demodulated data symbol before the real part extraction can be expressed as where Θ Δ,Δ  0 , 0 ;, is given in (A.1) (see Appendix A).If we set  =  0 and  =  0 , then the desired section of the data symbol can be expressed as Equation ( 9) shows that the desired section of the data symbol at the receiver is rotated by a phasor 2[Δ(Δ +  0  0 ) +  0  0 Δ], which is caused by the CFO and TO.To compensate the phase rotation, the demodulated data symbol â 0 , 0 is equalized with a single complex-tap.We can assume the equalizer to be a simple one-tap zero-forcing (ZF) equalizer at each subcarrier and each time instant, which eventually leads the equalized data symbol after real-parttaken operation to where    0  0 =  2[Δ(Δ+ 0  0 )+ 0  0 Δ] is the ZF equalizer coefficient at the ( 0 ,  0 ) subcarrier.It is then easy to separate the equalized data symbol into the useful part and ICI/ISI part; that is, Concerning the useful portion, the desired symbol is attenuated by Λ  ( max ,   , Δ, Δ) which is caused by the DD channel, TO, and CFO.

Mathematical Problems in Engineering
By taking into account the statistically independent property of the data symbols  , and the WSSUS assumption, the total energy of the equalized symbol can be calculated as (see Appendix A) Then, the desired signal energy can be calculated as The interference energy can be calculated as Hence, the SIR expression is defined by Equations ( 13)- (15) illustrate that the SIR depends on the prototype pulse through its ambiguity function when the channel scattering function, CFO, and TO are confirmed.Various prototype pulses, such as square-root raised-cosine (SRRC) [14] and isotropic orthogonal transform algorithm (IOTA) [15] pulses and Hermite pulse (HP) design based on the Hermite functions [16], have been studied.With SRRC pulse, one can get good frequency attenuation property, while IOTA and HP pulses lead to good and identical time and frequency localization.Additionally, the authors in [17] realized a Nyquist pulse in both time and frequency, which is referred to as optimal finite duration pulse (OFDP), by minimizing the out-of-band energy of a time-limited pulse.However, the above pulses are highly sensitive to desynchronization over DD channels.Therefore, in this paper, we search for the optimal prototype pulse based on SIR maximization criterion:

SIR Maximization Criterion Based on the Optimal Pulse.
The optimal pulse is researched as a linear combination of several Hermite functions which are maximally localized in both time and frequency [18].The th order Hermite function ℎ  () is given by where   () is the th order Hermite polynomial which can be computed recursively by  0 () = 1,  1 () = 2, and  +1 () = 2  () − 2 −1 ().We note that the Hermite function corresponding to  = 0 is the Gaussian pulse and provides the best time and frequency localization.It has been proved that ℎ  () is the eigenvector of its Fourier transformation: {ℎ  ()} =   ℎ  ().Moreover,   =  −/2 is the corresponding eigenvalue.Therefore, if we choose  = 2,  = 0, 1, 2, . .., both ℎ  () and its Fourier transformation are real and even functions.Finally, the proposed optimal pulse can be expressed as Equation (18) indicates that the  most concentrated Hermite functions are used.The coefficients   are assumed real.Substituting (18) into (3), we obtain where   (, ) = ∫  ℎ 2 ( + /2)ℎ 2 * ( − /2) −2  and they can be calculated according to [19].Using ( 13), ( 14), (15), and ( 19), after a few steps of straightforward manipulations, we obtain with To precisely evaluate the SIR, it is sufficient to consider only the energy perturbation from neighboring symbols with negligible performance loss by employing the well localized pulses.Therefore, we consider | −  0 | ≤ 4 and | −  0 | ≤ 4 in (22).Finally, the optimization task amounts to finding the coefficients   so that the objective function SIR( max ,   , Δ, Δ) is maximized.

Optimization Strategy and Algorithm
4.1.Optimization Problem Formulation.The optimization of the prototype pulse is achieved by the optimization of the coefficients   that maximize the SIR expression in (20).In order to get simplified, easier-to-manipulate expressions, let a denote the  × 1 column vector containing the coefficients to be optimized: .Then, we introduce the  2 ×1 column vector w = a⊗a, where ⊗ denotes the matrix Kronecker product.As a function of w, the SIR expression in ( 20) can be reorganized in the following matrix form: where S = (  ) and I = (  ) and 0 ≤ ,  ≤  2 − 1 are symmetric and positive definite (SPD) matrices, respectively, defined by   = (, ,   ,   ) and   = (, ,   ,   ), where  =  +  and  =    +   .The ambiguity function   (, ) can also be reorganized in the following matrix form: where A(, ) = (  ) and 0 ≤ ,  ≤  − 1 is defined by   =   (, ).In [16], it has been proposed that the orthogonality constraints (4) were applied at the origin position (, V) = 0 and other  − 1 nearest grid points around it.The operation relies on the good localization property of ().Finally, substituting (24) in ( 4), the coefficients a can be determined by solving the resulting  equations.
Since the ambiguity function is real-valued, the orthogonality conditions in (4) can be transformed to Substituting ( 19) into (25), we obtain where B(, V) = (  ) and 0 ≤ ,  ≤  2 − 1, and it is defined by , where  =  +  and  =    +   .Considering the effect of time and frequency dispersion in DD channels, the constraints on the nulls of the ambiguity function are not necessary [11].Therefore, we simply relax the constraint (26) to w  B(0, 0)w = 1.Besides the orthogonality constraints on w, there are some other inherent constraints due to the operation of matrix Kronecker product.After some computations, reported in Appendix B, we find that two kinds of equality constraints and one inequality constraint are introduced: Cw = 0, w  D  w = 0, and Gw ≥ 0. Eventually, the optimization problem can be formulated as follows:

Iterative Algorithm for SIR Maximization.
As the optimization problem (P) is a fractional programming problem, the objective function (OF) is nonconvex in this case.This kind of problem is NP-hard generally.In order to solve the problem, it is convenient to adopt Dinkelbach's parametric approach [20].By defining Ω as the set of feasible solutions of problem (P), the equivalent parametric problem associated with this problem can be expressed as If there exists  * ≥ 0 for which ( * ) = 0, it has been proved that w * = w( * ) is the optimal solution of problem (P) [20].Then, an effective iterative algorithm which is known as Dinkelbach's method is proposed in Algorithm 1.
By employing the above iterative algorithm, problem (P) can be transformed to a series of QCQP problems as follows: However, considering the variety of , the optimization problem in (29) is a nonconvex quadratic program problem.Fortunately, the semidefinite relaxation (SDR) method can be applied to generate the suboptimal solutions with quite encouraging results [21].It is crucial for the SDR method to observe that w  Hw = Tr(H(ww  )), where Tr(⋅) stands for the trace of a matrix.To get simplified expressions, we introduce a matrix variable M = ww  , which is equivalent to M being a rank-one and positive semidefinite (PSD) matrix.Therefore, we can obtain the following equivalent formulation of (29): Input:  outer (the maximum number of iterations) outer > 0 (the convergence tolerance) (1) Calculating the largest eigenvalue  of I −1 S and the corresponding eigenvector .If  ∈ Ω, terminate.The eigenvector  is the globally optimal solution in this case; otherwise go to Step (2) (2)  ← 0,  0 ← 0 (3) do while  <  outer and ( −1 , w *  ) >  outer (4)  ←  + 1 (5) Solving ( −1 ) to obtain the optimal solution w *  (inner loop) ( 6) end do (8) return Algorithm 1: Iterative algorithm for SIR maximization.
In the equation, M≻0 indicates that M is PSD.It is clear that except for the nonconvex rank-one constraint rank(M) = 1 in 5, the remaining problem is convex.Finally, by dropping the nonconvex constraint, we can obtain a relaxed semidefinite programming (SDP) problem, which can be efficiently solved using the general purpose CVX package [22].

Numerical Results
In this section, the SIR performance of the proposed pulseshaping algorithm is evaluated.The set of prototype pulses such as Hermite, IOTA, SRRC, and OFDP will serve as a basis for our comparison.The CP-OFDM with rectangular pulse is also included for comparison.The experiments are established through a 20 MHz channel for OFDM/OQAM systems where  = 256 subcarriers and 4-QAM are employed.The one-tap ZF equalizer is used at the output of the demodulator to compensate the phase rotation caused by CFO and TO.The WSSUS channel with exponential power delay profile and U-shape Doppler spectrum is adopted.The rectangular time-frequency lattice is set to adapt to the channel property, that is,  0 / 0 =  max /  .Therefore, the transmission pattern of the system is fixed.In addition, Δ and Δ are normalized to  0 and  0 , respectively.The time and frequency responses of the optimal pulse are shown in Figure 2 when the channel spread factor is fixed to 0.01 and both of the normalized CFO and TO are fixed to 0.05.In the figure, the curves of rectangular pulse and OFDP are included for comparison.It is clear that the OFDP yields a better result in frequency domain.On the other hand, the optimal pulse has a better SIR performance in the presence of desynchronization over DD channels which will be confirmed in the following simulations.
It is obvious that the SIR performance will be improved while increasing the number  of Hermite functions.Nevertheless, the improvement decreases with .Meanwhile, increasing  will induce computational complexity and even cause instability in the optimization process.Thus, we choose  = 5 as a compromise.
Figure 3 shows the SIR performance of the considered prototype pulses as a function of the channel spread factor  when the synchronization is perfect (i.e., Δ = Δ = 0).In the simulation, the optimization is redone for each set of .From the figure, it can be seen that the optimal design significantly outperforms other prototype pulses when  ≥ 0.006.However, the SIR is about only 27 dB when  = 0 due to the fact that the constraints on the nulls of the ambiguity function of the proposed pulse are abandoned.The Hermite pulse and IOTA pulse lead to approximately the same performance, so IOTA is not included in the following comparisons.
The SIR curves as a function of the normalized CFO are presented in Figures 4 and 5 for the case where timing is perfect.In the scenario, we set the channel spread factor  = 0.003 in Figure 4 and  = 0.018 in Figure 5.The optimal pulse is uniquely determined by a fixed CFO   which is set to 0.05 in the simulations.From the figures, an improvement of the robustness against the CFO can be addressed by using the optimal pulse.By varying   , we can get a series of optimal pulses and the corresponding maximum SIR.Then, the SIR upper bound can be obtained which shows the improvement margin for a given CFO.It can also be concluded that when the optimal pulse is employed, more SIR performance gain can be obtained for larger CFO and channel spread factor.
Figures 6 and 7 show the SIR results for the case when Δ = 0.The results are given as a function of the normalized TO for a fixed channel spread factor  = 0.003 in Figure 6 and  = 0.018 in Figure 7.In the simulations, we set the optimized point   = 0.03.As is shown in the figures, except for the rectangular pulse in OFDM, the optimal pulse has the highest SIR at the optimized point.CP-OFDM perfectly suppresses the interference caused by TO within cyclic prefix but performs poorly over DD channels.By varying   , we can also get a SIR upper bound.Similar results can be concluded that more SIR performance gain can be obtained for larger TO and channel spread factor.
We also note that, in Figures 4 and 6, where the small value of channel spread factor is used, the optimal pulse performs slightly inferior to other pulses in the small CFO and TO regions.However, the SIR level is still higher than 20 dB.confirms once again the efficiency of our method to acquire a prototype pulse robust to insufficient synchronization over DD channels.

Conclusion
In this paper, we have proposed a procedure for designing prototype pulses for OFDM/OQAM systems that are optimal with respect to SIR performance over doubly dispersive channels with insufficient synchronization.Based on SIR maximization criterion, the pulse design procedure has been formulated as a nonconvex fractional programming problem which was efficiently solved by employing an iterative algorithm.Numerical simulations demonstrated that the optimal pulse provides a significant performance gainover traditional pulses.The advantage is especially prominent in the case of seriously dispersive channels with large carrier frequency and timing offset.

Figure 2 :Figure 3 :
Figure 2: Pulse shapes in time and frequency domain.

Figure 8 Figure 4 :
Figure 8 gives the SIR comparison where the same increment step for CFO and TO is set.we set  = 0.01 and the optimized point   =   = 0.05.The result

Figure 5 :Figure 6 :
Figure 5: SIR as a function of normalized CFO for Δ = 0 over a DD channel with  = 0.018.

Figure 7 :Figure 8 :
Figure 7: SIR as a function of normalized TO for Δ = 0 over a DD channel with  = 0.018.