Linearization of Non-Linear Amplifiers

Abstract

A linearization device (380) is disclosed, which is configured to determine pre-distortion parameters associated with a plurality of non-linear amplifiers (331, 332, 333, 334). Each of the non-linear amplifiers is associated with one of a plurality of transmit antenna elements and with a non-linear transfer function defining an output of the non-linear amplifier based on an input of the non-linear amplifier and based on a reflection signal for the non-linear amplifier, resulting from mutual couplings among the plurality of transmit antenna elements. The linearization device comprises a first port (381), a second port (382), and determination circuitry (383). The first port is configured to receive a plurality of channel coefficients indicative of channel characteristics of a plurality of communication paths between the plurality of non-linear amplifiers and two or more transmit observation receivers (370, 371, 372). Each transmit observation receiver is configured to receive a sum of transmission signals generated by the plurality of non-linear amplifiers and transferred over the communication paths between the plurality of non-linear amplifiers and the transmit observation receiver. The second port is configured to receive the sums of transmission signals from the transmit observation receivers. The determination circuitry is configured to determine the pre-distortion parameters based on the received plurality of channel coefficients, the received sums of transmission signals, and a model of the non-linear transfer functions of the non-linear amplifiers.

Claims

1-20. (canceled)

21. A linearization device configured to determine pre-distortion parameters associated with a plurality of non-linear amplifiers, each associated with one of a plurality of transmit antenna elements and with a non-linear transfer function defining an output of the non-linear amplifier based on an input of the non-linear amplifier and based on a reflection signal for the non-linear amplifier, resulting from mutual couplings among the plurality of transmit antenna elements, the linearization device comprising: a first port configured to receive a plurality of channel coefficients indicative of channel characteristics of a plurality of communication paths between the plurality of non-linear amplifiers and two or more transmit observation receivers, wherein each transmit observation receiver is configured to receive a sum of transmission signals generated by the plurality of non-linear amplifiers and transferred over the communication paths between the plurality of non-linear amplifiers and the transmit observation receiver; a second port configured to receive the sums of transmission signals from the transmit observation receivers; and determination circuitry configured to determine the pre-distortion parameters based on the received plurality of channel coefficients, the received sums of transmission signals, and a model of the non-linear transfer functions of the non-linear amplifiers.

22. The linearization device of claim 21, wherein the model is represented, for each of the non-linear amplifiers, by a plurality of amplifying coefficients in a space spanned by a regression matrix of the inputs and the reflection signals.

23. The linearization device of claim 22, wherein each reflection signal is modeled as a linear function of the outputs of the plurality of non-linear amplifiers, the linear function being defined by reflection coefficients.

24. The linearization device of claim 23, wherein the determination circuitry is further configured to determine the reflection coefficients by: defining initial reflection coefficients; determining intermediate amplifying coefficients based on the initial reflection coefficients; determining intermediate reflection coefficients based on the intermediate amplifying coefficients; and recursively determining refined amplifying coefficients based on previous reflection coefficients, and refined reflection coefficients based on previous amplifying coefficients.

25. The linearization device of claim 22, wherein the linearization device is configured to determine the pre-distortion parameters based on the amplifying coefficients.

26. The linearization device of claim 25, wherein the linearization device is configured to determine the amplifying coefficients based on N samples of the sums of transmission signals received by the transmit observation receivers via L receive antenna ports by: determining, for each of the plurality of non-linear amplifiers, a Kronecker product between a column vector having as elements the corresponding estimated channel coefficients and the regression matrix; calculating a generalized pseudo-inverse of a matrix formed by concatenating the determined Kronecker products of the plurality of non-linear amplifiers; and determining a matrix product between the calculated generalized pseudo-inverse and a column vector having as elements the N samples of the sums of transmission signals from each of the L receive antenna ports.

27. The linearization device of claim 21, wherein the plurality of communication paths are radio communication paths between the plurality of non-linear amplifiers and the transmit observation receivers via a corresponding plurality of transmit antenna elements and at least two receive antenna elements.

28. An apparatus for a wireless transmitter node, the wireless transmitter node comprising at least two receive antenna ports, each connectable to a corresponding receive antenna element, and a plurality of transmit antenna ports, each connectable to a corresponding transmit antenna element of an active antenna array having a plurality of transmit antenna elements, each transmit antenna port being associated with a respective transmitter branch, wherein each transmitter branch comprises: a non-linear amplifier associated with one of the plurality of transmit antenna elements and with a non-linear transfer function defining an output of the non-linear amplifier based on an input of the non-linear amplifier and based on a reflection signal for the non-linear amplifier, resulting from mutual couplings among the plurality of transmit antenna elements; and digital pre-distortion circuitry configured to compensate the non-linear transfer function by pre-distorting the signal of the transmitter branch based on pre-distortion parameters; the apparatus comprising: a channel estimator configured to estimate a plurality of channel coefficients indicative of channel characteristics of a plurality of communication paths between the plurality of non-linear amplifiers and two or more transmit observation receivers; the two or more transmit observation receivers, each associated with one of the at least two receive antenna ports and each configured to receive a sum of transmission signals generated by the plurality of non-linear amplifiers and transferred over the communication paths between the plurality of non-linear amplifiers and the transmit observation receiver; and the linearization device of claim 21, wherein the first port is connected to the channel estimator and the second port is connected to the transmit observation receiver and wherein the linearization device is configured to provide the determined pre-distortion parameters to the digital pre-distortion circuitry of the transmitter branches.

