CN101825715A - 一种卫星导航接收设备波束形成的方法 - Google Patents
一种卫星导航接收设备波束形成的方法 Download PDFInfo
- Publication number
- CN101825715A CN101825715A CN 201010160077 CN201010160077A CN101825715A CN 101825715 A CN101825715 A CN 101825715A CN 201010160077 CN201010160077 CN 201010160077 CN 201010160077 A CN201010160077 A CN 201010160077A CN 101825715 A CN101825715 A CN 101825715A
- Authority
- CN
- China
- Prior art keywords
- array
- weights
- array element
- antenna
- group
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000015572 biosynthetic process Effects 0.000 claims description 29
- 238000012937 correction Methods 0.000 claims description 4
- 238000012423 maintenance Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种卫星导航接收设备的波束形成方法,它涉及基于阵列天线的卫星导航接收设备的波束形成的方法,所述方法包括:一、确定期望天线相位中心的相位值;二、在权值求解过程中加入保证天线相位中心恒定的约束条件;三、权值输出,在权值输出前应根据阵元的损坏情况对天线相位中心进行修正。通过本发明可以保证阵列天线相位中心的稳定,提高阵列天线相位中心的精度,并具有简单实用的优点。
Description
技术领域
本发明涉及卫星导航接收设备的波束形成方法,尤其涉及基于阵列天线的卫星导航接收设备的波束形成的方法。
背景技术
目前,阵列信号处理已经应用到了通信、雷达、声纳等系统中,具有处理灵活、抗干扰能力强等特点,波束形成方法是阵列信号处理中的一个重要技术。波束形成方法就是计算一组权值,阵列天线根据这组权值可以产生空间指向的波束,在设定的方向上形成主波束,接收到用户信号。
对于接收设备来说,波束形成的过程是把阵列天线接收到的信号经下变频变换成中频或基带信号后,进行A/D采样数字化,然后送到波束形成器中。在波束形成器里,波束合成单元根据波束形成单元计算出的波束形成权值,完成对各路信号的复加权处理,最后合成出所需的波束输出信号,其中波束形成单元所使用的波束形成方法是波束形成的核心。
目前,波束形成方法可以分为直接波束形成方法和自适应波束形成方法。直接波束形成方法就是根据所需的波束指向、旁瓣、零点位置等特性计算权值;自适应波束形成方法就是根据天线接收到的信号和一定的权值计算准则形成权值,常用的权值计算准则有最小均方误差(MMSE)准则、最大信号与干扰噪声比(SINR)准则、最小而乘(LS)准则等。
在卫星导航接收设备中,天线相位中心是其中很重要的一个指标参数,在接收卫星导航信号时,希望天线相位中心保持恒定,如果忽略天线相位中心对卫星导航定位有一定的影响。
现有的波束形成方法的主要问题是:在权值计算过程中,没有保证天线相位中心不变的约束条件,不能保证得到的权值不对天线相位中心造成影响。
发明内容
本发明的目的在于避免上述背景技术中的不足而提供一种应用于卫星导航接收设备的波束形成方法。本发明方法实现简单,能够避免权值影响天线相位中心,可以在此基础上形成相应的软件和硬件实现。
本发明一种卫星导航接收设备波束形成的方法,是把阵列天线接收到的信号经下变频变换成中频或基带信号后,进行A/D采样数字化,然后送到波束形成器中,在波束形成器里,波束合成单元根据波束形成单元计算出的波束形成权值,完成对各路信号的复加权处理并合路输出,其特征在于包括步骤:
①确定阵列天线相位中心的相位值为Φ;
②在权值求解过程中加入阵列天线相位中心约束条件,计算得到的一组权值;
③权值输出。
作为本发明的改进,所述步骤②中在权值求解过程中加入的阵列天线相位中心约束条件为:
计算得到的一组权值,使得整个波束指向空间的天线相位中心保持恒定为Φ;其中,N为阵元个数,wn是第n个阵元的权值;an(θ)是θ方向入射信号的方向向量在第n个阵元上的值;Θ是整个波束扫描空间;A表示增益;Φ为步骤①的相位值;j是虚数单位,
作为本发明的改进,步骤②中在权值求解过程中加入的阵列天线相位中心约束条件为:
计算得到一组权值,使得期望波束指向θ0上保持天线相位中心恒定为Φ;其中,θ0为期望的波束指向,N为阵元个数,wn是第n个阵元的权值;an(θ)是θ方向入射信号的方向向量在第n个阵元上的值;A表示增益;Φ为步骤①的相位值;j是虚数单位,此部分技术方案的改进,是针对在权值求解过程中加入天线相位中心约束条件为(θ∈Θ)时有可能造成无法求解权值的情况。当不能求解权值或者不需要在整个波束指向空间保持天线相位中心恒定时,可以采用一种简化的处理方式:首先求解出一组权值wn,然后对求解得到的权值采用约束条件(θ=θ0)修正,这样可保证可以得到一组权值,能够在期望波束指向θ0上保持天线相位中心恒定为Φ。
作为本发明的改进,当所述的阵列天线为偶数2M个阵元,并且布阵方式满足条件
a)阵列天线可以唯一的确定一个几何中心;
b)阵列天线中的所有阵元可以分为M组关于几何中心对称的阵元组,每个阵元组有两个阵元,两个阵元与几何中心在一条直线上,并且每个阵元到几何中心的距离相等;
所述步骤②中的天线相位中心约束条件可以进一步简化:对于M组关于几何中心对称的阵元组,每个阵元组的a、b两个阵元对应的权值为wa、wb,wa、wb满足
angle(wa)+angle(wb)=Φ
abs(wa)=abs(wb)
其中,angle(·)表示取一个复数的相位的运算,abs(·)表示取一个复数的模值的运算。
作为本发明的改进,当所述阵列天线为奇数2M+1个阵元,布阵方式满足条件:
a)阵列天线可以唯一的确定一个几何中心点;
b)阵列天线中有一个阵元位于几何中心,其他阵元可以分为M组关于几何中心对称的阵元组,每个阵元组有两个阵元,两个阵元与几何中心在一条直线上,并且每个阵元到几何中心的距离相等;
步骤②中的天线相位中心约束条件可以进一步简化:当阵元数为偶数2M+1时,位于几何中心的阵元o对应的权值为wo,wo满足:
angle(wo)=Φ
对于其余M组关于几何中心对称的阵元组,每个阵元组的a、b两个阵元对应的权值为wa、wb,wa、wb满足:
angle(wa)+angle(wb)=Φ
abs(wa)=abs(wb)
作为本发明的改进,步骤③中在权值输出时,需要根据阵元的损坏情况,修正权值,在期望波束指向θ0上天线相位中心保持恒定为Φ,则阵元权值为:
w′n=w′n·B·exp[j·β]
其中,w′n为完好阵元的权值,损坏的阵元权值将不用考虑;β为阵元损坏造成的在期望波束指向θ0上产生的天线相位中心偏移;B是对权值幅度修正的系数。
angle(·)表示取一个复数的相位的运算,∑wn是所有阵元权值的和,∑w′n是所有完好阵元权值的和。
本发明相比背景技术具有如下优点:
(i)本发明在权值求解过程中,加入了保持天线相位中心恒定的约束条件,具有提高天线相位中心的精度的优点。
(ii)在有阵元损坏时,本发明能够修正阵元权值无效带来的天线相位中心偏移,保证了天线相位中心的恒定,具有实用性强的特点。
(iii)在具有几何中心的阵列天线中,本发明提出了的简化的约束条件,能够简化波束求解过程,具有简单实用的特点。
附图说明
图1是本发明波束形成方法的流程图;
图2是两个阵元关于几何中心对称的示意图;
图3是根据本发明的实施例的波束形成方法的流程图;
图4实施例的阵列天线布阵形式。
具体实施方式
下面结合具体实施例和附图对本发明做进一步的描述:
本发明的波束形成方法的流程图如图1所示,这种波束形成方法具体包括步骤:
步骤101:确定阵列天线相位中心的相位值为Φ。
步骤102:在权值求解过程中加入天线相位中心约束条件
计算得到的一组权值,使得整个波束指向空间的天线相位中心保持恒定为Φ。其中,N为阵元个数,wn是第n个阵元的权值;an(θ)是θ方向入射信号的方向向量在第n个阵元上的值;Θ是整个波束扫描空间;A表示增益;Φ是步骤①的相位值;j是虚数单位,
在权值求解过程中加入天线相位中心约束条件,有可能造成无法求解权值的情况。当不能求解权值或者不需要在整个波束指向空间保持天线相位中心恒定时,可以采用一种简化的处理方式:首先求解出一组权值wn,然后对求解得到的权值采用约束条件
修正,这保证可以得到一组权值,能够在期望波束指向θ0上保持天线相位中心恒定为Φ。其中,θ0为期望的波束指向,N为阵元个数,wn是第n个阵元的权值;an(θ)是θ方向入射信号的方向向量在第n个阵元上的值;A表示增益;Φ为步骤①的相位值;j是虚数单位,
步骤103:权值输出。根据阵元的损坏情况,在权值输出时修正权值,使得在期望波束指向θ0上天线相位中心保持恒定为Φ,则阵元权值为:
w′n=w′n·B·exp[j·β]
其中,w′n为完好阵元的权值;β为阵元损坏造成的在期望波束指向θ0上产生的天线相位中心偏移;B是对权值幅度修正的系数。
angle(·)表示取一个复数的相位的运算,∑wn是所有阵元权值的和,∑w′n是所有完好阵元权值的和。
如图2所示,所述步骤102中,当阵列天线的布阵方式满足条件:
a)阵列天线可以唯一的确定一个几何中心;
b)当阵元数为偶数2M时,阵列天线中的所有阵元可以分为M组关于几何中心对称的阵元组,每个阵元组有两个阵元,两个阵元与几何中心在一条直线上,并且每个阵元到几何中心的距离相等;
c)当阵元数为奇数2M+1时,有一个阵元位于阵列天线的几何中心,其他阵元可以分为为M组关于几何中心对称的阵元组,每个阵元组有两个阵元,两个阵元与几何中心在一条直线上,并且每个阵元到几何中心的距离相等;
步骤102中的权值求解约束条件进一步简化为:
1)当阵元数为偶数2M时,对于M组关于几何中心对称的阵元组,每个阵元组的a、b两个阵元对应的权值为wa、wb,wa、wb满足:
angle(wa)+angle(wb)=Φ
abs(wa)=abs(wb)
其中,angle(·)表示取一个复数的相位的运算,abs(·)表示取一个复数的模值的运算。
2)当阵元数为偶数2M+1时,位于几何中心的阵元o对应的权值为wo,wo满足:
angle(wo)=Φ
对于其余M组关于几何中心对称的阵元组,每个阵元组的a、b两个阵元对应的权值为wa、wb,wa、wb满足:
angle(wa)+angle(wb)=Φ
abs(wa)=abs(wb)
图3是具体实施例的波束形成流程图,具体步骤包括
步骤301:确定期望的天线相位中心的相位值为Φ=0。
步骤302:具体实施例的权值求解采用均匀加权求解方式,求解过程中加入天线相位中心约束条件,保证整个波束指向空间的天线相位中心保持恒定为Φ:
具体实施例中阵列天线采用4阵元的线阵,如图4所示,布阵方式满足条件
a)阵列天线可以唯一的确定一个几何中心;
b)阵元天线有4个阵元,为偶数阵元,阵元间隔相等,阵列天线中的所有阵元可以分为2组关于几何中心对称的阵元组(阵元1、4和阵元2、3),每个阵元组有两个阵元,两个阵元与几何中心在一条直线上,并且每个阵元到几何中心的距离相等;
求解过程中加入的天线相位中心约束条件可以简化为:阵列天线中的4个阵元天线可以分为2组关于几何中心对称的阵元组,每个阵元组的a、b两个阵元对应的权值为wa、wb,wa、wb满足
angle(wa)+angle(wb)=Φ
abs(wa)=abs(wb)
其中,angle(·)表示取一个复数的相位的运算,abs(·)表示取一个复数的模值的运算。
实施例中能够得到满足天线相位中心约束条件的权值,求解得到的阵元权值为
步骤303:根据阵元的损坏情况,修正权值,在期望波束指向θ0上天线相位中心保持恒定为Φ,则阵元权值为:
w′n=w′n·B·exp[j·β]
其中,w′n为完好阵元的权值,损坏的阵元权值将不用考虑;β为阵元损坏造成的在期望波束指向θ0上产生的天线相位中心偏移;B是对阵元幅度修正的系数。
angle(·)表示取一个复数的相位的运算,∑wn是所有阵元权值的和,∑w′n是所有完好阵元权值的和。
步骤304:权值输出。
除上述实施例外,本发明还可以有其他实施方式。凡采用等同替换或等效变换形式的技术方案,均落在本发明要求的保护范围。
Claims (6)
1.一种卫星导航接收设备波束形成的方法,是把阵列天线接收到的信号经下变频变换成中频或基带信号后,进行A/D采样数字化,然后送到波束形成器中,在波束形成器里,波束合成单元根据波束形成单元计算出的波束形成权值,完成对各路信号的复加权处理并合路输出,其特征在于包括步骤:
①确定阵列天线相位中心的相位值为Φ;
②在权值求解过程中加入阵列天线相位中心约束条件,计算得到的一组权值;
③权值输出。
3.根据权利要求1所述的一种卫星导航接收设备波束形成的方法,其特征在于:步骤②中在权值求解过程中加入的阵列天线相位中心约束条件为:
计算得到一组权值,使得期望波束指向θ0上保持天线相位中心恒定为Φ;其中,θ0为期望的波束指向,N为阵元个数,wn是第n个阵元的权值;an(θ)是θ方向入射信号的方向向量在第n个阵元上的值;A表示增益;Φ为步骤①的相位值;j是虚数单位,
4.根据权利要求1所述的一种卫星导航接收设备波束形成的方法,其特征在于:所述的阵列天线为偶数2M个阵元,并且布阵方式满足条件
a)阵列天线可以唯一的确定一个几何中心;
b)阵列天线中的所有阵元可以分为M组关于几何中心对称的阵元组,每个阵元组有两个阵元,两个阵元与几何中心在一条直线上,并且每个阵元到几何中心的距离相等;
权利要求2所述步骤②中的天线相位中心约束条件可以进一步简化:对于M组关于几何中心对称的阵元组,每个阵元组的a、b两个阵元对应的权值为wa、wb,wa、wb满足
angle(wa)+angle(wb)=Φ
abs(wa)=abs(wb)
其中,angle(·)表示取一个复数的相位的运算,abs(·)表示取一个复数的模值的运算。
5.根据权利要求1所述的一种卫星导航接收设备波束形成的方法,其特征在于步骤②中,当阵列天线为奇数2M+1个阵元,布阵方式满足条件:
a)阵列天线可以唯一的确定一个几何中心点;
b)阵列天线中有一个阵元位于几何中心,其他阵元可以分为M组关于几何中心对称的阵元组,每个阵元组有两个阵元,两个阵元与几何中心在一条直线上,并且每个阵元到几何中心的距离相等;
权利要求2步骤②中的天线相位中心约束条件可以进一步简化:当阵元数为偶数2M+1时,位于几何中心的阵元o对应的权值为wo,wo满足:
angle(wo)=Φ
对于其余M组关于几何中心对称的阵元组,每个阵元组的a、b两个阵元对应的权值为wa、wb,wa、wb满足:
angle(wa)+angle(wb)=Φ
abs(wa)=abs(wb)
6.根据权利要求1至5任一所述的一种卫星导航接收设备波束形成的方法,其特征在于:步骤③中在权值输出时,需要根据阵元的损坏情况,修正权值,在期望波束指向θ0上天线相位中心保持恒定为Φ,则阵元权值为:
w′n=w′n·B·exp[j·β]
其中,w′n为完好阵元的权值,损坏的阵元权值将不用考虑;β为阵元损坏造成的在期望波束指向θ0上产生的天线相位中心偏移;B是对权值幅度修正的系数。
angle(·)表示取一个复数的相位的运算,∑wn是所有阵元权值的和,∑w′n是所有完好阵元权值的和。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010160077 CN101825715B (zh) | 2010-04-30 | 2010-04-30 | 一种卫星导航接收设备波束形成的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010160077 CN101825715B (zh) | 2010-04-30 | 2010-04-30 | 一种卫星导航接收设备波束形成的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101825715A true CN101825715A (zh) | 2010-09-08 |
CN101825715B CN101825715B (zh) | 2013-04-24 |
Family
ID=42689740
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010160077 Expired - Fee Related CN101825715B (zh) | 2010-04-30 | 2010-04-30 | 一种卫星导航接收设备波束形成的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101825715B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102721971A (zh) * | 2011-12-31 | 2012-10-10 | 武汉苍穹数码仪器有限公司 | 一种提高gnss接收机抗干扰性能的方法 |
CN103018713A (zh) * | 2012-12-24 | 2013-04-03 | 中国电子科技集团公司第五十四研究所 | 基于导航数字多波束接收阵列天线的卫星跟踪测角方法 |
CN103064089A (zh) * | 2012-12-28 | 2013-04-24 | 中国电子科技集团公司第五十四研究所 | 卫星导航数字多波束发射阵列天线相位中心的标定方法 |
CN103176191A (zh) * | 2013-03-11 | 2013-06-26 | 西北工业大学 | 一种基于多约束波束形成的伪卫星远近效应抑制方法 |
CN103248412A (zh) * | 2013-04-10 | 2013-08-14 | 中国电子科技集团公司第五十四研究所 | 一种卫星导航多波束形成的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030011514A1 (en) * | 2001-07-11 | 2003-01-16 | Richard Kirchofer | Interference rejection GPS antenna system |
CN1486006A (zh) * | 2002-09-24 | 2004-03-31 | 深圳市中兴通讯股份有限公司 | 智能天线上下行权值转换方法 |
CN1725556A (zh) * | 2005-07-07 | 2006-01-25 | 哈尔滨工业大学 | 基于遗传算法的智能天线复合波束形成方法 |
-
2010
- 2010-04-30 CN CN 201010160077 patent/CN101825715B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030011514A1 (en) * | 2001-07-11 | 2003-01-16 | Richard Kirchofer | Interference rejection GPS antenna system |
CN1486006A (zh) * | 2002-09-24 | 2004-03-31 | 深圳市中兴通讯股份有限公司 | 智能天线上下行权值转换方法 |
CN1725556A (zh) * | 2005-07-07 | 2006-01-25 | 哈尔滨工业大学 | 基于遗传算法的智能天线复合波束形成方法 |
Non-Patent Citations (4)
Title |
---|
《武汉理工大学学报》 20070831 李佳靖; 金荣洪; 耿军平 基于鲍威尔法的快速盲波束形成方法 第29卷, 第8期 * |
《系统仿真学报》 20090930 罗显志; 杨滕; 魏海涛; 李隽; 蔚保国 基于MATLAB的卫星导航信号仿真和验证平台 第21卷, 第18期 * |
《西安邮电学院学报》 20090531 张磊; 卢光跃 随机波束形成技术的仿真研究 第14卷, 第3期 * |
《飞行器测控学报》 20040930 蔚保国; 姚奇松 基于数字多波束天线的多星测控系统 第23卷, 第3期 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102721971A (zh) * | 2011-12-31 | 2012-10-10 | 武汉苍穹数码仪器有限公司 | 一种提高gnss接收机抗干扰性能的方法 |
CN103018713A (zh) * | 2012-12-24 | 2013-04-03 | 中国电子科技集团公司第五十四研究所 | 基于导航数字多波束接收阵列天线的卫星跟踪测角方法 |
CN103018713B (zh) * | 2012-12-24 | 2014-09-24 | 中国电子科技集团公司第五十四研究所 | 基于导航数字多波束接收阵列天线的卫星跟踪测角方法 |
CN103064089A (zh) * | 2012-12-28 | 2013-04-24 | 中国电子科技集团公司第五十四研究所 | 卫星导航数字多波束发射阵列天线相位中心的标定方法 |
CN103064089B (zh) * | 2012-12-28 | 2014-11-05 | 中国电子科技集团公司第五十四研究所 | 卫星导航数字多波束发射阵列天线相位中心的标定方法 |
CN103176191A (zh) * | 2013-03-11 | 2013-06-26 | 西北工业大学 | 一种基于多约束波束形成的伪卫星远近效应抑制方法 |
CN103248412A (zh) * | 2013-04-10 | 2013-08-14 | 中国电子科技集团公司第五十四研究所 | 一种卫星导航多波束形成的方法 |
CN103248412B (zh) * | 2013-04-10 | 2016-01-20 | 中国电子科技集团公司第五十四研究所 | 一种卫星导航多波束形成的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101825715B (zh) | 2013-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109407055A (zh) | 基于多径利用的波束形成方法 | |
CN101825715B (zh) | 一种卫星导航接收设备波束形成的方法 | |
CN105445709B (zh) | 一种稀布阵列近场无源定位幅相误差校正方法 | |
CN103984676A (zh) | 一种基于协方差矩阵重构的正交投影自适应波束形成方法 | |
CN108051782A (zh) | 基于子阵划分的大规模相控阵差波束形成系统 | |
CN105137409B (zh) | 基于幅相约束的目标信号稳健空时自适应处理方法 | |
CN107092007A (zh) | 一种虚拟二阶阵列扩展的波达方向估计方法 | |
CN107526077B (zh) | 一种便携式高频地波雷达浪高反演方法 | |
CN103037494B (zh) | 多波束天线功率控制方法 | |
CN103630887A (zh) | 一种多波束形成方法及使用该方法的多波束声纳 | |
CN112596022B (zh) | 低轨星载多波束正六边形相控阵天线的波达角估计方法 | |
CN101483280B (zh) | 一种稳健波束合成器权值求解方法 | |
CN110824414A (zh) | 一种到达角估计的装置及方法 | |
CN103902830A (zh) | 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法 | |
CN103248412B (zh) | 一种卫星导航多波束形成的方法 | |
CN104931937B (zh) | 基于协方差矩阵归一化的子阵级正交投影波束形成方法 | |
CN109521393A (zh) | 一种基于信号子空间旋转特性的波达方向估计算法 | |
CN111162878B (zh) | 一种基于子阵降维带约束的多域联合抗干扰方法 | |
CN109669172B (zh) | 基于主瓣内强干扰抑制的弱目标方位估计方法 | |
CN103605115A (zh) | 基于网络化雷达阵列合成的抗有源主瓣干扰方法 | |
CN116679306A (zh) | 一种基于系数修正的扩展孔径声纳成像方法 | |
Islam et al. | A computationally efficient near field broadband beamformer | |
Hossain et al. | Efficient robust broadband beamforming using circular antenna arrays | |
CN110554363A (zh) | 一种mimo雷达天线最优保护通道的生成方法 | |
CN116500550B (zh) | 一种星载sar距离模糊抑制方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130424 |