CN107612856A - A kind of digital pre-distortion processing method and device - Google Patents

A kind of digital pre-distortion processing method and device Download PDF

Info

Publication number
CN107612856A
CN107612856A CN201710934026.7A CN201710934026A CN107612856A CN 107612856 A CN107612856 A CN 107612856A CN 201710934026 A CN201710934026 A CN 201710934026A CN 107612856 A CN107612856 A CN 107612856A
Authority
CN
China
Prior art keywords
matrix
structural matrix
default
structural
3rd
Prior art date
Application number
CN201710934026.7A
Other languages
Chinese (zh)
Other versions
CN107612856B (en
Inventor
吴卓智
刁穗东
Original Assignee
京信通信系统(中国)有限公司
京信通信系统(广州)有限公司
京信通信技术(广州)有限公司
天津京信通信系统有限公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 京信通信系统(中国)有限公司, 京信通信系统(广州)有限公司, 京信通信技术(广州)有限公司, 天津京信通信系统有限公司 filed Critical 京信通信系统(中国)有限公司
Priority to CN201710934026.7A priority Critical patent/CN107612856B/en
Publication of CN107612856A publication Critical patent/CN107612856A/en
Application granted granted Critical
Publication of CN107612856B publication Critical patent/CN107612856B/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/38Synchronous or start-stop systems, e.g. for Baudot code
    • H04L25/40Transmitting circuits; Receiving circuits
    • H04L25/49Transmitting circuits; Receiving circuits using code conversion at the transmitter; using predistortion; using insertion of idle bits for obtaining a desired frequency spectrum; using three or more amplitude levels ; Baseband coding techniques specific to data transmission systems

Abstract

The present invention, which provides a kind of digital pre-distortion processing method and device, methods described, to be included:Obtain input signal, the input signal of power amplifier and the feedback signal of digital pre-distortion DPD module;First structural matrix F is determined according to the feedback signal and default predistortion memory polynomial1;According to the first structural matrix F1Determine the first matrix U1, however, it is determined that first matrix U1Conditional number be more than default first threshold, then to the first structural matrix F1It is adjusted, until according to the first structural matrix F after adjustment1The conditional number of determination is less than default Second Threshold;First matrix U1For the first structural matrix F1Conjugate matrices and the first structural matrix F1Product;According to the first structural matrix F after the adjustment1Corresponding digital pre-distortion coefficient matrices A carries out pre-distortion to the input signal.

Description

A kind of digital pre-distortion processing method and device

Technical field

The application is related to communication technical field, more particularly to a kind of digital pre-distortion processing method and device.

Background technology

With the development of radio communication, the fluctuating of its envelope is increasing with peak-to-average force ratio, and RF system circuit produces substantial amounts of Non-linear distortion, reduce communication quality.The non-linear distortion of signal be mainly derived from power amplifier (Power Amplifier, PA non-linear, electric memory effect and dynamic nonlinear effect especially caused by wide-band radio frequency power amplifier).It is right at present The linearization technique of power amplifier mainly includes, and manufacture linearizes higher power amplifier apparatus, method for suppressing peak to average ratio With digital pre-distortion (Digital Pre-Distortion, DPD) method etc..Wherein, high, the applicable band of DPD methods stability Wide scope is wide, precision is high, realizes that difficulty is low.It is main that a distortion with input signal is produced using predistorter in DPD methods The opposite signal of characteristic, so as to offset distortion component, is corrected to distorted signal, is then sent to the signal after correction Power amplifier is amplified output, so as to obtain undistorted signal output.

But in actual applications, as bandwidth is increasing, the digital pre-distortion model of power amplifier is also increasingly Complexity, in some application scenarios, it may appear that DPD solution has unstable phenomenon, can be caused when serious to power amplifier forever Property damage.

In summary, how to improve the stability of DPD solutions is a urgent problem to be solved.

The content of the invention

A kind of digital pre-distortion processing method of the application offer and device, the memory polynomial model based on DPD, according to The conditional number of the pre-distortion coefficient matrix of DPD memory polynomial, judge the stability of DPD solutions, realize the memory according to DPD The conditional number of polynomial pre-distortion coefficient matrix, determine that conditional number is less than the pre-distortion coefficient matrix of threshold value, reduce DPD's The unstable possibility of solution, ensure that the stability of power amplifier output signal, improves the overall performance of communication equipment.

The embodiment of the present application provides a kind of digital pre-distortion processing method, including:

Obtain input signal, the input signal of power amplifier and the feedback signal of digital pre-distortion DPD module;

First structural matrix F is determined according to the feedback signal and default predistortion memory polynomial1

According to the first structural matrix F1Determine the first matrix U1, however, it is determined that first matrix U1Conditional number be more than Default first threshold, then to the first structural matrix F1It is adjusted, until according to the first construction square after adjustment Battle array F1The conditional number of determination is less than default Second Threshold;First matrix U1For the first structural matrix F1Conjugate torque Battle array and the first structural matrix F1Product;

According to the first structural matrix F after the adjustment1Corresponding digital pre-distortion coefficient matrices A is to the input signal Carry out pre-distortion.

A kind of possible implementation, it is described if it is determined that first matrix U1Conditional number be more than default first threshold Value, then to the first structural matrix F1It is adjusted, until according to the first structural matrix F after adjustment1The condition of determination Number is less than default Second Threshold, including:

Maximum order K and maximal memory depth L corresponding to the default predistortion memory polynomial are adjusted to respectively Default first maximum order K2With default first maximal memory depth L2, wherein the K2<K, L2<L;

According to the K2And L2Obtain the second structural matrix F2, and the second structural matrix F2Corresponding second matrix U2's Conditional number is less than the default Second Threshold, second matrix U2For the second structural matrix F2Conjugate matrices and institute State the second structural matrix F2Product;

Element in matrix Q is sequentially added into the second structural matrix F2The 3rd structural matrix F of middle acquisition3, and calculate every Add the 3rd structural matrix F obtained during an element3Corresponding 3rd matrix U3Conditional number;3rd matrix U3For institute State the 3rd structural matrix F3Conjugate matrices and the 3rd structural matrix F3Product;The matrix Q is by appearing in described first Structural matrix F1And the second structural matrix F is not appeared in2Element composition;

By all 3rd matrix Us3Conditional number be ranked up, obtain minimum conditional number and the minimum condition 3rd structural matrix F corresponding to number3

Judge whether the minimum conditional number is less than the default Second Threshold, if so, then with the minimum bar 3rd structural matrix F corresponding to number of packages3As the second structural matrix F2, and recalculate the 3rd structural matrix F3With the 3rd Structural matrix F3Corresponding 3rd matrix U3Conditional number the step of;If it is not, then by the second structural matrix F2As described One structural matrix F1The matrix obtained after adjustment.

Optionally, the default Second Threshold is smaller 0.3dB than the default first threshold.

A kind of possible implementation, it is described if it is determined that first matrix U1Conditional number be more than default first threshold Value, then to the first structural matrix F1It is adjusted, until according to the first structural matrix F after adjustment1The condition of determination Number is less than default Second Threshold, including:

Maximum order K and maximal memory depth L corresponding to the default predistortion memory polynomial are adjusted to respectively Second default maximum order K3With the second default maximal memory depth L3, wherein the K3<K, the L3<L;

According to the K3And L3Obtain the 4th structural matrix F4, the 4th structural matrix F4For the first structural matrix F1 The matrix obtained after adjustment, and the 4th structural matrix F4Corresponding 4th matrix U4Conditional number be less than described default the Two threshold values;4th matrix U4For the 4th structural matrix F4Conjugate matrices and the 4th structural matrix F4Product.

Optionally, the first structural matrix F1By element F10To element FKLComposition, wherein for either element FklBy with Lower formula determines:

Fkl=y (n-l) | y (n-l) |k-1

Wherein, y (n-l) is the sample sequence of the feedback signal of the power amplifier, and n span is [0, M-1], M is the sampling sum of feedback signal, and K is maximum order, and k is polynomial order, and span, between [1, K], l is memory Depth, and span is [0, L], L is maximal memory depth.

First matrix U1Determined by below equation:

First matrix U1Conditional number C1Determined by below equation:

C1=cond (U1)=| | U1||||U1 -1||。

In specific implementation process, the pre-distortion, including:

According to the first structural matrix F and PA input signal z (n), pre-distortion coefficient matrix A is determined:A= (FHF)-1FHZ;

According to pre-distortion coefficient matrix A and the PA input signal, new PA input signal is generated.

The embodiment of the present application provides a kind of digital predistortion process apparatus, and described device includes:

Data acquisition unit, for obtaining input signal, the input signal of power amplifier of digital pre-distortion DPD module And feedback signal;

DPD coefficient matrix determining units, are used for:

First structural matrix F is determined according to the feedback signal and default predistortion memory polynomial1

According to the first structural matrix F1Determine the first matrix U1, however, it is determined that first matrix U1Conditional number be more than Default first threshold, then to the first structural matrix F1It is adjusted, until according to the first construction square after adjustment Battle array F1The conditional number of determination is less than default Second Threshold;First matrix U1For the first structural matrix F1Conjugate torque Battle array and the product of the first structural matrix F;

DPD processing units, for according to the first structural matrix F after the adjustment1Corresponding digital pre-distortion coefficient square Battle array A carries out pre-distortion to the input signal.

A kind of possible implementation, the DPD coefficient matrixes determining unit, is specifically used for:

Maximum order K and maximal memory depth L corresponding to the default predistortion memory polynomial are adjusted to respectively Default first maximum order K2With default first maximal memory depth L2, wherein the K2<K, L2<L;

According to the K2And L2Obtain the second structural matrix F2, and the second structural matrix F2Corresponding second matrix U2's Conditional number is less than the default Second Threshold;Second matrix U2For the second structural matrix F2Conjugate matrices and institute State the second structural matrix F2Product;

Element in matrix Q is sequentially added into the second structural matrix F2The 3rd structural matrix F of middle acquisition3, and calculate every Individual 3rd structural matrix F3Corresponding 3rd matrix U3Conditional number;3rd matrix U3For the 3rd structural matrix F3's Conjugate matrices and the 3rd structural matrix F3Product;The matrix Q is by appearing in the first structural matrix F1And do not occur In the second structural matrix F2Element composition;

By all 3rd matrix Us3Conditional number be ranked up, obtain minimum conditional number and the minimum condition 3rd structural matrix F corresponding to number3

Judge whether the minimum conditional number is less than the default Second Threshold, if so, then with the minimum bar 3rd structural matrix F corresponding to number of packages3As the second structural matrix F2, and recalculate the 3rd structural matrix F3With the 3rd Structural matrix F3Corresponding 3rd matrix U3Conditional number the step of;If it is not, then by the second structural matrix F2As described One structural matrix F1The matrix obtained after adjustment.

Optionally, the DPD coefficient matrixes determining unit, is specifically used for:The default Second Threshold is more default than described The small 0.3dB of first threshold.

A kind of possible implementation, the DPD coefficient matrixes determining unit, is specifically used for:

The first structural matrix F1By element F10To element FKLComposition, wherein for either element FklBy below equation It is determined that:

Fkl=y (n-l) | y (n-l) |k-1

Wherein, y (n-l) is the sample sequence of the feedback signal of the power amplifier, and n span is [0, M-1], M is the sampling sum of feedback signal, and K is maximum order, and the integer of k value between [1, K], l value is between [0, L] Integer, L is maximal memory depth.

First matrix U1Determined by below equation:

First matrix U1Conditional number C1Determined by below equation:

C1=cond (U1)=| | U1||||U1 -1||。

The DPD processing units, are specifically used for:

According to the first structural matrix F and PA input signal z (n), pre-distortion coefficient matrix A is determined:A= (FHF)-1FHZ;

In specific implementation process, according to pre-distortion coefficient matrix A and the input signal, the input for generating new PA is believed Number.

Using such scheme, due to the conditional number of the pre-distortion coefficient matrix of the memory polynomial according to DPD, DPD is judged The stability of solution, the conditional number of the pre-distortion coefficient matrix of memory polynomial according to DPD is realized, determines that conditional number is less than threshold The pre-distortion coefficient matrix of value, the unstable possibility of DPD solution is reduced, ensure that the stabilization of power amplifier output signal Property, improve the overall performance of communication equipment.

Brief description of the drawings

Fig. 1 is a kind of schematic architectural diagram for digital pre-distortion processing method that the embodiment of the present application provides;

Fig. 2 is a kind of digital pre-distortion processing method schematic flow sheet that the embodiment of the present application provides;

Fig. 3 is a kind of digital predistortion process apparatus structural representation that the embodiment of the present application provides.

Embodiment

The application is described in further detail below in conjunction with accompanying drawing.

Fig. 1 is a kind of structural representation for digital pre-distortion processing method that the embodiment of the present application provides.As shown in Figure 1 In the structural representation of digital pre-distortion method, mainly include DPD processing units 101, PA102, data acquisition unit 103, DPD Coefficient matrix determining unit 104.Input signal generates PA input signal, the input signal of the PA by DPD processing units Amplify generation PA output signal by PA, and generate feedback signal, data acquisition unit 103 is defeated by the PA's of acquisition Enter signal and feedback signal, be sent to DPD coefficient matrixes determining unit 104, DPD coefficient matrixes determining unit 104 is according to described PA input signal and feedback signal determines required DPD coefficient matrixes, and generation DPD coefficient matrixes are sent to DPD processing Unit 101.Fig. 1 is the rough schematic view illustrated, and can also include other equipment in system, not drawn in Fig. 1.

In order to more fully understand the application, the application is illustrated below with reference to accompanying drawing.

With reference to foregoing description, referring to Fig. 2, a kind of digital pre-distortion processing method flow provided for the embodiment of the present application is shown It is intended to.This method comprises the following steps:

Step 201:Obtain input signal, the input signal of power amplifier and the feedback letter of digital pre-distortion DPD module Number;

Step 202:First structural matrix F is determined according to the feedback signal and default predistortion memory polynomial;

Step 203:According to the first structural matrix F1Determine the first matrix U1, however, it is determined that first matrix U1Bar Number of packages is more than default first threshold, then to the first structural matrix F1It is adjusted, until according to described the after adjustment One structural matrix F1The conditional number of determination is less than default Second Threshold;First matrix U1For the first structural matrix F1 Conjugate matrices and the first structural matrix F1Product;

Step 204:According to the first structural matrix F after the adjustment1Corresponding digital pre-distortion coefficient matrices A is to described Input signal carries out pre-distortion.

With reference to Fig. 1, in step 201, the input signal of the PA is by input of the DPD processing units 101 to DPD module Signal obtains after carrying out pre-distortion, and the feedback signal is putting according to the PA102 PA exported input signal What big signal obtained.

Specifically, assume it is current there is input signal to input to DPD processing units 101, then DPD processing units 101, to input Input signal carry out digital pre-distortion processing after obtain PA input signal, and export to PA102.By PA102 to the defeated of PA Enter signal be amplified processing after obtain PA input signal amplified signal, by amplified signal divided by PA102 multiplication factor It can obtain feedback signal.

In step 202, it is assumed that the PA currently got input signal is z (n), n=0,1,2 ... M-1, and feedback signal is Y (n), n=0,1,2 ... M-1, wherein, M is PA input signal and the sampling total length of feedback signal, then with memory polynomial Exemplified by model, below equation is obtained:

Wherein, aklFor pre-distortion coefficients to be determined, z (n) is the sample sequence of PA input signal, and y (n) is feedback signal Signal sample sequence, k value is the integer in [1, K] section, and l value is the integer in [0, L] section.Maximum order K, maximal memory depth L is the pre-distortion parameters in the predistortion memory polynomial, can be set according to actual conditions.

The form of formula [1] matrix of being write as is:

Z=F1A [2]

Wherein:

Z=[z (0), z (1), z (2), z (3) ..., z (M-1)]T [3]

Wherein, the first structural matrix F1For:

F1=[F10,F20,...,FK0,F11,...,F1L,...,FKL] [4]

First structural matrix F1In by element F10To element FKLComposition, wherein for either element FklFor:

Fkl=[Fkl(0),Fkl(1),...,Fkl(M-1)]T [5]

Wherein, Fkl(n) it is the first structural matrix FklIn n coefficient value.

For example, the first structural matrix F1Element F10For:

F10=[F10(0),F10(1),...,F10(M-1)]T

Wherein, F10(n) it is the first structural matrix F10In n coefficient value.

Fkl(n)=y (n-l) | y (n-l) |k-1 [6]

In step 203, first matrix U1Determined by below equation:

Wherein, F1 HFor F1Conjugate matrices.

First matrix U1Conditional number C1Determined by below equation:

C1=cond (U1)=| | U1||||U1 -1|| [8]

Wherein, U1 -1For U1Inverse matrix, | | | | to take norm operator.

If the first matrix F1Conditional number C1It is more than default first threshold, then the default predistortion memory is multinomial Maximum order K corresponding to formula and maximal memory depth L is adjusted to default first maximum order K respectively2With default first most Big memory depth L2, wherein the K2<K, L2<L;

According to the K2And L2Generate the second structural matrix F2;The second structural matrix F2Determination method and described the One structural matrix F1Determination method it is identical, difference is the second structural matrix F2In maximum order K value adjustment For default first maximum order K2, maximal memory depth L value is adjusted to default first maximal memory depth L2, and institute State the second structural matrix F2Corresponding second matrix U2Conditional number be less than the default Second Threshold.

In a kind of possible implementation, the default first threshold is obtained by many experiments, when digital pre- When unusual fluctuations occurs in anamorphic system, PA input signal and predistortion feedback signal, and conditional number corresponding to determination are recorded, in advance If first threshold be according to multiple abnormal data determine conditional number average value.The default Second Threshold can be equal to Or less than the default first threshold.

In step 203, the second structural matrix F2For:

Wherein, the second structural matrix F2Second matrix U determined2Conditional number C2Less than described default Two threshold values.

A kind of possible implementation, to improve the precision of predistortion model, if second matrix U2Conditional number it is small In the default Second Threshold, the conditional number of second matrix after ensureing to adjust still meets to be less than default second threshold Under conditions of value, the second structural matrix F can be updated in the following manner2

Element in matrix Q is sequentially added into the second structural matrix F2The 3rd structural matrix F of middle acquisition3, and calculate every Add the 3rd structural matrix F obtained during an element3Corresponding 3rd matrix U3Conditional number;3rd matrix U3For institute State the 3rd structural matrix F3Conjugate matrices and the 3rd structural matrix F3Product;The matrix Q is by appearing in described first Structural matrix F1And the second structural matrix F is not appeared in2Element composition;

By all 3rd matrix Us3Conditional number be ranked up, obtain minimum conditional number and the minimum condition 3rd structural matrix F corresponding to number3

Judge whether the minimum conditional number is less than the default Second Threshold, if so, then with the minimum bar 3rd structural matrix F corresponding to number of packages3To the second structural matrix F2It is replaced, the second structural matrix F after being replaced2, and The 3rd structural matrix F will be recalculated3With the 3rd structural matrix F3Corresponding 3rd matrix U3Conditional number the step of;If it is not, then By the second structural matrix F2The matrix obtained after being adjusted as first structural matrix.

For example, however, it is determined that second matrix U2Conditional number C2Less than default Second Threshold, then the matrix Q For:

According to the second structural matrix F2With the matrix Q, N number of 3rd structural matrix F is obtained3, specifically, described Three structural matrix F3For:

Wherein, FiFor a coefficient in matrix Q;N is matrix Q number of coefficients.Each 3rd structural matrix F3From matrix The coefficient chosen in Q is different, therefore, the 3rd structural matrix F3There are N kinds may.

Each 3rd structural matrix F3The 3rd matrix U can be determined according to formula [7]3, each 3rd matrix U3Root again Conditional number C is determined according to formula [8]3, N number of conditional number C can be determined altogether3, determine N number of conditional number C3In minimum bar Number of packages.

If the minimal condition number is less than the default Second Threshold, by the second structural matrix F2It is updated to institute State the 3rd structural matrix F corresponding to minimal condition number3, will the second structural matrix F2Corresponding to the minimal condition number 3rd structural matrix F3Replace, obtain and realize to the second structural matrix F2Renewal.The second structural matrix F is updated2 Afterwards, return to and calculate the 3rd structural matrix F3With the 3rd structural matrix F3Corresponding 3rd matrix U3Conditional number the step of, after The second structural matrix F of continuous adjustment2

Accordingly, if the minimal condition number is more than or equal to the default Second Threshold, the institute can be determined State the second structural matrix F2For required structural matrix, so as to according to the second structural matrix F2Corresponding predistortion Coefficient carries out pre-distortion to the input signal.

The first structural matrix F after the adjustment obtained according to this method1In element, after the adjustment first construction Matrix F1Corresponding first matrix meets that the conditional number is less than under conditions of the Second Threshold, first structure after adjustment Make matrix F1The element of corresponding digital pre-distortion coefficient is enough, greatly improves the precision of digital pre-distortion model, and And it ensure that the stability of digital pre-distortion model.

In step 203, a kind of possible implementation, will be maximum corresponding to the default predistortion memory polynomial Exponent number K and maximal memory depth L is adjusted to the second default maximum order K respectively3With the second default maximal memory depth L3, wherein The K3<K, the L3<L;

According to the K3And L3Obtain the 4th structural matrix F4, the 4th structural matrix F4For the first structural matrix F1 The matrix obtained after adjustment, and the 4th structural matrix F4Corresponding 4th matrix U4Conditional number be less than described default the Two threshold values;4th matrix U4For the 4th structural matrix F4Conjugate matrices and the 4th structural matrix F4Product.

According to the first structural matrix F after the adjustment1Corresponding digital pre-distortion coefficient matrices A is to the input signal Carry out pre-distortion.

The embodiment of the present application provides a kind of digital pre-distortion processing method, comprises the following steps:

Step 1:Assuming that it is 4 that K values, which are 7, L values, then the first structural matrix F is determined according to formula [4]1For:

F1=[F10,F20,…,F70,F11,…,F71,…F14,…,F74]

Step 2:If the first structural matrix F1Corresponding first matrix U1Conditional number C1More than default first threshold, then Generate the second structural matrix F2

Assuming that K2For 6, L2For 3, then the second structural matrix F2For:

F2=[F10,F20,...,F60,F11,...,F13,...,F63]

Step 3:Second structural matrix F2Second matrix U determined2Conditional number C2Less than default Second Threshold, Then according to the second structural matrix F2Pre-distortion coefficient matrix is determined with the input signal of the PA, according to pre-distortion coefficients square Battle array carries out pre-distortion to the input signal.

The embodiment of the present application provides a kind of digital pre-distortion processing method, comprises the following steps:

Step 1:Assuming that it is 4 that K values, which are 7, L values, then the first structural matrix F is determined according to formula [4]1For:

F1=[F10,F20,…,F70,F11,…,F71,…F14,…,F74]

Step 2:If the first structural matrix F1Corresponding first matrix U1Conditional number C1More than default first threshold, then Generate the second structural matrix F2

Assuming that K2For 6, L2For 3, then the second structural matrix F2For:

F2=[F10,F20,...,F60,F11,...,F13,...,F63]

Step 3:If it is determined that second matrix U2Conditional number C2Less than default Second Threshold, then matrix Q is:

Q=[F70,F71,F72,F73,F74,F14,F24,F34,F44,F54,F64]

According to the second structural matrix F2With matrix Q, 11 the 3rd structural matrix F are obtained3, specifically, the 3rd structure Make matrix F3For:

F301=[F10,F20,...,F60,F11,...,F13,...,F63,F70]

F311=[F10,F20,...,F60,F11,...,F13,...,F63,F74]

Each 3rd structural matrix F3The 3rd matrix U can be determined according to formula [7]3, each 3rd matrix U3Root again Conditional number C is determined according to formula [8]3, 11 conditional number C can be determined altogether3, determine 11 conditional number C3In minimum Conditional number number C311

Step 4:If the minimal condition number C311Less than the default Second Threshold, then with the 3rd structural matrix F311 To the second structural matrix F2It is replaced, the second structural matrix F after being replaced2, and step 3 is back to, if minimal condition Number C311More than or equal to the default Second Threshold, then step 5 is gone to.Detailed process is as follows:

Now, the second structural matrix F2For:

F2=[F10,F20,...,F60,F11,...,F13,...,F63,F74]

Matrix Q is:

Q=[F70,F71,F72,F73,F14,F24,F34,F44,F54,F64]

According to the second structural matrix F2With the 3rd structural matrix F3, obtain 10 the 4th structural matrix F3, specifically It is as follows:

F301=[F10,F20,...,F60,F11,...,F13,...,F63,F74,F70]

F310=[F10,F20,...,F60,F11,...,F13,...,F63,F74,F64]

Each 3rd structural matrix F3The 3rd matrix U can be determined according to formula [7]3, each 3rd matrix U3Root again Conditional number C is determined according to formula [8]3, 10 conditional number C can be determined altogether3, determine 10 conditional number C3In minimum Conditional number.Assuming that F310Corresponding conditional number is minimal condition number C310

If minimal condition number C310More than or equal to the default Second Threshold, then according to the second structural matrix F2 (now, second structural matrix has been updated to F311), pre-distortion is carried out to the input signal;

If the minimal condition number C310Less than the default Second Threshold, then by the second structural matrix F2Renewal For the 3rd structural matrix F310, and return to step 3.Detailed process is same as described above, will not be repeated here.

Circulation step three and step 4, until the 3rd structural matrix F of generation3Corresponding 3rd matrix U3Minimal condition Number is more than or equal to the default Second Threshold, stops updating second structural matrix.

Step 5:According to the second structural matrix F2Pre-distortion is carried out to the input signal.

Optionally, to avoid frequently carrying out step 203, default Second Threshold could be arranged to be less than default first threshold Value, during practical application, default Second Threshold value can be smaller 0.3dB than default first threshold.

In step 204, pre-distortion coefficient matrix A is:

A=[a10,a20,...,aK0,a11,...,a1L,...,aKL]T[12]

The solution of formula [2] is equivalent to the solution of below equation group:

F1 HF1A=F1 HZ[13]

Wherein, F1 HFor F1Conjugate matrices, work as F1 HF1Can the inverse time, then can obtain:

A=(F1 HF1)-1F1 HZ[14]

Wherein, (F1 HF1)-1For F1 HF1Inverse matrix.

New PA input signal is generated according to pre-distortion coefficient matrix A and the input signal.Formula loses in advance in [14] The specific determination method of true coefficient matrix A, the embodiment of the present application are not limited this, will not be repeated here.

As shown in figure 3, the embodiment of the present application provides a kind of digital predistortion process apparatus, described device includes:

Data acquisition unit 301, the input of input signal, power amplifier for obtaining digital pre-distortion DPD module Signal and feedback signal;

DPD coefficient matrixes determining unit 302, is used for:

First structural matrix F is determined according to the feedback signal and default predistortion memory polynomial1

According to the first structural matrix F1Determine the first matrix U1, however, it is determined that first matrix U1Conditional number be more than Default first threshold, then to the first structural matrix F1It is adjusted, until according to the first construction square after adjustment Battle array F1The conditional number of determination is less than default Second Threshold;First matrix U1For the first structural matrix F1Conjugate torque Battle array and the product of the first structural matrix F;

DPD processing units 303, for according to the first structural matrix F after the adjustment1Corresponding digital pre-distortion coefficient Matrix A carries out pre-distortion to the input signal.

A kind of possible implementation, DPD coefficient matrixes determining unit 302, is specifically used for:

Maximum order K and maximal memory depth L corresponding to the default predistortion memory polynomial are adjusted to respectively Default first maximum order K2With default first maximal memory depth L2, wherein the K2<K, L2<L;

According to the K2And L2Obtain the second structural matrix F2, and the second structural matrix F2Corresponding second matrix U2's Conditional number is less than the default Second Threshold;Second matrix U2For the second structural matrix F2Conjugate matrices and institute State the second structural matrix F2Product;

Element in matrix Q is sequentially added into the second structural matrix F2The 3rd structural matrix F of middle acquisition3, and calculate every Individual 3rd structural matrix F3Corresponding 3rd matrix U3Conditional number;3rd matrix U3For the 3rd structural matrix F3's Conjugate matrices and the 3rd structural matrix F3Product;The matrix Q is by appearing in the first structural matrix F1And do not occur In the second structural matrix F2Element composition;

By all 3rd matrix Us3Conditional number be ranked up, obtain minimum conditional number and the minimum condition 3rd structural matrix F corresponding to number3

Judge whether the minimum conditional number is less than the default Second Threshold, if so, then with the minimum bar 3rd structural matrix F corresponding to number of packages3To the second structural matrix F2It is replaced, the second structural matrix F after being replaced2, and Recalculate the 3rd structural matrix F3With the 3rd structural matrix F3Corresponding 3rd matrix U3Conditional number the step of;If it is not, then will The second structural matrix F2The matrix obtained after being adjusted as first structural matrix.

Optionally, DPD coefficient matrixes determining unit 302, is specifically used for:The default Second Threshold is more default than described The small 0.3dB of first threshold.

A kind of possible implementation, DPD coefficient matrixes determining unit 302, is specifically used for:

The first structural matrix F1By element F10To element FKLComposition, wherein for either element FklBy below equation It is determined that:

Fkl=y (n-l) | y (n-l) |k-1

Wherein, y (n-l) is the sample sequence of the feedback signal of the power amplifier, and n span is [0, M-1], M is the sampling sum of feedback signal, and K is maximum order, and the integer of k value between [1, K], l value is between [0, L] Integer, L is maximal memory depth.

First matrix U1Determined by below equation:

First matrix U1Conditional number C1Determined by below equation:

C1=cond (U1)=| | U1||||U1 -1||。

DPD processing units 303, are specifically used for:

According to the first structural matrix F1With the input signal z (n) of the PA, pre-distortion coefficient matrix A is determined:A= (F1 HF1)-1F1 HZ;

In specific implementation process, according to pre-distortion coefficient matrix A and the input signal, the input for generating new PA is believed Number.

In the embodiment of the present application, by judging whether the conditional number of the first matrix is more than default first threshold, to judge The stability of DPD solutions, if the conditional number of the first matrix is more than default first threshold, generates the second structural matrix.If second The conditional number of second matrix corresponding to structural matrix is more than default Second Threshold, then determines predistortion according to the second structural matrix Coefficient matrix;If the conditional number of the second matrix corresponding to the second structural matrix is less than default Second Threshold, according to the first structure Matrix is made to select the coefficient of the second structural matrix.In prior art, according to the different maximum order of selection and most The value of big memory depth, it is determined that optimal DPD solution.The method that the embodiment of the present application uses is more flexible, avoids repetition meter Pre-distortion coefficient matrix is calculated, greatly reduces amount of calculation, more stable DPD solution can be obtained faster, reduce DPD Solution unstable possibility, ensure that the stability of power amplifier output signal, improve the overall performance of communication equipment.

It should be understood by those skilled in the art that, embodiments herein can be provided as method, system or determine machine program Product.Therefore, the application can use the reality in terms of complete hardware embodiment, complete software embodiment or combination software and hardware Apply the form of example.Moreover, the application can use the determination machine for wherein including determination machine usable program code in one or more The shape for the determination machine program product that usable storage medium is implemented on (including but is not limited to magnetic disk storage, optical memory etc.) Formula.

The application be with reference to according to the present processes, equipment (system) and determine machine program product flow chart and/or Block diagram describes.It should be understood that can by each flow in determination machine programmed instruction implementation process figure and/or block diagram and/or Square frame and the flow in flow chart and/or block diagram and/or the combination of square frame.These determination machine programmed instruction can be provided to arrive General determination machine, special determination machine, the processor of Embedded Processor or other programmable data processing devices are to produce one Machine so that produced by the instruction of determination machine or the computing device of other programmable data processing devices and flowed for realizing The device for the function of being specified in one flow of journey figure or multiple flows and/or one square frame of block diagram or multiple square frames.

These determination machine programmed instruction, which may be alternatively stored in, can guide determination machine or other programmable data processing devices with spy Determine in the determination machine readable memory that mode works so that the instruction being stored in the determination machine readable memory, which produces, to be included referring to Make the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one square frame of block diagram or The function of being specified in multiple square frames.

These determination machine programmed instruction can be also loaded into determination machine or other programmable data processing devices so that true Determine on machine or other programmable devices perform series of operation steps with produce determination machine realization processing, so as to it is determined that machine or The instruction performed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one The step of function of being specified in individual square frame or multiple square frames.

Obviously, those skilled in the art can carry out the essence of various changes and modification without departing from the application to the application God and scope.So, if these modifications and variations of the application belong to the scope of the application claim and its equivalent technologies Within, then the application is also intended to comprising including these changes and modification.

Claims (10)

1. a kind of digital pre-distortion processing method, it is characterised in that methods described includes:
Obtain input signal, power amplifier the PA input signal and feedback signal of digital pre-distortion DPD module;
First structural matrix F is determined according to the feedback signal and default predistortion memory polynomial1
According to the first structural matrix F1Determine the first matrix U1, however, it is determined that first matrix U1Conditional number be more than it is default First threshold, then to the first structural matrix F1It is adjusted, until according to the first structural matrix F after adjustment1Really Fixed conditional number is less than default Second Threshold;First matrix U1For the first structural matrix F1Conjugate matrices and institute State the first structural matrix F1Product;
According to the first structural matrix F after the adjustment1Input of the corresponding digital pre-distortion coefficient matrices A to the DPD module Signal carries out pre-distortion.
2. according to the method for claim 1, it is characterised in that described if it is determined that first matrix U1Conditional number be more than Default first threshold, then to the first structural matrix F1It is adjusted, until according to the first construction square after adjustment Battle array F1The conditional number of determination is less than default Second Threshold, including:
Maximum order K and maximal memory depth L corresponding to the default predistortion memory polynomial are adjusted to default respectively The first maximum order K2With default first maximal memory depth L2, wherein the K2<K, L2<L;
According to the K2And L2Obtain the second structural matrix F2, and the second structural matrix F2Corresponding second matrix U2Condition Number is less than default Second Threshold, second matrix U2For the second structural matrix F2Conjugate matrices and second structure Make matrix F2Product;
Element in matrix Q is sequentially added into the second structural matrix F2The 3rd structural matrix F of middle acquisition3, and calculate and often add The 3rd structural matrix F obtained during one element3Corresponding 3rd matrix U3Conditional number;3rd matrix U3For described Three structural matrix F3Conjugate matrices and the 3rd structural matrix F3Product;The matrix Q is by appearing in first construction Matrix F1And the second structural matrix F is not appeared in2Element composition;
By all 3rd matrix Us3Conditional number be ranked up, obtain minimum conditional number and the minimum conditional number pair The 3rd structural matrix F answered3
Judge whether the minimum conditional number is less than the default Second Threshold, if so, then by the minimum conditional number Corresponding 3rd structural matrix F3As the second structural matrix F2, and recalculate the 3rd structural matrix F3With the 3rd construction Matrix F3Corresponding 3rd matrix U3Conditional number the step of;If it is not, then by the second structural matrix F2As first structure Make matrix F1The matrix obtained after adjustment.
3. method according to claim 1 or 2, it is characterised in that the default Second Threshold is than described default The small 0.3dB of one threshold value.
4. according to the method for claim 1, it is characterised in that described if it is determined that first matrix U1Conditional number be more than Default first threshold, then to the first structural matrix F1It is adjusted, until according to the first construction square after adjustment Battle array F1The conditional number of determination is less than default Second Threshold, including:
Maximum order K and maximal memory depth L corresponding to the default predistortion memory polynomial are adjusted to default respectively The second maximum order K3With default second maximal memory depth L3, wherein the K3<K, the L3<L;
According to the K3And L3Obtain the 4th structural matrix F4, the 4th structural matrix F4For the first structural matrix F1Adjustment The matrix obtained afterwards, and the 4th structural matrix F4Corresponding 4th matrix U4Conditional number be less than default second threshold Value;4th matrix U4For the 4th structural matrix F4Conjugate matrices and the 4th structural matrix F4Product.
5. according to the method for claim 1, it is characterised in that the first structural matrix F1By element F10To element FKLGroup Into wherein for either element FklDetermined by below equation:
Fkl=y (n-l) | y (n-l) |k-1
Wherein, y (n-l) is the sample sequence of the feedback signal of the power amplifier, and n span is [0, M-1], and M is The sampling sum of feedback signal, K are maximum order, and k is polynomial order, and span, between [1, K], l is deep for memory Degree, and span is [0, L], L is maximal memory depth.
6. a kind of digital predistortion process apparatus, it is characterised in that described device includes:
Data acquisition unit, for obtaining the input signal of digital pre-distortion DPD module, the input signal of power amplifier and anti- Feedback signal;
DPD coefficient matrix determining units, are used for:
First structural matrix F is determined according to the feedback signal and default predistortion memory polynomial1
According to the first structural matrix F1Determine the first matrix U1, however, it is determined that first matrix U1Conditional number be more than it is default First threshold, then to the first structural matrix F1It is adjusted, until according to the first structural matrix F after adjustment1Really Fixed conditional number is less than default Second Threshold;First matrix U1For the first structural matrix F1Conjugate matrices and institute State the first structural matrix F product;
DPD processing units, for according to the first structural matrix F after the adjustment1Corresponding digital pre-distortion coefficient matrices A pair The input signal carries out pre-distortion.
7. device according to claim 6, it is characterised in that the DPD coefficient matrixes determining unit, be specifically used for:
Maximum order K and maximal memory depth L corresponding to the default predistortion memory polynomial are adjusted to default respectively The first maximum order K2With default first maximal memory depth L2, wherein the K2<K, L2<L;
According to the K2And L2Obtain the second structural matrix F2, and the second structural matrix F2Corresponding second matrix U2Condition Number is less than the default Second Threshold;Second matrix U2For the second structural matrix F2Conjugate matrices and described the Two structural matrix F2Product;
Element in matrix Q is sequentially added into the second structural matrix F2The 3rd structural matrix F of middle acquisition3, and calculate each Three structural matrix F3Corresponding 3rd matrix U3Conditional number;3rd matrix U3For the 3rd structural matrix F3Conjugation Matrix and the 3rd structural matrix F3Product;The matrix Q is by appearing in the first structural matrix F1And institute is not appeared in State the second structural matrix F2Element composition;
By all 3rd matrix Us3Conditional number be ranked up, obtain minimum conditional number and the minimum conditional number pair The 3rd structural matrix F answered3
Judge whether the minimum conditional number is less than the default Second Threshold, if so, then with the minimum conditional number Corresponding 3rd structural matrix F3As the second structural matrix F2, and recalculate the 3rd structural matrix F3With the 3rd construction Matrix F3Corresponding 3rd matrix U3Conditional number the step of;If it is not, then by the second structural matrix F2As first structure Make matrix F1The matrix obtained after adjustment.
8. device according to claim 6, it is characterised in that the DPD coefficient matrixes determining unit, be specifically used for:
The default Second Threshold is smaller 0.3dB than the default first threshold.
9. device according to claim 6, it is characterised in that the DPD coefficient matrixes determining unit, be specifically used for:
The first structural matrix F1In by element F10To element FKLComposition, wherein for either element FklIt is true by below equation It is fixed:
Fkl=y (n-l) | y (n-l) |k-1
Wherein, y (n-l) is the sample sequence of the feedback signal of the power amplifier, and n span is [0, M-1], M For the sampling sum of feedback signal, K is maximum order, and the integer of k value between [1, K], l value is between [0, L] Integer, L is maximal memory depth.
10. device according to claim 6, it is characterised in that the DPD processing units, be specifically used for:
According to the first structural matrix F and PA input signal z (n), pre-distortion coefficient matrix A is determined:A=(FHF)- 1FHZ;
According to pre-distortion coefficient matrix A and the PA input signal, new PA input signal is generated.
CN201710934026.7A 2017-10-10 2017-10-10 Digital pre-distortion processing method and device CN107612856B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710934026.7A CN107612856B (en) 2017-10-10 2017-10-10 Digital pre-distortion processing method and device

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201710934026.7A CN107612856B (en) 2017-10-10 2017-10-10 Digital pre-distortion processing method and device
PCT/CN2018/096530 WO2019071992A1 (en) 2017-10-10 2018-07-20 Digital pre-distortion processing method and apparatus

Publications (2)

Publication Number Publication Date
CN107612856A true CN107612856A (en) 2018-01-19
CN107612856B CN107612856B (en) 2020-01-14

Family

ID=61068296

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710934026.7A CN107612856B (en) 2017-10-10 2017-10-10 Digital pre-distortion processing method and device

Country Status (2)

Country Link
CN (1) CN107612856B (en)
WO (1) WO2019071992A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019071992A1 (en) * 2017-10-10 2019-04-18 京信通信系统(中国)有限公司 Digital pre-distortion processing method and apparatus

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1547801A (en) * 2002-05-14 2004-11-17 松下电器产业株式会社 Method and device for compensating hybrid distortion
CN101459636A (en) * 2007-12-12 2009-06-17 中兴通讯股份有限公司 Adaptive pre-distortion method
CN102651634A (en) * 2011-02-28 2012-08-29 英特尔移动通信有限公司 Joint adaptive bias point adjustment and digital pre-distortion for power amplifier
CN104168238A (en) * 2013-05-17 2014-11-26 上海无线通信研究中心 Method for lowering self-adaptive digital pre-distortion algorithm computation complexity
US20140355714A1 (en) * 2013-06-04 2014-12-04 Qualcomm Incorporated Configurable pre-emphasis component for transmission circuitry
CN105763503A (en) * 2016-04-25 2016-07-13 京信通信技术(广州)有限公司 Digital pre-distortion processing method and apparatus

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101610093B (en) * 2008-06-17 2012-09-05 富士通株式会社 Predistortion device, predistortion system and predistortion method
AU2010364182B2 (en) * 2010-11-16 2016-07-07 Telefonaktiebolaget L M Ericsson (Publ) Configurable basis-function generation for nonlinear modeling
CN102739586B (en) * 2012-06-14 2015-03-25 京信通信系统(中国)有限公司 Method, equipment and system for adjusting linear performance of predistortion
CN103618684B (en) * 2013-11-25 2017-09-29 京信通信系统(中国)有限公司 A kind of method and device of estimation DPD coefficients
CN103701738B (en) * 2013-12-09 2017-02-15 京信通信系统(中国)有限公司 Pre-distortion parameter determination method, equipment and digital pre-distortion processing system
CN103957179B (en) * 2014-04-04 2018-04-20 京信通信系统(中国)有限公司 DPD method and systems
CN107612856B (en) * 2017-10-10 2020-01-14 京信通信系统(中国)有限公司 Digital pre-distortion processing method and device

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1547801A (en) * 2002-05-14 2004-11-17 松下电器产业株式会社 Method and device for compensating hybrid distortion
CN101459636A (en) * 2007-12-12 2009-06-17 中兴通讯股份有限公司 Adaptive pre-distortion method
CN102651634A (en) * 2011-02-28 2012-08-29 英特尔移动通信有限公司 Joint adaptive bias point adjustment and digital pre-distortion for power amplifier
CN104168238A (en) * 2013-05-17 2014-11-26 上海无线通信研究中心 Method for lowering self-adaptive digital pre-distortion algorithm computation complexity
US20140355714A1 (en) * 2013-06-04 2014-12-04 Qualcomm Incorporated Configurable pre-emphasis component for transmission circuitry
CN105763503A (en) * 2016-04-25 2016-07-13 京信通信技术(广州)有限公司 Digital pre-distortion processing method and apparatus

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019071992A1 (en) * 2017-10-10 2019-04-18 京信通信系统(中国)有限公司 Digital pre-distortion processing method and apparatus

Also Published As

Publication number Publication date
WO2019071992A1 (en) 2019-04-18
CN107612856B (en) 2020-01-14

Similar Documents

Publication Publication Date Title
Bianco et al. High-order central schemes for hyperbolic systems of conservation laws
Ding et al. Memory polynomial predistorter based on the indirect learning architecture
KR100864196B1 (en) Digital predistortion system and method for linearizing an RF power amplifier with nonlinear gain characteristics and memory effects
CN100511974C (en) Self-adapting digital pre-distortion circuti, corresponding method and wireless communication system
JP5231218B2 (en) In-place transform with application to encoding and decoding of various code classes
US20050195919A1 (en) Digital predistortion system and method for high efficiency transmitters
US20090088093A1 (en) Signal predistortion in radio transmitter
Cavers Optimum table spacing in predistorting amplifier linearizers
CN100544203C (en) Nonlinear filter
Litsyn Peak power control in multicarrier communications
JP2008504721A (en) Method and apparatus for equalizing speakers
JP2011527860A (en) Method and system for reducing the PAPR of an OFDM signal
JP4685860B2 (en) Improvements in multi-carrier modulation systems
US9225501B2 (en) Non-linear modeling of a physical system using look-up table with polynomial interpolation
US8831133B2 (en) Recursive digital pre-distortion (DPD)
JP2008295089A (en) Distortion compensating amplifier apparatus
US7091779B2 (en) Non-linear modeling method
US6737918B2 (en) Distortion compensation apparatus
US8271859B2 (en) Segmented CRC design in high speed networks
JP4077449B2 (en) Table reference predistorter
CN1790920A (en) A transmitter with an envelope tracking power amplifier utilizing digital predistortion of the signal envelope
Wiechman et al. An improved sphere-packing bound for finite-length codes over symmetric memoryless channels
US8037399B2 (en) Techniques for segmented CRC design in high speed networks
US7729446B2 (en) System and method for digitally correcting a non-linear element using a multiply partitioned architecture for predistortion
US7822146B2 (en) System and method for digitally correcting a non-linear element

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant