CN102238126B - Method for reducing peak-to-average power ratio of OFDM (orthogonal frequency division multiplexing)/OQAM system based on selective sequence - Google Patents

Method for reducing peak-to-average power ratio of OFDM (orthogonal frequency division multiplexing)/OQAM system based on selective sequence Download PDF

Info

Publication number
CN102238126B
CN102238126B CN 201110169363 CN201110169363A CN102238126B CN 102238126 B CN102238126 B CN 102238126B CN 201110169363 CN201110169363 CN 201110169363 CN 201110169363 A CN201110169363 A CN 201110169363A CN 102238126 B CN102238126 B CN 102238126B
Authority
CN
China
Prior art keywords
sequence
ofdm
prime
signal
phase place
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
CN 201110169363
Other languages
Chinese (zh)
Other versions
CN102238126A (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.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and 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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN 201110169363 priority Critical patent/CN102238126B/en
Publication of CN102238126A publication Critical patent/CN102238126A/en
Application granted granted Critical
Publication of CN102238126B publication Critical patent/CN102238126B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a method for reducing a peak-to-average power ratio of an OFDM (orthogonal frequency division multiplexing)/OQAM (orthogonally multiplexed QAM) system based on a selective sequence. The method comprises the following steps of: coding, interleaving and modulating an input data stream, performing up-sampling on the obtained information source data packets, and then, enablingthe information source data packets to respectively pass through a filter bank; respectively selecting an optimized phase rotation sequence for a time domain signal of each obtained data block in order; and finally, selecting a signal, of which the peak-to-average power ratio is lowest, for transmission, and transmitting information of the phase rotation sequence as sideband information. The method provided by the invention can effectively lower the peak-to-average power ratio of the OFDM/OQAM system without influencing the power spectral density and the error rate performance of the system, improve work efficiency an HPA (high-power amplifier) and reduce energy consumption.

Description

Reduce the method for the peak-to-average power ratio of OFDM/OQAM system based on the selectivity sequence
Technical field
The invention belongs to the wireless and cable communicating technology field based on OFDM (OFDM/OQAM) system of offset quadrature amplitude modulation(PAM), be specifically related to adopt the selectivity sequence to reduce the method for the peak-to-average power ratio of OFDM/OQAM system.
Background technology
OFDM (OFDM/OQAM) technology based on the offset quadrature amplitude modulation(PAM) is a kind of promising multi-carrier modulation technology, because it has advantages such as high spectrum resolution and Large Volume Data transfer rate, has obtained extensive concern.Compare with traditional ofdm system, the part technical advantage of OFDM/OQAM system may be summarized as follows: 1) when the user is asynchronous, can carry out the distribution of subcarrier well; 2) meet the prototype filter pulse that condition is reinvented in (being similar to) perfection by design, under the prerequisite that does not need Cyclic Prefix, still can obtain littler intersymbol interference and interchannel interference; 3) by being combined with the OQAM modulated filter bank, can obtain maximized message transmission rate; 4) can finish high performance frequency spectrum perception; 5) can make the secondary lobe of frequency spectrum lower by the original shape Filter Design, the interference of the system that makes between can successive bands is lower; 6) frequency spectrum perception and transfer of data can be carried out synchronously.Therefore, for LTE-advanced and IMT-advanced etc., OFDM/OQAM is suggested as a good candidate scheme.For example, insert and cognitive radio for dynamic spectrum, OFDM/OQAM has been chosen as the up-and-coming technology of physical layer, wherein frequency spectrum perception and transfer of data can be carried out simultaneously, that is to say, same bank of filters not only can be used for frequency spectrum perception also can be used for transfer of data, and this has guaranteed the compatible performance of system.This is a key facility for opportunistic communication system efficiently.
But in the OFDM/OQAM system, some challenging problems still do not solve.Wherein one of problem of most critical is exactly the design of its high power amplifier (HPA), because the peak-to-average power ratio of OFDM/OQAM signal (PAPR) is too big.Generally, HPA is operated in the zone of saturation or near its power output efficient maximum of zone of saturation.Therefore, because the high PAPR of input signal, memoryless nonlinear distortion will be introduced in communication channel, cause out of band spectrum regeneration.Therefore, in the OFDM/OQAM wireless communication system, the high PAPR that reduces signal is extremely important.
People have proposed the PAPR that a lot of technology reduce traditional ofdm system.For example: coding techniques, multitone is reserved, companding technology etc.But all methods that reduce the PAPR of ofdm signal all can not directly apply in the OFDM/OQAM system.Because in traditional ofdm system, each OFDM symbol is independently in time domain, does not exist overlapping between the adjacent OFDM symbol.Therefore, the reduction of the PAPR of symbol before or after the reduction of the PAPR of an OFDM symbol can not influence.But in the OFDM/OQAM system, time-domain symbol exists overlapping.Therefore, the PAPR that we reduce each OFDM/OQAM symbol independently can impact other symbols, thereby the effect that causes final PAPR to reduce is bad.
Summary of the invention
Higher at signal peak-to-average power power ratio in the OFDM/OQAM system, and the method for existing control and reduction ofdm signal peak-to-average power ratio can not directly apply in the OFDM/OQAM system, the present invention proposes the method that reduces the peak-to-average power ratio of OFDM/OQAM system based on the selectivity sequence, this method is in the peak-to-average power ratio that reduces the OFDM/OQAM system, power spectral density and the bit error rate performance that can guarantee system are unaffected, and the operating efficiency of raising HPA, save power consumption.
Provided by the inventionly a kind ofly reduce the method for the peak-to-average power ratio of OFDM/OQAM system based on the selectivity sequence, the information source packet of establishing an OFDM/OQAM system comprises M data block, and each data block comprises N q system symbol; Wherein, q=2 l, l, M and N are positive integer, it is characterized in that, and this method comprises the steps:
(1) data flow of input obtains an information source packet X=[X through after encoding, interweaving, modulate 0, X 1..., X M-1], X wherein mBe m data block among the X, be expressed as X m=[X m(0), X m(1) ..., X m(N-1)] T, X wherein m(n) n data symbol of m data block of expression, and X m(n)=a m(n)+jb m(n), n=0,1 ..., N-1, wherein a m(n) and b m(n) be X respectively m(n) real part and imaginary part, j is the imaginary part of symbol, the transposition of T representing matrix;
(2) the input data symbol is carried out N times of up-sampling, X mIn the up-sampling that obtains of n data symbol after signal
Figure BDA0000070294050000031
Be designated as x m n ( t ) = X m ( n ) δ ( t - mT ) , Wherein T is a data symbol duration, and t is time variable, t 〉=0, impulse function δ ( t ) = 1 , t = 0 0 , else , Then the signal that obtains is passed through prototype filter h (t), and will
Figure BDA0000070294050000034
Imaginary part time-delay
Figure BDA0000070294050000035
X then mIn n data symbol by filter after and be modulated to mark signal on the carrier wave
Figure BDA0000070294050000036
Be expressed as
s m n ( t ) = { a m ( n ) h ( t - mT ) + j b m ( n ) h ( t - T 2 - mT ) } e jn ( 2 π T t + T 2 )
Wherein e jn ( 2 π T t + T 2 ) = sin ( n ( 2 π T t + T 2 ) ) + j cos ( n ( 2 π T t + T 2 ) ) ;
(3) adopt the selectivity sequence to reduce peak-to-average power ratio:
(3.1) produce U phase place rotatable sequence, wherein each phase rotation coefficient is selected from and defaultly gathers that { 1 ,-1} remembers that U phase place rotatable sequence is
Figure BDA0000070294050000039
U=0,1 ..., U-1, wherein
Figure BDA00000702940500000310
N=0,1 ..., N-1, note Expression acts on the phase place rotatable sequence of m data block, and makes P M, u=P u
(3.2) initialization m=0 multiply by U phase place rotatable sequence respectively to m data block and obtains U candidate signal
s ~ 0 , u ( t ) = Σ n = 0 N - 1 s 0 n ( t ) P 0 , u n , u = 0,1 , . . . , U - 1
And calculate this U candidate signal
Figure BDA00000702940500000313
PAPR 0 u = 10 log 10 max 0 ≤ t ≤ ( K - 1 ) T [ | s ~ 0 , u ( t ) | 2 ] E [ | s ~ 0 , u ( t ) | 2 ] , u = 0,1 , . . . , U - 1
E[ wherein] expression averages;
(3.3) note is minimum
Figure BDA00000702940500000315
For
Figure BDA00000702940500000316
Its corresponding block is
Figure BDA00000702940500000317
Corresponding phase place rotatable sequence sequence number note is S 0, and make m=m+1;
(3.4) m data block be multiply by U phase place rotatable sequence respectively and obtain U candidate signal
s ~ m , u ( t ) = Σ n = 0 N - 1 s m n ( t ) P m , u n , u = 0,1 , . . . , U - 1
And calculate this U candidate signal
PAPR m u = 10 log 10 max mT ≤ t ≤ ( m + K - 1 ) T [ | Σ p = 0 m - 1 s ~ p * ( t ) + s ~ m , u ( t ) | 2 ] E [ | Σ p = 0 m - 1 s ~ p * ( t ) + s ~ m , u ( t ) | 2 ] , u = 0,1 , . . . , U - 1 ;
(3.5) note is minimum
Figure BDA0000070294050000044
For
Figure BDA0000070294050000045
Its corresponding block is
Figure BDA0000070294050000046
Corresponding phase place rotatable sequence sequence number note is S m, and make m=m+1;
(3.6) if m<M forwards step (3.4) to; Otherwise, calculate final output signal
Figure BDA0000070294050000047
s ~ ( t ) = Σ m = 0 M - 1 s ~ m * ( t )
Forward step (4) to;
(4) the phase place rotatable sequence sequence number of establishing optimization is S=[S 0, S 1..., S M-1], with the signal of output
Figure BDA0000070294050000049
Send into channel, and sequence S is sent into channel as side information.
The present invention is a kind of method that reduces the peak-to-average power ratio of OFDM/OQAM signal based on the selectivity sequence, it not only can reduce the peak-to-average power ratio of OFDM/OQAM signal preferably, and can also keep the power spectral density of OFDM/OQAM signal constant, thereby obtain bit error rate performance preferably., through after encoding, interweaving, modulate the information source packet that obtains is carried out passing through bank of filters respectively after the up-sampling in data.Then, with the phase place rotatable sequence that the OFDM/OQAM signal that obtains is optimized according to the rule selection, finally select the signal of peak-to-average power ratio minimum to transmit, and the information of phase place rotatable sequence is sent to receiving terminal as side information.The present invention can reduce the peak-to-average power ratio of OFDM/OQAM signal effectively, simultaneously, can guarantee that the power spectral density of system and bit error rate performance are unaffected, improves the operating efficiency of HPA, saves power consumption.
Description of drawings
Fig. 1 is the system flow block diagram of the inventive method;
Fig. 2 adopts the selectivity sequence to reduce the FB(flow block) of peak-to-average power ratio among the present invention.
Embodiment
The present invention is further described below in conjunction with accompanying drawing and one group of concrete parameter:
Employing selectivity sequence provided by the invention reduce the OFDM/OQAM system peak-to-average power ratio method handling process as shown in Figure 1, wherein, adopt handling process that the selectivity sequence reduces peak-to-average power ratio as shown in Figure 2.
The duration of prototype filter h (t) is KT, and wherein K is a positive even numbers; General, K is more big, and secondary lobe is more little; But the K increase can cause higher computation complexity; Therefore, the K value of choosing in the reality is no more than 8 usually.
The step of the inventive method (1) and (2) are identical with the OFDM/OQAM system original treatment step of transmitting terminal.
Below by by embodiment the present invention being described in further detail, but following examples only are illustrative, and protection scope of the present invention is not subjected to the restriction of these embodiment.
Example:
Parameter declaration: adopt the QAM modulation of q=4, sub-carrier number N=64, the OFDM/OQAM system of prototype filter parameter K=4, the phase place rotatable sequence number that produces in the selectivity sequence is got U=16 and U=128 respectively.
Simulation result shows that the ability that the present invention reduces PAPR is better than the ability that traditional selectivity mapping method directly applies to the reduction PAPR of OFDM/OQAM system far away.Simultaneously, the present invention can not influence the power spectral density of OFDM/OQAM signal.And, owing to the reduction of PAPR, improved the operating efficiency of HPA, saved consumed power.
When U=128, at Pr{PAPR>PAPR 0}=10 -3The time, the reduction amount that the selectivity sequence method that adopts the present invention to propose can obtain PAPR is 4.7dB; The reduction amount that adopts original selectivity mapping method can obtain PAPR is 1.3dB.Obviously, when parameter was identical, the present invention had improved the ability that reduces PAPR greatly.Simultaneously, the present invention's power spectral density of reducing the OFDM/OQAM signal after the PAPR does not overlap fully with the ideal power spectrum density that the OFDM/OQAM primary signal is passed through under the HPA situation.When U=16, to compare with primary signal, the method for proposition has been saved 45% energy, and the energy that traditional selectivity mapping method is saved only is 10%.Simultaneously, the ability that the present invention reduces PAPR is better than the ability that traditional selectivity mapping method directly applies to the reduction PAPR of OFDM/OQAM system far away, has promoted the ability of saving energy simultaneously greatly.
The above is preferred embodiment of the present invention, but the present invention should not be confined to the disclosed content of this embodiment and accompanying drawing.So everyly do not break away from the equivalence of finishing under the spirit disclosed in this invention or revise, all fall into the scope of protection of the invention.

