CN110376613B - 一种基于无人机磁航向的gps欺骗干扰检测方法 - Google Patents

一种基于无人机磁航向的gps欺骗干扰检测方法 Download PDF

Info

Publication number
CN110376613B
CN110376613B CN201910622580.0A CN201910622580A CN110376613B CN 110376613 B CN110376613 B CN 110376613B CN 201910622580 A CN201910622580 A CN 201910622580A CN 110376613 B CN110376613 B CN 110376613B
Authority
CN
China
Prior art keywords
unmanned aerial
aerial vehicle
section
gps
detection
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.)
Active
Application number
CN201910622580.0A
Other languages
English (en)
Other versions
CN110376613A (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.)
Northwestern Polytechnical University
Xian Aisheng Technology Group Co Ltd
Original Assignee
Northwestern Polytechnical University
Xian Aisheng Technology Group Co Ltd
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 Northwestern Polytechnical University, Xian Aisheng Technology Group Co Ltd filed Critical Northwestern Polytechnical University
Priority to CN201910622580.0A priority Critical patent/CN110376613B/zh
Publication of CN110376613A publication Critical patent/CN110376613A/zh
Application granted granted Critical
Publication of CN110376613B publication Critical patent/CN110376613B/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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/21Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
    • G01S19/215Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service issues related to spoofing

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

本发明涉及一种基于无人机磁航向的GPS欺骗干扰检测方法,以无人机当前位置为起点A,当前航向
Figure DDA0002125999730000011
为方向,构建虚拟航线AB段、BC段、CD段和DA段,引导无人机沿AB段进行飞行并对飞行耗时Tf进行记录,若无人机偏离虚拟航线的距离Dd小于等于门限值Kd,每隔1秒计算一次无人机当前磁航向与角度
Figure DDA0002125999730000012
之差,记为γ,当无人机到达B点后,对此段时间内的角度差进行均值滤波,得到平均偏差,并引导无人机进入BC航段;若Dd大于Kd,则开始计时,若持续时间超过5秒,则判定为无人机GPS受到干扰并退出检测;若飞行耗时Tf>2t仍未判断到达B点,则判定为无人机GPS受到干扰并退出检测;依次进行BC段、CD段和DA段的判断。

Description

