CN104301271A - Frequency spectrum conversion method based on polarization rotational angular velocity discretization encoding - Google Patents

Frequency spectrum conversion method based on polarization rotational angular velocity discretization encoding Download PDF

Info

Publication number
CN104301271A
CN104301271A CN201410586362.3A CN201410586362A CN104301271A CN 104301271 A CN104301271 A CN 104301271A CN 201410586362 A CN201410586362 A CN 201410586362A CN 104301271 A CN104301271 A CN 104301271A
Authority
CN
China
Prior art keywords
angular speed
symbolic information
frequency
rotation angular
polarized rotation
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.)
Granted
Application number
CN201410586362.3A
Other languages
Chinese (zh)
Other versions
CN104301271B (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.)
China Academy of Space Technology CAST
Original Assignee
China Academy of Space Technology CAST
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 China Academy of Space Technology CAST filed Critical China Academy of Space Technology CAST
Priority to CN201410586362.3A priority Critical patent/CN104301271B/en
Publication of CN104301271A publication Critical patent/CN104301271A/en
Application granted granted Critical
Publication of CN104301271B publication Critical patent/CN104301271B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B14/00Transmission systems not characterised by the medium used for transmission
    • H04B14/002Transmission systems not characterised by the medium used for transmission characterised by the use of a carrier modulation
    • H04B14/008Polarisation modulation

Abstract

The invention relates to a frequency spectrum conversion method based on polarization rotational angular velocity discretization encoding. According to the method, different polarization rotational angular velocities are adopted to modulate signals, values of the polarization rotational angular velocities are periodic rectangular pulse functions, main energy of the modulated signals carrying symbol information is distributed at set frequency points, the energy of the modulation signals distributed outside the frequency points is extremely low and will not interfere in other communication systems, and therefore occupation of the signals on frequency spectrum resources is reduced. The different polarization rotational angular velocities are converted into angular velocity Doppler frequency shifts at received ends, so that signal energy is distributed at positions of frequency points deviating from a centre carrier frequency, interference signals not having polarization rotational angular velocities are still located at frequency points of the centre carrier frequency, and therefore the interference signals located at the frequency points of the centre carrier frequency and effective signals deviating from the frequency points of the centre carrier frequency can be effectively distinguished through a band-pass filter, and the interference resistance ability of a communication system can be improved.

Description

A kind of Spectrum Conversion method based on polarized rotation angular speed discretization coding
Technical field
The present invention relates to the communications field, particularly a kind of Spectrum Conversion method based on polarized rotation angular speed discretization coding.
Background technology
Under existing communication system, use electromagnetic spectrum feature representation raw information, various signal madulation mode is adopted to carry out information transmission, the modulation technique wherein adopted comprises amplitude modulation double side band modulation (AM), carrier-suppressed double sideband amplitude-modulating modulation (DSB), AMSSB (Amplitude Modulation Single-Side-Band) modulation (SSB), vestigial sideband amplitudemodulation modulation (VSB), frequency modulation(FM) (FM), phase-modulation (PM), amplitude keying (ASK), frequency keying (FSK), phase keying (PSK, DPSK, QPSK, QAM, MSK), pulse-ampli tude modulation (PAM), pulse-width modulation (PWM), pulse-position modulation (PPM), pulse-code modulation (PCM) and differential pulse code modulation (DPCM) etc.
When adopting above various modulation system to realize information transmission, inromation bandwidth is limited by channel spectral bandwidth, and in order to improve the transmission rate of information system, needs to increase channel width.But along with spacecraft quantity is increasing, spatial frequency spectrum resource day is becoming tight, and unrestrictedly can not increase channel width.
Summary of the invention
The object of the invention is to overcome the deficiencies in the prior art, a kind of Spectrum Conversion method based on polarized rotation angular speed discretization coding is provided, the method adopts different polarized rotation angular speed to modulate signal, and the value of polarized rotation angular speed is Periodic Rectangular impulse function, therefore the Energy distribution of the modulation signal of symbolic information is carried at the frequency place of setting, the energy level being distributed in the modulation signal outside this frequency is low, interference can not be caused to other communication systems, this signal taking frequency spectrum resource can be reduced.
Above-mentioned purpose of the present invention is achieved by following technical solution:
Based on a Spectrum Conversion method for polarized rotation angular speed discretization coding, comprise the following steps:
(1), according to system works centre frequency f 0generate the single frequency carrier signal S of linear polarization 0(t), wherein, S 0(t)=cos (ω 0t), ω 0=2 π f 0;
(2) the N kind symbolic information transmitted, as required determines N number of polarized rotation angular speed, wherein, and the n-th polarized rotation angular speed ‖ ω corresponding with n-th kind of symbolic information c,nt the computing formula of () ‖ is as follows:
ω c,n(t)=ω o 0≤t≤τ n
ω c,n(t)=0 τ n<t<τ n+T n
Wherein, Δ θ nfor setting the n-th polarized rotation step angle and 0 ° of < Δ θ n<360 °; T nfor the rotation stepping time interval corresponding to the n-th symbolic information of setting, T n>T 0; N-th polarized rotation angular speed ‖ ω c,nt () ‖ is ω c,nthe periodic extension function of (t), continuation cycle T n'=τ n+ T n; Wherein, n=1,2 ..., N;
(3), in the sending time slots of each symbolic information, the polarized rotation angular speed corresponding according to described symbolic information, adopt the polarized rotation angular speed of discretization to modulate carrier signal, wherein said symbolic information sending time slots length is T s, concrete modulations achievement is as follows:
If the symbolic information transmitted in described symbolic information sending time slots is n-th kind of symbolic information, then according to the n-th polarized rotation angular speed ‖ ω c,nt () ‖ modulates carrier signal, obtain the signal after modulating as follows:
S c,n(t)=cos(||ω c,n(t)||·t)·cos(ω o·t)+j·sin(||ω c,n(t)||·t)·cos(ω o·t);
(4), at symbolic information receiving terminal adopt reception antenna A and reception antenna B to carry out Signal reception respectively, obtain Received signal strength y (t)=y a(t)+jy b(t), wherein, y at () is the Received signal strength of reception antenna A, y bt () is the Received signal strength of reception antenna B, reception antenna A is horizontal linear polarization antenna, and reception antenna B is vertical line poliarizing antenna;
(5), in each receiving slot, Received signal strength y (t) obtained is received to step (4) and carries out Fourier transform, obtain Received signal strength frequency spectrum Y (w), and described signal spectrum Y (w) is searched for obtain frequency ω corresponding to Amplitude maxima maxif, described in or the symbolic information then received in described receiving slot is n-th kind of symbolic information.
In the above-mentioned Spectrum Conversion method based on polarized rotation angular speed discretization, in step (3), adopt the polarized rotation angular speed ‖ ω of discretization c,nt () ‖ is to carrier signal S 0(t)=cos (ω 0when t) modulating, produce short electromagnetic pulse interference S j,n(t), wherein, described impulse disturbances S j,nthe cycle T of (t) j,nn+ T n, pulse duration τ j,nn, pulse amplitude
In the above-mentioned Spectrum Conversion method based on polarized rotation angular speed discretization, corresponding to the rotation stepping time interval T of the n-th symbolic information n=k nt 0, wherein k nfor positive integer, n=1,2 ..., N.
In the above-mentioned Spectrum Conversion method based on polarized rotation angular speed discretization, described symbolic information sending time slots length is T smeet following condition: T s &GreaterEqual; 360 &Delta; &theta; n ( &tau; n + T n ) , n = 1,2 , . . . , N .
The present invention compared with prior art has following beneficial effect:
(1), the present invention adopts different polarized rotation angular speed to modulate signal, and the value of polarized rotation angular speed is Periodic Rectangular impulse function, therefore the main energetic carrying the modulation signal of symbolic information is distributed in the frequency place of setting, this signal taking frequency spectrum resource can be reduced, under the prerequisite that channel width is constant, improve the transmission rate of communication system, therefore can be used for the passback problem solving low rail high-resolution satellite magnanimity reconnaissance data;
(2), the present invention adopts different polarized rotation angular speed to modulate signal, different polarized rotation angular speed is converted to angular speed Doppler frequency shift at receiving terminal, signal energy is made to concentrate on the frequency place of offset carrier centre frequency, and the polarization characteristic of traditional electronic jamming signal does not change, its angular speed Doppler frequency shift is zero, interference signal is still positioned at centre carrier frequency frequency, therefore band pass filter can be utilized in the present invention effectively to distinguish the useful signal of interference signal and the offset carrier center frequency point being positioned at carrier center frequency, thus improve the antijamming capability of communication system.
Accompanying drawing explanation
Fig. 1 is the Spectrum Conversion method flow diagram based on angular velocity of rotation discretization coding of the present invention;
Fig. 2 is the n-th polarized rotation angular speed ‖ ω that in the present invention, n-th kind of symbolic information is corresponding c,nthe time domain waveform of (t) ‖;
Fig. 3 is the normalized power spectrum of signal after modulation in the present invention.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is described in further detail:
The present invention utilizes the polarization characteristic of electromagnetic signal to carry out information transmission, wherein often kind of symbolic information sets a polarized rotation angular speed function, namely the polarisation angles modulating rear signal is made to change according to the step angle of setting between 0 ~ 360 degree, thus the main energetic of signal after modulation is concentrated on carrier signal frequency place, and adopt two linear polarized antennas to receive at receiving terminal, different polarized rotation angular speed is converted to angular speed Doppler frequency shift at receiving terminal, make the Energy distribution of signal at the frequency place of offset carrier centre frequency, therefore symbolic information can be demodulated according to this values of Doppler frequency shift.
The inventive method flow chart as shown in Figure 1, the Spectrum Conversion method based on polarized rotation angular speed discretization coding of the present invention comprises the following steps:
(1), according to system works centre frequency f 0generate the single frequency carrier signal S of linear polarization 0(t), wherein, S 0(t)=cos (ω 0t), ω 0=2 π f 0;
(2) the N kind symbolic information transmitted, as required determines N number of polarized rotation angular speed, wherein, and the n-th polarized rotation angular speed ‖ ω corresponding with n-th kind of symbolic information c,nt the computing formula of () ‖ is as follows:
ω c,n(t)=ω o 0≤t≤τ n
ω c,n(t)=0 τ n<t<τ n+T n
Wherein, Δ θ nfor setting the n-th polarized rotation step angle and 0 ° of < Δ θ n<360 °; T nfor the rotation stepping time interval corresponding to the n-th symbolic information of setting, T n>T 0, can preferably this time interval be T n=k nt 0, wherein k nfor positive integer, optimal value is k n=1; N-th polarized rotation angular speed ‖ ω c,nt () ‖ is ω c,nthe periodic extension function of (t), continuation cycle T n'=τ n+ T n, this polarized rotation angular speed ‖ ω c,nthe time domain beamformer of (t) ‖ as shown in Figure 2, wherein, n=1,2 ..., N;
(3), in each symbol sending time slots, the polarized rotation angular speed corresponding according to described symbol, adopt the polarized rotation angular speed of discretization to modulate carrier signal, wherein said symbol sending time slots length is T s, this symbol sending time slots length is T smeet following condition: it is hereby ensured that the changing value of the polarisation angles of modulation signal in sending time slots interval is greater than 360 degree, namely polarisation angles can rotate a circle above;
Above-described concrete modulations achievement is as follows:
If the symbolic information transmitted in described symbol time slot is n-th kind of symbolic information, then according to the n-th polarized rotation angular speed ‖ ω c,nt () ‖ modulates carrier signal, obtain the signal after modulating as follows:
S c,n(t)=cos(||ω c,n(t)||·t)·cos(ω o·t)+j·sin(||ω c,n(t)||·t)·cos(ω o·t);
After modulation as shown in Figure 3, the normalized power spectrum of signal is known, signal S after this modulation c,nt the main energetic of () is distributed in center frequency point ω 0place.And in symbolic information sending time slots, signal S c,nt the polarisation angles of () is according to step angle Δ θ nchange, namely the value of this polarisation angles is respectively θ 0, θ 0+ Δ θ n, θ 0+ 2 Δ θ n..., wherein θ 0for initial polarization angle.
Short electromagnetic pulse interference S is created in above modulated process j,n(t), wherein, described impulse disturbances S j,nthe cycle T of (t) j,nn+ T n, pulse duration τ j,nn, pulse amplitude in order to reduce above short electromagnetic pulse interference to the impact of communication system, can according to pulse amplitude and step angle Δ θ nrelation, select suitable step angle Δ θ nvalue reduces pulse amplitude, makes this short electromagnetic pulse interference signal can not form interference to other communication systems.。
(4), at symbolic information receiving terminal adopt reception antenna A and reception antenna B to carry out Signal reception respectively, obtain Received signal strength y (t)=y a(t)+jy b(t), wherein, y at () is the Received signal strength of reception antenna A, y bt () is the Received signal strength of reception antenna B, reception antenna A is horizontal linear polarization antenna, and reception antenna B is vertical line poliarizing antenna;
Receiving terminal adopts two linear polarized antennas to receive, different polarized rotation angular speed is converted to different angular speed Doppler frequency shifts at receiving terminal, make the Energy distribution of signal at the frequency place of offset carrier centre frequency, therefore can demodulate symbolic information according to this values of Doppler frequency shift.
(5), in each receiving slot, Received signal strength y (t) obtained is received to step (4) and carries out Fourier transform, obtain Received signal strength frequency spectrum Y (w), and described signal spectrum Y (w) is searched for obtain frequency ω corresponding to Amplitude maxima maxif, described in or the symbolic information then received in described receiving slot is n-th kind of symbolic information.
Adopt Spectrum Conversion method of the present invention, the Spectrum Conversion method that multiple carrier frequency carries out above-described polarized rotation angular speed discretization can be chosen in communication band bandwidth, when can guarantee to adopt multiple carrier frequency parallel processing according to the performing step of above invention, can not spectral interference be caused, and can compatible existing communication system.The polarized rotation angle of m user should be greater than 2 π .m2 Δ f m, wherein 2 Δ f mit is the communication channel bandwidth of m user.
The above; be only the embodiment of the best of the present invention, but protection scope of the present invention is not limited thereto, is anyly familiar with those skilled in the art in the technical scope that the present invention discloses; the change that can expect easily or replacement, all should be encompassed within protection scope of the present invention.
The content be not described in detail in specification of the present invention belongs to the known technology of professional and technical personnel in the field.

Claims (4)

1., based on a Spectrum Conversion method for polarized rotation angular speed discretization coding, it is characterized in that comprising the following steps:
(1), according to system works centre frequency f 0generate the single frequency carrier signal S of linear polarization 0(t), wherein, S 0(t)=cos (ω 0t), ω 0=2 π f 0;
(2) the N kind symbolic information transmitted, as required determines N number of polarized rotation angular speed, wherein, the n-th polarized rotation angular speed corresponding with n-th kind of symbolic information || ω c,n(t) || computing formula as follows:
ω c,n(t)=ω o 0≤t≤τ n
ω c,n(t)=0 τ n<t<τ n+T n
Wherein, Δ θ nfor setting the n-th polarized rotation step angle and 0 ° of < Δ θ n<360 °; T nfor the rotation stepping time interval corresponding to the n-th symbolic information of setting, T n>T 0; N-th polarized rotation angular speed || ω c,n(t) || be ω c,nthe periodic extension function of (t), continuation cycle T n'=τ n+ T n; Wherein, n=1,2 ..., N;
(3), in the sending time slots of each symbolic information, the polarized rotation angular speed corresponding according to described symbolic information, adopt the polarized rotation angular speed of discretization to modulate carrier signal, wherein said symbolic information sending time slots length is T s, concrete modulations achievement is as follows:
If the symbolic information transmitted in described symbolic information sending time slots is n-th kind of symbolic information, then according to the n-th polarized rotation angular speed || ω c,n(t) || carrier signal is modulated, obtains the signal after modulating as follows:
S c,n(t)=cos(||ω c,n(t)||·t)·cos(ω o·t)+j·sin(||ω c,n(t)||·t)·cos(ω o·t);
(4), at symbolic information receiving terminal adopt reception antenna A and reception antenna B to carry out Signal reception respectively, obtain Received signal strength y (t)=y a(t)+jy b(t), wherein, y at () is the Received signal strength of reception antenna A, y bt () is the Received signal strength of reception antenna B, reception antenna A is horizontal linear polarization antenna, and reception antenna B is vertical line poliarizing antenna;
(5), in each receiving slot, Received signal strength y (t) obtained is received to step (4) and carries out Fourier transform, obtain Received signal strength frequency spectrum Y (w), and described signal spectrum Y (w) is searched for obtain frequency ω corresponding to Amplitude maxima maxif, described in &omega; max = &omega; 0 + 2 &pi;&tau; n ( &tau; n + T n ) T 0 Or &omega; max = &omega; 0 - 2 &pi;&tau; n ( &tau; n + T n ) T 0 , The symbolic information then received in described receiving slot is n-th kind of symbolic information.
2. a kind of Spectrum Conversion method based on polarized rotation angular speed discretization coding according to claim 1, is characterized in that, in step (3), adopts the polarized rotation angular speed of discretization || ω c,n(t) || to carrier signal S 0(t)=cos (ω 0when t) modulating, produce short electromagnetic pulse interference S j,n(t), wherein, described impulse disturbances S j,nthe cycle T of (t) j,nn+ T n, pulse duration τ j,nn, pulse amplitude | r &RightArrow; n | = 2 &CenterDot; 1 - cos ( &Delta;&theta; n ) .
3. a kind of Spectrum Conversion method based on polarized rotation angular speed discretization coding according to claim 1, is characterized in that, corresponding to the rotation stepping time interval T of the n-th symbolic information n=k nt 0, wherein k nfor positive integer, n=1,2 ..., N.
4. a kind of Spectrum Conversion method based on polarized rotation angular speed discretization coding according to claim 1, it is characterized in that, described symbolic information sending time slots length is T smeet following condition: T s &GreaterEqual; 360 &Delta; &theta; n ( &tau; n + T n ) , n = 1,2 , . . . , N .
CN201410586362.3A 2014-10-27 2014-10-27 A kind of Spectrum Conversion method based on polarized rotation angular speed discretization coding Expired - Fee Related CN104301271B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410586362.3A CN104301271B (en) 2014-10-27 2014-10-27 A kind of Spectrum Conversion method based on polarized rotation angular speed discretization coding

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410586362.3A CN104301271B (en) 2014-10-27 2014-10-27 A kind of Spectrum Conversion method based on polarized rotation angular speed discretization coding

