CN107727877A - 一种基于仪表着陆系统的地速测量方法 - Google Patents

一种基于仪表着陆系统的地速测量方法 Download PDF

Info

Publication number
CN107727877A
CN107727877A CN201710785737.2A CN201710785737A CN107727877A CN 107727877 A CN107727877 A CN 107727877A CN 201710785737 A CN201710785737 A CN 201710785737A CN 107727877 A CN107727877 A CN 107727877A
Authority
CN
China
Prior art keywords
mrow
msub
mtd
ground velocity
mtr
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.)
Pending
Application number
CN201710785737.2A
Other languages
English (en)
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.)
Luoyang Institute of Electro Optical Equipment AVIC
Original Assignee
Luoyang Institute of Electro Optical Equipment AVIC
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 Luoyang Institute of Electro Optical Equipment AVIC filed Critical Luoyang Institute of Electro Optical Equipment AVIC
Priority to CN201710785737.2A priority Critical patent/CN107727877A/zh
Publication of CN107727877A publication Critical patent/CN107727877A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P3/00Measuring linear or angular speed; Measuring differences of linear or angular speeds

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明涉及一种基于仪表着陆系统的地速测量方法,属于航电显示控制技术领域。该方法首先利用姿态航向基准系统输出的垂直加速度、大气数据计算机/无线电高度表输出的飞机高度、飞行管理系统输出的机场海拔高度,通过互补滤波计算飞机垂直速度;再利用几何关系计算基于仪表着陆系统的地速;最后利用沿跑道方向的加速度对地速进行滤波处理,获得平滑可用的地速。本发明的一种基于仪表着陆系统的地速测量方法,综合利用垂直加速度、飞机高度、机场海拔高度、下滑道偏差、参考下滑道角度计算姿态航向基准系统无法输出的地速的同时,通过加速度互补滤波,抑制了下滑道偏差引入的噪声,获得的地速平滑、完整性较高,可以支持平视飞行导引等高安全性应用。

Description

