CN105137395B - 一种设定 fda 雷达的频率差的方法及装置 - Google Patents

一种设定 fda 雷达的频率差的方法及装置 Download PDF

Info

Publication number
CN105137395B
CN105137395B CN201510451891.7A CN201510451891A CN105137395B CN 105137395 B CN105137395 B CN 105137395B CN 201510451891 A CN201510451891 A CN 201510451891A CN 105137395 B CN105137395 B CN 105137395B
Authority
CN
China
Prior art keywords
frequency difference
value
cost function
difference matrix
control variable
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
CN201510451891.7A
Other languages
English (en)
Other versions
CN105137395A (zh
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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201510451891.7A priority Critical patent/CN105137395B/zh
Publication of CN105137395A publication Critical patent/CN105137395A/zh
Application granted granted Critical
Publication of CN105137395B publication Critical patent/CN105137395B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种设定FDA雷达的频率差的方法及装置，该方法包括：构建步骤：根据FDA雷达的模糊函数构建包含有频率差的代价函数；以及搜索步骤：搜索频率差，将使得代价函数最小的频率差设定为FDA雷达的频率差。

Description

一种设定FDA雷达的频率差的方法及装置
技术领域
本发明涉及FDA雷达技术领域，尤其涉及一种设定FDA雷达的频率差的方法及装置。
背景技术
频率分级阵列雷达(FDA，Frequency Diverse Radar)在同一时间对不同的阵元施加了不同的频率差，即：每个阵元发射出去的信号频率不相同。这些发射出去的信号在空间上相互叠加，会使FDA雷达波束在空间上呈现出某些距离角度位置上加强，其他距离角度位置上减弱的特性。
现有的FDA雷达的频率差通常是线性设定的，也就是说，通常会选取一个远小于FDA雷达的中心频率f0的频率差Δfm，第一个阵元的发射信号的中心频率为f0，第n个阵元的发射信号的中心频率为f0+(n-1)Δfm。
现有的一些技术文献研究了如何优化FDA雷达的发射信号。但是，还没有披露通过调整频率差来优化FDA雷达发射信号的技术方案。
发明内容
本发明的目的在于提供一种设定FDA雷达的频率差的方法及装置，使得FDA雷达的发射信号得到优化。
本发明的一个实施例提供了一种设定FDA雷达的频率差的方法，包括：构建步骤：根据FDA雷达的模糊函数构建包含有频率差的代价函数；以及搜索步骤：搜索频率差，将使得代价函数最小的频率差设定为FDA雷达的频率差。
本发明的另一个实施例提供了一种设定FDA雷达的频率差的装置，包括：构建模块，用于根据FDA雷达的模糊函数构建包含有频率差的代价函数；以及搜索模块，用于搜索频率差，将使得代价函数最小的频率差设定为FDA雷达的频率差。
附图说明
通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。其中在附图中，参考数字之后的字母标记指示多个相同的部件，当泛指这些部件时，将省略其最后的字母标记。在附图中：
图1为本发明的设定FDA雷达的频率差的方法的一个实施例的流程图；
图2为图1所示方法中的步骤101的一个实施例的流程图；
图3为图1所示方法中的步骤102的一个实施例的流程图；
图4为本发明的设定FDA雷达的频率差的装置一个实施例的示意性框图。
在附图中，使用相同或类似的标号来指代相同或类似的元素。
具体实施方式
现在将参考附图来详细描述本发明的示例性实施方式。应当理解，附图中示出和描述的实施方式仅仅是示例性的，意在阐释本发明的原理和精神，而并非限制本发明的范围。
参考图1，图1为本发明的设定FDA雷达的频率差的方法的一个实施例100的流程图。图1所示的实施例100可以包含如下步骤101至102。
步骤101是构建步骤：根据FDA雷达的模糊函数构建包含有频率差的代价函数。
在本发明的一个实施例中，参考图2，步骤101可以包含如下子步骤201至203。
子步骤201是第一积分计算步骤：计算模糊函数的P范数在预先设定的最小距离角度范围内的第一积分。
在本发明的一个实施例中，模糊函数可以用如下公式(1)来表示：
其中，M表示FDA雷达中的阵元总数，exp表示以自然对数e为底的指数函数，T0表示一个固定的导向矢量时刻，τ表示延时，υ表示多普勒频移，θ表示第一方位角，r表示第一距离，θ′表示第二方位角，r′表示第二距离，c0表示光速，f0表示中心频率，Δfm表示第一频率差，Δfm′表示第二频率差。其中，*表示取共轭，fm＝f0+Δfm，φm(t)表示第m个阵元的发射信号的包络。
相应的，上述的模糊函数的P范数在预先设定的较小距离角度范围内的第一积分可以表示为：
其中，预先设定的较小距离角度范围是由rmin＜r＜rmax且θmin＜θ＜θmax所限定的范围，rmin表示较小距离范围中的距离最小值，rmax表示较小距离范围中的距离最大值，θmin表示较小角度范围中的角度最小值，θmax表示较小角度范围中的角度最大值。预先设定的较小距离角度范围可以是需要干扰或跟踪的目标所在的区域范围。
在本发明的一个实施例中，P范数中的P值可以等于1或者2。
子步骤202是第二积分计算步骤：计算模糊函数的P范数在预先设定的较大距离角度范围内的第二积分。
在本发明的一个实施例中，模糊函数的P范数在预先设定的较大距离角度范围内的第二积分可以表示为：
其中，预先设定的较大距离角度范围是由Rmin＜r＜Rmax且所限定的范围，Rmin表示较大距离范围中的距离最小值，Rmax表示较大距离范围中的距离最大值。预先设定的较大距离角度范围可以是FDA雷达能够干扰或跟踪的区域范围。
在本发明的一个实施例中，预先设定的较小距离角度范围是完全被包含在较大距离角度范围内的区域。
子步骤203是比值计算步骤：将第一积分与第二积分的比值作为代价函数。
在本发明的一个实施例中，代价函数可以用如下公式(2)来表示：
其中，C表示由各个阵元的频率差所组成的矩阵，即：频率差矩阵。
步骤102是搜索步骤：搜索频率差，将使得代价函数最小的频率差设定为FDA雷达的频率差。
在本发明的一个实施例中，可以在预先设定的较小距离角度范围内，搜索多组频率差，并根据这些频率差计算代价函数值，然后将代价函数值最小的频率差作为在该预先设定的较小距离角度范围内最优的频率差。
在本发明的一个实施例中，参考图3，步骤102可以包含如下子步骤301至308。
子步骤301是频率差矩阵产生步骤：在预先设定的元素取值范围内，随机选取多个互不相等的元素值以产生频率差矩阵并保存。
在本发明的一个实施例中，频率差矩阵C中的每一个元素的初始值可以从预先设定的范围{1,…,K}中随机选取。在本发明的一个实施例中，选取时还可以使得频率差矩阵中的任意两个元素的值均不相等。
子步骤302是第一代价函数值计算步骤：根据保存的频率差矩阵计算第一代价函数值。
在本发明的一个实施例中，可以将保存的频率差矩阵C代入上述的公式(2)以计算得到第一代价函数值Fp(C)。
需要说明的是，当第一次执行子步骤302时，保存的频率差矩阵指的是经由子步骤301产生的初始的频率差矩阵。当非第一次执行子步骤302时，保存的频率差矩阵指的是经由子步骤304或305保存下来的频率差矩阵。
子步骤303是元素值更改步骤：更改保存的频率差矩阵中的任一元素值。
在本发明的一个实施例中，可以在保存的频率差矩阵中随机地选择一个元素，并从预先设定的范围{1,…,K}中随机选取一个值。由于子步骤303可能会被多次循环执行，矩阵中的某一个元素值也就有可能被不止一次地被修改，在本发明的一个实施例中，更改后的元素值与该元素在之前任一次的值均不相同。
子步骤304是第二代价函数值计算步骤：根据更改后的频率差矩阵计算第二代价函数值。
如果经过子步骤303更改以后的频率差矩阵用C'表示，则可以将C'代入上述的公式(2)并计算得到第二代价函数值Fp(C′)。
子步骤305是第一更新步骤：当第一代价函数值大于第二代价函数值时，将更改后的频率差矩阵作为保存的频率差矩阵。
当第一代价函数值大于第二代价函数值时，说明更改后的频率差C'是更为优化的选择，则可以将更改后的频率差矩阵作为保存的频率差矩阵并舍弃原有的频率差矩阵C。
子步骤306是第二更新步骤：当第一代价函数值小于第二代价函数值时，根据第一代价函数值与第二代价函数值的差值和控制变量计算接受概率值并根据接受概率值将更改后的频率差矩阵作为保存的频率差矩阵。
当第一代价函数值小于第二代价函数值时，为了防止一直停留于局部极小点，可以以一定的概率值决定是将更改后的频率差矩阵C'作为保存的频率差矩阵还是将更改前的频率差矩阵C作为保存的频率差矩阵。具体来说，可以通过如下方法：首先根据第一代价函数值与第二代价函数值的差值和控制变量计算接受概率值，具体来说，可以通过如下公式来计算接受概率值：
p＝exp(-ΔF/k*Tmp)，其中，p表示接受概率值，exp()表示以自然对数e为底的指数函数，ΔF表示第一代价函数值与第二代价函数值的差值，k表示波尔兹曼常数，Tmp表示控制变量，控制变量可以有一个初始值，并且，控制变量的值可以通过下面的子步骤307调整。
然后产生0到1之间的随机数a，若p＞a则将更改后的频率差矩阵C'作为保存的频率差矩阵，否则将更改前的频率差矩阵C作为保存的频率差矩阵C。
子步骤307是控制变量调整步骤：当在控制变量不变的情况下频率差矩阵被更改的次数达到预先设定的次数时，调整控制变量。
在本发明的一个实施例中，当在控制变量不变的情况下频率差矩阵被更改的次数达到预先设定的次数时，可以按照预先设计的调整策略，降低控制变量Tmp的值并进入子步骤308。
在本发明的一个实施例中，当在控制变量不变的情况下频率差矩阵被更改的次数小于预先设定的次数时，则可以保持控制变量不变并返回子步骤302并执行子步骤302到子步骤306。
上述的预先设定的次数可以是在性能和运算量/运算时间之间进行权衡后设定的。
子步骤308是频率差输出步骤：当控制变量小于预先设定的下限值时，输出保存的频率差矩阵以设定FDA雷达的频率差。
当经过子步骤307调整后的控制变量小于预先设定的下限值时，可以认为对频率差的搜索已经达到可以接受的程度，因此，可以用保存的频率差矩阵中的元素值对FDA雷达的频率差进行设定。
该预先设定的下限值可以是在性能和运算量/运算时间之间进行权衡后设定的。
在本发明的一个实施例中，当经过子步骤307调整后的控制变量大于预先设定的下限值时，可以基于调整后的控制变量返回子步骤302。
至此描述了根据本发明实施例的设定FDA雷达的频率差的方法。
本发明提出的设定FDA雷达的频率差的方法，能够通过选取适当的频率差来对FDA雷达的发射信号进行优化，从而使得其模糊函数在主瓣值不降的前提下旁瓣值降低。
与该方法类似，本发明还提供了相应的设定FDA雷达的频率差的装置。
图4所示为本发明的设定FDA雷达的频率差的装置一个实施例400的示意性框图。
如图所示，装置400可以包括：构建模块401，用于根据FDA雷达的模糊函数构建包含有频率差的代价函数；以及搜索模块402，用于搜索频率差，将使得代价函数最小的频率差设定为FDA雷达的频率差。
在本发明的一个实施例中，构建模块401可以进一步包括：第一积分计算模块，用于计算模糊函数的P范数在预先设定的较小距离角度范围内的第一积分；第二积分计算模块，用于计算模糊函数的P范数在预先设定的较大距离角度范围内的第二积分；以及比值计算模块，用于将第一积分与第二积分的比值作为代价函数。
在本发明的一个实施例中，搜索模块402可以进一步包括：频率差矩阵产生模块，用于在预先设定的元素取值范围内，随机选取多个互不相等的元素值以产生频率差矩阵并保存；第一代价函数值计算模块，用于根据保存的频率差矩阵计算第一代价函数值；元素值更改模块，用于更改保存的频率差矩阵中的任一元素值；第二代价函数值计算模块，用于根据更改后的频率差矩阵计算第二代价函数值；第一更新模块，用于当第一代价函数值大于第二代价函数值时，将更改后的频率差矩阵作为保存的频率差矩阵；第二更新模块，用于当第一代价函数值小于第二代价函数值时，根据第一代价函数值与第二代价函数值的差值和控制变量计算接受概率值并根据接受概率值将更改后的频率差矩阵作为保存的频率差矩阵；控制变量调整模块，用于当在控制变量不变的情况下频率差矩阵被更改的次数达到预先设定的次数时，调整控制变量；以及频率差输出模块，用于当控制变量小于预先设定的下限值时，输出保存的频率差矩阵以设定FDA雷达的频率差。
至此描述了根据本发明实施例的设定FDA雷达的频率差的装置。
本发明提出的设定FDA雷达的频率差的装置，能够通过选取适当的频率差来对FDA雷达发射信号进行优化，从而使得其模糊函数在主瓣值不降的前提下旁瓣值降低。

Claims (8)

1.一种设定FDA雷达的频率差的方法，其特征是，包括：
构建步骤：根据所述FDA雷达的模糊函数构建包含有所述频率差的
代价函数；以及
搜索步骤：搜索所述频率差，将使得所述代价函数最小的频率差设定为所述FDA雷达的频率差；
其中，所述构建步骤进一步包括：
第一积分计算步骤：计算所述模糊函数的P范数在预先设定的较小距离角度范围内的第一积分；
第二积分计算步骤：计算所述模糊函数的P范数在预先设定的较大距离角度范围内的第二积分；以及
比值计算步骤：将所述第一积分与所述第二积分的比值作为所述代价函数。
2.根据权利要求1所述的方法，其特征是，所述搜索步骤进一步包括：
频率差矩阵产生步骤：在预先设定的元素取值范围内，随机选取多个互不相等的元素值以产生频率差矩阵并保存；
第一代价函数值计算步骤：根据保存的频率差矩阵计算第一代价函数值；
元素值更改步骤：更改保存的频率差矩阵中的任一元素值；
第二代价函数值计算步骤：根据更改后的频率差矩阵计算第二代价函数值；
第一更新步骤：当第一代价函数值大于第二代价函数值时，将更改后的频率差矩阵作为保存的频率差矩阵；
第二更新步骤：当第一代价函数值小于第二代价函数值时，根据第一代价函数值与第二代价函数值的差值和控制变量计算接受概率值并根据所述接受概率值将更改后的频率差矩阵作为保存的频率差矩阵，所述控制变量具有一个初始值；
控制变量调整步骤：当在所述控制变量不变的情况下所述频率差矩阵被更改的次数达到预先设定的次数时，调整所述控制变量；以及频率差输出步骤：当所述控制变量小于预先设定的下限值时，输出保存的频率差矩阵以设定FDA雷达的频率差。
3.根据权利要求2所述的方法，其特征是，所述搜索步骤进一步包括：当在所述控制变量不变的情况下所述频率差矩阵被更改的次数小于预先设定的次数时，返回所述第一代价函数值计算步骤。
4.根据权利要求2所述的方法，其特征是，所述搜索步骤进一步包括：当所述控制变量大于预先设定的下限值时，返回所述第一代价函数值计算步骤。
5.根据权利要求2所述的方法，其特征是，在所述元素值更改步骤中，更改后的元素值与该元素在之前任一次的值均不相同。
6.根据权利要求2所述的方法，其特征是，所述接受概率值根据公式p＝exp(-ΔF/k*Tmp)计算得到，其中，p表示所述接受概率值，exp()表示以自然对数e为底的指数函数，ΔF表示第一代价函数值与第二代价函数值的差值，k表示波尔兹曼常数，Tmp表示所述控制变量。
7.一种设定FDA雷达的频率差的装置，其特征是，包括：
构建模块，用于根据所述FDA雷达的模糊函数构建包含有所述频率差的代价函数；以及
搜索模块，用于搜索所述频率差，将使得所述代价函数最小的频率差设定为所述FDA雷达的频率差；
其中，所述构建模块进一步包括：
第一积分计算模块，用于计算所述模糊函数的P范数在预先设定的较小距离角度范围内的第一积分；
第二积分计算模块，用于计算所述模糊函数的P范数在预先设定的较大距离角度范围内的第二积分；以及
比值计算模块，用于将所述第一积分与所述第二积分的比值作为所述代价函数。
8.根据权利要求7所述的装置，其特征是，所述搜索模块进一步包括：
频率差矩阵产生模块，用于在预先设定的元素取值范围内，随机选取多个互不相等的元素值以产生频率差矩阵并保存；
第一代价函数值计算模块，用于根据保存的频率差矩阵计算第一代价函数值；
元素值更改模块，用于更改保存的频率差矩阵中的任一元素值；
第二代价函数值计算模块，用于根据更改后的频率差矩阵计算第二代价函数值；
第一更新模块，用于当第一代价函数值大于第二代价函数值时，将更改后的频率差矩阵作为保存的频率差矩阵；
第二更新模块，用于当第一代价函数值小于第二代价函数值时，根据第一代价函数值与第二代价函数值的差值和控制变量计算接受概率值并根据所述接受概率值将更改后的频率差矩阵作为保存的频率差矩阵，所述控制变量具有一个初始值；
控制变量调整模块，用于当在所述控制变量不变的情况下所述频率差矩阵被更改的次数达到预先设定的次数时，调整所述控制变量；以及
频率差输出模块，用于当所述控制变量小于预先设定的下限值时，输出保存的频率差矩阵以设定FDA雷达的频率差。
CN201510451891.7A 2015-07-29 2015-07-29 一种设定 fda 雷达的频率差的方法及装置 Expired - Fee Related CN105137395B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510451891.7A CN105137395B (zh) 2015-07-29 2015-07-29 一种设定 fda 雷达的频率差的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510451891.7A CN105137395B (zh) 2015-07-29 2015-07-29 一种设定 fda 雷达的频率差的方法及装置

Publications (2)

Publication Number Publication Date
CN105137395A CN105137395A (zh) 2015-12-09
CN105137395B true CN105137395B (zh) 2017-06-20

Family

ID=54722801

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510451891.7A Expired - Fee Related CN105137395B (zh) 2015-07-29 2015-07-29 一种设定 fda 雷达的频率差的方法及装置

Country Status (1)

Country Link
CN (1) CN105137395B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108761414B (zh) * 2018-05-25 2022-05-20 电子科技大学 一种基于频控阵的s形干扰波束测试方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103592635A (zh) * 2013-11-11 2014-02-19 电子科技大学 一种认知fda雷达的信号发射方法、装置及雷达

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7319427B2 (en) * 2005-01-12 2008-01-15 The United States Of America As Represented By The Secretary Of The Air Force Frequency diverse array with independent modulation of frequency, amplitude, and phase
US7535408B2 (en) * 2007-08-31 2009-05-19 Lockheed Martin Corporation Apparatus and methods for detection of multiple targets within radar resolution cell

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103592635A (zh) * 2013-11-11 2014-02-19 电子科技大学 一种认知fda雷达的信号发射方法、装置及雷达

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Frequency Divere Array Radar With Time-Dependnt Frequency Offset;Waseem Khan et.al;《IEEE ANTENNAS AND WIRELESS PROPAGATION LETIERS》;20140410;第758-761页 *
Frequency Diverse Array Radar With Logarithmically Increasing Frequency Offset;Waseem Khan et.al;《IEEE ANTENNAS AND WIRELESS PROGAGATION LETTERS》;20141110;第499-502页 *
Impacts of frequency increment errors on frequency diverse array beampattern;Kuandong Gao et.al;《EURASIP Journal on Advances in Signal Processing》;20150409;全文 *
Precisely Beam Steering for Frequency Diverse Arrays Based on Frequency Offset Selection;Long Zhuang et.al;《RADAR》;20100325;全文 *

Also Published As

Publication number Publication date
CN105137395A (zh) 2015-12-09

Similar Documents

Publication Publication Date Title
US11245441B2 (en) Radar device
CN109495189A (zh) 一种阵列天线校准方法及装置
CN108141256B (zh) 天线阵列的相位调整方法及装置
CN112511209B (zh) 基于频率分集的阵列天线的波束控制方法、系统及控制器
US20230036414A1 (en) Neural network acceleration circuit and method
CN108562897A (zh) 一种mimo穿墙雷达的结构稀疏成像方法和装置
CN106683130A (zh) 一种深度图像获得方法及装置
CN109633526A (zh) 基于方向函数的非均匀圆阵相位干涉仪测向解模糊方法
CN112382855B (zh) 基于稀疏阵列的有源相控阵天线波束控制方法
CN104515979A (zh) 车辆用探知传感器补正装置及其方法
CN105137395B (zh) 一种设定 fda 雷达的频率差的方法及装置
EP3836454B1 (en) Method and device for improving phase measurement accuracy
CN105737806A (zh) 基于星载微波辐射计反演海面盐度的方法和装置
JP2021005861A (ja) フェーズドアレイアンテナを用いた高速空間探索
CN106374224A (zh) 电磁波成像系统及天线阵列信号校正方法
CN108614235B (zh) 一种多鸽群信息交互的单快拍测向方法
CN110806578A (zh) 一种波束控制方法及装置、可读存储介质
WO2022040424A1 (en) Rf scene generation simulation with external maritime surface
CN105116383B (zh) 产生fda雷达的发射波束的方法、装置及fda雷达
CN110890932B (zh) 5g毫米波信道信息确定的方法、相关装置及存储介质
KR102360811B1 (ko) 영상 레이다의 prf를 선정하는 장치 및 방법
US11983888B2 (en) Image change detection device and image change detection method
JP2000171551A (ja) マルチビーム・レーダ装置
EP3879301A1 (en) Radar target simulation system and method
CN114980032B (zh) 基于感通一体的v2v波束跟踪方法

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: 20170620

Termination date: 20180729