LIU Yongjin,CHEN Xihong,and ZHAO Yu
Air and Missile Defense College,Air Force Engineering University,Xi’an 710051,China
Abstract:This paper investigates the problem of synchronization for offset quadrature amplitude modulation based orthogonal frequency division multiplexing(OFDM/OQAM)systems based on the genetic algorithm.In order to increase the spectrum efficiency,an improved preamble structure without guard symbols is derived at first.On this basis,instead of deriving the log likelihood function of power spectral density,joint estimation of the symbol timing offset and carrier frequency offset based on the preamble proposed is formulated into a bivariate optimization problem.After that,an improved genetic algorithm is used to find its global optimum solution.Conclusions can be drawn from simulation results that the proposed method has advantages in the joint estimation of synchronization.
Key words:offset quadrature amplitude modulation based orthogonal frequency division multiplexing(OFDM/OQAM),synchronization,joint estimation,genetic algorithm.
Due to the advantages in resistance to inter-symbol interference and inter-carrier interference,offset quadrature amplitude modulation based orthogonal frequency division multiplexing(OFDM/OQAM)has become the candidate of 5th-generation(5G)modulation[1].However,because there are no cyclic prefix inserted in OFDM/OQAM systems,symbol time offset(STO)will lead to malposition during fast Fourier transformation(FFT)and cause block interference,carrier frequency offset(CFO)will destroy the orthogonality between carriers and cause difficulty in correct demodulation of the systems[2].Furthermore,channel estimation and even the performance of the whole system is related to STO and CFO.Therefore,joint STO and CFO estimation for the OFDM/OQAM systems should be valued.
In recent years,the problem of synchronization for OFDM/OQAM systems has attracted more and more attention.All these methods can be divided into two classes:data-aided methods and blind methods.Regarding the blind methods,Fusco proved the conjugate-symmetry property of received signals in the frequency selective channel,which was then used to estimate CFO in[3].The weakness of this method is the low convergence speed.Similarly,the approximate conjugate-symmetry property was also utilized by Mattera and Tanda to estimate the synchronization of the systems blindly in[4].Moreover,signal characteristics were taken full advantage of in[5]to blindly estimate CFO of the doubly-selective fading channels.However,methods aforementioned can be further improved with respect to convergence speed and computational complexity.
Regarding the data-aided methods,Singh and Vasudevan adopted cross correlation and preamble respectively to estimate the coarse synchronization and the fine CFO in[6].In[7],the fine CFO estimation in the Rayleigh fading channels was also achieved based on the interference approximate method.Synchronization was estimated based on the zero autocorrelation code in[8],whose preamble included two symbol intervals.
The synchronization results of the aforementioned dataaided methods are closely related to the preamble structures.Therefore,how to improve the preamble structure is very important in the data-aided synchronization algorithms.In[9],guard symbols were placed at two sides,instead of the same side of reference preamble symbols to modify the preamble structure.However,the spectrum of the system would be reduced due to the guard symbols.For the purpose of overcoming this defect,pilot symbols were coded to carry information that need to be transmitted in[10]to save the system spectrum.However,there are two symbol intervals in the preamble which can be further improved.In[11],a preamble with less spectrum consumption was proposed.In order to increase the transmission ef-ficiency of the system,an efficient preamble structure was designed in this paper according to the thought proposed in[11].
Furthermore,data-aided synchronization algorithms mentioned above need at least two steps to estimate STO and CFO.There is so much preprocessing before the final estimation.It causes complex calculation.The genetic algorithm(GA)is able to solve this difficulty and has the ability to deal with both continuous and discontinuous functions[12].
There have been many researches on GA[13–19]recently.The main features of GA are the fast convergence speed and the optimal answer[20,21].As a computerized search,GA was used by Kakandikar and Nandedkar in[22]to predict and optimize the thinning in automotive sealing cover.In[23],Podlena and Hendtlass proposed an improved GA based on the history network to speed up the convergence.The weakness of this method lies in the existence of random operators.Parents are chosen randomly in the selection operator.Chromosomes change randomly in the mutation operator.The directed random search(DRS)which is used to narrow the search space also depends on the random change of an individual.In addition,updating the exemplar stored in the history network is by merging the old exemplar with the current input randomly.GA will converge slowly with all these random operators.In[24],Seyyed and Samane adopted the pruning operator to avoid the random assignment of the value.In this paper,the GA with the pruning operator is adopted to speed up the convergence and search for the optimal answer.
In this paper,joint estimation of STO and CFO based on the preamble is formulated into an optimization problem,whose objective is to minimize the mean square error(MSE)caused by CFO and STO.After that,the pruning operator is used to improve the GA and find the global solution of the optimization problem.
Contributions of this paper are two-fold:
(i)Preamble structure design
A novel structure of the preamble is designed to save the system spectrum.Zero symbols are generally inserted between preambles and data symbols.In this study,preambles without zero symbols are considered to increase the transmission efficiency.
(ii)Synchronization problem formulation
Generally,power spectral density(PSD)of the system is derived before the log likelihood function is adopted to estimate STO and CFO.Instead,in this paper,joint STO and CFO estimation is straightly formulated as bivariate optimization to reduce calculations.
This paper is organized as follows.In Section 2,the OFDM/OQAM system model is described and the preamble structure is designed.In Section 3,joint STO and CFO estimation is formulated as bivariate optimization and the GA with the pruning operator is used to solve this optimization problem.In Section 4,simulation results are presented.Finally,conclusions are drawn in Section 5.
Fig.1 shows the model of the OFDM/QAM system[25].C2R0represents converting complex signals into real signals.R2C0represents converting real signals into complex signals.F0(z)is the modulation filter function.Q0(z)is the demodulation filter function.H(z)represents the channel impulse response.η[k]represents the noise.The baseband equivalent of the discrete-time OFDM/OQAM transmitting signal can be written as
Fig.1O FDM/OQAM system model
whereM,n,m∈{0,1,...,M?1},cm,nanddm,nrepresent the number of subcarriers,the time index,the frequency index,the input symbols and the real valued OQAM symbols,respectively.ym,n=θm,ndm,n,θm,n=jm+nand filterfm[k]satisfies
whereG[k]represents the prototype filter whose length isLp.
The perfect reconstruction condition of the signalχm,n[k]satisfies
where?{·}is the operation of taking the real part,(·)?represents the complex conjugate,δm,m?denotes the Dirac’s delta function,and
The preamble structure proposed in[26]is shown in Fig.2.It can be seen from Fig.2 that there arewcolumns of zero symbols.
Fig.2Preamble structure proposed in[6]
According to[27],the continuous-time signals(k)can be described as
The ambiguity function ofg(m)is defined as
whereνdenotes the subcarrier spacing andτdenotes the symbol duration.
The relationship of the discrete-time signals[k]and the continuous-time signals(k)issR(kTs)Δ=sR[k],whereTsis the sampling period.
DefineCp,q=jp+q+pq+2pn0Ag(?qτ,pν),and Ω?Δm0,Δn0=ΩΔm0,Δn0?(0,0),where ΩΔm0,Δn0denotes the neighborhood of position(m0,n0).Cp,qgradually approaches to zero along with the increasing of|p|and|q|.For example,for the prototype filter of the isotropic orthogonal transform algorithm(IOTA),when(p,q)/∈Ω1,1,we have[29]
Conclusion can be drawn from(9)that one column of guard symbols is adequate to prevent the pilot symbol from the interference.If there is only one column of guard symbols in the preamble,its channel estimation(CE)performance will be worse than that of the interference approximation method(IAM)[30].As shown in Fig.3,two columns of zero symbols are inserted in IAM.
Fig.3Structure of IAM
According to[11],a higher power of the pseudo pilot means a better performance of CE.The power of pseudo piloty?for the IAM is
Obviously,there are two columns of zero symbols,which leads to the low spectrum efficiency.In order to increase the pseudo pilot power and spectrum efficiency,a preamble without zero symbols is proposed in this paper.The proposed structure of the preamble is shown in Fig.4.Three-tap interferences are taken into account in this paper due to sufficient accuracy.
Fig.4The proposed structure for OQAM system
The power of pseudo piloty?for the proposed structure in Fig.4 can be written as
Conclusion can be drawn from(11)that the structure proposed in Fig.4 can increase the power of the pseudo pilot.It means that the structure proposed in Fig.4 has the ability to decrease the interference caused by the adjacent unknown data symbol.
Considering the additive white Gaussian noise(AWGN)channel with a zero-meanη[k]which is statistically independent ofs[k],when the transmitted signals pass through this channel,the received signalr[k][6]is written as
whereτ0is STO corresponding to one half of a symbol duration andμ0is normalized CFO corresponding to the sub-channel frequency spacing,μ0∈[?Δ,Δ].Assume thatΔis known to the receiver.
According to the previous researches,the second-order effects of STO and CFO are negligible,which therefore can be neglected in the OFDM/OQAM system model[31].
The demodulation signal at preamble position(m0,n0)is
Equation(14)implies that the demodulation of the preamble signal is composed of two parts:the first term presents the interference caused by the surrounding symbols,the second term is the transmitted signalym0,0plus attenuation.It means that CFO and STO will introduce interference into the demodulation process.In this section,the MSE is derived to describe the influence of the interference caused by CFO and STO.
Define
Rewrite(14)in vector notation as
whereζis the noise vector,and
Considering zero-mean AWGNη[k]with varianceδ2,the covariance matrix ofζis
where
According to[32],the MSE can be written as
Therefore,the joint CFO and STO are formulated into bivariate optimization as follows:
GA is a powerful tool to deal with complicated issues and optimization problems.There are many advantages of GA over other techniques,such as the ability to obtain the global optimal answer rather than local optima.The global optimal answers are all contained in the search space which may be expanded by gorges in GA.Furthermore,the space that does not contain answers may cause redundant search,which will slow down the convergence speed of GA.
In[23],Podlena and Hendtlass utilized the history network to store the record of spaces that already have been searched.In this way,the redundant repetition can be avoided in the same generation.For the next generation,the history network will be inquired before iteration about whether the space has been searched.This function is called Baldwin effect.It can be seen that Baldwin effect is able to avoid the redundant search in invalid space.However,Podlena and Hendtlass have not attached importance to the random operators in GA,which probably causes early termination of iteration and even the wrong solution.In[24],Seyyed and Samane applied the pruning operator to avoid the random assignment of the value.There is much in their method that can be used.In order to speed up the convergence,the pruning operator can be adopted to improve the method proposed by Podlena and Hendtlass.In this paper,the pruning operator and the history network are combined with each other to improve the GA.
The time order of the pruned initialization is shown in
whereNcrepresents the number of chromosomes of each generation,Ngrepresents the number of genomes in each chromosome,andNpis determined by the pruning function.It is shown in
whereηis the reduced ratio of the search space.
The time order of the pruned search loop is shown in(22),whereNsrepresents the solution space.
In this section,the CE performance based on the proposed preamble structure is simulated.Truncation of an IOTA filter with length 4T0[33]is chosen as the prototype filter in this simulation.The fundamental parameters are listed in Table 1 where ITU is the International Telecommunications Union.
Table 1Fundamental parameters of CE
Fig.5 illustrates the CE performances of two methods for different kinds of constellation mapping.For convergence,the proposed structure in this paper is named as A and the IAM is named as B.For 16-QAM modulation,a higher signal to noise ratio(SNR)means a better CE performance and B shows worse CE performance than A when SNR<15 dB.The reason is that there are imaginary preamble symbols in the proposed structure to increase the pseudo pilot power.Additionally,because there exists a performance platform,curves of B become flatten when SNR>15 dB.For 4-QAM modulation,the similar phenomenon occurs when SNR>20 dB.
Fig.5Normalized MSE(NMSE)performance of the proposed structure A and the IAM structure B for 4-QAM and 16-QAM in CE
In Fig.6,the CE performance of two methods are compared with each other when the number of subcarriers are respectively 256 and 512.For structure A,the curve with 512 subcarriers is lower than that with 256 subcarriers.Curves of structure B show the same results.It implies that a larger number of subcarriers means a better CE performance.Furthermore,it can be seen from Fig.6 that when BER=10?3,structure A with 256 subcarriers and 512 subcarriers are respectively better than that of structure B about 0.4 dB and 1 dB,which proves that structure A has better CE performance than structure B.
Fig.6Bit error ratio(BER)performance of the proposed structure A and the IAM structure B for different numbers of subcarriers in CE
The fundamental parameters are listed in Table 2.
Table 2Fundamental parameters of synchronization estimation
Fig.7,Fig.8 and Fig.9 respectively show the performance of the STO estimation on different channels for different numbers of subcarriers.The purpose of analyzing the coarse estimation is to be compared with the proposed STO estimation method.In Fig.7,it is obvious that the performance of the proposed method is superior to the coarse procedure on the ITU-A and ITU-B channels.The proposed method achieves a better performance on channel ITU-A than the coarse estimation on the AWGN channel.
Fig.7Performance of the STO estimation in different channels for M=1 024
Fig.8Performance of the STO estimation on different channels for M=2 048
Fig.9Performance of the STO estimation on different channels for M=4 096
For the same channel,a larger value of the threshold means a worse performance.Moreover,the convergence speed of the proposed method on channel ITU-A is slower than that on ITU-B and AWGN.The reason is that the ITU-A channel model is with mobility,time variation of the channel will be enhanced,and Doppler spread and Doppler shift will become obvious with the increased mobility.The lines in Fig.8 have the same tendency as that in Fig.7.Comparing Fig.9 with Fig.7 and Fig.8,the effect of the number of subcarriers on the performance of the proposed estimation method can be analyzed.A larger number of subcarriers means a better estimation performance.
The CFO estimation performances of three methods in different simulation environments are compared with each in Fig.10,Fig.11 and Fig.12.For convenience,methods proposed in [9],[19]and this paper are respectively named as method A,method B and method C.The computer simulations in Fig.11 show that the performance of CFO estimation in the context of the AWGN channel is the best while that in the context of channel ITU-B is the worst,no matter which method is used.For one certain method,the performance gap between different channels is more obvious when SNR>10 dB due to the multipath effects.In the context of the AWGN channel,method C is the best one to estimate CFO while the rest two methods show the similar performance,especially when SNR>14 dB.Things are different when it comes to the rest two channels,method A becomes the best one and there exist intersections between different lines.The performance gap between different lines is more obvious when SNR>20 dB.Furthermore,a larger SNR always means a better CFO estimation performance for whichever simulation environment.In conclusion,the method proposed in this paper is effective in different simulation environments and has outperformance in the CFO estimation domain.Similar conclusion can be drawn from Fig.10 and Fig.12.
Fig.10Performance of the three CFO estimators on different channels for M=1 024
Fig.11Performance of the three CFO estimators on different channels for M=2 048
Fig.12Performance of the three CFO estimators on different channels for M=4 096
The sensitivity of the OFDM/OQAM system to the STO and CFO directly results in the poor BER performance.In order to better verify the advantages of the proposed synchronization estimation method,BER performance of the proposed synchronization estimation method is compared with the improved data-aided joint CFO and time offset estimation method in[35].Fig.13 shows that the proposed method has a better BER performance than that in[35].When SNR=10 dB,the method proposed in this paper is better than the method in[35]by about 1.3 dB.The gap becomes larger with the increase of SNR.It can be seen from Fig.4 that pilot symbols+1,–1,j and–j are alternately set in the preamble structure proposed in this paper.In this structure,the interference around pilot symbols will be cancelled by each other.Therefore,the proposed structure can reduce the inter-symbol interference and increase the pseudo pilot power.Furthermore,the accuracy of the synchronization estimation will be increased and the BER performance will be better than other methods.
Fig.13BER performance of the methods proposed in this paper and[35]
This paper investigates the joint estimation of STO and CFO based on an improved GA for OFDM/OQAM systems.Specially,a new preamble structure without zero symbols is proposed at first to increase the spectral efficiency.On this basis,the MSE caused by CFO and STO is derived.The joint CFO and STO estimation is formulated into a bivariate optimization problem,and minimizing the MSE is treated as the objective.Furthermore,the pruning operator is used to improved GA and find the global solution of this optimization problem.Simulations results show that STO and CFO can be jointly estimated by the proposed method with acceptable performance.Because the preamble proposed in this paper carries no information and GA belongs to the evolutionary algorithm,how to code the pilot to improve the data transmission speed and utilize other evolutionary algorithms in joint synchronization estimation are the directions for future improvement.
Journal of Systems Engineering and Electronics2020年4期