CN103163504A - 无线电测高雷达实时回波模拟方法及系统 - Google Patents
无线电测高雷达实时回波模拟方法及系统 Download PDFInfo
- Publication number
- CN103163504A CN103163504A CN2013100952626A CN201310095262A CN103163504A CN 103163504 A CN103163504 A CN 103163504A CN 2013100952626 A CN2013100952626 A CN 2013100952626A CN 201310095262 A CN201310095262 A CN 201310095262A CN 103163504 A CN103163504 A CN 103163504A
- Authority
- CN
- China
- Prior art keywords
- radar
- scattering point
- distance
- real
- expression
- 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
Abstract
本发明提出一种无线电测高雷达实时回波模拟方法及系统。其中,方法包括:通过雷达获取与雷达位于预定范围内的三维地图数据;计算雷达的当前发射脉冲重复周期时雷达的天线坐标和覆盖范围;根据雷达接收的回波信号和覆盖范围获取与覆盖范围对应的三维地图数据,以得到每个散射点到天线的距离;根据距离和对应的三维地图数据计算每个散射点对应的延迟值、相位值和幅度值;根据延迟值、相位值和幅度值得到对应散射点的冲激响应,并根据冲激响应得到响应函数;以及对响应函数和雷达的当前发射脉冲进行处理得到对应的回波信号。根据本发明实施例的方法,通过将雷达和三维地图数据结合进行处理得到回波信号,实现了强实时和低延时,并且提高了灵活性。
Description
技术领域
本发明涉及雷达回波模拟技术领域,特别涉及一种无线电测高雷达实时回波模拟方法及系统。
背景技术
传统的InSAR回波模拟是采用非实时的方法。首先,通过全数字仿真产生InSAR回波数据并对其进行存储,然后采用回放的方式实现对InSAR回波的模拟。相对于回放式回波模拟方式,采用实时方式能够满足制导装置在回路闭环仿真的实时性要求。回放式回波模拟方式,首先不能根据发射波形的变化而实时变化回波波形,而且也不能根据航迹和航姿的变化实时确定波束照射范围。
发明内容
本发明的目的旨在至少解决上述的技术缺陷之一。
为达到上述目的,本发明一方面的实施例提出一种无线电测高雷达实时回波模拟方法,包括以下步骤:通过雷达获取与所述雷达位于预定范围内的三维地图数据;计算所述雷达的当前发射脉冲重复周期时所述雷达的天线坐标和覆盖范围;根据所述雷达接收的回波信号和所述覆盖范围获取与所述覆盖范围对应的三维地图数据,以得到每个散射点到天线的距离;根据所述距离和对应的三维地图数据计算每个散射点对应的延迟值、相位值和幅度值;根据所述延迟值、相位值和幅度值得到对应散射点的冲激响应,并根据所述冲激响应得到响应函数;以及对所述响应函数和所述雷达的当前发射脉冲进行处理得到对应的回波信号。
根据本发明实施例的方法,通过将雷达和三维地图数据结合进行处理得到回波信号,实现了强实时和低延时,并且提高了灵活性。
在本发明的一个实施例中,对所述相应函数和雷达的发射脉冲信号的处理方法采用卷积、滤波和变频。
在本发明的一个实施例中,所述雷达到单个散射点的距离通过如下公式计算, 其中,(xna,nr,yna,nr,zna,nr)表示第(na,nr)个散射点的三维坐标,(x0,y0,z0)表示雷达的坐标,Rna,nr表示第(na,nr)个散射点到所述雷达的距离,na和nr分别表示散射点在方位向和距离向上的序号。
在本发明的一个实施例中,所述延迟值dbin通过如下公式计算,其中,Rna,nr表示第(na,nr)个散射点到所述雷达的距离,fsnin表示一个采样点所代表的距离,na和nr分别表示散射点在方位向和距离向上的序号。
为达到上述目的,本发明另一方面的实施例提出一种无线电测高雷达实时回波模拟系统,包括:第一获取模块,用于通过雷达获取与所述雷达位于预定范围内的三维地图数据;第一计算模块,用于计算所述雷达的当前发射脉冲重复周期时所述雷达的天线坐标和覆盖范围;第二获取模块,用于根据所述雷达接收的回波信号和所述覆盖范围获取与所述覆盖范围对应的三维地图数据,以得到每个散射点到天线的距离;第二计算模块,用于根据所述距离和对应的三维地图数据计算每个散射点对应的延迟值、相位值和幅度值;响应模块,用于根据所述延迟值、相位值和幅度值得到对应散射点的冲激响应,并根据所述冲激响应得到响应函数;以及处理模块,用于对所述响应函数和所述雷达的当前发射脉冲进行处理得到对应的回波信号。
根据本发明实施例的系统,通过将雷达和三维地图数据结合进行处理得到回波信号,实现了强实时和低延时,并且提高了灵活性。
在本发明的一个实施例中,对所述相应函数和雷达的发射脉冲信号的处理方法采用卷积、滤波和变频。
在本发明的一个实施例中,所述雷达到单个散射点的距离通过如下公式计算, 其中,(xna,nr,yna,nr,zna,nr)表示第(na,nr)个散射点的三维坐标,(x0,y0,z0)表示雷达的坐标,Rna,nr表示第(na,nr)个散射点到所述雷达的距离,na和nr分别表示散射点在方位向和距离向上的序号。
在本发明的一个实施例中,所述延迟值dbin通过如下公式计算,其中,Rna,nr表示第(na,nr)个散射点到所述雷达的距离,fsnin表示一个采样点所代表的距离,na和nr分别表示散射点在方位向和距离向上的序号。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明一个实施例的无线电测高雷达实时回波模拟方法的流程图;以及
图2为根据本发明一个实施例的无线电测高雷达实时回波模拟系统的框架图。
具体实施方式
下面详细描述本发明的实施例,实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
在本发明的描述中,需要理解的是,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中,“多个”的含义是两个或两个以上,除非另有明确具体的限定。
图1为本发明实施例的无线电测高雷达实时回波模拟方法的流程图。如图1所示,根据本发明实施例的无线电测高雷达实时回波模拟方法,包括以下步骤:
步骤S101,通过雷达获取与雷达位于预定范围内的三维地图数据。
具体地,通过雷达将数字三维地图数据(DEM)下载到大容量高速缓存中,并将两个天线的三维方向图下载到大容量高速缓存中。
步骤S102,计算雷达的当前发射脉冲重复周期时雷达的天线坐标和覆盖范围。
具体地,计算当前脉冲重复周期时刻雷达两个天线的坐标、姿态,以及波束覆盖范围。第一天线坐标为(x0,y0,z0),方位角为θA0,俯仰角为θP0,高度值为H0。第二天线坐标为(x1,y1,z1),方位角为θA1,俯仰角为θP1,高度值为H1。在FPGA中,通过利用内嵌CPU进行计算获得相应的数据。从地图文件中取出当前脉冲重复周期两个天线波束照射的数字高程地图数据。地图数据存放于外存,且应为同步SRAM介质。
步骤S103,根据雷达接收的回波信号和覆盖范围获取与覆盖范围对应的三维地图数据,以得到每个散射点到天线的距离。
具体地,计算每个散射点到雷达天线的距离,其中,利用定制可编程处理器实现。每个节点的处理器根据单个散射点坐标位置和雷达坐标位置计算之间距离,为保证回波相位精度,需要采用双精度格式进行计算。
单个散射点到雷达的距离通过如下公式计算, 其中,(xna,nr,yna,nr,zna,nr)表示第(na,nr)个散射点的三维坐标,(x0,y0,z0)表示雷达的坐标,Rna,nr表示第(na,nr)个散射点到雷达的距离,na和nr分别表示散射点在方位向和距离向上的序号。
步骤S104,根据距离和对应的三维地图数据计算每个散射点对应的延迟值、相位值和幅度值。
具体地,根据距离值,计算每个散射点对应的延迟值和相位值。其中,延迟值及相位值可采用查表法获取。
从片外RAM取出当前周期的数字高程地图数据,确定幅度值σna,nr。
步骤S105,根据延迟值、相位值和幅度值得到对应散射点的冲激响应,并根据冲激响应得到响应函数
具体地,计算每个散射点和雷达天线中心的角度值,从而查表得到在天线方向图中的功率Fna,nr。
根据每个散射点的延迟值、相位值和幅度值得到每个散射点的冲激响应,并利用天线方向图进行对齐加权。
步骤S106,对响应函数和雷达的当前发射脉冲进行处理得到对应的回波信号。
具体地,在零中频通过发射信号与系统响应函数进行卷积,得到当前脉冲重复周期的回波函数。为满足强实时、低延时的仿真要求,应采用时域卷积法进行计算,以最大限度提高运算效率。然后在对零中频回波信号内插滤波、上变频,得到两个天线各自对应的中频或射频回波信号实现实时回波模拟。
根据本发明实施例的方法,通过将雷达和三维地图数据结合进行处理得到回波信号,实现了强实时和低延时,并且提高了灵活性。
图2为根据本发明一个实施例的无线电测高雷达实时回波模拟系统的框架图。如图2所示,根据本发明实施例的无线电测高雷达实时回波模拟系统包括第一获取模块100、第一计算模块200、第二获取模块300、第二计算模块400、响应模块500和处理模块600。
第一获取模块100用于通过雷达获取与雷达位于预定范围内的三维地图数据。
具体地,通过雷达将数字三维地图数据(DEM)下载到大容量高速缓存中,并将两个天线的三维方向图下载到大容量高速缓存中。
第一计算模块200用于计算雷达的当前发射脉冲重复周期时雷达的天线坐标和覆盖范围。
具体地,计算当前脉冲重复周期时刻雷达两个天线的坐标、姿态,以及波束覆盖范围。第一天线坐标为(x0,y0,z0),方位角为θA0,俯仰角为θP0,高度值为H0。第二天线坐标为(x1,y1,z1),方位角为θA1,俯仰角为θP1,高度值为H1。在FPGA中,通过利用内嵌CPU进行计算获得相应的数据。从地图文件中取出当前脉冲重复周期两个天线波束照射的数字高程地图数据。地图数据存放于外存,且应为同步SRAM介质。
第二获取模块300用于根据雷达接收的回波信号和覆盖范围获取与覆盖范围对应的三维地图数据,以得到每个散射点到天线的距离。
具体地,计算每个散射点到雷达天线的距离,其中,利用定制可编程处理器实现。每个节点的处理器根据单个散射点坐标位置和雷达坐标位置计算之间距离,为保证回波相位精度,需要采用双精度格式进行计算。
单个散射点到雷达的距离通过如下公式计算, 其中,(xna,nr,yna,nr,zna,nr)表示第(na,nr)个散射点的三维坐标,(x0,y0,z0)表示雷达的坐标,Rna,nr表示第(na,nr)个散射点到雷达的距离,na和nr分别表示散射点在方位向和距离向上的序号。
第二计算模块400用于根据距离和对应的三维地图数据计算每个散射点对应的延迟值、相位值和幅度值。
具体地,根据距离值,计算每个散射点对应的延迟值和相位值。其中,延迟值及相位值可采用查表法获取。
从片外RAM取出当前周期的数字高程地图数据,确定幅度值σna,nr。
响应模块500用于根据延迟值、相位值和幅度值得到对应散射点的冲激响应,并根据冲激响应得到响应函数。
具体地,计算每个散射点和雷达天线中心的角度值,从而查表得到在天线方向图中的功率Fna,nr。
根据每个散射点的延迟值、相位值和幅度值得到每个散射点的冲激响应,并利用天线方向图进行对齐加权。
处理模块600用于对响应函数和雷达的当前发射脉冲进行处理得到对应的回波信号。
具体地,在零中频通过发射信号与系统响应函数进行卷积,得到当前脉冲重复周期的回波函数。为满足强实时、低延时的仿真要求,应采用时域卷积法进行计算,以最大限度提高运算效率。然后在对零中频回波信号内插滤波、上变频,得到两个天线各自对应的中频或射频回波信号实现实时回波模拟。
根据本发明实施例的系统,通过将雷达和三维地图数据结合进行处理得到回波信号,实现了强实时和低延时,并且提高了灵活性。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (10)
1.一种无线电测高雷达实时回波模拟方法,其特征在于,包括以下步骤:
通过雷达获取与所述雷达位于预定范围内的三维地图数据;
计算所述雷达的当前发射脉冲重复周期时所述雷达的天线坐标和覆盖范围;
根据所述雷达接收的回波信号和所述覆盖范围获取与所述覆盖范围对应的三维地图数据,以得到每个散射点到天线的距离;
根据所述距离和对应的三维地图数据计算每个散射点对应的延迟值、相位值和幅度值;
根据所述延迟值、相位值和幅度值得到对应散射点的冲激响应,并根据所述冲激响应得到响应函数;以及
对所述响应函数和所述雷达的当前发射脉冲进行处理得到对应的回波信号。
2.如权利要求1所述的无线电测高雷达实时回波模拟方法,其特征在于,对所述相应函数和雷达的发射脉冲信号的处理方法采用卷积、滤波和变频。
3.如权利要求1所述的无线电测高雷达实时回波模拟方法,其特征在于,所述雷达到单个散射点的距离通过如下公式计算,
其中,(xna,nr,yna,nr,zna,nr)表示第(na,nr)个散射点的三维坐标,(x0,y0,z0)表示雷达的坐标,Rna,nr表示第(na,nr)个散射点到所述雷达的距离,na和nr分别表示散射点在方位向和距离向上的序号。
4.如权利要求1所述的无线电测高雷达实时回波模拟方法,其特征在于,所述延迟值dxin通过如下公式计算,
其中,Rna,nr表示第(na,nr)个散射点到所述雷达的距离,fsnin表示一个采样点所代表的距离,na和nr分别表示散射点在方位向和距离向上的序号。
6.一种无线电测高雷达实时回波模拟系统,其特征在于,包括:
第一获取模块,用于通过雷达获取与所述雷达位于预定范围内的三维地图数据;
第一计算模块,用于计算所述雷达的当前发射脉冲重复周期时所述雷达的天线坐标和覆盖范围;
第二获取模块,用于根据所述雷达接收的回波信号和所述覆盖范围获取与所述覆盖范围对应的三维地图数据,以得到每个散射点到天线的距离;
第二计算模块,用于根据所述距离和对应的三维地图数据计算每个散射点对应的延迟值、相位值和幅度值;
响应模块,用于根据所述延迟值、相位值和幅度值得到对应散射点的冲激响应,并根据所述冲激响应得到响应函数;以及
处理模块,用于对所述响应函数和所述雷达的当前发射脉冲进行处理得到对应的回波信号。
7.如权利要求6所述的无线电测高雷达实时回波模拟系统,其特征在于,所述处理模块对所述相应函数和雷达的发射脉冲信号的处理方法采用卷积、滤波和变频。
8.如权利要求6所述的无线电测高雷达实时回波模拟系统,其特征在于,所述雷达到单个散射点的距离通过如下公式计算,
其中,(xna,nr,yna,nr,zna,nr)表示第(na,nr)个散射点的三维坐标,(x0,y0,z0)表示雷达的坐标,Rna,nr表示第(na,nr)个散射点到所述雷达的距离,na和nr分别表示散射点在方位向和距离向上的序号。
9.如权利要求6所述的无线电测高雷达实时回波模拟系统,其特征在于,所述延迟值dbin通过如下公式计算,
其中,Rna,nr表示第(na,nr)个散射点到所述雷达的距离,fsnin表示一个采样点所代表的距离,na和nr分别表示散射点在方位向和距离向上的序号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310095262.6A CN103163504B (zh) | 2013-03-22 | 2013-03-22 | 无线电测高雷达实时回波模拟方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310095262.6A CN103163504B (zh) | 2013-03-22 | 2013-03-22 | 无线电测高雷达实时回波模拟方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103163504A true CN103163504A (zh) | 2013-06-19 |
CN103163504B CN103163504B (zh) | 2014-07-16 |
Family
ID=48586750
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310095262.6A Active CN103163504B (zh) | 2013-03-22 | 2013-03-22 | 无线电测高雷达实时回波模拟方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103163504B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105093203A (zh) * | 2015-08-31 | 2015-11-25 | 中国科学院电子学研究所 | 星体次表层遥感探测雷达回波模拟及参数反演方法 |
CN106707254A (zh) * | 2017-03-31 | 2017-05-24 | 北京润科通用技术有限公司 | 一种雷达回波信号的模拟方法及装置 |
CN108614248A (zh) * | 2018-04-11 | 2018-10-02 | 南京理工大学 | 基于分形理论的调频探测目标回波模拟方法 |
CN109239684A (zh) * | 2018-10-15 | 2019-01-18 | 中国航天科工集团八五研究所 | 基于电磁散射模型的雷达目标回波模拟系统及其模拟方法 |
CN110208760A (zh) * | 2019-05-27 | 2019-09-06 | 西安空间无线电技术研究所 | 一种基于时域上采样的雷达回波仿真方法 |
CN111551905A (zh) * | 2020-04-20 | 2020-08-18 | 北京无线电计量测试研究所 | 一种散射点延时参数校准方法和系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102778241A (zh) * | 2012-07-31 | 2012-11-14 | 清华大学 | 通用机载无线电高度表回波模拟器、模拟回波生成方法 |
CN102778673A (zh) * | 2012-07-24 | 2012-11-14 | 清华大学 | 雷达回波信号高精度距离模拟方法、装置及目标模拟器 |
CN102830390A (zh) * | 2012-08-24 | 2012-12-19 | 电子科技大学 | 一种双基地sar回波模拟方法 |
CN103176179A (zh) * | 2013-02-22 | 2013-06-26 | 清华大学 | Sar雷达回波实时模拟方法 |
-
2013
- 2013-03-22 CN CN201310095262.6A patent/CN103163504B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102778673A (zh) * | 2012-07-24 | 2012-11-14 | 清华大学 | 雷达回波信号高精度距离模拟方法、装置及目标模拟器 |
CN102778241A (zh) * | 2012-07-31 | 2012-11-14 | 清华大学 | 通用机载无线电高度表回波模拟器、模拟回波生成方法 |
CN102830390A (zh) * | 2012-08-24 | 2012-12-19 | 电子科技大学 | 一种双基地sar回波模拟方法 |
CN103176179A (zh) * | 2013-02-22 | 2013-06-26 | 清华大学 | Sar雷达回波实时模拟方法 |
Non-Patent Citations (1)
Title |
---|
孟庆虎等: "一种基于FPGA的通用雷达回波实时模拟器", 《电子技术应用》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105093203A (zh) * | 2015-08-31 | 2015-11-25 | 中国科学院电子学研究所 | 星体次表层遥感探测雷达回波模拟及参数反演方法 |
CN106707254A (zh) * | 2017-03-31 | 2017-05-24 | 北京润科通用技术有限公司 | 一种雷达回波信号的模拟方法及装置 |
CN106707254B (zh) * | 2017-03-31 | 2019-11-19 | 北京润科通用技术有限公司 | 一种雷达回波信号的模拟方法及装置 |
CN108614248A (zh) * | 2018-04-11 | 2018-10-02 | 南京理工大学 | 基于分形理论的调频探测目标回波模拟方法 |
CN108614248B (zh) * | 2018-04-11 | 2022-04-08 | 南京理工大学 | 基于分形理论的调频探测目标回波模拟方法 |
CN109239684A (zh) * | 2018-10-15 | 2019-01-18 | 中国航天科工集团八五研究所 | 基于电磁散射模型的雷达目标回波模拟系统及其模拟方法 |
CN109239684B (zh) * | 2018-10-15 | 2020-09-11 | 中国航天科工集团八五一一研究所 | 基于电磁散射模型的雷达目标回波模拟系统及其模拟方法 |
CN110208760A (zh) * | 2019-05-27 | 2019-09-06 | 西安空间无线电技术研究所 | 一种基于时域上采样的雷达回波仿真方法 |
CN110208760B (zh) * | 2019-05-27 | 2021-07-13 | 西安空间无线电技术研究所 | 一种基于时域上采样的雷达回波仿真方法 |
CN111551905A (zh) * | 2020-04-20 | 2020-08-18 | 北京无线电计量测试研究所 | 一种散射点延时参数校准方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103163504B (zh) | 2014-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103163504B (zh) | 无线电测高雷达实时回波模拟方法及系统 | |
CN103176179B (zh) | Sar雷达回波实时模拟方法 | |
CN105259552A (zh) | 一种基于非线性调频信号的合成孔径雷达成像方法和装置 | |
CN105242243B (zh) | 基于去斜处理和两次延时的宽带接收数字波束形成方法 | |
CN105204004B (zh) | 基于数字延时和相位补偿的发射数字波束形成方法 | |
CN103064080B (zh) | 一种连续波目标引导雷达 | |
CN110515050B (zh) | 一种基于gpu的星载sar实时回波模拟器 | |
CN109188434B (zh) | 基于调频连续波体制的sar系统及其处理方法 | |
CN102680975B (zh) | 一种天基步进频分时测角雷达空间非合作目标成像方法 | |
CN103095220B (zh) | 基于快行fir滤波器的微型sar数字下变频器设计方法 | |
CN103308917B (zh) | 用于实时测量超低空段风的声探测系统 | |
RU2546999C1 (ru) | Короткоимпульсный радиолокатор с электронным сканированием в двух плоскостях и с высокоточным измерением координат и скорости объектов | |
CN115336420B (zh) | 底视合成孔径与近底视双刈幅干涉成像的一体化探测装置 | |
CN111649803B (zh) | 基于垂直线性阵列的三维雷达物位计及其设计方法 | |
CN105676190B (zh) | 一种校正合成孔径雷达回波数据的方法和装置 | |
CN102914775A (zh) | 一种改进的双基地合成孔径雷达nlcs成像算法 | |
CN103872463A (zh) | 一种具有高辐射效率的小型s波段无源相控阵天线 | |
CN110988862A (zh) | 一种基于极近距离毫米波雷达感知方法及系统 | |
CN103630905A (zh) | 阵列天线sar极坐标交叠子孔径成像方法 | |
CN108254728B (zh) | 基于局部线性误差模型的调频连续波sar运动补偿方法 | |
CN104237886A (zh) | 一种高精度合成孔径雷达成像方法 | |
CN104035094A (zh) | 接收距离向多孔径宽幅星载sar回波的方法和装置 | |
CN106772326A (zh) | 一种多子阵合成孔径声纳相位误差分析方法 | |
CN103675781A (zh) | 一种准确的近场回波获取方法 | |
CN107390210B (zh) | 一种物位测量中差拍信号的数字处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |