CN102201832B - 基于时分复用的并行匹配滤波器 - Google Patents

基于时分复用的并行匹配滤波器 Download PDF

Info

Publication number
CN102201832B
CN102201832B CN 201110123462 CN201110123462A CN102201832B CN 102201832 B CN102201832 B CN 102201832B CN 201110123462 CN201110123462 CN 201110123462 CN 201110123462 A CN201110123462 A CN 201110123462A CN 102201832 B CN102201832 B CN 102201832B
Authority
CN
China
Prior art keywords
accumulator
register
output
input
data
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
Application number
CN 201110123462
Other languages
English (en)
Other versions
CN102201832A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN 201110123462 priority Critical patent/CN102201832B/zh
Publication of CN102201832A publication Critical patent/CN102201832A/zh
Application granted granted Critical
Publication of CN102201832B publication Critical patent/CN102201832B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

基于时分复用的并行匹配滤波器。它涉及扩频通信领域,它针对传统的匹配滤波器硬件消耗大的缺点。它包括M个移位寄存器、M个乘法器、M+1个数据选择器、一个移位累加器、一个M级的FIFO结构、N个寄存器、一个第一累加器和一个第二累加器;把输入序列与本地整个周期扩频码的相关运算分段进行,最后叠加起来形成全周期的相关值。对每段的相关运算使用相同的寄存器、乘法器和加法单元,并把每段的相关值存入FIFO结构中,最后把各段数据相加形成整个周期的相关值。在保持扩频增益不变的情况下,由于每段数据跟本地数据的相关值计算采用了相同的逻辑资源,因而硬件消耗得到大幅度降低。这适用于长周期伪码的快速捕获。

Description

