CN109765520A - 一种基于时间调制阵列谐波波束扫描的波达角估计方法 - Google Patents
一种基于时间调制阵列谐波波束扫描的波达角估计方法 Download PDFInfo
- Publication number
- CN109765520A CN109765520A CN201811600163.8A CN201811600163A CN109765520A CN 109765520 A CN109765520 A CN 109765520A CN 201811600163 A CN201811600163 A CN 201811600163A CN 109765520 A CN109765520 A CN 109765520A
- Authority
- CN
- China
- Prior art keywords
- array
- time
- power meter
- output
- switch
- 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
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明属于电子侦察领域,涉及一种基于时间调制阵列的波束扫描波达方向估计方法,配置接收端天线阵列数量为N,共N路射频天线;每个阵元均与单刀单掷开关模块相连,开关的通断由数字信号处理模块控制,调制后的信号经过功率分配器合为一路,合路后的信号经过功率放大器进行放大经过下变频模块,将+2次谐波分量下变频到基带,经过低通滤波器,提取+2次谐波分量;功率计实时检测低通滤波器输出的能量大小,持续记录功率计的输出,经过一次Tm的记录后,根据功率计最大输出当时阵列方向图的最大增益方向,查表得到入射波方向θs。本发明省略了高速数字采样模块,降低了系统对数字信号处理芯片的要求,降低了系统对数字处理模块的要求。
Description
技术领域
本发明属于电子侦察领域,涉及一种基于时间调制阵列的波束扫描波达方向估计方法。
背景技术
波达角(DirectionofArrival,DOA)估计算法研究是阵列信号处理中的关键问题之一,主要研究任务是如何从背景噪声中估计信号的波达方向。常见的测向方法有比幅法、干涉仪测向法、到达时间差测向法、多普勒测向、空间谱测向等等。
幅度比较式测向体制的工作原理是依据电波在行进中,利用测向天线阵或测向天线的方向特性,对不同方向来波接收信号幅度的不同,测定来波方向。干涉仪测向体制的测向原理是依据电波在行进中,从不同方向来的电波到达测向天线阵时,在空间上各测向天线单元接收的相位不同,因而相互间的相位差也不同,通过测定来波相位和相位差,即可确定来波方向。到达时间差(Time Difference of Arrival,TDOA)测向体制的测向原理依据电波在行进中,通过测量电波到达测向天线阵各个测向天线单元时间上的差别,确定电波到来的方向。多普勒测向是依据电波在传播中,遇到与它相对运动的测向天线时,被接收的电波信号产生多普勒效应,测定多普勒效应产生的频移,可以确定来波的方向。比幅法和多普勒法结构简单,但测向精度较差。干涉仪方法具有中等精度和系统复杂性,在实际的测向系统中最常使用。空间谱估计方法和TDOA方法的精度非常好,而它们的硬件和算法复杂度非常高。上世纪五十年代,Bartlett提出了一种波束形成法,这是最早的DOA估计法,但由于阵列物理孔径的限制,该算法对DOA估计的精度很低且抗噪性很差,为了提高算法的估计准确率,研究者们开始了新的研究。
Burg等人提出了最大熵法将DOA估计的研究从线性预测转变到了谱估计领域,开启了现代基于谱估计法研究的时代。Capon等人提出了最小方差法,该方法是基于空间最优滤波的思想,使噪声和来自非目标方向的干扰最小,同时保持目标方向上的信号不变,提高了最小方差法的分辨率,但其运算量很大。Schmidt提出了多重信号分类(Multiple SignalClassification,MUSIC)算法,这一方法的提出开创了空间谱估计研究的新时代,该算法的提出使得DOA估计的研究有了突飞猛进的发展,其利用子空间的特性进行估计的思想也推动了子空间类估计算法的发展。随后科研人员对多重信号分类算法进行了改进和优化。基于子空间分解的思想,Roy和Paulraj等人提出了旋转不变子空间的参数估计法,该算法使用子空间之间旋转不变的特性来估计DOA,为了提高该算法的性能,研究人员提出了多个其衍生版本。
上述DOA估计算法,需要大量射频通道,硬件复杂度较高;此外需要进行复杂的数学运算,占用大量的运算资源。本专利提出了基于时间调制阵列谐波波束扫描的DOA估计方法,避免了大量的数学运算,降低了数字运算模块复杂度及算法的复杂度。
发明内容
本发明提供一种复杂度更低、实现简单、精度较高的DOA估计方法,旨在对接收到的电磁波信号的进行实时高效的达波方向估计。
本发明提出了一种基于时间调制阵列的波束扫描波达方向估计方法,包括以下步骤:
步骤1:配置接收端天线阵列数量为N,共N路射频天线,阵列天线为半波长间距的均匀直线阵,时间调制阵列在接收状态;
步骤2:每个阵元均与单刀单掷开关模块相连,开关的通断由数字信号处理模块控制,通过开关通断,实现对接收信号的调制;假设开关周期为T,则每路开关的导通时间为τ=T/N;
步骤3:调制后的信号经过功率分配器合为一路,合路后的信号经过功率放大器进行放大,功率放大器输出的信号经过下变频模块,将+2次谐波分量下变频到基带,经过低通滤波器,提取+2次谐波分量;
步骤4:经过步骤3后的信号接通到功率计,功率计实时检测低通滤波器输出的能量大小,并将结果传送至数字信号处理模块进行记录;
步骤5:保持各开关导通时间长度τ不变,周期性调整各开关通断的相对关系,周期为Tm,实现阵列+2谐波上方向图的周期性调整,从而实现波束扫描;当方向图的最大增益方向指向入射波方向时,此时功率计输出的能量最大;
步骤6:步骤4中数字信号处理模块,持续记录功率计的输出,经过一次Tm的记录后,根据记录的数据,找到功率计最大输出;
步骤7:根据所记录的功率计最大输出时刻,根据当时阵列方向图的最大增益方向,查表得到入射波方向θs。
本发明提出一种基于时间调制阵列的波束扫描波达方向估计方法,其优点在于:
(1)本发明省略了高速数字采样模块,仅需要记录功率计输出的+2次谐波信号瞬时功率信息,降低了系统对数字信号处理芯片的要求。
(2)本发明不需要复杂的数学运算,对于计算及存储单元的要求很低,进一步降低了系统对数字处理模块的要求。
附图说明
图1为本发明系统流程图;
图2为实施例zΔt=0T时的调制时序;
图3为实施例zΔt=0.7T时的调制时序;
图4为实施例+2次谐波分量方向图;
图5为实施例经过一次扫描周期得到的+2次谐波功率;
图6为实施例z=z'时阵列+2次谐波的方向图。
具体实施方式
下面结合附图对本发明提供的技术方案进行详细说明。
如图1所示,一种基于时间调制阵列的波束扫描波达方向估计方法,包括以下步骤:
步骤1:配置接收端天线阵列数量为N,共N路射频天线,阵列天线为半波长间距的均匀直线阵,时间调制阵列在接收状态。
第n路天线接收到的信号可以表示为:
其中,A0为接收信号的幅度,Fc为载波频率,θ为入射角,λ为波长。
步骤2:每个阵元均与单刀单掷开关模块相连,开关的通断由数字信号处理模块控制,通过开关通断,实现对接收信号的调制,将产生新的频率分量。假设开关周期为T,则每路开关的导通时间为τ=T/N。第n个开关的控制信号Un(t)可以表示为:
其中τn,ON为第n个开关的导通时刻,k为功率计在进行功率估计时需要用到的调制周期T的个数;在扫描的初始,其通断时序如图2和图3所示。为了避免谐波分量混叠,调制信号频率Fm=1/T应当大于待测信号带宽。
步骤3:调制后的信号经过功率分配器合为一路,合路后的信号经过功率放大器进行放大,功率放大器的输出经过下变频模块,将+2次谐波分量下变频到基带。
不同于传统的测向体制将信号的Fc线性搬移至基带,此处需要将功分器输出信号中Fc+2*Fm频率处的信号搬移至基带,并通过低通滤波器滤除其他谐波分量,提取+2次谐波分量。
步骤4:经过步骤3后的信号接通到功率计入口。功率计显示当前输入功率的瞬时值,经过k个调制周期T后,功率计得到了此时信号的+2次谐波分量能量的大小,并将结果传送至数字信号处理模块进行记录。
+2次谐波分量方向图如图4所示。
步骤5:保持各开关导通时间长度τ不变,以每k个调制周期T为单位,调整开关通断时序,此时开关控制信号可表示为:
其中每个天线导通起始时间可表示为0<z≤Zmax。当进行到第Zmax次调整时,即为波束扫描过程的最后一次调制时序的调整,当进行k个周期后,功率计输出Zmax+1个有效数值,此时一次波束扫描结束。经过以上过程,实现阵列+2次谐波上方向图的周期性调整,从而经过Tm=k(Zmax+1)T后,完成一个波束扫描过程。当方向图的最大增益方向指向入射波方向时,此时功率计输出的能量最大。
步骤6:步骤4中数字信号处理模块,持续记录功率计的输出,经过一次Tm的记录后,得到Zmax+1个功率信息,根据记录的数据,找到功率计最大输出所对应的z'。
步骤7:由于z次调整的第二次谐波方向图最大增益方向不断改变,方向图的最大增益方向与z为一一对应关系,在数字信号处理单元中,可以将每z次调整所对应的方向θ0,θ1,θ2,θ3,…,存储在芯片中。根据所记录的功率计最大输出z',查找列表,即可找到对应的方向θs,并输出一次结果。
结果表示:
假设入射信号的DOA在-30度,经过一次扫描周期Tm内,功率计的输出会不断改变,并由数字信号处理模块记录,如图5所示;一次扫描完成后,找到最大值出现的位置z=z',根据该时刻+2次谐波方向图最大值所对应的角度,可以得到目标信号的DOA在-30度,z=z'时的+2次谐波方向图如图6所示,其中方向图最大增益方向为-30度。
Claims (2)
1.一种基于时间调制阵列的波束扫描波达方向估计方法,其特征在于,包括以下步骤:
步骤1:配置接收端天线阵列数量为N,共N路射频天线,阵列天线为半波长间距的均匀直线阵,时间调制阵列在接收状态;
步骤2:每个阵元均与单刀单掷开关模块相连,开关的通断由数字信号处理模块控制,通过开关通断,实现对接收信号的调制;假设开关周期为T,则每路开关的导通时间为τ=T/N;
步骤3:调制后的信号经过功率分配器合为一路,合路后的信号经过功率放大器进行放大,功率放大器输出的信号经过下变频模块,将+2次谐波分量下变频到基带,经过低通滤波器,提取+2次谐波分量;
步骤4:经过步骤3后的信号接通到功率计,功率计实时检测低通滤波器输出的能量大小,并将结果传送至数字信号处理模块进行记录;
步骤5:保持各开关导通时间长度τ不变,周期性调整各开关通断的相对关系,周期为Tm,实现阵列+2谐波上方向图的周期性调整,从而实现波束扫描;
步骤6:步骤4中数字信号处理模块,持续记录功率计的输出,经过一次Tm的记录后,根据记录的数据,找到功率计最大输出;
步骤7:根据所记录的功率计最大输出时刻,根据当时阵列方向图的最大增益方向,查表得到入射波方向θs。
2.根据权利要求1所述的一种基于时间调制阵列的波束扫描波达方向估计方法,其特征在于,所述的步骤5中,当方向图的最大增益方向指向入射波方向时,此时功率计输出的能量最大。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811600163.8A CN109765520B (zh) | 2018-12-26 | 2018-12-26 | 一种基于时间调制阵列谐波波束扫描的波达角估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811600163.8A CN109765520B (zh) | 2018-12-26 | 2018-12-26 | 一种基于时间调制阵列谐波波束扫描的波达角估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109765520A true CN109765520A (zh) | 2019-05-17 |
CN109765520B CN109765520B (zh) | 2023-02-24 |
Family
ID=66450977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811600163.8A Active CN109765520B (zh) | 2018-12-26 | 2018-12-26 | 一种基于时间调制阵列谐波波束扫描的波达角估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109765520B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110221242A (zh) * | 2019-05-20 | 2019-09-10 | 北京航空航天大学 | 一种基于时间调制阵列的无人机侦测方法 |
CN110290457A (zh) * | 2019-05-20 | 2019-09-27 | 北京航空航天大学 | 一种应用于无人机及飞控源探测定位的时间选择阵列系统 |
CN110336627A (zh) * | 2019-07-10 | 2019-10-15 | 电子科技大学 | 基于时间调制的阵列天线幅相调控系统及其实现方法 |
CN111537947A (zh) * | 2020-05-12 | 2020-08-14 | 上海交通大学 | 单射频通道空间谱估计测向系统与方法 |
CN111693932A (zh) * | 2020-06-08 | 2020-09-22 | 西安电子科技大学 | 基于微波光子的多微波信号波达方向角估计方法 |
CN111948641A (zh) * | 2020-08-14 | 2020-11-17 | 上海交通大学 | 基于时间调制技术的调频连续波探测系统及方法 |
CN112014793A (zh) * | 2020-09-15 | 2020-12-01 | 上海交通大学 | 基于方向图可重构天线阵列的二维测向系统及方法 |
CN113572509A (zh) * | 2021-06-23 | 2021-10-29 | 南京理工大学 | 基于时间调制阵列的mimo雷达系统及波束形成方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030025629A1 (en) * | 2001-07-26 | 2003-02-06 | Barrick Donald E. | Multi-station HF FMCW radar frequency sharing with GPS time modulation multiplexing |
CN102246597A (zh) * | 2008-12-12 | 2011-11-16 | 塞瑞斯逻辑公司 | 具有时分环境光反馈响应的基于发光二极管的照明系统 |
CN102856665A (zh) * | 2012-09-06 | 2013-01-02 | 中国铁路通信信号股份有限公司 | 一种频率和波束方向可重构的时间调制阵列天线系统 |
CN102868435A (zh) * | 2012-09-06 | 2013-01-09 | 中国铁路通信信号股份有限公司 | 一种能补偿多普勒效应的自适应阵列天线系统 |
US20140152477A1 (en) * | 2012-12-05 | 2014-06-05 | Crest Semiconductors, Inc. | Time Interleaving Analog-to-Digital Converter |
CN106329153A (zh) * | 2016-08-31 | 2017-01-11 | 电子科技大学 | 一种用于大型异构四维天线阵综合的联合优化方法 |
CN106526532A (zh) * | 2016-10-31 | 2017-03-22 | 电子科技大学 | 一种基于四维天线阵的多普勒测向装置 |
CN106961019A (zh) * | 2016-10-18 | 2017-07-18 | 电子科技大学 | 一种基于多时间调制频率的四维天线阵边带抑制方法 |
CN108196229A (zh) * | 2018-02-11 | 2018-06-22 | 哈尔滨工业大学 | 一种基于时间调制阵列的通信雷达一体化设计方法 |
CN108414966A (zh) * | 2018-01-09 | 2018-08-17 | 上海交通大学 | 一种基于时间调制的宽带线性调频信号测向系统及方法 |
US20180310272A1 (en) * | 2017-03-16 | 2018-10-25 | Lonprox Corporation | Systems and methods for indoor positioning using wireless positioning nodes |
-
2018
- 2018-12-26 CN CN201811600163.8A patent/CN109765520B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030025629A1 (en) * | 2001-07-26 | 2003-02-06 | Barrick Donald E. | Multi-station HF FMCW radar frequency sharing with GPS time modulation multiplexing |
CN102246597A (zh) * | 2008-12-12 | 2011-11-16 | 塞瑞斯逻辑公司 | 具有时分环境光反馈响应的基于发光二极管的照明系统 |
CN102856665A (zh) * | 2012-09-06 | 2013-01-02 | 中国铁路通信信号股份有限公司 | 一种频率和波束方向可重构的时间调制阵列天线系统 |
CN102868435A (zh) * | 2012-09-06 | 2013-01-09 | 中国铁路通信信号股份有限公司 | 一种能补偿多普勒效应的自适应阵列天线系统 |
US20140152477A1 (en) * | 2012-12-05 | 2014-06-05 | Crest Semiconductors, Inc. | Time Interleaving Analog-to-Digital Converter |
CN106329153A (zh) * | 2016-08-31 | 2017-01-11 | 电子科技大学 | 一种用于大型异构四维天线阵综合的联合优化方法 |
CN106961019A (zh) * | 2016-10-18 | 2017-07-18 | 电子科技大学 | 一种基于多时间调制频率的四维天线阵边带抑制方法 |
CN106526532A (zh) * | 2016-10-31 | 2017-03-22 | 电子科技大学 | 一种基于四维天线阵的多普勒测向装置 |
US20180310272A1 (en) * | 2017-03-16 | 2018-10-25 | Lonprox Corporation | Systems and methods for indoor positioning using wireless positioning nodes |
CN108414966A (zh) * | 2018-01-09 | 2018-08-17 | 上海交通大学 | 一种基于时间调制的宽带线性调频信号测向系统及方法 |
CN108196229A (zh) * | 2018-02-11 | 2018-06-22 | 哈尔滨工业大学 | 一种基于时间调制阵列的通信雷达一体化设计方法 |
Non-Patent Citations (6)
Title |
---|
CHONG HE: "Direction Finding by Time-Modulated Array With Harmonic Characteristic Analysis", 《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》 * |
QUN CHEN: "A Novel SDMA Uplink Method Based on Time-Modulated Array for UAV Communications", 《20185 IEEE INTERNATIONAL CONFERENCE ON COMMUNICATION SYSTEMS(ICCS)》 * |
曹岸杰等: "一种用于双通道跳频技术的时间调制阵列系统", 《电波科学学报》 * |
杜海明等: "基于时间调制阵列的单通道测向估计方法", 《探测与控制学报》 * |
贺冲: "时间调制阵列理论与应用研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
高飞等: "基于谐波分析的两阵元时间调制阵列测向方法", 《电子科技》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110221242B (zh) * | 2019-05-20 | 2021-07-02 | 北京航空航天大学 | 一种基于时间调制阵列的无人机侦测方法 |
CN110290457A (zh) * | 2019-05-20 | 2019-09-27 | 北京航空航天大学 | 一种应用于无人机及飞控源探测定位的时间选择阵列系统 |
CN110221242A (zh) * | 2019-05-20 | 2019-09-10 | 北京航空航天大学 | 一种基于时间调制阵列的无人机侦测方法 |
CN110336627A (zh) * | 2019-07-10 | 2019-10-15 | 电子科技大学 | 基于时间调制的阵列天线幅相调控系统及其实现方法 |
CN110336627B (zh) * | 2019-07-10 | 2020-09-11 | 电子科技大学 | 基于时间调制的阵列天线幅相调控系统及其实现方法 |
CN111537947A (zh) * | 2020-05-12 | 2020-08-14 | 上海交通大学 | 单射频通道空间谱估计测向系统与方法 |
CN111693932A (zh) * | 2020-06-08 | 2020-09-22 | 西安电子科技大学 | 基于微波光子的多微波信号波达方向角估计方法 |
CN111948641A (zh) * | 2020-08-14 | 2020-11-17 | 上海交通大学 | 基于时间调制技术的调频连续波探测系统及方法 |
CN111948641B (zh) * | 2020-08-14 | 2022-08-26 | 上海交通大学 | 基于时间调制技术的调频连续波探测系统及方法 |
CN112014793A (zh) * | 2020-09-15 | 2020-12-01 | 上海交通大学 | 基于方向图可重构天线阵列的二维测向系统及方法 |
CN112014793B (zh) * | 2020-09-15 | 2023-11-24 | 上海交通大学 | 基于方向图可重构天线阵列的二维测向系统及方法 |
CN113572509A (zh) * | 2021-06-23 | 2021-10-29 | 南京理工大学 | 基于时间调制阵列的mimo雷达系统及波束形成方法 |
CN113572509B (zh) * | 2021-06-23 | 2024-04-05 | 南京理工大学 | 基于时间调制阵列的mimo雷达系统及波束形成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109765520B (zh) | 2023-02-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109765520A (zh) | 一种基于时间调制阵列谐波波束扫描的波达角估计方法 | |
CN108414966B (zh) | 一种基于时间调制的宽带线性调频信号测向系统及方法 | |
CN110221242B (zh) | 一种基于时间调制阵列的无人机侦测方法 | |
CN110515033A (zh) | 一种基于Toeplitz矩阵恢复的欠通道测向系统及方法 | |
Bialer et al. | Super resolution wide aperture automotive radar | |
CN106950528A (zh) | 一种基于线性调频信号的波达方向估计方法 | |
CN103913645A (zh) | 光纤传感器阵列及天线方向图测量装置、测量方法 | |
Zhou et al. | Accurate DOA estimation with adjacent angle power difference for indoor localization | |
CN110290457A (zh) | 一种应用于无人机及飞控源探测定位的时间选择阵列系统 | |
Wang et al. | Design and implementation of a FPGA and DSP based MIMO radar imaging system | |
CN107144811B (zh) | 一种单通道接收阵列信号的循环子空间测向方法 | |
CN115856767A (zh) | 一种可重构智能超表面辅助的波到达方向估计方法 | |
CN109799475A (zh) | 一种基于时间调制阵列谐波能量检测的无线电测向方法 | |
Lin et al. | A super-resolution-assisted fingerprinting method based on channel impulse response measurement for indoor positioning | |
CN111537947A (zh) | 单射频通道空间谱估计测向系统与方法 | |
CN110531309A (zh) | 存在幅相误差时基于原子范数的相关信号到达角估计方法 | |
RU2517365C2 (ru) | Способ обнаружения и пеленгации источников радиоизлучения на одной частоте | |
CN112986976B (zh) | 一种水流测速雷达 | |
CN114563756A (zh) | 双通道相关干涉仪测向样本线性插值方法 | |
Kim et al. | Cascade AOA Estimation Using Uniform Rectangular Array Antenna | |
Gu et al. | High speed moving target detection using distributed OFDM-MIMO phased radar | |
Chang et al. | Subspace based joint delay and direction of arrival estimation for GNSS multipath signals | |
RU2292650C1 (ru) | Способ обнаружения радиоэлектронных средств | |
Ghobadzadeh et al. | Low-complexity 2D root-MUSIC pairing for an L-shaped array | |
Xin et al. | The realization of digital beamforming based on FPGA and DSP |
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 |