CN110161492B - 舰船航向航速提取方法 - Google Patents

舰船航向航速提取方法 Download PDF

Info

Publication number
CN110161492B
CN110161492B CN201910066973.8A CN201910066973A CN110161492B CN 110161492 B CN110161492 B CN 110161492B CN 201910066973 A CN201910066973 A CN 201910066973A CN 110161492 B CN110161492 B CN 110161492B
Authority
CN
China
Prior art keywords
radar
radars
speed
doppler frequency
measured
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910066973.8A
Other languages
English (en)
Other versions
CN110161492A (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.)
Beijing Research Institute of Mechanical and Electrical Technology
Original Assignee
Beijing Research Institute of Mechanical and Electrical Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Research Institute of Mechanical and Electrical Technology filed Critical Beijing Research Institute of Mechanical and Electrical Technology
Priority to CN201910066973.8A priority Critical patent/CN110161492B/zh
Publication of CN110161492A publication Critical patent/CN110161492A/zh
Application granted granted Critical
Publication of CN110161492B publication Critical patent/CN110161492B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/588Velocity or trajectory determination systems; Sense-of-movement determination systems deriving the velocity value from the range measurement
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/62Sense-of-movement determination

Landscapes

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

Abstract

本发明涉及精确制导技术领域,公开了一种舰船航向航速提取方法。其中,该方法包括:基于双基地分布的两个雷达分别距舰船的距离RT、RR和两个雷达之间的距离L计算双基地角β,两个雷达包括第一雷达和第二雷达;基于两个雷达的工作波长和测得的多普勒频率以及双基地角β计算航速预估值
Figure DDA0001956015470000011
和航向预估值φ0;计算航速预估值
Figure DDA0001956015470000012
和航向预估值φ0点处的雅可比矩阵G和多普勒频率预估值向量F0;基于雅可比矩阵G和多普勒频率预估值向量F0计算航速真实值vg和航向真实值φ。由此,可以实现舰船目标航向、航速的快速准确提取。

Description

舰船航向航速提取方法
技术领域
本发明涉及精确制导技术领域,尤其涉及一种舰船航向航速提取方法。
背景技术
对海作战环境下,航向、航速是舰船目标的重要特征,快速准确的目标航向、航速特征提取对舰船编队识别、战场态势构建等具有重要的意义。传统基于目标的位置及其变化信息的航向、航速舰船目标特征提取,受雷达测量精度特别是测角精度影响较大,而且往往需要相对较长的时间进行滤波来提高精度,在现代反舰导弹对海作战中的应用受到限制。
发明内容
本发明的目的在于克服现有技术不足,提供了一种舰船航向航速提取方法,能够解决现有技术中精度低且时间长的问题。
本发明的技术解决方案:一种舰船航向航速提取方法,其中,该方法包括:
基于双基地分布的两个雷达分别距舰船的距离RT、RR和两个雷达之间的距离L计算双基地角β,两个雷达包括第一雷达和第二雷达;
基于两个雷达的工作波长和测得的多普勒频率以及双基地角β计算航速预估值
Figure BDA0001956015450000011
和航向预估值φ0
计算航速预估值
Figure BDA0001956015450000012
和航向预估值φ0点处的雅可比矩阵G和多普勒频率预估值向量F0
基于雅可比矩阵G和多普勒频率预估值向量F0计算航速真实值vg和航向真实值φ。
优选地,通过下式基于双基地分布的两个雷达分别距舰船的距离RT、RR和两个雷达之间的距离L计算双基地角β:
Figure BDA0001956015450000021
优选地,通过下式中的任意两个方程基于两个雷达的工作波长和测得的多普勒频率以及双基地角β计算航速预估值
Figure BDA0001956015450000022
和航向预估值φ0
Figure BDA0001956015450000023
其中,fd1、λ1分别为所述第一雷达测得的多普勒频率和工作波长,fd2、λ2分别为第二雷达测得的多普勒频率和工作波长,fd12、λ12分别为测得的双基地多普勒频率和双基地雷达工作波长。
优选地,通过下式计算航速预估值
Figure BDA0001956015450000024
和航向预估值φ0点处的雅可比矩阵G和多普勒频率预估值向量F0
Figure BDA0001956015450000025
Figure BDA0001956015450000031
优选地,通过下式基于雅可比矩阵G和多普勒频率预估值向量F0计算航速真实值vg和航向真实值φ:
Figure BDA0001956015450000032
其中,
Figure BDA0001956015450000033
F为测得的多普勒频率向量且F=(fd1,fd2,fd12)T,N为雷达多普勒测量误差协方差矩阵。
通过上述技术方案,可以根据双基地分布的两个雷达分别距舰船的距离和雷达之间的距离计算双基地角,进而可以根据两个雷达的工作波长及雷达测得的多普勒频率以及计算得到的双基地角得到航速vg、航向φ的预估值
Figure BDA0001956015450000034
φ0,从而可以计算得到航速和航向预估值点处的雅可比矩阵G和多普勒频率预估值向量F0,最终可以基于雅可比矩阵G和多普勒频率预估值向量F0计算得到航速真实值vg和航向真实值φ。由此,可以实现舰船目标航向、航速的快速准确提取。
附图说明
所包括的附图用来提供对本发明实施例的进一步的理解,其构成了说明书的一部分,用于例示本发明的实施例,并与文字描述一起来阐释本发明的原理。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种舰船航向航速提取方法的流程图;
图2为本发明实施例中双基地雷达位置关系示意图。
具体实施方式
下面将结合附图对本发明的具体实施例进行详细说明。在下面的描述中,出于解释而非限制性的目的,阐述了具体细节,以帮助全面地理解本发明。然而,对本领域技术人员来说显而易见的是,也可以在脱离了这些具体细节的其它实施例中实践本发明。
在此需要说明的是,为了避免因不必要的细节而模糊了本发明,在附图中仅仅示出了与根据本发明的方案密切相关的设备结构和/或处理步骤,而省略了与本发明关系不大的其他细节。
图1为本发明实施例提供的一种舰船航向航速提取方法的流程图。
在本发明一种实施例中,可以通过双基地分布式的两部雷达测量的关于舰船目标距离及多普勒信息解算出舰船目标的航向、航速信息。两部雷达双基地分布可以如图2所示。
图2为本发明实施例中双基地雷达位置关系示意图。
如图2所示,可以假设第一雷达(如图2中Tx)和第二雷达2(如图2中Rx)到舰船目标(如图2所示Tgt)的距离分别为RT和RR,而两雷达之间距离即双基地基线为L,则关于舰船目标的双基地角可以为
Figure BDA0001956015450000041
若舰船目标的航速为vg、航向(航向角)为φ,则有:
Figure BDA0001956015450000042
其中:fd1、λ1分别为第一雷达测量的多普勒频率和工作波长;fd2、λ2分别为第二雷达测量的多普勒频率和工作波长;fd12、λ12分别为测得的双基地多普勒频率和工作波长;εk,k=1,2,3分别为多普勒测量误差。
此外,图2中的θT和θR均为北向夹角。
根据上述雷达测量的多普勒信息并利用最小二乘法可实现对舰船目标的航向、航速信息的快速准确提取。
下面将结合图1具体描述本发明实施例提供的一种舰船目标的航向、航速信息提取方法。
如图1所示,本发明实施例提供了一种舰船航向航速提取方法,其中,该方法可以包括:
S100,基于双基地分布的两个雷达分别距舰船的距离RT、RR和两个雷达之间的距离L计算双基地角β,两个雷达包括第一雷达和第二雷达;
也就是,可以根据两部雷达测量的舰船目标距离和基线长度,即RT、RR及L计算双基地角β。
S102,基于两个雷达的工作波长和测得的多普勒频率以及双基地角β计算航速预估值
Figure BDA0001956015450000051
和航向预估值φ0
其中,两个雷达测得的多普勒频率例如可以表示为:
F=(fd1,fd2,fd12)T, (2)
S104,计算航速预估值
Figure BDA0001956015450000052
和航向预估值φ0点处(即,
Figure BDA0001956015450000053
处)的雅可比矩阵G和多普勒频率预估值向量F0
S106,基于雅可比矩阵G和多普勒频率预估值向量F0计算航速真实值vg和航向真实值φ。
举例来讲,可以根据泰勒公式并利用最小二乘方法基于雅可比矩阵G和多普勒频率预估值向量F0计算舰船目标的航速真实值vg和航向真实值φ。
通过上述技术方案,可以根据双基地分布的两个雷达分别距舰船的距离和雷达之间的距离计算双基地角,进而可以根据两个雷达的工作波长及雷达测得的多普勒频率以及计算得到的双基地角得到航速vg、航向φ的预估值
Figure BDA0001956015450000061
φ0,从而可以计算得到航速和航向预估值点处的雅可比矩阵G和多普勒频率预估值向量F0,最终可以基于雅可比矩阵G和多普勒频率预估值向量F0计算得到航速真实值vg和航向真实值φ。由此,可以实现舰船目标航向、航速的快速准确提取。
根据本发明一种实施例,可以通过下式基于双基地分布的两个雷达分别距舰船的距离RT、RR和两个雷达之间的距离L计算双基地角β:
Figure BDA0001956015450000062
根据本发明一种实施例,可以通过下式中的任意两个方程基于两个雷达的工作波长和测得的多普勒频率以及双基地角β计算航速预估值
Figure BDA0001956015450000063
和航向预估值φ0
Figure BDA0001956015450000064
其中,fd1、λ1分别为所述第一雷达测得的多普勒频率和工作波长,fd2、λ2分别为第二雷达测得的多普勒频率和工作波长,fd12、λ12分别为测得的双基地多普勒频率和双基地雷达工作波长。
也就是,可以根据式(4)中的任意两个方程求解航速预估值
Figure BDA0001956015450000065
和航向预估值φ0。通过式(4)得到的预估值是方程解出来的近似预估值。
根据本发明一种实施例,通过下式计算航速预估值
Figure BDA0001956015450000066
和航向预估值φ0点处的雅可比矩阵G和多普勒频率预估值向量F0
Figure BDA0001956015450000071
Figure BDA0001956015450000072
同理,通过式(6)得到的多普勒频率预估值是方程解出来的近似预估值。
根据本发明一种实施例,可以通过下式基于雅可比矩阵G和多普勒频率预估值向量F0计算航速真实值vg和航向真实值φ:
Figure BDA0001956015450000073
其中,
Figure BDA0001956015450000074
F为测得的多普勒频率向量且F=(fd1,fd2,fd12)T,N为雷达多普勒测量误差协方差矩阵。
举例来讲,雷达多普勒测量误差协方差矩阵N可以根据雷达统计结果预先确定。
下面结合实例对本发明实施例所述的方法进行描述。
在本实例中,描述了如何运用本发明上述实施例所述的方法通过对舰船目标的多普勒测量来提取舰船目标的航向、航速信息。
假设两个雷达到目标(例如,舰船)的距离分别为RT=10km和RR=10km,而两雷达之间距离即基线长度L=6km,则通过式(3)得到相应的双基地角
Figure BDA0001956015450000075
假设舰船目标的航向、航速分别为vg=30m/s、φ=π/3=60°,若两个雷达工作频率均为f=5GHz,则相应的目标真实的多普勒频率(即,不存在误差的时候,对应于式(1)中的ε1、ε2和ε3均为0的情况)为:
Figure BDA0001956015450000081
假设雷达多普勒频率测量误差协方差矩阵为:
Figure BDA0001956015450000082
而由于雷达多普勒测量误差导致雷达测得的多普勒频率(即,存在误差的情况)为:
Figure BDA0001956015450000083
利用数值法解非线性方程组(即,解式(4)中任意两个方程):
Figure BDA0001956015450000084
求解得到航速真实值vg和航向真实值φ的预估值
Figure BDA0001956015450000085
φ0=58.87°。
然后,计算在
Figure BDA0001956015450000086
处的雅可比矩阵G和多普勒频率预估值向量F0
Figure BDA0001956015450000087
Figure BDA0001956015450000091
最后,根据泰勒公式并利用最小二乘方法可以计算得到解舰船目标的航速真实值、航向真实值:
Figure BDA0001956015450000092
由此可见,利用本发明实施例所述的方法,可以实现舰船目标航向、航速的快速准确提取。
如上针对一种实施例描述和/或示出的特征可以以相同或类似的方式在一个或更多个其它实施例中使用,和/或与其它实施例中的特征相结合或替代其它实施例中的特征使用。
应该强调,术语“包括/包含”在本文使用时指特征、整件、步骤或组件的存在,但并不排除一个或更多个其它特征、整件、步骤、组件或其组合的存在或附加。
本发明以上的方法可以由硬件实现,也可以由硬件结合软件实现。本发明涉及这样的计算机可读程序,当该程序被逻辑部件所执行时,能够使该逻辑部件实现上文所述的装置或构成部件,或使该逻辑部件实现上文所述的各种方法或步骤。本发明还涉及用于存储以上程序的存储介质,如硬盘、磁盘、光盘、DVD、flash存储器等。
这些实施例的许多特征和优点根据该详细描述是清楚的,因此所附权利要求旨在覆盖这些实施例的落入其真实精神和范围内的所有这些特征和优点。此外,由于本领域的技术人员容易想到很多修改和改变,因此不是要将本发明的实施例限于所例示和描述的精确结构和操作,而是可以涵盖落入其范围内的所有合适修改和等同物。
本发明未详细说明部分为本领域技术人员公知技术。

Claims (5)

1.一种舰船航向航速提取方法,其特征在于,该方法包括:
基于双基地分布的两个雷达分别距舰船的距离RT、RR和两个雷达之间的距离L计算双基地角β,两个雷达包括第一雷达和第二雷达;
基于两个雷达的工作波长和测得的多普勒频率以及双基地角β计算航速预估值
Figure FDA0002685765880000011
和航向预估值φ0
计算航速预估值
Figure FDA0002685765880000012
和航向预估值φ0点处的雅可比矩阵G和多普勒频率预估值向量F0
基于雅可比矩阵G和多普勒频率预估值向量F0计算航速真实值vg和航向真实值φ,
其中,通过下式计算多普勒频率预估值向量F0
Figure FDA0002685765880000013
其中,λ1为所述第一雷达测得的工作波长,λ2为第二雷达测得的工作波长,λ12为测得的双基地雷达工作波长。
2.根据权利要求1所述的方法,其特征在于,通过下式基于双基地分布的两个雷达分别距舰船的距离RT、RR和两个雷达之间的距离L计算双基地角β:
Figure FDA0002685765880000014
3.根据权利要求2所述的方法,其特征在于,通过下式中的任意两个方程基于两个雷达的工作波长和测得的多普勒频率以及双基地角β计算航速预估值
Figure FDA0002685765880000021
和航向预估值φ0
Figure FDA0002685765880000022
其中,fd1、λ1分别为所述第一雷达测得的多普勒频率和工作波长,fd2、λ2分别为第二雷达测得的多普勒频率和工作波长,fd12、λ12分别为测得的双基地多普勒频率和双基地雷达工作波长。
4.根据权利要求3所述的方法,其特征在于,通过下式计算航速预估值
Figure FDA0002685765880000023
和航向预估值φ0点处的雅可比矩阵G:
Figure FDA0002685765880000024
5.根据权利要求4所述的方法,其特征在于,通过下式基于雅可比矩阵G和多普勒频率预估值向量F0计算航速真实值vg和航向真实值φ:
Figure FDA0002685765880000031
其中,
Figure FDA0002685765880000032
F为测得的多普勒频率向量且F=(fd1,fd2,fd12)T,N为雷达多普勒测量误差协方差矩阵。
CN201910066973.8A 2019-01-24 2019-01-24 舰船航向航速提取方法 Active CN110161492B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910066973.8A CN110161492B (zh) 2019-01-24 2019-01-24 舰船航向航速提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910066973.8A CN110161492B (zh) 2019-01-24 2019-01-24 舰船航向航速提取方法

Publications (2)

Publication Number Publication Date
CN110161492A CN110161492A (zh) 2019-08-23
CN110161492B true CN110161492B (zh) 2020-12-08

Family

ID=67644816

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910066973.8A Active CN110161492B (zh) 2019-01-24 2019-01-24 舰船航向航速提取方法

Country Status (1)

Country Link
CN (1) CN110161492B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103698759A (zh) * 2013-12-27 2014-04-02 武汉大学 一种基于单频网的外辐射源雷达系统及其信号处理方法
CN106597428A (zh) * 2016-12-20 2017-04-26 中国航空工业集团公司雷华电子技术研究所 一种海面目标航向航速估算方法
CN108089183A (zh) * 2017-11-28 2018-05-29 西安电子科技大学 一种用于异步多基地雷达系统的检测跟踪一体化方法
EP3330732A1 (de) * 2016-12-02 2018-06-06 Thales Deutschland GmbH Verfahren und verarbeitungseinheit einer anordnung zur flugsicherung und/oder flugleitung von luftfahrzeugen
CN108363054A (zh) * 2018-02-07 2018-08-03 电子科技大学 用于单频网络和多路径传播的被动雷达多目标跟踪方法
CN108398678A (zh) * 2018-02-08 2018-08-14 电子科技大学 一种集中式mimo雷达快速自适应目标跟踪方法
CN108919231A (zh) * 2018-05-25 2018-11-30 长江大学 联合误差条件下双基地mimo雷达角度估计方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3301470A3 (en) * 2016-09-29 2018-06-20 Panasonic Corporation Multi-radar system
CN107783078B (zh) * 2017-09-11 2020-09-22 西北大学 一种波束-多普勒酉esprit多目标角度估计方法
CN108226920B (zh) * 2018-01-09 2021-07-06 电子科技大学 一种基于预测值处理多普勒量测的机动目标跟踪系统及方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103698759A (zh) * 2013-12-27 2014-04-02 武汉大学 一种基于单频网的外辐射源雷达系统及其信号处理方法
EP3330732A1 (de) * 2016-12-02 2018-06-06 Thales Deutschland GmbH Verfahren und verarbeitungseinheit einer anordnung zur flugsicherung und/oder flugleitung von luftfahrzeugen
CN106597428A (zh) * 2016-12-20 2017-04-26 中国航空工业集团公司雷华电子技术研究所 一种海面目标航向航速估算方法
CN108089183A (zh) * 2017-11-28 2018-05-29 西安电子科技大学 一种用于异步多基地雷达系统的检测跟踪一体化方法
CN108363054A (zh) * 2018-02-07 2018-08-03 电子科技大学 用于单频网络和多路径传播的被动雷达多目标跟踪方法
CN108398678A (zh) * 2018-02-08 2018-08-14 电子科技大学 一种集中式mimo雷达快速自适应目标跟踪方法
CN108919231A (zh) * 2018-05-25 2018-11-30 长江大学 联合误差条件下双基地mimo雷达角度估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Direction finding for bistatic MIMO radar with uniform circular array";Yunhe Cao等;《Electrical and electronic engineering》;20131231;第1-6页 *
"双基地雷达目标速度计算的FPGA实现";宋万均等;《电子技术应用》;20140226;第40卷(第1期);第47-49页 *

Also Published As

Publication number Publication date
CN110161492A (zh) 2019-08-23

Similar Documents

Publication Publication Date Title
CN110554376A (zh) 用于运载工具的雷达测程法
JP6597517B2 (ja) 物標検出装置
CN107102303B (zh) 机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法
US9846229B1 (en) Radar velocity determination using direction of arrival measurements
CN106990395B (zh) 一种雷达目标的校准方法及装置
CN107024687B (zh) 一种离线快速实现pos/激光雷达安装误差标定的方法
CN110361020B (zh) 用于确定车道线坐标的方法及系统
CN114755661A (zh) 一种移动激光扫描系统参数标定方法及标定装置
CN114063622B (zh) 无人船自主停泊定位方法、装置及相关组件
CN110006432B (zh) 一种基于几何先验信息下的室内机器人快速重定位的方法
CN110161492B (zh) 舰船航向航速提取方法
KR101173954B1 (ko) Sar 영상 형성을 위한 도플러 중심 주파수 정보 추정방법 및 그를 실행시키기 위한 프로그램을 기록한 컴퓨터 판독가능한 매체
CN102508217B (zh) 建立雷达测量误差标定模型的方法
CN112154429A (zh) 高精度地图定位方法、系统、平台及计算机可读存储介质
CN105301593A (zh) 一种基于双频rcs信息融合的飞行目标高度估计方法
RU2638177C1 (ru) Способ определения координат источника радиоизлучений с борта летательного аппарата по двум азимутальным пеленгам
KR102198298B1 (ko) 지상 또는 해상의 표적을 추적하는 항공기 탑재 레이다 장치 및 그 동작 방법
CN109375159B (zh) 纯方位加权约束总体最小二乘定位方法
WO2017168588A1 (ja) 測定装置、測定方法、及び、プログラム
CN111781591B (zh) 基于分数阶傅里叶变换的雷达目标速度快速估计方法
US10037035B2 (en) Methods and apparatus for positioning aircraft based on images of mobile targets
JP2017181195A (ja) 測定装置、測定方法、及び、プログラム
CN108872973B (zh) 一种弹道导弹目标定轨的ekf滤波方法
CN113538699A (zh) 基于三维点云的定位方法、装置、设备及存储介质
RU2684710C1 (ru) Система коррекции ошибок инс летательного аппарата по дорожной карте местности

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant