CN103344246B - 火星动力下降段减弱动力学系统误差的两步滤波方法 - Google Patents
火星动力下降段减弱动力学系统误差的两步滤波方法 Download PDFInfo
- Publication number
- CN103344246B CN103344246B CN201310287572.8A CN201310287572A CN103344246B CN 103344246 B CN103344246 B CN 103344246B CN 201310287572 A CN201310287572 A CN 201310287572A CN 103344246 B CN103344246 B CN 103344246B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- msubsup
- overbar
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000001914 filtration Methods 0.000 title claims abstract description 46
- 238000005183 dynamical system Methods 0.000 title claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims description 45
- 230000001133 acceleration Effects 0.000 claims description 31
- 238000005259 measurement Methods 0.000 claims description 30
- 230000005484 gravity Effects 0.000 claims description 14
- 230000003313 weakening effect Effects 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000007704 transition Effects 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 241000651994 Curio Species 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Landscapes
- Indication And Recording Devices For Special Purposes And Tariff Metering Devices (AREA)
- Navigation (AREA)
Abstract
一种火星动力下降段减弱动力学系统误差的两步滤波方法,它包括以下步骤:一、建立工程实际方程;二、给定初始值;三、状态量滤波;四、动力学系统偏差滤波;五、更新相关系数、校正状态估计和动力学偏差估计;六、令k=k+1,返回步骤三往下进行,直到k等于火星动力下降段时间截止对应的时刻T时,即火星着陆器着陆为止,至此完成火星动力下降段减弱动力学系统误差的两步滤波方法。本方法统筹考虑了火星动力下降过程中,非线性、非高斯随机系统在动力学系统误差条件下的航天器位置速度估计问题。通过该方法在计算过程中引入对动力学系统误差进行估计和补偿,减弱了动力学系统误差对滤波引起的导航误差,有效保证航天器在火星动力下降段的位置速度估计。
Description
技术领域
本发明涉及火星动力下降段减弱动力学系统误差的两步滤波方法。属于航天导航技术领域。
背景技术
Kalman方法是非常常见的一种确定航天器位置速度方法。它要求动力学系统和量测系统的精确已知。在实际工程中动力学系统难精确得到,存在一定的不确定性和动力学系统误差。在火星动力下降段中,由于飞行器动力学系统采用加速度计和陀螺仪量(IMU)的惯性导航进行航迹递推。而加速度计测量的为非惯性力,即火星重力加速度无法测得。通常只能利用事先估计的重力带入计算,这样将产生较大的误差。如2012年美国的好奇号在其着陆点附近存在着一定的重力加速度误差。未来火星任务将在一些具有科学意义的地点着陆,这些地点的地貌复杂,其重力加速度将在火星着陆器着陆前更加难以精确测定,存在的重力加速度误差将影响着高精度的着陆。另外,由于初始安装在着陆器上的测量仪器存在一定的安装误差,这些误差将使得动力学系统和量测系统存在一定的未知系统误差。过大的误差会导致位置速度误差的增大甚至发散,引起航天器导航误差、降低导航精度。为了方便描述,下文以火星重力加速度为例,对火星重力加速度进行有效的估计和校正从而进一步提高导航精度。
现有技术中,可以用于确定航天器位置速度的方法有多种。
现有技术一,基于泰勒展开的扩展Kalman滤波估计方法。该方法忽略了动力学系统中的误差,如火星重力加速度误差。该方法给出了非线性动力学方程和非线性测量方程的泰勒展开加权融合的估算公式。
现有技术二,基于sigma点集(为正态分布采样策略)的无迹Kalman滤波方法。先根据正态分布的均值和方差计算出sigma点集,并确定出各点的权值,再通过动力学方程计算出航天器的位置速度,然后通过量测方程得到的量测数据对航天器的位置速度进行调整修正。它同样不考虑火星重力加速度误差对动力学系统的影响。
现有技术一适用于动力学系统和量测系统精确可知条件下。在火星动力下降过程中,火星重力加速度具有一定的不可估计的误差,因此不太适用于火星动力下降段。
现有技术二在测量手段有限,测量数据少,难以动力学系统误差估算和减弱。因此动力学系统中的误差将影响航天器的位置速度估计,难以进行调整修正,因此不太适用于火星动力下降段。
现有技术一二对于未知的火星环境,其动力学系统存在误差,如火星重力加速度误差。即使事先对着陆地点的重力加速度进行估计也还存在一定的误差。特别是在一些地形地貌复杂的具有科学探测的着陆点附近,火星重力加速度将难以进行事先估计校正。因此不太适用于火星动力下降段。
发明内容
1、目的:本发明的目的是提供一种火星动力下降段减弱动力学系统误差的两步滤波方法,以减小航天器位置速度误差,提高其精度。
2、技术方案:本发明的目的是通过以下技术方案来实现的。
本发明一种火星动力下降段减弱动力学系统误差的两步滤波方法,它包括以下步骤:
步骤一、建立工程实际方程:离散的动力学系统和量测系统
zk=h(xk)+vk (3)
其中xk表示系统状态量,zk是测量系统测量值,bgk是未知的动力学系统误差(即火星重力加速度误差)。非线性方程f(·)和h(·)分别是状态转移方程和量测方程并且关于可xk微。矩阵Bk为未知的动力学系统误差(即火星重力加速度误差)对动力系统的驱动矩阵,矩阵Ak为未知的动力学系统误差(即火星重力加速度误差)转移矩阵,矩阵Bk,Ak具有恰当的维数。是动力学系统误差的噪声,和vk分别是动力学系统噪声和量测噪声,它们是不相关的高斯白噪声满足以下式子。
步骤二、给定初始值:
为初始状态的估计值,为初始状态估计均方误差,为初始误差的估计值通常取0,为初始误差的估计均方误差,为初始状态和初始误差的估计均方误差。将这些以上初始值通过以下式子计算得到滤波初始条件。
其中,为滤波初始状态的估计值,为滤波初始状态估计均方误差,为滤波初始误差的估计值通常取0,为滤波初始误差的估计均方误差阵,V0为滤波初始状态和初始误差的相关系数。
步骤三、状态量滤波
其中
式中:为tk-1时刻的状态量,uk-1为tk-1时刻的控制输入量,为tk-1时刻的滤波校正输入量。为状态的滤波一步预测。为tk-1时刻的滤波状态估计均方误差,Φk-1为tk-1时刻到tk时刻的一步转移矩阵;为系统的噪声的滤波校正方差阵,为一步预测均方误差。Hk为量测阵,为状态增益。I为单位阵,为tk时刻的滤波状态估计均方误差。为量测新息,为滤波状态估计。为状态的一步预测由式(17)计算得到。
步骤四、动力学系统偏差滤波
其中步骤三步骤四中的对应式子为
式中,Ak-1为tk-1时刻动力学系统误差的转移矩阵,为tk-1时刻动力学系统误差(火星重力加速度误差)的滤波估计,为动力学系统误差(火星重力加速度误差)的一步预测,为tk-1时刻动力学系统误差的估计均方误差阵,为动力学系统误差的一步预测均方误差,为动力学系统误差的噪声的方差阵,为动力学系统误差状态增益,Rk为量测噪声的方差阵。为tk时刻动力学系统误差的估计均方误差阵,为动力学系统误差(火星重力加速度误差)的量测新息,为动力学系统误差的状态估计。Nk为tk时刻动力学系统偏差对量测系统校正量测阵。为系统噪声方差阵。Uk,Vk为状态估计均方误差阵U-V分解的对应阵。
步骤五、更新相关系数、校正状态估计和动力学偏差估计:情况如下
为tk时刻校正后的状态量,为tk时刻校正后的状态估计均方误差,为tk时刻校正后的状态量,为tk时刻校正后的状态估计均方误差。
步骤六、令k=k+1,返回步骤三往下进行。直到k等于火星动力下降段时间截止对应的时刻T时,即火星着陆器着陆为止。至此完成火星动力下降段减弱动力学系统误差的两步滤波方法。
其中,步骤一中要根据实际情况合理建立相应的矩阵Bk,Ak。
其中,步骤二中根据实际情况估计初值,为火星动力下降段之前的伞降段末端状态估计值和估计均方误差。且动力学系统偏差和测量系统中的未知测量系统误差为不相关量。
其中,在步骤一中所述的建立工程实际方程,其步骤如下:
a、分析动力学不确定性之间的相互关系,并进行相应的数值计算分析;
b、得到这些不确定性因素对动力学模型的影响主要是通过哪些量对动力学模型产生影响从而引起误差的传播;
c、将动力学系统xk+1=f(xk,uk)改写为考虑动力学系统偏差的动力学系统
3、优点和功效:
本发明统筹考虑了火星动力下降过程中,非线性、非高斯随机系统在动力学系统误差条件下的航天器位置速度估计问题。通过火星动力下降段减弱火星重力加速度误差的两步滤波方法在计算过程中引入了对动力学系统误差进行估计和补偿,减弱了动力学系统误差对滤波引起的导航误差。因而本发明提出的算法可以有效保证航天器在火星动力下降段的位置速度估计。
附图说明
图1为三轴位移误差在在本文方法和传统EKF方法下的计算结果
图2为三轴速度误差在在本文方法和传统EKF方法下的计算结果
图3上图为高度方向位移误差,下图为高度方向速度误差(利用本文方法)
图4上图为火星重力加速偏差值和估计值,下图为火星重力加速度估计误差
图5为本发明所述方法流程图
图中的代号、符号说明如下:
TEKF为火星动力下降段减弱动力学系统误差的两步滤波方法
EKF为火星动力下降段的扩展Kalman滤波方法
具体实施方式
本发明涉及火星动力下降段减弱动力学系统误差的两步滤波方法,其航天器在火星动力下降过程中,其依靠IMU的火星着陆动力学系统如下方程:
其中为动力学系统的状态量,r=[r1,r2,r3]T为相对火星着陆点坐标系下的三轴位移,v=[v1,v2,v3]T为火星着陆器相对火星着陆点坐标系下的三轴速度,e=[φ,θ,ψ]T为三轴Euler姿态角,φ,θ,ψ分别为滚转角、俯仰角、航向角。为加速度计测量值,为陀螺仪测量值,ba为加速度计偏差,bω为陀螺仪偏差,g为火星重力加速度,bg为火星重力加速度偏差。ξa为加速度计测量白噪声,ξω为陀螺仪测量白噪声。为测量仪器坐标系到火星着陆点坐标系的转换矩阵。K为姿态运动学矩阵。矩阵Bk=[01×5 1 01×9]T,矩阵Ak=1。
目前火星动力下降段,目前美国喷气推力实验室(JPL)开发了一种集高度计和速度计为一体的小体积、低能耗的微型传感器(MCAV),它采用激光测距测速来实现高度、速度信息的同时输出。其对应的测量系统方程可以如下式来表示:
其中,r3为高度方向,为火星着陆点坐标系到测量仪器坐标系的转换矩阵,v=[v1,v2,v3]T为火星着陆器相对火星着陆点坐标系下的三轴速度,即为测量仪器测得速度数据。
本发明一种火星动力下降段减弱动力学系统误差的两步滤波方法,见图5所示,它包括以下步骤:
步骤一:建立工程实际方程:火星动力下降段对应的离散动力学系统可以改写为如下形式:
xk+1=f(xk,wk)+Bkbgk (22)
其中x为火星动力下降段动力学系统中的状态量具体包含式(19)中的左边的各分量,即矩阵Bk=[01×5 1 01×9]T。其对应的离散量测方程为
zk=h(xk)+vk (23)
其中
式中,测量高度h=r3,vb=[vx vy vz]T是加速度计在体坐标系下的测量值,可以通过式(21)进行换算。vk为量测噪声。
步骤二、给定初始值:
初始值为飞行器在伞降段末端状态估计得到,如表一
表一 火星动力下降段的真实值及初始估计值
其中本文中火星动力下降段的真实值及初始估计值为仿真选取值。也可以是另外两组不相同的数值。初始状态估计均方误差 初始状态和初始误差的估计均方误差初始误差的估计均方误差初始误差的估计
步骤三、状态量滤波:
按照公式(4)-(10)进行滤波,对火星动力下降段的状态进行估计。其中动力学系统噪声方差阵为 量测噪声方差阵为 动力学系统误差的噪声的方差阵tk-1时刻到tk时刻的一步转移矩阵Φk-1=I+Fk-1·Δt,tk时刻的量测阵 其中,I3×3为3阶单位矩阵, Δt为测量时间间隔。按照公式(16)校正后再代入量测阵Hk中。
步骤四、动力学系统偏差滤波:
按照式(11)-式(15)对动力学系统偏差进行滤波估计,其中根据实际动力学系统分析选取相应矩阵Ak-1=1。
步骤五、更新相关系数、校正状态估计和动力学偏差估计:
按照式(16)-式(18)对步骤三、四中的状态及其均误差进行校正。
步骤六、令k=k+1,返回步骤三往下进行。直到k等于火星动力下降段时间截止对应的时刻T时。至此完成火星动力下降段减弱动力学系统误差的两步滤波方法。
其中截止时间主要取决于飞行器与火星表面的高度及速度,能否满足软着陆的条件要求。
通过火星动力下降段减弱动力学系统误差的两步滤波方法得到着陆器三轴位移误差在在本文方法和传统EKF方法下的计算结果见图1,图2为三轴速度误差在在本文方法和传统EKF方法下的计算结果;图3上图为高度方向位移误差,下图为高度方向速度误差(利用本文方法)。图4上图为火星重力加速偏差值和估计值,下图为火星重力加速度估计误差。通过以上步骤,火星动力下降段减弱动力学系统误差的两步滤波方法能够成功消除动力学模型中的偏差,有效地提高导航的精度。
以上所述仅为本发明较佳的实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化和替换都应涵盖在本发明的保护范围之内,另外本发明提供的方法可以集成到火星动力下降的航天器位置速度估计软件中。
Claims (2)
1.一种火星动力下降段减弱动力学系统误差的两步滤波方法,其特征在于:它包括以下步骤:
步骤一、建立工程实际方程:离散的动力学系统和量测系统
zk=h(xk)+vk (3)
其中xk表示系统状态量,zk是测量系统测量值,bgk是未知的动力学系统误差即火星重力加速度误差;非线性方程f(·)和h(·)分别为状态转移方程和量测方程,且关于xk可导;矩阵Bk为未知的动力学系统误差即火星重力加速度误差对动力系统的驱动矩阵,矩阵Ak为未知的动力学系统误差即火星重力加速度误差转移矩阵,矩阵Bk,Ak具有恰当的维数;是动力学系统误差的噪声,和vk分别是动力学系统噪声和量测噪声,它们是不相关的高斯白噪声满足以下式子:
其中Rk≥0,δkj是δ函数,当k=j时δkj=1当k≠j时δkj=0.
步骤二、给定初始值:
为初始状态的估计值,为初始状态估计均方误差,为初始误差的估计值通常取0,为初始误差的估计均方误差,为初始状态和初始误差的估计均方误差;将这些以上初始值通过以下式子计算得到滤波初始条件:
其中,为滤波初始状态的估计值,为滤波初始状态估计均方误差,为滤波初始误差的估计值通常取0,为滤波初始误差的估计均方误差阵,V0为滤波初始状态和初始误差的相关系数;
步骤三、状态量滤波
其中
式中:为tk-1时刻的状态量,uk-1为tk-1时刻的控制输入量,为tk-1时刻的滤波校正输入量;为状态的滤波一步预测;为tk-1时刻的滤波状态估计均方误差,Φk-1为tk-1时刻到tk时刻的一步转移矩阵;为系统的噪声的滤波校正方差阵,为一步预测均方误差;Hk为量测阵,为状态增益;I为单位阵,为tk时刻的滤波状态估计均方误差;为量测新息,为滤波状态估计;为状态的一步预测由式(17)计算得到;
步骤四、动力学系统偏差滤波
其中步骤三步骤四中的对应式子为
式中,Ak-1为tk-1时刻动力学系统误差的转移矩阵,为tk-1时刻动力学系统误差即火星重力加速度误差的滤波估计,为动力学系统误差即火星重力加速度误差的一步预测,为tk-1时刻动力学系统误差的估计均方误差阵,为动力学系统误差的一步预测均方误差,为动力学系统误差的噪声的方差阵,为动力学系统误差状态增益,Rk为量测噪声的方差阵;为tk时刻动力学系统误差的估计均方误差阵,为动力学系统误差即火星重力加速度误差的量测新息,为动力学系统误差的状态估计;Nk为tk时刻动力学系统偏差对量测系统校正量测阵;为系统噪声方差阵;Uk,Vk为状态估计均方误差阵U-V分解的对应阵;
步骤五、更新相关系数、校正状态估计和动力学偏差估计:情况如下
为tk时刻校正后的状态量,为tk时刻校正后的状态估计均方误差,为tk时刻校正后的状态量,为tk时刻校正后的状态估计均方误差;
步骤六、令k=k+1,返回步骤三往下进行,直到k等于火星动力下降段时间截止对应的时刻T时,即火星着陆器着陆为止;至此完成火星动力下降段减弱动力学系统误差的两步滤波方法。
2.根据权利要求1所述的一种火星动力下降段减弱动力学系统误差的两步滤波方法,其特征在于:在步骤一中所述的建立工程实际方程,其步骤如下:
a、分析动力学不确定性之间的相互关系,并进行相应的数值计算分析;
b、得到这些不确定性因素对动力学模型的影响是通过哪些量对动力学模型产生影响从而引起误差的传播;
c、将动力学系统xk+1=f(xk,uk)改写为考虑动力学系统偏差的动力学系统
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310287572.8A CN103344246B (zh) | 2013-07-10 | 2013-07-10 | 火星动力下降段减弱动力学系统误差的两步滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310287572.8A CN103344246B (zh) | 2013-07-10 | 2013-07-10 | 火星动力下降段减弱动力学系统误差的两步滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103344246A CN103344246A (zh) | 2013-10-09 |
CN103344246B true CN103344246B (zh) | 2015-08-26 |
Family
ID=49279062
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310287572.8A Expired - Fee Related CN103344246B (zh) | 2013-07-10 | 2013-07-10 | 火星动力下降段减弱动力学系统误差的两步滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103344246B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103884333B (zh) * | 2014-03-31 | 2017-03-15 | 北京控制工程研究所 | 一种深空探测自主导航初始基准捕获方法 |
CN105151327B (zh) * | 2015-09-11 | 2017-02-01 | 北京理工大学 | 一种多约束下着陆器转移能力确定方法 |
CN105947238B (zh) * | 2016-05-30 | 2017-05-10 | 北京航空航天大学 | 一种转动惯量不确定的火星着陆器抗干扰姿态控制方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101270993A (zh) * | 2007-12-12 | 2008-09-24 | 北京航空航天大学 | 一种远程高精度自主组合导航定位方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6882964B2 (en) * | 2002-03-06 | 2005-04-19 | California Institute Of Technology | High accuracy inertial sensors from inexpensive components |
-
2013
- 2013-07-10 CN CN201310287572.8A patent/CN103344246B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101270993A (zh) * | 2007-12-12 | 2008-09-24 | 北京航空航天大学 | 一种远程高精度自主组合导航定位方法 |
Non-Patent Citations (3)
Title |
---|
着陆器在火星大气进入段的自主导航方法研究;王轶博;《中国优秀硕士论文全文数据库》;20120515;全文 * |
自适应扩展增量Kalman滤波方法;傅惠民等;《航空动力学报》;20120831;第27卷(第8期);第1734-1737页 * |
自适应无迹增量滤波方法;傅惠民等;《航空动力学报》;20130228;第28卷(第2期);第259-263页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103344246A (zh) | 2013-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107525503B (zh) | 基于双天线gps和mimu组合的自适应级联卡尔曼滤波方法 | |
CN101846510B (zh) | 一种基于星敏感器和陀螺的高精度卫星姿态确定方法 | |
CN103245359B (zh) | 一种惯性导航系统中惯性传感器固定误差实时标定方法 | |
CN104764467B (zh) | 空天飞行器惯性传感器误差在线自适应标定方法 | |
CN102519470B (zh) | 多级嵌入式组合导航系统及导航方法 | |
CN101963513B (zh) | 消除水下运载体捷联惯导系统杆臂效应误差的对准方法 | |
CN106500693A (zh) | 一种基于自适应扩展卡尔曼滤波的ahrs算法 | |
CN113340298B (zh) | 一种惯导和双天线gnss外参标定方法 | |
CN102645223B (zh) | 一种基于比力观测的捷联惯导真空滤波修正方法 | |
CN105300387B (zh) | 一种火星大气进入段非线性非高斯秩滤波方法 | |
CN103344259A (zh) | 一种基于杆臂估计的ins/gps组合导航系统反馈校正方法 | |
CN103884340B (zh) | 一种深空探测定点软着陆过程的信息融合导航方法 | |
CN104215244B (zh) | 基于发射惯性坐标系的空天飞行器组合导航鲁棒滤波方法 | |
CN103363991A (zh) | 一种适应月面崎岖地形的imu与测距敏感器融合方法 | |
CN103344260A (zh) | 基于rbckf的捷联惯导系统大方位失准角初始对准方法 | |
CN102680000A (zh) | 应用零速/航向修正的光纤捷联惯组在线标定方法 | |
CN103218482B (zh) | 一种动力学系统中不确定参数的估计方法 | |
CN105157724A (zh) | 一种基于速度加姿态匹配的传递对准时间延迟估计与补偿方法 | |
CN103708045A (zh) | 一种探月飞船跳跃式再入的在线参数辨识方法 | |
CN106441291A (zh) | 一种基于强跟踪sdre滤波的组合导航系统及导航方法 | |
CN104359496A (zh) | 基于垂线偏差补偿的高精度姿态修正方法 | |
CN103344246B (zh) | 火星动力下降段减弱动力学系统误差的两步滤波方法 | |
CN111190207B (zh) | 基于pstcsdref算法的无人机ins bds组合导航方法 | |
CN105865432B (zh) | 一种针对陀螺仪多源噪声的混合滤波方法与测试平台 | |
CN103411627B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150826 Termination date: 20170710 |
|
CF01 | Termination of patent right due to non-payment of annual fee |