[go: up one dir, main page]

CN115097432A - 基于数字可编程超表面的近距离目标探测方法 - Google Patents

基于数字可编程超表面的近距离目标探测方法 Download PDF

Info

Publication number
CN115097432A
CN115097432A CN202210712830.1A CN202210712830A CN115097432A CN 115097432 A CN115097432 A CN 115097432A CN 202210712830 A CN202210712830 A CN 202210712830A CN 115097432 A CN115097432 A CN 115097432A
Authority
CN
China
Prior art keywords
signal
target
sum
difference
metasurface
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
Application number
CN202210712830.1A
Other languages
English (en)
Other versions
CN115097432B (zh
Inventor
万向
黄梓艾
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN202210712830.1A priority Critical patent/CN115097432B/zh
Publication of CN115097432A publication Critical patent/CN115097432A/zh
Application granted granted Critical
Publication of CN115097432B publication Critical patent/CN115097432B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/42Simultaneous measurement of distance and other co-ordinates
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/66Radar-tracking systems; Analogous systems
    • G01S13/72Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar

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

本发明公开了一种基于数字可编程超表面的近距离目标探测方法。探测信息包括目标水平位移、径向距离和径向速度。其中,目标水平位移估计依靠超表面生成的和差波束完成,通过差波束与和波束的比值得出差和比进行估计,进一步地,可通过不同指向的和差波束实现波束跟踪;对于径向距离和速度估计,应用基于和波束上的二相编码的距离多普勒算法实现,其中二相编码信号由超表面进行两种相同辐射模式拨通辐射相位的口面编码切换生成。

Description

基于数字可编程超表面的近距离目标探测方法
技术领域
本发明属于无线通讯技术领域,尤其涉及一种基于数字可编程超表面的近距 离目标探测方法。
背景技术
电磁感知被视作6G无线通信的重点,它包括目标定位、成像、跟踪和分类, 所述目标定位是指目标距离、速度和角度探测,这些研究目前主要基于传统相控 阵天线,由于相控阵系统需要给每一个天线安装收发器,其系统架构复杂度高, 成本高。数字可编程超表面在电磁波调控方面具有卓越的性能,与传统相控阵相 比,超表面只需要一个收发器,可以大大降低成本和硬件复杂度;另外,超表面 可以通过动态调控口面的单元状态来调控电磁波,降低了在时域、频域和空域上 进行信号处理的复杂度,对超表面在电磁感知的应用带来极大的便利。
发明内容
本发明的目的在于提出一种基于数字可编程超表面的近距离目标探测方法。 探测信息包括目标水平位移、径向距离和径向速度。其中,目标水平位移估计依 靠超表面生成的和差波束完成,通过差波束与和波束的比值得出差和比进行估计, 进一步地,可通过不同指向的和差波束实现波束跟踪;对于径向距离和速度估计, 应用基于和波束上的二相编码的距离多普勒算法实现,其中二相编码信号由超表 面进行两种相同辐射模式拨通辐射相位的口面编码切换生成。
本发明的基于数字可编程超表面的近距离目标探测方法,包括如下步骤:
步骤1,超表面发射波束对动目标进行照射,通过设计口面编码,使超表面 生成和差波束对动目标的水平位移进行探测;
对于和波束,辐射模式用定向波描述,定向波为
Figure BDA0003707523710000011
其中, θ和
Figure BDA0003707523710000012
分别代表定向波束的俯仰角和方位角。
对于差波束,辐射模式由两个定向波束的差值构成:
Figure BDA0003707523710000013
步骤2,计算差和比R(d),如下:
Figure BDA0003707523710000014
其中,FΔ(d)和F(d)分别表示差波束与和波束在目标移动迹线处的场强。
步骤3,对于近距离的在与超表面平行的水平线上的动目标,通过差和比获 得鉴距离线,即获得了差和比R(d)和水平位移d的关系,从而获得水平位移。
本发明的基于数字可编程超表面的近距离目标探测方法,还包括如下步骤:
步骤1,假设超表面发射带二进制辐射相位的定向波束,发射信号st(t)和接 收信号sr(t)由下列式子表示:
Figure BDA0003707523710000021
Figure BDA0003707523710000022
其中,a(t)和
Figure BDA0003707523710000023
分别是发射信号的幅度和相位。对于二相编码信号,
Figure BDA0003707523710000024
是0 或π,fc表示信号的载波频率,tr表示信号从发射到被接收经过的延时;
步骤2,脉冲压缩技术通过匹配滤波获得瞬时信噪比的信号scr(t),目标的 径向距离通过脉冲压缩技术进行估计,下式给出脉冲压缩过程:
Figure BDA0003707523710000025
其中,
Figure BDA0003707523710000026
Figure BDA0003707523710000027
分别为st(t)和sr(t)的傅里叶变换形式,信号scr(t)在t=tr时刻取得峰值;
步骤3,对于动目标,回波信号包括径向速度的信息。假设目标以速度vr靠 近超表面,目标与超表面的初始径向距离为R0,则回波信号sr(t)可以表示为:
Figure BDA0003707523710000028
其中,
Figure BDA0003707523710000029
从式(6)可知sr(t)包含延时信息tr和多普勒频率fd,延时信 息tr和多普勒频率fd分别与目标的径向距离和径向速度相关。
信号srd(t,f)在t=tr,f=fd时取得峰值,利用
Figure BDA00037075237100000210
以及d=c*tr能够实现 对动目标的径向距离和速度的估计。
进一步的,fd的值可通过多普勒滤波收集多个周期的回波信号并进行傅里叶 变换压缩回波信号进行估计,如下式所述:
srd(t,f)=FFT[scr(t+nTr)],n=0,1,…,N (7)
其中,t∈(0,Tr),
Figure BDA00037075237100000211
Tr是信号周期,N是信号累积周期数。
有益效果:超表面通过动态调控口面的单元状态来调控电磁波,实现了目标定位,降低了在时域、频域和空域上进行信号处理的复杂度,对超表面在电磁感知 的应用带来极大的便利。
附图说明
图1是本发明方法所基于数字可编程超表面的目标探测模型图;
图2是超表面的单元结构特性图;
图3是三组和差波束的全波仿真结果图;
图4是超表面和差波束和鉴距离线的全波仿真结果图;
图5是基于超表面对单目标场景和多目标场景估计的仿真结果;
图6是目标水平位移估计和跟踪实验结果图;
图7是对目标径向距离和速度进行估计的仿真图。
具体实施方式
本发明的基于数字可编程超表面的近距离目标探测方法,基于如图1所示 的目标探测模型,透射式超表面包含10×10个单元,由辐射28GHz单频波的喇 叭天线作为馈源垂直照射。图2展示了所提出的超表面的单元结构特性。图2(a) 展示了单元的结构组成。图2(b)是单元在二进制状态下的透射幅度图;图2
(c)是单元在二进制状态下的透射相位图。
超表面发射波束对动目标进行照射,通过设计口面编码,使超表面生成和 差波束对动目标的水平位移进行探测。另外,超表面通过二相调控在和波束上 生成二相编码信号,因此目标的径向距离和速度可以通过距离多普勒算法获得。
和差波束通常被应用于单脉冲雷达系统进行目标跟踪,因为和差波束的回 波携带着目标空间角度信息。与传统的单脉冲雷达天线相比,超表面可以通过 对口面编码进行重构实现和差波束,而不需要和差网络。对于和波束,辐射模 式可以用定向波描述,定向波为
Figure BDA0003707523710000031
其中,θ和
Figure BDA0003707523710000032
分别代表定 向波束的俯仰角和方位角。对于差波束,辐射模式可以由两个定向波束的差值 构成:
Figure BDA0003707523710000033
当两个定向波束彼此非常接近时,相减会导致辐射方向图在指向处生成凹陷。 经过定义足够的和差波束编码可以实现和差波束转向。对应的和差波束的口面编 码可以通过几何光学原理生成。
图3给出了三组和差波束的全波仿真结果。图3中的(a)-(c)依次给出了和 波束指向-8.53°,0°和8.53°的口面编码图,图3中的(d)-(f)依次给出了 差波束指向-8.53°,0°和8.53°的口面编码图,图3中的(g)-(i)依次给出 了和波束和差波束指向-8.53°,0°和8.53°的口面编码对应的辐射模式图。 仿真结果验证了基于超表面生成和差波束的可行性。
对于近距离的在与超表面平行的水平线上的动目标,其水平位移的求法类似 于利用差波束与和波束的比值得出单脉冲鉴角曲线的方法,可以通过差和比获得 鉴距离线,即差和比R(d)如下:
Figure BDA0003707523710000041
其中,FΔ(d)和F(d)分别表示差波束和和波束在目标移动迹线处的场强。应 用图3的和差波束编码,图4给出了超表面和差波束和鉴距离线的全波仿真结果, 图4中的(a)-(c)给出了距离超表面12mm处的平行线上,和差波束指向-8.53°, 0°和8.53°的仿真场强的结果,图4中的(d)-(f)展示了鉴距离线,即差和比R和 水平位移d的关系。这些线的单调性可以用来对动目标的水平位移进行估计。
二相编码信号由于具有高时宽带宽积、脉冲压缩和抗干扰特性,在雷达探测 领域被广泛使用。通过切换适当的口面编码,超表面可以在一个固定的模式上进 行相位调制,因此二相编码信号可以直接由超表面生成。假设超表面发射带二进 制辐射相位的定向波束,发射信号st(t)和接收信号sr(t)由下列式子表示:
Figure BDA0003707523710000042
Figure BDA0003707523710000043
其中,a(t)和
Figure BDA0003707523710000044
分别是发射信号的幅度和相位。对于二相编码信号,
Figure BDA0003707523710000045
是 0或π,fc表示信号的载波频率,tr表示信号从发射到被接收经过的延时。
脉冲压缩技术通过匹配滤波获得峰瞬时信噪比的信号,目标的径向距离可以 通过脉冲压缩技术进行估计,下式给出脉冲压缩过程:
Figure BDA0003707523710000046
其中,
Figure BDA0003707523710000047
Figure BDA0003707523710000048
分别为st(t)和sr(t)的傅里叶变换形式。信号scr(t)在 t=tr时刻取得峰值,此峰值便可用于对目标径向距离进行估计。
对于动目标,回波信号包括径向速度的信息。假设目标以速度vr靠近超表面, 目标与超表面的初始径向距离为R0,则回波信号sr(t)可以表示为:
Figure BDA0003707523710000051
其中,
Figure BDA0003707523710000052
可以看出sr(t)包含延时信息tr和多普勒频率fd,它们分 别与目标的径向距离和径向速度相关。fd的值可通过多普勒滤波收集多个周期的 回波信号并进行傅里叶变换压缩回波信号进行估计,如下式所述:
srd(t,f)=FFT[scr(t+nTr)],n=0,1,…,N (7)
其中,t∈(0,Tr),
Figure BDA0003707523710000053
Tr是信号周期,N是信号累积周期数。 信号srd(t,f)在(t=tr,f=fd)取得峰值,利用
Figure BDA0003707523710000054
以及d=c*tr能够实现 对动目标的径向距离和速度的估计。
举例验证理论,分别针对单目标和多目标的场景进行了仿真。仿真中255bit 的m序列作为二相编码信号,信号的符号速率为200MHz,对应超表面口面编码 的切换速率。在单目标场景中,目标被设定与超表面间距2.1m,并以1m/s的速 度靠近超表面,即(2.1m,1m/s);在多目标场景中,设定三个目标在同一指向上, 参数分别为(2m,-1m/s)、(4m,1m/s)和(4m,-1m/s)。图5中的(a)展示了单目 标场景下的估计结果,峰值位于(2.1m,1.0463m/s),与设定的参数非常接近。 图5中的(b)给出了三个目标场景下的估计结果,峰值位于(2.025m,-1.0463m/s), (4.05m,1.0463m/s)和(4.05m,-1.0463m/s),也与设定参数相近。仿真结果验 证了1超表面进行径向距离和速度估计的有效性。
为了估计水平位移,测试了距离超表面1m的平行线上的近场。图6中的 (a)-(c)依次展示了和差波束指向-8.53°,0°和8.53°的近场测量结果,明显 看出和差波束导致了不同的近场分布,这个分布可以用来对目标水平位移进行估 计。由于超表面和馈电喇叭之间存在多次反射,测量的近场值会出现一些抖动。 图6中的(d)-(f)是基于式(2)得出的和差波束结果。实线通过最小二乘算法拟合 而成。这些线的单调性证实了目标水平位移估计的有效性。
为了对目标径向距离和速度进行估计,超表面配置了两组口面编码,这两组 编码生成指向0°的和波束但是相位信息相反。图7给出了测量的近场信息,图 7中的(a)展示了距离超表面1m的水平线上的近场幅度和相位的测量结果,可以 看出二者幅度相近但相位相反,为生成二相编码信号奠定了基础。根据m序列切 换超表面的口面编码,在超表面前方的目标会被二相编码信号照射,图7中的(b) 是距离多普勒图像的结果,估计的参数为(1.05m,1.0463m/s),非常接近设定的 参数(1m,1m/s)。
本发明成功建立了基于数字可编程超表面的近距离目标探测,包括目标水平 位移、径向距离和径向速度。依靠超表面和差波束的比值和不同指向实现了目标 水平位移估计和跟踪,依靠二相编码信号和距离多普勒算法实现目标径向距离和 径向速度估计。仿真和实验结果都证明了实验的有效性。本发明展示了数字可编 程超表面在近场探测的应用潜力,并希望进一步开拓可编程超表面在智能电磁感 知领域的应用。

Claims (3)

1.基于数字可编程超表面的近距离目标探测方法,其特征在于,包括如下步骤:
步骤1,超表面发射波束对动目标进行照射,通过设计口面编码,使超表面生成和差波束对动目标的水平位移进行探测;
对于和波束,辐射模式用定向波描述,定向波为
Figure FDA0003707523700000011
其中,θ和
Figure FDA0003707523700000012
分别代表定向波束的俯仰角和方位角;
对于差波束,辐射模式由两个定向波束的差值构成:
Figure FDA0003707523700000013
步骤2,计算差和比R(d),如下:
Figure FDA0003707523700000014
其中,FΔ(d)和F(d)分别表示差波束与和波束在目标移动迹线处的场强;
步骤3,对于近距离的在与超表面平行的水平线上的动目标,通过差和比获得鉴距离线,即获得了差和比R(d)和水平位移d的关系,从而获得水平位移。
2.基于数字可编程超表面的近距离目标探测方法,其特征在于,包括如下步骤:
步骤1,假设超表面发射带二进制辐射相位的定向波束,发射信号st(t)和接收信号sr(t)由下列式子表示:
Figure FDA0003707523700000015
Figure FDA0003707523700000016
其中,a(t)和
Figure FDA0003707523700000017
分别是发射信号的幅度和相位;对于二相编码信号,
Figure FDA0003707523700000018
是0或π,fc表示信号的载波频率,tr表示信号从发射到被接收经过的延时;
步骤2,脉冲压缩技术通过匹配滤波获得瞬时信噪比的信号scr(t),目标的径向距离通过脉冲压缩技术进行估计,下式给出脉冲压缩过程:
Figure FDA0003707523700000019
其中,
Figure FDA00037075237000000110
Figure FDA00037075237000000111
分别为st(t)和sr(t)的傅里叶变换形式,信号scr(t)在t=tr时刻取得峰值;
步骤3,对于动目标,回波信号包括径向速度的信息;假设目标以速度vr靠近超表面,目标与超表面的初始径向距离为R0,则回波信号sr(t)可以表示为:
Figure FDA00037075237000000112
其中,
Figure FDA00037075237000000113
从式(6)可知sr(t)包含延时信息tr和多普勒频率fd,延时信息tr和多普勒频率fd分别与目标的径向距离和径向速度相关;
信号srd(t,f)在t=tr,f=fd时取得峰值,利用
Figure FDA0003707523700000021
以及d=c*tr能够实现对动目标的径向距离和速度的估计。
3.根据权利要求2所述基于数字可编程超表面的近距离目标探测方法,其特征在于,
fd的值可通过多普勒滤波收集多个周期的回波信号并进行傅里叶变换压缩回波信号进行估计,如下式所述:
srd(t,f)=FFT[scr(t+nTr)],n=0,1,...,N (7)
其中,t∈(0,Tr),
Figure FDA0003707523700000022
Tr是信号周期,N是信号累积周期数。
CN202210712830.1A 2022-06-22 2022-06-22 基于数字可编程超表面的近距离目标探测方法 Active CN115097432B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210712830.1A CN115097432B (zh) 2022-06-22 2022-06-22 基于数字可编程超表面的近距离目标探测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210712830.1A CN115097432B (zh) 2022-06-22 2022-06-22 基于数字可编程超表面的近距离目标探测方法

Publications (2)

Publication Number Publication Date
CN115097432A true CN115097432A (zh) 2022-09-23
CN115097432B CN115097432B (zh) 2025-02-11

Family

ID=83293699

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210712830.1A Active CN115097432B (zh) 2022-06-22 2022-06-22 基于数字可编程超表面的近距离目标探测方法

Country Status (1)

Country Link
CN (1) CN115097432B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2574957A1 (fr) * 2011-09-30 2013-04-03 Thales Procédé d'estimation de la fréquence doppler non ambigue d'une cible mobile, notamment marine, et radar mettant en oeuvre le procédé
CN103616679A (zh) * 2013-11-19 2014-03-05 北京航空航天大学 基于差波束调制和波形分析的pd雷达测距测角方法
RU2574079C1 (ru) * 2014-10-30 2016-02-10 Федеральное государственное казенное военное образовательное учреждение высшего профессионального образования "Военная академия войсковой противовоздушной обороны Вооруженных Сил Российской Федерации имени Маршала Советского Союза А.М. Василевского" Министерства Обороны Российской Федерации Способ однозначного измерения радиальной скорости цели в когерентно-импульсной радиолокационной станции
CN112147593A (zh) * 2020-09-02 2020-12-29 西安电子科技大学 一种高速密集爆炸物破片目标的四维参数估计方法
CN113109811A (zh) * 2021-04-15 2021-07-13 东南大学 一种基于可编程超表面的二相编码一维距离成像方法
CN113820699A (zh) * 2021-10-19 2021-12-21 中安锐达(北京)电子科技有限公司 一种不对称波束情形下的一维同时多波束雷达测角方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2574957A1 (fr) * 2011-09-30 2013-04-03 Thales Procédé d'estimation de la fréquence doppler non ambigue d'une cible mobile, notamment marine, et radar mettant en oeuvre le procédé
CN103616679A (zh) * 2013-11-19 2014-03-05 北京航空航天大学 基于差波束调制和波形分析的pd雷达测距测角方法
RU2574079C1 (ru) * 2014-10-30 2016-02-10 Федеральное государственное казенное военное образовательное учреждение высшего профессионального образования "Военная академия войсковой противовоздушной обороны Вооруженных Сил Российской Федерации имени Маршала Советского Союза А.М. Василевского" Министерства Обороны Российской Федерации Способ однозначного измерения радиальной скорости цели в когерентно-импульсной радиолокационной станции
CN112147593A (zh) * 2020-09-02 2020-12-29 西安电子科技大学 一种高速密集爆炸物破片目标的四维参数估计方法
CN113109811A (zh) * 2021-04-15 2021-07-13 东南大学 一种基于可编程超表面的二相编码一维距离成像方法
CN113820699A (zh) * 2021-10-19 2021-12-21 中安锐达(北京)电子科技有限公司 一种不对称波束情形下的一维同时多波束雷达测角方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
黄梓艾: "基于可编程超表面的微波计算成像技术研究", 《信息科技》, 31 December 2023 (2023-12-31) *

Also Published As

Publication number Publication date
CN115097432B (zh) 2025-02-11

Similar Documents

Publication Publication Date Title
CN113567928A (zh) Mimo雷达装置和mimo雷达方法
CN102680975B (zh) 一种天基步进频分时测角雷达空间非合作目标成像方法
CN102590812B (zh) 基于调频连续波的sar实时成像方法
CN103969629A (zh) 一种基于主瓣杂波配准的机载雷达杂波自适应抑制方法
CN109283497B (zh) 双基地fda-mimo雷达距离欺骗干扰识别方法
CN109507669A (zh) 基于相参累积的地面微弱运动目标参数估计方法
CN105487060A (zh) 一种双通道四斜率调制的多目标提取方法
CN104730498A (zh) 基于Keystone和加权旋转FFT的目标检测方法
CN106291489B (zh) 适用于多种发射信号波形的合成孔径雷达回波仿真方法
Choi et al. High angular resolution estimation methods for vehicle FMCW radar
Hu et al. Generalized ambiguity function properties of ground-based wideband MIMO imaging radar
CN101793960B (zh) 具有啁啾特性的高精度动态门限目标监测方法
US20230305103A1 (en) Mimo radar apparatus and mimo radar method
Hussain Ambiguity functions for monostatic and bistatic radar systems using UWB throb signal
CN115097432A (zh) 基于数字可编程超表面的近距离目标探测方法
CN107942307B (zh) 调频连续波着陆雷达波束的地面入射角估计方法
CN104360333A (zh) 一种同时校正一阶和二阶距离徙动的相参积累检测方法
CN112068112B (zh) 一种宽窄脉冲测角方法
Zhao et al. Using sky-wave echoes information to extend HFSWR's maximum detection range
CN115629378A (zh) 非视距高分辨率微波强度相干凝视成像装置及方法
Hongtu et al. Efficient Raw Signal Generation Based on Equivalent Scatterer and Subaperture Processing for SAR with Arbitrary Motion.
Wang et al. Bistatic ISAR imaging based on BP algorithm
Bayat et al. Target detection by passive coherent FM based bistatic radar
Hussain et al. Space-time array processing based on ultrawideband throb signal
CN101509973A (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