Disclosure of Invention
In view of this, an object of the present application is to provide a method, an apparatus, an electronic device, and a medium for extracting DPD parameters, which can continuously update the DPD parameters according to various data under the condition of a low sampling rate, so that out-of-band distortion of PA signals is smaller and smaller, and finally, after the DPD parameters are updated, the same DPD effect as that under the condition of a high sampling rate is achieved, thereby reducing the limitation of the DPD when the DPD is applied to a wideband.
In a first aspect, an embodiment of the present application provides a DPD parameter extraction method applied to 5G, where the method includes:
inputting the received 5G new air interface 5G NR signal into a digital predistortion model, and performing digital predistortion processing to obtain a predistortion signal;
inputting the predistortion signal into a determined power amplifier model, and performing power amplification treatment to obtain a power amplifier PA signal with a complete frequency spectrum;
sequentially performing band-limited filtering and low-speed analog-to-digital converter (ADC) sampling under the limited frequency on the PA signal to obtain a band-limited frequency spectrum in a truncation interval corresponding to the limited frequency in the complete frequency spectrum;
determining a spectrum extrapolation signal with an approximately complete spectrum based on a mapping operator for performing time domain space projection on the truncation interval, an in-band signal corresponding to the band-limited spectrum, a data matrix formed by a basis function of the power amplifier model and a vector corresponding to the 5G NR signal, and the PA signal;
according to the iterative formula betat+1=βt-μ(XHX)-1XH(y' -x), iterating the Digital Predistortion (DPD) parameters in the digital predistortion model until the out-of-band distortion of the PA signal is smaller than a preset threshold value, and taking the finally obtained numerical value as the DPD parameters of the digital predistortion model; wherein, betatIs the DPD parameter after the t-th iteration, X is the data matrix, XHAnd a conjugate transpose matrix representing the data matrix, wherein x is the 5G NR signal, y' is the spectrum extrapolation signal, and mu is a preset parameter.
In a possible embodiment, determining a spectral extrapolation signal having an approximately complete spectrum based on a mapping operator for performing time-domain spatial projection on the truncated interval, an in-band signal corresponding to the band-limited spectrum, a data matrix formed by basis functions of the power amplifier model and vectors corresponding to the 5G NR signal, and the PA signal includes:
determining a first projection of an in-band signal in the time domain based on the mapping operator and the in-band signal;
determining an extrapolation operator based on the mapping operator, the data matrix, and the PA signal;
determining the spectral extrapolation signal based on the extrapolation operator and the first projection.
In a possible implementation, the spectral extrapolation module, when determining an extrapolation operator based on the mapping operator, the data matrix, and the PA signal, includes:
determining a second projection of the in-band signal in the time domain based on the mapping operator and the PA signal;
determining power amplifier parameters of the power amplifier model based on the second projection, the mapping operator and the data matrix;
and determining the extrapolation operator based on the power amplifier parameters.
In one possible embodiment, the formula for calculating the mapping operator is:
wherein, W
MIs a discrete fourier matrix of dimension M,
is W
MConjugate matrix of W
NIs a discrete Fourier matrix of dimension N, H
i=[0I
M0]∈R
M×NIs a selection matrix for determining the extent in the frequency domain, I, of the inband signal
MIs an M matrix, where M is the number of in-band frequency points.
In a second aspect, an embodiment of the present application provides a DPD parameter extracting apparatus applied to 5G, where the DPD parameter extracting apparatus includes:
the first processing module is used for inputting the received 5G new air interface 5G NR signal into a digital predistortion model and carrying out digital predistortion processing to obtain a predistortion signal;
the second processing module is used for inputting the predistortion signal into a determined power amplifier model and carrying out power amplification processing to obtain a power amplifier PA signal with a complete frequency spectrum;
the sampling module is used for sequentially carrying out band-limited filtering and low-speed analog/digital converter (ADC) sampling under the limited frequency on the PA signal to obtain a band-limited frequency spectrum in a truncation interval corresponding to the limited frequency in the complete frequency spectrum;
the spectrum extrapolation module is used for determining a spectrum extrapolation signal with an approximate complete spectrum based on a mapping operator for performing time-domain space projection on the truncation interval, the in-band signal corresponding to the band-limited spectrum, a data matrix formed by a basis function of the power amplifier model and a vector corresponding to the 5GNR signal, and the PA signal;
a parameter extraction module for extracting the parameters according to an iterative formula betat+1=βt-μ(XHX)-1XH(y' -x), iterating the Digital Predistortion (DPD) parameters in the digital predistortion model until the out-of-band distortion of the PA signal is smaller than a preset threshold value, and taking the finally obtained numerical value as the DPD parameters of the digital predistortion model; wherein, betatIs the DPD parameter after the t-th iteration, X is the data matrix, XHAnd a conjugate transpose matrix representing the data matrix, wherein x is the 5G NR signal, y' is the spectrum extrapolation signal, and mu is a preset parameter.
In a possible implementation, the spectrum extrapolation module is specifically configured to:
determining a first projection of an in-band signal in the time domain based on the mapping operator and the in-band signal;
determining an extrapolation operator based on the mapping operator, the data matrix, and the PA signal;
determining the spectral extrapolation signal based on the extrapolation operator and the first projection.
In a possible embodiment, said determining an extrapolation operator based on said mapping operator, said data matrix and said PA signal comprises:
determining a second projection of the in-band signal in the time domain based on the mapping operator and the PA signal;
determining power amplifier parameters of the power amplifier model based on the second projection, the mapping operator and the data matrix;
and determining the extrapolation operator based on the power amplifier parameters.
In one possible embodiment, the formula for calculating the mapping operator is:
wherein, W
MIs a discrete fourier matrix of dimension M,
is W
MConjugate matrix of W
NIs a discrete Fourier matrix of dimension N, H
i=[0I
M0]∈R
M×NIs a selection matrix for determining the extent in the frequency domain, I, of the inband signal
MIs an M matrix, where M is the number of in-band frequency points.
In a third aspect, an embodiment of the present application provides an electronic device, including: the DPD parameter extraction method comprises a processor, a storage medium and a bus, wherein the storage medium stores machine-readable instructions executable by the processor, when an electronic device runs, the processor and the storage medium communicate through the bus, and the processor executes the machine-readable instructions to execute the steps of the DPD parameter extraction method applied to 5G according to any one of the first aspects.
In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the steps of the DPD parameter extraction method applied to 5G according to any one of the first aspects are performed.
The method and the device for extracting the DPD parameters applied to the 5G, the electronic device and the medium can continuously update the DPD parameters according to various data under the condition of low sampling rate, so that the out-of-band distortion of PA signals is smaller and smaller, finally, after the DPD parameters are changed, the same DPD effect as that under the condition of high sampling rate is achieved, and the limit of the DPD to broadband application is reduced.
Detailed Description
In order to make the purpose, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it should be understood that the drawings in the present application are for illustrative and descriptive purposes only and are not used to limit the scope of protection of the present application. Additionally, it should be understood that the schematic drawings are not necessarily drawn to scale. The flowcharts used in this application illustrate operations implemented according to some embodiments of the present application. It should be understood that the operations of the flow diagrams may be performed out of order, and steps without logical context may be performed in reverse order or simultaneously. One skilled in the art, under the guidance of this application, may add one or more other operations to, or remove one or more operations from, the flowchart.
In addition, the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. The components of the embodiments of the present application, generally described and illustrated in the figures herein, can be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present application, presented in the accompanying drawings, is not intended to limit the scope of the claimed application, but is merely representative of selected embodiments of the application. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present application without making any creative effort, shall fall within the protection scope of the present application.
It should be noted that in the embodiments of the present application, the term "comprising" is used to indicate the presence of the features stated hereinafter, but does not exclude the addition of further features.
The DPD parameter extraction method applied to the 5G realizes the DPD algorithm realization for reducing the ADC sampling rate, can meet the application in a 5G broadband system, and greatly reduces the operation fund and equipment cost.
4G and previous signals have been subjected to a spectral extrapolation method, and the application is applied to 5G NR signals, and the 5G NR signals are wider in bandwidth and higher in complexity and cost requirements.
Referring to fig. 1, in the system model provided in the embodiment of the present application, a 5G New air interface (NR, New Radio) signal x passes through a Digital Predistortion (DPD) model and a Power Amplifier (PA) model in sequence to obtain a signal y, and due to nonlinearity of the power amplifier model, a frequency spectrum of y is very wide. After the band-limited sampling of the ADC, the obtained output digital signal is yBDue to yBIs limited by the sampling rate of the ADC, and the sampling bandwidth is relatively narrow.
Referring to fig. 2, a flowchart of a DPD parameter extraction method applied to 5G according to an embodiment of the present application is provided, where the method includes:
s201, inputting the received 5G new air interface 5G NR signal into a digital predistortion model, and performing digital predistortion processing to obtain a predistortion signal.
In the embodiment of the present application, the maximum bandwidth of the applicable 5G NR signal can reach 300 megabits per second (Mbps).
Inputting 5G NR signal x into digital predistortion model
The digital pre-distortion processing is carried out to pre-correct the signal distortion generated after the signal passes through the power amplifier model. Wherein, the digital predistortion model is modeled by adopting a memory polynomial fitting mode, and beta
kqIs the DPD parameter, K is the order of the polynomial and Q is the memory depth of the polynomial.
S202, inputting the predistortion signal into a determined power amplifier model, and performing power amplification treatment to obtain a power amplifier PA signal with a complete frequency spectrum.
And inputting the predistortion signal obtained after the digital predistortion model is preprocessed into the determined power amplifier model for power amplification processing to obtain a PA signal y with a complete frequency spectrum. The power amplifier adopts a filter and polynomial modeling to realize a power amplifier model of a memory polynomial.
S203, performing band-limited filtering and low-speed analog-to-digital converter (ADC) sampling under the limited frequency on the PA signal in sequence to obtain a band-limited frequency spectrum in a truncation interval corresponding to the limited frequency in the complete frequency spectrum.
Sequentially carrying out band-limited filtering and low-speed sampling under the limited frequency on the PA signal to obtain a frequency spectrum signal in a filter band, namely an in-band signal; since the sampling frequency is limited and the spectrum of the PA signal y is wide, the entire spectrum cannot be intercepted, but only a part of the spectrum.
S204, determining a spectrum extrapolation signal with an approximate complete spectrum based on a mapping operator for performing time domain space projection on the truncation interval, the in-band signal corresponding to the band-limited spectrum, a data matrix formed by a basis function of the power amplifier model and a vector corresponding to the 5G NR signal, and the PA signal.
Mapping operator P based on time domain space projection aiming at truncation intervaliIn-band signal y corresponding to band-limited spectrumBAnd a data matrix X, PA signal y formed by a basis function of the power amplifier model and a vector corresponding to the 5G NR signal reduces the band-limited spectrum to obtain an approximate complete spectrum, thereby determining a spectrum extrapolation signal y' corresponding to the approximate complete spectrum.
Since the solution of the spectrum extrapolation signal y' is related to the power amplifier model and the received 5G NR signal, the influence of the nonlinear characteristic of the power amplifier model and the influence of the memory effect of the power amplifier model under different 5G NR signal inputs are considered.
S205, according to an iterative formula betat+1=βt-μ(XHX)-1XH(y' -x), iterating the Digital Predistortion (DPD) parameters in the digital predistortion model until the out-of-band distortion of the PA signal is smaller than a preset threshold value, and taking the finally obtained numerical value as the DPD parameters of the digital predistortion model; wherein, betatIs the DPD parameter after the t-th iteration, X is the data matrix, XHAnd a conjugate transpose matrix representing the data matrix, wherein x is the 5G NR signal, y' is the spectrum extrapolation signal, and mu is a preset parameter.
In order to weaken signal distortion caused by a power amplifier model, based on a direct learning structure, beta can be deduced from a Gauss-Newton (Gauss-Newton) iterative formulat+1=βt-μ(XHX)-1XH(y' -x), digital predistortion parameters for the digital predistortion model (i.e., β in the equation of step S201)kq) Iteration is carried out until the out-of-band distortion of the PA signal y is less than a certain preset threshold value so as to obtain the number finallyThe values are used as DPD parameters of the digital predistortion model; wherein, betatIs the DPD parameter after the t-th iteration, X is a data matrix, XHAnd the matrix represents a conjugate transpose matrix of the data matrix, x is a 5G NR signal, y' is a spectrum extrapolation signal, and mu is a preset parameter of an adjustable integer value.
The iterative formula derivation process is as follows: under the low-speed sampling condition under the band-limited filtering and limiting frequency, based on the overall error solving mode, the objective function can be written as:
wherein β is a DPD parameter.
The first derivative of the objective function is solved to obtain a Jacobian matrix which is:
the second derivative of the objective function is solved to obtain a Hessian matrix which is:
the iterative formula for obtaining the DPD parameters by the Gaussian-Newton iterative formula is as follows:
in the iterative formula of the DPD model parameters, the solution is performed based on the total error of direct learning, and compared with a method based on the segmentation error solution, the calculation of Discrete Fourier Transform (DFT) and Inverse Discrete Fourier Transform (IDFT) is avoided when the Hessian matrix is calculated, so that the calculation is simpler.
Referring to fig. 3, in a possible embodiment, the determining a spectral extrapolation signal having an approximately complete spectrum based on a mapping operator performing time-domain spatial projection on the truncated interval, an in-band signal corresponding to the band-limited spectrum, a data matrix formed by basis functions of the power amplifier model and a vector corresponding to the 5G NR signal, and the PA signal includes:
s301, determining a first projection of the in-band signal in a time domain based on the mapping operator and the in-band signal.
S302, determining an extrapolation operator based on the mapping operator, the data matrix and the PA signal.
S303, determining the frequency spectrum extrapolation signal based on the extrapolation operator and the first projection.
With reference to steps S301 to S303, a specific derivation method of the spectrum extrapolation signal is:
based on a mapping operator P for band-limited spectra (i.e. truncation intervals)iFor in-band signal y corresponding to band-limited spectrumBPerforming projection on a time domain function space to obtain a first projection of an in-band signal in a time domain, wherein the first projection is as follows:
determining an extrapolation operator xi for carrying out spectrum extrapolation reduction on the in-band signal based on a data matrix, a PA signal and a mapping operator which are formed by a basis function of the power amplifier model and a vector of the 5G NR signal; the specific derivation method of the extrapolation operator is shown in steps S401-S403.
In-band signal with incomplete spectrum based on extrapolation operator xi
And performing full-spectrum extrapolation reduction to obtain a spectrum extrapolation signal:
the extrapolation operator can extrapolate the in-band signal under the limited sampling frequency to a signal with complete frequency spectrum which can be obtained under the expected sampling rate, and meanwhile, the solving mode of the extrapolation operator depends on the selection of the power amplifier model and the vector corresponding to the received 5G NR signal, so that the influence of how to balance the static nonlinear characteristic of the power amplifier model is considered, and the influence of how to deal with the memory effect of the power amplifier model is also considered.
Referring to fig. 4, in a possible implementation, the determining an extrapolation operator based on the mapping operator, the data matrix, and the PA signal includes:
s401, determining a second projection of the in-band signal in a time domain based on the projection operator and the PA signal.
S402, determining power amplifier parameters of the power amplifier model based on the second projection, the mapping operator and the data matrix.
And S403, determining the extrapolation operator based on the power amplifier parameters.
With reference to steps S401 to S403, a specific derivation method of the extrapolation operator is as follows:
because the power amplifier model is selected, the data matrix formed by the basis function of the power amplifier model and the 5G NR signal vector is X, and the power amplifier parameter of the corresponding power amplifier model is a, then the PA signal y can be expressed as:
y=Xa
multiplying the mapping operator P on both sides of the above equationiThen, obtaining:
Piy=PiXa
because ADC for sampling the PA signal is close to perfect sampling, the band-limited spectrum obtained after sampling does not generate signal distortion compared with the spectrum corresponding to the truncation interval in the PA signal (the mapping operator only performs time domain projection on the band-limited spectrum part, and does not perform projection on the part except the truncation interval), therefore, the second projection of the in-band signal obtained by the product of the mapping operator and the PA signal t in the time domain and the first projection of the in-band signal obtained by the product of the mapping operator and the in-band signal in the previous step in the time domain
Are the same. The above equation thus becomes:
therefore, according to the formula, the power amplifier parameters under the current sampling frequency are obtained by the least square method:
estimating a frequency spectrum extrapolation signal approximate to the PA signal according to the obtained power amplifier parameter a as follows:
and target equation
And comparing, wherein the extrapolation operator to be solved is as follows:
the extrapolation operator xi calculated by the method can extrapolate the in-band signal under the limited sampling frequency to the spectrum extrapolation signal obtained under the expected sampling rate, and meanwhile, the solving method of the extrapolation operator depends on the selection of the power amplifier model, so that the influence of how to balance the static nonlinear characteristic of the power amplifier model and how to deal with the memory effect of the power amplifier model are considered.
In one possible embodiment, the formula for calculating the mapping operator is:
wherein, W
MIs a discrete fourier matrix of dimension M,
is W
MConjugate matrix of W
NIs a discrete Fourier matrix of dimension N, H
i=[0I
M0]∈R
M×NIs a selection matrix for determining the extent in the frequency domain, I, of the inband signal
MIs an M matrix, where M is the number of in-band frequency points.
Referring to fig. 5, after the parameters in the digital predistortion are iterated, in the present application, when the input signal is a 5G signal with a signal bandwidth of 160Mbps and the feedback path employs an ADC with a sampling rate of 250Mbps, the Adjacent Channel Leakage Ratio (ACLR) can reach-60 dB, and the effect of the high-speed ADC processing with a sampling rate of 900Mbps is equivalent to that of upper ACLR0 being an upper direct Adjacent Channel (where lower ACLR0 is an upper direct Adjacent Channel, upper ACLR1 is an upper direct Adjacent Channel with a distance BW (one bandwidth), and lower ACLR1 is a lower Adjacent Channel with a distance BW).
Referring to fig. 6, the input Signal is a 5G NR Signal with a bandwidth of 160Mbps, and the input Signal is a spectrum diagram corresponding to a Perfect Signal (Perfect Signal), a Signal without DPD (PA Signal) and a Signal (DPD + PA Signal) with DPD parameters determined by the method of the present application (the ADC of the feedback path uses a sampling rate of 250 Mbps).
Based on the same inventive concept, the embodiment of the present application further provides a DPD parameter extraction device applied to 5G corresponding to the DPD parameter extraction method applied to 5G, and since the principle of the device in the embodiment of the present application to solve the problem is similar to the DPD parameter extraction method applied to 5G described above in the embodiment of the present application, the implementation of the device may refer to the implementation of the method, and repeated details are not repeated.
Referring to fig. 7, a schematic diagram of a DPD parameter extracting apparatus applied to 5G according to an embodiment of the present application is shown, where the DPD parameter extracting apparatus includes:
a first processing module 701, configured to input the received 5G new air interface 5G NR signal into a digital predistortion model, and perform digital predistortion processing to obtain a predistortion signal;
a second processing module 702, configured to input the predistortion signal into a determined power amplifier model, and perform power amplification processing to obtain a power amplifier PA signal having a complete frequency spectrum;
a sampling module 703, configured to perform band-limited filtering and low-speed analog-to-digital converter ADC sampling on the PA signal sequentially under a limited frequency to obtain a band-limited spectrum in a truncation interval corresponding to the limited frequency in the complete spectrum;
a spectrum extrapolation module 704, configured to determine a spectrum extrapolation signal having an approximately complete spectrum based on a mapping operator performing time-domain spatial projection on the truncated interval, an in-band signal corresponding to the band-limited spectrum, a data matrix formed by a basis function of the power amplifier model and a vector corresponding to the 5G NR signal, and the PA signal;
a parameter extraction module 705 for extracting the parameter according to the iterative formula βt+1=βt-μ(XHX)-1XH(y' -x), iterating the Digital Predistortion (DPD) parameters in the digital predistortion model until the out-of-band distortion of the PA signal is smaller than a preset threshold value, and taking the finally obtained numerical value as the DPD parameters of the digital predistortion model; wherein, betatIs the DPD parameter after the t-th iteration, X is the data matrix, XHAnd a conjugate transpose matrix representing the data matrix, wherein x is the 5G NR signal, y' is the spectrum extrapolation signal, and mu is a preset parameter.
In a possible implementation, the spectrum extrapolation module 704 is specifically configured to:
determining a first projection of an in-band signal in the time domain based on the mapping operator and the in-band signal;
determining an extrapolation operator based on the mapping operator, the data matrix, and the PA signal;
determining the spectral extrapolation signal based on the extrapolation operator and the first projection.
In a possible implementation, the spectrum extrapolation module 705, when determining an extrapolation operator based on the mapping operator, the data matrix, and the PA signal, includes:
determining a second projection of the in-band signal in the time domain based on the mapping operator and the PA signal;
determining power amplifier parameters of the power amplifier model based on the second projection, the mapping operator and the data matrix;
and determining the extrapolation operator based on the power amplifier parameters.
In one possible embodiment, the formula for calculating the mapping operator is:
wherein, W
MIs a discrete fourier matrix of dimension M,
is W
MConjugate matrix of W
NIs a discrete Fourier matrix of dimension N, H
i=[0I
M0]∈R
M×NIs a selection matrix for determining the extent in the frequency domain, I, of the inband signal
MIs an M matrix, where M is the number of in-band frequency points.
Referring to fig. 8, an electronic device 800 according to an embodiment of the present application includes: a processor 801, a memory 802 and a bus, wherein the memory 802 stores machine-readable instructions executable by the processor 801, when the electronic device is operated, the processor 801 communicates with the memory 802 through the bus, and the processor 801 executes the machine-readable instructions to execute the steps of the DPD parameter extraction method applied to 5G as described above.
Specifically, the memory 802 and the processor 801 can be general-purpose memories and processors, which are not specifically limited herein, and the DPD parameter extracting method applied to 5G can be performed when the processor 801 runs a computer program stored in the memory 802.
Corresponding to the method, an embodiment of the present application further provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and the computer program is executed by a processor to perform the steps of the DPD parameter extraction method applied to 5G.
It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working processes of the system and the apparatus described above may refer to corresponding processes in the method embodiments, and are not described in detail in this application. In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. The above-described apparatus embodiments are merely illustrative, and for example, the division of the modules is merely a logical division, and there may be other divisions in actual implementation, and for example, a plurality of modules or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection of devices or modules through some communication interfaces, and may be in an electrical, mechanical or other form.
The modules described as separate parts may or may not be physically separate, and parts displayed as modules may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a non-volatile computer-readable storage medium executable by a processor. Based on such understanding, the technical solution of the present application or portions thereof that substantially contribute to the prior art may be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: various media capable of storing program codes, such as a U disk, a removable hard disk, a ROM, a RAM, a magnetic disk, or an optical disk.
The above description is only for the specific embodiments of the present application, but the scope of the present application is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present application, and shall be covered by the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.