29. The apparatus of claim 28 further comprising the at least two receive antenna ports and the at least two receive antenna elements.

30. A wireless transmitter node comprising the apparatus of claim 28 and further comprising a linearization device, the linearization device comprising: a first port configured to receive a plurality of channel coefficients indicative of channel characteristics of a plurality of communication paths between the plurality of non-linear amplifiers and two or more transmit observation receivers, wherein each transmit observation receiver is configured to receive a sum of transmission signals generated by the plurality of non-linear amplifiers and transferred over the communication paths between the plurality of non-linear amplifiers and the transmit observation receiver; a second port configured to receive the sums of transmission signals from the transmit observation receivers; and determination circuitry configured to determine the pre-distortion parameters based on the received plurality of channel coefficients, the received sums of transmission signals, and a model of the non-linear transfer functions of the non-linear amplifiers.

31. A cloud based server node comprising the linearization device of claim 21, wherein the cloud based server node is configured to provide the determined pre-distortion parameters to a wireless transmitter node.

32. A method for a linearization device for determining pre-distortion parameters associated with a plurality of non-linear amplifiers, each associated with one of a plurality of transmit antenna elements and with a non-linear transfer function defining an output of the non-linear amplifier based on an input of the non-linear amplifier and based on a reflection signal for the non-linear amplifier, resulting from mutual couplings among the plurality of transmit antenna elements, the method comprising: receiving, via a first port of the linearization device, a plurality of channel coefficients indicative of channel characteristics of a plurality of communication paths between the plurality of non-linear amplifiers and two or more transmit observation receivers, wherein each transmit observation receiver is configured to receive a sum of transmission signals generated by the plurality of non-linear amplifiers and transferred over the communication paths between the plurality of non-linear amplifiers and the transmit observation receiver; receiving, via a second port of the linearization device, the sums of transmission signals from the transmit observation receiver; and determining the pre-distortion parameters based on the received plurality of channel coefficients, the received sums of transmission signals, and a model of the non-linear transfer functions of the non-linear amplifiers.

33. The method of claim 32, wherein the model is represented, for each of the non-linear amplifiers, by a plurality of amplifying coefficients in a space spanned by a regression matrix of the inputs and the reflection signals.

34. The method of claim 33, wherein each reflection signal is modelled as a linear function of the outputs of the plurality of non-linear amplifiers, the linear function being defined by reflection coefficients.

35. The method of claim 34, further comprising determining the reflection coefficients by: defining initial reflection coefficients; determining intermediate amplifying coefficients based on the initial reflection coefficients; determining intermediate reflection coefficients based on the intermediate amplifying coefficients; and recursively determining refined amplifying coefficients based on previous reflection coefficients, and refined reflection coefficients based on previous amplifying coefficients.

36. The method of claim 33, wherein the pre-distortion parameters are determined based on the amplifying coefficients.

37. The method of claim 36, wherein determining the amplifying coefficients based on N samples of the sums of transmission signals received by the transmit observation receiver via L receive antenna ports comprises: determining, for each of the plurality of non-linear amplifiers, a Kronecker product between a column vector having as elements the corresponding estimated channel coefficients and the regression matrix; calculating a generalized pseudo-inverse of a matrix formed by concatenating the determined Kronecker products of the plurality of non-linear amplifiers; and determining a matrix product between the calculated generalized pseudo-inverse and a column vector having as elements the N samples of the sums of transmission signals from each of the L receive antenna ports.

38. The method of claim 32, wherein the plurality of communication paths are radio communication paths between the plurality of non-linear amplifiers and the transmit observation receiver via a corresponding plurality of transmit antenna elements and at least two receive antenna elements.

39. The method of claim 32 further comprising estimating the plurality of channel coefficients.

40. A computer program product comprising a non-transitory computer readable medium, having thereon a computer program comprising program instructions, the computer program being loadable into a data processing unit and configured to cause execution of the method according to claim 32 when the computer program is run by the data processing unit.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0046] Further objects, features and advantages will appear from the following detailed description of embodiments, with reference being made to the accompanying drawings. The drawings are not necessarily to scale, emphasis instead being placed upon illustrating the example embodiments.

[0047] FIG. 1 is a schematic block diagram illustrating an example arrangement according to the prior art;

[0048] FIG. 2 is a schematic block diagram illustrating an example arrangement according to the prior art;

[0049] FIG. 3 is a schematic block diagram illustrating an example arrangement according to some embodiments;

[0050] FIG. 4 is a flowchart illustrating example method steps according to some embodiments;

[0051] and

[0052] FIG. 5 is a schematic drawing illustrating an example computer readable medium according to some embodiments.

DETAILED DESCRIPTION

[0053] As already mentioned above, it should be emphasized that the term “comprises/comprising” when used in this specification is taken to specify the presence of stated features, integers, steps, or components, but does not preclude the presence or addition of one or more other features, integers, steps, components, or groups thereof. As used herein, the singular forms “a”, “an” and “the” are intended to include the plural forms as well, unless the context clearly indicates otherwise.

[0054] Embodiments of the present disclosure will be described and exemplified more fully hereinafter with reference to the accompanying drawings. The solutions disclosed herein can, however, be realized in many different forms and should not be construed as being limited to the embodiments set forth herein.

[0055] In the following, embodiments will be described where two or more TOR:s can be used in connection to determination of pre-distortion parameters for a plurality of transmitter branches without the need for directional couplers or switching networks.

[0056] Embodiment presented herein also provides for handling of the issue encountered for conventional pre-distortion (DPD) approaches; the effects emanating from mutual coupling between antennas. To this end, some embodiments accommodate a need for isolation between transmit branches to cope with the effects of mutual coupling. It may be noted that existing solutions for reduced analog hardware complexity in pre-distortion (DPD) approaches take mutual coupling between antennas into consideration.

[0057] This is achieved by (simultaneous) observation of at least two sums of transmission signals, wherein each transmission signal is generated by a respective one of a plurality of non-linear amplifiers and is transferred over the communication path between the non-linear amplifier and the corresponding transmit observation receiver. The channel coefficients indicative of channel characteristics of the plurality of communication paths between the plurality of non-linear amplifiers and the transmit observation receivers may be estimated by any suitable method or may be otherwise known.

[0058] The pre-distortion parameters are then determined for the plurality of transmitter branches based on the channel coefficients, the received sums of transmission signals, and a model of non-linear transfer functions of the non-linear amplifiers.

[0059] Some embodiments thus provide simplified pre-distortion for antenna arrays with mutual coupling. Embodiments presented herein uses few—typically much fewer than the number of transmission branches—TOR:s (each connected to an observation receiver) to simultaneously identify the pre-distortion parameters for a larger number of different transmit branches, thereby saving cost and/or implementation complexity. Thus, according to some embodiments, the amount of measurement receivers needed to compensate for power amplifier distortion in a large array is reduced compared to prior art approaches.

[0060] Furthermore, the mutual coupling between transmit antenna elements is taken into account in the pre-distortion parameter determination as exemplified herein.

[0061] FIG. 3 schematically illustrates an example arrangement 300 for a wireless transmitter node according to some embodiments. The wireless transmitter node comprises at least two receive antenna ports 355, 356, 357, each connectable to a corresponding receive antenna element 350, 351, 352, and a plurality (exemplified as four in FIG. 3 for illustration) of transmit antenna ports 345, 346, 347, 348, each connectable to a corresponding transmit antenna element of an active antenna array having a plurality of transmit antenna elements 341, 342, 343, 344.

[0062] Each transmit antenna port is associated with a respective transmitter branch, wherein each transmitter branch comprises a non-linear power amplifier (PA) 331, 332, 333, 334 associated with a non-linear transfer function and configured to amplify a signal of the transmitter branch.

[0063] The non-linear transfer function defines an output of the non-linear amplifier based on an input of the non-linear amplifier and based on a reflection signal for the non-linear amplifier, resulting from mutual couplings among the plurality of transmit antenna elements.

[0064] Each transmitter branch also comprises digital pre-distortion circuitry (DPD) 321, 322, 323, 324 configured to compensate the non-linear transfer function by pre-distorting the signal of the transmitter branch based on pre-distortion parameters.

[0065] The pre-distortion circuitry may, for example, comprise respective pre-distortion circuitry for each of the transmitter branches as illustrated in FIG. 3 or common pre-distortion circuitry shared for all of the transmitter branches (see e.g. 320 of FIG. 3). Generally, the pre-distortion circuitry may comprise one or more actuators.

[0066] In one typical example, a multi-antenna transmitter DPD (a common pre-distortion circuitry) 320 may comprise one dual-input DPD per transmission branch (the dual inputs being a signal input and a crosstalk input), and one shared crosstalk mismatch model (CTMM) 325 taking the signal input of each DPD as inputs and providing the crosstalk inputs (estimates of the reflected signals) to each DPD. The CTMM may include the coupling (reflection) coefficients for all M transmit paths, the m.sup.th dual-input DPD may include DPD coefficients (parameters) for the corresponding PA, and all coefficients are identified (and updated if there is an adaptive algorithm used) by the linearization device 380.

[0067] A signal to be transmitted is divided into respective signals to be transmitted by each transmitter branch by a transmit pre-coder (TPC) 310 according to any suitable approach. It may be preferable that the signal to be transmitted comprises a wide range of signal properties (e.g. different amplitudes and/or different frequencies), since such diverse signal properties provides for proper linearization of the non-linear amplifiers.

[0068] Typically, but not necessarily, the arrangement 300 is implemented in a different hardware unit than the transmitter branches.

[0069] The arrangement 300 comprises at least two transmit observation receivers (TOR) 370, 371, 372, each associated with one of the receive antenna ports 355, 356, 357 and configured to receive a sum of transmission signals generated by the plurality of non-linear amplifiers and transferred over a plurality of communication paths 391, 392, 393, 394 (e.g. radio communication paths) between the plurality of non-linear amplifiers and the transmit observation receiver. In FIG. 3, only the communication paths 391, 392, 393, 394 between the plurality of non-linear amplifiers and the transmit observation receiver 370 have been labeled with reference numbers. The communication paths between the plurality of non-linear amplifiers and the other transmit observation receivers 371, 372 are illustrated ad dotted lines without reference number labels.

[0070] If the reflection coefficients (thus, at least implicitly, the mutual couplings among the plurality of transmit antenna elements is known) are known, already estimated, or otherwise in no need for being estimated, the arrangement comprises at least two (for example exactly two) transmit observation receivers (TOR).

[0071] If the reflection coefficients (and implicitly the mutual couplings among the plurality of transmit antenna elements) are not known, not estimated, or otherwise in need for being estimated, the arrangement comprises at least three (for example exactly three) transmit observation receivers (TOR).

[0072] Thus, the wireless transmitter node may comprise two or more receive antenna ports (e.g. as many receive antenna ports as it comprises transmit antenna ports). Typically, a convergence time of the linearization process to be described herein decreases as the number of receive antenna ports increases.

