CN103516642B - The method and apparatus of Combined estimator interference signal physical parameter - Google Patents
The method and apparatus of Combined estimator interference signal physical parameter Download PDFInfo
- Publication number
- CN103516642B CN103516642B CN201310405521.0A CN201310405521A CN103516642B CN 103516642 B CN103516642 B CN 103516642B CN 201310405521 A CN201310405521 A CN 201310405521A CN 103516642 B CN103516642 B CN 103516642B
- Authority
- CN
- China
- Prior art keywords
- signal
- compensation
- frequency
- time delay
- estimating
- 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
Links
Landscapes
- Noise Elimination (AREA)
- Radio Relay Systems (AREA)
Abstract
The method of a kind of Combined estimator interference signal physical parameter, the method comprises the steps: to receive a mixed signal；The sampling of described mixed signal is obtained discrete observation sample of signal, and wherein, described discrete observation sample of signal includes an interference signal and a useful signal；Carry the identical local signal of information to one with described interference signal and carry out delay compensation；Estimate the frequency deviation of described interference signal；According to the frequency deviation estimated, the described local signal after delay compensation is carried out frequency compensation；According to carrying out the local signal after delay compensation and frequency compensation and the phase shift of described discrete observation sample estimated disturbance signal；According to the phase shift estimated, the local signal after carrying out delay compensation and frequency compensation is carried out phase compensation；And the amplitude of described interference signal is estimated according to the local signal after carrying out delay compensation, frequency compensation and phase shift compensation and discrete observation sample.The invention still further relates to the device of a kind of Combined estimator interference signal physical parameter.
Description
Technical Field
The invention relates to an interference elimination technology after two paths of signals are mixed and overlapped in a communication system, in particular to a method and a device for jointly estimating physical parameters of interference signals.
Background
In some application scenarios of wireless communication, two stations cannot directly communicate with each other due to the limitation of distance or geographic environment, and a relay station must be used for forwarding and amplifying signals. In order to improve the spectrum utilization and the level of information interception resistance, it is also usually necessary that the signals transmitted by two stations communicating with each other completely overlap in time domain, frequency domain, or even code space.
During communication, the two stations respectively send uplink signals with the same frequency to the relay station at the same time, the uplink signals are superposed on the relay station, and after down-conversion and power amplification of the relay station, the two stations can both receive downlink mixed signals formed by local signal components (namely interference signals) and opposite signal components (namely useful signals). Since each station knows exactly the information transmitted by itself and the signal processing procedure of local transceiving, it can make a more efficient estimation of the local signal component (interference signal) in the mixed signal. After the part of signals are eliminated from the mixed downlink signal as interference signals, reliable demodulation of the opposite signal component (useful signal) can be realized. Thus, even if the signals of the two parties are completely overlapped on the time domain, the frequency domain and the code space, the utilization rate of the channel resources can be doubled when the links are symmetrical.
In order to eliminate the interference signal as much as possible, the receiver used by the station needs to add a series of physical parameter estimation modules at the front end of the signal demodulator. The physical parameters to be estimated include frequency offset, time delay, phase shift, amplitude, etc. of the interfering signal. Whether the estimation of the physical parameters is accurate or not directly affects whether the useful signals can be demodulated correctly, and therefore, the method becomes the most key technology in the communication application scene. In order to separate the mixed signals, a generalized correlation method is generally used to estimate physical parameters of the interference signals, then the interference signals are reconstructed on the basis, and finally the interference signals are subtracted from the mixed signals to obtain useful signals required by the receiving end.
However, the conventional generalized correlation method assumes that the rest of the parameters are known parameters or does not consider one or more of the parameters before estimation, and only one parameter value can be estimated at a time. Such an assumption is not only not realistic, but also does not take into account the interplay between the physical parameters to be estimated. The physical parameters estimated by the method are not only inaccurate, but also inconvenient for real-time processing.
Disclosure of Invention
In view of the above problems, it is necessary to provide a method for jointly estimating physical parameters of an interference signal, which can accurately estimate various physical parameters of the interference signal.
In addition, it is necessary to provide a device for jointly estimating physical parameters of an interference signal, which can accurately estimate various physical parameters of the interference signal.
A method of jointly estimating a physical parameter of an interfering signal, the method comprising the steps of:
receiving a mixed signal;
sampling the mixed signal to obtain a discrete observation signal sample, wherein the discrete observation signal sample comprises an interference signal and a useful signal;
performing time delay compensation on a local signal which carries the same information as the interference signal;
estimating a frequency offset of the interfering signal;
performing frequency compensation on the local signal after the time delay compensation according to the estimated frequency offset;
estimating the phase shift of an interference signal according to the local signal after time delay compensation and frequency compensation and the discrete observation sample;
performing phase compensation on the local signal subjected to time delay compensation and frequency compensation according to the estimated phase shift; and
and estimating the amplitude of the interference signal according to the local signal and the discrete observation sample after the time delay compensation, the frequency compensation and the phase shift compensation are carried out.
An apparatus for jointly estimating a physical parameter of an interfering signal, the apparatus comprising:
the signal receiving module is used for receiving a mixed signal;
the sampling module is used for sampling the mixed signal to obtain a discrete observation signal sample, wherein the discrete observation signal sample comprises an interference signal and a useful signal;
the time delay compensation module is used for carrying out time delay compensation on a local signal which carries the same information as the interference signal;
the frequency offset estimation module is used for estimating the frequency offset of the interference signal;
the frequency compensation module is used for carrying out frequency compensation on the local signal after the time delay compensation according to the estimated frequency offset;
the phase shift estimation module is used for estimating the phase shift of the interference signal according to the local signal after time delay compensation and frequency compensation and the discrete observation sample;
the phase compensation module is used for carrying out phase compensation on the local signal subjected to time delay compensation and frequency compensation according to the estimated phase shift; and
and the amplitude estimation module is used for estimating the amplitude of the interference signal according to the local signal and the discrete observation sample after the time delay compensation, the frequency compensation and the phase shift compensation are carried out.
The method and the device for jointly estimating the physical parameters of the interference signal estimate the physical parameters of the interference signal such as time delay, frequency offset, phase shift and the like in sequence. The previous estimated parameter is compensated (aligned) in advance before estimating the next parameter to be estimated, and the overall performance of the whole device is optimized. Compared with the traditional frequency offset estimation method based on frequency subdivision and searching in a certain range, the method and the device for jointly estimating the physical parameters of the interference signals can estimate the frequency offset of the interference signals more quickly and accurately and are more suitable for real-time processing.
Drawings
Fig. 1 is a diagram of a communication system employing the method and apparatus for jointly estimating physical parameters of interfering signals according to a preferred embodiment of the present invention.
Fig. 2 is a functional block diagram of an apparatus for jointly estimating physical parameters of interfering signals according to a preferred embodiment of the present invention.
Fig. 3 is a flowchart of a method for jointly estimating physical parameters of interfering signals according to a preferred embodiment of the present invention.
Description of the main elements
| Wireless communication system | 100 |
| Site | A、B |
| Relay station | C |
| Device for jointly estimating physical parameters of interference signals | 10 |
| Signal receiving module | 11 |
| Sampling module | 12 |
| Time delay compensation module | 13 |
| Correlation operation module | 14 |
| Autocorrelation parameter calculation module | 15 |
| Frequency offset estimation module | 16 |
| Frequency compensation module | 17 |
| Phase shift estimation module | 18 |
| Phase compensation module | 19 |
| Amplitude estimation module | 20 |
The following detailed description will further illustrate the invention in conjunction with the above-described figures.
Detailed Description
The method and apparatus for jointly estimating physical parameters of an interference signal according to the preferred embodiment of the present invention are applied to a wireless communication system, and are used for jointly estimating various physical parameters of the interference signal.
First, various symbols, functions, formulas, and the like used in the present invention are defined and explained.
For any one of plural numbers，RepresentsThe real part of (a) is,representsThe imaginary part of (a) is,is an imaginary unit of the complex number.Representing and gettingIs conjugated, i.e.。Representative pairPerforming modulo arithmetic, i.e.。
Any one of pluralAre all aligned with one in a complex plane with originThe starting point is the point of the process,vectors that are endpoints correspond one-to-one. The argument of the complex number is the vector starting from the positive half axis of the x-axisIn the ray (starting point isAt an end point of) Corner being terminal edge. Therefore, a plurality ofCan also be expressed asWherein。
Any one of pluralHas an infinite number of values and the values differ from each otherInteger multiples of. Is suitable forSpoke angle ofValue of (A) isThe main value of (A) is recorded as。
Referring to fig. 1, the method and apparatus for jointly estimating physical parameters of interfering signals according to the preferred embodiment of the present invention are applied to a wireless communication system 100, which includes a station a, a station B and a relay station C. The relay station C is used for forwarding and amplifying signals sent by the station a and the station B. Specifically, when communication is required between the station a and the station B, the station a and the station B transmit signals of the same frequency to the relay station C at the same time, respectively. For example, let station A output a signal of(ii) a The signal output by station B is. The relay station C superposes the signals sent by the station A and the station B to form a mixed signalAnd for the mixed signalAnd performing frequency conversion and power amplification, and outputting the frequency-converted and power-amplified signals to a site A and a site B respectively. That is, the relay station C only processes the mixed signalFrequency conversion and power amplification, i.e. performing transparent forwarding functions, without mixing signalsAnd performing other signal processing. Thus, not only will site a receive the signal from site B that is forwarded via relay CAt this timeThe station a receives a useful signal from the station a itself and a signal forwarded by the relay station CAt this timeIs an interfering signal with respect to site a. Similarly, station B will not only receive the signal from station a that is forwarded via relay CAt this timeThe useful signal relative to the station B is received by the signal transmitted by the relay station C from the station BAt this timeIs an interfering signal with respect to site B.
No matter which one of the station a or the station B is used as the observation station, if the information of the useful signal needs to be demodulated correctly locally, the physical parameters of the interference signal transmitted locally and forwarded by the relay station C, including time delay, frequency offset, phase shift, amplitude, etc., need to be estimated first.
Referring to fig. 2, an apparatus 10 for jointly estimating physical parameters of an interference signal according to a preferred embodiment of the present invention includes a signal receiving module 11, a sampling module 12, a delay compensation module 13, a correlation operation module 14, an autocorrelation parameter calculation module 15, a frequency offset estimation module 16, a frequency compensation module 17, a phase shift estimation module 18, a phase compensation module 19, and an amplitude estimation module 20. The function of the modules of the apparatus 10 for jointly estimating the physical parameter of the interfering signal will be described in detail in fig. 3.
Referring to fig. 3, a method and an apparatus for jointly estimating physical parameters of an interfering signal according to the present invention are described in the following with reference to a station a as an example.
Step S1: the signal receiving module 11 receives the mixed signal sent by the relay station C。The complex baseband model of (a) is expressed as:
wherein,andwaveforms of an interference signal and a useful signal relative to a site A are respectively, the mean value is 0, the variance is 1, and statistics are independent;andthe amplitudes of the interfering signal and the useful signal, respectively;is complex white gaussian noise.In order to delay the transmission of the interfering signal,is a carrier frequency offset and is,is the carrier phase shift.
Step S2: the sampling module 12 is used for receiving the mixed signalSampling to obtain discrete observation samples. In this embodiment, the sampling is done at the symbol rate, i.e. the sampling rate，Is the duration of a unit symbol. The sampling module 12 first applies the received mixed signalMatched filtering and then mixing the signalsAnd (5) extracting. Assuming transmission delay of interfering signalsThe exact value of (A) can be obtained, corresponding to the point of unsigned crosstalk, the instant of decimationThen obtaining discrete observation samplesIt is expressed as:
wherein,、andrespectively representing the interfering signalsUseful signalAnd noiseThe resulting samples are sampled at the symbol rate.Transmission delay for interference signalsAnd after sampling, corresponding discrete values.Is from 1 to，To observe the length.
Step S3: the delay compensation module 13 is used for sampling the one-sum interference signalAnd (3) carrying local signals with the same information to perform time delay compensation, namely:
the local signal is stored at site a,namely the local signal after time delay compensation according to the symbol synchronization method.
Step S4: the correlation operation module 14 compares each obtained discrete observation sample with the local signal after the completion of the delay compensationPerforming a correlation operation to obtain a correlation signal：
Wherein, for two complex symbolsAndthe correlation operation meansMultiplication byIs conjugated, i.e.. Therefore, the temperature of the molten metal is controlled,to representAndby correlation, i.e.。
Step S5: the autocorrelation parameter calculation module 15 calculates the correlation signalIs/are as followsIndividual autocorrelation parameter：
Wherein,is from 1 to，Is the number of symbol intervals, satisfies。Representing an integer number up (plus infinity).
Step S6: the frequency offset estimation module 16 is based onIndividual autocorrelation parameterEstimating interfering signalsFrequency deviation of：
The larger the value of (A), the frequency offsetThe more accurate the estimated value. In the present embodiment, it is considered thatThe limitation of the complexity of the calculation,value of (2) is taken to maximum。
Step S7: the frequency compensation module 17 estimates the frequency offsetTime-delay compensated local signalAnd performing frequency compensation to align the frequency and update the local signal. The updated local signal is：
That is, the signalThe local signal is the local signal after time delay compensation and frequency compensation are made.
Step S8: the phase shift estimation module 18 compensates the local signal according to the time delay and frequencyAnd discrete observation samplesEstimating interfering signalsPhase shift of：
Step S9: the phase compensation module 19 is based on the estimated phase shiftFor the local signal after time delay compensation and frequency compensationPhase compensation is performed to align the phases and the local signal is updated again. The local signal after being updated again is：
That is, the signalThe local signal is subjected to time delay compensation, frequency compensation and phase shift compensation.
Step S10: the amplitude estimation module 20 performs time delay compensation, frequency compensation and phase shift compensation according to the local signalAnd discrete observation samplesEstimating interfering signalsAmplitude of (2)：
The method and the device for jointly estimating the physical parameters of the interference signal estimate the physical parameters of the interference signal such as time delay, frequency offset, phase shift and the like in sequence. The previous estimated parameter is compensated (aligned) in advance before estimating the next parameter to be estimated, and the overall performance of the whole device is optimized. Compared with the traditional frequency offset estimation method based on frequency subdivision and searching in a certain range, the method and the device for jointly estimating the physical parameters of the interference signals can estimate the frequency offset of the interference signals more quickly and accurately and are more suitable for real-time processing.
Claims (10)
1. A method of jointly estimating a physical parameter of an interfering signal, the method comprising the steps of:
(a) receiving a mixed signal;
(b) sampling the mixed signal to obtain a discrete observation sample, wherein the discrete observation sample comprises an interference signal and a useful signal;
(c) performing time delay compensation on a local signal which carries the same information as the interference signal;
(d) estimating a frequency offset of the interfering signal;
(e) performing frequency compensation on the local signal after the time delay compensation according to the estimated frequency offset;
(f) estimating the phase shift of an interference signal according to the local signal after time delay compensation and frequency compensation and the discrete observation sample;
(g) performing phase compensation on the local signal subjected to time delay compensation and frequency compensation according to the estimated phase shift; and
(h) and estimating the amplitude of the interference signal according to the local signal and the discrete observation sample after time delay compensation, frequency compensation and phase compensation are carried out.
2. The method of claim 1, wherein the step (d) comprises the following steps:
performing correlation operation on each obtained discrete observation sample and the local signal after time delay compensation is completed to obtain a correlation signal zk：Wherein r iskRepresenting the discrete observation samples in question,to complete the delay compensated local signal,to representThe value of k is from 1 to N, and N is the observation length;
calculating L autocorrelation parameters p of the correlation signall：Wherein, the value of L is from 1 to L, L is a symbol interval number, L is less than or equal to ceil (N/2), ceil represents an integer taken upwards (in a positive infinite direction); and
based on L autocorrelation parameters plEstimating a frequency offset of the interfering signal
Wherein, TsIs the duration of a unit symbol.
3. The method for jointly estimating physical parameters of an interfering signal according to claim 2, wherein in step (e), the frequency compensation is performed on the delay-compensated local signal by the following formula:
wherein,the local signal after time delay compensation and frequency compensation is shown.
4. A method for jointly estimating a physical parameter of an interfering signal according to claim 3, characterized in that: in step (f), the estimated phase shift of the interference signalComprises the following steps:and in step (g), the local signals after time delay compensation, frequency compensation and phase compensation are carried outComprises the following steps:
5. the method of jointly estimating physical parameters of an interfering signal according to claim 4, characterized in that: in step (h), the estimated amplitude of the interference signalComprises the following steps:
6. an apparatus for jointly estimating a physical parameter of an interfering signal, the apparatus comprising:
the signal receiving module is used for receiving a mixed signal;
the sampling module is used for sampling the mixed signal to obtain a discrete observation sample, wherein the discrete observation sample comprises an interference signal and a useful signal;
the time delay compensation module is used for carrying out time delay compensation on a local signal which carries the same information as the interference signal;
the frequency offset estimation module is used for estimating the frequency offset of the interference signal;
the frequency compensation module is used for carrying out frequency compensation on the local signal after the time delay compensation according to the estimated frequency offset;
the phase shift estimation module is used for estimating the phase shift of the interference signal according to the local signal after time delay compensation and frequency compensation and the discrete observation sample;
the phase compensation module is used for carrying out phase compensation on the local signal subjected to time delay compensation and frequency compensation according to the estimated phase shift; and
and the amplitude estimation module is used for estimating the amplitude of the interference signal according to the local signal and the discrete observation sample after time delay compensation, frequency compensation and phase compensation are carried out.
7. The apparatus for jointly estimating physical parameters of an interfering signal as claimed in claim 6, characterized in that the apparatus further comprises:
a correlation operation module for performing correlation operation on each discrete observation sample and the local signal after completion of the time delay compensation to obtain a correlation signal zk：Wherein r iskRepresenting the discrete observation samples in question,to complete the delay compensated local signal,to representThe value of k is from 1 to N, and N is the observation length;
an autocorrelation parameter calculation module for calculating L autocorrelation parameters p of the correlation signall：
Wherein, the value of L is from 1 to L, L is a symbol interval number, L is less than or equal to ceil (N/2), ceil represents an upward (positive infinite direction) integer, and N is an observation length;
the frequency deviation estimation module is used for estimating the frequency deviation according to L autocorrelation parameters plEstimating a frequency offset of the interfering signal
Wherein, TsIs the duration of a unit symbol.
8. The apparatus for jointly estimating physical parameters of an interfering signal according to claim 7, characterized in that: the frequency compensation module performs frequency compensation on the local signal after the time delay compensation according to the following formula:
wherein,the local signal after time delay compensation and frequency compensation is shown.
9. The apparatus for jointly estimating physical parameters of an interfering signal according to claim 8, characterized in that: the phase shift of the interference signal estimated by the phase shift estimation moduleComprises the following steps:and the local signal after time delay compensation, frequency compensation and phase compensation is madeComprises the following steps:
10. the apparatus for jointly estimating physical parameters of an interfering signal according to claim 9, characterized in that: the amplitude of the interference signal estimated by the amplitude estimation moduleComprises the following steps:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310405521.0A CN103516642B (en) | 2013-09-09 | 2013-09-09 | The method and apparatus of Combined estimator interference signal physical parameter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310405521.0A CN103516642B (en) | 2013-09-09 | 2013-09-09 | The method and apparatus of Combined estimator interference signal physical parameter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103516642A CN103516642A (en) | 2014-01-15 |
| CN103516642B true CN103516642B (en) | 2017-01-04 |
Family
ID=49898696
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310405521.0A Expired - Fee Related CN103516642B (en) | 2013-09-09 | 2013-09-09 | The method and apparatus of Combined estimator interference signal physical parameter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103516642B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111585926B (en) * | 2020-04-28 | 2021-07-27 | 北京升哲科技有限公司 | Multi-antenna diversity receiver |
| CN112039613B (en) * | 2020-08-07 | 2022-07-29 | 中国卫通集团股份有限公司 | Processing method and device for asymmetric PCMA (pulse code division multiple Access) mixed signal |
| CN113258984B (en) * | 2021-04-29 | 2022-08-09 | 东方红卫星移动通信有限公司 | Multi-user self-adaptive frequency offset elimination method and device, storage medium and low-orbit satellite communication system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6791964B1 (en) * | 1998-12-31 | 2004-09-14 | Samsung Electronics Co., Ltd. | Multicarrier CDMA rake receiver apparatus |
| CN101741776A (en) * | 2009-11-10 | 2010-06-16 | 华为技术有限公司 | Method and device for eliminating interference signals |
| CN102932041A (en) * | 2012-11-21 | 2013-02-13 | 西安电子科技大学 | Method for encoding and decoding asynchronous space-time code for collaborative multi-point transmission |
-
2013
- 2013-09-09 CN CN201310405521.0A patent/CN103516642B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6791964B1 (en) * | 1998-12-31 | 2004-09-14 | Samsung Electronics Co., Ltd. | Multicarrier CDMA rake receiver apparatus |
| CN101741776A (en) * | 2009-11-10 | 2010-06-16 | 华为技术有限公司 | Method and device for eliminating interference signals |
| CN102932041A (en) * | 2012-11-21 | 2013-02-13 | 西安电子科技大学 | Method for encoding and decoding asynchronous space-time code for collaborative multi-point transmission |
Non-Patent Citations (1)
| Title |
|---|
| 多天线信号联合接收的合成技术研究;沈彩耀;《中国博士学位论文电子期刊网》;20120715;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103516642A (en) | 2014-01-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106461749B (en) | For determining the method, apparatus and computer program of arrival time | |
| CN102111205B (en) | Channel estimation for communication system with multiple transmitting antennas | |
| US7561643B1 (en) | DC offset correction in mobile communication system | |
| EP2194665B1 (en) | CFO and I/Q imbalance estimation and compensation | |
| CN107911329A (en) | A kind of signal analyzer ofdm signal demodulation method | |
| CN105490973A (en) | IQ signal calibration method and apparatus | |
| Marey et al. | EM-based joint channel estimation and IQ imbalances for OFDM systems | |
| CN103516642B (en) | The method and apparatus of Combined estimator interference signal physical parameter | |
| CN102263710B (en) | Channel estimation method and device as well as frequency offset estimation method and device | |
| US20110206105A1 (en) | Method for determining hybrid domain compensation parameters for analog loss in ofdm communication systems and compensating for the same | |
| CN106375251A (en) | Noise power estimating method and apparatus, and user equipment | |
| TW201406097A (en) | Method and device for estimating channel and noise | |
| CN103457630B (en) | Interference signal amplitude method of estimation under a kind of complex field and device | |
| CN103188067B (en) | A kind of chip clock frequency departure estimation error of spread spectrum system and the method for correction | |
| CN106130573B (en) | The pilot signal method for generation of NB-IoT pilot signal generating device | |
| CN105162533B (en) | Transmitter amplitude imbalance and phase imbalance measuring method | |
| Yao et al. | Low-complexity timing synchronization for decode-and-forward cooperative communication systems with multiple relays | |
| EP3035622B1 (en) | Maximum likelihood sequence detection in the phase domain | |
| Rice et al. | On frequency offset estimation using the iNET preamble in frequency selective fading | |
| CN105075131B (en) | Device and method for frequency offset estimation | |
| CN104052690A (en) | Frequency synchronization method for cooperative communication system | |
| CN100373792C (en) | Decoding method of received data in short range radio network | |
| Guan et al. | Iterative phase noise suppression for full-duplex communication systems | |
| US20020009162A1 (en) | Method and apparatus for correcting a signal | |
| KR100764356B1 (en) | A timing estimator in a oqpsk demodulator |
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 | ||
| 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: 20170104 Termination date: 20210909 |