一种基于无人机磁航向的GPS欺骗干扰检测方法
技术领域
本发明涉及无人机导航领域,更具体地说,是指一种基于无人机磁航向的GPS欺骗干扰检测方法。
背景技术
目前大多数无人机的导航都是建立在GPS接收机基础上,地面操纵人员根据GPS接收机输出的位置信息确定无人机位置、速度及航迹角,控制无人机飞往指定区域,当地面操纵人员失去对无人机的监控链路后,无人机的自主返航过程也是依靠GPS接收机输出的位置信息来导航至预定的区域,因此,GPS接收机输出数据的正确性对无人机导航至关重要。
GPS干扰分为压制式干扰和欺骗式干扰。欺骗式干扰会使GPS接收机在无意识的情况下接收到干扰信号,从而解算出错误的定位数据,进而导致采用该定位数据进行导航的无人机受到影响,不能正确进行导航。这种干扰在无人机受控情况下很容易能被操纵人员觉察从而采取应对措施,而在无人机自主返航这种应急模式下,必须研究能够检测GPS欺骗干扰的方法,保证无人机不受影响,顺利返航。
磁航向传感器是根据地球磁场信息来确定方位的设备,对于采用磁航向传感器获取航向的无人机,其航向信息不会受到GPS欺骗干扰的影响,因此可以在无人机自主返航时利用该信息来进行欺骗干扰检测,从而及时采取应对措施,避免无人机被欺骗、被干扰。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种基于无人机磁航向的GPS欺骗干扰检测方法,实现无人机在失去控制应急返航这种模式下对GPS欺骗干扰的检测,及时为自主导航控制提供欺骗警示,从而降低无人机受欺骗风险。
技术方案
一种基于无人机磁航向的GPS欺骗干扰检测方法,其特征在于步骤如下:
步骤1:无人机失控后,自动进入应急返航模式,无人机按预定航线进行返航,此时启动无人机GPS欺骗干扰检测,每T分钟进行一次判断,进入应急返航模式后立即进行一次检测;
步骤2:以无人机当前位置(x0,y0)为起点A,当前航向
Figure BDA0002125999710000023
为方向,构建虚拟航线AB段、BC段、CD段和DA段,其中B点位置为(XB,YB),C点位置为(XC,YC),D点位置为(XD,YD),各点位置按以下公式进行计算:
Figure BDA0002125999710000021
式中,d为虚拟航线段的长度,v为无人机巡航速度,t为单条检测航线段飞行时长;
步骤3:引导无人机沿AB段进行飞行并对飞行耗时Tf进行记录,若无人机偏离虚拟航线的距离Dd小于等于门限值Kd,每隔1秒计算一次无人机当前磁航向与角度
Figure BDA0002125999710000024
之差,记为γ,当无人机到达B点后,对此段时间内的角度差进行均值滤波,得到平均偏差
Figure BDA0002125999710000022
并引导无人机进入BC航段;若Dd大于Kd,则开始计时,若持续时间超过5秒,则判定为无人机GPS受到干扰并退出检测;若飞行耗时Tf>2t仍未判断到达B点,则判定为无人机GPS受到干扰并退出检测;
步骤4:引导无人机沿BC段进行飞行,按照步骤3的方法继续进行判断,得到在BC段的平均偏差
Figure BDA0002125999710000031
引导无人机进入CD航段或判定为受到干扰并退出检测;
步骤5:引导无人机沿CD段进行飞行,按照步骤3的方法继续进行判断,得到在CD段的平均偏差
Figure BDA0002125999710000032
引导无人机进入DA航段或判定为受到干扰并退出检测;
步骤6:引导无人机沿DA段进行飞行,按照步骤3的方法继续进行判断,得到在DA段的平均偏差
Figure BDA0002125999710000033
到达A点后,计算无人机当前位置(x1,y1)与(x0,y0)的距离
Figure BDA0002125999710000034
Figure BDA0002125999710000035
若该距离大于门限值Ks,则判定为无人机GPS受到干扰并退出检测;若在每个航段均未检测到受干扰,则计算AB段与CD段平均偏差之和γ1=γABCD以及BC段与DA段平均偏差之和γ2=γBCDA,若γ1和γ2任意一个大于门限值Kγ,则判定为无人机GPS受到干扰并退出检测;
步骤7:当判定为无人机GPS未受到干扰,则引导无人机继续沿原返航航线进行返航;当判定为无人机GPS受到干扰,则不再使用GPS位置信息进行飞行。
步骤1中的T取无人机飞行10km~20km所需时长。
步骤2中的t取值为1~2分钟。
步骤3中的Kd取值为100m~500m。
步骤6中的Ks取值为100m~1000m。
步骤6中的Kγ取值为1度~5度。
有益效果
本发明提出的一种基于无人机磁航向的GPS欺骗干扰检测方法,具有以下3点优点:
1、借助于不受GPS信号干扰的磁航向传感器作为检测基础,依据无人机自身控制策略,不受时间空间因素限制且能保证欺骗检测的可信度。
2、通过构建虚拟的正四边形检测航线,能有效降低因空中风所带来的航向偏差影响,提高检测的准确性。
3、检测过程中多个门限值可调,使用灵活,可针对多种无人机进行适配。
附图说明
图1构建虚拟检测航线示意图
图2检测流程图
具体实施方式
现结合实施例、附图对本发明作进一步描述:
本发明所涉方法步骤如下:
步骤1:无人机失控后,自动进入应急返航模式,无人机按预定航线进行返航,此时启动无人机GPS欺骗干扰检测,每T分钟进行一次判断,该时间取无人机飞行10km~20km所需时长,进入应急返航模式后立即进行一次检测;
步骤2:以无人机当前位置(x0,y0)为起点A,当前航向
Figure BDA0002125999710000042
为方向,构建虚拟航线AB段,BC段,CD段和DA段,其中B点位置为(XB,YB),C点位置为(XC,YC),D点位置为(XD,YD),各点位置按以下公式进行计算:
Figure BDA0002125999710000041
上式中d为虚拟航线段的长度,v为无人机巡航速度,t为单条检测航线段飞行时长,t应保证无人机进入航线后至少能稳定飞行30秒,一般取值为1~2分钟。
步骤3:引导无人机沿AB段进行飞行并对飞行耗时Tf进行记录,若无人机偏离虚拟航线的距离Dd小于等于门限值Kd,该门限值根据无人机的抗风能力来设定,一般取值为100m~500m,每隔1秒计算一次无人机当前磁航向与角度
Figure BDA0002125999710000043
之差,记为γ,当无人机到达B点后,对此段时间内的角度差进行均值滤波,得到平均偏差
Figure BDA0002125999710000051
并引导无人机进入BC航段;若Dd大于Kd,则开始计时,若持续时间超过5秒,则判定为无人机GPS受到干扰并退出检测;若飞行耗时Tf>2t仍未判断到达B点,则判定为无人机GPS受到干扰并退出检测。
步骤4:引导无人机沿BC段进行飞行,按照步骤3的方法继续进行判断,得到在BC段的平均偏差
Figure BDA0002125999710000052
引导无人机进入CD航段或判定为受到干扰并退出检测。
步骤5:引导无人机沿CD段进行飞行,按照步骤3的方法继续进行判断,得到在CD段的平均偏差
Figure BDA0002125999710000053
引导无人机进入DA航段或判定为受到干扰并退出检测。
步骤6:引导无人机沿DA段进行飞行,按照步骤3的方法继续进行判断,得到在DA段的平均偏差
Figure BDA0002125999710000054
到达A点后,计算无人机当前位置(x1,y1)与(x0,y0)的距离
Figure BDA0002125999710000055
Figure BDA0002125999710000056
若该距离大于门限值Ks,该门限值根据无人机判断航点切换的距离决定,一般取值为100m~1000m,则判定为无人机GPS受到干扰并退出检测;若在每个航段均未检测到受干扰,则计算AB段与CD段平均偏差之和γ1=γABCD以及BC段与DA段平均偏差之和γ2=γBCDA,若γ1和γ2任意一个大于门限值Kγ,该门限值根据无人机航向控制精度来确定,一般取值为1度~5度,则判定为无人机GPS受到干扰并退出检测。
步骤7:当判定为无人机GPS未受到干扰,则引导无人机继续沿原返航航线进行返航;当判定为无人机GPS受到干扰,则不再使用GPS位置信息进行飞行。
实施例:
步骤1:无人机失控后,自动进入应急返航模式,飞行速度按巡航速度140km/h进行控制,无人机按预定航线进行返航,此时启动无人机GPS欺骗干扰检测,每8分钟进行一次判断,进入应急返航模式后立即进行一次检测;
步骤2:以无人机当前位置(x0,y0)为起点A,当前航向
Figure BDA0002125999710000061
为方向,构建虚拟航线AB段,BC段,CD段和DA段,其中B点位置为(XB,YB),C点位置为(XC,YC),D点位置为(XD,YD),各点位置按以下公式进行计算:
Figure BDA0002125999710000062
上式中d为单条虚拟航线段的长度,取值2400米,t为单条检测航线段飞行时长,取值1分钟。
步骤3:引导无人机沿AB段进行飞行并对飞行耗时Tf进行记录,若无人机偏离虚拟航线的距离Dd小于等于门限值100米,每隔1秒计算一次无人机当前磁航向与角度
Figure BDA0002125999710000063
之差,记为γ,当无人机到达B点后,对此段时间内的角度差进行均值滤波,得到平均偏差
Figure BDA0002125999710000064
并引导无人机进入BC航段;若Dd大于100米,则开始计时,若持续时间超过5秒,则判定为无人机GPS受到干扰并退出检测;若飞行耗时Tf>2分钟仍未判断到达B点,则判定为无人机GPS受到干扰并退出检测。
步骤4:引导无人机沿BC段进行飞行,按照步骤3的方法继续进行判断,得到在BC段的平均偏差
Figure BDA0002125999710000065
引导无人机进入CD航段或判定为受到干扰并退出检测。
步骤5:引导无人机沿CD段进行飞行,按照步骤3的方法继续进行判断,得到在CD段的平均偏差
Figure BDA0002125999710000066
引导无人机进入DA航段或判定为受到干扰并退出检测。
步骤6:引导无人机沿DA段进行飞行,按照步骤3的方法继续进行判断,得到在DA段的平均偏差
Figure BDA0002125999710000071
到达A点后,计算无人机当前位置(x1,y1)与(x0,y0)的距离
Figure BDA0002125999710000072
Figure BDA0002125999710000073
若该距离大于门限值300米,则判定为无人机GPS受到干扰并退出检测;若在每个航段均未检测到受干扰,则计算AB段与CD段平均偏差之和γ1=γABCD以及BC段与DA段平均偏差之和γ2=γBCDA,若γ1和γ2任意一个大于门限值5度,则判定为无人机GPS受到干扰并退出检测。
步骤7:当判定为无人机GPS未受到干扰,则引导无人机继续沿原返航航线进行返航;当判定为无人机GPS受到干扰,则不再使用GPS位置信息进行飞行。
从而实现无人机GPS诱骗干扰的检测。