Claims (1)

1. one kind is reduced the method for the peak-to-average power ratio of OFDM/OQAM system based on the selectivity sequence, and the information source packet of establishing an OFDM/OQAM system comprises M data block, and each data block comprises N q system symbol; Wherein, q=2 l, l, M and N are positive integer, it is characterized in that, and this method comprises the steps:
(1) data flow of input obtains an information source packet X=[X through after encoding, interweaving, modulate 0, X 1..., X M-1], X wherein mBe m data block among the X, be expressed as X m=[X m(0), X m(1) ..., X m(N-1)] T, X wherein m(n) n data symbol of m data block of expression, and X m(n)=a m(n)+jb m(n), n=0,1 ..., N-1, wherein a m(n) and b m(n) be X respectively m(n) real part and imaginary part, j is the imaginary part of symbol, TThe transposition of representing matrix;
(2) the input data symbol is carried out N times of up-sampling, X mIn the up-sampling that obtains of n data symbol after signal
Figure FDA00002911320400011
Be designated as Wherein T ' is a data symbol duration, and t is time variable, t 〉=0, impulse function δ ( t ) = 1 , t = 0 0 , else , Then the signal behind the up-sampling that obtains is passed through prototype filter h (t), and will Imaginary part time-delay
Figure FDA00002911320400015
X then mIn n data symbol by filter after and be modulated to mark signal on the carrier wave
Figure FDA00002911320400016
Be expressed as
s m n ( t ) = { a m ( n ) h ( t - m T ′ ) + j b m ( n ) h ( t - T ′ 2 - m T ′ ) } e jn ( 2 π T ′ t + T ′ 2 )
Wherein, e jn ( 2 π T ′ t + T ′ 2 ) = sin ( n ( 2 π T ′ t + T ′ 2 ) ) + j cos ( n ( 2 π T ′ t + T ′ 2 ) ) ;
(3) adopt the selectivity sequence to reduce peak-to-average power ratio:
(3.1) produce U phase place rotatable sequence, wherein each phase rotation coefficient is selected from and defaultly gathers that { 1 ,-1} remembers that U phase place rotatable sequence is
Figure FDA00002911320400019
U=0,1 ..., U-1, wherein , n=0,1 ..., N-1, note
Figure FDA000029113204000111
Expression acts on the phase place rotatable sequence of m data block, and makes P M, u=P u
(3.2) initialization m=0 is to m data block X mIn n data symbol by filter after and be modulated to mark signal on the carrier wave
Figure FDA00002911320400021
Multiply by U phase place rotatable sequence respectively and obtain U candidate signal
Figure FDA00002911320400022
s ~ m , u ( t ) = Σ n = 0 N - 1 s m n ( t ) P m , u n , u = 0,1 , . . . , U - 1
And calculate this U candidate signal
PAPR m u = 10 log 10 max 0 ≤ t ≤ ( K - 1 ) T [ | s ~ m , u ( t ) | 2 ] E [ | s ~ m , u ( t ) | 2 ] , u = 0,1 , . . . , U - 1
E[wherein] expression averages;
(3.3) obtain in the note step (3.2)
Figure FDA00002911320400026
Middle minimum value is designated as
Figure FDA00002911320400027
, its corresponding block is
Figure FDA00002911320400028
, corresponding phase place rotatable sequence sequence number note is S 0, and make m=m+1;
(3.4) m data block be multiply by U phase place rotatable sequence respectively and obtain U candidate signal
s ~ m , u ( t ) = Σ n = 0 N - 1 s m n ( t ) P m , u n , u = 0,1 , . . . , U - 1
And calculate this U candidate signal
Figure FDA000029113204000210
PAPR m u = 10 log 10 max mT ≤ t ≤ ( m + K - 1 ) T [ | Σ p = 0 m - 1 s ~ p * ( t ) + s ~ m , u ( t ) | 2 ] E [ | Σ p = 0 m - 1 s ~ p * ( t ) s ~ m , u ( t ) | 2 ] , u = 0,1 , . . . , U - 1 ;
(3.5) note is minimum
Figure FDA000029113204000212
For
Figure FDA000029113204000213
, its corresponding block is
Figure FDA000029113204000214
, corresponding phase place rotatable sequence sequence number note is S m, and make m=m+1;
(3.6) if m<M forwards step (3.4) to; Otherwise, calculate final output signal
Figure FDA000029113204000215
s ~ ( t ) = Σ m = 0 M - 1 s ~ m * ( t )
Forward step (4) to;
(4) the phase place rotatable sequence sequence number of establishing optimization is S=[S 0, S 1..., S M-1], with the signal of output
Figure FDA000029113204000217
Send into channel, and sequence S is sent into channel as side information.
CN 201110169363 2011-06-22 2011-06-22 Method for reducing peak-to-average power ratio of OFDM (orthogonal frequency division multiplexing)/OQAM system based on selective sequence Active CN102238126B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110169363 CN102238126B (en) 2011-06-22 2011-06-22 Method for reducing peak-to-average power ratio of OFDM (orthogonal frequency division multiplexing)/OQAM system based on selective sequence

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110169363 CN102238126B (en) 2011-06-22 2011-06-22 Method for reducing peak-to-average power ratio of OFDM (orthogonal frequency division multiplexing)/OQAM system based on selective sequence