Publications (2)

Publication Number Publication Date
CN104301271A true CN104301271A (en) 2015-01-21
CN104301271B CN104301271B (en) 2017-11-28

Family

ID=52320843

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410586362.3A Expired - Fee Related CN104301271B (en) 2014-10-27 2014-10-27 A kind of Spectrum Conversion method based on polarized rotation angular speed discretization coding

Country Status (1)

Country Link
CN (1) CN104301271B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111525960A (en) * 2019-02-01 2020-08-11 华为技术有限公司 Quantum communication method, device and system
CN112310640A (en) * 2020-09-21 2021-02-02 西安电子科技大学 high-Q antenna assembly and signal modulation method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101296208A (en) * 2008-04-21 2008-10-29 上海大学 Two-bit-per-symbol small linear frequency-modulation key modulation communication method
CN101741805A (en) * 2010-01-29 2010-06-16 华为技术有限公司 Signal modulating method and device as well as signal modulating system
CN102457459A (en) * 2010-10-16 2012-05-16 中兴通讯股份有限公司 Complex carrier modulation method, device and system
EP2549706A1 (en) * 2010-10-16 2013-01-23 ZTE Corporation Complex carrier demodulation method and device

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101296208A (en) * 2008-04-21 2008-10-29 上海大学 Two-bit-per-symbol small linear frequency-modulation key modulation communication method
CN101741805A (en) * 2010-01-29 2010-06-16 华为技术有限公司 Signal modulating method and device as well as signal modulating system
CN102457459A (en) * 2010-10-16 2012-05-16 中兴通讯股份有限公司 Complex carrier modulation method, device and system
EP2549706A1 (en) * 2010-10-16 2013-01-23 ZTE Corporation Complex carrier demodulation method and device

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘海洋: "《一种新的卫星地面干扰源查找方法》", 《中国无线电》 *
张柳: "《基于相位载波技术的全光纤分布式传感系统及其稳定性研究》", 《中国优秀硕士学位论文全文数据库(电子期刊)》 *
李安平: "《基于GEO卫星通信中EIRP值的无线电测量及分析方法》", 《电波科学学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111525960A (en) * 2019-02-01 2020-08-11 华为技术有限公司 Quantum communication method, device and system
CN111525960B (en) * 2019-02-01 2022-01-14 华为技术有限公司 Quantum communication method, device and system
US11476950B2 (en) 2019-02-01 2022-10-18 Huawei Technologies Co., Ltd. Quantum communication method, apparatus, and system
CN112310640A (en) * 2020-09-21 2021-02-02 西安电子科技大学 high-Q antenna assembly and signal modulation method

