1. Field of the Invention
The present invention relates to a carrier synchronization and channel equalization method, and in particular to a joint carrier synchronization and channel equalization method for OFDM systems.
2. The Prior Arts
In the prior art, a device configuration of a baseband equivalent model for Orthogonal Frequency Division Multiplexer (OFDM) is shown in
wherein, T=Tu/N represents a sample interval.
In the structure mentioned above, a channel impulse response of a multi-path fading channel 16 is represented by h(t)=Σihi(t)·δ(t−τi), wherein, hi (t) and τi represent, respectively, an attenuation and a delay spread of the ith path, and then the time domain signal samples xn,l are serially sent into a Digital-to-Analog converter (DAC) 14, and then they are transmitted into channel 16, such that channel noises n(t) exists in channel 16. As such, the channel output can be expressed by y(t)=Σihi (t)·x(t−τi)+n(t), wherein, n(t) is a white Gaussian noise with its expectation value as zero-mean. In this framework, the timing frequency offset between a Digital-to-Analog converter (DAC) 14 and an Analog-to-Digital converter (ADC) 22 is assumed to be ideal. The output signal of the channel is damaged by a carrier frequency offset (CFO) effect, thus upon being sampled by an Analog-to-Digital converter (ADC) 22, the nth reception signal sample of the lth OFDM symbol can be expressed by the following formula:
y
n,l
=y(t)·ej2πΔft|t=l(N+Ng)T+NgT+NT, (2)
formula (2) explains that a carrier frequency offset (CFO) Δf induces a linear increment of phase offsets in the time domain signals.
Subsequently, upon removing the guard interval from the reception signal samples, the remaining reception signal samples are sent into DFT for demodulation processing. Therefore, the kth sub-channel signal of the ith OFDM symbol can be expressed by the following formula:
Y
k,l
=S
k,l
+I
k,l
+N
k,l (3)
wherein, Sk,l, Ik,l and Nk,l represent a signal term, inter-carrier interference (ICI), and white Gaussian noise on the kth sub-channel respectively. Moreover, Sk,l and Ik,l can be derived as follows:
wherein Hk,l is a channel response of the kth sub-channel, and it must satisfy the stationary property in an OFDM symbol. In addition, a local subcarrier frequency offset φqk and an attenuation factor si(πφqk) can be expressed by the following formula:
wherein, ε=ΔfNT is a normalized CFO, and it represents the amount of Δf transferred from the time domain to the frequency domain through DFT. It is evident that the attenuation factor and ICI term are proportional to ε, as shown in formulae (5), (6.1) and (6.2).
In view of the fact that ε in formulae (6.1) and (6.2) is extremely small, while the system enters into a tracking stage, then si(πφkk) is very close to 1 and si(πφqk) almost approaching zero. As such, si(πφkk) in formula (4) can be ignored, and the ICI term in formula (3) can be eliminated. Upon performing the simplification mentioned above, an equivalent channel response {tilde over (H)}k,l on the kth sub-channel can be expressed in polar coordinate as follows:
wherein, θH
In addition, Gk,l=GH
Finally, formula (3) can be rearranged as
Y
k,l
={tilde over (H)}
k,l
·X
k,l
+N
k,l (8)
In order to resolve the adverse effects caused by the CFO and the channel distortion to a received signal, thus a carrier synchronization and a channel equalization techniques are proposed to overcome these problems. In a framework of the prior arts, the carrier synchronization technique is realized by a frequency control loop having an individual frequency detector. However, such a frequency estimation mechanism is not a perfect process. In practice, the carrier frequency jitter will not be zero. Therefore, such a phenomenon will result in the constellation rotation on each sub-channel in an OFDM system, hereby further degrading the system performance. In a practical OFDM system, a carrier phase compensation on each sub-channel is necessary to overcome the constellation rotation.
In addition, in OFDM transmission system, in general, a channel estimation method utilized on each sub-channel is based on the least square (LS) algorithm. However, this method is not very accurate. The residual CFO will destroy the accuracy of the channel estimation on each sub-channel since the residual CFO has not been fully removed. As such, when the carrier synchronization process enters into a tracking stage, the channel information on each sub-channel has to be updated by the least-mean square (LMS) algorithm to track the channel variations.
For the reasons mentioned above, in general, the carrier frequency synchronization and the channel equalization are restrained based on an individual cost function. Furthermore, the mutual interference will occur between the carrier frequency synchronization and the channel equalization to degrade the system performance, namely increase Bit Error Rate (BER) for an OFDM system.
As such, presently, the performance of the carrier frequency synchronization and the channel equalization techniques of the prior art is still not quite satisfactory, and it has much room for improvements.
In view of the problems and shortcomings of the prior arts, the present invention provides a joint carrier synchronization and channel equalization method for OFDM systems, so as to solve the problems of the prior arts.
A major objective of the present invention is to provide a joint carrier synchronization and channel equalization method for OFDM systems, wherein, a single cost function mechanism is employed to minimize the error power on each sub-channel and to further maximize the sub-channel SNR concurrently. The joint method is used to obtain the phase error and the gain error of each sub-channel signal simultaneously, hereby to solve the problems caused by the CFO and the channel distortion, further to increase the system performance of an OFDM system, and also to minimize Bit-Error-Rate (BER) for the OFDM system.
In order to achieve the above-mentioned objectives, the present invention provides a joint carrier synchronization and channel equalization method, that is suitable for use in a receiver of an OFDM system, including the following steps: firstly, receiving a reception signal sample of an OFDM symbol, then outputting a sub-channel signal in the frequency domain, while obtaining a phase error and a gain error of each sub-channel simultaneously. Next, obtaining an execution phase compensation factor based on phase error on each sub-channel, and obtaining an execution carrier frequency offset factor after calculating an average value of the summed-up phase error for all sub-channels, then obtaining an execution gain compensation factor based on a gain error on each sub-channel. Subsequently, through first multiplying an execution carrier frequency offset factor by a reception signal sample of a next OFDM symbol, thus eliminating the phase offset in a reception signal sample in time domain. Finally, in a next step, the magnitude distortion and the phase distortion on each sub-channel signal in the frequency domain are compensated, respectively, by using the sub-channel signal generated by multiplying an execution gain compensation factor and an execution phase compensation factor by a reception signal sample of a next OFDM symbol.
Further scope of the applicability of the present invention will become apparent from the detailed description given hereinafter. However, it should be understood that the detailed description and specific examples, while indicating preferred embodiments of the present invention, are given by way of illustration only, since various changes and modifications within the spirit and scope of the present invention will become apparent to those skilled in the art from this detailed descriptions.
The related drawings in connection with the detailed description of the present invention to be made later are described briefly as follows, in which:
The purpose, construction, features, functions and advantages of the present invention can be appreciated and understood more thoroughly through the following detailed descriptions with reference to the attached drawings.
Firstly, referring to
In the following descriptions referring to
Furthermore, the outer loop includes an average phase error calculator 36, a second loop filter 38, an adder 40, a second numerically-controlled oscillator (NCO) 42, a de-rotator 44, a discrete Fourier transformation 46. Wherein, all the elements included in an outer loop are connected in the above-mentioned sequence, and an average phase error calculator 36 is connected to a phase error detector 28 on each sub-channel in the frequency domain.
In addition, a gain equalization loop is also constructed on each sub-channel in the frequency domain. The gain equalization loop includes a signal slicer 24, a decision error element 26, a gain error detector 48, a third loop filter 50, and a multiplier used as a gain equalization compensator 52. Wherein, the signal slicer 24 and the decision error element 26 are shared by an inner loop, and signal slicer 24 are connected in parallel with decision error element 26; and the signal slicer 24 are connected to the following devices in a sequential manner: a gain error detector 48, a third loop filter 50, a gain equalization compensator 52, and an inner loop compensator 34; finally, the inner loop compensator 34 is connected to signal slicer 24. Moreover, the gain equalization compensator 52 is connected to the discrete Fourier transformation 46 of the outer loop.
Through the application of the present invention, the problems of the residual carrier frequency offset (CFO) and the channel equalization can be solved simultaneously. From the system performance point of view, as long as the power of decision error signal on each sub-channel can be minimized, then the system performance can be enhanced. Therefore, based on the minimum mean square error (MMSE), the cost function J(•) of the joint carrier synchronization and channel equalization method is presented to minimize the power of the decision error signal on each sub-channel, as represented by the following formula:
wherein, E[•] is an expectation operator. k represents the sub-channel index in the frequency domain. Ek and Ek represent respectively a decision error value of a decision error signal and the power of decision error signal on the kth sub-channel, and Ŷk is a equalized sub-channel signal of a sub-channel signal Yk after the gain compensation and the phase compensation as performed by the gain equalization compensator 52 and the inner loop compensator 34. {circumflex over (X)}k is the sub-channel decision signal as performed by a signal slicer 24. Ĝk and {circumflex over (θ)}k are the compensated gain and the compensated phase on the kth sub-channel in the frequency domain respectively, and the derivations of the signals and parameters mentioned above will be described in more detail in the following. Conceptually, the signal-to-noise ratio (SNR) on the kth sub-channel can be expressed as
wherein, Sk is the signal power on the kth sub-channel. In formulae (9) and (10), it is described that, the accurate phase and gain factors are able to make the power of the decision error signal on the kth sub-channel minimum as well as the signal-to-noise ratio (SNR) on the kth sub-channel maximum simultaneously. Therefore, the closed-loop control technique is used to obtain the accurate phase compensation factor and gain compensation factor on the kth sub-channel, hereby realizing a method of joint carrier synchronization and channel equalization.
Subsequently, referring to
Furthermore, the initial gain and phase factors on a kth sub-channel can be acquired based on a training symbol, such as the long preamble in IEEE 802.11a, and represented as
where Xk and Yk are the transmitted and the received training symbols on the kth sub-channel respectively.
Therefore, upon receiving signals at the first time, as shown in Step S10, the de-rotator 44 is used to receive a first reception signal sample of a first OFDM symbol with CFO effect, namely, at this time, the yn in
Subsequently, upon completing the execution of step S16, then proceeding with the execution of step S18. At this time, a signal slicer 24 is used to receive a first received sub-channel equalization signal from the output of the phase compensator of the inner loop 34, and perform a slicing processing for a first received sub-channel equalization signal, hereby outputting a sub-channel decision signal {circumflex over (X)}k. Then, as shown in step S20, a decision error element 26 is used to receive a sub-channel equalization signal and a first received sub-channel decision signal, and output a decision error signal after subtracting a first received sub-channel equalization signal from the sub-channel decision signal. Furthermore, as shown in step S22, a phase error detector 28, and a gain error detector 48 are utilized separately to detect and obtain simultaneously a phase error and a gain error of each sub-channel signal in the frequency domain based on formulae (3) and (4) by using a decision error signal and a sub-channel decision signal, and outputting the phase error and the gain error thus obtained.
εk,p=Im({circumflex over (X)}kEk*), (13)
εk,g=Re({circumflex over (X)}kEk*), (14)
wherein, Re(.), Im(.) and * represent respectively a real-part, an imaginary-part, and a complex-conjugate operators, {circumflex over (X)}k is a decision signal on the kth sub-channel, and Ek is a decision error on the kth sub-channel corresponding to the decision signal {circumflex over (X)}k and the equalized signal Ŷk.
Moreover, upon completing the execution of step S22, proceeding with the execution of step S24. At this time, proceeding with obtaining an execution carrier frequency offset factor Δ{circumflex over (f)}={circumflex over (f)}initial+{circumflex over (f)}r, an execution phase compensation factor e−j{circumflex over (θ)}
A first loop filter 30 of an inner loop is used to receive a phase error output by a phase error detector 28, and then output a compensation phase {circumflex over (θ)}k thus obtained. Then, a first numerically-controlled oscillator (NCO) 32 is used to receive this compensation phase, and calculate to obtain an execution phase compensation factor e−j{circumflex over (θ)}
An average phase error calculator 36 of an outer loop is used to receive a phase error output by a phase error detector 28 of each sub-channel in the frequency domain, and the phase errors for all sub-channels are summed up according to formula (15), then the result of this summation is an average phase error.
wherein, k represents the sub-channel index in the frequency domain, εk,p is a phase error on the kth sub-channel,
Subsequently, a second loop filter 38 is used to receive this average phase error, and output an estimated residual CFO thus obtained, namely {circumflex over (f)}r. Then, an adder 40 is used to receive this estimated residual CFO {circumflex over (f)}r and an initial CFO {circumflex over (f)}initial, and add them together to obtain a compensation CFO, namely Δ{circumflex over (f)}={circumflex over (f)}initial+{circumflex over (f)}r. Finally, a second numerically-controlled oscillator (NCO) 42 is used to receive this compensation CFO, and calculate to obtain an execution carrier frequency offset factor e−jφ
In addition, a third loop filter 50 of gain equalization loop is used to receive a gain error output by a gain error detector 48, and then output a execution gain compensation factor Ĝk of a real number thus obtained.
Upon obtaining all the execution gain compensation factor, the execution phase compensation factor, and the execution carrier frequency offset factor, then proceeding with execution of step S26. At this time, a de-rotator 44 will receive again a reception signal sample of an OFDM symbol in the second time. Herein, the signal is defined as a second reception signal sample of a second OFDM symbol, and at this time, the yn in
ŷn,l=yn,le−jφ
Then, upon completing the execution of step S26, then proceeding with the execution of step S28. At this time, discrete Fourier transformation 46 is used to receive a second reception signal sample from de-rotator 44, and perform discrete Fourier transformation on the second reception signal sample, such that outputting a second received sub-channel signal on each sub-channel in the frequency domain, namely, at this time, the Yk in
Subsequently, as shown in step S30, a gain equalization compensator 52 is used to receive a second received sub-channel signal and an execution gain compensation factor, then multiply the signal and the gain compensation factor, thus compensating the magnitude distortion of a second received sub-channel signal on each sub-channel in the frequency domain.
Afterwards, proceeding with execution of step S32, during which, an inner loop compensator 34 is used to receive an execution phase compensation factor, and receive a second received sub-channel signal having its magnitude distortion compensated from a gain equalization compensator 52, multiply the signal and the received execution phase compensation factor, and then output a second received sub-channel equalization signal, namely, at this time, Ŷk in
Upon receiving again a reception signal sample of a next OFDM symbol, the de-rotator 44 may proceed with phase and gain compensations for a reception signal sample of a next OFDM symbol in a same manner as utilizing a first reception signal sample in compensating a second reception signal sample. Namely, proceeding with compensation by using the carrier frequency offset factor e−j2πΔ{circumflex over (f)}nT, phase compensation factor e−j{circumflex over (θ)}
In order to derive the phase variation on each sub-channel, the phase transfer function of each loop must be derived as described below, meanwhile, referring to
In the inner loop, the first loop filter 30 is a proportional controller with gain κp. Based on this result, the inner loop is constructed by utilizing a type-1 digital phase-locked loop (DPLL). Actually, the inner loop is realized on each sub-channel in the frequency domain. Wherein, the “type” means the number of integrators within a closed-loop. The phase transfer function Hi(z) and phase error transfer function Ei(z) of the inner loop can be derived as follows:
wherein θ(z) is a phase of a sub-channel signal output by a discrete Fourier transformation 46, {circumflex over (θ)}(z) is a phase of a sub-channel signal output by a first numerically-controlled oscillator (NCO) 32, and θe(z)=θ(z)−{circumflex over (θ)}(z). κl=κdκloκp is an open loop gain. κd is a gain of a phase error detector 28. κlo is a gain of a first numerically-controlled oscillator (NCO) 32 of the inner loop. The subscripts k as shown in Hi(z) and Ei(z) are omitted since κl is the same for all sub-channels. The stability condition of the inner loop Hi (z) must satisfy the condition 0<κi<2, since in this condition the pole of Hi (z) is located within the unit circle.
Since the gain equalization loop on each sub-channel is a first order equalization loop, and its loop equation can be expressed as
Ĝ
k(l)=Ĝk(l−1)+μgεk,g(l−1) (18)
wherein, l is a symbol index. μg is an open loop gain step-size. εk,g is a gain error on the kth sub-channel, and gain compensation factor Ĝk is a real number.
For a gain equalization loop on each sub-channel, the gain compensation factor on a gain equalization loop must be updated for every OFDM symbol. As such, in a tracking stage, the gain error will gradually decreases to a very small value. The closed-loop transfer function of the first order loop can be expressed as
wherein, the subscript k is omitted since the open loop gain μg is the same for all sub-channels. Evidently, the stability condition of Hg(z) must satisfy the condition 0<μg<2, because in this condition the pole is located within the unit circle.
In the following description, please refer to
The outer loop is operated over two-rate regions, namely, sample-rate, and symbol-rate. In addition, K inner loops are located in the outer loop, and the carrier phase error caused by CFO is identical for all sub-channels. Therefore, K inner loops can be simplified into a single loop to derive the outer loop transfer function. In addition, the average phase error can be ignored in the signal flow graph as shown in
According to a domain transformation as shown in
wherein, {circumflex over (φ)}(z) is the phase information of the output of a second numerically-controlled oscillator (NCO) 42 in the time domain, φe(z) is the phase error of the output of de-rotator 144 in the time domain, and κ2=κdκoκi is an open loop gain of the outer loop. κo is a gain of a second numerically-controlled oscillator (NCO) 42 of the outer loop. Based on F(z), the phase transfer function Ho(z) and the phase error transfer function Eo (z) of the outer loop can be derived as
wherein, φ(z) is a phase information of each signal sample of an received OFDM symbol in the time domain.
Moreover, the dual-loop carrier synchronization device can not only track CFO in time domain by using the outer loop, but it can also recover the residual CFO jitter and the sub-channel phase distortion through utilizing the inner loop on each sub-channel in the frequency domain. Therefore, the phase transfer function Hd(z) and the phase error transfer function Ed(z) of the dual-loop carrier synchronization loop can be expressed as
Subsequently, referring to
Finally, referring to
Number | Date | Country | Kind |
---|---|---|---|
098109273 | Mar 2009 | TW | national |