[0073] The arrangement 300 also comprises a channel estimator (CE) 360 configured to estimate a plurality of channel coefficients indicative of channel characteristics of the plurality of communication paths.

[0074] In various embodiments, the arrangement may also comprise some or all of: the at least two receive antenna ports 355, 356, 357, the at least one receive antenna element 350, 351, 352, the plurality of transmit antenna ports 345, 345, 347, 348, the active antenna array comprising the transmit antenna elements 341, 342, 343, 344, and the transmitter branches.

[0075] The arrangement 300 also comprises a linearization device (LIN; e.g. linearization circuitry) 380. The linearization device 380 is configured to determine the pre-distortion parameters associated with the non-linear amplifiers 331, 332, 333, 334 to be used by the digital pre-distortion circuitry 321, 322, 323, 324. To this end, the linearization device comprises a first port 381, a second port 382 and determination circuitry (DET) 383.

[0076] Generally, when a linearization device is referred to herein, a linearization device may be defined as a device which is configured to compensate for a distortion in a signal introduced by one or more non-linear components in a radio system. A non-linear component is a component which has an output that is not a linear function of an input of the component. An example of a non-linear component is a non-linear power amplifier, such as the power amplifiers discussed in connection to FIG. 3. As mentioned above, each of the non-linear power amplifiers is associated with a non-linear transfer function defining an output of the non-linear amplifier based on an input of the non-linear amplifier and based on a reflection signal for the non-linear amplifier, resulting from mutual couplings among the plurality of transmit antenna elements. Thus, the mutual coupling among the plurality of transmit antenna elements is considered generally herein as subject to the linearization.

[0077] The first port 381 is configured to receive a plurality of channel coefficients indicative of channel characteristics of the plurality of communication paths 391, 392, 393, 394 between the plurality of non-linear amplifiers and the transmit observation receivers 370, 371, 372. In the embodiment of FIG. 3, the channel coefficients are estimated by and received from the channel estimator 360. It should be noted that the channel coefficients may be determined or estimated in any suitable manner. In some embodiments, the channel coefficients may even be known beforehand.

[0078] Generally, when channel coefficients are referred to herein, the channel coefficients may be defined according to any suitable, known or future, approach. For example, each channel coefficient may refer to a complex value which describes a difference (in terms of amplitude and phase) between a transmitted signal and a part of the corresponding signal received at a particular time, wherein different channel coefficients refer to different particular times to, e.g. to describe a multipath channel model which in this case is a model of the plurality of communication paths 391, 392, 393, 394 between the plurality of non-linear amplifiers and the transmit observation receivers 370, 371, 372. Thus, the channel estimation referred to herein is not for access purposes (as is typically the case for channel estimation), but is an estimation of the channel comprising the plurality of communication paths 391, 392, 393, 394 between the plurality of non-linear amplifiers and the transmit observation receiver 370.

[0079] Also generally, the communication paths may be wireless (e.g. radio) communication paths and the linearization may be based on over-the-air measurements.

[0080] The second port 382 linearization device is configured to receive the sums of transmission signals from the transmit observation receivers 370, 371, 372, and the determination circuitry 383 is configured to determine the pre-distortion parameters based on the received plurality of channel coefficients, the received sums of transmission signals, and a model of the non-linear transfer functions of the non-linear amplifiers.

[0081] In the embodiment of FIG. 3, the first port 381 is connected to the channel estimator and the second port 382 is connected to the transmit observation receivers. The linearization device is configured to provide the determined pre-distortion parameters, and the determined reflection coefficients (if applicable), to the digital pre-distortion circuitry of the transmitter branches as illustrated by 384 in FIG. 3.

[0082] In FIG. 3, the linearization device 380 (and the arrangement 300) has been illustrated as comprised in a wireless transmitter node. It should be noted that, in other embodiments, the linearization device 380 may be comprised in another node, for example a cloud based server node. In such embodiments, the other node may be configured to provide the determined pre-distortion parameters to a wireless transmitter node for application therein.

[0083] Alternatively or additionally, the channel estimator 360 may be comprised in another node than the wireless transmitter node (which other node may be the same or different than the node comprising the linearization device). In such embodiments, the node comprising the channel estimator may be configured to provide the channel coefficients to the node comprising the linearization device for application therein. In fact, the channel coefficients may not even be provided by a channel estimator according to some embodiments, but may be acquired in some other way (e.g. known beforehand).

[0084] The determination, by the linearization device, of the pre-distortion parameters will be exemplified in the following. In this context it should be noted that some embodiments provide an approach wherein it is possible to determine the pre-distortion parameters from the sum of transmission signals; i.e. it may not be necessary, or even possible, to separate the individual transmission signals of each respective transmitter branch.

[0085] In the following example of the determination of the pre-distortion parameters, it is assumed that there are M transmitter branches, each with a respective non-linear amplifier, the index m identifying the relevant transmitter branch and amplifier, and that there are L receiver ports, each with a respective TOR, the index l identifying the relevant receiver port.

[0086] The output from the m.sup.th amplifier may be modeled using a non-linear transfer function which describes the input-output relation of the non-linear amplifier. The non-linear transfer function for the m.sup.th amplifier defines an output b.sub.2m, 302, of the non-linear amplifier based on an input a.sub.1m, 301, of the non-linear amplifier and based on a reflection signal a.sub.2m, 303, for the non-linear amplifier (the reflected wave injected via mutual coupling), wherein the reflection signal results from mutual couplings among the plurality of transmit antenna elements. The non-linear transfer function for the m.sup.th amplifier may be written as: b.sub.2m=ƒ.sub.m(a.sub.1m, a.sub.2m).

