CN109116297A - 一种被动雷达空间谱估计与合成波束的联合测向方法 - Google Patents

一种被动雷达空间谱估计与合成波束的联合测向方法 Download PDF

Info

Publication number
CN109116297A
CN109116297A CN201710487934.6A CN201710487934A CN109116297A CN 109116297 A CN109116297 A CN 109116297A CN 201710487934 A CN201710487934 A CN 201710487934A CN 109116297 A CN109116297 A CN 109116297A
Authority
CN
China
Prior art keywords
module
signal
estimation
spatial spectrum
airspace
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
CN201710487934.6A
Other languages
English (en)
Other versions
CN109116297B (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 Institute of Remote Sensing Equipment
Original Assignee
Beijing Institute of Remote Sensing Equipment
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 Institute of Remote Sensing Equipment filed Critical Beijing Institute of Remote Sensing Equipment
Priority to CN201710487934.6A priority Critical patent/CN109116297B/zh
Publication of CN109116297A publication Critical patent/CN109116297A/zh
Application granted granted Critical
Publication of CN109116297B publication Critical patent/CN109116297B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/80Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
    • G01S3/802Systems for determining direction or deviation from predetermined direction

Landscapes

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

Abstract

本发明公开了一种被动雷达空间谱估计与合成波束的联合测向方法,通过构建被动雷达空间谱估计与合成波束的联合测向系统,包括宽带接收天线阵(1)、接收通道(2)、信号处理单元(3)、信号提取模块、信号采集解算模块(5)、空域扫描模块(6)和空间谱估计模块(7);信号提取模块置于接受通道内,信号采集解算模块(5)、空域扫描模块(6)和空间谱估计模块(7)置于信号处理单元(3)内;宽带接收天线阵(1)接收目标辐射源的微波信号,信号经由接收通道(2)传递给信号处理单元(3),信号处理单元(3)完成对目标辐射源的测向。本发明比单纯采用谱估计进行全空域扫描的计算量降低了40%以上,能够在1ms内实现大空域搜索功能。

Description

一种被动雷达空间谱估计与合成波束的联合测向方法
技术领域
本发明涉及一种联合测向方法,特别是一种被动雷达空间谱估计与合成波束的联合测向方法。
背景技术
目前典型的微波被动导引头常用测向方法有比幅法测向、干涉仪测向、空间谱估计。上述测向方法均基于对电场分布的测量,从而计算出电波到达方向。其中比幅法测向是基于对目标幅度的测量,干涉仪测向是基于对目标相位差的测量,而空间谱估计是基于对目标信号电磁特征的多次测量。空间谱估计是利用天线阵列实现超分辨测向的一种测向体制,但其计算量大,在进行目标大空域搜索时,难以满足测向实时性的要求。
发明内容
本发明目的在于提供一种被动雷达空间谱估计与合成波束的联合测向方法,解决微波被动导引头大空域搜索的实时性问题。
一种被动雷达空间谱估计与合成波束的联合测向方法的具体步骤为:
第一步构建被动雷达空间谱估计与合成波束的联合测向系统
被动雷达空间谱估计与合成波束的联合测向系统,包括:宽带接收天线阵、接收通道、信号处理单元、信号提取模块、信号采集解算模块、空域扫描模块和空间谱估计模块;所述信号提取模块置于接受通道内,所述信号采集解算模块、空域扫描模块和空间谱估计模块置于信号处理单元内。
所述宽带接收天线阵包含N个天线单元,以宽带接收天线阵中心为坐标原点建立坐标系,以第i个天线为例,则第i个天线坐标表示为(xi,yi,zi);宽带接收天线阵与接收通道用微波电缆连接;接收通道的两个输出端口与信号处理单元的射频电缆和控制电缆连接;信号处理单元的通信接口与弹上惯导设备连接。
信号提取模块的功能为:提取宽带接收天线阵接收信号的幅度和相位;
信号采集解算模块的功能为:采集并解算接收通道的幅度和相位;
空域扫描模块的功能为:进行数字波束合成并进行空域扫描;
空间谱估计模块的功能为:进行空间谱估计计算,获取信号的入射方向。
第二步信号提取模块提取宽带接收天线阵的幅度和相位
信号提取模块根据宽带接收天线阵接收目标辐射源辐射的微波信号,经由接收通道下变频,频率为ω,由此得到宽带接收天线阵接收到的信号频率f0及各天线的幅度和相位。
第三步信号采集解算模块采集并解算接收通道的幅度和相位
信号处理单元采集接收通道输出的模拟信号,并根据计算得到通道的幅度和相位信息。第i个通道的幅度为Ei,相位为φi
第四步空域扫描模块进行数字波束合成并进行空域扫描
空域扫描模块进行数字波束合成和空域扫描处理,设来波方位角为α,来波俯仰角为β,则第i个通道与宽带接收天线阵中心的相位差为:
ψi=2πf0(xi sinβ+yi sinαcosβ+zi cosαcosβ)/c。
式中:c为光速。
又根据预先的通道功率补偿值,第i个通道预先的通道功率补偿值表示为 Ei',得到一组向量如下:
A=[E11,E22,…,Eii,…,ENN]T
根据得到的通道的幅度和相位信息,得到一组向量如下:
B=[E1φ1,E2φ2,…,Eiφi,…,ENφN]T
分别改变α、β的取值,使α、β的取值在全视场范围内以一定步长扫描,计算空域扫描值:
C=ABT
选取C值最大点对应的α、β的取值,记为αm、βm
第五步空间谱估计模块进行空间谱估计计算并获取信号的入射方向
空间谱估计模块根据得到的通道的幅度和相位信息,得到数据矩阵X如下,其中N为通道数,M为数据点数,
由X得到数据协方差矩阵R,
R=XXH
对R进行特征值分解得到R的噪声子空间导向矢量Un
用A与Un采用MUSIC算法对方位角αm、俯仰角βm附近空域进行搜索,搜索范围按实际情况确定,公式为:
式中:P为谱峰值。
找出极大值点对应的角度即信号的入射方向,完成测向。
至此,完成了被动雷达空间谱估计与合成波束联合测向。
本发明采用合成波束进行空域粗扫,空间谱估计在粗扫的结果上进行精扫,比单纯采用谱估计进行全空域扫描的计算量降低了40%以上,能够在1ms内实现大空域搜索功能。
附图说明
图1一种被动雷达空间谱估计与合成波束联合测向方法所述系统组成示意图;
图2一种被动雷达空间谱估计与合成波束联合测向方法所述宽带接收天线阵示意图。
1.宽带接收天线阵 2.接收通道 3.信号处理单元 4.信号提取模块 5.信号采集解算模块 6.空域扫描模块 7.空间谱估计模块
具体实施方式
一种被动雷达空间谱估计与合成波束的联合测向方法的具体步骤为:
第一步构建被动雷达空间谱估计与合成波束的联合测向系统
被动雷达空间谱估计与合成波束的联合测向系统,包括:宽带接收天线阵1、接收通道2、信号处理单元3、信号提取模块、信号采集解算模块5、空域扫描模块6和空间谱估计模块7;所述信号提取模块置于接受通道内,所述信号采集解算模块5、空域扫描模块6和空间谱估计模块7置于信号处理单元3内。
所述宽带接收天线阵1包含N个天线单元,以宽带接收天线阵1中心为坐标原点建立坐标系,以第i个天线为例,则第i个天线坐标表示为xi,yi,zi;宽带接收天线阵1与接收通道2用微波电缆连接;接收通道2的两个输出端口与信号处理单元3的射频电缆和控制电缆连接;信号处理单元3的通信接口与弹上惯导设备连接。
信号提取模块的功能为:提取宽带接收天线阵1接收信号的幅度和相位;
信号采集解算模块5的功能为:采集并解算接收通道2的幅度和相位;
空域扫描模块6的功能为:进行数字波束合成并进行空域扫描;
空间谱估计模块7的功能为:进行空间谱估计计算,获取信号的入射方向。
第二步信号提取模块4提取宽带接收天线阵1的幅度和相位
信号提取模块根据宽带接收天线阵1接收目标辐射源辐射的微波信号,经由接收通道2下变频,频率为ω,由此得到宽带接收天线阵1接收到的信号频率 f0及各天线的幅度和相位。
第三步信号采集解算模块5采集并解算接收通道2的幅度和相位
信号处理单元3采集接收通道2输出的模拟信号,并根据计算得到通道的幅度和相位信息。第i个通道的幅度为Ei,相位为φi
第四步空域扫描模块6进行数字波束合成并进行空域扫描
空域扫描模块6进行数字波束合成和空域扫描处理,设来波方位角为α,来波俯仰角为β,则第i个通道与宽带接收天线阵1中心的相位差为:
ψi=2πf0(xi sinβ+yi sinαcosβ+zi cosαcosβ)/c。
式中:c为光速。
又根据预先的通道功率补偿值,第i个通道预先的通道功率补偿值表示为 Ei',得到一组向量如下:
A=[E11,E22,…,Eii,…,ENN]T
根据得到的通道的幅度和相位信息,得到一组向量如下:
B=[E1φ1,E2φ2,…,Eiφi,…,ENφN]T
分别改变α、β的取值,使α、β的取值在全视场范围内以一定步长扫描,计算空域扫描值:
C=ABT
选取C值最大点对应的α、β的取值,记为αm、βm
第五步空间谱估计模块7进行空间谱估计计算并获取信号的入射方向
空间谱估计模块7根据得到的通道的幅度和相位信息,得到数据矩阵X如下,其中N为通道数,M为数据点数,
由X得到数据协方差矩阵R,
R=XXH
对R进行特征值分解得到R的噪声子空间导向矢量Un
用A与Un采用MUSIC算法对方位角αm、俯仰角βm附近空域进行搜索,搜索范围按实际情况确定,公式为:
式中:P为谱峰值。
找出极大值点对应的角度即信号的入射方向,完成测向。
至此,完成了被动雷达空间谱估计与合成波束联合测向。

Claims (3)

1.一种被动雷达空间谱估计与合成波束的联合测向方法,其特征在于具体步骤为:
第一步 构建被动雷达空间谱估计与合成波束的联合测向系统
被动雷达空间谱估计与合成波束的联合测向系统,包括:宽带接收天线阵(1)、接收通道(2)、信号处理单元(3)、信号提取模块、信号采集解算模块(5)、空域扫描模块(6)和空间谱估计模块(7);所述信号提取模块置于接受通道内,所述信号采集解算模块(5)、空域扫描模块(6)和空间谱估计模块(7)置于信号处理单元(3)内;
所述宽带接收天线阵(1)包含N个天线单元,以宽带接收天线阵(1)中心为坐标原点建立坐标系,以第i个天线为例,则第i个天线坐标表示为(xi,yi,zi);宽带接收天线阵(1)与接收通道(2)用微波电缆连接;接收通道(2)的两个输出端口与信号处理单元(3)的射频电缆和控制电缆连接;信号处理单元(3)的通信接口与弹上惯导设备连接;
信号提取模块的功能为:提取宽带接收天线阵(1)接收信号的幅度和相位;
信号采集解算模块(5)的功能为:采集并解算接收通道(2)的幅度和相位;
空域扫描模块(6)的功能为:进行数字波束合成并进行空域扫描;
空间谱估计模块(7)的功能为:进行空间谱估计计算,获取信号的入射方向;
第二步 信号提取模块(4)提取宽带接收天线阵(1)的幅度和相位
信号提取模块根据宽带接收天线阵(1)接收目标辐射源辐射的微波信号,经由接收通道(2)下变频,频率为ω,由此得到宽带接收天线阵(1)接收到的信号频率f0及各天线的幅度和相位;
第三步 信号采集解算模块(5)采集并解算接收通道(2)的幅度和相位
信号处理单元(3)采集接收通道(2)输出的模拟信号,并根据计算得到通道的幅度和相位信息;第i个通道的幅度为Ei,相位为φi
第四步 空域扫描模块(6)进行数字波束合成并进行空域扫描
空域扫描模块(6)进行数字波束合成和空域扫描处理,设来波方位角为α,来波俯仰角为β,根据预先的通道功率补偿值以及通道的幅度和相位信息,得到两组向量A和B,分别改变α、β的取值,使α、β的取值在全视场范围内以一定步长扫描,计算空域扫描值:
C=ABT
选取C值最大点对应的α、β的取值,记为αm、βm
第五步 空间谱估计模块(7)进行空间谱估计计算并获取信号的入射方向
根据得到的通道的幅度和相位信息,得到数据矩阵X,由X得到数据协方差矩阵R,
R=XXH
对R进行特征值分解得到R的噪声子空间导向矢量Un
用A与Un采用MUSIC算法对方位角αm、俯仰角βm附近空域进行搜索,搜索范围按实际情况确定,公式为:
式中:P为谱峰值;
找出极大值点对应的角度即信号的入射方向,完成测向;
至此,完成了被动雷达空间谱估计与合成波束联合测向。
2.如权利要求1所述的被动雷达空间谱估计与合成波束的联合测向方法,其特征在于:所述空域扫描模块(6)进行数字波束合成和空域扫描处理的具体过程为:
首先获取第i个通道与宽带接收天线阵中心的相位差ψi为:
ψi=2πf0(xisinβ+yisinαcosβ+zicosαcosβ)/c;
式中:c为光速;
又根据预先的通道功率补偿值,第i个通道预先的通道功率补偿值表示为Ei',得到一组向量为:
A=[E11,E22,…,Eii,…,ENN]T
根据得到的通道的幅度和相位信息,得到一组向量为:
B=[E1φ1,E2φ2,…,Eiφi,…,ENφN]T
分别改变α、β的取值,使α、β的取值在全视场范围内以一定步长扫描,计算空域扫描值为:
C=ABT
选取C值最大点对应的α、β的取值,记为αm、βm
3.如权利要求1或2所述的被动雷达空间谱估计与合成波束的联合测向方法,其特征在于:所述空间谱估计模块(7)进行空间谱估计计算并获取信号的入射方向时,数据矩阵X为:
其中N为通道数,M为数据点数。
CN201710487934.6A 2017-06-23 2017-06-23 一种被动雷达空间谱估计与合成波束的联合测向方法 Active CN109116297B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710487934.6A CN109116297B (zh) 2017-06-23 2017-06-23 一种被动雷达空间谱估计与合成波束的联合测向方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710487934.6A CN109116297B (zh) 2017-06-23 2017-06-23 一种被动雷达空间谱估计与合成波束的联合测向方法

Publications (2)

Publication Number Publication Date
CN109116297A true CN109116297A (zh) 2019-01-01
CN109116297B CN109116297B (zh) 2022-09-27

Family

ID=64732089

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710487934.6A Active CN109116297B (zh) 2017-06-23 2017-06-23 一种被动雷达空间谱估计与合成波束的联合测向方法

Country Status (1)

Country Link
CN (1) CN109116297B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109828232A (zh) * 2019-03-26 2019-05-31 北京遥感设备研究所 一种阵列天线测角模糊置信度计算装置及方法
CN110018362A (zh) * 2019-04-16 2019-07-16 中国电子科技集团公司第二十九研究所 一种主波束对称的宽带天线的相位中心测量方法
CN112068096A (zh) * 2020-06-09 2020-12-11 周口师范学院 一种被动雷达复合扫描三维空间测向方法及系统
CN113093167A (zh) * 2021-04-19 2021-07-09 成都安智杰科技有限公司 一种目标测向方法及装置
CN113203980A (zh) * 2021-04-20 2021-08-03 北京通广龙电子科技有限公司 高精度快速无线电测向方法及系统
CN117289037A (zh) * 2023-11-23 2023-12-26 南京华成微波技术有限公司 一种大功率相控阵天线平面近场测试方法及系统

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004088347A1 (fr) * 2003-03-31 2004-10-14 Da Tang Mobile Communications Equipment Co., Ltd. Procede d'estimation de la direction de l'arrivee d'un faisceau de fixation dans l'espace
CN1752771A (zh) * 2005-10-20 2006-03-29 武汉大学 单通道相控阵接收信号重构及空间信号处理方法
JP2008026035A (ja) * 2006-07-18 2008-02-07 Murata Mfg Co Ltd レーダ
CN101470187A (zh) * 2007-12-26 2009-07-01 中国科学院声学研究所 一种用于线列阵的高精度测向方法
CN103383452A (zh) * 2013-06-26 2013-11-06 西安电子科技大学 分布式阵列目标到达角估计方法
CN104459680A (zh) * 2014-12-03 2015-03-25 西安电子科技大学 Mimo雷达估计目标方向的快速方法
CN104749555A (zh) * 2014-12-19 2015-07-01 中国航天科技集团公司第五研究院第五一三研究所 一种组合相位差测向和空间谱测向的测向定位系统
US20150362580A1 (en) * 2014-06-11 2015-12-17 Agency For Defense Development Method for estimating angle of arrival of multi-target moving at high speed
CN105242264A (zh) * 2015-11-26 2016-01-13 上海无线电设备研究所 一种平面阵列多目标角度高分辨实现方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004088347A1 (fr) * 2003-03-31 2004-10-14 Da Tang Mobile Communications Equipment Co., Ltd. Procede d'estimation de la direction de l'arrivee d'un faisceau de fixation dans l'espace
CN1752771A (zh) * 2005-10-20 2006-03-29 武汉大学 单通道相控阵接收信号重构及空间信号处理方法
JP2008026035A (ja) * 2006-07-18 2008-02-07 Murata Mfg Co Ltd レーダ
CN101470187A (zh) * 2007-12-26 2009-07-01 中国科学院声学研究所 一种用于线列阵的高精度测向方法
CN103383452A (zh) * 2013-06-26 2013-11-06 西安电子科技大学 分布式阵列目标到达角估计方法
US20150362580A1 (en) * 2014-06-11 2015-12-17 Agency For Defense Development Method for estimating angle of arrival of multi-target moving at high speed
CN104459680A (zh) * 2014-12-03 2015-03-25 西安电子科技大学 Mimo雷达估计目标方向的快速方法
CN104749555A (zh) * 2014-12-19 2015-07-01 中国航天科技集团公司第五研究院第五一三研究所 一种组合相位差测向和空间谱测向的测向定位系统
CN105242264A (zh) * 2015-11-26 2016-01-13 上海无线电设备研究所 一种平面阵列多目标角度高分辨实现方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A. MONAKOV等: "Direction finding for an extended target with possibly non-symmetric spatial spectrum", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 *
H.H. NGUYEN等: "A direct band-pass implementation of subspace-based techniques for direction finding and digital beamforming applications", 《2003 IEEE WIRELESS COMMUNICATIONS AND NETWORKING》 *
张正鸿: "数字波束合成及测向技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑(月刊)》 *
甄雪娇: "基于空间谱估计的短波测向技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑(月刊)》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109828232A (zh) * 2019-03-26 2019-05-31 北京遥感设备研究所 一种阵列天线测角模糊置信度计算装置及方法
CN110018362A (zh) * 2019-04-16 2019-07-16 中国电子科技集团公司第二十九研究所 一种主波束对称的宽带天线的相位中心测量方法
CN110018362B (zh) * 2019-04-16 2021-05-04 中国电子科技集团公司第二十九研究所 一种主波束对称的宽带天线的相位中心测量方法
CN112068096A (zh) * 2020-06-09 2020-12-11 周口师范学院 一种被动雷达复合扫描三维空间测向方法及系统
CN112068096B (zh) * 2020-06-09 2024-04-19 周口师范学院 一种被动雷达复合扫描三维空间测向方法及系统
CN113093167A (zh) * 2021-04-19 2021-07-09 成都安智杰科技有限公司 一种目标测向方法及装置
CN113203980A (zh) * 2021-04-20 2021-08-03 北京通广龙电子科技有限公司 高精度快速无线电测向方法及系统
CN117289037A (zh) * 2023-11-23 2023-12-26 南京华成微波技术有限公司 一种大功率相控阵天线平面近场测试方法及系统
CN117289037B (zh) * 2023-11-23 2024-02-09 南京华成微波技术有限公司 一种大功率相控阵天线平面近场测试方法及系统

Also Published As

Publication number Publication date
CN109116297B (zh) 2022-09-27

Similar Documents

Publication Publication Date Title
CN109116297A (zh) 一种被动雷达空间谱估计与合成波束的联合测向方法
EP3254133B1 (en) Direction finding using signal power
CN111458711A (zh) 星载双波段sar系统和舰船目标的探测方法
EP2437079A1 (en) Radar system
KR101807499B1 (ko) 상호간섭효과를 보상한 원형 배열 안테나 기반의 도래각 추정 장치 및 추정 방법
CN111649803B (zh) 基于垂直线性阵列的三维雷达物位计及其设计方法
CN107526077B (zh) 一种便携式高频地波雷达浪高反演方法
CN104777467B (zh) 基于频率扫描天线的目标检测方法
CN109597076B (zh) 用于地基合成孔径雷达的数据处理方法及装置
CN106329152A (zh) 一种针对半球覆盖波束成形的阵列设计方法
CN108089161B (zh) 一种基于功率反馈的天线阵列合成波束空域选择方法
RU2546330C1 (ru) Способ поляризационно-чувствительного радиоконтроля подвижных объектов
Sun et al. Four-dimensional high-resolution automotive radar imaging exploiting joint sparse-frequency and sparse-array design
CN105807252B (zh) 基于矩形天线阵列的波达角估计方法
CN113960523A (zh) 基于fpga的通用超宽带校正测向方法及系统
RU2546329C1 (ru) Способ поляризационно-чувствительного обнаружения подвижных объектов
CN109828232A (zh) 一种阵列天线测角模糊置信度计算装置及方法
Chen et al. Multichannel radar forward looking superresolution imaging via atomic norm minimization
CN114222242B (zh) 基于数字相控阵体制的通信定位一体化方法及系统
Hoffmann et al. Rx beamforming for long baseline multistatic radar networks
CN114994676A (zh) 一种一站固定式双站低频超宽带sar运动目标成像方法
RU2557250C1 (ru) Способ скрытной радиолокации подвижных объектов
CN105652234A (zh) 一种循环空间谱测向方法
Yu et al. W-Band Co-Prime Array FMCW MIMO Radar With IPCE Algorithm for High-Angular Resolution
CN108107424B (zh) 一种基于伺服驻留的弹载空域搜索方法

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