CN113504513A - 一种时域非线性调频信号产生方法 - Google Patents
一种时域非线性调频信号产生方法 Download PDFInfo
- Publication number
- CN113504513A CN113504513A CN202110737898.0A CN202110737898A CN113504513A CN 113504513 A CN113504513 A CN 113504513A CN 202110737898 A CN202110737898 A CN 202110737898A CN 113504513 A CN113504513 A CN 113504513A
- Authority
- CN
- China
- Prior art keywords
- frequency
- sequence
- control word
- address
- frequency control
- 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
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
- 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/28—Details of pulse systems
- G01S7/282—Transmitters
-
- 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种时域非线性调频信号产生方法,包括:确定雷达信号的调频频率带宽范围以及频率步进生成频率控制字查找表;先得到非线性函数的采样值序列;对采样值序列进行折线拟合处理,生成采样值序列的拟合系数序列;对拟合系数序列转换为二进制得到步进系数二进制序列并输出至FPGA;FPGA根据步进系数二进制序列生成读取频率控制字查找表的地址;根据地址在频率控制字查找表中查找出对应的频率控制字至直接数字式频率合成DDS;DDS根据频率控制字生成时域非线性调频信号。本发明实现信号频率的非线性变化,适用于雷达应用场景。
Description
技术领域
本发明涉及雷达信号产生技术,特别涉及一种基于时域非线性的FPGA+DDS非线性调频信号产生方法。
背景技术
非线性调频信号作为雷达发射信号时,在接收端直接进行匹配滤波就可以得到高主副瓣比和较窄主瓣宽度的信号,具有较强的目标识别和抗干扰能力。因此非线性调频信号的产生,对雷达模拟测试有着重要作用。传统的非线性调频方法是通过控制器将信号波形参数从波形存储器中读取出来,再传输给D/A转换得到期望的波形。波形存储库通常需要采用多级存储结构,采用非易失性存储器如flash作为上位机波形数据的临时存储器,而采用读写速率较快的SDRAM作为波形查找表。存储空间需求也比较大,否则产生的信号精度很低。若要产生不同参数的波形信号,需要重新产生波形数据并更新到波形存储库。控制器需要同时实现与上位机进行指令交互,存储器的读写控制逻辑以及D/A的控制逻辑,因此控制模块也比较复杂。同时,这种方法要求极高的参考时钟,在实际的工程中很难满足。
另一种方法是基于DDS的波形产生方法,只需向DDS配置相应的幅度、相位以及频率控制字就能产生对应波形。在采用DDS产生非线性调频时,一种普遍的方法是通过上位机将非线性调频信号的频率控制字实时发送到DDS中。虽然实现了非线性频率信号的输出,但对于快速跳频系统,往往频率切换时间要求达到ns量级。此种方法的非线性频率控制字的传送速度显然不能满足快速跳频系统的需求。
发明内容
本发明所要解决的技术问题是,针对传统非线性调频D/A转换法结构复杂与时钟要求极高以及DDS的直接控频方法频率切换时间精度较小的问题,提供一种能快速跳频的非线性调频信号的产生方法。
本发明为解上述技术问题所采用的技术方案是,一种时域非线性调频信号产生方法,包括以下步骤:
1)非线性函数的步进系数二进制序列生成步骤:
1-1)根据精度需求设置非线性曲线的采样点数,得到非线性函数的采样值序列;
1-2)对采样值序列进行折线拟合处理,生成采样值序列的拟合系数序列;
1-3)对拟合系数序列转换为二进制得到步进系数二进制序列并输出至FPGA;
2)FPGA中频率控制字的非线性读取方法:
2-1)确定雷达信号的调频频率带宽范围以及频率步进生成频率控制字查找表;
2-2)频率控制字按频率由大到小依次存放在FPGA的查找表中,地址大小为0~(设置的频率点数-1),FPGA根据步进系数二进制序列生成读取频率控制字查找表的地址,其中当二进制序列值为0时地址保持不变,当二进制序列值为1时且非线性函数采样点间的频率点数与每个频率点脉冲个数达到上限时,则根据当前地址加1或者减1操作;初始地址为用户输入频率下限对应的地址;
2-3)根据地址在频率控制字查找表中查找出对应的频率控制字至直接数字式频率合成DDS;
3)DDS根据频率控制字生成时域非线性调频信号。
以脉间非线性调频为例,设置信号带宽范围内的固定步进线性递增频率控制字查找表,并设置调频上限和调频上限。根据比例系数计算修正的上限值,利用修正后的频率上限计算带宽,以带宽与固定步进计算出带宽范围内频率点数。根据计算得到的调制带宽内频率点数与非线性函数采样点数,计算得到每个非线性采样阶段的频率点数。根据非线性调频周期、带宽内频率点数以及脉冲重复周期计算得到每个频率保持的脉冲个数。
非线性脉间频率调制的FPGA的控制参数包含脉冲重复周期、脉冲宽度、输入的起始频率对应的频率控制字的起始地址传递、每个非线性函数采样值阶段的频率点数以及每个频率点的脉冲个数。在FPGA控制逻辑中,在脉冲下降沿对频率点脉冲个数进行加1操作,计数值达到计数上限后,控制每个非线性阶段频率点数的加1操作,以控制非线性函数变化规律序列的访问。同时,频点脉冲个数计数到上限值后,判断此时非线性函数值的二进制序列值是否为1来确定是否在下个脉冲需要切换频率控制字。若满足上述条件,频率查找表的地址加1操作,更新频率控制字值,并传输到DDS内核实现对应频率信号的输出。当二进制序列地址值达到上限时,地址将会进行减1操作,实现频率值在峰值之间的震荡。
本发明的有益效果是,只需存储线性频率步进查找表,且只需用户给出如调制周期、脉冲重复周期以及脉冲宽度等简单参数配置,实现了频率控制字查找表为线性的条件下,根据频率值保持时间的非线性特征,能快速收实现线性频率步进下的非线性调频信号的产生,实现信号频率的非线性变化,适用于雷达应用场景。
附图说明
图1为本发明提出的一种基于时域非线性的FPGA+DDS非线性调频信号产生方法示意图;
图2为本发明提出的非线性函数二进制序列生成原理示意图;
图3为本发明提出的非线性调频信号频率变化仿真图;
图4为本发明提出的FPGA频率控制字的非线性读取时序图。
具体实施方式
下面结合附图和具体实施例,对本发明进行详细的描述。
本发明提出一种基于时域非线性的FPGA+DDS非线性调频信号产生方法,方法步骤框图如图1所示,包括非线性函数二进制序列生成方法,FPGA中频率控制字的非线性读取方法。
根据提出的频率上限与频率下限,以及决定频率精度的频率点数,计算得到频率步进。频率带宽下限以频率步进递增,得到原始频率列表,再通过式(1)生成频率控制字查找表。
FTW为频率控制字,round为四舍五入取整函数,fout为输出频率,fsysclk为DDS参考时钟。
在FPGA中,若需要变化输出频率,只需通过控制频率控制字索引变化,从而从查找表输出频率值。根据需求设置非线性函数采样点数,得到函数值采样序列。结合频率变化应该由频率下限开始,选取的非线性函数应相位变化,使得函数值从最小值时开始变化。如图2所示,对标准余弦函数做32点采样得到采样值序列,再对函数采样值序列的数值进行折线拟合,准余弦函数中前半个周期的函数采样值序列的16个采样值分别为1、0.98、0.92、0.83、0.71、0.56、0.38、0.2、0、0.38、0.56、0.71、0.83、0.92、0.98,得到拟合反映非线性函数变化规律的步进系数二进制序列,其中图2中的16个采样值的拟合值分别为-0.2*5、-0.2*5、-0.2*4、-0.2*4、-0.2*3、-0.2*3、-0.2*2、-0.2*1、0、0.2*1、0.2*2、0.2*3、0.2*3、0.2*4、0.2*4、0.2*5,拟合系数分别为5、5、4、4、3、3、2、1、0、1、2、3、3、4、4、5,16个拟合系数按照下一个步进系数变化对应的二进制序列为1,其余情况均为0的原则生成的步进系数二进制序列为0、1、0、1、0、1、1、1、1、1、1、0、1、0、1、0;按照此方法最终生成32位采样点的余弦函数的步进系数二进制序列。
步进系数二进制序列将作为查找表频率控制字索引变化的标志,当二进制序列值为1时输出可以跳变(这个方法就是在设置的频率上下限内以非线性函数规律进行输出,频率上下限内,查找表的每个频率都应该被输出,不然会出现频率的丢失,因此地址只会加1),为零时输出保持不变,即以时序非线性化的方式实现频率变化的非线性化。
若采用初始设置的上限频率来计算频率变化的点数,并且根据非线性函数序列来进行频率变化。此时,频率跳变次数变小且步进不变,变化后的频率上限值将不能达到实际的上限值。因此,需要对上限频率进行修正,实现频率根据二进制序列跳变次数减少时,也能变化到实际的上限频率。修正的频率上限的计算方法如下:
得出修正后上限频率为:
所述fmin表示频率下限值,fmax表示频率上限值,Δf表示频率步进,fm′ax表示修正后的频率上限值,p表示非线性函数二进制序列中1的个数,q表示非线性函数的采样点数。
根据修正后的上限频率计算出调频周期内的频率点数N为:
那么每个余弦序列阶段所对应的频率点数fre_point_num为:
最后由用户输入的调制周期与重频周期计算出每个频率点的脉冲个数fre_pulse_num为:
其中Tcos为余弦函数周期,T为重频周期。
利用MATLAB对算法仿真,为了能更好的与非线性函数做对比,设置上限频率为1,下限频率为0,频率点数设置了一个任意值为147。根据上述计算参数与二进制余弦序列拟合的按32采样点余弦函数规律变化的非线性频率输出如图3所示。每个余弦函数采样点之间的间隔为实际输出的频率点数,可以看到频率变化规律与标准非线性函数非常逼近。
在FPGA中的设计控制逻辑中,当fre_pulse_num计满上限值是将控制fre_point_num进行加1操作。当fre_point_num计数到达其上限值时,将会控制非线性函数的二进制序列索引(二进制序列索引是指二进制序列的元素索引地址,因为需要按序读取二进制序列值来判断序列值是否为1)的加1操作。当访问到的二进制序列值为1且fre_pulse_num计数满上限值时,对查找表地址进行加1,即输出频率发生变化;否则当前访问查找表地址不变,输出频率保持不变。控制时序如图4所示,假设根据当前设置的上下限频率,计算得到的每个频率点的脉冲点数fre_pulse_num_cnt与每个正弦阶段包含的频率点数fre_point_num_cnt分别为4和2,当前顺序读取到二进制序列索引为2的序列值为0,频率控制字查找表地址为addr0。当非线性调频信号模块使能EN为高,从图4的第一个osk下降沿开始,每当且DDS外部输出开关信号osk的下降沿与时钟clk上升沿同时有效时,fre_pulse_num_cnt计满清零,fre_point_num_cnt也计满清零。fre_pulse_num_cnt和fre_point_num_cnt计满清零后,fre_pulse_num_cnt从0开始计数,期间fre_point_num_cnt保持不变,当fre_pulse_num_cnt计数值为3(计满4个clk周期)时,fre_pulse_num_cnt清零。在下一个时钟上升沿与osk下降沿重新计数,且fre_point_num_cnt开始加1,直到再次计满4个clk周期,在下一个时钟上升沿与osk下降沿同时来到时,重新开始计数。二进制序列索引index则以一个fre_point_num_cnt计数周期进行加1,直到计满15后减1,循环检测二进制序列值输出周期性按余弦规律变化的频率控制字。二进制序列索引index加1指向下一个二进制序列值。若下一个二进制序列值为1,拉高频率变化指示信号keep_fre_flag。在下一个时钟上升沿,检测到keep_fre_flga为高,频率控制器查找表地址加1,经过一个时钟周期的读取延迟读出此地址对应的频率控制字。频率控制字的索引起始地址是频率下限对应的地址,因此频率控制字变化规律是按照非线性序列变化输出到DDS中,实现了对频率变化的一种非线性函数逼近输出。DDS最终根据频率控制字输出频率按非线性规律变化调频信号。
本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其他各种具体变形和组合,这些变形和组合仍然在本发明的保护范围之内。
Claims (1)
1.一种时域非线性调频信号产生方法,其特征在于,包括以下步骤:
1)非线性函数的步进系数二进制序列生成步骤:
1-1)根据精度需求设置非线性曲线的采样点数,得到非线性函数的采样值序列;
1-2)对采样值序列进行折线拟合处理,生成采样值序列的拟合系数序列;
1-3)对拟合系数序列转换为二进制得到步进系数二进制序列并输出至FPGA;
2)FPGA中频率控制字的非线性读取方法:
2-1)确定雷达信号的调频频率带宽范围以及频率步进生成频率控制字查找表;
2-2)频率控制字按频率由大到小依次存放在FPGA的查找表中,地址大小为0~(设置的频率点数-1),FPGA根据步进系数二进制序列生成读取频率控制字查找表的地址,其中当二进制序列值为0时地址保持不变,当二进制序列值为1时且非线性函数采样点间的频率点数与每个频率点脉冲个数达到上限时,则根据当前地址加1或者减1操作;初始地址为用户输入频率下限对应的地址;
2-3)根据地址在频率控制字查找表中查找出对应的频率控制字至直接数字式频率合成DDS;
3)DDS根据频率控制字生成时域非线性调频信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110737898.0A CN113504513B (zh) | 2021-06-30 | 2021-06-30 | 一种时域非线性调频信号产生方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110737898.0A CN113504513B (zh) | 2021-06-30 | 2021-06-30 | 一种时域非线性调频信号产生方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113504513A true CN113504513A (zh) | 2021-10-15 |
CN113504513B CN113504513B (zh) | 2023-05-23 |
Family
ID=78011451
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110737898.0A Active CN113504513B (zh) | 2021-06-30 | 2021-06-30 | 一种时域非线性调频信号产生方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113504513B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115166640A (zh) * | 2022-06-08 | 2022-10-11 | 中国计量大学 | 一种时域Frank码信号产生方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1645163A (zh) * | 2005-01-28 | 2005-07-27 | 武汉大学 | 线性数字调频信号的产生方法 |
CN102073032A (zh) * | 2010-11-02 | 2011-05-25 | 中国兵器工业第二○六研究所 | 多波形雷达信号的模块化产生方法 |
CN102780490A (zh) * | 2012-08-14 | 2012-11-14 | 武汉滨湖电子有限责任公司 | 一种dds直接形成式的超宽带调频信号产生电路及方法 |
CN103780523A (zh) * | 2012-10-24 | 2014-05-07 | 中兴通讯股份有限公司 | 数字预失真数据的处理的方法及装置 |
WO2014094106A1 (en) * | 2012-12-21 | 2014-06-26 | Raytheon Canada Limited | Methods and apparatus for a radar having windfarm mitigation |
CN104267385A (zh) * | 2014-10-16 | 2015-01-07 | 中国科学院电子学研究所 | 具有预失真功能的线性调频信号发生器 |
CN204832855U (zh) * | 2015-06-30 | 2015-12-02 | 安徽四创电子股份有限公司 | 一种基于ad9914设计的高采样率频率分集脉冲信号发生器 |
CN108594214A (zh) * | 2018-04-17 | 2018-09-28 | 西安电子科技大学 | 基于fpga的参数可调的线性调频信号产生装置及其产生方法 |
WO2019023529A1 (en) * | 2017-07-27 | 2019-01-31 | Blackmore Sensors and Analytics Inc. | METHOD AND SYSTEM FOR USING SQUARE WAVE DIGITAL FLUCTUATION SIGNAL FOR OPTICAL FLUCTUATION RANGE DETECTION |
CN109521992A (zh) * | 2018-11-14 | 2019-03-26 | 桂林电子科技大学 | 一种无乘法器的基于cordic算法的线性调频信号生成方法 |
US20190346556A1 (en) * | 2018-10-31 | 2019-11-14 | Institute Of Electronics, Chinese Academy Of Sciences | Method for range ambiguity suppression based on multi-degree-of-freedom frequency modulation signal |
-
2021
- 2021-06-30 CN CN202110737898.0A patent/CN113504513B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1645163A (zh) * | 2005-01-28 | 2005-07-27 | 武汉大学 | 线性数字调频信号的产生方法 |
CN102073032A (zh) * | 2010-11-02 | 2011-05-25 | 中国兵器工业第二○六研究所 | 多波形雷达信号的模块化产生方法 |
CN102780490A (zh) * | 2012-08-14 | 2012-11-14 | 武汉滨湖电子有限责任公司 | 一种dds直接形成式的超宽带调频信号产生电路及方法 |
CN103780523A (zh) * | 2012-10-24 | 2014-05-07 | 中兴通讯股份有限公司 | 数字预失真数据的处理的方法及装置 |
WO2014094106A1 (en) * | 2012-12-21 | 2014-06-26 | Raytheon Canada Limited | Methods and apparatus for a radar having windfarm mitigation |
CN104267385A (zh) * | 2014-10-16 | 2015-01-07 | 中国科学院电子学研究所 | 具有预失真功能的线性调频信号发生器 |
CN204832855U (zh) * | 2015-06-30 | 2015-12-02 | 安徽四创电子股份有限公司 | 一种基于ad9914设计的高采样率频率分集脉冲信号发生器 |
WO2019023529A1 (en) * | 2017-07-27 | 2019-01-31 | Blackmore Sensors and Analytics Inc. | METHOD AND SYSTEM FOR USING SQUARE WAVE DIGITAL FLUCTUATION SIGNAL FOR OPTICAL FLUCTUATION RANGE DETECTION |
CN110869798A (zh) * | 2017-07-27 | 2020-03-06 | 布莱克莫尔传感器和分析有限责任公司 | 用于使用方波数字线性调频信号进行光学线性调频距离检测的方法和系统 |
CN108594214A (zh) * | 2018-04-17 | 2018-09-28 | 西安电子科技大学 | 基于fpga的参数可调的线性调频信号产生装置及其产生方法 |
US20190346556A1 (en) * | 2018-10-31 | 2019-11-14 | Institute Of Electronics, Chinese Academy Of Sciences | Method for range ambiguity suppression based on multi-degree-of-freedom frequency modulation signal |
CN109521992A (zh) * | 2018-11-14 | 2019-03-26 | 桂林电子科技大学 | 一种无乘法器的基于cordic算法的线性调频信号生成方法 |
Non-Patent Citations (5)
Title |
---|
YU MI: "Design of High Performance Linear Frequency Modulation Source for Terahertz Imaging", 《2018 11TH UK-EUROPE-CHINA WORKSHOP ON MILLIMETER WAVES AND TERAHERTZ TECHNOLOGIES (UCMMT)》 * |
叶忠彬: "基于FPGA的数字频率合成系统的设计与实现", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
姚琳: "基于直接数字频率合成的多模式信号源研究与设计", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
范欢欢: "可重构非线性调频信号的设计与实现", 《信息通信》 * |
薛文: "调频信号的数字合成方法", 《信息技术与信息化》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115166640A (zh) * | 2022-06-08 | 2022-10-11 | 中国计量大学 | 一种时域Frank码信号产生方法 |
CN115166640B (zh) * | 2022-06-08 | 2024-05-07 | 中国计量大学 | 一种时域Frank码信号产生方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113504513B (zh) | 2023-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4153026B2 (ja) | Ad変換器およびad変換方法 | |
CN105718404B (zh) | 一种基于fpga的方波发生器及方法 | |
CN102882517B (zh) | 一种生成低失真度的低频正弦信号的装置及方法 | |
CN109142820B (zh) | 一种脉冲波发生装置 | |
CN108710404B (zh) | 一种混合信号发生器 | |
CN110488228B (zh) | 线性调频信号生成方法、装置及存储介质 | |
US20190196999A1 (en) | Fpga-based square-wave generator and square-wave generation method | |
WO2004019496A1 (ja) | 波形発生方法、波形発生プログラム、波形発生回路、及びレーダ装置 | |
CN109521992B (zh) | 一种无乘法器的基于cordic算法的线性调频信号生成方法 | |
CN113504513A (zh) | 一种时域非线性调频信号产生方法 | |
CN108227828B (zh) | 一种序列信号发生器及序列信号产生方法 | |
US20210004041A1 (en) | Sequence signal generator and sequence signal generation method | |
CN103178843A (zh) | 一种具有扫频功能的信号源及其生成方法 | |
CN109709851B (zh) | 一种复杂调制信号高速生成装置和方法 | |
Sharma et al. | Design and implementation of a re-configurable versatile direct digital synthesis-based pulse generator | |
CN103178784A (zh) | 一种振荡键控调制方法、装置和函数信号发生器 | |
CN113376585B (zh) | 一种高分辨率脉冲信号合成装置 | |
CN107797441B (zh) | 低抖动时钟信号的时间幅度调制数字域表示方法 | |
CN104660218A (zh) | 一种任意波形合成器 | |
Liu et al. | Direct digital frequency synthesizer based on curve approximation | |
Zhang et al. | An FPGA based real-time radar target simulator with high spur suppression | |
CN111641490A (zh) | 一种采样时钟高精度相位校准与时间基准确定方法 | |
CN111756357B (zh) | 周期信号的参数确定方法及电子设备 | |
JP2017163525A (ja) | 増幅回路とアンテナモジュール及び無線通信装置 | |
CN110808935B (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 |