[0087] The above expression of the non-linear transfer function for the m.sup.th amplifier may be written in linear matrix-vector form as: b.sub.2m=G(a.sub.1m, a.sub.2m)θ.sub.m, where θ.sub.m is a vector of amplifying coefficients θ.sub.jm in a (possibly polynomial) model for the m.sup.th amplifier, and G(a.sub.1m, a.sub.2m) is the regression matrix of the m.sup.th amplifier (in many cases the regression matrix may be considered to be equal for all amplifiers).

[0088] The regression matrix represents a model of the non-linear amplifier. The model is represented by a plurality of amplifying coefficients θ.sub.jm in a space spanned by the regression matrix of the inputs and the reflection signals.

[0089] The reflection signal a.sub.2m may be expressed as a function of the M different outputs from the power amplifiers: a.sub.2m=h.sub.m(b.sub.21, . . . , b.sub.2M), where h.sub.m is a function describing the mutual coupling.

[0090] The reflection signal for the m.sup.th amplifier may be expressed as a scalar product: a.sub.2m=b.sub.2.sup.Tλ.sub.m, where b.sub.2=[b.sub.21, b.sub.2M].sup.T and Δ.sub.m=[λ.sub.m1, . . . , λ.sub.mM].sup.T is a vector of reflection coefficients modeling the mutual coupling between the m.sup.th transmit antenna element and the other transmit antenna elements. Thus, each reflection signal may be modelled as a linear function (e.g. linear combination) of the outputs of the plurality of non-linear amplifiers, the linear function being defined by reflection coefficients.

[0091] The received signal at the l.sup.th receiver port is measured as the sum: r.sub.l=Σ.sub.m=1.sup.Mη.sub.lmb.sub.2m, where η.sub.lm is the channel coefficient describing the wireless channel, 391, 392, 393, 304 between the m.sup.th transmitter and the l.sup.th TOR. Given the assumptions above, the received signal may be modelled as: r.sub.l=Σ.sub.m=1.sup.Mη.sub.lmG(a.sub.1m, a.sub.2m)θ.sub.m.

[0092] Using a generalized system model (compare with FIG. 3) with L measurement receivers (i.e. L receiver antenna ports) and M transmitters, the measurement vector may be written as:

[00001]$\left[\begin{array}{c}{r}_1\\ \mathrm{.Math.}\\ r_L\end{array}\right]=(\begin{array}{ccc}[{\eta }_1.Math.G\left(a_{11},a_{21}\right)& \mathrm{.Math.}& {\eta }_M.Math.G\left(a_{1M},a_{2M}\right)])\end{array}\left[\begin{array}{c}{\theta }_1\\ \mathrm{.Math.}\\ {\theta }_M\end{array}\right]=\hspace{1em}\left[\begin{array}{ccc}{\eta }_{11}G\left(a_{11},a_{21}\right)& \mathrm{.Math.}& {\eta }_{1M}G\left(a_{1M},a_{2M}\right)\\ \mathrm{.Math.}& \ddots & \mathrm{.Math.}\\ {\eta }_{\mathrm{L1}}G\left(a_{11},a_{21}\right)& \mathrm{.Math.}& {\eta }_{\mathrm{LM}}G\left(a_{1M},a_{2M}\right)\end{array}\right]\left[\begin{array}{c}{\theta }_1\\ \mathrm{.Math.}\\ {\theta }_M\end{array}\right]$

[0093] where r.sub.l is the vector of N samples received at the l.sup.th receive antenna port, η represents the channel coefficients, and is the Kronecker-product.

[0094] Thus, [η.sub.1G(a.sub.11, a.sub.21) . . . η.sub.MG(a.sub.1M, a.sub.2M)] represents a concatenation of the Kronecker products (of the plurality of non-linear amplifiers) between a column vector having as elements the corresponding estimated channel coefficients and a regression matrix, and r represents N samples of the sum of transmission signals received by the transmit observation receiver via L receive antenna ports.

[0095] Provided that the reflection coefficients modeling the mutual coupling are known, procedure may be used to in identify the amplifier coefficients as:

[00002]$={\left[\begin{array}{ccc}{\eta }_{11}G\left(a_{11},a_{21}\right)& \mathrm{.Math.}& {\eta }_{1M}G\left(a_{1M},a_{2M}\right)\\ \mathrm{.Math.}& \ddots & \mathrm{.Math.}\\ {\eta }_{\mathrm{L1}}G\left(a_{11},a_{21}\right)& \mathrm{.Math.}& {\eta }_{\mathrm{LM}}G\left(a_{1M},a_{2M}\right)\end{array}\right]}^{†}\left[\begin{array}{c}{r}_1\\ \mathrm{.Math.}\\ r_L\end{array}\right]$

[0096] wherein [⋅].sup.† denotes the generalized pseudo-invers, and denotes the least-square estimate of the coefficients.

[0097] The estimated amplifier coefficients thus provided may be used in a next step of an iterative procedure further refining the estimate of the amplifier coefficients.

[0098] An example of such an iterative procedure for estimating the amplifier coefficients can be illustrated by a pseudo-code function using known reflection coefficients and two over-the-air measurements (i.e. two transmit observation receivers, L=2) as follows, where NMSE denotes the normalized mean square error, and {tilde over (B)}.sub.2M.sup.(0) is an augmented N×M matrix consisting of all N×1-dimensional b.sub.2m vectors for all M transmit branches:

TABLE-US-00001 known: input signals a.sub.1m reflection coefficients λ.sub.m channel coefficients η.sub.1m, η.sub.2m measured: received signals r.sub.1, r.sub.2 NMSE.sub.des = X, [define desired accuracy] NMSE = ∞, {tilde over (B)}.sub.2M.sup.(0) = 0, [set initial value of PA model] i = 0, while NMSE ≥ NMSE.sub.des do [iterate to desired accuracy] i = i + 1, For all M compute: ã.sub.2m.sup.(i) = {tilde over (B)}.sub.2M.sup.(i-1)λ.sub.m, [00003]$\left[\begin{array}{c}{\tilde{\theta }}_1^{\left(i\right)}\\ \mathrm{.Math.}\\ {\tilde{\theta }}_M^{\left(i\right)}\end{array}\right]={\left[\begin{array}{ccc}{\eta }_{11}G\left(a_{11},{\tilde{a}}_{21}^{\left(i\right)}\right)& \mathrm{.Math.}& {\eta }_{1M}G\left(a_{1M},{\tilde{a}}_{2M}^{\left(i\right)}\right)\\ {\eta }_{21}G\left(a_{11},{\tilde{a}}_{21}^{\left(i\right)}\right)& \mathrm{.Math.}& {\eta }_{2M}G\left(a_{1M},{\tilde{a}}_{2M}^{\left(i\right)}\right)\end{array}\right]}^{†}\left[\begin{array}{c}{r}_1\\ r_2\end{array}\right],$ {tilde over (B)}.sub.2M.sup.(i) = [G(a.sub.11, ã.sub.21.sup.(i)){tilde over (θ)}.sub.1.sup.(i) . . . G(a.sub.1M, ã.sub.2M.sup.(i)){tilde over (θ)}.sub.M.sup.(i)], [{tilde over (r)}.sub.1.sup.(i) {tilde over (r)}.sub.2.sup.(i)] = {tilde over (B)}.sub.2M.sup.(i)[η.sub.1 η.sub.2], NMSE = max.sub.l=1,2 {NMSE(r.sub.l, {tilde over (r)}.sub.l.sup.(i))}, end while θ.sub.m = {tilde over (θ)}.sub.m.sup.(i).

[0099] If the reflection coefficients modeling the mutual coupling (e.g. due to cross talk and/or mismatch) are not known, they may be determined (i.e. estimated) using a recursive approach in which initial reflection coefficients are defined as an initial assumption. Such initial reflection coefficients may be selected using any suitable approach—For example, it may be initially assumed that there is no mutual coupling between the transmit antenna elements. Then, intermediate amplifying coefficients may be determined based on the initial reflection coefficients and intermediate reflection coefficients may be determined based on the intermediate amplifying coefficients, and so on until some stopping criterion is met (e.g. that the mean square error between subsequent recursions falls below a threshold value and/or that a maximum number of allowable recursions have been carried out).

[0100] An example of such a recursive approach for estimating the amplifier coefficients and the reflection coefficients can be illustrated by a pseudo-code function using three over-the-air measurements (i.e. three transmit observation receivers, L=3) as follows, where r is a 3N×1-dimensional vector, [⋅]* denotes the complex conjugate, P is the dimensionality of the polynomial model (indexed by p), and ƒ.sub.m.sup.(⋅) is a vector of N samples from all the corresponding ƒ.sub.m.sup.(⋅):

TABLE-US-00002 known: input signals a.sub.1m channel coefficients η.sub.1m, η.sub.2m, η.sub.3m measured: received signals r = [r.sub.1.sup.T, r.sub.2.sup.T, r.sub.3.sup.T].sup.T NMSE.sub.des = X,  [define desired accuracy] NMSE = ∞, {tilde over (λ)}.sub.m.sup.(0) = 1, ã.sub.2m.sup.(0) = 0  [set initial values] i = 0, while NMSE ≥ NMSE.sub.des do  [iterate to desired accuracy] i = i + 1, STEP 1:  [find PA model coefficients] [00004]$\left[\begin{array}{c}{\tilde{\theta }}_1^{\left(i\right)}\\ \mathrm{.Math.}\\ {\tilde{\theta }}_M^{\left(i\right)}\end{array}\right]={\left[\begin{array}{ccc}{\eta }_{11}G\left(a_{11},{\tilde{a}}_{21}^{\left(i­1\right)}\right)& \mathrm{.Math.}& {\eta }_{1M}G\left(a_{1M},{\tilde{a}}_{2M}^{\left(i­1\right)}\right)\\ {\eta }_{21}G\left(a_{11},{\tilde{a}}_{21}^{\left(i­1\right)}\right)& \mathrm{.Math.}& {\eta }_{2M}G\left(a_{1M},{\tilde{a}}_{2M}^{\left(i­1\right)}\right)\\ {\eta }_{31}G\left(a_{11},{\tilde{a}}_{21}^{\left(i­1\right)}\right)& \mathrm{.Math.}& {\eta }_{3M}G\left(a_{1M},{\tilde{a}}_{2M}^{\left(i­1\right)}\right)\end{array}\right]}^{†}\left[\begin{array}{c}{r}_1\\ r_2\\ r_3\end{array}\right],$ if i > 1 then STEP 2:  [find reflection coefficients]  f.sub.m.sup.(0) = Σ.sub.p=0.sup.(P-1)/2 {tilde over (α)}.sub.mp.sup.(i)a.sub.1m.sup.p+1a.sub.1m.sup.*p,  f.sub.m.sup.(1) = Σ.sub.p=0.sup.(P-1)/2 {tilde over (β)}.sub.mp.sup.(i)a.sub.1m.sup.pa.sub.1m.sup.*p,  f.sub.m.sup.(2) = Σ.sub.p=0.sup.(P-1)/2 {tilde over (γ)}.sub.mp.sup.(i)a.sub.1m.sup.p+1a.sub.1m.sup.*p−1,  F.sub.m.sup.(1) = diag{f.sub.m.sup.(1)}{tilde over (β)}.sub.2M.sup.(i-1),  F.sub.m.sup.(2) = diag{f.sub.m.sup.(2)}{tilde over (β)}.sub.2M.sup.(i-1)*,  [00005]$f^{\left(0\right)}=\left[\begin{array}{c}{\eta }_{11}{f}_1^{\left(0\right)}+\mathrm{.Math.}+{\eta }_{1M}{f}_M^{\left(0\right)}\\ {\eta }_{21}{f}_1^{\left(0\right)}+\mathrm{.Math.}+{\eta }_{2M}{f}_M^{\left(0\right)}\\ {\eta }_{31}{f}_1^{\left(0\right)}+\mathrm{.Math.}+{\eta }_{3M}{f}_M^{\left(0\right)}\end{array}\right],$  [00006]$F^{\left(1\right)}=\left[\begin{array}{ccc}{\eta }_{11}{F}_1^{\left(1\right)}& \mathrm{.Math.}& {\eta }_{1M}{F}_M^{\left(1\right)}\\ {\eta }_{21}{F}_1^{\left(1\right)}& \mathrm{.Math.}& {\eta }_{2M}{F}_M^{\left(1\right)}\\ {\eta }_{31}{F}_1^{\left(1\right)}& \mathrm{.Math.}& {\eta }_{3M}{F}_M^{\left(1\right)}\end{array}\right],$  [00007]$F^{\left(2\right)}=\left[\begin{array}{ccc}{\eta }_{11}{F}_1^{\left(2\right)}& \text{.Math.}& {\eta }_{1M}{F}_M^{\left(2\right)}\\ {\eta }_{21}{F}_1^{\left(2\right)}& \text{.Math.}& {\eta }_{2M}{F}_M^{\left(2\right)}\\ {\eta }_{31}{F}_1^{\left(2\right)}& \text{.Math.}& {\eta }_{3M}{F}_M^{\left(2\right)}\end{array}\right],$  [00008]$\left[\begin{array}{c}\mathrm{Re}\left\{{\tilde{\lambda }}^{\left(1\right)}\right\}\\ \mathrm{Im}\left\{{\tilde{\lambda }}^{\left(1\right)}\right\}\end{array}\right]={\left[\begin{array}{cc}\mathrm{Re}\left\{F^{\left(1\right)}+F^{\left(2\right)}\right\}& \mathrm{Im}\left\{-F^{\left(1\right)}+F^{\left(2\right)}\right\}\\ \mathrm{Im}\left\{F^{\left(1\right)}+F^{\left(2\right)}\right\}& \mathrm{Re}\left\{F^{\left(1\right)}-F^{\left(2\right)}\right\}\end{array}\right]}^{†}\left[\begin{array}{c}\mathrm{Re}\left\{r-f^{\left(0\right)}\right\}\\ \mathrm{Im}\left\{r-f^{\left(0\right)}\right\}\end{array}\right],$  [where{tilde over (λ)} =[{tilde over (λ)}.sub.1.sup.T, {tilde over (λ)}.sub.2.sup.T, {tilde over (λ)}.sub.3.sup.T].sup.T]  {tilde over (λ)}.sub.m.sup.(i) = Re{{tilde over (λ)}.sub.m.sup.(i)} + jIm{{tilde over (λ)}.sub.m.sup.(i)},  {tilde over (λ)}.sub.m.sup.(i) = {tilde over (λ)}.sub.m.sup.(i)/max.sub.q=1 . . . M;q≠m {tilde over (λ)}.sub.mq.sup.(i), {tilde over (λ)}.sub.mm.sup.(i) = 1,  [normalization to avoid numerical problems]  ã.sub.2m.sup.(i) = {tilde over (B)}.sub.2m.sup.(i-1){tilde over (λ)}.sub.m.sup.(i),  [calculate new estimate of reflection coefficients] end if {tilde over (B)}.sub.2m.sup.(i) = [G(a.sub.11, ã.sub.21.sup.(i)){tilde over (θ)}.sub.1.sup.(i) . . . G(a.sub.1M, ã.sub.2M.sup.(i)){tilde over (θ)}.sub.M.sup.(i)],  [calculate new estimate of PA output signals] [{tilde over (r)}.sub.1.sup.(i) {tilde over (r)}.sub.2.sup.(i) {tilde over (r)}.sub.3.sup.(i)] = {tilde over (B)}.sub.2M.sup.(i)[η.sub.1 η.sub.2 η.sub.3],  [calculate new estimate of received signals] NMSE = max.sub.l=1,2,3 {NMSE(r.sub.l, {tilde over (r)}.sub.l.sup.(i))},  [evaluate accuracy] end while

[0101] Once the amplifier coefficients and the reflection coefficients are estimated or known, the model of the non-linear amplifier characteristics is complete and the digital pre-distortion parameters may be determined using any suitable approach. Examples of such approaches include the MILA (Model-based Indirect Learning Algorithm, described in P. Landin, A. Mayer, and T. Eriksson, “MILA—A Noise Mitigation Technique for RF Power Amplifier Linearization,” in International Multi-Conference on Systems, Signals & Devices, Conference on Communication & Signal Processing, 2014) and other DLA (Direct Learning Algorithm) based methods (which may involve model inversion).