基于时分复用的并行匹配滤波器
技术领域
本发明涉及扩频通信领域。
背景技术
扩频通信技术由于具有抗干扰能力强,截获概率低和多址接入能力强等优点,因而在军事和民用上都有广泛的应用,例如全球定位系统(GPS)和第三代移动通信技术。然而,要充分发挥扩频通信的优点,接收机必须首先捕获接收扩频码序列,否则将会导致后续的解扩解调过程难以进行。伪随机码的捕获是扩频通信的关键技术之一。扩频接收中的同步捕获技术是决定扩频系统能否正常工作的关键。传统的捕获方法有:串行捕获法、顺序估计捕获法、匹配滤波器同步捕获法等。
下面分析传统的串行捕获方法和匹配滤波器同步捕获法,不失一般性,采用的扩频码长度为255,匹配滤波器的输入数据的位宽为12bits,每个码片采集一个点。
传统的串行捕获方法如图1所示,本地PN码序列发生器产生的扩频码C1,C2,…,C255与输入数据di+1,di+2,…分别相乘,送入到一个积分清零器。当积分清零器的累加数据达到255个时,积分的结果即为255个输入数据和本地255个扩频码的相关值R(i),如公式一所示,
Figure GDA00002948657700011
(公式一)。然后累加器清零,并通过判决模块对R(i)进行判断,如果R(i)大于阈值η,则认为捕获成功,并发送信号到PN码跟踪模块,否则产生相位调整信号,调整本地PN码序列的相位,在下一个相位上重新进行相关运算。由此可见,在最坏的情况下,需要搜索255个相位才能完成捕获。因此,串行捕获的捕获时间比较长,不适合快速捕获的要求,也发挥不出硬件实现的特点。
传统的并行匹配滤波器的结构如图2所示。其由一组移位寄存器组成,di+1~di+255是存储在移位寄存器中的输入数据。C1~C255是本地一个周期的扩频码,每个寄存器中的数据跟本地的扩频码进行相乘累加得到输入序列与本地扩频码的相关值R(i),并作为匹配滤波器的输出,如公式一所示。当匹配滤波器的输出超过一定的阈值,表明捕获成功。采用这种结构,输入序列在一个扩频码周期内就完成捕获。但是,当扩频码的长度较长时,会导致匹配滤波器占用很多的寄存器、乘法器等逻辑资源。在图2的结构中,移位寄存器占用了255×12=3060个寄存器,系数相乘需要消耗255个乘法器。匹配滤波器所消耗的逻辑资源也会随着扩频码长度的增大而急剧增加。因此,扩频码的长度受到硬件资源的限制,降低硬件消耗是实现匹配滤波技术的关键。因此从上述论述中可以得出并行运算匹配滤波器同步捕获法的捕获速度最快,只需要一个扩频码周期就可以完成捕获,但是其代价是硬件消耗大,当扩频码长度很长时,巨大的硬件资源消耗往往导致其难以实现。
发明内容
本发明针对传统的匹配滤波器硬件消耗大的缺点,提出了基于时分复用的并行匹配滤波器。
本发明的基于时分复用的并行匹配滤波器包括M个移位寄存器、M个乘法器、M+1个数据选择器、一个移位累加器、一个M级的先进先出(FIFO)结构、N个寄存器、一个第一累加器和一个第二累加器;其中,N为输入数据与本地整个周期的扩频码的相关函数R(i)所分成的计算段数,本地整个周期的扩频码周期除以所述的段数进位取整为M,其中N和M为整数,M个移位寄存器、M个乘法器和M个数据选择器组成M组运算器,每一组运算器均包括一个移位寄存器、一个乘法器和一个数据选择器,所述的移位寄存器的输出端与所述的乘法器的一个输入端连接,所述的乘法器的另一个输入端与所述的数据选择器的输出端连接,所述的数据选择器的N个输入端分别输入第i个扩频码、第i+M×1个扩频码至第i+M×(N-1)个扩频码,并且每一组运算器中的移位寄存器、乘法器和数据选择器的连接关系相同,所述的M个移位寄存器为串行输入并行输出,所述的M个乘法器的输出端均与第一累加器的输入端连接,所述的M个数据选择器的受控端与移位累加器的控制端连接,所述的移位累加器的控制端还与第M+1个数据选择器的受控端连接,第一累加器的输出端与第M+1个数据选择器的输入端连接,第M+1个数据选择器的N个输出端分别与M级的先进先出(FIFO)结构的输入端和第N个寄存器的输入端连接,M级的先进先出(FIFO)结构的输出端分别与N-1个寄存器的输入端连接,N-1个寄存器的输出端和第N个寄存器的输出端同时与第二累加器的输入端连接,第二累加器的输出端为匹配滤波输出端。
本发明的基于时分复用的并行匹配滤波器,把输入序列与本地整个周期扩频码的相关运算分段进行,最后叠加起来形成全周期的相关值。对每段的相关运算使用相同的寄存器、乘法器和加法单元,并把每段的相关值存入FIFO结构中,最后把各段数据相加形成整个周期的相关值。在保持扩频增益不变的情况下,由于每段数据跟本地数据的相关值计算采用了相同的逻辑资源,因而硬件消耗得到大幅度降低。同时随着复用次数的增加,硬件消耗能够进一步降低。这对于长周期伪码的快速捕获具有重要的实用价值。
附图说明
图1是传统串行捕获方法的滤波器的结构示意图,图2是传统并行匹配滤波器的结构示意图,图3是具体实施方式二的结构示意图,图4是具体实施方式三的结构示意图。
具体实施方式
具体实施方式一:结合图3和图4说明本实施方式,本实施方式包括M个移位寄存器1、M个乘法器2、M+1个数据选择器3、一个移位累加器4、一个M级的FIFO结构(First In First Out先入先出结构)5、N个寄存器Reg、一个第一累加器6和一个第二累加器7;
其中,N为输入数据与本地整个周期的扩频码的相关函数R(i)所分成的计算段数,本地整个周期的扩频码的周期除以所述的段数进位取整为M,其中N为整数,M为整数,
M个移位寄存器1、M个乘法器2和M个数据选择器3组成M组运算器,每一组运算器均包括一个移位寄存器1、一个乘法器2和一个数据选择器3,所述的移位寄存器1的输出端与所述的乘法器2的一个输入端连接,所述的乘法器2的另一个输入端与所述的数据选择器3的输出端连接,所述的数据选择器3的N个输入端分别输入第i个扩频码、第i+M×1个扩频码至第i+M×(N-1)个扩频码,并且每一组运算器中的移位寄存器1、乘法器2和数据选择器3的连接关系相同,
所述的M个移位寄存器1为串行输入并行输出,
所述的M个乘法器2的输出端均与第一累加器6的输入端连接,
所述的M个数据选择器3的受控端与移位累加器4的控制端连接,
所述的移位累加器4的控制端还与第M+1个数据选择器3的受控端连接,
第一累加器6的输出端与第M+1个数据选择器3的输入端连接,
第M+1个数据选择器3的N个输出端分别与M级的FIFO结构5的输入端和第N个寄存器Reg的输入端连接,
M级的FIFO结构5的输出端与第N-1个寄存器Reg的输入端连接,
第N-1个寄存器Reg的输出端和第N个寄存器Reg的输出端同时与第二累加器7的输入端连接,
第二累加器7的输出端为匹配滤波输出端。
具体实施方式二:结合图3说明本实施方式,本实施方式与具体实施方式一不同点在于输入数据与本地整个周期的扩频码的相关函数R(i)分成两段来计算,记N=2,则包括128个移位寄存器1、128个乘法器2、129个数据选择器3、1个移位累加器4、1个FIFO结构5、2个寄存器Reg、1个第一累加器6和1个第二累加器7;其连接关系与具体实施方式相同,由此可以得到本实施方式的基于时分复用的并行匹配滤波器的结构,如图3所示。本实施方式的基于时分复用的并行匹配滤波器的工作过程如下:
如公式二所示
R ( i ) = Σ m = 1 255 C m d i + m = Σ m = 1 128 C m d i + m + Σ m = 129 255 C m d i + m    公式二
令n=m-128,则有
R ( i ) = Σ m = 1 128 C m d i + m + Σ n = 1 127 C n + 128 d i + n + 128    公式三
在此引入两个变量来表示这两段的计算结果,如公式四和公式五所示:
R 1 ( i ) = Σ m = 1 128 C m d i + m    公式四
R 2 ( i ) = Σ n = 1 127 C n + 128 d i + n    公式五
则由公式五可得
R 2 ( i + 128 ) = Σ n = 1 127 C n + 128 d i + n + 128    公式六
则对于先输入的128个数据di+1,di+2…di+128,其与本地前128个扩频码C1,C2…C128的相关值为R1(i),而di+1,di+2…di+127与本地后127个扩频码C129,C130…C255的相关值为R2(i)。那么由公式三、公式四和公式六可得输入序列与本地整个周期扩频码的相关函数R(i)可以分解成R1(i)与R2(i+128)之和,如公式七所示。
R(i)=R1(i)+R2(i+128)   公式七
本实施方式的基于时分复用的并行匹配滤波器的输入移位寄存器个数截短为传统匹配滤波器的一半,也就是128个。滤波器的工作时钟为clk,其频率为输入数据率的两倍。每输入一个数据,基于时分复用的并行匹配滤波器需要两个clk才能计算出输入的255个数据和本地整个周期扩频码的相关值:第一个clk,计算出di+1,di+2,…,di+128和C1,C2,…,C128的相关值R1(i)并存入到一个128级的FIFO结构5中,同时FIFO寄存器5把其最前面的数据R1(i-128)弹出,存入到第二寄存器Reg2中,并且将移位寄存器1的抽头系数切换为C129,C130,…,C255,0;第二个clk,计算出di+1,di+2,…,di+128和C129,C130,…,C255,0的相关值R2(i)并且存入第一寄存器Reg1中,同时将移位寄存器1的抽头系数切换为C1,C2,…,C128,并且更新输入数据为di+2,di+3,…,di+129。这样第一寄存器Reg1和第二寄存器Reg2之和R1(i-128)+R2(i)则为di-128,di-127,…,di+126,di+127和C1,C2,…,C255的相关值R(i-128),并作为匹配滤波器的输出。
表1列出了基于时分复用的并行匹配滤波器中移位寄存器1、FIFO结构5、第一寄存器Reg1、第二寄存器Reg2和最终输出的关系,从该表中可以看出,第一寄存器Reg1存储着当前寄存器中的127个输入数据与C129,C130,…,C255的相关和,而第二寄存器Reg2存放的是刚从移位寄存器1中移出的128个数据与C1,C2,…,C128的相关和,两者之和恰好为255个输入数据与C1,C2,...,C255的相关值R(i-128)。
表1基于时分复用的并行匹配滤波器中的数据变化情况表
下面分析基于时分复用的并行匹配滤波器的硬件消耗:输入数据移位寄存器1使用12×128=1536个寄存器;系数相乘消耗128个乘法器;FIFO寄存器5的长度为128,每个存储单元为19位,共使用128×19=2432bit RAM。可见基于时分复用的并行匹配滤波器所需的寄存器和乘加单元是传统匹配滤波器长度的一半。
具体实施方式三:结合图4说明本实施方式,本实施方式与具体实施方式一不同点在于输入数据与本地整个周期的扩频码的相关函数R(i)分成四段来计算,记N=4,则包括64个移位寄存器1、64个乘法器2、65个数据选择器3、1个移位累加器4、3个FIFO结构5、4个寄存器Reg、1个第一累加器6和1个第二累加器7;其连接关系与具体实施方式相同,由此可以得到本实施方式的基于时分复用的并行匹配滤波器的结构,如图4所示。本实施方式的基于时分复用的并行匹配滤波器的结构为提升结构,该结构可以进一步降低硬件消耗。本实施方式的基于时分复用的并行匹配滤波器的工作过程如下:
如公式八所示:
R ( i ) = Σ m = 1 255 C m d i + m = Σ m = 1 64 C m d i + m + Σ m = 65 128 C m d i + m + Σ m = 129 192 C m d i + m + Σ m = 193 255 C m d i + m   公式八
令n=m-64,k=m-128,l=m-172,做变量代换得公式九
R ( i ) = Σ m = 1 64 C m d i + m + Σ n = 1 64 C n + 64 d i + n + 64 + Σ k = 1 64 C k + 128 d i + k + 128 + Σ l = 1 63 C l + 192 d i + l + 192   公式九
相应地把本地的255个扩频码分成4段:C1,C2,…,C64;C65,C66,…,C128;C129,C130,…,C192和C193,C194,…,C255。定义输入的64个数据di+1,di+2,…,di+64与这四段扩频码的相关值分别为R1′(i),R2′(i),R3′(i)和R4′(i),如公式十、公式十一、公式十二和公式十三所示。
R 1 ′ ( i ) = Σ m = 1 64 C m d i + m      公式十
R 2 ′ ( i ) = Σ n = 1 64 C n + 64 d i + n      公式十一
R 3 ′ ( i ) = Σ k = 1 64 C k + 128 d i + k      公式十二
R 4 ′ ( i ) = Σ l = 1 63 C l + 192 d i + l      公式十三
因此R(i)分解成R1′(i),R2′(i+64),R3′(i+128)和R4′(i+192)之和,如公式十四所示。
R(i)=R′1(i)+R′2(i+64)+R′3(i+128)+R′4(i+192)公式十四
匹配滤波器的长度可以进一步截短为传统匹配滤波器的1/4,也就是64级。输入数据时钟为clk,其频率为输入数据率的4倍。在连续的3个clk中,其分别计算出R1′(i),R2′(i),R3′(i)并存入在3个FIFO结构5中,同时这3个FIFO结构5的最前面的数据弹出,分别存储在第一寄存器Reg1,第二寄存器Reg2和第三寄存器Reg3中。由于这3个FIFO寄存器5的级数分别为192,128和64,则这3个寄存器Reg中存储的数据分别为R1′(i-192),R2′(i-128),R3′(i-64)。在第4个clk2时钟求出R4′(i),并存入第四寄存器Reg4中。第一寄存器Reg1,第二寄存器Reg2、第三寄存器Reg3和第四寄存器Reg4中的数据相加,得到R1′(i-192)+R2′(i-128)+R3′(i-64)+R4′(i)=R(i-192),并作为匹配滤波器的输出。
下面分析其硬件消耗:移位寄存器1使用12×64=768个寄存器,系数相乘需要消耗64个乘法器2,FIFO寄存器5的总长度为192+128+64=384,每个存储单元18位,共使用384×18=6912bits RAM。由于移位寄存器1的长度为传统匹配滤波器长度的1/4,并且FIFO寄存器5使用的是RAM资源,因此该滤波器的硬件消耗降低为传统匹配滤波器的1/4。在芯片允许时钟范围内,通过加大复用的次数,可以显著的降低匹配滤波器的硬件资源的消耗。
以Altera公司的Stratix II系列FPGA芯片EP2S60F1020C5为平台,利用VHDL语言编写传统匹配滤波器模块和基于图3的基于时分复用的并行匹配滤波器,编译后两者的资源使用情况如表2所示。从表中可以看出,基于时分复用的并行匹配滤波器相对于传统匹配滤波器能节省将近一半的寄存器和组合逻辑消耗,虽然增加了RAM的消耗,但是由于FPGA的RAM的资源远大于寄存器和组合逻辑的资源,因此从总体上大大降低了硬件的开销。在FPGA中,RAM的资源通常远远大于寄存器的资源,因此通过增加RAM来减少寄存器使总体硬件开销变小。
表2传统匹配滤波器和基于时分复用的并行匹配滤波器的硬件资源消耗对比
占用芯片资源/资源总量 传统匹配滤波器 时分复用并行匹配滤波器
寄存器 5397/48352 2776/48352
组合逻辑 4472/48352 2799/48352
RAM 0/2544192 2432/2544192
本发明内容不仅限于上述各实施方式的内容,其中一个或几个具体实施方式的组合同样也可以实现发明的目的。

