CN114185030A - 一种雷达信号脉内调制识别方法及系统 - Google Patents
一种雷达信号脉内调制识别方法及系统 Download PDFInfo
- Publication number
- CN114185030A CN114185030A CN202111478263.XA CN202111478263A CN114185030A CN 114185030 A CN114185030 A CN 114185030A CN 202111478263 A CN202111478263 A CN 202111478263A CN 114185030 A CN114185030 A CN 114185030A
- Authority
- CN
- China
- Prior art keywords
- signal
- bandwidth
- frequency
- modulation signal
- radar signal
- 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.)
- Pending
Links
Images
Classifications
-
- 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
- G01S13/00—Systems 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
-
- 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/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
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
本发明公开了一种雷达信号脉内调制识别方法及系统,属于电子对抗技术领域,包括:获取雷达信号的3dB带宽;响应于检测到3dB带宽大于预设带宽阈值,将雷达信号输入预建立的调频信号识别模型,输出调频信号类型,所述调频信号类型为脉内调制识别结果;响应于检测到3dB带宽不大于预设带宽阈值,将雷达信号输入预建立的调相信号识别模型,输出调相信号类型,所述调相信号类型为脉内调制识别结果;对常用的LPI雷达信号提供了一种系统性的脉内调制识别方法,提高识别精确度。
Description
技术领域
本发明涉及一种雷达信号脉内调制识别方法及系统,属于电子对抗技术领域。
背景技术
目前,宽频带和复杂波形体制的雷达在军事领域中得到了广泛的应用,传统的基于信号载频、脉宽、到达方向、到达角度、信号幅度五参数识别方法,已很难对其进行有效的识别;为了更加有效的检测雷达信号,对雷达信号进行脉内特征的识别已经成为一项紧迫而严峻的任务,脉内特征分析的主要目的是实现不同脉内调制类型信号的自动区分与识别,并检测出相应的脉内调制参数(起始频率、结束频率、调频斜率、码元个数、码元宽度、码元序列等);通过脉内调制特征分析,不仅有利于更全面地描述雷达信号,提高信号分选与识别的准确率,而且能够了解对方雷达的用途与性能,还可以向我方干扰机提供更准确的敌方雷达参数,从而提高对敌干扰的效果;对于电子战应用来说,要求雷达侦察设备具有实时的脉内分析能力已成为一种必然趋势。
采用脉内频率调制和相位调制的主要信号形式是脉冲压缩信号大时宽带宽信号;常见的LPI雷达信号有线性调频(LFM)、非线性调频(NLFM)、二相编码(BPSK)、四相编码(QPSK)、频率编码(FSK)等;在实际系统中,人们研究不同雷达辐射源信号的脉内调制特征,并提出了多种脉内调制特征的识别方法,但这些方法都只适用于少数二、三种信号、没有充分考虑系统性和实时性的识别问题,如何确保实时性和复杂环境下调制特征识别的可靠性,这是非常困难的任务;因此亟须一种实时性好可靠性高系统性高的调制特征识别方法。
发明内容
本发明的目的在于提供一种雷达信号脉内调制识别方法及系统,对常用的LPI雷达信号提供了一种系统性的脉内调制识别方法,提高识别精确度。
为实现以上目的,本发明是采用下述技术方案实现的:
第一方面,本发明提供了一种雷达信号脉内调制识别方法,包括:
获取雷达信号的3dB带宽;
响应于检测到3dB带宽大于预设带宽阈值,将雷达信号输入预建立的调频信号识别模型,输出调频信号类型,所述调频信号类型为脉内调制识别结果;
响应于检测到3dB带宽不大于预设带宽阈值,将雷达信号输入预建立的调相信号识别模型,输出调相信号类型,所述调相信号类型为脉内调制识别结果。
结合第一方面,进一步的,所述雷达信号的3dB带宽通过以下方法得到:
超过半功率点的最大频点减去最小频点得到的信号带宽为3dB带宽。
结合第一方面,进一步的,将雷达信号输入预建立的调频信号识别模型,输出调频信号类型:
获取雷达信号的有效带宽,将3dB带宽除以有效带宽得到带宽比值;
响应于检测到带宽比值大于预设比值阈值,输出调频信号类型为FSK信号;
响应于检测到带宽比值不大于预设比值阈值,对雷达信号进行短时傅里叶变换,搜索短时傅里叶变换中的最大值得到时频脊线,根据时频脊线计算频率变化率,若频率变化率小于预设变化率阈值,输出调频信号类型为LFM信号,否则输出调频信号类型为NLFM信号。
结合第一方面,进一步的,所述有效带宽为雷达信号一个波束内超过半功率点的信号带宽。
结合第一方面,进一步的,所述短时傅里叶变换计算公式如下:
其中,s(τ)为输入的雷达信号,h*(τ-t)为分析窗函数取共轭。
结合第一方面,进一步的,根据时频脊线计算频率变化率:
对时频脊线进行二次多项式曲线拟合,利用时频脊线线性插值后的起始时刻、中间时刻及末端时刻这三点将整段雷达信号分为两段,起始时刻至中间时刻为前一段,中间时刻至末端时刻为后一段,分别计算前一段和后一段的频率变化率并做差值,得到最终的频率变化率。
结合第一方面,进一步的,将雷达信号输入预建立的调相信号识别模型,输出调相信号类型:
获取雷达信号经傅里叶变换后的谱峰个数,对雷达信号进行离散Haar小波变换得到离散Haar小波变换峰值;
响应于检测到谱峰个数小于等于3且离散Haar小波变换峰值小于预设峰值阈值,输出调相信号类型为点频信号;
响应于检测到谱峰个数小于等于3和离散Haar小波变换峰值小于预设峰值阈值中任一条件不满足,将雷达信号平方后进行傅里叶变换,获取平方后的谱峰个数,若平方后的谱峰个数小于等于2,输出调相信号类型为BPSK信号,否则将雷达信号四次方后进行傅里叶变换,获取四次方后的谱峰个数,若四次方后的谱峰个数小于等于2,输出调相信号类型为QPSK信号,否则输出调相信号类型为其他信号。
结合第一方面,进一步的,所述离散Haar小波变换的计算公式如下:
其中,α是离散Haar小波变换尺度,n是数据长度,s(k)是原始时域数据。
第二方面,本发明还提供了一种雷达信号脉内调制识别系统,包括:
带宽获取模块:用于获取雷达信号的3dB带宽;
调频信号识别模块:用于响应于检测到3dB带宽大于预设带宽阈值,将雷达信号输入预建立的调频信号识别模型,输出调频信号类型,所述调频信号类型为脉内调制识别结果;
调相信号识别模块:用于响应于检测到3dB带宽不大于预设带宽阈值,将雷达信号输入预建立的调相信号识别模型,输出调相信号类型,所述调相信号类型为脉内调制识别结果。
与现有技术相比,本发明所达到的有益效果是:
本发明提供的一种雷达信号脉内调制识别方法及系统,通过对雷达信号的3dB带宽的粗识别,判断雷达信号是调频信号还是调相信号;若雷达信号是调频信号,则将雷达信号输入预建立的调频信号识别模型中进行进一步的调制特征识别,精确判断雷达信号属于哪一种调频信号;若雷达信号是调相信号,则将雷达信号输入预建立的调相信号识别模型中进行进一步的调制特征识别,精确判断雷达信号属于哪一种调相信号;综上所述,本发明提供的一种雷达信号脉内调制识别方法及系统,对常用的LPI雷达信号提供了一种系统性的脉内调制识别方法,提高识别精确度;
本发明提供的一种雷达信号脉内调制识别方法及系统,通过短时傅里叶变换,搜索短时傅里叶变换中的最大值得到时频脊线,根据时频脊线计算频率变化率,根据频率变化率对LFM信号和NLFM信号进行区分识别,提高识别精确度;
本发明提供的一种雷达信号脉内调制识别方法及系统,根据傅里叶变换的谱峰个数和离散Haar小波变换,对常规点频信号和其他调相信号进行可靠性识别,提高识别精确度。
附图说明
图1是本发明实施例提供的一种雷达信号脉内调制识别方法的流程图。
具体实施方式
下面结合附图对本发明作进一步描述,以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
实施例一
如图1所示,本发明实施例提供了一种雷达信号脉内调制识别方法,首先将雷达信号进行傅里叶变换(FFT),得到雷达信号的3dB带宽,雷达信号经FFT处理后信号频谱主要为单根谱线。
对X(k)数据进行M点数据平滑,平滑取的点数为40,根据X(k)数据进行3dB带宽和有效带宽的估计,3dB带宽定义为超过半功率点的最大频点减去最小频点得到的信号带宽,有效带宽定义为雷达信号一个波束内超过半功率点的信号带宽。
对雷达信号进行初次识别(粗识别),将3dB带宽与预设带宽阈值T3dB比较,T3dB为10M;若3dB带宽大于预设带宽阈值,则雷达信号为调频信号,否则雷达信号为调相信号。
若雷达信号经初次识别后结果为调频信号,对雷达信号进行调频信号识别(细识别),即响应于检测到3dB带宽大于预设带宽阈值,将雷达信号输入预建立的调频信号识别模型,输出调频信号类型,所述调频信号类型为脉内调制识别结果。
获取雷达信号的有效带宽,将3dB带宽除以有效带宽得到带宽比值,将所述带宽比值和预设比值阈值Tr比较,Tr为0.3;若所述带宽比值大于预设比值阈值,则输出雷达信号的调频信号类型为FSK信号;若所述带宽比值小于等于预设比值阈值,对雷达信号进行短时傅里叶变换,搜索短时傅里叶变换中的最大值,得到时频脊线,通过对时频脊线的特征提取来进行调频信号的细分,根据时频脊线计算频率变化率,若频率变化率小于预设变化率阈值Tf,Tf为0.15;则雷达信号为LFM信号,否则雷达信号为NLFM信号。
所述短时傅里叶变换计算公式如下:
其中,s(τ)为输入的雷达信号,h*(τ-t)为分析窗函数取共轭。
频率变化率的计算方法为:对时频脊线进行二次多项式曲线拟合,利用时频脊线线性插值后的起始时刻、中间时刻及末端时刻这三点将整段雷达信号分为两段,起始时刻至中间时刻为前一段,中间时刻至末端时刻为后一段,分别计算前一段和后一段的频率变化率并做差值,得到最终的频率变化率。
若雷达信号经初次识别后结果为调相信号,对雷达信号进行调相信号识别(细识别),即响应于检测到3dB带宽不大于预设带宽阈值,将雷达信号输入预建立的调相信号识别模型,输出调相信号类型,所述调相信号类型为脉内调制识别结果。
统计X(k)的谱峰个数m(即雷达信号经傅里叶变换后的谱峰个数),对雷达信号进行离散Haar小波变换得到离散Haar小波变换峰值;若谱峰个数小于等于3且离散Haar小波变换峰值小于预设峰值阈值Th,则输出调相信号类型为点频信号;若谱峰个数小于等于3和离散Haar小波变换峰值小于预设峰值阈值中任一条件不满足,将雷达信号平方后进行傅里叶变换(FFT),获取平方后的谱峰个数m psk,若平方后的谱峰个数小于等于2,则输出调相信号类型为BPSK信号,若平方后的谱峰个数大于2,将雷达信号四次方后进行傅里叶变换,获取四次方后的谱峰个数m psk,若四次方后的谱峰个数小于等于2,则输出调相信号类型为QPSK信号,否则输出调相信号类型为其他信号。
离散Haar小波变换的计算公式如下:
其中,α是离散Haar小波变换尺度,n是数据长度,s(k)是原始时域数据。
Haar小波变换尺度α计算方法为α=fs/Bv,fs为采样率,Bv为有效带宽。
雷达信号平方后,3dB带宽内谱峰只有1到2根谱线,我们将其判为二相编码信号(BPSK);雷达信号四次方处理后,3dB带宽内谱峰只有1到2根谱线,我们将其判为四相编码信号(QPSK)。
本发明实施例提供的一种雷达信号脉内调制识别方法,识别出的几类典型的雷达辐射源信号:常规雷达信号(CW)、线性调频信号(LFM)、非线性调频信号(NLFM)、二相编码信号(BPSK)、四相编码信号(QPSK)、频率编码信号(FSK),其中频率编码信号(FSK)最大支持16个频点的FSK信号。
实施例二
本发明实施例提供了一种雷达信号脉内调制识别系统,包括:
带宽获取模块:用于获取雷达信号的3dB带宽;
调频信号识别模块:用于响应于检测到3dB带宽大于预设带宽阈值,将雷达信号输入预建立的调频信号识别模型,输出调频信号类型,所述调频信号类型为脉内调制识别结果;
调相信号识别模块:用于响应于检测到3dB带宽不大于预设带宽阈值,将雷达信号输入预建立的调相信号识别模型,输出调相信号类型,所述调相信号类型为脉内调制识别结果。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (10)
1.一种雷达信号脉内调制识别方法,其特征在于,包括:
获取雷达信号的3dB带宽;
响应于检测到3dB带宽大于预设带宽阈值,将雷达信号输入预建立的调频信号识别模型,输出调频信号类型,所述调频信号类型为脉内调制识别结果;
响应于检测到3dB带宽不大于预设带宽阈值,将雷达信号输入预建立的调相信号识别模型,输出调相信号类型,所述调相信号类型为脉内调制识别结果。
2.根据权利要求1所述的一种雷达信号脉内调制识别方法,其特征在于,所述雷达信号的3dB带宽通过以下方法得到:
超过半功率点的最大频点减去最小频点得到的信号带宽为3dB带宽。
3.根据权利要求1所述的一种雷达信号脉内调制识别方法,其特征在于,将雷达信号输入预建立的调频信号识别模型,输出调频信号类型:
获取雷达信号的有效带宽,将3dB带宽除以有效带宽得到带宽比值;
响应于检测到带宽比值大于预设比值阈值,输出调频信号类型为FSK信号;
响应于检测到带宽比值不大于预设比值阈值,对雷达信号进行短时傅里叶变换,搜索短时傅里叶变换中的最大值得到时频脊线,根据时频脊线计算频率变化率,若频率变化率小于预设变化率阈值,输出调频信号类型为LFM信号,否则输出调频信号类型为NLFM信号。
4.根据权利要求3所述的一种雷达信号脉内调制识别方法,其特征在于,所述有效带宽为雷达信号一个波束内超过半功率点的信号带宽。
6.根据权利要求3所述的一种雷达信号脉内调制识别方法,其特征在于,根据时频脊线计算频率变化率:
对时频脊线进行二次多项式曲线拟合,利用时频脊线线性插值后的起始时刻、中间时刻及末端时刻这三点将整段雷达信号分为两段,起始时刻至中间时刻为前一段,中间时刻至末端时刻为后一段,分别计算前一段和后一段的频率变化率并做差值,得到最终的频率变化率。
7.根据权利要求1所述的一种雷达信号脉内调制识别方法,其特征在于,将雷达信号输入预建立的调相信号识别模型,输出调相信号类型:
获取雷达信号经傅里叶变换后的谱峰个数,对雷达信号进行离散Haar小波变换得到离散Haar小波变换峰值;
响应于检测到谱峰个数小于等于3且离散Haar小波变换峰值小于预设峰值阈值,输出调相信号类型为点频信号;
响应于检测到谱峰个数小于等于3和离散Haar小波变换峰值小于预设峰值阈值中任一条件不满足,将雷达信号平方后进行傅里叶变换,获取平方后的谱峰个数,若平方后的谱峰个数小于等于2,输出调相信号类型为BPSK信号,否则将雷达信号四次方后进行傅里叶变换,获取四次方后的谱峰个数,若四次方后的谱峰个数小于等于2,输出调相信号类型为QPSK信号,否则输出调相信号类型为其他信号。
10.一种雷达信号脉内调制识别系统,其特征在于,包括:
带宽获取模块:用于获取雷达信号的3dB带宽;
调频信号识别模块:用于响应于检测到3dB带宽大于预设带宽阈值,将雷达信号输入预建立的调频信号识别模型,输出调频信号类型,所述调频信号类型为脉内调制识别结果;
调相信号识别模块:用于响应于检测到3dB带宽不大于预设带宽阈值,将雷达信号输入预建立的调相信号识别模型,输出调相信号类型,所述调相信号类型为脉内调制识别结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111478263.XA CN114185030A (zh) | 2021-12-06 | 2021-12-06 | 一种雷达信号脉内调制识别方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111478263.XA CN114185030A (zh) | 2021-12-06 | 2021-12-06 | 一种雷达信号脉内调制识别方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114185030A true CN114185030A (zh) | 2022-03-15 |
Family
ID=80603470
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111478263.XA Pending CN114185030A (zh) | 2021-12-06 | 2021-12-06 | 一种雷达信号脉内调制识别方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114185030A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115616490A (zh) * | 2022-10-13 | 2023-01-17 | 南京中新赛克科技有限责任公司 | 一种线性和非线性调频雷达信号的实时识别方法 |
CN117289236A (zh) * | 2023-11-27 | 2023-12-26 | 成都立思方信息技术有限公司 | 短时雷达信号脉内调制类型识别方法、装置、设备及介质 |
-
2021
- 2021-12-06 CN CN202111478263.XA patent/CN114185030A/zh active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115616490A (zh) * | 2022-10-13 | 2023-01-17 | 南京中新赛克科技有限责任公司 | 一种线性和非线性调频雷达信号的实时识别方法 |
CN115616490B (zh) * | 2022-10-13 | 2023-06-09 | 南京中新赛克科技有限责任公司 | 一种线性和非线性调频雷达信号的实时识别方法 |
CN117289236A (zh) * | 2023-11-27 | 2023-12-26 | 成都立思方信息技术有限公司 | 短时雷达信号脉内调制类型识别方法、装置、设备及介质 |
CN117289236B (zh) * | 2023-11-27 | 2024-02-09 | 成都立思方信息技术有限公司 | 短时雷达信号脉内调制类型识别方法、装置、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114185030A (zh) | 一种雷达信号脉内调制识别方法及系统 | |
KR101235059B1 (ko) | 레이더 신호 변조 형태 식별 장치 및 방법 | |
CN110133632B (zh) | 一种基于cwd时频分析的复合调制信号识别方法 | |
CN111693944A (zh) | 雷达有源干扰信号参数提取与干扰样式识别方法及装置 | |
CN111060878B (zh) | 适用于单脉冲的lfm雷达工作模式实时分类方法及装置 | |
CN111505598B (zh) | 一种基于frft域三特征联合检测装置及方法 | |
CN112751629A (zh) | 基于时频图像处理的宽带特定信号检测方法 | |
CN116797796A (zh) | Drfm间歇采样下基于时频分析与深度学习的信号识别方法 | |
Yang et al. | Radar emitter signal recognition based on time-frequency analysis | |
Khan et al. | Robust sparse reconstruction of signals with gapped missing samples from multi-sensor recordings | |
CN106842141B (zh) | 一种高阶多次共轭滞后积脉内分析方法 | |
Berezowski et al. | Identification of fast-changing signals by means of adaptive chaotic transformations | |
CN110632563B (zh) | 一种基于短时傅里叶变换的脉内频率编码信号参数测量方法 | |
CN115795302B (zh) | 一种无线电跳频信号识别方法、系统、终端及介质 | |
CN115219991A (zh) | 一种基于希尔伯特变换的二相编码调制信号识别方法 | |
CN112505649B (zh) | 一种对雷达辐射源信号进行可视化分析的方法和系统 | |
CN111308426B (zh) | 一种适用于单天线接收机的低信噪比周期调频信号检测与分离方法 | |
CN112327260B (zh) | 一种sar回波数据中脉冲式干扰信号的抑制方法和装置 | |
Ren et al. | Specific emitter identification based on intrinsic time-scale-decomposition and image texture feature | |
JP2626579B2 (ja) | 合成開口レーダの画像データ処理方法及び装置 | |
CN109581423B (zh) | 一种运用子空间投影的线性扫频干扰抑制方法及系统 | |
Ahmad et al. | Analysis and Classification of Airborne Radar Signal Types Using Time-Frequency Analysis | |
CN114019236A (zh) | 一种电网谐波单通道混叠目标信号检测方法和装置 | |
GB2588500A (en) | Computer implemented method of decoding a signal | |
EP3712626A1 (en) | High-rate dft-based data manipulator and data manipulation method for high performance and robust signal processing |
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 |