一种基于仪表着陆系统的地速测量方法
技术领域
本发明涉及一种基于仪表着陆系统的地速测量方法,属于航电显示控制技术领域。
背景技术
飞机的地速是平视飞行导引功能中进近导引、滑跑导引、起飞导引重要的输入参数,该参数的精度和完整性直接影响到平视飞行导引的性能和可用性。现有的姿态航向基准系统无法输出地速,GPS地速虽然在精度方面满足平视飞行导引功能的需求,但GPS信号容易受到干扰,信号的完整性无法支持平视飞行导引这种高安全应用,需要对GPS地速信号的完整性进行监控。
发明内容
本发明的目的是提供一种基于仪表着陆系统的地速测量方法,以支持平视飞行导引功能等高安全性应用。
本发明为解决上述技术问题提供了一种基于仪表着陆系统的地速测量方法,该方法包括以下步骤:
步骤1:采集无线电高度、气压高度、机场海拔高度,计算飞机离地高度,再采集飞机垂直加速度,利用互补滤波方法获得飞机的垂直速度;
步骤2:采集飞机相对标准下滑道的垂直偏差、设定的参考下滑道角度,结合飞机垂直速度,通过几何关系计算基于仪表着陆系统的地速;
步骤3:根据飞机北东地坐标系下的三轴加速度、跑道航向、基于仪表着陆系统的地速,利用互补滤波方法获得平滑处理后的地速。
进一步的优选方案,所述一种基于仪表着陆系统的地速测量方法,其特征在于:步骤1中飞机垂直速度的计算过程为:
步骤1.1:根据气压高度hbaro,无线电高度hradalt,机场海拔高度hR,计算飞机离地高度:
hagl=(hbaro-hR)N+hradalt(1-N)
其中,N为高度转换过渡时间常数,飞机进近着陆过程中,当hradalt>90ft时,N=1,当hradalt≤90ft时,N在Δt的时间内从1线性地减小到0;
步骤1.2:根据离地高度hagl,垂直加速度az_AHRS,利用互补滤波方法获得飞机的垂直速度:
其中,Tv为垂直速度互补滤波时间常数,s为拉普拉斯算子。
进一步的优选方案,所述一种基于仪表着陆系统的地速测量方法,其特征在于:步骤2中基于仪表着陆系统的地速的计算过程为:
根据下滑道偏差εGS、设定的参考下滑道角度γREF、垂直速度Vver,通过几何关系计算基于仪表着陆系统的地速VGS
进一步的优选方案,所述一种基于仪表着陆系统的地速测量方法,其特征在于:步骤3中地速平滑处理的计算过程为:
步骤3.1:根据飞机北东地坐标系下的三轴加速度ax_NED、ay_NED、az_NED,跑道航向ψR,根据以下公式计算沿跑道方向的加速度ax_RW
步骤3.2:根据沿跑道方向的加速度ax_RW,基于仪表着陆系统的地速VGS,利用互补滤波方法获得平滑处理后的地速VGS_smooth
其中,Ta为地速平滑互补滤波时间常数。
有益效果
本发明的有益效果是:本发明首先采集无线电高度、气压高度、机场海拔高度,计算飞机离地高度,再采集飞机垂直加速度,利用互补滤波方法获得飞机的垂直速度;再基于飞机精密进近过程中在进近区域紧密跟踪下滑路径这种假设,利用几何关系计算基于仪表着陆系统的地速;最后利用沿跑道方向的加速度进行互补滤波,以获得平滑处理后的地速。本发明的一种基于仪表着陆系统的地速测量方法,综合利用垂直加速度、飞机高度、机场海拔高度、下滑道偏差、参考下滑道角度,计算姿态航向基准系统无法输出的地速,同时通过加速度互补滤波,抑制了下滑道偏引入的噪声,获得的地速平滑、完整性较高,可以支持平视飞行导引等高安全性应用。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明基于仪表着陆系统的地速测量方法的原理示意图;
图2是本发明实施例中基于仪表着陆系统的地速测量方法的流程图。
具体实施方式
下面详细描述本发明的实施例,所述实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
本发明基于仪表着陆系统的地速测量方法的原理如图1所示,首先利用姿态航向基准系统输出的垂直加速度、大气数据计算机/无线电高度表输出的飞机高度、飞行管理系统输出的机场海拔高度,通过互补滤波计算飞机垂直速度;再基于飞机精密进近过程中精确跟踪下滑路径这种假设,利用几何关系计算基于仪表着陆系统的地速;最后利用沿跑道方向的加速度对地速进行滤波处理,获得平滑可用的地速。
该方法的具体实施过程如图2所示,具体实施例步骤如下:
1.解算飞机垂直速度
本步骤需要采集飞机气压高度、无线电高度、机场海拔高度、垂直加速度。
本实施例是从大气数据计算机获取气压高度hbaro,从无线电高度表获取无线电高度hradalt,从飞行管理系统获取机场海拔高度hR,从姿态航向基准系统获取飞机垂直加速度az_AHRS
首先利用气压高度hbaro、无线电高度hradalt、机场海拔高度hR,计算飞机离地高度hagl
hagl=(hbaro-hR)N+hradalt(1-N)
其中,N为高度转换过渡时间常数。飞机进近着陆过程中,当hradalt>90ft时,N=1,当hradalt≤90ft时,N在Δt的时间内,从1线性地减小到0。
利用飞机垂直加速度az_AHRS,结合飞机离地高度hagl,通过互补滤波获得飞机的垂直速度Vver
其中,Tv为垂直速度互补滤波时间常数,s为拉普拉斯算子。
2.解算基于仪表着陆系统的地速
本步骤需要采集下滑道偏差、设定的参考下滑道角度。
本实施例是从多模接收机获取下滑道偏差εGS、从飞行管理系统获取设定的参考下滑道角度γREF
基于飞机精密进近过程中在进近区域紧密跟踪下滑路径这种假设,利用几何关系计算基于仪表着陆系统的地速VGS
3.地速平滑处理
本步骤需要采集北东地坐标系下的三轴加速度、跑道航向。
本实施例是从姿态航向基准系统采集北东地坐标系下的三轴加速度ax_NED、ay_NED、az_NED,从飞行管理系统采集跑道航向ψR
利用北东地坐标系和跑道坐标系之间的映射关系,计算沿跑道方向的加速度ax_RW
利用沿跑道方向的加速度ax_RW,通过互补滤波方法,抑制下滑道偏差中的噪声,获得平滑处理后的地速VGS_smooth
其中,当Ta为地速平滑互补滤波时间常数。
本实施例中基于仪表着陆系统的地速测量方法,综合利用垂直加速度、飞机气压高度、无线电高度、机场海拔高度、下滑道偏差、参考下滑道角度,计算姿态航向基准系统无法输出的地速,同时通过加速度互补滤波,抑制了下滑道偏引入的噪声,获得的地速平滑、完整性较高,可以支持平视飞行导引等高安全性应用。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (4)

1.一种基于仪表着陆系统的地速测量方法,其特征在于:包括以下步骤:
步骤1:采集无线电高度、气压高度、机场海拔高度,计算飞机离地高度,再采集飞机垂直加速度,利用互补滤波方法获得飞机的垂直速度;
步骤2:采集飞机相对标准下滑道的垂直偏差、设定的参考下滑道角度,结合飞机垂直速度,通过几何关系计算基于仪表着陆系统的地速;
步骤3:根据飞机北东地坐标系下的三轴加速度、跑道航向、基于仪表着陆系统的地速,利用互补滤波方法获得平滑处理后的地速。
2.根据权利要求1所述一种基于仪表着陆系统的地速测量方法,其特征在于:步骤1中飞机垂直速度的计算过程为:
步骤1.1:根据气压高度hbaro,无线电高度hradalt,机场海拔高度hR,计算飞机离地高度:
hagl=(hbaro-hR)N+hradalt(1-N)
其中,N为高度转换过渡时间常数,飞机进近着陆过程中,当hradalt>90ft时,N=1,当hradalt≤90ft时,N在Δt的时间内从1线性地减小到0;
步骤1.2:根据离地高度hagl,垂直加速度az_AHRS,利用互补滤波方法获得飞机的垂直速度:
<mrow> <msub> <mi>V</mi> <mrow> <mi>v</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mi>s</mi> <mrow> <msub> <mi>T</mi> <mi>v</mi> </msub> <mi>s</mi> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <msub> <mi>h</mi> <mrow> <mi>a</mi> <mi>g</mi> <mi>l</mi> </mrow> </msub> <mo>+</mo> <mfrac> <msub> <mi>T</mi> <mi>v</mi> </msub> <mrow> <msub> <mi>T</mi> <mi>v</mi> </msub> <mi>s</mi> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <msub> <mi>a</mi> <mrow> <mi>z</mi> <mo>_</mo> <mi>A</mi> <mi>H</mi> <mi>R</mi> <mi>S</mi> </mrow> </msub> </mrow>
其中,Tv为垂直速度互补滤波时间常数,s为拉普拉斯算子。
3.根据权利要求1所述一种基于仪表着陆系统的地速测量方法,其特征在于:步骤2中基于仪表着陆系统的地速的计算过程为:
根据下滑道偏差εGS、设定的参考下滑道角度γREF、垂直速度Vver,通过几何关系计算基于仪表着陆系统的地速VGS
<mrow> <msub> <mi>V</mi> <mrow> <mi>G</mi> <mi>S</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>V</mi> <mrow> <mi>v</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> <mrow> <mi>t</mi> <mi>a</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;epsiv;</mi> <mrow> <mi>G</mi> <mi>S</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>R</mi> <mi>E</mi> <mi>F</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>.</mo> </mrow>
4.根据权利要求1所述一种基于仪表着陆系统的地速测量方法,其特征在于:步骤3中地速平滑处理的计算过程为:
步骤3.1:根据飞机北东地坐标系下的三轴加速度ax_NED、ay_NED、az_NED,跑道航向ψR,根据以下公式计算沿跑道方向的加速度ax_RW
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>a</mi> <mrow> <mi>x</mi> <mo>_</mo> <mi>R</mi> <mi>W</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mrow> <mi>y</mi> <mo>_</mo> <mi>R</mi> <mi>W</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mrow> <mi>z</mi> <mo>_</mo> <mi>R</mi> <mi>W</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>cos&amp;psi;</mi> <mi>R</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>sin&amp;psi;</mi> <mi>R</mi> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>sin&amp;psi;</mi> <mi>R</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>cos&amp;psi;</mi> <mi>R</mi> </msub> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>a</mi> <mrow> <mi>x</mi> <mo>_</mo> <mi>N</mi> <mi>E</mi> <mi>D</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mrow> <mi>y</mi> <mo>_</mo> <mi>N</mi> <mi>E</mi> <mi>D</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mrow> <mi>z</mi> <mo>_</mo> <mi>N</mi> <mi>E</mi> <mi>D</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>
步骤3.2:根据沿跑道方向的加速度ax_RW,基于仪表着陆系统的地速VGS,利用互补滤波方法获得平滑处理后的地速VGS_smooth
<mrow> <msub> <mi>V</mi> <mrow> <mi>G</mi> <mi>S</mi> <mo>_</mo> <mi>s</mi> <mi>m</mi> <mi>o</mi> <mi>o</mi> <mi>t</mi> <mi>h</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msub> <mi>T</mi> <mi>a</mi> </msub> <mi>s</mi> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <msub> <mi>V</mi> <mrow> <mi>G</mi> <mi>S</mi> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mi>a</mi> </msub> <mi>s</mi> </mrow> <mrow> <msub> <mi>T</mi> <mi>a</mi> </msub> <mi>s</mi> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <msub> <mi>a</mi> <mrow> <mi>x</mi> <mo>_</mo> <mi>R</mi> <mi>W</mi> </mrow> </msub> </mrow>
其中,Ta为地速平滑互补滤波时间常数。
CN201710785737.2A 2017-09-04 2017-09-04 一种基于仪表着陆系统的地速测量方法 Pending CN107727877A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710785737.2A CN107727877A (zh) 2017-09-04 2017-09-04 一种基于仪表着陆系统的地速测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710785737.2A CN107727877A (zh) 2017-09-04 2017-09-04 一种基于仪表着陆系统的地速测量方法

Publications (1)

Publication Number Publication Date
CN107727877A true CN107727877A (zh) 2018-02-23

Family

ID=61205671

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710785737.2A Pending CN107727877A (zh) 2017-09-04 2017-09-04 一种基于仪表着陆系统的地速测量方法

Country Status (1)

Country Link
CN (1) CN107727877A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110488864A (zh) * 2019-08-15 2019-11-22 中国商用飞机有限责任公司 用于修正飞机的飞行控制系统中的loc信号的方法及系统
CN112230587A (zh) * 2020-10-29 2021-01-15 中国民用航空总局第二研究所 一种仪表着陆系统同频异呼运行监测系统及方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103197682A (zh) * 2013-03-22 2013-07-10 北京航空航天大学 一种无人机进场着陆下滑通道设计与下滑轨迹调整方法
CN103941742A (zh) * 2014-04-29 2014-07-23 中国科学院自动化研究所 一种无人机滑跑纠偏控制装置及方法
CN104238580A (zh) * 2014-09-30 2014-12-24 中国航天空气动力技术研究院 一种应用于无人机航空物探的低空飞行控制方法
CN105235911A (zh) * 2014-07-04 2016-01-13 空中客车运营简化股份公司 用于飞行器的着陆辅助方法和着陆辅助设备
CN105352495A (zh) * 2015-11-17 2016-02-24 天津大学 加速度与光流传感器数据融合无人机水平速度控制方法
US20160347467A1 (en) * 2015-05-28 2016-12-01 Airbus Helicopters Method of determining the longitudinal air speed and the longitudinal ground speed of a rotary wing aircraft depending on its exposure to the wind
CN106586026A (zh) * 2016-12-06 2017-04-26 中国航空工业集团公司洛阳电光设备研究所 一种飞机相对跑道侧向偏差速率的测量方法
CN106774376A (zh) * 2017-01-25 2017-05-31 上海拓攻机器人有限公司 一种无人机仿地飞行控制方法及系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103197682A (zh) * 2013-03-22 2013-07-10 北京航空航天大学 一种无人机进场着陆下滑通道设计与下滑轨迹调整方法
CN103941742A (zh) * 2014-04-29 2014-07-23 中国科学院自动化研究所 一种无人机滑跑纠偏控制装置及方法
CN105235911A (zh) * 2014-07-04 2016-01-13 空中客车运营简化股份公司 用于飞行器的着陆辅助方法和着陆辅助设备
CN104238580A (zh) * 2014-09-30 2014-12-24 中国航天空气动力技术研究院 一种应用于无人机航空物探的低空飞行控制方法
US20160347467A1 (en) * 2015-05-28 2016-12-01 Airbus Helicopters Method of determining the longitudinal air speed and the longitudinal ground speed of a rotary wing aircraft depending on its exposure to the wind
CN105352495A (zh) * 2015-11-17 2016-02-24 天津大学 加速度与光流传感器数据融合无人机水平速度控制方法
CN106586026A (zh) * 2016-12-06 2017-04-26 中国航空工业集团公司洛阳电光设备研究所 一种飞机相对跑道侧向偏差速率的测量方法
CN106774376A (zh) * 2017-01-25 2017-05-31 上海拓攻机器人有限公司 一种无人机仿地飞行控制方法及系统

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110488864A (zh) * 2019-08-15 2019-11-22 中国商用飞机有限责任公司 用于修正飞机的飞行控制系统中的loc信号的方法及系统
CN110488864B (zh) * 2019-08-15 2021-12-03 中国商用飞机有限责任公司 用于修正飞机的飞行控制系统中的loc信号的方法及系统
CN112230587A (zh) * 2020-10-29 2021-01-15 中国民用航空总局第二研究所 一种仪表着陆系统同频异呼运行监测系统及方法
CN112230587B (zh) * 2020-10-29 2021-08-10 中国民用航空总局第二研究所 一种仪表着陆系统同频异呼运行监测系统及方法

Similar Documents

Publication Publication Date Title
CN105352495B (zh) 加速度与光流传感器数据融合无人机水平速度控制方法
CN106586026B (zh) 一种飞机相对跑道侧向偏差速率的测量方法
US9128109B1 (en) Method and system for detecting errors in indicated air speed
CN111123966B (zh) 基于机载近地告警系统判断飞行阶段的方法
CN108008434B (zh) 具有低频仪表着陆系统定位信标异常检测的飞行控制系统及使用方法
CN107240063A (zh) 一种面向移动平台的旋翼无人机自主起降方法
CN104808231B (zh) 基于gps与光流传感器数据融合的无人机定位方法
KR20150043392A (ko) 바람 벡터 추산
CA2978836C (en) Flight control system with synthetic intertial glideslope deviation and method of use
EP3315909B1 (en) Flight control system with synthetic inertial localizer deviation and method of use
CN101246078A (zh) 应用于低速无人机上总静压受感器位置误差的修正方法
CN107727877A (zh) 一种基于仪表着陆系统的地速测量方法
CN109725649A (zh) 一种基于气压计/imu/gps多传感器融合的旋翼无人机定高算法
CN112201083A (zh) 管理飞行器在目视进近跑道阶段的飞行的方法、电子系统及计算机程序
CN109798918A (zh) 一种下滑天线下滑角的测试方法
CN106372307B (zh) 一种基于气动模型的民机气流角估计方法
CN103808309A (zh) 一种无人机三维航空摄影测树方法
CN110375710A (zh) 飞行器高度的确定方法、装置和飞行器系统
CN106651220B (zh) 一种基于曲线拟合的飞行员操作质量评估方法
CN117553773A (zh) 一种仪表着陆系统下滑道扩展方法
TWI805141B (zh) 用於無人機的定位方法和設備
CN114778887A (zh) 一种基于改进三角矢量模型的无人机测风方法和装置
RU2570339C1 (ru) Способ определения координат центра масс самолета в полете и устройство для его осуществления
US11046451B1 (en) Angle of attack estimation system and method
CN113514049A (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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180223