CN111669132B - Method for improving index accuracy of short wave digital predistortion processing - Google Patents

Method for improving index accuracy of short wave digital predistortion processing Download PDF

Info

Publication number
CN111669132B
CN111669132B CN202010534060.7A CN202010534060A CN111669132B CN 111669132 B CN111669132 B CN 111669132B CN 202010534060 A CN202010534060 A CN 202010534060A CN 111669132 B CN111669132 B CN 111669132B
Authority
CN
China
Prior art keywords
lookup table
input signal
power amplifier
power
predistortion
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.)
Active
Application number
CN202010534060.7A
Other languages
Chinese (zh)
Other versions
CN111669132A (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.)
Shaanxi Fenghuo Communication Group Co Ltd
Original Assignee
Shaanxi Fenghuo Communication Group Co Ltd
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 Shaanxi Fenghuo Communication Group Co Ltd filed Critical Shaanxi Fenghuo Communication Group Co Ltd
Priority to CN202010534060.7A priority Critical patent/CN111669132B/en
Publication of CN111669132A publication Critical patent/CN111669132A/en
Application granted granted Critical
Publication of CN111669132B publication Critical patent/CN111669132B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03FAMPLIFIERS
    • H03F1/00Details of amplifiers with only discharge tubes, only semiconductor devices or only unspecified devices as amplifying elements
    • H03F1/32Modifications of amplifiers to reduce non-linear distortion
    • H03F1/3241Modifications of amplifiers to reduce non-linear distortion using predistortion circuits
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Abstract

The application belongs to the technical field of adaptive digital predistortion, and discloses a method for improving the index precision of short wave digital predistortion processing, which comprises the following steps: step 1, amplifying an original input signal by a power amplifier to obtain an amplified output signal; step 2, processing the original input signal and the output signal after the power amplification by adopting a predistortion algorithm module to obtain the lookup table parameters of the inverse characteristics of the power amplification distortion; the predistortion algorithm module adopts an integer integrated small-scale table index structure; step 3, correcting an input signal of the power amplifier by adopting a lookup table parameter of the distortion inverse characteristic of the power amplifier to obtain a signal after predistortion treatment; the method designs an integral small-scale table index structure aiming at the characteristics of a short wave signal system, so that the table index precision can be improved, the performance of a digital predistortion system is improved, and compared with the integral table index structure, the third-order intermodulation index (2-30 MHz) of a transmitter can be improved by more than 5 dB.

Description

Method for improving index accuracy of short wave digital predistortion processing
Technical Field
The application relates to the technical field of self-adaptive digital predistortion, in particular to a method for improving the index precision of short-wave digital predistortion processing, which is suitable for improving the efficiency of a short-wave radio frequency power amplifier and the performance of a digital predistortion system.
Background
The frequency of the wireless electromagnetic wave with the frequency range of 1.6 MHz-30 MHz is generally called a short wave frequency band, and the short wave frequency is utilized for carrying out worldwide broadcast transmission for unidirectional communication, which is generally called short wave broadcast; a radio capable of receiving the certain frequency is called a short wave radio station. Because short wave communication mainly relies on back and forth reflection and refraction between the ionosphere and the ground for propagation, short waves can propagate far no matter in daytime or at night. With the continuous development of computer, microelectronics and wireless communication technologies, the short-wave communication technology has breakthrough progress, and the short-wave communication technology is used for emergency communication and disaster-resistant communication, and particularly plays an important and wide role in the aspect of land, sea and air unified communication command required by military.
Currently, short-wave radio stations have become digitalized, and the working frequency band of the short-wave radio stations is not limited to the original short-wave frequency band range, so that the short-wave radio stations have the characteristics of multiple bands and multiple channels. In the digitizing process of the short-wave radio station, the performance of the short-wave power amplifier is required to be higher and higher, namely, the power amplifier has higher efficiency on the premise of meeting higher linearity requirements. To meet this requirement, the amplifier is made to be both linear and efficient, various means are adopted to achieve high efficiency and high linearity of the amplifier, and the predistortion technology can also utilize the adaptive principle to track and compensate the error of the power amplifier caused by the change of environmental factors such as temperature, humidity and the like. In a word, the predistortion technology not only can improve the efficiency of the transmitter and reduce the cost and the volume, but also can effectively increase the linearity of the transmitter to improve the system efficiency and the communication quality, and has great practical significance for the development and the realization of future high-efficiency short-wave radio stations.
For ease of implementation, short wave digital predistortion architecture employs a look-up Table (LUT) method, as shown in fig. 1, by means of which predistortion of this form works. The LUT is retrieved as a function of the amplitude of the signal, or the input amplitude, and then corrects the amplitude and phase of the signal applied to the power amplifier input. The predistortion system architecture includes two channels: a loop path for data training and a predistortion path. The data training channel is a loop structure, the core part of the data training channel is a predistortion algorithm module, and the module processes a feedback signal (obtained by power amplifier output coupling) after power amplification and an original input signal to extract the distortion characteristic of the power amplifier, thereby obtaining the LUT parameter of the distortion inverse characteristic of the power amplifier. When the power amplifier characteristic changes along with the change of time or external environment, the predistortion anti-characteristic LUT parameter can be updated through the adaptive predistortion algorithm.
FIG. 2 is a schematic diagram of a basic predistortion parameter extraction method; where X is the input signal and Y is the output (feedback) signal. Firstly, determining output power according to requirements, when the system runs for the first time, directly connecting the system, obtaining a first group of output data after LUT parameters are all 1 and LUT (X) =X and signals pass through a power amplifier, feeding back to carry out predistortion, obtaining the LUT parameters by solving min|LUT (Y) -LUT (X) |, and continuously correcting the LUT parameters in an iterative process until the output signal Y obtains a satisfactory effect, and extracting the LUT parameters at the moment to be used as a predistortion LUT under the output power for an open loop system.
In the prior predistortion model table index, an integer index structure is adopted, namely, equal interval integer division is adopted for the amplitude of an input signal, so that a table index address is obtained. Therefore, the index accuracy of the table is reduced, and under the condition of limited table length, the predistortion effect of the table can generate 'burrs' in the time domain due to the problem of insufficient index accuracy, which can further influence other short-wave system processing links such as capturing and the like.
Disclosure of Invention
Aiming at the problems existing in the prior art, the application aims to provide a method for improving the index precision of short-wave digital predistortion processing, and the method designs a memory polynomial model of an integral small-number table index structure aiming at the characteristics of a short-wave signal system, so that the index precision of a table can be improved, and the performance of a digital predistortion system is improved; and compared with an integer table index structure, the third-order intermodulation index (2 MHz-30 MHz) of the transmitter can be improved by more than 5 dB.
In order to achieve the above purpose, the present application is realized by the following technical scheme.
A method for improving index accuracy of short wave digital predistortion processing comprises the following steps:
step 1, amplifying an original input signal X (n) by a power amplifier to obtain an amplified output signal; wherein the input signal of the power amplifier is equal to the original input signal X (n);
step 2, processing the original input signal and the output signal after the power amplification by adopting a predistortion algorithm module, and extracting the distortion characteristic of the power amplification so as to obtain the lookup table parameters of the distortion inverse characteristic of the power amplification; the predistortion algorithm module adopts an integer integrated small-scale table index structure;
and 3, correcting the amplitude and the phase of the input signal X (n) of the power amplifier by adopting the lookup table parameters of the power amplifier distortion inverse characteristic to obtain a signal Z (n) after predistortion treatment.
Further, the step 2 specifically includes the following sub-steps:
2.1, performing power calculation on an original input signal X (n) to obtain the power of the input signal;
and 2.2, dividing the power of the input signal into an integer part and a decimal part, and determining the corresponding lookup table parameters of the power amplifier distortion inverse characteristic according to the integer part and the decimal part.
Further, substep 2.2 is specifically:
determining a corresponding table index address according to the integer part, and then finding a corresponding integer part lookup table and an adjacent lookup table which is 1 larger than the integer part according to the table index address, so that integer part lookup table parameters corresponding to the integer part lookup table and adjacent lookup table parameters corresponding to the adjacent lookup table can be obtained;
performing difference operation on the adjacent lookup table parameters and the integer part lookup table parameters to obtain a difference value; multiplying the difference value with the decimal part to obtain a product;
and carrying out summation operation on the product and the integer part lookup table parameters to obtain a summation value, wherein the summation value is the lookup table parameter of the power amplifier distortion inverse characteristic.
Further, in the substep 2.2, the calculation formula of the lookup table parameter of the power amplifier distortion inverse characteristic is:
LUT(A.B)=LUT(A)+(A.B-A)*(LUT(A+1)-LUT(A))
wherein a.b is the power of the input signal, a is the integer part of the power of the input signal, B is the fractional part of the power of the input signal, LUT (a) is the integer part lookup table parameter, (LUT (a+1) is the adjacent lookup table parameter, and LUT (a.b) is the lookup table parameter of the power amplifier distortion inverse characteristic.
Further, the step 3 specifically includes:
and multiplying the sum value with the input signal X (n) of the power amplifier to obtain a signal Z (n) after predistortion treatment.
Compared with the prior art, the application has the beneficial effects that:
because the characteristics of the power amplifier are different along with the characteristics of the input signals, the application designs a memory polynomial model of an integral small-number table index structure aiming at the short-wave power amplifier digital predistortion linearization according to the factors such as actual predistortion linearization performance, realization complexity and the like, thereby improving the table index precision and further improving the performance of a digital predistortion system; and compared with the integer form index structure, the integer integrated small form index structure can eliminate burrs generated by the integer index structure in the time domain, and meanwhile, the third-order intermodulation index (2 MHz-30 MHz) of the transmitter can be improved by more than 5 dB.
Drawings
The application will now be described in further detail with reference to the drawings and to specific examples.
FIG. 1 is a diagram of the overall architecture of a basic predistortion architecture;
FIG. 2 is a diagram illustrating the extraction of a pre-distortion table according to the prior art;
FIG. 3 is a diagram of a prior art predistortion model table index block;
FIG. 4 is a diagram of the prior art predistortion model integer table index time domain simulation result; wherein, (a) is an overall simulation result diagram; (b) is a partial enlarged view at a in fig. (a);
FIG. 5 is a diagram of an integral integer form index structure of a predistortion model in accordance with the present application;
FIG. 6 is a graph of the results of the predistortion model integer integrated wavelet table index time domain simulation of the present application; (a) is a whole simulation result graph; (b) is a partial enlarged view at a in fig. (a);
FIG. 7 is a graph comparing the frequency domain simulation results of the predistortion model integer index and the integer integral fraction index; wherein; (a) is a comparison graph of the overall simulation result; (b) is a partial enlarged view at A in the drawing (a).
Detailed Description
Embodiments of the present application will be described in detail below with reference to examples, but it will be understood by those skilled in the art that the following examples are only for illustrating the present application and should not be construed as limiting the scope of the present application.
Example 1
A method for improving index accuracy of short wave digital predistortion processing comprises the following steps:
step 1, amplifying an original input signal X (n) by a power amplifier to obtain an amplified output signal; wherein the input signal of the power amplifier is equal to the original input signal X (n).
Step 2, processing the original input signal and the output signal after the power amplification by adopting a predistortion algorithm module, and extracting the distortion characteristic of the power amplification so as to obtain the lookup table parameters of the distortion inverse characteristic of the power amplification; the predistortion algorithm module adopts an integer integrated small-scale table index structure.
Specifically, step 2 comprises the following sub-steps:
2.1, performing power calculation on an original input signal X (n) to obtain the power of the input signal;
sub-step 2.2, dividing the power of the input signal into an integer part and a decimal part; and determining the lookup table parameters of the corresponding power amplifier distortion inverse characteristic according to the integer part and the decimal part, wherein the lookup table parameters are specifically as follows:
determining a corresponding table index address according to the integer part, and then finding a corresponding integer part lookup table and an adjacent lookup table which is 1 larger than the integer part according to the table index address, so that integer part lookup table parameters corresponding to the integer part lookup table and adjacent lookup table parameters corresponding to the adjacent lookup table can be obtained;
performing difference operation on the adjacent lookup table parameters and the integer part lookup table parameters to obtain a difference value; multiplying the difference value with the decimal part to obtain a product;
and carrying out summation operation on the product and the integer part lookup table parameter to obtain a summation value, wherein the summation value is the lookup table parameter of the power amplifier distortion inverse characteristic, and the specific calculation formula is as follows:
LUT(A.B)=LUT(A)+(A.B-A)*(LUT(A+1)-LUT(A))
wherein a.b is the power of the input signal, a is the integer part of the power of the input signal, B is the fractional part of the power of the input signal, LUT (a) is the integer part lookup table parameter, (LUT (a+1) is the adjacent lookup table parameter, and LUT (a.b) is the lookup table parameter of the power amplifier distortion inverse characteristic.
And 3, correcting the amplitude and the phase of the input signal X (n) of the power amplifier by adopting the lookup table parameters of the power amplifier distortion inverse characteristic to obtain a signal Z (n) after predistortion treatment.
The method comprises the following steps: and multiplying the sum value with the input signal X (n) of the power amplifier to obtain a signal Z (n) after predistortion treatment.
The predistortion model integer integrated small table index structure of the application as shown in fig. 5, when the calculated power of the input signal is 36.6, adopts the integer integrated small table index structure, and obtains the lookup table parameter LUT (36.6) of the inverse characteristic of the power amplifier distortion by the following formula:
LUT(36.6)=LUT(36)+(36.6-36)*(LUT(37)-LUT(36))
and then multiplying the input signal X (n) of the power amplifier by using the lookup table parameter LUT (36.6) of the obtained power amplifier distortion inverse characteristic to obtain a signal Z (n) after predistortion treatment.
FIG. 6 is a graph of the results of a predistortion model integer integrated wavelet table index time domain simulation; as shown in fig. 6, when the integer-integrated wavelet table index is adopted, the predistortion effect is improved due to the improvement of the index precision, and the "burr" generated by the integer index structure is eliminated in the time domain.
Comparative example 1
Fig. 3 is a diagram showing a table index structure of a conventional predistortion model. The input signal X (n) is subjected to power calculation to obtain a corresponding table index address P (X (n)), then a corresponding LUT parameter, namely LUT (P (X (n)), is obtained according to the table index address P (X (n)), and then the amplitude and the phase of the signal X (n) applied to the input end of the power amplifier are corrected by adopting the LUT (P (X (n)), so that the Z (n) after predistortion processing is obtained.
In the table index, an integer index structure is adopted, that is, equal interval integer division is adopted for the amplitude of the input signal, so that a table index address is obtained. As shown in fig. 4, when the calculated power is 36.6, the method of rounding is generally adopted to directly round 36.6 into 37, and then the table index lookup table parameters corresponding to 37 are adopted to correct the amplitude and phase of the signal X (n) applied to the input end of the power amplifier. As can be seen from fig. 4, the index accuracy of the table is reduced by using the integer table index, and under the condition of limited table length, the predistortion effect can generate a "burr" phenomenon in the time domain due to the problem of insufficient index accuracy, which further affects other short wave system processing links such as capturing.
The predistortion model integer-integrated decimal index structure of example 1 and the integer index structure of comparative example 1 were simulated in the frequency domain, and the results are shown in fig. 7. As can be seen from fig. 7, after the integral decimal index structure is adopted, the third-order intermodulation index (2 MHz-30 MHz) of the transmitter can be improved by more than 5dB, and the spectrum quality of the transmitter is also greatly improved.
While the application has been described in detail in this specification with reference to the general description and the specific embodiments thereof, it will be apparent to one skilled in the art that modifications and improvements can be made thereto. Accordingly, such modifications or improvements may be made without departing from the spirit of the application and are intended to be within the scope of the application as claimed.

Claims (1)

1. The method for improving the index precision of the short wave digital predistortion processing is characterized by comprising the following steps of:
step 1, amplifying an original input signal X (n) by a power amplifier to obtain an amplified output signal; wherein the input signal of the power amplifier is equal to the original input signal X (n);
step 2, processing the original input signal X (n) and the output signal after the power amplification by adopting a predistortion algorithm module, and extracting the distortion characteristic of the power amplification so as to obtain the lookup table parameters of the distortion inverse characteristic of the power amplification; the predistortion algorithm module adopts an integer integrated small-scale table index structure;
the method specifically comprises the following substeps:
2.1, performing power calculation on an original input signal X (n) to obtain the power of the input signal;
sub-step 2.2, dividing the power of the input signal into an integer part and a decimal part, and determining the corresponding lookup table parameters of the power amplifier distortion inverse characteristic according to the integer part and the decimal part;
the method comprises the following steps:
determining a corresponding table index address according to the integer part, and then finding a corresponding integer part lookup table and an adjacent lookup table which is 1 larger than the integer part according to the table index address, so that integer part lookup table parameters corresponding to the integer part lookup table and adjacent lookup table parameters corresponding to the adjacent lookup table can be obtained;
performing difference operation on the adjacent lookup table parameters and the integer part lookup table parameters to obtain a difference value; multiplying the difference value with the decimal part to obtain a product;
carrying out summation operation on the product and the integer part lookup table parameters to obtain a summation value, wherein the summation value is the lookup table parameter of the power amplifier distortion inverse characteristic;
the calculation formula of the lookup table parameters of the power amplifier distortion inverse characteristic is as follows:
LUT(A.B)=LUT(A)+(A.B-A)*(LUT(A+1)-LUT(A))
wherein, A.B is the power of the input signal, A is the integer part of the power of the input signal, B is the decimal part of the power of the input signal, LUT (A) is the integer part lookup table parameter, (LUT (A+1) is the adjacent lookup table parameter, LUT (A.B) is the lookup table parameter of the power amplifier distortion inverse characteristic;
step 3, correcting the amplitude and phase of the input signal X (n) of the power amplifier by adopting the lookup table parameters of the power amplifier distortion inverse characteristic to obtain a signal Z (n) after predistortion treatment;
the method comprises the following steps:
and multiplying the sum value with the input signal X (n) of the power amplifier to obtain a signal Z (n) after predistortion treatment.
CN202010534060.7A 2020-06-12 2020-06-12 Method for improving index accuracy of short wave digital predistortion processing Active CN111669132B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010534060.7A CN111669132B (en) 2020-06-12 2020-06-12 Method for improving index accuracy of short wave digital predistortion processing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010534060.7A CN111669132B (en) 2020-06-12 2020-06-12 Method for improving index accuracy of short wave digital predistortion processing

Publications (2)

Publication Number Publication Date
CN111669132A CN111669132A (en) 2020-09-15
CN111669132B true CN111669132B (en) 2023-11-07

Family

ID=72387180

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010534060.7A Active CN111669132B (en) 2020-06-12 2020-06-12 Method for improving index accuracy of short wave digital predistortion processing

Country Status (1)

Country Link
CN (1) CN111669132B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1610249A (en) * 2003-10-22 2005-04-27 安德鲁公司 Amplifier pre-distortion processing based on composite look-up tables
CN1689295A (en) * 2002-10-31 2005-10-26 中兴通讯股份有限公司 Broadband predistortion linearization method and system
CN101459636A (en) * 2007-12-12 2009-06-17 中兴通讯股份有限公司 Adaptive pre-distortion method
CN102404263A (en) * 2010-09-19 2012-04-04 中兴通讯股份有限公司 Method and device for processing digital pre-distortion
CN108920748A (en) * 2018-05-22 2018-11-30 复旦大学 A kind of non-linear arithmetic circuit structure with high flexibility and precision

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1689295A (en) * 2002-10-31 2005-10-26 中兴通讯股份有限公司 Broadband predistortion linearization method and system
CN1610249A (en) * 2003-10-22 2005-04-27 安德鲁公司 Amplifier pre-distortion processing based on composite look-up tables
CN101459636A (en) * 2007-12-12 2009-06-17 中兴通讯股份有限公司 Adaptive pre-distortion method
CN102404263A (en) * 2010-09-19 2012-04-04 中兴通讯股份有限公司 Method and device for processing digital pre-distortion
CN108920748A (en) * 2018-05-22 2018-11-30 复旦大学 A kind of non-linear arithmetic circuit structure with high flexibility and precision

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Digital Predistorter for Short-Wave Power Amplifier with Improving Index Accuracy of Lookup Table Based on FPGA;Ren, Jijun;《IEEE ACCESS 》;全文 *
一种有效的自适应预失真算法;武占生,葛建华,韩乔;无线通信技术(01);全文 *
新型查找表数字预失真器的设计与实现;席重;《中国优秀硕士学位论文全文数据库 (信息科技辑)》;全文 *

Also Published As

Publication number Publication date
CN111669132A (en) 2020-09-15

Similar Documents

Publication Publication Date Title
CN101175061B (en) Self-adapting digital predistortion method and apparatus for OFDM transmitter
US10075324B2 (en) Predistortion processing apparatus and method
CN102143107B (en) Method and device for realizing pre-distortion of digital baseband
US8737937B2 (en) Distortion compensation apparatus, transmitter, and distortion compensation method
CN102075469B (en) Estimation method for signal delay time of digital pre-distortion system
CN111064439B (en) System and method for improving short wave digital predistortion performance
CN107359864B (en) Self-adaptive agile digital predistortion method for frequency agile power amplifier
CN107592083A (en) The digital pre-distortion system and its method of millimeter wave broadband power amplifier
CN101860326A (en) Indirect learning predistortion linearized system based on Volterra series
CN103296978B (en) A kind of digital pre-distortion method and device
CN107659273B (en) Reconfigurable digital predistortion system and method for millimeter wave broadband power amplifier
CN102111361A (en) Self-adaptive predistorter design based on table look-up method of amplifier estimator
CN114400979A (en) System and method for improving digital predistortion of short wave interaction and harmonic distortion
CN105262447A (en) Pre-distortion method and device for power amplifier and radio frequency system
CN103929136A (en) Broadband mobile communication LTE system transmitter power amplifier predistortion processing method
CN111669132B (en) Method for improving index accuracy of short wave digital predistortion processing
CN111082829B (en) Signal processing device, radio remote unit and base station
CN102904846B (en) A kind of digital pre-distortion processing method adapting to fast changed signal
CN110086438B (en) Digital predistortion system and method for passive multi-beam transmitter
Li et al. Iterative learning control assisted neural network for digital predistortion of MIMO power amplifier
CN106850493B (en) A kind of pouring-in peak-to-average power ratio of composite non-linear inhibits and digital pre-distortion method
CN111901264B (en) Short wave digital predistortion method based on non-uniform quantization lookup table
CN108134584B (en) In-band and out-of-band combined digital predistortion system and method for broadband radio frequency power amplifier
CN107592084B (en) Distortion compensation device of broadband power amplifier based on dynamic modeling
Long et al. Indirect learning hybrid memory predistorter based on polynomial and look-up-table

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