CN107180259B - 一种基于系统辨识的stap训练样本选择方法 - Google Patents
一种基于系统辨识的stap训练样本选择方法 Download PDFInfo
- Publication number
- CN107180259B CN107180259B CN201710337764.3A CN201710337764A CN107180259B CN 107180259 B CN107180259 B CN 107180259B CN 201710337764 A CN201710337764 A CN 201710337764A CN 107180259 B CN107180259 B CN 107180259B
- Authority
- CN
- China
- Prior art keywords
- clutter
- cut
- model
- sample
- training
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- 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
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
-
- 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
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/417—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section involving the use of neural networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Artificial Intelligence (AREA)
- Computer Networks & Wireless Communication (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于系统辨识的STAP训练样本选择方法,属于雷达技术领域。本发明首先给出了基于系统辨识的样本选择模型,在此基础上,进一步提出用神经网络辨识CUT的杂波模型,然后利用辨识出的模型对其他距离单元进行滤波,最后根据滤波后输出的方差大小选择出与CUT杂波协方差矩阵相似的样本。本发明解决了STAP中训练样本不足的问题,估计的杂波协方差矩阵更准确,提高了杂波的抑制性能。
Description
技术领域
本发明属于雷达技术领域,具体涉及到基于系统辨识的STAP训练样本选择方法。
背景技术
空时自适应处理(STAP)是雷达和通信信号处理中的一项关键技术,被广泛用于雷达及通信信号处理之中。STAP要求设计最优权向量时,使得输出信杂噪比 (SCNR)最大。其中,s是目标信号空时导向向量,RCUT为待检测距离单元的杂波协方差矩阵。然而实际中RCUT是未知的,需通过选取的训练样本来估计。设训练样本的杂波协方差矩阵是RTS,则选取的训练样本的要求是:
(1)应满足RTS=RCUT。
(2)训练样本要足够多,因为用于估计协方差矩阵的训练样本数至少大于系统自由度的 2倍,才能使STAP的输出信噪比损失小于3dB。
目前训练样本选择方法主要分为三类。第一类是功率选择训练算法,其基本原理是挑选杂波功率大的样本,从而加深杂波凹口深度。第二类是广义内积(GIP)算法,其基本原理是通过广义内积统计量来剔除掉训练样本中与待检测单元杂波背景统计特性不同的样本。第三类是基于波形相似的方法。
其中,基于波形相似的方法尤其受到关注。该类方法通常是在时域或频域上选取与CUT (待检测距离单元)杂波波形相似性较高的样本。文献《Zhang X,Yang Q,DengW.Weak target detection within the nonhomogeneous ionospheric clutterbackground of HFSWR based on STAP[J].International Journal of Antennas andPropagation,2013,2013》中提出根据时域波形相似性选取样本,计算待选样本与CUT杂波的时域波形相关系数,选择相关系数大于所设阈值的样本作为训练样本。但该方法直接舍弃了相似性较低的样本,导致样本的利用率较低。针对该问题,文献《Xin Z,Yanhua S,Qiang Y,et al.Space-time adaptive processing-based algorithm for meteor trailsuppression in high-frequency surface wave radar[J].IET Radar,Sonar&Navigation,2015,9(4):429-436》进一步提出基于相关系数加权的方法估计协方差矩阵。相关系数小的样本被赋予小的权值,从而使得相似性较低的样本也能被用来估计RCUT,从而提高了样本利用率。文献《Yifeng W,Tong W,Jianxin W,et al.Robust training samplesselection algorithm based on spectral similarity for space–time adaptiveprocessing in heterogeneous interference environments[J].IET Radar,Sonar&Navigation,2015,9(7):778-782》和《Wu Y, Wang T,Wu J,et al.Training sampleselection for space-time adaptive processing in heterogeneous environments[J].IEEE Geoscience and Remote Sensing Letters,2015,12(4):691-695》则研究了根据频域波形相似性选择样本,该方法选择与CUT杂波频谱相似的样本作为训练样本,这些方法有效提高了CUT杂波协方差矩阵估计的准确度。
训练样本选择问题的本质是寻找与待检测距离单元杂波具有相同协方差矩阵的样本。而在两个信号的协方差矩阵相同时,这两个信号的波形可能完全不相似,因此选取的训练样本中可能存在波形完全不相似的样本,传统基于波形相似性的样本选择方法容易遗漏大量相似性低的可用样本。
发明内容
本发明针对STAP中基于传统训练样本选择方法的样本漏选问题,提出了一种新型的基于系统辨识的训练样本选择方法。
对于包括N个阵元的雷达天线,用M表示一个相干处理间隔(CPI)内的脉冲数,则第k 个距离单元的回波信号xk可以如下表示:
xk=ξks+ck+nk (1)
其中ck是杂波信号,nk是接收的噪声,s是目标信号空时导向向量,ξk是目标的增益系数, xk,ck,nk,s均为MN维的复向量。用xCUT表示待检测距离单元的回波信号。
STAP的最优权向量w可以通过求解如下优化问题得到:
其中RCUT为待检测距离单元的杂波协方差矩阵,符号“(·)H”表示共轭转置。计算可求得最优权向量为:
然而在实际情况中,由于地形、地貌的空间变化,强散射点以及阵列等原因,导致训练样本的统计特性偏离待检测单元的统计特性,即难以满足。此时,应选择与待检测单元杂波协方差矩阵尽可能相似的样本,并且所选训练样本数量应尽可能多。
现有STAP训练样本的选择方式是挑选时域或频域波形与CUT相似的样本。然而协方差矩阵相同不等价于波形相似,即,
具有相同杂波协方差矩阵的两个距离单元,其波形不一定相似。
对上述的结论证明如下:
设时间序列分别为x=(x1,x2,…,xl,…,xn)T,y=(y1,y2,…,yl,…,yn)T的2个波形满足 xxH=yyH,其中,而中的e表示自然底数,i表示虚数单位,θl,βl表示对应的相位幅角。由于xxH=yyH,则它们的杂波协方差矩阵相同。为简便起见,此处只写出xxH的表达式:
由于xxH=yyH,x与y的模与辐角应分别满足式(6)与(7):
|y1|2=|x1|2,|y2|2=|x2|2,…|yn|2=|xn|2 (6)
β1-β2=θ1-θ2,β1-β3=θ1-θ3,…,β1-βn=θ1-θn (7)
记α=β1-θ1=β2-θ2=…=βn-θn,则x与y的相关系数γx,y为:
显然只有当α=0时,x与y才线性相关;当α=π/2时,则x与y完全不相关。因此,具有相同杂波协方差矩阵的两个距离单元,其波形不一定相似。
因此,在现有STAP训练样本的选择方式会漏选波形不相似、杂波协方差矩阵却相似的样本。其训练样本选择的范围变小了,一些性能更好的样本可能被遗漏了。在地貌、地形变化很快的环境下,此类方法可获得的训练样本将更少,协方差矩阵估计的误差将更大,导致 STAP的性能严重下降。
如果直接选取杂波协方差矩阵与CUT相同的训练样本,将可以避免上述基于波形相似性产生的问题。
本发明中,通过滤波输出方差来判断两个向量的协方差矩阵是否相同。其原理分析如下:
设有两时间序列{x(n)}与{y(n)},记作向量x=(x(1),x(2),…,x(n))T, y=(y(1),y(2),…,y(n))T,Rxx与Ryy分别是各自的协方差矩阵。设时间序列{x(n)}的白化滤波器为H,这里的H可以看作是{x(n)}的模型。滤波输出为白噪声{ex(n)},其方差为用H对{y(n)}滤波,滤波输出是方差为的白噪声{ey(n)}。
因上述模型H可能是线性,也可能是非线性系统,因此下面分两种情况论证。
(A)当模型H是线性系统时。
一个线性时间序列通常可用自回归滑动平均模型(ARMA)建模,{x(n)}的ARMA模型为:
则上述向量x的自相关函数rx(m)为:
(B)当模型H是非线性系统时。
在实际处理中,信号可能是一个复杂的非线性时间序列。双线性模型是一种一般性的非线性时间序列模型,理论上可以证明,任何连续的因果泛函都可以用一个双线性系统来逼近。此时非线性时间序列{x(n)}可采用双线性模型对其建模。Grander和Anderson定义的阶数为 p,q,r,s的双线性时间序列模型满足式(12)所示差分方程:
上式中,ai,bj,ck,l为系统的参数,{ex(n)}是方差为σx 2的白噪声序列。显然双线性模型是线性ARMA模型的直接推广。
上述向量x的自相关函数rx(m)典型表达式为:
其中D为积分域[-π,π],α(ω)、β(ω)、γ0(ω)及Z的表达式如式(14)(15)所示,
多项式α(ω)、β(ω)、γ0(ω)只与该双线性系统的模型参数有关,对于确定的模型,Z只与有关。因此,式(13)中自相关函数rx(m)只取决于白噪声的方差于是,当模型H是非线性系统时,如果那么Rxx=Ryy。
在选择STAP训练样本时,如果把CUT杂波视作上述时间序列{x(n)},待选样本xk视作上述时间序列{y(n)},H为CUT杂波模型。当其对应的滤波输出方差相等时,则可以认为其协方差矩阵相同,则样本可以选作训练样本。于是,样本选择问题转变成了CUT杂波的系统模型H的辨识问题。
因H既可能是线性系统,也可能是复杂的非线性系统。而神经网络不仅能对线性系统进行建模,而且还具有很强的非线性映射能力。基于该考虑,本发明采用神经网络对CUT杂波进行建模,辨识出其模型H。此外,考虑到回声状态网络(ESN)是一种新型递归网络,相较于传统的神经网络具有更好的稳定性及精度,因此优选回声状态神经网络(如雷达科学与技术2015年8月第4期公开的“基于回声状态网络的OTHR海杂波抑制方法”中所提及的回声状态神经网络)辨识CUT的杂波模型H,即基于CUT进行多次杂波预测模型(如回声状态神经网络、或者其他神经网络(如径向基神经网络)、Volterra模型等)训练,通过训练不断调整自身的参数,以近似逼近杂波模型H,取滤波输出方差最小的杂波预测模型作为CUT杂波模型H。辨识出CUT杂波模型H后,使用H对待选样本进行滤波最后根据滤波输出的方差大小选择训练样本。即本发明的基于系统辨识的STAP训练样本选择方法,具体包括下列步骤:
由于杂波具有混沌特性,理论上CUT单元的杂波经模型H后的滤波输出为0。因此,如果待选样本xk与CUT的杂波协方差矩阵相同,其滤波后的输出方差也为0。然而由于雷达回波信号中的噪声等影响,实际上是一个不为0的很小的值,因此可通过设置合适的阈值来选择那些接近于0的待选样本xk作为训练样本,从而实现选择与CUT具有相同的杂波协方差矩阵的待选样本的目的。
本发明适用于STAP中估计协方差矩阵的训练样本选择问题,即基于本发明的训练样本的选择结果,估计CUT杂波的协方差矩阵RCUT,计算STAP的最优权向量w,得出滤波器的输出SCNRout。从而有效提升现有STAP的性能。
与传统基于波形相似选择训练样本的方法相比,本发明的有益效果为:
(1)能获得更多的有效样本,可有效解决样本的短缺问题;
(2)估计的杂波协方差矩阵更准确,杂波的抑制性能更好。
附图说明
图1为本发明具体实施方式的流程图。
图2为第435距离单元雷达回波信号的频谱。
图3为各距离单元与CUT的相关系数。
图4为各距离单元滤波后的归一化输出方差。
图5(a)和图5(b)为两种方法的杂波抑制结果,其中图5(a)为基于相似性的方法,图5(b)为本发明所提方法。
图6为不同输入SCNR条件下的杂波抑制性能。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合实施方式和附图,对本发明作进一步地详细描述。
将本发明用于STAP中的估计协方差矩阵,具体实现步骤如图1所示:
步骤S3:设置合适的阈值(取值接近0的正数),如果小于或等于该阈值,则当前待选样本xk选作训练样本否则,舍弃当前待选样本xk。其原因是:由于杂波具有混沌特性,理论上CUT单元的杂波经H后的滤波输出为0。因此,如果待选样本xk与 CUT的杂波协方差矩阵相同,其滤波后的输出方差也为0。然而由于雷达回波信号中的噪声等影响,实际上是一个不为0的很小的值。
实施例
以天波雷达工作频率f0=18.3MHz,脉冲重复周期T=12ms,脉冲积累个数M=512,相干积累时间CIT=6.144s。回波数据中,已知待检测的第435距离单元有一个多普勒频率为 -5.859的目标,其频谱如图2所示。图3是各距离单元数据与待检测单元之间的相似性,图4 表示各距离单元经训练好的神经网络滤波后的归一化输出方差。
如果基于相似性来选择训练样本,理论上应尽可能选取相关系数接近1的样本。然而受制于可选样本数量的限制,此处根据参考文献《Zhang X,Yang Q,Deng W.Weak targetdetection within the nonhomogeneous ionospheric clutter background of HFSWRbased on STAP[J]. International Journal of Antennas and Propagation,2013,2013》,选取相关系数大于0.7的样本作为参考样本,总共可获得25个参考样本。
若采用本发明所提参考样本选取办法,选取归一化预测误差方差小于0.165的样本作为参考样本,最终能获得的56个有效的训练样本。表1统计了这56个训练样本在各个相关系数区间的数量。从表1可以看出,与CUT相关系数小于0.1的训练样本数量为16,它们与CUT几乎完全不相关。然而,这些样本与CUT杂波却具有相似的协方差矩阵,因此可以选做训练样本。
表1所选训练样本在各相关系数区间的数量
传统的基于相似性的样本选择方法漏选了与CUT波形不相似、杂波协方差矩阵却相似的样本。本发明所提方法的选择目标是与CUT具有相似杂波协方差矩阵的单元,因此能获得更多的训练样本。
分别采用两种方法对第435距离单元进行杂波抑制,仿真结果如图5(a)、5(b)所示。图 5(a)是利用相似性选择样本的方法进行海杂波抑制后的各频率归一化输出SCNR结果。目标 -5.859Hz多普勒频率处的输出SCNR最高,比第二高峰高出11.202dB。图5(b)是本发明所提方法的仿真结果。目标信号的多普勒频率处有一个更高的突起,杂波的输出归一化SCNR全被抑制在-13.34dB以下。相比于图5(a)中基于相似性选择样本的方法,本发明所提方法对杂波的抑制效果更好,更容易检测出目标信号。
为了比较本文所提方法与基于相似性选择样本的方法在不同输入SCNR条件下的杂波抑制性能,本次仿真实验使用另一雷达回波数据,在第144距离单元模拟添加多普勒频率为 fd=1Hz的目标,输入SCNR的取值从-30dB到-10dB取值。图6为两种方法在不同输入SCNR 条件下的杂波抑制性能。纵坐标为目标频率处输出SCNR与其余频率范围内最大输出SCNR 的差值。该差值表征了目标多普勒频率处的输出SCNR的凸起程度,差值越大,说明杂波抑制效果越好,目标更容易被检测到。经计算,本文所提的方法其输出SCNR差值比使用基于相似性的方法平均高出2.71dB,杂波的抑制性能更好。
通过以上仿真,验证了本发明的两点优势:(1)能获得更多的有效样本,可有效解决样本的短缺问题;(2)估计的杂波协方差矩阵更准确,杂波的抑制性能更好。
Claims (2)
1.一种基于系统辨识的STAP训练样本选择方法,其特征在于,包括下列步骤:
步骤1:采用神经网络对CUT杂波进行建模,辨识CUT杂波模型:
其中,杂波预测模型用于白化滤波,所述杂波预测模型为回声状态神经网络、径向基神经网络或Volterra模型;
所述杂波预测模型包括线性和非线性模式;
若为线性模式,则采用自回归滑动平均模型对CUT杂波的时间序列进行建模;
若为非线性模式,则采用双线性模型对CUT杂波的时间序列进行建模;
步骤3:基于两个信号经过滤波后的方差相等,则对应的协方差矩阵相等的原则,选择与CUT具有相似的杂波协方差矩阵的待选样本作为STAP训练样本:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710337764.3A CN107180259B (zh) | 2017-05-15 | 2017-05-15 | 一种基于系统辨识的stap训练样本选择方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710337764.3A CN107180259B (zh) | 2017-05-15 | 2017-05-15 | 一种基于系统辨识的stap训练样本选择方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107180259A CN107180259A (zh) | 2017-09-19 |
CN107180259B true CN107180259B (zh) | 2021-03-30 |
Family
ID=59832556
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710337764.3A Active CN107180259B (zh) | 2017-05-15 | 2017-05-15 | 一种基于系统辨识的stap训练样本选择方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107180259B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109061598B (zh) * | 2018-08-28 | 2022-10-14 | 电子科技大学 | 一种stap杂波协方差矩阵估计方法 |
CN109061599B (zh) * | 2018-08-28 | 2022-06-03 | 电子科技大学 | 一种基于循环平稳和对称先验知识的stap方法 |
US12044797B2 (en) | 2018-12-31 | 2024-07-23 | Atai Labs Pvt Ltd. | Object classification using machine learning |
CN111654264B (zh) * | 2020-05-27 | 2023-04-28 | 桂林长海发展有限责任公司 | 一种信号数据模拟器生成信号脉冲序列方法及系统 |
CN111796518B (zh) * | 2020-06-09 | 2021-08-20 | 吉林大学 | 磁控形状记忆合金执行器位移控制方法 |
CN112906476B (zh) * | 2021-01-20 | 2022-11-04 | 电子科技大学中山学院 | 一种基于信杂噪比损失的机载雷达训练样本选择方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102879767B (zh) * | 2012-06-28 | 2014-02-19 | 北京理工大学 | 一种用于空时自适应处理的干扰目标检测方法 |
CN104215939B (zh) * | 2014-10-10 | 2017-02-15 | 北京航空航天大学 | 一种融合广义对称结构信息的知识辅助空时自适应处理方法 |
CN106375065A (zh) * | 2016-08-31 | 2017-02-01 | 电子科技大学 | 基于存在频偏的空间调制系统的迭代mmse软检测方法 |
-
2017
- 2017-05-15 CN CN201710337764.3A patent/CN107180259B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107180259A (zh) | 2017-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107180259B (zh) | 一种基于系统辨识的stap训练样本选择方法 | |
Sira et al. | Waveform-agile sensing for tracking | |
CN104155650A (zh) | 一种基于熵权值法点迹质量评估的目标跟踪方法 | |
CN106483515A (zh) | 均匀k分布杂波下的最优自适应检测方法 | |
CN107229040B (zh) | 基于稀疏恢复空时谱估计的高频雷达目标检测方法 | |
Han et al. | Censoring outliers in radar data: an approximate ML approach and its analysis | |
Li et al. | A training samples selection method based on system identification for STAP | |
CN110515052A (zh) | 一种基于时间反演的超宽带频域非等间隔采样目标检测方法 | |
CN108333568B (zh) | 冲击噪声环境下基于Sigmoid变换的宽带回波Doppler和时延估计方法 | |
CN112255608A (zh) | 一种基于正交投影的雷达杂波自适应抑制方法 | |
CN106154241B (zh) | 脉冲噪声环境下的韧性平行因子分析算法 | |
CN108957416B (zh) | 脉冲噪声环境下的线性调频信号参数估计方法 | |
CN110426701A (zh) | 一种基于速度维稀疏表示的捷变频雷达目标速度测量方法 | |
CN104101868B (zh) | 基于干扰子空间重构的雷达多假目标干扰抑制方法 | |
Sun et al. | Waveform fusion in sonar signal processing | |
Jiang et al. | Frequency hopping signal parameter estimation algorithm based on time-frequency point correlation | |
Zhang et al. | Impulsive noise excision using robust smoothing | |
De Maio et al. | KB-GLRT: exploiting knowledge of the clutter ridge in airborne radar | |
Sira et al. | A subspace-based approach to sea clutter suppression for improved target detection | |
Roy et al. | A GLRT detector in partially correlated texture based compound-Gaussian clutter | |
Chen et al. | On probing waveforms and adaptive receivers for active sonar | |
Watts | Modelling of coherent detectors in sea clutter | |
Ashok et al. | Performance comparison of matched filter, wavelet denoising and wiener filter technique in communication receivers | |
Eraslan et al. | Turbo Parametric Spectral Estimation Method of Clutter Profile for Adaptive Radar Detection | |
CN113138372B (zh) | 一种基于改进型m/n检测器的多径环境下雷达目标检测方法 |
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 |