Publications (2)

Publication Number Publication Date
CN102238126A CN102238126A (en) 2011-11-09
CN102238126B true CN102238126B (en) 2013-09-18

Family

ID=44888350

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110169363 Active CN102238126B (en) 2011-06-22 2011-06-22 Method for reducing peak-to-average power ratio of OFDM (orthogonal frequency division multiplexing)/OQAM system based on selective sequence

Country Status (1)

Country Link
CN (1) CN102238126B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103973395A (en) * 2013-01-28 2014-08-06 华为技术有限公司 Data processing method, device and system in wireless communication
US9419770B2 (en) * 2014-03-31 2016-08-16 Huawei Technologies Co., Ltd. Method and apparatus for asynchronous OFDMA/SC-FDMA
CN106572042B (en) * 2015-10-09 2020-03-31 华为技术有限公司 Method and device for transmitting data
CN106789828B (en) * 2016-12-26 2019-10-08 中南大学 A method of being fed back based on peak time tracking reduces FBMC-OQAM system peak-to-average power ratio
CN108270708B (en) 2016-12-30 2020-07-07 维沃移动通信有限公司 Information transmission method, sending end and receiving end for reducing PAPR
CN108322295B (en) * 2017-01-17 2021-12-24 维沃移动通信有限公司 Sending method, receiving method, sending end and receiving end of sideband information
CN106850496B (en) * 2017-03-02 2019-08-20 北京邮电大学 A kind of quadrature amplitude modulation signal phase recovery method and device
CN109194603B (en) * 2018-11-21 2019-06-14 济南大学 Reduce the guiding deinterleaving method of MIMO-OFDM system peak-to-average power ratio

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101778071A (en) * 2010-01-29 2010-07-14 华中科技大学 Method for reducing peak-to-average power ratio of OFDM signal

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080137767A1 (en) * 2006-06-20 2008-06-12 Alcatel Lucent Ofdm clipping using sidebands and up-sampling

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101778071A (en) * 2010-01-29 2010-07-14 华中科技大学 Method for reducing peak-to-average power ratio of OFDM signal

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
蒙文武,朱光喜,胡修林等.《OFDM系统中低峰均功率比的SLM算法》.《华中科技大学学报》.2008, *

Also Published As

Publication number Publication date
CN102238126A (en) 2011-11-09

Similar Documents

Publication Publication Date Title
CN102238126B (en) Method for reducing peak-to-average power ratio of OFDM (orthogonal frequency division multiplexing)/OQAM system based on selective sequence
US9813276B2 (en) Method and apparatus for transmitting a signal
CN100556012C (en) The frequency domain equalization of single-carrier signal
CN101783781B (en) Information transmission method for lowering peak to average power ratio of OFDM system signal
CN102025681B (en) Sideband information transmission method of SLM (selective mapping method) and PTS (partial transmit sequences) in OFDM (orthogonal frequency division multiplexing) system reduce peak average ratio technology
CN105530217A (en) Signal emission and reception method of GFDM system based on weight fraction Fourier transform
AU2003251169B2 (en) Transmitting/receiving apparatus and method for reducing PAPR in an OFDM mobile communication system
CN101005474B (en) New compression extension method for reducing OFDM signal peak to average ratio
CN101340417A (en) Improved iterative PTS method for lowering peak-average-ratio in OFDM system
CN101155168A (en) Method for reducing PAPR of OFDM system
CN101848184B (en) Selective mapping scrambling method for reducing peak-to-average power ratio of orthogonal frequency division multiplexing system
CN102497350B (en) OFDM (Orthogonal Frequency Division Multiplexing) peak-to-average power ratio lowering method based on constellation linear expansion
CN102075483A (en) Method for reducing peak to average power ratio of OFDM signal
CN101753500A (en) Method and system for efficiently restraining large peak-to-average power ratio (PAPR) of OFDM system based on companding
CN105553909B (en) A kind of method that iteration limit filtration reduces ofdm signal peak-to-average force ratio
CN101534282A (en) Low-complexity OFDM signal limitation companding method
Wang et al. A low-complexity peak-to-average power ratio reduction technique for OFDM systems
CN105429911A (en) Multicarrier communication method and device having no time and frequency protection intervals
CN103457896A (en) OFDM peak-to-average ratio restraining method
Lin et al. Modified selective mapping technique for PAPR reduction in OFDM systems
CN106209717A (en) A kind of adaptive polarization QAM combined modulation method of hoisting work exoergic effect
CN111711593B (en) OFDM/OQAM-based modulation method for visible light communication system
CN101729134B (en) Generalized multiple carrier frequency division multiple access transmission equipment based on orthogonal transformation processing
CN101951354A (en) Power line carrier baseband signal modulator-demodulator based on COFDM
CN107592280A (en) A kind of method that combination channel estimation avoids PTS transmitted sideband information

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant