CN117784040A - 一种时频域抗自卫式复合干扰方法 - Google Patents
一种时频域抗自卫式复合干扰方法 Download PDFInfo
- Publication number
- CN117784040A CN117784040A CN202410202388.7A CN202410202388A CN117784040A CN 117784040 A CN117784040 A CN 117784040A CN 202410202388 A CN202410202388 A CN 202410202388A CN 117784040 A CN117784040 A CN 117784040A
- Authority
- CN
- China
- Prior art keywords
- interference
- data
- time
- frequency
- target
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000002131 composite material Substances 0.000 title claims abstract description 12
- 230000006835 compression Effects 0.000 claims abstract description 11
- 238000007906 compression Methods 0.000 claims abstract description 11
- 238000001228 spectrum Methods 0.000 claims description 21
- 239000006185 dispersion Substances 0.000 claims description 16
- 238000001914 filtration Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 7
- 230000000903 blocking effect Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 4
- 230000003595 spectral effect Effects 0.000 claims description 4
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 230000001629 suppression Effects 0.000 claims description 3
- 241000209094 Oryza Species 0.000 description 4
- 235000007164 Oryza sativa Nutrition 0.000 description 4
- 235000009566 rice Nutrition 0.000 description 4
- 210000000582 semen Anatomy 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种时频域抗自卫式复合干扰方法,利用雷达的和通道提取受干扰的脉冲数据,再利用雷达的发射信号构造一个时频域二维的滤波器;利用对多普勒频率的搜索确定目标的多普勒频率,再重构出抑制干扰后的目标信号,进而通过脉冲压缩实现目标的检测;最后得到目标的距离信息;由于发射信号的调频率是已知的,不仅可以避免了估计干扰调频率的复杂过程,从而可以采用先期估计滤波器,并进行存贮。
Description
技术领域
本发明涉及各类发射线性调频信号技术领域,尤其涉及一种时频域抗自卫式复合干扰方法。
背景技术
由于近年来,由于电子对抗技术飞速发展,对抗雷达的各处干扰信号被开发出来,并得到广泛的运用。其中切片干扰和频谱弥散干扰是最常用的自卫式干扰,它们也是目前最有效的对抗线性调频信号的干扰。这类干扰充分利用了线性调频信号的时频特性,产生大量的虚假目标从而压制真实目标,对雷达抗干扰带来了很大的难度。当前,雷达最有效的对抗干扰的方法就是频率捷变,但由于雷达工作需要毫秒级的驻留时间,而切片转发式干扰可以在微秒量级跟上雷达的频率,从而导致雷达的这一对抗措施失效。
目前,对抗切片干扰最常用的是时域匿影方法,但该方法只适用于只存在切片转发式干扰的场景,且目标信号的能量会有严重损失,所以其使用条件受限。另外,对抗频谱弥散干扰则采用频域旋转的方法,再通过低通过滤波或带通滤波的方法来抑制干扰,这种方法需要首先估计出调制的干扰斜率。但在自卫式切片和频谱弥散混合干扰的场景下,由于切片干扰通常和噪声进行卷积处理,导致无法估计出干扰的调制斜率,同时由于频谱弥散干扰和切片干扰不同步,导致无法从时域找出切片干扰,这就导致现有的只能抗切片干扰或者抗频谱弥散干扰的方法无法适用于二种干扰均存在的混合干扰场景。
发明内容
本发明的目的在于针对上述现有技术的不足,提供了一种时频域抗自卫式复合干扰方法,通过构造时频域二维滤波器、通过滤波、搜索估计再重构的方式实现主瓣自卫式切片和频谱弥散混合干扰的抑制。
为实现上述目的,本发明采用了如下技术方案:
本发明提供了一种时频域抗自卫式复合干扰方法,包括以下步骤:
S1、利用雷达的和通道提取受干扰的脉冲数据,再利用雷达的发射信号构造一个时频域二维的滤波器;
S2、利用对多普勒频率的搜索确定目标的多普勒频率,再重构出抑制干扰后的目标信号,进而通过脉冲压缩实现目标的检测;
S3、最后得到目标的距离信息;
所述S1具体为:
S101、接收雷达的和通道的数据进行分析,从数据中提取有干扰的脉冲数据段:
(1);
其中, 为目标信号的回波幅度;为雷达的发射信号;为虚数;为数
据采样时间;为目标的多普勒频率;为自卫式干扰发射的切片干扰;为自
卫式干扰发射的频谱弥散干扰;为噪声;
雷达的发射信号的信号波形为:
(2);
雷达的发射信号能平分成份:
(3);
自卫式干扰发射的切片干扰的波形为:
(4);
自卫式干扰发射的频谱弥散干扰的波形为:
(5);
其中,为发射信号的起始频率;为线性调频信号的斜率;为发射脉冲的
周期;分别为二种干扰的幅度;为雷达发射信号的数据分块;为扰调制
噪声由干扰机按时序设置;为数据分块的采样时间;为频谱弥散信号的数据分块;为干扰机将干扰复制了份;为频谱弥散干扰的调频斜率;符号表示卷积运算;
S102、对雷达的发射信号和自卫式干扰发射的频谱弥散干扰进行短
时傅立叶变换:
;
(6);
对二维矩阵数据取绝对值,并进行归一化处理:
(7);
对二维矩阵数据进行归一化处理:
(8);
其中,、和均为二维矩阵,维数为,为时
间维采样的长度,为频率维采样的长度,为已知信息;为门限;为短时
傅立叶变换;为所得到的滤波器;
所述S3具体为:对所述时域数据进行检测,最大值对应的距离门就是目标所在的距离门,将目标的距离信息输出。
进一步,所述S2具体为:
S201、根据已知目标的跟踪信息确定多普勒搜索范围和搜索步长,搜索范围为:
(9);
搜索步长为:
(10);
其中,多普勒估计值为;为航迹的径向角,为目标的速度,为工作波长,为多普勒估计误差;其中,、、和均为已知信息;为搜
索的次数;
对提取的数据进行次的时频域滤波:
(11);
其中,为点乘,为短时傅立叶变换的反变换,,;
S202、对组时频域滤波的数据取绝对值求和,得到个数据,通过
搜索找出其中的最大值,并得到最大值对应的多普勒频率,则重构的数据为:
(12);
S203、对重构的数据进行特异点剔除,即:
(13);
其中,为门限;
S204、对滤波后的数据进行脉冲压缩,得到时域数据为:
(14);
其中,为脉冲压缩。
进一步,所述多普勒估计误差能够按整数倍来设置。
进一步,所述最大值对应的多普勒频率能直接对最大值周围数据进行幅度
拟合得到。
本发明的有益效果为:(1)由于发射信号的调频率是已知的,不仅可以避免了估计干扰调频率的复杂过程,从而可以采用先期估计滤波器,并进行存贮;
(2)由于发射信号已知,其时频域滤波可以构造得很理想,则很容易抑制大部分的带外干扰,而不需要通过设计专门的滤波器,实现简单方便;
(3)通过带外的滤波之后只剩余少量的带内干扰,此时带内干扰在时域上占的区域比例少且幅度远大于信号,所以通过时域方法可以很容易滤除剩余的频谱弥散干扰;
(4)该方法可用于改造现有雷达的信号处理系统,不需要额外增加处理通道和设备,只需对数字化雷达的进行算法升级即可。
附图说明
图1 为本发明一种时频域抗自卫式复合干扰方法的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,下面结合附图,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
假设雷达发射线性调频信号,干扰机发射切片转发和频谱弥散混合式干扰,实施例中信号的归一化的调频率,干扰的归一化的调频率/>,/>Hz,,/>,/>Hz,/>米/秒,/>米,/>度,/>。下面结合实施例说明一下整个发明的详细步骤:
请参阅图1,一种时频域抗自卫式复合干扰方法,包括以下步骤:
S1、利用雷达的和通道提取受干扰的脉冲数据,再利用雷达的发射信号构造一个时频域二维的滤波器;
S2、利用对多普勒频率的搜索确定目标的多普勒频率,再重构出抑制干扰后的目标信号,进而通过脉冲压缩实现目标的检测;
S3、最后得到目标的距离信息。
所述S1具体为:
S101、接收雷达的和通道的数据进行分析,从数据中提取有干扰的脉冲数据段:
(1);
其中,为目标信号的回波幅度;为雷达的发射信号;为虚数;为数据
采样时间;为目标的多普勒频率;为自卫式干扰发射的切片干扰;为自卫
式干扰发射的频谱弥散干扰;为噪声;
雷达的发射信号的信号波形为:
(2);
雷达的发射信号能平分成份:
(3);
自卫式干扰发射的切片干扰的波形为:
(4);
自卫式干扰发射的频谱弥散干扰的波形为:
(5);
其中,为发射信号的起始频率;为线性调频信号的斜率;为发射脉冲的
周期;分别为二种干扰的幅度;为雷达发射信号的数据分块;为扰调制
噪声由干扰机按时序设置;为数据分块的采样时间;为频谱弥散信号的数据分块;为干扰机将干扰复制了份;为频谱弥散干扰的调频斜率;符号表示卷积运算
在一种具体的实施例中,,切片干扰/>,频谱弥散干扰/>,说明整个波形被切片干扰分成了50份进行调制,每一份卷积一个随机噪声,再构成切片干扰。频谱弥散干扰/>,分成了15份进行调制,每一份是相同的。将二种干扰波形相加则形成干扰的数据。
S102、对雷达的发射信号和自卫式干扰发射的频谱弥散干扰进行短
时傅立叶变换:
;
(6);
对二维矩阵数据取绝对值,并进行归一化处理:
(7);
对二维矩阵数据进行归一化处理:
(8);
其中,、和均为二维矩阵,维数为,为时
间维采样的长度,为频率维采样的长度,为已知信息;为门限;为短时
傅立叶变换;为所得到的滤波器。
在一种具体的实施例中,假设,,,即大于0.5的被置
为1,其它情况均为0。另外,、、、的维数均为。
所述S2具体为:
S201、根据已知目标的跟踪信息确定多普勒搜索范围和搜索步长,
搜索范围为:
(9);
搜索步长为:
(10);
其中,多普勒估计值为;为航迹的径向角,为目标的速度,为工作波长,为多普勒估计误差;其中,、、和均为已知信息;为搜
索的次数;
对提取的数据进行次的时频域滤波:
(11);
其中,为点乘,为短时傅立叶变换的反变换,,;
在一种具体的实施例中,Hz,/>米/秒,/>米,/>度,则/>Hz,所以/>的搜索范围为/>,假设/>,则/>Hz,说明在280至320之间按1Hz步长进行搜索。
S202、对组时频域滤波的数据取绝对值求和,得到个数据,通过
搜索找出其中的最大值,并得到最大值对应的多普勒频率,则重构的数据为:
(12);
在一种具体的实施例中,共对41组数据进行取绝对值求和,并找出41个数据中最大值,由这个最大值确定多普勒频率,由于多普勒频率为310Hz,则第31个数据的幅度最大,对应的多普勒频率为310Hz,再利用这个值对数据进行重构,公式如下
;
S203、对重构的数据进行特异点剔除,即:
(13);
其中,为门限;
在一种具体的实施方式中,为/>所有数据的平均值的10倍,即大于平均值10倍的值被置零,其于值保留原值。
S204、对滤波后的数据进行脉冲压缩,得到时域数据为:
(14);
其中,为脉冲压缩。
所述S3具体为:对所述时域数据进行检测,最大值对应的距离门就是目标
所在的距离门,将目标的距离信息输出;其中,检测可以单元平均CFAR实现,也可以采用选
大、选小、排序等CFAR算法实现;实施例中采用单元平均CFAR方法。
所述多普勒估计误差能够按整数倍来设置。能够进一步增加多普勒的容差
范围;实施例中,是按1倍设置的,如容错范围需扩大至2倍,则可以设计为2;
所述最大值对应的多普勒频率能直接对最大值周围数据进行幅度拟合得
到,能够进一步增加估计的精度,实施例中,采用幅度最大的方式进行,没有采用拟合的方
式,如果需要提高精度,可以采用拟合的方式。
其中,S204中,脉冲压缩可以采用时域方法,也可以用频域方法。实施例中采用时域方法进行卷积处理得到。
以上所述实施例仅表达了本发明的实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求。
Claims (4)
1.一种时频域抗自卫式复合干扰方法,其特征在于,包括以下步骤:
S1、利用雷达的和通道提取受干扰的脉冲数据,再利用雷达的发射信号构造一个时频域二维的滤波器;
S2、利用对多普勒频率的搜索确定目标的多普勒频率,再重构出抑制干扰后的目标信号,进而通过脉冲压缩实现目标的检测;
S3、最后得到目标的距离信息;
所述S1具体为:
S101、接收雷达的和通道的数据进行分析,从数据中提取有干扰的脉冲数据段:
(1);
其中, 为目标信号的回波幅度;/>为雷达的发射信号;/>为虚数;/>为数据采样时间;/>为目标的多普勒频率;/>为自卫式干扰发射的切片干扰;/>为自卫式干扰发射的频谱弥散干扰;/>为噪声;
雷达的发射信号的信号波形为:
(2);
雷达的发射信号能平分成/>份:
(3);
自卫式干扰发射的切片干扰的波形为:
(4);
自卫式干扰发射的频谱弥散干扰的波形为:
(5);
其中,为发射信号的起始频率;/>为线性调频信号的斜率;/>为发射脉冲的周期;分别为二种干扰的幅度;/>为雷达发射信号的数据分块;/>为扰调制噪声由干扰机按时序设置;/>为数据分块的采样时间;/>为频谱弥散信号的数据分块;/>为干扰机将干扰复制了/>份;/>为频谱弥散干扰的调频斜率;符号/>表示卷积运算;
S102、对雷达的发射信号和自卫式干扰发射的频谱弥散干扰/>进行短时傅立叶变换:
;
(6);
对二维矩阵数据取绝对值,并进行归一化处理:
(7);
对二维矩阵数据进行归一化处理:
(8);
其中,、/>和/>均为二维矩阵,维数为/>,/>为时间维采样的长度,/>为频率维采样的长度,/>为已知信息;/>为门限;/>为短时傅立叶变换;/>为所得到的滤波器;
所述S3具体为:对所述时域数据进行检测,最大值对应的距离门就是目标所在的距离门,将目标的距离信息输出。
2.根据权利要求1所述的一种时频域抗自卫式复合干扰方法,其特征在于,所述S2具体为:
S201、根据已知目标的跟踪信息确定多普勒搜索范围和搜索步长,搜索范围为:
(9);
搜索步长为:
(10);
其中,多普勒估计值为;/>为航迹的径向角,/>为目标的速度,/>为工作波长,/>为多普勒估计误差;其中,/>、/>、/>和/>均为已知信息;/>为搜索的次数;
对提取的数据进行次的时频域滤波:
(11);
其中,为点乘,/>为短时傅立叶变换的反变换,,/>;
S202、对组时频域滤波/>的数据取绝对值求和,得到/>个数据,通过搜索找出其中的最大值,并得到最大值对应的多普勒频率/>,则重构的数据为:
(12);
S203、对重构的数据进行特异点剔除,即:
(13);
其中,为门限;
S204、对滤波后的数据进行脉冲压缩,得到时域数据为:
(14);
其中,为脉冲压缩。
3.根据权利要求2所述的一种时频域抗自卫式复合干扰方法,其特征在于,所述多普勒估计误差能够按整数倍来设置。
4.根据权利要求2所述的一种时频域抗自卫式复合干扰方法,其特征在于,所述最大值对应的多普勒频率能直接对最大值周围数据进行幅度拟合得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410202388.7A CN117784040B (zh) | 2024-02-23 | 2024-02-23 | 一种时频域抗自卫式复合干扰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410202388.7A CN117784040B (zh) | 2024-02-23 | 2024-02-23 | 一种时频域抗自卫式复合干扰方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117784040A true CN117784040A (zh) | 2024-03-29 |
CN117784040B CN117784040B (zh) | 2024-05-10 |
Family
ID=90391219
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410202388.7A Active CN117784040B (zh) | 2024-02-23 | 2024-02-23 | 一种时频域抗自卫式复合干扰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117784040B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5790475A (en) * | 1996-10-28 | 1998-08-04 | Multispec Corporation | Process and apparatus for improved interference suppression in echo-location and imaging systems |
UA89945U (ru) * | 2013-10-14 | 2014-05-12 | Запорізький Національний Технічний Університет | Способ помехозащиты когерентно-импульсной рлс от комбинированных помех |
CN113176542A (zh) * | 2021-03-24 | 2021-07-27 | 中国人民解放军空军预警学院 | 一种时域自适应抗主瓣自卫式干扰方法及系统 |
US20210231789A1 (en) * | 2018-06-28 | 2021-07-29 | Richwave Technology Corp. | Doppler signal processing device and method thereof for interference suppression |
CN113702921A (zh) * | 2021-08-17 | 2021-11-26 | 西安电子科技大学 | 基于时频分析的干扰抑制方法 |
CN115453469A (zh) * | 2022-09-21 | 2022-12-09 | 南京理工大学 | 针对fmcw毫米波雷达切片重构干扰的抑制方法及系统 |
CN116359854A (zh) * | 2023-03-23 | 2023-06-30 | 哈尔滨工程大学 | 一种基于YOLOv5的防空警戒雷达复合干扰参数估计方法 |
CN117008064A (zh) * | 2023-08-07 | 2023-11-07 | 青岛九维华盾科技研究院有限公司 | 一种基于改进型YOLOv5的雷达干扰自适应抑制方法 |
-
2024
- 2024-02-23 CN CN202410202388.7A patent/CN117784040B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5790475A (en) * | 1996-10-28 | 1998-08-04 | Multispec Corporation | Process and apparatus for improved interference suppression in echo-location and imaging systems |
UA89945U (ru) * | 2013-10-14 | 2014-05-12 | Запорізький Національний Технічний Університет | Способ помехозащиты когерентно-импульсной рлс от комбинированных помех |
US20210231789A1 (en) * | 2018-06-28 | 2021-07-29 | Richwave Technology Corp. | Doppler signal processing device and method thereof for interference suppression |
CN113176542A (zh) * | 2021-03-24 | 2021-07-27 | 中国人民解放军空军预警学院 | 一种时域自适应抗主瓣自卫式干扰方法及系统 |
CN113702921A (zh) * | 2021-08-17 | 2021-11-26 | 西安电子科技大学 | 基于时频分析的干扰抑制方法 |
CN115453469A (zh) * | 2022-09-21 | 2022-12-09 | 南京理工大学 | 针对fmcw毫米波雷达切片重构干扰的抑制方法及系统 |
CN116359854A (zh) * | 2023-03-23 | 2023-06-30 | 哈尔滨工程大学 | 一种基于YOLOv5的防空警戒雷达复合干扰参数估计方法 |
CN117008064A (zh) * | 2023-08-07 | 2023-11-07 | 青岛九维华盾科技研究院有限公司 | 一种基于改进型YOLOv5的雷达干扰自适应抑制方法 |
Non-Patent Citations (3)
Title |
---|
WANG XIAOGE 等: "Radar anti-false target jamming method based on phase modulation", SYSTEMS ENGINEERING AND ELECTRONICS, vol. 43, no. 9, 30 September 2021 (2021-09-30), pages 2476 - 2483 * |
周必雷 等: "基于盲分离的空时联合处理抗复合干扰方法", 系统工程与电子技术, no. 11, 30 October 2018 (2018-10-30), pages 6 - 15 * |
邓斌: "某型雷达性能参数测试仿真训练系统设计与实现", 2020中国系统仿真与虚拟现实技术高层论坛论文集, 31 December 2022 (2022-12-31), pages 98 - 102 * |
Also Published As
Publication number | Publication date |
---|---|
CN117784040B (zh) | 2024-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5376939A (en) | Dual-frequency, complementary-sequence pulse radar | |
CN109490851B (zh) | 一种基于脉间伪随机码的抗干扰方法 | |
CN106019243B (zh) | 一种基于三次初相和frft的drfm干扰检测和对抗方法 | |
CN112363121B (zh) | 基于回波预处理和相参积累的c&i干扰抑制方法 | |
CN113376601B (zh) | 基于clean算法的捷变频雷达旁瓣抑制方法 | |
CN110133632B (zh) | 一种基于cwd时频分析的复合调制信号识别方法 | |
CN114578296B (zh) | 基于相位编码信号的间歇采样干扰抑制方法 | |
Gulum et al. | Parameter extraction of FMCW modulated radar signals using Wigner-Hough Transform | |
CN114966572A (zh) | 基于脉内分段lfm波形的间歇采样转发干扰对抗方法 | |
CN113176542B (zh) | 一种时域自适应抗主瓣自卫式干扰方法及系统 | |
CN109061626B (zh) | 一种步进频相参处理检测低信杂比动目标的方法 | |
CN113176541B (zh) | 一种自适应抗频谱弥散干扰方法及系统 | |
CN117784040B (zh) | 一种时频域抗自卫式复合干扰方法 | |
CN116400303A (zh) | 一种脉内联合脉间时域编码抗isrj的波形设计方法 | |
Li et al. | A Radar Waveform Design of MCPC Method for Interrupted Sampling Repeater Jamming Suppression via Fractional Fourier Transform. | |
Li et al. | A study on blanket noise jamming to LFM pulse compression radar | |
Zhang et al. | Enhanced Signal Processing Method for Integrated Detection and Jamming System Considering the Complex Environment | |
CN112285656B (zh) | 一种基于航迹特征的抗时敏目标欺骗干扰方法 | |
Tian et al. | An adaptive scheme for chopping and interleaving jamming suppression based on the mismatched filter | |
Rodriguez et al. | A simple clutter suppression approach for OFDM-based passive radar exploiting reciprocal filter | |
Taboada et al. | Intercept receiver signal processing techniques to detect low probability of intercept radar signals | |
CN118244210B (zh) | 一种间歇采样转发干扰主动对抗方法、装置及系统 | |
Fan et al. | A New CCBFM/P4 RF Stealth Signal Waveform Design | |
Tian et al. | An Improved Frequency Domain Interference Suppressor for DBT Processing of High Order BOC Signals | |
CN116774163A (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 |