CN103809192A - 一种gnss接收机的动态校正算法 - Google Patents
一种gnss接收机的动态校正算法 Download PDFInfo
- Publication number
- CN103809192A CN103809192A CN201410064414.0A CN201410064414A CN103809192A CN 103809192 A CN103809192 A CN 103809192A CN 201410064414 A CN201410064414 A CN 201410064414A CN 103809192 A CN103809192 A CN 103809192A
- Authority
- CN
- China
- Prior art keywords
- gps signal
- frequency
- output matrix
- dbzp
- natural number
- 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
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/246—Acquisition or tracking or demodulation of signals transmitted by the system involving long acquisition integration times, extended snapshots of signals or methods specifically directed towards weak signal acquisition
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种GNSS接收机的动态校正算法,其采用差分相干技术,无需估计导航电文翻转,大大降低了运算开销;在相干累加之前加入动态校正环节:设计DBZP差分相干捕获算法输出矩阵,通过行DFT变换,提取与距离走动和距离弯曲相关的项,然后,设计Keystone变换算法去掉距离走动相关的项,然后依据相关函数在频域的相位信息,构造相位补偿函数,补偿距离弯曲,进一步削弱了多普勒频率及其变化率的影响,使得信号相干累加后的能量更加集中;而且,动态校正过程中没有用到中间估计量,可适用于较低信噪比的环境,最终,接收机可在低信噪比且具有较高动态的应用环境下保持正常工作。
Description
技术领域
本发明属于导航定位技术领域,具体涉及一种GNSS接收机的动态校正算法。
背景技术
卫星导航定位系统(GNSS)是一种以卫星为基础的无线电导航系统,能为陆、海、空的各类载体提供全天候、不间断、高精度、实时导航定位服务,已经应用于国民经济与日常生活的各个领域,如地面交通监管、飞机与船舶导航、精密受时、大地测量等。目前,全球范围内研发最早、应用最早的卫星定位系统GPS系统在我国已得到广泛应用,我国也参与了将于近两年建成的Galileo(伽利略)系统的建设,并正在自主研发全球卫星定位系统Compass(北斗二代),该系统2012年底已经在我国及其周边地区提供定位服务。因此,研究卫星导航及其接收机技术必将成为国内未来一段时间内的研究重点。
GNSS信号到达地面接收机时已相当微弱,比接收机内部热噪声低20~30dB,而且,多数接收机应用时一般处于运动状态。因此,要提高弱信号环境下捕获灵敏度,主要方法是增加相干累加时间,如全比特算法以及David M Lin等在2000年提出的二倍分组块补零(Double Block Zero Padding,DBZP)的GPS信号高灵敏度捕获算法,累加时间超过了导航电文比特长度的限制;但高动态会引起相干累加损耗,累加时间较长时,必须考虑运动的影响:当存在速度及加速度时,伪码相关峰会随着累加时间的变化而移动,形成所谓的距离走动和弯曲;此外利用快速傅氏变换做循环卷积时存在相关功率损失,当接收机和GNSS卫星相对速度较大时,接收机接收的GNSS信号将产生较大的多普勒频移,当多普勒频移较大且积累时间较长时,将造成码片速率发生较大变化,对码周期产生较大影响,从而会引起相干累加损耗,造成伪码相关峰的包络展宽、峰值降低。因此,当累加时间较长时,必须考虑动态的影响,需进行动态校正以克服这一问题。
针对动态校正算法,目前已经出现许多研究成果,如互相关法、谱峰跟踪法以及Yang J G、Huang X T等在2011年提出的使用Keystone变换校正距离弯曲的方法,但这些方法都只能用于高信噪比场合,在低信噪比情况下效果不好。
发明内容
针对现有技术所存在的上述技术问题,本发明提供了一种GNSS接收机的动态校正算法,能够克服弱信号环境下接收机运动对GNSS信号捕获的影响,补偿运动引起的相干累加损耗、伪码相关峰的包络展宽和峰值降低值,提高弱信号高动态应用环境下GNSS接收机捕获和跟踪电路的性能;当信噪比较低、并且接收机处于运动状态时,使用本发明动态校正算法的GNSS接收机也可稳定的给出定位结果。
本发明是在DBZP算法的基础上进行改进,在相干累加之前加入动态校正环节。
一种GNSS接收机的动态校正算法,包括如下步骤:
(1)对DBZP输出矩阵进行相关运算,得到相关器输出矩阵ψn;
(2)对相关器输出矩阵ψn做行DFT变换(离散傅里叶变换),得到对应的功率谱Ps(l,k);
(5)使相位δk替换相位θk代入功率谱Ps(l,k)中,并对新的功率谱进行DFT反变换,得到动态校正后的相关器输出矩阵进而对该相关器输出矩阵做列FFT变换(快速傅里叶变换),即可得到相干累加输出。
所述的相关器输出矩阵ψn的表达式如下:
其中:A和分别为GPS信号的幅值和初始相位,L为DBZP输出矩阵的行数,Sa()抽样函数,fd和fa分别为GPS信号的多普勒频率和多普勒频率变化率,N为DBZP输出矩阵列数,Ts为GPS信号的采样间隔,ξ=fd/f0,ζ=fa/f0,f0为GPS信号的载波频率,j为虚数单位,R()为DBZP输出矩阵的自相关函数,Δτ为伪码相位时延的估计误差,Tk=kNTs,Tn=nTs,k为自然数且0≤k≤L-1,n为自然数且0≤n≤N-1。
所述的功率谱Ps(l,k)的表达式如下:
其中:Pc(l)为伪码的功率谱密度函数,l为频点,且i为自然数,j为虚数单位,τ为伪码相位时延,fl为距离相关频率且fl=fs*l/N,fs为GPS信号的采样频率,N为DBZP输出矩阵列数,ξ=fd/f0,ζ=fa/f0,fd和fa分别为GPS信号的多普勒频率和多普勒频率变化率,f0为GPS信号的载波频率,Tk=kNTs,Ts为GPS信号的采样间隔,k为自然数且0≤k≤L-1,L为DBZP输出矩阵的行数。
所述的相位的表达式如下:
其中:fl为距离相关频率且fl=fs*l/N,l为频点,且i为自然数,fs为GPS信号的采样频率,N为DBZP输出矩阵列数,fd和fa分别为GPS信号的多普勒频率和多普勒频率变化率,f0为GPS信号的载波频率,Tk=kNTs,Ts为GPS信号的采样间隔,k为自然数且0≤k≤L-1,L为DBZP输出矩阵的行数。
其中:fl为距离相关频率且fl=fs*l/N,l为频点,且i为自然数,fs为GPS信号的采样频率,N为DBZP输出矩阵列数,fa为GPS信号的多普勒频率变化率,f0为GPS信号的载波频率,Tk=kNTs,Ts为GPS信号的采样间隔,k为自然数且0≤k≤L-1,L为DBZP输出矩阵的行数。
其中:fl为距离相关频率且fl=fs*l/N,l为频点,且i为自然数,fs为GPS信号的采样频率,N为DBZP输出矩阵列数,fd为GPS信号的多普勒频率,f0为GPS信号的载波频率,Tk=kNTs,Ts为GPS信号的采样间隔,k为自然数且0≤k≤L-1,L为DBZP输出矩阵的行数。
本发明的有益效果如下:
(1)本发明可以在低信噪比环境且存在高动态运动时保持正常工作,通过采用动态校正先进技术,提高GNSS定位的可靠性和定位精度等性能。
(2)DBZP捕获算法在导航电文跃变前将GNSS输入信号分成多段,有助于降低捕获过程中由大多普勒频移引起码片速率变化而造成的相关功率损失;然而在捕获过程中,在每个预检测积分时间内都需要对最可靠的数据位组合进行估计,并利用它去掉先前的数据位,这样需要较大的运算开销。为此,本发明采用差分相干技术,无需估计导航电文翻转,大大降低了运算开销。
(3)本发明利用DBZP捕获算法输出矩阵,设计算法估计载波多普勒频差并做补偿,进一步削弱了多普勒频率的影响,提高了捕获算法的精度。
(4)高动态会引起相干累加损耗,累加时间较长时,必须考虑运动的影响:当存在速度及加速度时,伪码相关峰会随着累加时间的变化而移动,形成所谓的距离走动和弯曲,此外,加速度还会带来多普勒扩展。可见,如相干累加时间较长,需进行动态校正克服运动会造成伪码相关峰的包络展宽、峰值降低等问题。故本发明在DBZP算法的基础上进行改进,在相干累加之前加入动态校正环节:首先采用Keystone变换校正多普勒频移,然后依据相关函数在频域的相位信息,构造相位补偿函数,补偿距离弯曲,从而使得信号相干累加后的能量更加集中,且动态校正过程中没有用到中间估计量,因此适用于低信噪比情况。
附图说明
图1为本发明GNSS接收机的结构示意图。
图2为输入信号与本地码信号的分块示意图。
图3为DBZP输出矩阵的示意图。
图4为分块扩组合与相关计算的示意图。
图5为本发明基于动态校正DBZP捕获算法的流程示意图。
具体实施方式
为了更为具体地描述本发明,下面结合附图及具体实施方式对本发明的技术方案进行详细说明。
本实施方式采用的GNSS接收机组成结构如图1所示,GNSS接收机由射频前端、基带信号处理模块、导航定位解算模块构成。基带信号处理模块包括信号捕获、跟踪、译码与导航电文提取等子模块组成。
GNSS信号捕获粗略估计伪码相位和载波多普勒频率,信号跟踪模块实现对伪码相位和载波多普勒频率的精确估计,以便实现GNSS信号的解扩和解调。译码与电文提取模块通过维特比译码获取导航电文,获取当前时刻下的卫星星历信息和伪距测量信息,导航定位解算模块用星历信息和测得的伪距和伪距率信息,实现导航定位解算。。
可见,捕获和跟踪部分实现信号的解扩、解调,其效果将直接影响接收机定位性能。本发明主要针对捕获算法部分的改进,以提高接收机的捕获速度和灵敏度。
接收机采用后的中频(IF)信号经数字下变频为基带信号后,采用DBZP算法进行捕获,DBZP捕获算法的步骤流程如图5所示,设预检测积分时间为TI毫秒,算法具体步骤如下:
步骤1:将TI毫秒输入信号分成Nb个长度为Ns的子块,同时将TI毫秒的本地C/A码分成Nb个长度为Ns的子块,如图2所示。
设相干积分时间为TI,信号采样间隔为Ts,GNSS接收信号下变频和采样后得到数字中频信号,则第k个相干时间段的数字中频信号表示为:
式中:fd为中频信号标称频率,fd和fa为多普勒频率和多普勒变化率,τ为伪码相位时延,ξ=fd/f0为伪码速率偏移量,ζ=fa/f0为伪码速率偏移量变化率,A为信号幅度。C()为伪随机码,伪码速率为fc。
第k个时间段本地产生的伪码信号模型为:
步骤2:将输入信号相邻的两个子块组合成一个长度为2×Ns的双块,最后一个分组块与大小为Ns的下一TI毫秒采样数据的第一采样块进行组合。在每个本地的C/A码子块后用零元素拓展成一个长度为2×Ns的双块。见图4所示。
对式(1)、(2)中信号进行双块零拓展得:
步骤3:将对应的双块用FFT循环卷积做相关运算,并将相关结果的第一个子块保存,即保存最开始的Ns个采样点,如图4所示。所保存的采样点将排列成一个大小为Nb×Ns的DBZP输出矩阵Mc,矩阵中每个列包含位于每个分组块中的下标相同的采样点。
将式(3)、(4)中的对应块利用FFT循环卷积进行相关运算,有:
式中:为下的不同码时延的相关结果,Rc()为C/A码循环卷积相关值,相关损耗为Sa(πfdNTs),为无关相位项。Tk=kNTs,Tn=nTs,k为自然数且0≤k≤L-1,n为自然数且0≤n≤N-1。
将式(5)中的第一个有用信息子块保存,即:
步骤4:为覆盖所有不确定的码时延,重复步骤2、3共L次(L为1ms码周期被分割的块数)。在每一次循环过程中,将复制码的分组块右循环移位一块,见图4所示。每次循环后,同上所述将所循环卷积结果最开始的Ns个采样点按顺序添加到矩阵Mc中,最后形成大小为Nb×(Ns×L)的矩阵Mc,如图3所示。
步骤5:对上述DBZP输出矩阵进行相干累加,得到:
步骤6:校正距离走动。
对式(7)中求和号内的信号进行离散傅里叶变换:
式中,l=-N/2,-N/2+1,",N/2,fl=fsl/N为距离相关频率,Pc(l)为伪码的功率谱密度函数。提取上式中最后一项的相位如下:
等号右边的两项包含有距离频率fl和列时间Tk的混合项,分别对应距离走动和距离弯曲。采用Keystone变换解除fl与Tk的耦合以消除距离走动。令k=f0k'/(f0+fl),k'=0,1,2,",L-1,带入式(9)后经推导得:
步骤7:校正距离弯曲和多普勒扩展。
可见,虽然距离走动已经校正,但距离弯曲和多普勒扩展仍然存在。为此,构造相位补偿函数:
式中:fa为多普勒变化率补偿值,且用k代替式k',式(10)减式(11),可得:
式中:△fa=fa-fa'为残留多普勒变化率,由泰勒展开可得f0/(f0+f1)≈(f0-f1)/f0,带入上式可得:
再将(11b)式带入(8)式得:
可见,距离走动已经消除,尽管距离弯曲和多普勒扩展仍然存在,但影响已经大大减小了。
步骤8:差分累加。对式(10)中保存的信息进行K次差分累加得:
利用式(13)中估计的多普勒频率偏差,校正差分相干捕获得到的载波频率,修正后的多普勒频率为:
修正后的结果进减小了残余多普勒频率的影响,提高了捕获频率精度。
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410064414.0A CN103809192A (zh) | 2014-02-25 | 2014-02-25 | 一种gnss接收机的动态校正算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410064414.0A CN103809192A (zh) | 2014-02-25 | 2014-02-25 | 一种gnss接收机的动态校正算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103809192A true CN103809192A (zh) | 2014-05-21 |
Family
ID=50706243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410064414.0A Pending CN103809192A (zh) | 2014-02-25 | 2014-02-25 | 一种gnss接收机的动态校正算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103809192A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107807374A (zh) * | 2017-10-26 | 2018-03-16 | 慧众行知科技(北京)有限公司 | 一种时变参数估计方法及系统 |
CN109471101A (zh) * | 2018-10-19 | 2019-03-15 | 湖北航天技术研究院总体设计所 | 一种弹载sar图像校正方法 |
CN109541531A (zh) * | 2018-11-22 | 2019-03-29 | 中电科仪器仪表有限公司 | 一种降低频差对时间延迟估计影响的方法 |
CN111337955A (zh) * | 2020-03-04 | 2020-06-26 | 中国人民解放军军事科学院国防科技创新研究院 | 星载单星测频无源定位接收机试验方法 |
CN112953517A (zh) * | 2021-02-24 | 2021-06-11 | 中国电子科技集团公司第十研究所 | 一种动态多普勒频偏校准方法 |
CN113009522A (zh) * | 2021-02-08 | 2021-06-22 | 浙江理工大学 | 多普勒频率残差校正的长时间相干积分捕获算法模块 |
CN113009523A (zh) * | 2021-02-22 | 2021-06-22 | 浙江理工大学 | 长时间相干积分捕获的多普勒频率估计与补偿方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080143598A1 (en) * | 2006-09-27 | 2008-06-19 | Purdue Research Foundation | Methods for improving computational efficiency in a global positioning satellite receiver |
CN101408616A (zh) * | 2008-11-24 | 2009-04-15 | 江南大学 | 适用于低信噪比数据的逆合成孔径雷达成像距离对齐方法 |
US20110037652A1 (en) * | 2006-11-08 | 2011-02-17 | Stafford Stephen J | Cross-PRN Phase Correction for GPS Processing |
CN102680988A (zh) * | 2012-05-23 | 2012-09-19 | 桂林电子科技大学 | 高灵敏度导航卫星信号非线性捕获方法及装置 |
CN102901956A (zh) * | 2012-09-27 | 2013-01-30 | 电子科技大学 | 一种雷达微弱目标检测方法 |
-
2014
- 2014-02-25 CN CN201410064414.0A patent/CN103809192A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080143598A1 (en) * | 2006-09-27 | 2008-06-19 | Purdue Research Foundation | Methods for improving computational efficiency in a global positioning satellite receiver |
US20110037652A1 (en) * | 2006-11-08 | 2011-02-17 | Stafford Stephen J | Cross-PRN Phase Correction for GPS Processing |
CN101408616A (zh) * | 2008-11-24 | 2009-04-15 | 江南大学 | 适用于低信噪比数据的逆合成孔径雷达成像距离对齐方法 |
CN102680988A (zh) * | 2012-05-23 | 2012-09-19 | 桂林电子科技大学 | 高灵敏度导航卫星信号非线性捕获方法及装置 |
CN102901956A (zh) * | 2012-09-27 | 2013-01-30 | 电子科技大学 | 一种雷达微弱目标检测方法 |
Non-Patent Citations (3)
Title |
---|
罗海坤: "一种基于运动补偿的高动态微弱直扩信号捕获算法", 《宇航学报》, vol. 34, no. 4, 30 April 2013 (2013-04-30), pages 541 * |
罗海坤等: "一种基于运动补偿的高动态微弱直扩信号捕获算法", 《宇航学报》, vol. 34, no. 4, 30 April 2013 (2013-04-30), pages 539 - 545 * |
莫建文等: "基于DBZP差分相干的GPS信号捕获算法", 《系统工程与电子技术》, vol. 34, no. 6, 30 June 2012 (2012-06-30), pages 1085 - 1089 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107807374A (zh) * | 2017-10-26 | 2018-03-16 | 慧众行知科技(北京)有限公司 | 一种时变参数估计方法及系统 |
CN109471101A (zh) * | 2018-10-19 | 2019-03-15 | 湖北航天技术研究院总体设计所 | 一种弹载sar图像校正方法 |
CN109541531A (zh) * | 2018-11-22 | 2019-03-29 | 中电科仪器仪表有限公司 | 一种降低频差对时间延迟估计影响的方法 |
CN111337955A (zh) * | 2020-03-04 | 2020-06-26 | 中国人民解放军军事科学院国防科技创新研究院 | 星载单星测频无源定位接收机试验方法 |
CN113009522A (zh) * | 2021-02-08 | 2021-06-22 | 浙江理工大学 | 多普勒频率残差校正的长时间相干积分捕获算法模块 |
CN113009522B (zh) * | 2021-02-08 | 2023-03-24 | 浙江理工大学 | 多普勒频率残差校正的长时间相干积分捕获算法模块 |
CN113009523A (zh) * | 2021-02-22 | 2021-06-22 | 浙江理工大学 | 长时间相干积分捕获的多普勒频率估计与补偿方法及系统 |
CN113009523B (zh) * | 2021-02-22 | 2023-09-29 | 浙江理工大学 | 长时间相干积分捕获的多普勒频率估计与补偿方法及系统 |
CN112953517A (zh) * | 2021-02-24 | 2021-06-11 | 中国电子科技集团公司第十研究所 | 一种动态多普勒频偏校准方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103809192A (zh) | 一种gnss接收机的动态校正算法 | |
CN103033828B (zh) | 一种高灵敏度北斗辅助授时装置和授时接收机及授时方法 | |
CN103645483B (zh) | 一种弱信号环境下北斗信号捕获方法 | |
CN104407323B (zh) | 一种高动态低信噪比扩频信号伪码时延测量方法 | |
CN102435999B (zh) | Gps接收机基带模块和gps信号捕获及跟踪方法 | |
CN102636798B (zh) | 基于环路状态自检测的sins/gps深组合导航方法 | |
CN107255822A (zh) | 多径环境下gnss接收机信号参数估计方法 | |
CN104049528B (zh) | 一种北斗授时方法及卫星导航接收机 | |
CN104483691B (zh) | 一种gnss组合精密单点定位方法 | |
CN104597471A (zh) | 面向时钟同步多天线gnss接收机的定向测姿方法 | |
CN113009524B (zh) | 长时间相干积分捕获的导航电文比特翻转估计方法及系统 | |
CN204101732U (zh) | 一种gnss新体制信号捕获装置 | |
CN107193023B (zh) | 一种具有闭式解的高精度北斗卫星系统单点定位方法 | |
CN104536016A (zh) | 一种gnss新体制信号捕获装置及方法 | |
CN102033236A (zh) | 一种卫星导航位置速度联合估计方法 | |
CN102890280A (zh) | 一种多模gnss组合接收机兼容捕获方法 | |
CN105425258A (zh) | 一种惯导系统辅助的高动态微弱信号gps捕获方法 | |
CN103543456A (zh) | 一种基于分段相关结合fft运算的大频偏gnss信号捕获方法 | |
CN113009522B (zh) | 多普勒频率残差校正的长时间相干积分捕获算法模块 | |
Wang et al. | A scheme for weak GPS signal acquisition aided by SINS information | |
CN104215979A (zh) | 一种基于分段相关结合fft运算的大频偏gnss信号捕获方法 | |
CN103645484A (zh) | 一种复合二进制偏移载波调制信号的处理方法 | |
Chen et al. | High dynamic GPS signal tracking based on UKF and carrier aiding technology | |
JP4777353B2 (ja) | Gps測位方法及びgps測位装置 | |
CN106154300B (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140521 |