[go: up one dir, main page]

CN113114138B - DPD parameter extraction method and device applied to 5G, electronic equipment and medium - Google Patents

DPD parameter extraction method and device applied to 5G, electronic equipment and medium Download PDF

Info

Publication number
CN113114138B
CN113114138B CN202110434196.5A CN202110434196A CN113114138B CN 113114138 B CN113114138 B CN 113114138B CN 202110434196 A CN202110434196 A CN 202110434196A CN 113114138 B CN113114138 B CN 113114138B
Authority
CN
China
Prior art keywords
signal
band
spectrum
dpd
power amplifier
Prior art date
Legal status (The legal status 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 status listed.)
Expired - Fee Related
Application number
CN202110434196.5A
Other languages
Chinese (zh)
Other versions
CN113114138A (en
Inventor
尹宇芳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
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 Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202110434196.5A priority Critical patent/CN113114138B/en
Publication of CN113114138A publication Critical patent/CN113114138A/en
Application granted granted Critical
Publication of CN113114138B publication Critical patent/CN113114138B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03FAMPLIFIERS
    • H03F3/00Amplifiers with only discharge tubes or only semiconductor devices as amplifying elements
    • H03F3/20Power amplifiers, e.g. Class B amplifiers, Class C amplifiers

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Amplifiers (AREA)

Abstract

本申请提供了一种应用于5G的DPD参数提取方法、装置、电子设备及介质,方法包括:将接收到的5G NR信号输入到DPD模型中得到预失真信号;将预失真信号输入到功放模型中得到频谱完整的功放PA信号;对PA信号进行带限滤波及低速采样得到完整频谱中处于限制频率对应的截断区间中的带限频谱;基于映射算子、带限频谱对应的带内信号、由功放模型的基函数和5G NR信号的矢量构成的数据矩阵、PA信号,确定拥有近似完整频谱的频谱外推信号;根据迭代公式βt+1=βt‑μ(XHX)‑1XH(y′‑x)对DPD模型中的DPD参数进行迭代直至收敛。本申请在对PA信号的采样频率较低的情况下,使用迭代收敛后的DPD参数,可达到与高采样速率情况时同样的DPD效果,从而大幅降低了设备成本和运营资金。

Figure 202110434196

The present application provides a DPD parameter extraction method, device, electronic device and medium applied to 5G. The method includes: inputting a received 5G NR signal into a DPD model to obtain a predistorted signal; inputting the predistorted signal into a power amplifier model Obtain the PA signal of the power amplifier with complete spectrum; perform band-limit filtering and low-speed sampling on the PA signal to obtain the band-limited spectrum in the truncation interval corresponding to the limiting frequency in the complete spectrum; based on the mapping operator, the in-band signal corresponding to the band-limited spectrum, The data matrix and PA signal composed of the basis function of the power amplifier model and the vector of the 5G NR signal determine the spectral extrapolation signal with an approximately complete spectrum; according to the iterative formula β t+1t ‑μ(X H X) ‑1 X H (y′‑x) iterates the DPD parameters in the DPD model until convergence. In the present application, when the sampling frequency of the PA signal is low, the DPD parameters after iterative convergence can be used to achieve the same DPD effect as in the case of a high sampling rate, thereby greatly reducing equipment costs and operating funds.

Figure 202110434196

Description

DPD parameter extraction method and device applied to 5G, electronic equipment and medium
Technical Field
The present application relates to the field of wireless communication technologies, and in particular, to a method and an apparatus for extracting DPD parameters for 5G, an electronic device, and a medium.
Background
In order to obtain larger output signal power, a Power Amplifier (PA) generally operates in an operating range close to a saturation point, so that the nonlinear region of the PA causes signal distortion.
Digital Pre-Distortion (DPD) is one of the most basic building blocks in the current wireless communication system, and distorts the signal generated when the pa operates in its non-linear region by distorting the signal in the opposite direction of the rf pa.
The traditional DPD parameter extraction algorithm relies on a sampling rate not less than 5 times the bandwidth of the original signal to obtain the feedback signal. However, with the arrival of 5G, the carrier bandwidth is wider and wider. This places extremely high demands on Analog-to-Digital converters (ADCs). A 5G DPD with 160M bandwidth, the ADC of the feedback path needs 900M sampling rate, and the power consumption and cost of the high-speed and high-precision ADC are very high, so that the DPD is greatly limited when being applied to a broadband.
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:
Figure GDA0003058040910000031
wherein, WMIs a discrete fourier matrix of dimension M,
Figure GDA0003058040910000032
is WMConjugate matrix of WNIs a discrete Fourier matrix of dimension N, Hi=[0IM0]∈RM×NIs a selection matrix for determining the extent in the frequency domain, I, of the inband signalMIs 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:
Figure GDA0003058040910000041
wherein, WMIs a discrete fourier matrix of dimension M,
Figure GDA0003058040910000042
is WMConjugate matrix of WNIs a discrete Fourier matrix of dimension N, Hi=[0IM0]∈RM×NIs a selection matrix for determining the extent in the frequency domain, I, of the inband signalMIs 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.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained from the drawings without inventive effort.
FIG. 1 illustrates a system model diagram provided by an embodiment of the present application;
fig. 2 shows a flowchart of a DPD parameter extraction method applied to 5G according to an embodiment of the present application;
fig. 3 shows a flowchart of a DPD parameter extraction method applied to 5G according to an embodiment of the present application;
fig. 4 shows a flowchart of a DPD parameter extraction method applied to 5G according to an embodiment of the present application;
FIG. 5 is a table diagram of simulation result data provided by an embodiment of the present application;
FIG. 6 is a diagram illustrating simulation results provided by an embodiment of the present application;
fig. 7 is a schematic structural diagram illustrating a DPD parameter extraction apparatus applied to 5G according to an embodiment of the present application;
fig. 8 shows a schematic structural diagram of an electronic device provided in an embodiment of the present application.
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
Figure GDA0003058040910000071
Figure GDA0003058040910000072
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 betakqIs 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:
Figure GDA0003058040910000091
wherein β is a DPD parameter.
The first derivative of the objective function is solved to obtain a Jacobian matrix which is:
Figure GDA0003058040910000092
the second derivative of the objective function is solved to obtain a Hessian matrix which is:
Figure GDA0003058040910000093
the iterative formula for obtaining the DPD parameters by the Gaussian-Newton iterative formula is as follows:
Figure GDA0003058040910000094
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:
Figure GDA0003058040910000101
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
Figure GDA0003058040910000102
And performing full-spectrum extrapolation reduction to obtain a spectrum extrapolation signal:
Figure GDA0003058040910000103
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
Figure GDA0003058040910000111
Are the same. The above equation thus becomes:
Figure GDA0003058040910000112
therefore, according to the formula, the power amplifier parameters under the current sampling frequency are obtained by the least square method:
Figure GDA0003058040910000113
estimating a frequency spectrum extrapolation signal approximate to the PA signal according to the obtained power amplifier parameter a as follows:
Figure GDA0003058040910000114
and target equation
Figure GDA0003058040910000121
And comparing, wherein the extrapolation operator to be solved is as follows:
Figure GDA0003058040910000122
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:
Figure GDA0003058040910000123
wherein, WMIs a discrete fourier matrix of dimension M,
Figure GDA0003058040910000124
is WMConjugate matrix of WNIs a discrete Fourier matrix of dimension N, Hi=[0IM0]∈RM×NIs a selection matrix for determining the extent in the frequency domain, I, of the inband signalMIs 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:
Figure GDA0003058040910000141
wherein, WMIs a discrete fourier matrix of dimension M,
Figure GDA0003058040910000142
is WMConjugate matrix of WNIs a discrete Fourier matrix of dimension N, Hi=[0IM0]∈RM×NIs a selection matrix for determining the extent in the frequency domain, I, of the inband signalMIs 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.

Claims (10)

1.一种应用于5G的DPD参数提取方法,其特征在于,所述方法包括:1. a DPD parameter extraction method applied to 5G, is characterized in that, described method comprises: 将接收到的5G新空口5G NR信号输入到数字预失真模型中,进行数字预失真处理,得到预失真信号;Input the received 5G NR 5G NR signal into the digital pre-distortion model, perform digital pre-distortion processing, and obtain a pre-distorted signal; 将所述预失真信号输入到确定的功放模型中,进行功率放大处理,得到拥有完整频谱的功放PA信号;Inputting the predistorted signal into the determined power amplifier model, and performing power amplification processing to obtain a power amplifier PA signal with a complete spectrum; 对所述PA信号依次进行带限滤波以及限制频率下的低速模拟/数字转换器ADC采样,得到所述完整频谱中处于所述限制频率对应的截断区间中的带限频谱;Perform band-limit filtering and low-speed analog/digital converter ADC sampling on the PA signal in turn to obtain a band-limited spectrum in the truncation interval corresponding to the limiting frequency in the complete spectrum; 基于针对所述截断区间进行时域空间投影的映射算子、所述带限频谱对应的带内信号、由所述功放模型的基函数和所述5G NR信号对应的矢量构成的数据矩阵、所述PA信号,对所述带限频谱进行还原,得到还原后的带限频谱对应的频谱外推信号;Based on the mapping operator that performs time-domain space projection for the truncation interval, the in-band signal corresponding to the band-limited spectrum, the data matrix composed of the basis function of the power amplifier model and the vector corresponding to the 5G NR signal, the the PA signal, restore the band-limited spectrum, and obtain a spectrum extrapolation signal corresponding to the restored band-limited spectrum; 根据迭代公式βt+1=βt-μ(XHX)-1XH(y′-x),对所述数字预失真模型中的数字预失真DPD参数进行迭代,直至所述PA信号的带外失真小于预设阈值,以将最终得到的数值作为所述数字预失真模型的DPD参数;其中,βt为第t次迭代后的所述DPD参数,X为所述数据矩阵,XH表示所述数据矩阵的共轭转置矩阵,x为所述5G NR信号,y′为所述频谱外推信号,μ为预设参量。According to the iterative formula β t+1t -μ(X H X) -1 X H (y'-x), iterate the digital predistortion DPD parameters in the digital predistortion model until the PA signal The out-of-band distortion of β is less than the preset threshold, so that the final value obtained is used as the DPD parameter of the digital predistortion model; wherein, β t is the DPD parameter after the t-th iteration, X is the data matrix, X H represents the conjugate transpose matrix of the data matrix, x is the 5G NR signal, y' is the spectrum extrapolation signal, and μ is a preset parameter. 2.根据权利要求1所述的应用于5G的DPD参数提取方法,其特征在于,所述基于针对所述截断区间进行时域空间投影的映射算子、所述带限频谱对应的带内信号、由所述功放模型的基函数和所述5G NR信号对应的矢量构成的数据矩阵、所述PA信号,对所述带限频谱进行还原,得到还原后的带限频谱对应的频谱外推信号,包括:2 . The DPD parameter extraction method applied to 5G according to claim 1 , wherein the mapping operator based on the time-domain space projection for the truncation interval, the in-band signal corresponding to the band-limited spectrum , the data matrix formed by the basis function of the power amplifier model and the vector corresponding to the 5G NR signal, the PA signal, restore the band-limited spectrum, and obtain the spectrum extrapolation signal corresponding to the restored band-limited spectrum ,include: 基于所述映射算子和所述带内信号,确定带内信号在时域里的第一投影;determining a first projection of the in-band signal in the time domain based on the mapping operator and the in-band signal; 基于所述映射算子、所述数据矩阵和所述PA信号,确定外推算子;determining an extrapolation operator based on the mapping operator, the data matrix and the PA signal; 基于所述外推算子和所述第一投影,确定所述频谱外推信号。Based on the extrapolation operator and the first projection, the spectrally extrapolated signal is determined. 3.根据权利要求2所述的应用于5G的DPD参数提取方法,其特征在于,所述基于所述映射算子、所述数据矩阵和所述PA信号,确定外推算子,包括:3. The method for extracting DPD parameters applied to 5G according to claim 2, wherein, determining an extrapolation operator based on the mapping operator, the data matrix and the PA signal, comprising: 基于所述映射算子和所述PA信号,确定带内信号在时域里的第二投影;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; 基于所述功放参数,确定所述外推算子。The extrapolation operator is determined based on the power amplifier parameters. 4.根据权利要求1至3中任一项所述的应用于5G的DPD参数提取方法,其特征在于,计算所述映射算子的公式为:4. The DPD parameter extraction method applied to 5G according to any one of claims 1 to 3, wherein the formula for calculating the mapping operator is:
Figure FDA0003513021580000021
Figure FDA0003513021580000021
其中,WM为维度为M的离散傅里叶矩阵,
Figure FDA0003513021580000022
是WM的共轭矩阵,WN为维度为N的离散傅里叶矩阵,Hi=[0 IM 0]∈RM×N是选择矩阵,用于确定所述带内信号在频域中的范围,IM为M×M矩阵,M是带内频率点数。
where W M is a discrete Fourier matrix of dimension M,
Figure FDA0003513021580000022
is the conjugate matrix of W M , W N is the discrete Fourier matrix of dimension N, H i =[0 I M 0]∈R M×N is the selection matrix, used to determine the in-band signal in the frequency domain In the range, IM is an M×M matrix, and M is the number of in-band frequency points.
5.一种应用于5G的DPD参数提取装置,其特征在于,所述装置包括:5. A DPD parameter extraction device applied to 5G, wherein the device comprises: 第一处理模块,用于将接收到的5G新空口5G NR信号输入到数字预失真模型中,进行数字预失真处理,得到预失真信号;The first processing module is used to input the received 5G new air interface 5G NR signal into the digital pre-distortion model, perform digital pre-distortion processing, and obtain a pre-distorted signal; 第二处理模块,用于将所述预失真信号输入到确定的功放模型中,进行功率放大处理,得到拥有完整频谱的功放PA信号;a second processing module, configured to input the predistorted signal into the determined power amplifier model, perform power amplification processing, and obtain a power amplifier PA signal with a complete spectrum; 采样模块,用于对所述PA信号依次进行带限滤波以及限制频率下的低速模拟/数字转换器ADC采样,得到所述完整频谱中处于所述限制频率对应的截断区间中的带限频谱;a sampling module, configured to sequentially perform band-limit filtering and low-speed analog/digital converter ADC sampling under the limiting frequency on the PA signal, to obtain a band-limited spectrum in the truncation interval corresponding to the limiting frequency in the complete spectrum; 频谱外推模块,用于基于针对所述截断区间进行时域空间投影的映射算子、所述带限频谱对应的带内信号、由所述功放模型的基函数和所述5G信号对应的矢量构成的数据矩阵、所述PA信号,对所述带限频谱进行还原,得到还原后的带限频谱对应的频谱外推信号;A spectrum extrapolation module for performing a mapping operator for time-domain space projection based on the truncation interval, an in-band signal corresponding to the band-limited spectrum, a vector corresponding to the basis function of the power amplifier model and the 5G signal The formed data matrix and the PA signal restore the band-limited spectrum to obtain a spectrum extrapolation signal corresponding to the restored band-limited spectrum; 参数提取模块,用于根据迭代公式βt+1=βt-μ(XHX)-1XH(y′-x),对所述数字预失真模型中的数字预失真DPD参数进行迭代,直至所述PA信号的带外失真小于预设阈值,以将最终得到的数值作为所述数字预失真模型的DPD参数;其中,βt为第t次迭代后的所述DPD参数,X为所述数据矩阵,XH表示所述数据矩阵的共轭转置矩阵,x为所述5G NR信号,y′为所述频谱外推信号,μ为预设参量。A parameter extraction module, configured to iterate the digital predistortion DPD parameters in the digital predistortion model according to the iterative formula β t+1t -μ(X H X) -1 X H (y'-x) , until the out-of-band distortion of the PA signal is less than the preset threshold, to use the final value as the DPD parameter of the digital predistortion model; where β t is the DPD parameter after the t-th iteration, and X is the In the data matrix, X H represents the conjugate transpose matrix of the data matrix, x is the 5G NR signal, y' is the spectrum extrapolation signal, and μ is a preset parameter. 6.根据权利要求5所述的应用于5G的DPD参数提取装置,其特征在于,所述频谱外推模块,具体用于:6. The DPD parameter extraction device applied to 5G according to claim 5, wherein the spectrum extrapolation module is specifically used for: 基于所述映射算子和所述带内信号,确定带内信号在时域里的第一投影;determining a first projection of the in-band signal in the time domain based on the mapping operator and the in-band signal; 基于所述映射算子、所述数据矩阵和所述PA信号,确定外推算子;determining an extrapolation operator based on the mapping operator, the data matrix and the PA signal; 基于所述外推算子和所述第一投影,确定所述频谱外推信号。Based on the extrapolation operator and the first projection, the spectrally extrapolated signal is determined. 7.根据权利要求6所述的应用于5G的DPD参数提取装置,其特征在于,所述频谱外推模块用于基于所述映射算子、所述数据矩阵和所述PA信号,确定外推算子时,包括:7 . The DPD parameter extraction device applied to 5G according to claim 6 , wherein the spectrum extrapolation module is configured to determine an extrapolation calculation based on the mapping operator, the data matrix and the PA signal. 8 . child hours, including: 基于所述映射算子和所述PA信号,确定带内信号在时域里的第二投影;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; 基于所述功放参数,确定所述外推算子。The extrapolation operator is determined based on the power amplifier parameters. 8.根据权利要求5至7中任一项所述的应用于5G的DPD参数提取装置,其特征在于,计算所述映射算子的公式为:8. The DPD parameter extraction device applied to 5G according to any one of claims 5 to 7, wherein the formula for calculating the mapping operator is:
Figure FDA0003513021580000041
Figure FDA0003513021580000041
其中,WM为维度为M的离散傅里叶矩阵,
Figure FDA0003513021580000042
是WM的共轭矩阵,WN为维度为N的离散傅里叶矩阵,Hi=[0 IM 0]∈RM×N是选择矩阵,用于确定所述带内信号在频域中的范围,IM为M×M矩阵,M是带内频率点数。
where W M is a discrete Fourier matrix of dimension M,
Figure FDA0003513021580000042
is the conjugate matrix of W M , W N is the discrete Fourier matrix of dimension N, H i =[0 I M 0]∈R M×N is the selection matrix, used to determine the in-band signal in the frequency domain In the range, IM is an M×M matrix, and M is the number of in-band frequency points.
9.一种电子设备,其特征在于,包括:处理器、存储介质和总线,所述存储介质存储有所述处理器可执行的机器可读指令,当电子设备运行时,所述处理器与所述存储介质之间通过总线通信,所述处理器执行所述机器可读指令,以执行如权利要求1至4任一项所述的应用于5G的DPD参数提取方法的步骤。9. An electronic device, comprising: a processor, a storage medium, and a bus, wherein the storage medium stores machine-readable instructions executable by the processor, and when the electronic device runs, the processor and the The storage media communicate through a bus, and the processor executes the machine-readable instructions to execute the steps of the method for extracting DPD parameters for 5G according to any one of claims 1 to 4. 10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器运行时执行如权利要求1至4任一项所述的应用于5G的DPD参数提取方法的步骤。10. A computer-readable storage medium, wherein a computer program is stored on the computer-readable storage medium, and the computer program executes the application according to any one of claims 1 to 4 when the computer program is run by a processor The steps of the DPD parameter extraction method for 5G.
CN202110434196.5A 2021-04-22 2021-04-22 DPD parameter extraction method and device applied to 5G, electronic equipment and medium Expired - Fee Related CN113114138B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110434196.5A CN113114138B (en) 2021-04-22 2021-04-22 DPD parameter extraction method and device applied to 5G, electronic equipment and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110434196.5A CN113114138B (en) 2021-04-22 2021-04-22 DPD parameter extraction method and device applied to 5G, electronic equipment and medium

Publications (2)

Publication Number Publication Date
CN113114138A CN113114138A (en) 2021-07-13
CN113114138B true CN113114138B (en) 2022-03-25

Family

ID=76719266

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110434196.5A Expired - Fee Related CN113114138B (en) 2021-04-22 2021-04-22 DPD parameter extraction method and device applied to 5G, electronic equipment and medium

Country Status (1)

Country Link
CN (1) CN113114138B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2025155041A1 (en) * 2024-01-15 2025-07-24 삼성전자주식회사 Electronic device and method for digital pre-distortion

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118282814A (en) * 2024-03-06 2024-07-02 密卡思(深圳)电讯有限公司 Undersampled digital predistortion method and system based on spectrum extrapolation

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3005554B1 (en) * 2013-06-05 2021-01-27 Telefonaktiebolaget LM Ericsson (publ) Digital predistortion of wideband power amplifiers with reduced observation bandwidth
US9608676B2 (en) * 2013-07-12 2017-03-28 Analog Devices Global Digital pre-distortion systems in transmitters
CN107276546B (en) * 2016-04-08 2020-05-08 大唐移动通信设备有限公司 Digital pre-distortion processing method and device
CN105958952B (en) * 2016-04-25 2019-06-21 中国科学技术大学 A single-bit digital predistortion method for power amplifier
CN106877825B (en) * 2017-01-25 2020-04-24 东南大学 Digital predistortion device and method based on band-limited simplified nonlinear filter
US11082013B2 (en) * 2017-05-12 2021-08-03 Mediatek Inc. Method of reducing memory effect of power amplifier
US10320340B1 (en) * 2018-01-11 2019-06-11 Analog Devices Global Unlimited Company Frequency-shaped digital predistortion
CN109996205B (en) * 2019-04-12 2021-12-07 成都工业学院 Sensor data fusion method and device, electronic equipment and storage medium
CN110266276B (en) * 2019-05-17 2022-12-02 杭州电子科技大学 Low-speed digital predistortion method for 5G ultra-wideband power amplifier
CN111030954B (en) * 2019-12-04 2022-02-25 辽宁工程技术大学 A predistortion method for wideband power amplifier based on multi-channel sampling based on compressed sensing
CN111082757B (en) * 2019-12-27 2023-05-23 东方红卫星移动通信有限公司 Method for applying adaptive band-limited digital predistortion in broadband satellite communication system

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2025155041A1 (en) * 2024-01-15 2025-07-24 삼성전자주식회사 Electronic device and method for digital pre-distortion

Also Published As

Publication number Publication date
CN113114138A (en) 2021-07-13

Similar Documents

Publication Publication Date Title
CN113114138B (en) DPD parameter extraction method and device applied to 5G, electronic equipment and medium
CN102487367B (en) Adaptive power amplifier digital baseband predistortion method
WO2015096735A1 (en) Digital pre-distortion parameter obtaining method and pre-distortion system
CN104796091B (en) Power amplifier modeling and digital pre-distortion method based on segmentation memory polynomial
CN113949350B (en) Digital predistortion method and system based on baseband-RF joint optimization
CN105393452B (en) Baseband equivalence Volterra series for digital predistorter in multiband power amplifier
CN106453172A (en) Memory polynomial digital pre-distortion method based on piecewise linear function
CN106877825A (en) Digital predistortion device and method based on band-limited simplified nonlinear filter
JP2013542696A (en) Nonlinear model with tap output normalization
CN103888396A (en) Predistortion method based on self-adaptive variable step size iteration
Liu et al. Digital predistortion linearization with deep neural networks for 5G power amplifiers
Pan et al. A predistortion algorithm based on accurately solving the reverse function of memory polynomial model
CN112202695A (en) An Undersampling Digital Predistortion Method and System Based on Landweber Iterative Algorithm
CN102142814B (en) Power amplifier related device, power amplifier predistortion system and modeling method
CN107786174B (en) Circuit and method for predistortion
CN103685110A (en) Predistortion processing method and system and predistortion factor arithmetic unit
CN109740225B (en) A Method for Calculation and Evaluation of Broadband Traveling Wave Tubes
CN103179074A (en) Orthogonal polynomial based adaptive predistortion system and method
CN114123991B (en) A method for implementing a digital predistorter for 5G wireless networks
CN102611661B (en) Predistortion device and method based on precise inverse solution memory polynomial model equation
CN115296619A (en) Memory polynomial digital predistortion method based on multi-dimensional amplitude switching
CN111988002B (en) Digital predistortion method, device, equipment and storage medium for MIMO power amplifier
CN106411271B (en) A kind of pre-distortion device and its parameter selection method for power amplifier
CN110266276B (en) Low-speed digital predistortion method for 5G ultra-wideband power amplifier
CN106105032B (en) System and method for adaptive filter

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
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220325