Also Published As

Publication number Publication date
CN104301271B (en) 2017-11-28

Similar Documents

Publication Publication Date Title
US8965290B2 (en) Amplitude enhanced frequency modulation
CN101729471A (en) Composite transmission communication method of analogue signal and digital signal of broadcast communication system
CN101588332B (en) Cosine signal and Chirp signal combined signal modulating and demodulating method, and signal transmitting and receiving method based on the method
CN102368758A (en) New improvement scheme about gaussian minimum shift keying (GMSK) modulation technology
CN104301271B (en) A kind of Spectrum Conversion method based on polarized rotation angular speed discretization coding
CN104717165A (en) Biorthogonal bipolar binary shift pulse keying modulation and demodulation method
CN104333525A (en) GMSK (Gaussian minimum shift keying) modulating system synchronization method
CN102014092B (en) Cascade-mode-based quaternary minimum-shift keying ( MSK) modulation method and device
CN105991508A (en) 8QAM modulation method and system thereof
CN103812810B (en) Quaternary Continuous Phase Modulation demodulation method
CN108983155A (en) A kind of radar-communication integration waveform design method
CN101026602A (en) Orthogonal modulated chaos communication method
CN108123724A (en) A kind of communication system based on shortwave narrowband waveform
CN101729195B (en) Analogue signal and digital information combined transmission communication method
CN101753495A (en) Very minimum chirp keying ultra-narrow bandwidth broadcast communication modulation method
CN101710888B (en) Compound signal transmission and communication method based on equal-amplitude equal-period modulated carrier wave technology
CN203014854U (en) Offset quadrature phase shift keying signal transmitter
CN103618688A (en) Multi-carrier-phase modulation and relevant demodulation method
CN103346856B (en) A kind of bpsk signal transmission method based on spacing wave
US20220078068A1 (en) Dual-modulation transmission in a wireless communication system
CN105119855B (en) A kind of modulation and demodulation method and device
Udawant et al. Digital image processing by using GMSK
Ikwuazom et al. Simulation Analysis of Information Transmitting Mechanism Based on Amplitude Modulation Technique
CN103001903B (en) Modulation and demodulation method of wave form synergy signal based on frequency domain and fractional Fourier domain orthogonal basis function
CN103795435B (en) A kind of image frequency suppressing method and device

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171128

Termination date: 20181027