[0102] Typically, the amplifier coefficients and the reflection coefficients are updated at a similar (or the same) frequency as the pre-distortion parameters. Typically, the process of determining the pre-distortion parameters may be continuously performed in an iterative manner or they may be updated responsive to an triggering event (e.g. detection of that the channel has changed, regular time intervals, detection that the DPD does not fulfil some performance requirement, etc.).

[0103] FIG. 4 illustrates an example method 400 according to some embodiments. The example method may, for example, be performed by a linearization device 380 or an arrangement 300 as described and exemplified in connection to FIG. 3, and any features described in connection with FIG. 3 may be equally applicable to the example method 400 according to various embodiments.

[0104] The example method 400 is for determining pre-distortion parameters associated with a plurality of non-linear amplifiers, each associated with a non-linear transfer function defining an output of the non-linear amplifier based on an input of the non-linear amplifier and based on a reflection signal for the non-linear amplifier, as described above.

[0105] In step 410, a plurality of channel coefficients are received (e.g. via a first port of the linearization device). The plurality of channel coefficients are indicative of channel characteristics of a plurality of communication paths between the plurality of non-linear amplifiers and two or more transmit observation receivers, wherein each transmit observation receiver is configured to receive a sum of transmission signals generated by the plurality of non-linear amplifiers and transferred over the communication paths between the plurality of non-linear amplifiers and the transmit observation receiver.

[0106] In step 420, the sums of transmission signals is received from the transmit observation receiver (e.g. via a second port of the linearization device).

[0107] In step 440, the pre-distortion parameters are determined based on the received plurality of channel coefficients, the received sums of transmission signals, and a model of the non-linear transfer functions of the non-linear amplifiers.

[0108] As described and exemplified above, the method may further comprise estimating the plurality of channel coefficients. Typically, a recursive procedure may be applied where intermediate amplifying coefficients are determined based on initially defined reflection coefficients, intermediate reflection coefficients are determined based on the intermediate amplifying coefficients, refined amplifying coefficients are determined recursively based on previous reflection coefficients, and refined reflection coefficients are determined recursively based on previous amplifying coefficients. Such recursion may, for example, continue a predefined number of times and/or until some stopping criterion is met.

[0109] The described embodiments and their equivalents may be realized in software or hardware or a combination thereof. The embodiments may be performed by general purpose circuitry. Examples of general purpose circuitry include digital signal processors (DSP), central processing units (CPU), co-processor units, field programmable gate arrays (FPGA) and other programmable hardware. Alternatively or additionally, the embodiments may be performed by specialized circuitry, such as application specific integrated circuits (ASIC). The general purpose circuitry and/or the specialized circuitry may, for example, be associated with or comprised in an apparatus such as a wireless transmitter node (e.g. a network node) or a cloud based server node.

[0110] Embodiments may appear within an electronic apparatus (such as a wireless transmitter node or a cloud based server node) comprising arrangements, circuitry, and/or logic according to any of the embodiments described herein. Alternatively or additionally, an electronic apparatus (such as a wireless transmitter node or a cloud based server node) may be configured to perform methods according to any of the embodiments described herein.

[0111] According to some embodiments, a computer program product comprises a computer readable medium such as, for example a universal serial bus (USB) memory, a plug-in card, an embedded drive or a read only memory (ROM). FIG. 5 illustrates an example computer readable medium in the form of a compact disc (CD) ROM 500. The computer readable medium has stored thereon a computer program comprising program instructions. The computer program is loadable into a data processor (PROC) 520, which may, for example, be comprised in a wireless transmitter node or a cloud based server node. When loaded into the data processing unit, the computer program may be stored in a memory (MEM) 530 associated with or comprised in the data-processing unit. According to some embodiments, the computer program may, when loaded into and run by the data processing unit, cause execution of method steps according to, for example, the method illustrated in FIG. 4 or otherwise described herein.

[0112] Generally, all terms used herein are to be interpreted according to their ordinary meaning in the relevant technical field, unless a different meaning is clearly given and/or is implied from the context in which it is used.

[0113] Reference has been made herein to various embodiments. However, a person skilled in the art would recognize numerous variations to the described embodiments that would still fall within the scope of the claims.

[0114] For example, the method embodiments described herein discloses example methods through steps being performed in a certain order. However, it is recognized that these sequences of events may take place in another order without departing from the scope of the claims. Furthermore, some method steps may be performed in parallel even though they have been described as being performed in sequence. Thus, the steps of any methods disclosed herein do not have to be performed in the exact order disclosed, unless a step is explicitly described as following or preceding another step and/or where it is implicit that a step must follow or precede another step.

[0115] In the same manner, it should be noted that in the description of embodiments, the partition of functional blocks into particular units is by no means intended as limiting. Contrarily, these partitions are merely examples. Functional blocks described herein as one unit may be split into two or more units. Furthermore, functional blocks described herein as being implemented as two or more units may be merged into fewer (e.g. a single) unit.

[0116] Any feature of any of the embodiments disclosed herein may be applied to any other embodiment, wherever suitable. Likewise, any advantage of any of the embodiments may apply to any other embodiments, and vice versa.

[0117] Hence, it should be understood that the details of the described embodiments are merely examples brought forward for illustrative purposes, and that all variations that fall within the scope of the claims are intended to be embraced therein.