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 PDF

Info

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
interference signal
frequency
local
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
CN201310405521.0A
Other languages
Chinese (zh)
Other versions
CN103516642A (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.)
Shenzhen Research Institute Tsinghua University
Original Assignee
Shenzhen Research Institute Tsinghua University
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 Shenzhen Research Institute Tsinghua University filed Critical Shenzhen Research Institute Tsinghua University
Priority to CN201310405521.0A priority Critical patent/CN103516642B/en
Publication of CN103516642A publication Critical patent/CN103516642A/en
Application granted granted Critical
Publication of CN103516642B publication Critical patent/CN103516642B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

The method and apparatus of Combined estimator interference signal physical parameter
Technical field
The present invention relates to the interference cancellation techniques after two paths of signals aliasing in communication system, particularly relate to a kind of Combined estimator The method and apparatus of interference signal physical parameter.
Background technology
In some application scenarios of radio communication, cannot be straight due to the restriction of distance or geographical environment between two websites Connect and communicate, it is necessary to utilize relay station to carry out forwarding and the amplification of signal.In order to improve the anti-of the availability of frequency spectrum and information Intercepting and capturing level, the signal that two websites generally also needing to be in communication with each other are launched is complete on time domain, frequency domain even code space Overlapping.
During communication, both sides' website sends the upward signal of same frequency in the identical time to relay station respectively, on relaying Superposition, after the down coversion and power amplification of relaying, two websites all can receive and (i.e. be disturbed letter by local signal composition Number) downmix chan-nel that forms with the other side's signal component (i.e. useful signal).Owing to each website knows for sure oneself institute The information sent and the signal processing of local transceiver, therefore, it is possible to the local signal composition (interference in mixed signal Signal) make and the most effectively estimating.After this part signal is eliminated from mixing downstream signal as interference signal, Realize the reliable demodulation to the other side's signal component (useful signal).Such that make both sides' signal on time domain, frequency domain, code space Completely overlapped, channel resource utilization rate during link symmetry can also double.
In order to eliminated as much as possible by interference signal, the receiver that website uses is accomplished by before demodulator of PM signal PM End increases a series of physical parameter estimation module.Physical parameter to be estimated include disturb the frequency deviation of signal, time delay, phase shift, Amplitude etc..Can the most whether the estimation of these physical parameters will directly influence useful signal be correctly demodulated, and therefore becomes For the most critical technology under this class communications applications scene.In order to realize the separation of mixed signal, the most first use generalized correlation Method estimates the physical parameter of interference signal, reconstructs interference signal the most on this basis, and signal will be disturbed the most again from mixed Close in signal and deduct the useful signal obtained required for receiving terminal.
But, traditional generalized correlation method is before estimation or assumes that remaining parameter is all known parameters, or does not considers Some or certain several parameters therein, and a parameter value can only be estimated every time.Such hypothesis does not only meet reality Border, does not accounts for influencing each other between physical parameter to be estimated yet.Adopt the physical parameter estimated in this way not But inaccurate, also inconvenience processes in real time.
Summary of the invention
For the problems referred to above, it is necessary to the method providing a kind of Combined estimator interference signal physical parameter, it can accurately be estimated Count out every physical parameter of interference signal.
It addition, there is a need to provide the device of a kind of Combined estimator interference signal physical parameter, it can accurately estimate dry Disturb every physical parameter of signal.
The method of a kind of Combined estimator interference signal physical parameter, the method comprises the steps:
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 One 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;
Interference is estimated according to the local signal after carrying out delay compensation and frequency compensation and described discrete observation sample The phase shift of signal;
According to the phase shift estimated, the local signal after carrying out delay compensation and frequency compensation is carried out phase compensation;With And
Estimate according to the local signal after carrying out delay compensation, frequency compensation and phase shift compensation and discrete observation sample The amplitude of described interference signal.
A kind of device of Combined estimator interference signal physical parameter, this device includes:
Signal receiving module, for receiving a mixed signal;
Sampling module, for obtaining discrete observation sample of signal, wherein, described discrete observation to the sampling of described mixed signal Sample of signal includes an interference signal and a useful signal;
Delay compensation module, carries out time delay benefit for carrying the identical local signal of information to one with described interference signal Repay;
Frequency deviation estimating modules, for estimating the frequency deviation of described interference signal;
Frequency compensation module, for carrying out frequency benefit according to the frequency deviation estimated to the described local signal after delay compensation Repay;
Phase shift estimation module, for according to carrying out the local signal after delay compensation and frequency compensation and described discrete The phase shift of observation sample estimated disturbance signal;
Phase compensation block, for believing this locality after carrying out delay compensation and frequency compensation according to the phase shift estimated Number carry out phase compensation;And
Amplitude Estimation module, for according to carry out the local signal after delay compensation, frequency compensation and phase shift compensation with And the amplitude of the discrete observation sample described interference signal of estimation.
The method and device of described Combined estimator interference signal physical parameter by estimation of the order go out to disturb signal time Prolong, the physical parameter such as frequency deviation and phase shift.Before estimating next one parameter to be estimated, upper one of pre-compensation (alignment) has estimated ginseng Number, optimizes the overall performance of whole device.Segment based on frequency relative to traditional, and scan in certain scope Frequency deviation estimating method, use Combined estimator provided by the present invention interference signal physical parameter method and device can be more Estimate the frequency deviation of interference signal rapidly and accurately, be more applicable for processing in real time.
Accompanying drawing explanation
Fig. 1 is the logical of the method and device of the Combined estimator interference signal physical parameter using better embodiment of the present invention The schematic diagram of communication system.
Fig. 2 is the functional block diagram of the device of the Combined estimator interference signal physical parameter of better embodiment of the present invention.
Fig. 3 is the flow chart of the method for the Combined estimator interference signal physical parameter of better embodiment of the present invention.
Main element symbol description
Wireless communication system 100
Website A、B
Relay station C
The device of Combined estimator interference signal physical parameter 10
Signal receiving module 11
Sampling module 12
Delay compensation module 13
Related operation module 14
Auto-correlation Parameters Calculation module 15
Frequency deviation estimating modules 16
Frequency compensation module 17
Phase shift estimation module 18
Phase compensation block 19
Amplitude Estimation module 20
Following detailed description of the invention will further illustrate the present invention in conjunction with above-mentioned accompanying drawing.
Detailed description of the invention
It is wireless that the method and device of the Combined estimator interference signal physical parameter of better embodiment of the present invention is applied to one Communication system, for carrying out Combined estimator to every physical parameter of interference signal.
First, various symbols, function and the formula etc. that use the present invention are defined and illustrate.
For any one plural number,RepresentReal part,RepresentImaginary part,Imaginary unit for plural number.Representative takesConjugation, i.e.Represent RightCarry out modulus computing, i.e.
Any one plural numberAll with one in complex plane with initial pointFor initial point,Vector one a pair for terminal Should.Argument of a complex number is with the positive axis of x-axis for an initial line, vector(starting point is the ray at place, terminal is) it is The angle on whole limit.Therefore, plural numberIt is also denoted as, wherein
Any one plural numberArgument have between unlimited number of value, and these values differIntegral multiple.Being suitable forArgumentValue, be calledThe main value of argument, be denoted as
Referring to Fig. 1, the method and device of the Combined estimator interference signal physical parameter of better embodiment of the present invention should For a wireless communication system 100, it includes website A, website B and relay station C.Relay station C is used for forwarding and amplify website A And the signal that website B sends.Specifically, when needing to communicate between website A with website B, website A and website B is respectively identical Time to relay station C send same frequency signal.Such as, if the signal of website A output is;The letter of website B output Number be.The signal that website A and website B sends is overlapped by relay station C, forms a mixed signal, and to described Mixed signalCarry out frequency conversion and power amplification, then export respectively to website A and website B.It is to say, relay station C Only to mixed signalCarry out frequency conversion and power amplification, i.e. complete transparent forwarding function, and not to mixed signal Do other signal processing.So, website A not only can receive the signal forwarded from website B through relay station C, nowIt is useful signal relative to website A, also can receive the signal forwarded through relay station C from website A self , nowIt is interference signal relative to website A.Similarly, website B not only can receive and turn through relay station C from website A The signal sent out, nowIt is useful signal relative to website B, also can receive and relay from passing through of website B self The signal that the C that stands forwards, nowIt is interference signal relative to website B.
No matter website A or website B which side as observation station, if the letter of useful signal correctly need to be demodulated in this locality Breath, then need first to estimate the physical parameter of the local interference signal forwarded through relay station C self launched, including Time delay, frequency deviation, phase shift, amplitude etc..
Incorporated by reference to refering to Fig. 2, the device 10 of the Combined estimator interference signal physical parameter of better embodiment of the present invention wraps Include signal receiving module 11, sampling module 12, delay compensation module 13, related operation module 14, auto-correlation Parameters Calculation module 15, frequency deviation estimating modules 16, frequency compensation module 17, phase shift estimation module 18, phase compensation block 19 and amplitude Estimation mould Block 20.The function of each module of the device 10 of Combined estimator interference signal physical parameter will be described in detail in figure 3.
Refer to Fig. 3, as a example by website A, the Combined estimator of the present invention is disturbed method and the dress of signal physical parameter below Put and illustrate.
Step S1: signal receiving module 11 receives the mixed signal that relay station C sendsComplex baseband model table It is shown as:
Wherein,WithBeing respectively the interference signal relative to website A and the waveform of useful signal, average is 0, side Difference is 1, and statistical iteration;WithIt is respectively interference signal and the amplitude of useful signal;For white complex gaussian noise.For The propagation delay time of interference signal,For carrier wave frequency deviation,For carrier phase.
Step S2: the sampling module 12 mixed signal to receivingSample, to obtain discrete observation sample。 In the present embodiment, sample according to character rate, i.e. sampling rate,During for unit code element lasting Between.Sampling module 12 is first to the mixed signal receivedCarry out matched filtering, then to mixed signalTake out Take.Assume to disturb the propagation delay time of signalExact value can obtain, corresponding to without intersymbol interference point, the extraction moment, then discrete observation sample can be obtained, it is expressed as:
Wherein,AndRepresent respectively interference signal, useful signalAnd noiseAccording to Character rate carries out gained sample of sampling.For disturbing the propagation delay time of signalCentrifugal pump corresponding after sampling.'s Value from 1 to,For observed length.
Step S3: delay compensation module 13 is to one and disturbs sample of signalThe local signal that information of carrying is identical is carried out Delay compensation, it may be assumed that
Described local signal is stored in website A,I.e. for carrying out the local letter after delay compensation according to symbol timing synchronization method Number.
Step S4: related operation module 14 each discrete observation sample to obtaining with complete this locality after delay compensation SignalCarry out related calculation, obtain a coherent signal:
Wherein, for two complex symbolsWith, related operation refers toIt is multiplied byConjugation, i.e..Therefore,RepresentWithRelated operation, i.e.
Step S5: auto-correlation Parameters Calculation module 15 calculates coherent signal'sIndividual auto-correlation parameter:
Wherein,Value from 1 to,For mark space number, meetRepresent upwards (just without Poor direction) round numbers.
Step S6: frequency deviation estimating modules 16 basisIndividual auto-correlation parameterEstimated disturbance signalFrequency deviation:
Value the biggest, then frequency deviationThe value estimated is the most accurate.In the present embodiment, it is contemplated that computation complexity Limit,Value maximum take
Step S7: frequency compensation module 17 is according to the frequency deviation estimatedTo the local signal after delay compensationCarry out Frequency compensation, with frequency of aliging, updates local signal.Local signal after renewal is:
It is to say, signalIt it is the local signal after carrying out delay compensation and frequency compensation.
Step S8: phase shift estimation module 18 is according to carrying out the local signal after delay compensation and frequency compensationAnd Discrete observation sampleEstimated disturbance signalPhase shift:
Step S9: phase compensation block 19 is according to the phase shift estimatedAfter carrying out delay compensation and frequency compensation Local signalCarry out phase compensation, with alignment phase, again update local signal.Local signal after again updating is:
It is to say, signalIt it is the local signal after carrying out delay compensation, frequency compensation and phase shift compensation.
Step S10: amplitude Estimation module 20 is believed according to this locality after carrying out delay compensation, frequency compensation and phase shift compensation NumberAnd discrete observation sampleEstimated disturbance signalAmplitude:
The method and device of described Combined estimator interference signal physical parameter by estimation of the order go out to disturb signal time Prolong, the physical parameter such as frequency deviation and phase shift.Before estimating next one parameter to be estimated, upper one of pre-compensation (alignment) has estimated ginseng Number, optimizes the overall performance of whole device.Segment based on frequency relative to traditional, and scan in certain scope Frequency deviation estimating method, use Combined estimator provided by the present invention interference signal physical parameter method and device can be more Estimate the frequency deviation of interference signal rapidly and accurately, be more applicable for processing in real time.

