A NEW INTERPOLATION ALGORITHM FOR PILOT-BASED CHANNEL ESTIMATION IN OFDM SYSTEMS
Field of the invention
[001 ] The invention relates to the field of multicarrier communications systems and more particularly to Orthogonal Frequency-Division Multiplexing (OFDM) systems, including wireless OFDM systems.
Background of the invention
[002] Orthogonal Frequency-Division Multiplexing (OFDM) is also referred to as
"multicarrier modulation" (MCM) or "discrete multi-tone modulation" (DMTM). OFDM splits up and encodes high-speed incoming serial data, modulating it over a plurality of different carrier frequencies (called "subcarriers") within a communication channel to transmit the data from one user to another. The serial information is broken up into a plurality of sub-signals that are transmitted simultaneously (in parallel) over the subcarriers.
[003] By spacing the subcarrier frequencies at intervals of the frequency of the symbols to transmit, the peak power of each modulated subcarrier lines up exactly with zero power components of the other modulated subcarriers, thereby providing orthogonality (independence and separability) of the individual subcarriers. This allows a good spectral efficiency (close to optimal) and minimal inter-channel interference (ICI), i.e. interference between the subcarriers.
[004] For these reasons, OFDM is used in many applications. Many digital
transmission systems have adopted OFDM as the modulation technique such as digital broadcasting terrestrial TV (DVB-T), digital audio broadcasting (DAB), terrestrial integrated services digital broadcasting (ISDB-T), digital subscriber line
(xDSL), WLAN systems, e.g. based on the IEEE 802.1 1 a/g standards, cable TV systems, etc.
[005] In an OFDM system, the transmitter modulates the message bit sequence into PSK/QAM symbols, performs IDFT (Inverse Discrete Fourier Transform) on the symbols to convert them into time-domain signals, and then sends them out through a transmission channel (in general a wireless channel).
[006] The received signal is usually distorted by the channel characteristics. In order to recover the transmitted bits, the channel effect must be estimated and compensated in the receiver; this is called channel equalization.
[007] As known by the OFDM theory, each subcarrier (i.e. each input of the Inverse DFT in the transmitter) can be regarded as an independent channel, as long as no ICI (Inter-Carrier Interference) occurs, and thus preserving the orthogonality among subcarriers. The orthogonality allows each subcarrier component of the received signal to be expressed as the product of the transmitted signal and channel frequency response (transfer function) at the subcarrier. Thus, the transmitted signal can be recovered by estimating the channel response just at each subcarrier.
[008] In general, the channel can be estimated by using a preamble symbols or pilot subcarriers known to both transmitter and receiver, which employ various interpolation techniques to estimate the channel response of the subcarriers between pilot tones.
[009] Preamble symbols are training data where all subcarriers are used for channel estimation. The preamble is usually transmitted at the beginning of each OFDM packet.
[010] Pilot subcarriers are dedicated subcarriers for channel estimation. The other subcarriers can be used for data transmission; Pilot subcarriers are transmitted in the same time as OFDM symbols.
[01 1 ] In general, a data signal or a training signal, or both, can be used for channel estimation. In order to choose the channel estimation technique for the OFDM system under consideration, many different aspects of implementations, including the required performance, computational complexity and time-variation of the channel must be taken into account.
[012] Several standards specify the use of channel estimation using pilot subcarriers.
For instance the WLAN 802.1 1 a/g standard is based on 4 pilot subcarriers that are sent with each OFDM symbol. The pilot subcarriers are spaced in frequency by having them transmitted by subcarriers -21 , -7, +7 and +21 .
[013] The data transmitted over each of the pilot subcarriers is known to the receiver.
In each symbol received by the receiver, the outputs of the DFT (Discrete Fourier Transform) corresponding to the subcarrier frequencies -21 , -7, +7 and +21 have phases and amplitudes that are known to the receiver (in particular, the amplitude is set to a "unit amplitude" and the phase is either 0 or 180 degrees).
[014] Based on these pilot subcarriers, the channel estimation can be performed in the frequency domain by the receiver, i.e. after the DFT module. For each subcarrier, the channel equation is: Y(k)=C(k).X(k)+Z(k)
[015] Where
• k is the subcarrier (or sub-channel) frequency index,
• X(k) is the subcarrier signal of kth index at the input to the
Inverse Digital Fourier Transform (IDFT) module of the
transmitter (i.e. the transmitted OFDM symbol k in the
frequency domain)
• Y(k) is the corresponding signal output from the Digital Fourier Transform (DFT) module of the receiver (i.e. the received OFDM symbol k in the frequency domain),
• C(k) is the frequency response (or transfer function) at the frequency of the kth subcarrier (i.e. kth sample of the channel transfer function), and,
• Z(k) is the noise frequency spectrum. If the DFT input noise is white, the output noise Z(k) is also white.
[01 6] For channel estimation, the components X(k)s are known pilot subcarriers with unit amplitude and phase at π or -π (i.e. a BPSK signal is sent over each pilot subcarrier).
The subcarrier response C (k) can be directly estimated as an interpolation of the product X'(kp ). Y (kp) , where kp corresponds to pilot indices. In the example of a WLAN 802.1 1 a/g based implementation, these pilot indices can take as values: -21 , -7, +7 and +21 .
[01 7] In order to find the other C"(k) , an interpolation (and an extrapolation) can be made on the basis of the know C"(k ). This estimates the channel frequency response p in the pass-band of the OFDM transmitted signal. In a WLAN
802.1 1 a/g, this pass-band corresponds to 52 samples spaced by 31 2.5 kHz.
[01 8] Once the transfer function C (k) is determined, it can be used for channel
equalization. This channel equalization process can consist in dividing the previously described equation Y(k)=C(k).X(k)+Z(k) by C(k).
[019] We can then write:
Y(k) X(k) Z(k)
= X(k) + noise
[020] However this technique is not very satisfactory, as the "noise" is not easy to
estimate and the interpolation/extrapolation error rate is often too high. There is a therefore a need for a solution improving the situation.
Statement of invention
[021 ] Accordingly the present invention provides an OFDM receiver to receive a
plurality of signals transmitted in parallel via multiple subcarrier channels, wherein
at least one subcarrier is a pilot subcarrier to receive a pilot signal and at least one subcarrier carries a data subcarrier to carry a data signal comprising: a pilot subcarrier transfer function estimation module arranged to receive the pilot signal and to estimate a transfer function for the pilot subcarrier response;
an extrapolation module responsive to the estimated pilot channel transfer function to extrapolate a transfer function for the channel response of the data subcarrier; and
an equaliser module responsive to the extrapolated data subcarrier transfer function to equalise a data signal received via the data sub carrier;
characterised in that the pilot sub carrier estimation module comprises a
Cartesian to polar converter to convert the received pilot signal from Cartesian to polar form before the pilot subcarrier transfer function estimation module estimates the pilot subcarrier transfer function.
[022] The direction of phase interpolation may be determined by a module which
determines a ratio of the difference in the phase of the pilot signals and the
difference in frequency. The required direction of phase interpolation is indicated according to whether the ratio is positive or negative.
[023] According to a second aspect of the present invention there is provided a method of equalising an OFDM signal received by an OFDM receiver comprising the steps of:
simultaneously receiving a pilot signal transmitted via a pilot sub carrier and a data signal transmitted via a data subcarrier;
converting the pilot signal from Cartesian to polar coordinates;
calculating a pilot channel transfer function for the pilot channel;
extrapolating a data channel transfer function from the pilot channel transfer function; and
equalising the data channel signal using the extrapolated data channel transfer function.
[024] The receiver may comprise an antenna and other conventional components to communicate the received pilot and data signals to a general purpose processor running machine readable code to implement the method and receiver of the invention. Accordingly the invention may be embodied in a package of machine readable code arranged for download into a device in which the steps of the receiver method and the components of the receiver are to be implemented by execution of the code. Further the invention may be embodied in a recording of the package of machine readable code.
[025] Converting the received frequency bins to the complex polar domain facilitates separate interpolation of the phase and amplitude for each subcarrier.
BRIEF DESCRI PTION OF THE DRAWI NGS
[026] Embodiments of a receiver arranged to implement the method of the present invention will now be described, by way of example only, with reference to the accompanying figures, wherein:
Fig. 1 is a diagram of an OFDM system comprising a transmitter and receiver.
DETAI LED DESCRI PTION OF EM BODIM ENTS OF THE I NVENTION
[027] FIG.1 illustrates a simplified and high-level block diagram of an OFDM system comprising a transmitter and a receiver (1 ) connected by a communication channel TC.
[028] The OFDM transmitter TX transmits several symbols (or signals) in parallel by assigning each of them to a different subcarrier, each carrier being orthogonal to the others. The number of carriers (or subcarriers) depends on the total bandwidth and on the duration of a symbol (i.e. the size of the time window in which the receiver can capture the transmitted symbols). The transmitter is substantially conventional.
[029] The symbols to be transmitted are represented as X0-XNC-I are initially in the frequency domain.
[030] The symbols X1 [k] to be transmitted at emitter EMT side are first appropriately encoded and then sent to a serial-to-parallel transformer (not depicted). This function block aims in extracting from a symbol, a number N of individual values X1 [k], where "k" ranges from 1 to N-1 and "1 " represents the ordinal number of the OFDM symbol to be transmitted.
[031 ] This number N is equal to the number of subcarriers (s(n)) used for the
transmission over the communication channel TC. The number of subcarriers
used is generally selected depending on the standard to which the transmitter receiver systems complies. For instance according to WLAN 802.1 1 a (on which relies the embodiment depicted on figure 5a), 52 subcarriers are used (48 data subcarriers + 4 pilot subcarriers), and 12 (6 on the "left" and 6 on the "right") subcarriers are not used (zero subcarriers).
[032] Here and in the following description, the variable k will be used to refer to the frequency domain and the variable n for time domain.
[033] The symbols Xl[k] which are output by the serial-to-parallel convertor are
provided to a frequency-to-time conversion module IDFT. This frequency-to-time convertor can perform an inverse discrete Fourier Transform of the symbols Xl[k] to generate symbols xl[n] in the time domain. These symbols xl[n] are then multiplexed by the parallel-to-serial transformer (MUX), to produce a signal transmitted over the transmission channel TC. Some cyclic prefixes are also added at the input of the transformer MUX. NCp is the number of cyclic prefixes.
[034] The multiplexer MUX modulates the orthogonal subcarriers by the symbols to be transmitted.
[035] A conventional multiplexer (MUX), transposes the signals into the time domain and modulates the subcarriers (s(n) for transmission over the communication channel TC.
[036] The receiver 1 includes a demultiplexer (DEMUX) 2 which transposes the
received signals back to the frequency domain to extract the transmitted symbols and noise.
[037] The transmitted OFDM symbol s(n) signal spectrum is the sum in the frequency domain of the orthogonal subcarrier sine functions that are superposed over each other. The individual symbols can be independently modulated by using different
types of modulations techniques, like for instance QAM (Quadrature Amplitude Modulation) or PSK (Phase-Shift Keying).
[038] The transmission channel TC is associated with a transfer function C(k) and an impulse response c(n). The transfer function represents the channel multipath behaviour, delay, and complex attenuation. The transfer function may vary over time and is therefore indexed by the number n of the transmitted symbol.
[039] Noise z(n) can also be introduced during the transmission over the
communication channel TC.
[040] The receiver 1 comprises receiving circuits, particularly a demultiplexer 2
(DEMUX) and a time to frequency conversion module 3 (DFT) which generates demodulated symbols Y1 [k] from the received symbols stream. The demodulated symbols can then be passed over to other circuitry of the receiver RCV (not depicted).
[041 ] The receiving circuitry may perform different steps which are usual in OFDM systems but may be subject to variations and different implementations. The receiver 1 receives the signal and performs a serial to parallel transformation in the serial-to-parallel transformer or demultiplexer 2 (DEMUX), which parallelizes each received symbol to a set of time symbols Y1 (n). It can also remove the cyclic prefixes.
[042] A time to frequency conversion is performed by the time-to-frequency conversion module 3 (DFT) performs a discrete Fourier Transform of the time symbols y1 (n) to generate a set of received frequency bins Yi (k) in the frequency domain, corresponding to determined subcarriers.
[043] At a fourth step a parallel-to-serial transformer (not depicted) produce symbols Yi(k) made of these received frequency bins Yi(k). These symbols can then be transmitted to other circuits of the receiver RCV (also not depicted in the figure). [044] In order to perform these steps, the receiver 1 shares the same subcarrier
frequencies, as well as the number of subcarriers N with the transmitter. They may have been communicated before the transmission, or they may be set or tuned previously in both the transmitter and receiver according to standards specifications.
[045] In addition to the receiving circuitry, the receiver 1 comprises an estimation circuit 4. [046] The estimation circuit 4 may comprise several functional modules. The functional modules are herein described separately for clarity but may be implemented on a single device or possibly even in software.
[047] A Cartesian to polar module 5 (M1 ) receives the frequency bins vectors Y0 -YNC-I in Cartesian form from the DFT 3 and performs the step of converting the frequency bins vectors to complex polar (cordic) values. A channel transfer function calculator module 6 (M2) receives the pilot subcarriers Y(-21 ), Y(-7), Y(7) and Y(21 ) from the Cartesian to polar module and computes a channel transfer function (C(-21 ), C(-7), C(7) and C(21 ) for the example of WLAN 802.1 1 a/g reception. The pilot subcarriers may be different for other applications. In this case the transfer function module 6 the transfer functions as follows: a. C (-21 ) = *(-21 ). V(-21 ) b C -7) = X'{-7). Y {-7) c. C 7) = X'{7). Y {7)
d. C (21 ) = X*(21 ). V (21 )
[048] These values are then provided to an extrapolation module 7 which performs polar interpolation (and extrapolation) to find all (or part of) the other channel transfer function samples, i.e. the values of C (k) of the other values of k ie ke[-21 , -7,7,21 ].
[049] Then the values of C (k) can be output by the extrapolation module 6 and
provided to a channel equalizer (not depicted). This channel equalizer can be conventional.
[050] By simulation and experience, it has been determined by the inventor that
performing the interpolation in the complex polar domain provides better results than the conventional process of interpolation in the complex Cartesian domain. The process of interpolation in the complex polar domain produces, relative to interpolation in the complex Cartesian domain, a reduced error rate and better channel estimation accuracy.
[051 ] The channel distorts the transmitted signal by distorting its amplitude (channel attenuation) and its phase (channel delay). Accordingly channel distortion occurs in the Polar form of the signal and not in its Cartesian form.
[052] In addition, according to an embodiment of the invention, further functional modules can be put in place in order to compute an interpolation of the phase. These modules M4, M5, M6 avoid doing phase interpolation errors caused by channel frequency response phase wrapping.
[053] The quantity X*(k)Y(k) can be interpolated in the complex polar domain by
separately interpolating the amplitude and the phase, and then forming the complex vector C(k) from the interpolated amplitude and interpolated phase. Interpolating the amplitude has been described previously.
[054] A phase extraction module 7 (M4) receives the frequency bins Yi(k) converted into polar coordinates from the Cartesian to polar module 5 and extract from them the phase A9(k).
[055] The phase of each frequency bin vector (Y1 k) is then provided to a module M5 in charge of computing the values A6(k)/Ak.
[056] These values are then provided to a channel phase wrapping detector 8 (M6).
This module is in charge of detecting and localizing each sub-channel at which the phase of the channel frequency response wraps from -π to +π or vice versa. Once each sub-channel is localized, the phase interpolation takes that into consideration in order to do the phase interpolation in a correct way.
[057] If A0(k)/Ak between pilots is positive, a clockwise phase interpolation is to be done.
[058] If A0(k)/Ak between pilots is negative, an anti-clockwise phase interpolation is to be done
[059] Channel phase wrapping is a phenomenon that happens in multipath channels where a Rayleigh effect is provoked. In this case the channel impulse response is complex and not real anymore. Being complex, and because the phase is a value that can wrap, while interpolating the phase between two consecutive pilot subcarriers we do not know if the interpolation must be done in the clockwise or anticlockwise direction. By this invention, doing polar domain interpolation becomes possible thanks to the idea of the channel wrapping detector.
[060] When the channel transfer functions are interpolated/extrapolated they are
applied to an equaliser module 1 1 to equalise the subcarriers.
The invention has been described with reference to preferred embodiments. However, many variations are possible within the scope of the invention.