Claims (1)

1.基于时分复用的并行匹配滤波器,其特征在于它包括M个移位寄存器(1)、M个乘法器(2)、M+1个数据选择器(3)、一个移位累加器(4)、一个M级的先进先出(FIFO)结构(5)、N个寄存器(Reg)、一个第一累加器(6)和一个第二累加器(7);其中,N为输入数据与本地整个周期的扩频码的相关函数R(i)所分成的计算段数,本地整个周期的扩频码周期除以所述的段数进位取整为M,其中N和M为整数,M个移位寄存器(1)、M个乘法器(2)和M个数据选择器(3)组成M组运算器,每一组运算器均包括一个移位寄存器(1)、一个乘法器(2)和一个数据选择器(3),所述的移位寄存器(1)的输出端与所述的乘法器(2)的一个输入端连接,所述的乘法器(2)的另一个输入端与所述的数据选择器(3)的输出端连接,所述的数据选择器(3)的N个输入端分别输入第i个扩频码、第i+M×1个扩频码至第i+M×(N-1)个扩频码,并且每一组运算器中的移位寄存器(1)、乘法器(2)和数据选择器(3)的连接关系相同,所述的M个移位寄存器(1)为串行输入并行输出,所述的M个乘法器(2)的输出端均与第一累加器(6)的输入端连接,所述的M个数据选择器(3)的受控端与移位累加器(4)的控制端连接,所述的移位累加器(4)的控制端还与第M+1个数据选择器(3)的受控端连接,第一累加器(6)的输出端与第M+1个数据选择器(3)的输入端连接,第M+1个数据选择器(3)的N个输出端分别与M级的先进先出(FIFO)结构(5)的输入端和第N个寄存器(Reg)的输入端连接,M级的先进先出(FIFO)结构(5)的输出端分别与N-1个寄存器(Reg)的输入端连接,N-1个寄存器(Reg)的输出端和第N个寄存器(Reg)的输出端同时与第二累加器(7)的输入端连接,第二累加器(7)的输出端为匹配滤波输出端。
CN 201110123462 2011-05-13 2011-05-13 基于时分复用的并行匹配滤波器 Expired - Fee Related CN102201832B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110123462 CN102201832B (zh) 2011-05-13 2011-05-13 基于时分复用的并行匹配滤波器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110123462 CN102201832B (zh) 2011-05-13 2011-05-13 基于时分复用的并行匹配滤波器

Publications (2)

Publication Number Publication Date
CN102201832A CN102201832A (zh) 2011-09-28
CN102201832B true CN102201832B (zh) 2013-07-03

Family

ID=44662253

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110123462 Expired - Fee Related CN102201832B (zh) 2011-05-13 2011-05-13 基于时分复用的并行匹配滤波器

Country Status (1)

Country Link
CN (1) CN102201832B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104765054A (zh) * 2015-04-20 2015-07-08 和芯星通科技(北京)有限公司 一种gnss中伪随机码的捕获方法及装置
CN105388501B (zh) * 2015-11-05 2017-12-01 天津津航计算技术研究所 一种星载北斗通信系统快速捕获方法
CN111510110B (zh) * 2020-04-30 2023-10-20 上海大学 一种并行处理的插值匹配滤波方法及滤波器

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1274841A (zh) * 1999-05-21 2000-11-29 浙江省中医药研究院 正品桑白皮的红外光谱鉴别法
US6775684B1 (en) * 1999-06-03 2004-08-10 Sharp Kabushiki Kaisha Digital matched filter

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1274841A (zh) * 1999-05-21 2000-11-29 浙江省中医药研究院 正品桑白皮的红外光谱鉴别法
US6775684B1 (en) * 1999-06-03 2004-08-10 Sharp Kabushiki Kaisha Digital matched filter

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李向东,郭树旭,赵蔚,贺蓉.《CDMA并行匹配滤波器的CPLD设计》.《吉林大学学报(理工版)》.2002,第41卷(第2期),全文. *

Also Published As

Publication number Publication date
CN102201832A (zh) 2011-09-28

Similar Documents

Publication Publication Date Title
EP2294518B1 (en) Adaptive correlation
CN111694027B (zh) 超大动态扩频信号捕获方法与装置
CN105791184B (zh) 基于相似度的非周期长码直扩信号扩频序列盲估计方法
CN102540217A (zh) 一种在fpga中实现glonass卫星信号的快速捕获系统
CN102201832B (zh) 基于时分复用的并行匹配滤波器
CN103885071A (zh) 一种卫星信号捕获的方法和设备
CN102025391A (zh) 实现多种扩频码长的高速滤波装置及其方法
CN101174849B (zh) 无线传感网节点的扩频码片同步捕获和跟踪方法及其装置
CN101494506B (zh) 一种流水结构的水声扩频通信带通信号捕获电路
CN103837878A (zh) 一种gnss卫星信号捕捉方法
CN101384992A (zh) 用来执行乘加运算的方法和装置
CN101494505A (zh) 一种面积优化的流水结构水声扩频通信带通信号捕获电路
CN108919200A (zh) 一种基于fpga的任意雷达波形在线设计方法
CN102710283A (zh) 直接序列扩频伪码捕获方法、捕获装置及通信系统
CN102185628A (zh) 基于自适应功率积累的扩频码相位捕获装置及其捕获方法
CN105388501B (zh) 一种星载北斗通信系统快速捕获方法
CN102323604B (zh) 一种加法器和相关器一体装置
US6601078B1 (en) Time-efficient real-time correlator
CN101762803B (zh) 一种应用于全球导航卫星系统接收机的差分相关器
CN104460444A (zh) 基于广义相关系数的fpga运算电路
CN102723967B (zh) 基于匹配滤波的相关运算电路
CN100438360C (zh) 处理扩频信号的方法、设备和接收扩频信号的接收机
Wang et al. Optimized fpga realization of digital matched filter in spread spectrum communication systems
Sun et al. Multi-channel signal acquisition algorithm for multi-beam satellite systems in high dynamic environment
CN106199656B (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: 20130703

Termination date: 20150513

EXPY Termination of patent right or utility model