Claims (6)

1.一种基于无人机磁航向的GPS欺骗干扰检测方法,其特征在于步骤如下:
步骤1:无人机失控后,自动进入应急返航模式,无人机按预定航线进行返航,此时启动无人机GPS欺骗干扰检测,每T分钟进行一次判断,进入应急返航模式后立即进行一次检测;
步骤2:以无人机当前位置(x0,y0)为起点A,当前航向
Figure FDA0003887852170000011
为方向,构建虚拟航线AB段、BC段、CD段和DA段,其中B点位置为(XB,YB),C点位置为(XC,YC),D点位置为(XD,YD),各点位置按以下公式进行计算:
Figure FDA0003887852170000012
式中,d为虚拟航线段的长度,v为无人机巡航速度,t为单条检测航线段飞行时长;
步骤3:引导无人机沿AB段进行飞行并对飞行耗时Tf进行记录,若无人机偏离虚拟航线的距离Dd小于等于门限值Kd,每隔1秒计算一次无人机当前磁航向与角度
Figure FDA0003887852170000013
之差,记为γ,当无人机到达B点后,对此段时间内的角度差进行均值滤波,得到平均偏差
Figure FDA0003887852170000014
并引导无人机进入BC航段;若Dd大于Kd,则开始计时,若持续时间超过5秒,则判定为无人机GPS受到干扰并退出检测;若飞行耗时Tf>2t仍未判断到达B点,则判定为无人机GPS受到干扰并退出检测;
步骤4:引导无人机沿BC段进行飞行,按照步骤3的方法继续进行判断,得到在BC段的平均偏差
Figure FDA0003887852170000015
引导无人机进入CD航段或判定为受到干扰并退出检测;
步骤5:引导无人机沿CD段进行飞行,按照步骤3的方法继续进行判断,得到在CD段的平均偏差
Figure FDA0003887852170000021
引导无人机进入DA航段或判定为受到干扰并退出检测;
步骤6:引导无人机沿DA段进行飞行,按照步骤3的方法继续进行判断,得到在DA段的平均偏差
Figure FDA0003887852170000022
到达A点后,计算无人机当前位置(x1,y1)与(x0,y0)的距离
Figure FDA0003887852170000023
Figure FDA0003887852170000024
若该距离大于门限值Ks,则判定为无人机GPS受到干扰并退出检测;若在每个航段均未检测到受干扰,则计算AB段与CD段平均偏差之和γ1=γABCD以及BC段与DA段平均偏差之和γ2=γBCDA,若γ1和γ2任意一个大于门限值Kθ,则判定为无人机GPS受到干扰并退出检测;
步骤7:当判定为无人机GPS未受到干扰,则引导无人机继续沿原返航航线进行返航;当判定为无人机GPS受到干扰,则不再使用GPS位置信息进行飞行。
2.根据权利要求1所述的一种基于无人机磁航向的GPS欺骗干扰检测方法,其特征在于步骤1中的T取无人机飞行10km~20km所需时长。
3.根据权利要求1所述的一种基于无人机磁航向的GPS欺骗干扰检测方法,其特征在于步骤2中的t取值为1~2分钟。
4.根据权利要求1所述的一种基于无人机磁航向的GPS欺骗干扰检测方法,其特征在于步骤3中的Kd取值为100m~500m。
5.根据权利要求1所述的一种基于无人机磁航向的GPS欺骗干扰检测方法,其特征在于步骤6中的Ks取值为100m~1000m。
6.根据权利要求1所述的一种基于无人机磁航向的GPS欺骗干扰检测方法,其特征在于步骤6中的Kθ取值为1度~5度。
CN201910622580.0A 2019-07-11 2019-07-11 一种基于无人机磁航向的gps欺骗干扰检测方法 Active CN110376613B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910622580.0A CN110376613B (zh) 2019-07-11 2019-07-11 一种基于无人机磁航向的gps欺骗干扰检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910622580.0A CN110376613B (zh) 2019-07-11 2019-07-11 一种基于无人机磁航向的gps欺骗干扰检测方法