Claims (10)

1. a method for Combined estimator interference signal physical parameter, the method comprises the steps:
A () receives a mixed signal;
B () obtains discrete observation sample to the sampling of described mixed signal, wherein, described discrete observation sample includes an interference signal And a useful signal;
C () is carried the identical local signal of information to one with described interference signal and is carried out delay compensation;
D () estimates the frequency deviation of described interference signal;
E () carries out frequency compensation according to the frequency deviation estimated to the described local signal after delay compensation;
F () estimates interference letter according to the local signal after carrying out delay compensation and frequency compensation and described discrete observation sample Number phase shift;
G () carries out phase compensation according to the phase shift estimated to the local signal after carrying out delay compensation and frequency compensation;With And
H () is estimated according to the local signal after carrying out delay compensation, frequency compensation and phase compensation and discrete observation sample The amplitude of described interference signal.
2. the method for Combined estimator interference signal physical parameter as claimed in claim 1, it is characterised in that step (d) is concrete Comprise the steps:
To each the discrete observation sample obtained with complete delay compensation after local signal carry out related calculation, obtain a phase OFF signal zk:Wherein, rkRepresent described discrete observation sample,Believe for this locality after completing delay compensation Number,RepresentConjugation, the value of k is from 1 to N, and N is observed length;
Calculate L auto-correlation parameter p of described coherent signall:Wherein, the value of l from 1 to L, L For mark space number, meet L≤ceil (N/2), ceil and represent upwards (positive infinity) round numbers;And
According to L auto-correlation parameter plEstimate the frequency deviation of described interference signal
Wherein, TsPersistent period for unit code element.
3. the method for Combined estimator interference signal physical parameter as claimed in claim 2, it is characterised in that in step (e), By equation below the local signal after delay compensation carried out frequency compensation:
Wherein,Represent the local letter after carrying out delay compensation and frequency compensation Number.
4. the method for Combined estimator interference signal physical parameter as claimed in claim 3, it is characterised in that: in step (f), The phase shift of the described interference signal estimatedFor:And in step (g), carry out time delay Local signal after compensation, frequency compensation and phase compensationFor:
5. the method for Combined estimator interference signal physical parameter as claimed in claim 4, it is characterised in that: in step (h), The amplitude of the described interference signal estimatedFor:
6. the device of a Combined estimator interference signal physical parameter, it is characterised in that this device includes:
Signal receiving module, for receiving a mixed signal;
Sampling module, for the sampling of described mixed signal is obtained discrete observation sample, wherein, described discrete observation sample includes One interference signal and a useful signal;
Delay compensation module, carries out delay compensation for carrying the identical local signal of information to one with described interference signal;
Frequency deviation estimating modules, for estimating the frequency deviation of described interference signal;
Frequency compensation module, for carrying out frequency compensation according to the frequency deviation estimated to the described local signal after delay compensation;
Phase shift estimation module, the local signal after basis carries out delay compensation and frequency compensation and described discrete observation The phase shift of sample estimated disturbance signal;
Phase compensation block, for entering the local signal after carrying out delay compensation and frequency compensation according to the phase shift estimated Line phase compensates;And
Amplitude Estimation module, for according to carry out the local signal after delay compensation, frequency compensation and phase compensation and from Dissipate observation sample and estimate the amplitude of described interference signal.
7. the device of Combined estimator interference signal physical parameter as claimed in claim 6, it is characterised in that this device also wraps Include:
Related operation module, for each the discrete observation sample obtained with complete delay compensation after local signal do phase Close computing, obtain a coherent signal zk:Wherein, rkRepresent described discrete observation sample,During for completing Prolong the local signal after compensation,RepresentConjugation, the value of k is from 1 to N, and N is observed length;
Auto-correlation Parameters Calculation module, for calculating L auto-correlation parameter p of described coherent signall:
Wherein, the value of l is from 1 to L, and L is mark space number, meets L≤ceil (N/2), ceil Representing upwards (positive infinity) round numbers, N is observed length;
Described frequency deviation estimating modules is according to L auto-correlation parameter plEstimate the frequency deviation of described interference signal
Wherein, TsPersistent period for unit code element.
8. the device of Combined estimator interference signal physical parameter as claimed in claim 7, it is characterised in that: described frequency compensation Module carries out frequency compensation according to equation below to the local signal after delay compensation:
Wherein,Represent the local letter after carrying out delay compensation and frequency compensation Number.
9. the device of Combined estimator interference signal physical parameter as claimed in claim 8, it is characterised in that: phase shift estimation module The phase shift of the described interference signal estimatedFor:And carry out delay compensation, frequency compensation And the local signal after phase compensationFor:
10. the device of Combined estimator interference signal physical parameter as claimed in claim 9, it is characterised in that: amplitude Estimation mould The amplitude of the described interference signal that block estimatesFor:
CN201310405521.0A 2013-09-09 2013-09-09 The method and apparatus of Combined estimator interference signal physical parameter Expired - Fee Related CN103516642B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
多天线信号联合接收的合成技术研究;沈彩耀;《中国博士学位论文电子期刊网》;20120715;全文 *

Also Published As

Publication number Publication date
CN103516642A (en) 2014-01-15

Similar Documents

Publication Publication Date Title
CN103929394A (en) High-precision frequency offset estimation method based on iteration algorithm
CN103516642B (en) The method and apparatus of Combined estimator interference signal physical parameter
WO2013121958A1 (en) Receiving device, receiving method and receiving program
TWI728028B (en) Chip, user equipment and method for enhanced channel estimation, manufacturing method of chip
WO2019076210A1 (en) Time domain correlation estimation method, apparatus, and device
CN104537257A (en) Distributed self-adaptation direct positioning method based on time difference
CN108541061A (en) The method and apparatus of enhancing reference signal time difference for LTE positioning
CN106936741A (en) A kind of mimo channel method of estimation based on Kalman filter
CN107466054A (en) Disturb localization method, device and the base station of cell-of-origin
TWI484787B (en) Method and device for estimating channel and noise
CN107342960A (en) A kind of unbound nucleus frequency deviation estimating method of suitable Amplitude phase shift keying
CN103188067B (en) A kind of chip clock frequency departure estimation error of spread spectrum system and the method for correction
JP6561471B2 (en) COMMUNICATION DEVICE, COMMUNICATION METHOD, AND COMMUNICATION PROGRAM
CN106878213B (en) LTE uplink frequency offset estimation method
CN103457630B (en) Interference signal amplitude method of estimation under a kind of complex field and device
Coutinho et al. 5G cascaded channel estimation using convolutional neural networks
Wang et al. Carrier phase tracking architecture for positioning in LTE networks under channel fading conditions
US10277440B1 (en) Determining common phase error
CN106375251A (en) Noise power estimating method and apparatus, and user equipment
CN105721075B (en) Data are aided in frequency reuse satellite interference signal intensity measuring method and device
CN105075131B (en) Device and method for frequency offset estimation
CN105072063A (en) Method and device for suppressing interference signal
CN106936481A (en) A kind of iterative channel estimation method based on Kalman filter
CN104779993A (en) Deep space measurement and control system and method on basis of frequency-domain equalization
CN108353066A (en) Joint non-coherent demodulation based on nonlinear filtering and offset correction of carrier frequency

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