The present invention will be better understood by reading the description of embodiments given purely by way of illustration and in no way limitative respectively with reference to the appended drawings in which:
Identical, similar or equivalent parts of the different figures described hereinafter carry the same digital references so as to facilitate moving between one figure and the next.
Reference is made first of all to
In the example in
The signals considered here are complex and sampled in the time or frequency domain.
Let s be a signal intended to be amplified by the power amplifier 2 and emerging from the amplifier 2 in the form s1 before being transmitted in a transmission channel 4. In the time domain, we have:
s
1(t)=g(∥s(t)∥).s(t) (1)
where g is the amplifier gain.
g is a function of the amplitude of the input signal s. This gain g may modify the amplitude and/or the phase of the input signal s in a non-linear way, introducing nonlinearities into the communication system 300 when the AM-AM characteristic, in other words the curve of the output signal amplitude as a function of the input signal amplitude, and/or the AM-PM characteristic, in other words the curve of the output signal phase as a function of the input signal amplitude, of the gain element g, here the amplifier 2, is not linear.
In the frequency domain, S1 is the result of the convolution between the gain G of the amplifier 2 and the incoming signal S in the amplifier:
S
1(f)=G(S(f))*S(f) (2)
The amplified signal S1 is then transmitted by one or more antennae, and passes into the transmission channel 4. At the output of the transmission channel 4 a signal S2 is obtained such that:
S
2(f)=H(f).S1(f)+η(f) (3)
where H stands for the transmission characteristic of the channel 4 in the frequency domain and η an additive white gaussian noise.
In this embodiment, prior to the transmission of the information signal S, at least one first learning sequence, or pilot, is used to estimate the characteristic of the channel 4. To do this, P1 follows the same transmission path as the information signals to be transmitted, in other words being firstly subject to an amplification stage by the power amplifier 2, then a passing through of the transmission channel 4. On reception, in other words at the output of the transmission channel 4, P1 is subject, in this embodiment, to a fast Fourier transform (FFT) (reference 6 in
After this FFT, the characteristic of the channel 4 is then subject to estimation 8. This estimation is implemented using channel characteristic estimation means 104. The amplitude of P1 is chosen preferably to be substantially constant or of small dynamic, in other words having a low PAPR (Peak to Average Power Ratio). Only the P1 phase changes. In this way, the nonlinearities present in the communication system 300 do not bias the estimation of the channel 4, since in that way, they have no or little influence over P1. The channel 4 is estimated by measuring, frequency by frequency, the difference between the sequence received at the output of the channel 4 and the initial learning sequence P1. For example, in the case of an OFDM system, the characteristic of the channel 4 is obtained directly from the ratio between the received sequence and the initial learning sequence P1, this sequence being for example an alternate set of symbols of values +1 and −1.
The nonlinearities present in the communication system 300 are subject to estimation 10, after the channel 4 estimation is obtained. In the example in
To do this, at least one pilot P2 is first amplified by the power amplifier 2 then sent through the transmission channel 4. Given that the communication system is of the OFDM type, the pilot P2 is subject to an FFT 6 by the FFT means 102, as previously for P1 during the estimation of the channel 4. This estimation 10 is implemented using nonlinearities estimation means 106.
To implement the estimation 10 of the nonlinearities, an orthonormal base of projection polynomials {p0,p1, . . . ,pn} is first calculated. This base is obtained from moments of the amplitude of the learning sequences or pilots transmitted, in other words from moments of the amplitude of initial P2 before it is sent through the power amplifier 2 and the channel 4.
After P2 has been subject to the FFT 6, coefficients {cA0,cA1, . . . ,cAn} of the AM-AM characteristic and/or coefficients {cφ0,cφ1, . . . ,cφn} of the AM-PM characteristic can be determined by projecting P2 on the orthonormal base {p0,p1, . . . ,pn}. These coefficients allow nonlinearities approximation to be obtained such that:
ĝ(∥x∥)=(cA0p0(∥x∥)+cA1p1(∥x∥)+. . . cAnpn(∥x∥)exp(j*(cφ0p0(∥x∥)+cφ1p1(∥x∥)+ . . . cφnpn(∥x∥))) (4)
This projection allows a good estimation of the whole nonlinearities characteristic to be obtained despite the presence of noise in the channel 4.
Lastly, from this approximation, the interpolation points of coordinates bi=ĝ(ai) of the nonlinearities characteristic ĝ are calculated.
The different amplitudes of the symbols forming the learning sequence P2 will be selected so as to cover as far as possible the whole dynamic of the amplitude of the signals transmitted subsequently by the communication system 300, so as to take into account all the nonlinearities sustained by the signals transmitted through the communication system 300.
The estimation of the AM-AM characteristic will now be specified.
Prior to calculating the orthonormal base of projection polynomials {p0,p1,pn}, it is appropriate to define a scalar product adapted to the vector space formed by this orthonormal base.
In the time domain, the scalar product associated with this orthonormalisation is the expected product over a certain observation time T, such that:
U,V
=E[u(t)v(t)]ε[0,T] (5)
For a system sampled at
In the frequency domain, the scalar product associated with this orthonormalisation is the expected convolution product over a certain observation frequency band B, such that:
U,V
=E[U(f)*V(f)*V(f)]fε[0,B] (7)
For a system sampled at Fe=F*N:
This scalar product selected is in fact a bilinear, symmetrical and defined positive function.
An orthonormal base of projection polynomials {p0,p1, . . . ,pn} can be obtained by the Gram-Schmidt method for the orthonormalisation of a base {1,x,x2,x3, . . . ,xn}.
Here for example are the first three polynomials obtained for the orthonormal base in the time domain:
With χ1=E[p], χ2=E[p2], χ3=E[p3], χ4=E[p4] the moments of p, in other words the moments of the amplitude of initial P2 before it is sent through the amplifier 2.
For an estimation in the frequency domain, an orhonormal base is constructed in the same way taking the formulas (7) and (8) for the scalar product calculations.
The number of polynomials forming this orthonormal base is selected as a function of the precision required in respect of estimating nonlinearities, and of the complexity of the algorithm generated by a large number of polynomials.
The channel estimation 8 implemented previously is then used to obtain y(t), the amplitude of P2 corrected by the effects of the transmission channel 4. In this way, an amplitude y(t) of the pilot P2 affected solely by the nonlinearities of the amplifier 2 is obtained.
The nonlinearities characteristic is then estimated by projecting y(t) on the orthonormal time or frequency base, previously calculated. The coefficients {cA0,cA1, . . . ,cAn} of the AM-AM characteristic are obtained by:
cAk=∥y∥,pk (9)
The calculation of the nonlinearities approximation for the AM-AM characteristic may therefore be:
In order to be able to use this approximation, a certain number of interpolation points of the nonlinearities characteristic ĝ is calculated. 5 interpolation points may for example be selected.
Let {a1,a2, . . . ,am} be the abscissas and {b1,b2, . . . ,bm} the ordinates of the interpolation points of the AM-AM characteristic.
In order to obtain a good description of the nonlinearities of the power amplifier 2, the ai may be chosen so as to have a majority of the interpolation points corresponding to the low amplitudes of the AM-AM characteristic, for example by distributing them exponentially.
In this embodiment, the ordinates of these interpolation points are given by the equation:
The Pk(x) and Pk(ai) may be calculated previously from the pilot P2 used, and the abscissas of the interpolation points ai, a multiplication-accumulation then allowing the bi to be obtained.
The interpolation points are therefore obtained at the cost of low complexity. This projection makes it possible to use all available samples of y(t), thereby obtaining a robust to noise estimation of nonlinearities.
In order to obtain an estimation of the AM-PM characteristic, by choosing the same orthonormal time base, the phase difference φ(t) is projected between the learning sequences, or pilots, received and transmitted, on this orthonormal base.
In a way similar to the AM-AM characteristic, the calculation of the nonlinearities approximation for the AM-PM characteristic may therefore be:
Let {c1,c2, . . . ,cn} be the abscissas and {d1,d2, . . . ,dn} the ordinates of the interpolation points of the AM-PM characteristic.
In this embodiment, the ordinates of these interpolation points are given by the equation:
In order to obtain a good description of the nonlinearities of the power amplifier, the abscissas ci of the in-phase characteristic may also be chosen so as to have more points for the low amplitudes, for example by distributing them exponentially.
To obtain both the AM-AM and AM-PM characteristics, it is possible either to project amplitude and phase of the received pilot successively according to the equations (10) and (12) or to use an orthonormal base of the space of the complex polynomials for the norm considered.
According to this variant, the Gram-Schmidt orthonormalisation is then carried out from a base {1,x,x*,x2,xx*, . . . }. The complex pilot received can then be projected directly on this base.
The nonlinearities approximation is therefore in this case:
The interpolation points obtained are then also complex.
Since the nonlinearities are estimated, it is now possible to implement compensation 12 for the nonlinearities sustained by an information signal s transmitted in the communications system 300.
The nonlinearities can be compensated immediately after the analogue-to-digital conversion of the signal in the communication system 300, after the transmission of signal s in the channel 4. This signal is denoted s2(t) in
It is possible to implement a compensation for the AM-AM characteristic, and then a compensation for the AM-PM characteristic.
In this case, an interpolating polynomial, for example a Lagrange polynomial, is calculated, using the interpolation points (ai,bi) to invert characteristic in amplitude of the amplifier 2. The following signal is then obtained:
The abscissas and ordinates of the interpolation points have been inverted. In this way the inverse nonlinearities characteristic is directly interpolated from inverse abscissa bi and ordinate ai interpolation points. The Lagrange polynomials are advantageous on account of their reliability in describing a bijective amplifier characteristic and their adapted implementation. Other interpolating polynomials, other than those of Lagrange, can also be selected (splines, Newton, etc.).
Using this compensation method, it is not necessary to have a high order Lagrange polynomial in order to get effective compensation for nonlinearities.
Since the AM-AM characteristic is compensated, phase distortion is then evaluated as a function of the amplitude by the equation:
This phase then just needs to be subtracted from the signal phase to correct the AM-PM characteristic and obtain at the output some compensation 12 of nonlinearities:
s
3(t)=s3′(t).exp(−j*angle(ĝ(∥s3′∥))) (17)
Instead of firstly compensating the AM-AM characteristic and then the AM-PM characteristic, it is also possible to implement a joint compensation of the AM-AM and AM-PM characteristics.
In this case, a complex Lagrange polynomial is used to interpolate the in amplitude and in phase characteristic of the amplifier:
where the pairs (ai,bi) are obtained after projection on a base of complex polynomials.
The gain in terms of bit error rate is shown in
The signal whereof the nonlinearities have been compensated can then be passed into the frequency domain by an FFT 6 (
One embodiment of the invention is compatible with a method for restricting or clipping the power to be transmitted upstream of the power amplifier 2, and a method for compensating the distortions caused by such clipping at the receiver end. To do this, the maximum amplitude of the learning sequence or pilot is defined as being lower than the clipping level.
The method according to one embodiment of the invention can for example be implemented in communication systems of the OFDM, WCDMA, EDGE type, or more generally in all radio, wire, optical communication systems etc.
Number | Date | Country | Kind |
---|---|---|---|
06 07147 | Aug 2006 | FR | national |