Publications (2)

Publication Number Publication Date
CN110376613A CN110376613A (zh) 2019-10-25
CN110376613B true CN110376613B (zh) 2022-12-23

Family

ID=68250895

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910622580.0A Active CN110376613B (zh) 2019-07-11 2019-07-11 一种基于无人机磁航向的gps欺骗干扰检测方法

Country Status (1)

Country Link
CN (1) CN110376613B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112198534A (zh) * 2020-09-11 2021-01-08 中国航空工业集团公司成都飞机设计研究所 一种基于飞机控制响应的位置导航信号防欺骗方法
CN112859013B (zh) * 2021-01-22 2023-06-30 西安电子科技大学 基于航迹的无人机gnss欺骗检测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103852077A (zh) * 2014-03-19 2014-06-11 北京航空航天大学 链路失效时无人机定位信息防欺骗自主判断方法
CN105929417A (zh) * 2016-06-06 2016-09-07 华东师范大学 一种捕获无人机的方法
CN107728629A (zh) * 2017-09-19 2018-02-23 富平县韦加无人机科技有限公司 无人机磁异常检测系统及方法
KR20180020384A (ko) * 2016-08-18 2018-02-28 한화시스템 주식회사 위성 항법 기만에 의한 무인기 유도 기만 장치 및 방법
JP2018147467A (ja) * 2017-03-03 2018-09-20 アルパイン株式会社 無人航空機の飛行制御装置および飛行制御方法
CN109088695A (zh) * 2018-07-24 2018-12-25 广州海格亚华防务科技有限公司 一种无人机侦测防御系统及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8922427B2 (en) * 2011-06-29 2014-12-30 Bae Systems Information And Electronic Systems Integration Inc. Methods and systems for detecting GPS spoofing attacks

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103852077A (zh) * 2014-03-19 2014-06-11 北京航空航天大学 链路失效时无人机定位信息防欺骗自主判断方法
CN105929417A (zh) * 2016-06-06 2016-09-07 华东师范大学 一种捕获无人机的方法
KR20180020384A (ko) * 2016-08-18 2018-02-28 한화시스템 주식회사 위성 항법 기만에 의한 무인기 유도 기만 장치 및 방법
JP2018147467A (ja) * 2017-03-03 2018-09-20 アルパイン株式会社 無人航空機の飛行制御装置および飛行制御方法
CN107728629A (zh) * 2017-09-19 2018-02-23 富平县韦加无人机科技有限公司 无人机磁异常检测系统及方法
CN109088695A (zh) * 2018-07-24 2018-12-25 广州海格亚华防务科技有限公司 一种无人机侦测防御系统及方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Covert Spoofing Algorithm of UAV Based on GPS/INS-Integrated Navigation;Yan Guo Et al.;《IEEE Transactions on Vehicular Technology》;20190531;第68卷(第7期);第6557-6564页 *
无人机G PS欺骗干扰方法及诱导模型的研究;戴博文等;《通信技术》;20170331;第50卷(第3期);第496-501页 *
无人机中磁航向自修正应用;刘程等;《压电与声光》;20140630;第36卷(第3期);第389-392页 *

Also Published As

Publication number Publication date
CN110376613A (zh) 2019-10-25

Similar Documents

Publication Publication Date Title
CN109992006B (zh) 一种电力巡检无人机的精准回收方法和系统
US9939819B2 (en) System and methods for automatically landing aircraft
US4700307A (en) Feature navigation system and method
US7839322B2 (en) System for detecting obstacles in the vicinity of a touchdown point
EP3410249B1 (en) System to estimate wind direction and strength using constant bank angle turn
US6792363B1 (en) System and method for trajectory optimization using adaptive navigation performance estimation
CN110376613B (zh) 一种基于无人机磁航向的gps欺骗干扰检测方法
CN107783106B (zh) 无人机与障碍物之间的数据融合方法
CN105390028B (zh) 船舶航行轨迹的纠正方法和系统
KR101136120B1 (ko) 무인항공기 이착륙 유도 시스템
KR20150019771A (ko) 무인 항공기의 착륙 방법 및 시스템
CN107728618A (zh) 一种无人船的避障方法
CN106352872A (zh) 一种无人机自主导航系统及其导航方法
EP1947542B1 (en) Method and system for automatically guiding an unmanned vehicle
US20170162066A1 (en) Method and system for aiding landing of an aircraft
US20130158749A1 (en) Methods, systems, and apparatuses for measuring fluid velocity
CN107861521A (zh) 无人靶机地面站系统
EP3680617A1 (en) Restoring navigational performance for a navigational system
CN114995520A (zh) 一种无人机协同防御控制方法、电子设备及介质
CN106043277A (zh) 汽车自动跟车控制系统及方法、汽车自动跟车系统及方法,控制雷达转向方法
CN114924581B (zh) 一种单余度无人机俯仰角失效的判定方法
CN110307837B (zh) 一种基于图像识别的无人机导航系统及方法
EP2843364A1 (en) Method and system for estimating aircraft course
CN116430901A (zh) 一种基于移动停机坪的无人机返航控制方法及系统
RU2475802C1 (ru) Способ дистанционного управления полетом беспилотного летательного аппарата и беспилотная авиационная система

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