CN111505580B - 基于方位角和多普勒信息的多平台协同目标定位方法 - Google Patents

基于方位角和多普勒信息的多平台协同目标定位方法 Download PDF

Info

Publication number
CN111505580B
CN111505580B CN202010292402.9A CN202010292402A CN111505580B CN 111505580 B CN111505580 B CN 111505580B CN 202010292402 A CN202010292402 A CN 202010292402A CN 111505580 B CN111505580 B CN 111505580B
Authority
CN
China
Prior art keywords
target
platform
doppler
azimuth
spectrum
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
CN202010292402.9A
Other languages
English (en)
Other versions
CN111505580A (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering 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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN202010292402.9A priority Critical patent/CN111505580B/zh
Publication of CN111505580A publication Critical patent/CN111505580A/zh
Application granted granted Critical
Publication of CN111505580B publication Critical patent/CN111505580B/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
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/18Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

基于方位角和多普勒信息的多平台协同目标定位方法,涉及水声目标定位跟踪领域。本发明是为了解决现有缺少通过获得目标运动速度和航向来准确定位目标的方法的问题。本发明利用多个平台获取的目标方位和多普勒频移,被动实现运动目标初始位置、运动速度和航向的估计实现对目标的定位。它用于定位目标。

Description

基于方位角和多普勒信息的多平台协同目标定位方法
技术领域
本发明涉及目标定位方法。属于水声目标定位跟踪领域。
背景技术
目标位置信息是海洋水声目标的重要参数,亦是被动声呐设备关心的重要对象。方位交汇作为一种常用的被动目标定位方法,其定位误差会随着目标活动区域的扩大而迅速增大,其使用区域范围有所限制。
现有的文献中曾出现基于目标方位角的被动定位方法,现总结如下:
文献1:公开号为CN107797091A,主题名称为一种基于子空间的新型纯方位目标定位方法的公开文件中介绍了,利用“信号”与“噪声”子空间正交的思想,在经典纯方位目标定位的基础上,构建了新型纯方位定位迭代方程来提高目标定位精度。其研究是依然属于纯方位定位技术,得到的目标参数仅限于目标位置。
文献2(杨长虹.基于目标方位角测量的目标定位方法研究[D].东南大学硕士论文,2003.):利用纯方位角度测量值及传感器位置实现对目标初始位置和速度的估计。其研究仅利用了目标方位信息,得到的目标参数则限于目标位置和速度。
然而现有缺少通过获得目标运动速度和航向来准确定位目标的方法。
发明内容
本发明是为了解决现有缺少通过获得目标运动速度和航向来准确定位目标的方法的问题。现提供基于方位角和多普勒信息的多平台协同目标定位方法。
基于方位角和多普勒信息的多平台协同目标定位方法,所述方法包括以下步骤:
步骤一、每个平台将每个时刻采集到的声学信号分别进行低通滤波和快速傅里叶变换,经过低通滤波获得环境噪声连续谱,经过快速傅里叶变换获得信号谱,在每个时刻下将每个平台得到的信号谱和环境噪声连续谱做差,将差值高于6dB的信号谱保留,否则舍弃,从保留的信号谱中选取多个平台在同一时刻下均存在的信号谱,从选取的每个信号谱中得到多个线谱分量;
步骤二、在同一时刻下,将步骤一中任意一个平台的一个线谱分量与其他平台中的任意一个线谱分量均做差,对获得的所有差值中同时小于等于设定的多普勒变化范围的线谱分量进行保留,
步骤三、对步骤二筛选出的同一时刻每个平台的线谱分量,将同一频率范围内的两个线谱分量做差,所述两个线谱分量分别位于两个平台内,得到每个时刻下每两个平台间的多普勒频移;
步骤四、利用步骤三中得到的多普勒频移和步骤二中保留的每个线谱分量对应的方位角,构建方位角-多普勒超定方程组;
步骤五、对步骤四中的方位角-多普勒超定方程组进行最小二乘估计,得到目标初始位置坐标和目标速度;
步骤六、根据目标速度,得到目标径向运动速度和航向角,从而根据得到的目标初始位置坐标、目标径向运动速度和航向角,实现对目标的定位。
本发明的有益效果为:
本发明涉及UUV、滑翔机等水下移动平台利用方位角和多普勒信息进行多平台协同的目标定位方法。本发明利用多个平台获取的目标方位和多普勒频移,被动实现运动目标初始位置、运动速度和航向的估计。
基于方位角和多普勒信息的多平台协同目标定位方法不同于文献1和文献2,本申请在利用目标方位角同时,融合了平台间的多普勒频移,进而利用最小二乘优化思想将目标定位问题转化为超定线性方程组求解问题来提升目标初始位置、运动速度和航向的估计稳健性。
附图说明
图1为具体实施方式一所述的基于方位角和多普勒信息的多平台协同目标定位方法的流程图;
图2为三个平台协同定位目标位置的结构图;
图3为仿真中平台位置与目标位置的关系图;
图4为目标预测轨迹和目标真实运动轨迹的对比曲线图。
具体实施方式
具体实施方式一:参照图1具体说明本实施方式,本实施方式所述的基于方位角和多普勒信息的多平台协同目标定位方法,所述方法包括以下步骤:
步骤一、每个平台将每个时刻采集到的声学信号分别进行低通滤波和快速傅里叶变换,经过低通滤波获得环境噪声连续谱,经过快速傅里叶变换获得信号谱,在每个时刻下将每个平台得到的信号谱和环境噪声连续谱做差,将差值高于6dB的信号谱保留,否则舍弃,从保留的信号谱中选取多个平台在同一时刻下均存在的信号谱,从选取的每个信号谱中得到多个线谱分量;
步骤二、在同一时刻下,将步骤一中任意一个平台的一个线谱分量与其他平台中的任意一个线谱分量均做差,对获得的所有差值中同时小于等于设定的多普勒变化范围的线谱分量进行保留,
步骤三、对步骤二筛选出的同一时刻每个平台的线谱分量,将同一频率范围内的两个线谱分量做差,所述两个线谱分量分别位于两个平台内,得到每个时刻下每两个平台间的多普勒频移;
步骤四、利用步骤三中得到的多普勒频移和步骤二中保留的每个线谱分量对应的方位角,构建方位角-多普勒超定方程组;
步骤五、对步骤四中的方位角-多普勒超定方程组进行最小二乘估计,得到目标初始位置坐标和目标速度;
步骤六、根据目标速度,得到目标径向运动速度和航向角,从而根据得到的目标初始位置坐标、目标径向运动速度和航向角,实现对目标的定位。
本实施方式中,设定的多普勒变化范围为δf,假设平台的数量为3个,三个平台在第一时刻均得到一个线谱分量,分别为:f1j1、f2j1、f3j1,若|f1m1-f2n1|≤δf和|f1m1-f3n1|≤δf同时成立,则认为该线谱存在。
在数据融合中心,计算第一时刻时平台1和平台2间的多普勒频移Δf12j1、平台1和平台3间的多普勒频移Δf13j1、平台2和平台3间的多普勒频移Δf23j1,即:
Δf12j1=f1j1-f2j1
Δf13j1=f1j1-f3j1
Δf23j1=f2j1-f3j1
本实施方式中,步骤二中,将一个平台中的一个线谱分量分别与其他平台中的一个线谱分量做差,判断差值的绝对值是否同时满足小于等于设定的多普勒变化范围,该内容具体实现的例子为:
假设平台包括1号平台、2号平台和3号平台,经过步骤一后1号平台得到两个线谱分量值分别为1和2,
2号平台得到两个线谱分量值分别为3和4,
3号平台得到两个线谱分量值分别为5和6,
要比较1-3的差值绝对值和1-5的差值绝对值是否同时满足小于等于设定的多普勒变化范围,如果同时满足,则留下值为1、3和5的线谱分量,如果不满足,则舍弃1、3和5;按照这种方式在比较2-3的差值绝对值与2-5的差值绝对值是否同时满足小于等于设定的多普勒变化范围,如果满足则留下值为2、3和5的线谱分量,如果不满足,则舍弃2、3和5;通过这种比较方式,最终从三个平台中选出满足范围且数量相同的线谱分量。
本实时方式中,步骤三中的频率范围为预设的频率范围。低通滤波可以采用双向α滤波器实现。
本申请中一共有N个时刻,相邻两个时刻之间的时间间隔为△T,经过步骤一至步骤三得到N个△T时间段的线谱方位及平台间的多普勒频移观测序列。如果有3个平台的话,其中的多普勒频移观测序列指N个时刻的平台1和2间多普勒频移历程△f12j1、△f12j2、△f12j3、……△f12jN,平台1和3间多普勒频移历程△f13j1、△f13j2、△f13j3、……△f13jN,平台2和3间多普勒频移历程△f23j1、△f23j2、△f23j3、……△f23jN
本申请的每个平台可以采用四元十字阵实现。
具体实施方式二:本实施方式是对具体实施方式一所述的基于方位角和多普勒信息的多平台协同目标定位方法作进一步说明,本实施方式中,步骤一中,得到信号谱的过程为:
利用快速傅里叶变换将每个平台采集的4路声信号pi1(t)、pi2(t)、pi3(t)和pi4(t)转换成对应的频域信号Pi1(f)、Pi2(f)、Pi3(f)和Pi4(f),利用公式1,对频域信号Pi1(f)、Pi2(f)、Pi3(f)和Pi4(f)进行计算,得到信号谱Pi(f):
Pi(f)=|Pi1(f)|2+|Pi2(f)|2+|Pi3(f)|2+|Pi4(f)|2 公式1。
具体实施方式三:本实施方式是对具体实施方式一所述的基于方位角和多普勒信息的多平台协同目标定位方法作进一步说明,本实施方式中,步骤四中所述方位角-多普勒超定方程组为:
Figure BDA0002450921040000041
式中,p=(x0,y0,vx,vy),x0为目标初始位置的横坐标,y0为目标初始位置的纵坐标,vx为目标运动速度的x向分量,vy为目标的运动速度的y向分量,
Figure BDA0002450921040000042
Figure BDA0002450921040000043
分别表示方位角-多普勒伪线性表达式的测量矩阵和测量向量,
Figure BDA0002450921040000044
Figure BDA0002450921040000045
具体表达式分别为:
Figure BDA0002450921040000051
其中,
Figure BDA0002450921040000052
Figure BDA0002450921040000053
Figure BDA0002450921040000054
Figure BDA0002450921040000055
Figure BDA0002450921040000056
Figure BDA0002450921040000057
Figure BDA0002450921040000061
Figure BDA0002450921040000062
θajk表示第a个平台在第k时刻接收到的第j个线谱分量相对于大地北的角度,k=1,2,…,N,xa和ya分别表示第a个平台的x轴和y轴的坐标位置,ΔfabjN表示第a平台和第b平台在第N时刻的多普勒频移,a,b表示两个平台的编号,a,b=1,2,3,且a≠b,j=1,2,3,j表示线谱编号,△T为相邻两个时刻之间的时间间隔;c为水中的声速;f为信号谱的中心频率。
具体实施方式四:本实施方式是对具体实施方式一所述的基于方位角和多普勒信息的多平台协同目标定位方法作进一步说明,本实施方式中,步骤五中,对步骤四中的方位角-多普勒超定方程组进行最小二乘估计,得到p的估计值
Figure BDA0002450921040000063
Figure BDA0002450921040000064
根据公式3,得到x0,y0和vx,vy
具体实施方式五:本实施方式是对具体实施方式四所述的基于方位角和多普勒信息的多平台协同目标定位方法作进一步说明,本实施方式中,根据公式4,得到目标径向运动速度v:
Figure BDA0002450921040000071
具体实施方式六:本实施方式是对具体实施方式五所述的基于方位角和多普勒信息的多平台协同目标定位方法作进一步说明,本实施方式中,根据公式5,得到目标航向角
Figure BDA0002450921040000073
Figure BDA0002450921040000072
通过仿真实验验证用本申请的方法得到目标位置的准确性:
仿真中,每个平台中安装有正四面体,且正四面体的基元间距为35cm;目标辐射线谱信号频率为311Hz;信噪比为20dB;以10kHz采样率连续采集。
如3给出了平台位置与目标位置设定,三个平台的位置分别位于S1=[300,200]T,S2=[-200,100]T,S3=[400,-500]T;未知目标的初始位置为(x0,y0)=[400,500]T,以5.5m/s的速度向着北偏东120°方向匀速航行,即p=[400,500,4.7631,-2.75]T。图4给出了仿真目标运动轨迹与本发明估计得到的目标估计对比,其中,实线表示目标真实的运动轨迹,“◇”表示本发明根据估计结果预测的目标运动轨迹。这里估计得到的目标速度为5.4784m/s,航向角度为120.2052°,与实际的目标速度和运动方向相吻合。因此,本申请能够准确的定位目标。

Claims (5)

1.基于方位角和多普勒信息的多平台协同目标定位方法,其特征在于,所述方法包括以下步骤:
步骤一、每个平台将每个时刻采集到的声学信号分别进行低通滤波和快速傅里叶变换,经过低通滤波获得环境噪声连续谱,经过快速傅里叶变换获得信号谱,在每个时刻下将每个平台得到的信号谱和环境噪声连续谱做差,将差值高于6dB的信号谱保留,否则舍弃,从保留的信号谱中选取多个平台在同一时刻下均存在的信号谱,从选取的每个信号谱中得到多个线谱分量;
步骤二、在同一时刻下,将步骤一中任意一个平台的一个线谱分量与其他平台中的任意一个线谱分量均做差,对获得的所有差值中同时小于等于设定的多普勒变化范围的线谱分量进行保留;
步骤三、对步骤二筛选出的同一时刻每个平台的线谱分量,将同一频率范围内的两个线谱分量做差,所述两个线谱分量分别位于两个平台内,得到每个时刻下每两个平台间的多普勒频移;
步骤四、利用步骤三中得到的多普勒频移和步骤二中保留的每个线谱分量对应的方位角,构建方位角-多普勒超定方程组;
步骤五、对步骤四中的方位角-多普勒超定方程组进行最小二乘估计,得到目标初始位置坐标和目标速度;
步骤六、根据目标速度,得到目标径向运动速度和航向角,从而根据得到的目标初始位置坐标、目标径向运动速度和航向角,实现对目标的定位;
步骤四中所述方位角-多普勒超定方程组为:
Figure FDA0003406807190000011
式中,p=(x0,y0,vx,vy),x0为目标初始位置的横坐标,y0为目标初始位置的纵坐标,vx为目标运动速度的x向分量,vy为目标的运动速度的y向分量,
Figure FDA0003406807190000012
Figure FDA0003406807190000013
分别表示方位角-多普勒伪线性表达式的测量矩阵和测量向量,
Figure FDA0003406807190000014
Figure FDA0003406807190000015
具体表达式分别为:
Figure FDA0003406807190000021
其中,
Figure FDA0003406807190000022
Figure FDA0003406807190000023
Figure FDA0003406807190000024
Figure FDA0003406807190000025
Figure FDA0003406807190000026
Figure FDA0003406807190000027
Figure FDA0003406807190000031
Figure FDA0003406807190000032
θajk表示第a个平台在第k时刻接收到的第j个线谱分量相对于大地北的角度,k=1,2,…,N,xa和ya分别表示第a个平台的x轴和y轴的坐标位置,ΔfabjN表示第a平台和第b平台在第N时刻的多普勒频移,a,b表示两个平台的编号,a,b=1,2,3,且a≠b,j=1,2,3,j表示线谱编号,△T为相邻两个时刻之间的时间间隔;c为水中的声速;f为信号谱的中心频率。
2.根据权利要求1所述基于方位角和多普勒信息的多平台协同目标定位方法,其特征在于,步骤一中,得到信号谱的过程为:
利用快速傅里叶变换将每个平台采集的4路声信号pi1(t)、pi2(t)、pi3(t)和pi4(t)转换成对应的频域信号Pi1(f)、Pi2(f)、Pi3(f)和Pi4(f),利用公式1,对频域信号Pi1(f)、Pi2(f)、Pi3(f)和Pi4(f)进行计算,得到信号谱Pi(f):
Pi(f)=|Pi1(f)|2+|Pi2(f)|2+|Pi3(f)|2+|Pi4(f)|2 公式1。
3.根据权利要求1所述基于方位角和多普勒信息的多平台协同目标定位方法,其特征在于,步骤五中,对步骤四中的方位角-多普勒超定方程组进行最小二乘估计,得到p的估计值
Figure FDA0003406807190000033
Figure FDA0003406807190000041
根据公式3,得到x0,y0和vx,vy
4.根据权利要求3所述基于方位角和多普勒信息的多平台协同目标定位方法,其特征在于,根据公式4,得到目标径向运动速度v:
Figure FDA0003406807190000042
5.根据权利要求4所述基于方位角和多普勒信息的多平台协同目标定位方法,其特征在于,根据公式5,得到目标航向角
Figure FDA0003406807190000043
Figure FDA0003406807190000044
CN202010292402.9A 2020-04-14 2020-04-14 基于方位角和多普勒信息的多平台协同目标定位方法 Active CN111505580B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010292402.9A CN111505580B (zh) 2020-04-14 2020-04-14 基于方位角和多普勒信息的多平台协同目标定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010292402.9A CN111505580B (zh) 2020-04-14 2020-04-14 基于方位角和多普勒信息的多平台协同目标定位方法

Publications (2)

Publication Number Publication Date
CN111505580A CN111505580A (zh) 2020-08-07
CN111505580B true CN111505580B (zh) 2022-04-15

Family

ID=71872732

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010292402.9A Active CN111505580B (zh) 2020-04-14 2020-04-14 基于方位角和多普勒信息的多平台协同目标定位方法

Country Status (1)

Country Link
CN (1) CN111505580B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112763980B (zh) * 2020-12-28 2022-08-05 哈尔滨工程大学 一种基于方位角及其变化率的目标运动分析方法
CN114577213B (zh) * 2022-03-01 2022-11-18 哈尔滨工程大学 一种分布式多平台水下多目标关联与被动定位方法
CN115390077B (zh) * 2022-07-28 2023-06-06 哈尔滨工程大学 一种双基元条件下水下机动目标瞬时航速解析求解方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003001233A1 (de) * 2001-06-20 2003-01-03 Atlas Elektronik Gmbh Verfahren zum passiven bestimmen von zieldaten
CN105589066A (zh) * 2015-12-14 2016-05-18 西北工业大学 一种利用垂直矢量阵估计水下匀速运动航行器参数的方法
EP3088913A1 (fr) * 2015-04-29 2016-11-02 Zodiac Data Systems Procédé de suivi d'une cible mobile au moyen d'un modèle a priori
CN109540172A (zh) * 2018-12-27 2019-03-29 中国船舶重工集团公司第七0研究所 一种用于水雷平台的目标运动参数估计方法
CN110441769A (zh) * 2018-05-03 2019-11-12 北京航空航天大学 基于sar序贯图像的目标定位方法、装置及存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003001233A1 (de) * 2001-06-20 2003-01-03 Atlas Elektronik Gmbh Verfahren zum passiven bestimmen von zieldaten
EP3088913A1 (fr) * 2015-04-29 2016-11-02 Zodiac Data Systems Procédé de suivi d'une cible mobile au moyen d'un modèle a priori
CN105589066A (zh) * 2015-12-14 2016-05-18 西北工业大学 一种利用垂直矢量阵估计水下匀速运动航行器参数的方法
CN110441769A (zh) * 2018-05-03 2019-11-12 北京航空航天大学 基于sar序贯图像的目标定位方法、装置及存储介质
CN109540172A (zh) * 2018-12-27 2019-03-29 中国船舶重工集团公司第七0研究所 一种用于水雷平台的目标运动参数估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
单站无源定位技术研究;邱汇;《电子科技》;20140115;第78-80页 *
基于纯方位估计的远程被动定位技术研究;刘圣松;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20150115;第7-25页 *

Also Published As

Publication number Publication date
CN111505580A (zh) 2020-08-07

Similar Documents

Publication Publication Date Title
CN111505580B (zh) 基于方位角和多普勒信息的多平台协同目标定位方法
CN110749891B (zh) 一种可估计未知有效声速的自适应水下单信标定位方法
CN109669160B (zh) 一种水下瞬态声信号的检测方法
CN109489796A (zh) 一种基于单元辐射法的水下复杂结构辐射噪声源定位识别与声辐射预报方法
CN110738275A (zh) 基于ut-phd的多传感器序贯融合跟踪方法
CN110703205B (zh) 基于自适应无迹卡尔曼滤波的超短基线定位方法
CN114280541A (zh) 一种基于深海分布式垂直线列阵的目标被动定位方法
CN110196407B (zh) 一种基于频率预估的单矢量水听器信号来波方向估计方法
CN115828074A (zh) 一种辐射源目标定位参数联合估计方法
Han et al. Maneuvering target tracking using retrospective-cost input estimation
CN111781590B (zh) 基于高效frft的雷达目标参数估计方法
CN112666519A (zh) 一种基于广义二阶时延差的水下目标高精度定位方法
CN111856400B (zh) 一种水下目标声源定位方法及系统
CN114994608B (zh) 基于深度学习的多设备自组织麦克风阵列声源定位方法
RU2260197C2 (ru) Способ автоматического сопровождения маневрирующей цели в режиме активной локации гидроакустического или радиолокационного комплекса
CN110728026A (zh) 一种基于角速度量测的末端弹道目标被动跟踪方法
CN113219406B (zh) 一种基于扩展卡尔曼滤波的直接跟踪方法
CN110736988B (zh) 双基地pfa运动目标参数估计和成像方法
CN108981707B (zh) 基于时差量测箱粒子phd的被动跟踪多目标方法
CN109270344B (zh) 脉冲丢失下的相参脉冲信号频率估计方法
CN111220158A (zh) 一种基于时间方位历程图的线谱目标运动参数估计方法
CN116106910A (zh) 一种基于静止被动声纳多普勒的目标距离及轨迹估计方法
CN117031397B (zh) 一种运动物体噪声源定位和评估的快速计算方法
CN117908034B (zh) 基于自适应波束跟踪水下目标高稳健模基doa估计方法
Zhang et al. Research on trajectory smoothing optimization algorithm of underwater long baseline positioning system

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