CN101800719B - 离散导频同步方法与装置及其通信系统 - Google Patents

离散导频同步方法与装置及其通信系统 Download PDF

Info

Publication number
CN101800719B
CN101800719B CN 200910006341 CN200910006341A CN101800719B CN 101800719 B CN101800719 B CN 101800719B CN 200910006341 CN200910006341 CN 200910006341 CN 200910006341 A CN200910006341 A CN 200910006341A CN 101800719 B CN101800719 B CN 101800719B
Authority
CN
China
Prior art keywords
impulse response
channel impulse
scattered pilot
channel
circuit
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 200910006341
Other languages
English (en)
Other versions
CN101800719A (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.)
Ali Corp
Original Assignee
Ali Corp
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 Ali Corp filed Critical Ali Corp
Priority to CN 200910006341 priority Critical patent/CN101800719B/zh
Publication of CN101800719A publication Critical patent/CN101800719A/zh
Application granted granted Critical
Publication of CN101800719B publication Critical patent/CN101800719B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Mobile Radio Communication Systems (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

本发明的示范实施例提供了一种离散导频同步方法与装置及其通信系统。离散导频同步方法可用于数字视频广播系统,并具有以下步骤:(1)分别萃取频域符号的第一至第x位置集合的多个子载波,其中,x为大于1的整数;(2)获得第一至第x通道脉冲响应,其中,第i通道脉冲响应是根据第i位置的子载波所获得第,且i为1至x的整数;(3)根据第一至第x通道脉冲响应决定频域符号的多个离散导频的正确位置。

Description

离散导频同步方法与装置及其通信系统
技术领域
本发明是有关于一种离散导频同步方法,且特别是有可应用于数字视频广播(Digital Video Broadcast,简称为DVB)系统中的基于时域的一种离散导频同步方法与装置及其通信系统。
背景技术
近年来,使用DVB系统的电子装置已经逐渐普及。在无线传输通道中,由于多路径与时变效应的关系,接收端必须知道无线传输通道的通道信息才能顺利地将接收到的信号进行解调并正确地获得传输端所传送的信息。DVB系统的传送端在传送的过程中会在每一个符号的频域上插入离散导频(Scattered Pilot),因此,在DVB系统的接收端进行解调时,能够依赖每一个符号的频域上的离散导频来获得无线传输通道信息。
由于离散导频插入的位置在各符号内不完全相同,所以在接收端进行解调前,接收端必须知道离散导频的位置,其中,获得离散导频的位置的方法又称为离散导频同步方法。传统的离散导频同步方法使用接收端对传输参数信令(Transmission Parameter Signaling,简称为TPS)进行解调,并藉此得知离散导频的位置。但是,在DVB系统中,TPS仅会在每一个帧(fame)中出现一次,而每一个帧中包括了68个符号。因此,传统的离散导频同步方法必须耗费至少68个符号的时间才能得知离散导频的位置。
请参照图1,图1是DVB系统中的符号的离散导频在频域上的位置的示意图。如同图1所示,符号中的离散导频在频域上的位置不完全相同,每一个符号在频域上的第0个位置与第Kmax个位置上都有离散导频。第(4×i)个符号在第0个位置之后与第Kmax个位置之前的离散导频位于第(12+12×p)个位置,第(1+4×i)个符号在第0个位置之后与第Kmax个位置之前的离散导频位于第(3+12×p)个位置,第(2+4×i)个符号在第0个位置之后与第Kmax个位置之前的离散导频位于第(6+12×p)个位置,而第(3+4×i)个符号在第0个位置之后与第Kmax个位置之前的离散导频位于第(9+12×p)个位置。其中,i为任意整数,而p等于0至pmax的任意整数。
另外,若DVB系统是采用2K模式时,则Kmax等于1704,第(1+4×i)个符号、第(2+4×i)个符号与第(3+4×i)个符号的pmax为141,而第(4×i)个符号的pmax为140;若DVB系统是采用4K模式时,则Kmax等于3408,第(1+4×i)个符号、第(2+4×i)个符号与第(3+4×i)个符号的pmax为283,而第(4×i)个符号的pmax为282;若DVB系统是采用8K模式时,则Kmax等于6816,第(1+4×i)个符号、第(2+4×i)个符号与第(3+4×i)个符号的pmax为567,而第(4×i)个符号的pmax为566。
再参照图1,DVB标准规定离散导频的起始位置以四个符号为周期,离散导频的能量为数据的能量的4/3倍,位于同一个频域位置上的不同符号的离散导频彼此相同,而符号内的每一个离散导频彼此不相同。
为了减少离散导频同步方法的处理时间,美国专利公开申请案US2006/0088133A1提出的一种基于能量的离散导频同步方法的示意图,此基于能量的离散导频同步方法会在一个符号内提取出离散导频的四种可能位置的能量。由于离散导频的能量比数据的能量高,因此,通过比较能量的方式,便可以在四种可能位置中找到离散导频的正确位置。
请参照图2,图2是基于能量的离散导频同步方法的估计能量的示意图。其中,纵轴表示时间,横轴表示频率。基于能量的离散导频同步方法在对第x个符号进行离散导频同步时,会将计算频域上第12、24、36、...、12+12×px个位置的子载波的能量总和E1(x)、频域上第3、15、27、...、3+12×px个位置的子载波的能量总和E2(x)、频域上第6、18、30、...、6+12×px个位置的子载波的能量总和E3(x)与频域上第9、21、33、...、9+12×px个位置的子载波的能量总和E4(x)。能量总和E1(x)、E2(x)、E3(x)与E4(x)的计算方式分别表示如下:
E 1 ( x ) = | Σ p = 0 p x S ( x , 12 × p + 12 ) · S * ( x , 12 × p + 12 ) | ,
E 2 ( x ) = | Σ p = 0 p x S ( x , 12 × p + 3 ) · S * ( x , 12 × p + 3 ) | ,
E 3 ( x ) = | Σ p = 0 p x S ( x , 12 × p + 6 ) · S * ( x , 12 × p + 6 ) | ,
E 4 ( x ) = | Σ p = 0 p x S ( x , 12 × p + 9 ) · S * ( x , 12 × p + 9 ) | .
其中,px为大于0且小于pmax的整数,S(x,i)表示第x个符号中频域上第i个子载波的信号值,S*(x,i)表示第x个符号中频域上第i个子载波的共轭(conjugated)信号值。接着,基于能量的离散导频同步方法会比较能量总和E1(x)、E2(x)、E3(x)与E4(x),而能量总和最大者所表示的位置即为离散导频的正确位置。根据以上所述,基于能量的离散导频同步方法仅需要一个符号时间就能完成离散导频同步。
接着,请参照图3,图3是基于相关性(correlation)的离散导频同步方法的计算相关性的示意图。如同图3所示,基于相关性的离散导频同步方法在对第x个符号进行离散导频同步时,会计算第x个与第(x-4)个符号的频域上第12、24、36、...、12+12×px个位置的子载波的相关性总和C1(x)、第x个与第(x-4)个符号的频域上第3、15、27、...、3+12×px个位置的子载波的相关性总和C2(x)、第x个与第(x-4)个符号的频域上第6、18、30、...、6+12×px个位置的子载波的相关性总和C3(x)与第x个与第(x-4)个符号的频域上第9、21、33、...、9+12×px个位置的子载波的相关性总和C4(x)。相关性总和C1(x)、C2(x)、C3(x)与C4(x)的表示如下:
C 1 ( x ) = | Σ p = 0 p x S ( x , 12 × p + 12 ) · S * ( x - 4 , 12 × p + 12 ) | ,
C 2 ( x ) = | Σ p = 0 p x S ( x , 12 × p + 3 ) · S * ( x - 4 , 12 × p + 3 ) | ,
C 3 ( x ) = | Σ p = 0 p x S ( x , 12 × p + 6 ) · S * ( x - 4 , 12 × p + 6 ) | ,
C 4 ( x ) = | Σ p = 0 p x S ( x , 12 × p + 9 ) · S * ( x - 4 , 12 × p + 9 ) | .
接着,基于相关性的离散导频同步方法会比较相关性总和C1(x)、C2(x)、C3(x)与C4(x),而相关性总和最大者所表示的位置即为离散导频的正确位置。根据以上所述,基于相关性的离散导频同步方法仅需要五个符号时间就能完成离散导频同步。
接着,请参照图4,图4是时频相关性的离散导频同步方法计算时频相关性的示意图。如同图4所示,基于时频相关性的离散导频同步方法在对第x个符号进行离散导频同步时,会计算第x个符号的频域上第12、24、36、...、12+12×px个位置与第(x-1)个符号的频域上第9、21、33、...、9+12×px个位置的子载波的时频相关性总和TFC1(x)、第x个符号的频域上第3、15、27、...、3+12×px个位置与第(x-1)个符号的频域上第0、12、24、...、12×px个位置的子载波的时频相关性总和TFC2(x)、第x个符号的频域上第6、18、30、...、6+12×px个位置与第(x-1)个符号的频域上第3、15、27、...、3+12×px个位置的子载波的时频相关性总和TFC3(x)与第x个符号的频域上第9、21、33、...、9+12×px个位置与第(x-1)个符号的频域上第6、18、30、...、6+12×px个位置的子载波的时频相关性总和TFC4(x)。相关性总和TFC1(x)、TFC2(x)、TFC3(x)与TFC4(x)的表示如下:
TFC 1 ( x ) = | Σ p = 0 p x ( S ( x , 12 × p + 12 ) · P * ( 12 × p + 12 ) ) · S * ( x - 1,12 × p + 9 ) · P ( 12 × p + 9 ) | ,
TFC 2 ( x ) = | Σ p = 0 p x ( S ( x , 12 × p + 3 ) · P * ( 12 × p + 3 ) ) · ( S * ( x - 1,12 × p ) · P ( 12 × p ) ) | ,
TFC 3 ( x ) = | Σ p = 0 p x ( S ( x , 12 × p + 6 ) · P * ( 12 × p + 6 ) ) · ( S * ( x - 1,12 × p + 3 ) · P ( 12 × p + 3 ) ) | ,
TFC 4 ( x ) = | Σ p = 0 p x ( S ( x , 12 × p + 9 ) · P * ( 12 × p + 9 ) ) · ( S * ( x - 1,12 × p + 9 ) · P ( 12 × p + 9 ) ) | .
其中,P(x)表示接收端已知的频域上第x个位置的离散导频的信号值,P*(x)表示接收端已知的频域上第x个位置的离散导频的共轭信号值。接着,基于时频相关性的离散导频同步方法会比较时频相关性总和TFC1(x)、TFC2(x)、TFC3(x)与TFC4(x),而时频相关性总和最大者所表示的位置即为离散导频的正确位置。根据以上所述,基于时频相关性的离散导频同步方法仅需要两个符号时间就能完成离散导频同步。
据此,虽然已知技术已提供了多种快速的离散导频同步方法,然而,上述的离散导频同步方法却无法在多种不同的通道环境都能够具有很好的效能。例如,在加性白色高斯噪声(Additive White Gaussian Noise,简称为AWGN)通道下,基于能量的离散导频同步方法可能会得到错误的位置;在时变(timevaring)通道下,基于相关性的离散导频同步方法可能会得到错误的位置;而在多径(multi path)通道下,基于时频相关性的离散导频同步方法可能会得到错误的位置。
发明内容
本发明的示范实施例提供了一种离散导频同步方法。此离散导频同步方法用于一通信系统,且此基于时域的离散导频同步方法包括以下步骤:(1)接收频域符号;(2)分别萃取所述频域符号的第一至第x位置集合的多个子载波,其中,x为大于1的整数;(3)获得第一至第x通道脉冲响应,其中,第i通道脉冲响应是根据第i位置的子载波所获得,i为1至x的整数;(4)根据第一至第x通道脉冲响应决定频域符号的多个离散导频的位置。
本发明的示范实施例提供了一种离散导频同步装置,此离散导频同步装置用于通信系统,其包括符号寄存器、至少一通道脉冲响应计算电路与位置决定电路。其中,通道脉冲响应计算电路耦接于符号寄存器,位置决定电路耦接于通道脉冲响应计算电路。符号寄存器用以接收频域符号,通道脉冲响应计算电路用以萃取频域符号的第一至第x位置集合的多个子载波,且通道脉冲响应计算电路更用以获得第一至第x通道脉冲响应。其中,x为大于1的整数,第i通道脉冲响应是根据第i位置的子载波计算得来,i为1至x的整数。位置决定电路根据第一至第x通道脉冲响应决定频域符号的多个离散导频的位置。
本发明的示范实施例提供了一种通信系统,此通信系统包括接收端以及发射端。接收端包括基于时域的离散导频同步装置,此基于时域的离散导频同步装置包括符号寄存器、至少一通道脉冲响应计算电路与位置决定电路。其中,通道脉冲响应计算电路耦接于符号寄存器,位置决定电路耦接于通道脉冲响应计算电路。符号寄存器用以接收频域符号,通道脉冲响应计算电路用以萃取频域符号的第一至第x位置集合的多个子载波且通道脉冲响应计算电路更用以获得第一至第x通道脉冲响应。其中,x为大于1的整数,第i通道脉冲响应是根据第i位置的子载波计算得来,i为1至x的整数。位置决定电路根据第一至第x通道脉冲响应决定频域符号的多个离散导频的位置。
基于所述,本发明的示范实施例所提供的基于时域的离散导频同步方法、装置及其通信系统在各种不同的通道环境皆能够具有良好的效能,并能够准确地检测出离散导频的正确位置。
附图说明
图1是DVB系统中的符号的离散导频在频域上的位置的示意图。
图2是基于能量的离散导频同步方法的估计能量的示意图。
图3是基于相关性(correlation)的离散导频同步方法的计算相关性的示意图。
图4是时频相关性的离散导频同步方法计算时频相关性的示意图。
图5A~图5D分别是通道脉冲响应
Figure GDA0000131549670000071
Figure GDA0000131549670000072
的通道脉冲响应示意图。
图6是本发明的示范实施例所提供的基于时域的离散导频同步装置的电路方块图。
图7是本发明的示范实施例所提供的基于时域的离散导频同步方法的流程图。
图8是中国移动多媒体广播系统中的符号的离散导频在频域上的位置的示意图。
图9是基于能量的离散导频同步方法、基于相关性的离散导频同步方法、基于时频相关性的离散导频同步方法与基于时域的离散导频同步方法在AWGN通道下的效能示意图。
图10是基于能量的离散导频同步方法、基于相关性的离散导频同步方法、基于时频相关性的离散导频同步方法与基于时域的离散导频同步方法在时变通道下的效能示意图。
图11是基于能量的离散导频同步方法、基于相关性的离散导频同步方法、基于时频相关性的离散导频同步方法与基于时域的离散导频同步方法在多径通道下的效能示意图。
附图标号
60:基于时域的离散导频同步装置
61:第一通道脉冲响应估计电路
62:第二通道脉冲响应估计电路
63:第三通道脉冲响应估计电路
64:第四通道脉冲响应估计电路
65:位置决定电路
66:符号寄存器
611:第一萃取器
621:第二萃取器
623:第三萃取器
624:第四萃取器
612、622、632、642:通道频率响应计算电路
613、623、633、643:IFFT电路
651~654:数据处理电路
656:比较电路
S70~S75:步骤
S801~804、S901~S904、S1001~S1004:曲线
具体实施方式
为让本发明的上述特征和优点能更明显易懂,下文特举实施例,并配合所附附图作详细说明如下。
本发明的示范实施例提供了一种离散导频同步装置,此离散导频同步装置用于DVB系统。DVB系统包含了发射端与接收端,而上述的离散导频同步系统配置于DVB系统的接收端。另外,DVB系统不限定是DVB-T(Terrestrial)或DVB-H(Handhold)规格的DVB系统。
由于离散导频在每一个符号内为等间隔(间隔为12)分布,并且在离散导频上所传输的子载波的信号值对接收端来说为已知的信号值,因此,上述的多个离散导频会与通道频率响应Hx的取样(sampling)有关。通道频率响应Hx包含多个取样点,在此用Hx,k来表示第x个符号频域上第k个位置所估计出的通道频率响应的取样点,其中,k表示离散导频的子载波的位置。第x个符号频域上第k个位置所估计出的通道频率响应的取样点Hx,k可以表示为:
H x , k = Y x , k / P x , k = Y x , k · P x , k * / | P x , k | 2 .
其中,Px,k表示第x个符号频域上第k个位置已知的离散导频的子载波的信号值。
接着,将利用离散导频所计算出的通道频率响应Hx进行补零(zero padding)的动作,以使得补零后的通道频率响应Hx具有2N个点,其中,N为正整数。接着,对补零后的通道频率响应Hx进行反快速傅立叶转换(Inversed FastFourier Transform,IFFT)以获得第x个符号的通道脉冲响应hx。第x个符号的通道脉冲响应hx包括2N个点,第x个符号的第n个取样时间的通道脉冲响应取样点hx,n表示为:
h x , n = IFFT ( H x , k ) = 1 2 N × Σ k = 0 2 N - 1 H x , k · e j 2 πnk 2 N
其中,n为0至2N-1的整数。
在一个符号内取四种可能位置的离散导频,通过上述的步骤可以得到四种不同的通道脉冲响应,然而仅有正确位置的离散导频所得到的通道脉冲响应可以完整且明显地看出通道的多径信息,而其他错误位置的离散导频所得到的通道脉冲响应看起来则像是白色噪声(white noise)。
请回去参照图1,假设对第3个符号进行离散导频同步以获得离散导频的正确位置,则本发明的示范实施例所提出的基于时域的离散导频同步方法会取第3个符号频域上的第12、24、...、12×pmax+12上的子载波来计算通道脉冲响应
Figure GDA0000131549670000092
同理,本发明的示范实施例所提出的基于时域的离散导频同步方法亦会同时取第3个符号频域上的第3、15、...、12×pmax+3上的子载波来计算通道脉冲响应
Figure GDA0000131549670000093
取第3个符号频域上的第6、18、...、12×pmax+6上的子载波来计算通道脉冲响应与取第3个符号频域上的第9、21、...、12×pmax+9上的子载波来计算通道脉冲响应
Figure GDA0000131549670000095
其中,通道脉冲响应
Figure GDA0000131549670000096
(i等于1至4的整数)的第n个取样时间的通道脉冲响应取样点
Figure GDA0000131549670000097
的计算方式如同前面所述。
接着,请参照图5A~图5D,图5A~图5D分别是通道脉冲响应
Figure GDA0000131549670000101
Figure GDA0000131549670000102
Figure GDA0000131549670000103
的通道脉冲响应示意图。如同前面所述,仅有正确位置的离散导频所得到的通道脉冲响应可以完整且明显地看出通道的多径信息,而其他错误位置的离散导频所得到的通道脉冲响应看起来则像是白色噪声。因此,在这个例子中,仅有图5C的通道脉冲响应
Figure GDA0000131549670000104
可以看得到完整的多径信息,而图5A、图5B、图5D的通道响应
Figure GDA0000131549670000105
Figure GDA0000131549670000106
看起来则像是白色噪声。因此,可以知道图1的第3个符号的频散导频位于频域上的第6、18、30、...、12×pmax+6与Kmax个位置上。
因为,仅有正确位置的离散导频所得到的通道脉冲响应可以完整且明显地看出通道的多径信息,而其他错误位置的离散导频所得到的通道脉冲响应看起来则像是白色噪声,所以,本发明的示范实施例所提供的基于时域的离散导频同步方法可以比较通道脉冲响应
Figure GDA0000131549670000107
Figure GDA0000131549670000108
的方差Cov(1)、Cov(2)、Cov(3)与Cov(4),其中,方差最大者表示具有多径信息,而在此例中,Cov(1)最大。Cov(u)的计算方式可以表示如下:
Cov ( u ) = Σ n = 0 2 N - 1 ( h 3 , n u ) 2 - 1 2 N ( Σ n = 0 2 N - 1 h 3 , n u ) 2
其中,u为1至4的整数。除此之外,上述计算方差的公式仅是一种实施方式,并非用以限定本发明。
值得一提的是,虽然上述的实施方式是比较通道脉冲响应
Figure GDA00001315496700001010
Figure GDA00001315496700001011
的方差来获得离散导频的正确位置,但本发明并非用以限定于此种实施方式。因为正确位置的离散导频所获得的通道脉冲响应会具有多径信息,因此可以通过比较通道脉冲响应
Figure GDA00001315496700001012
Figure GDA00001315496700001013
的峰值平均功率比(Peak toAverage Power Ratio,简称为PAPR)PAPR(1)、PAPR(2)、PAPR(3)与PAPR(4)来获得离散导频的正确位置。PAPR(u)的计算方式如下:
PAPR ( u ) = max arg n [ ( h 3 , n u ) 2 ] / mean _ power ( h 3 u )
其中,表示通道脉冲响应
Figure GDA0000131549670000113
的最大峰值功率,而
Figure GDA0000131549670000114
通道脉冲响应
Figure GDA0000131549670000115
的平均功率。除此之外,上述计算PAPR的公式仅是一种实施方式,并非用以限定本发明。
除了计算方差与PAPR的方式之外,比较通道脉冲响应
Figure GDA0000131549670000117
中最大径的值也可以是一种实施方式。另外,比较通道脉冲响应
Figure GDA0000131549670000118
Figure GDA0000131549670000119
中超过某一门槛值(Threshold)的径数,或比较减去最大径之后的通道脉冲响应
Figure GDA00001315496700001110
Figure GDA00001315496700001111
的能量也可以是一种实施方式。
接着,请参照图6,图6是本发明的示范实施例所提供的基于时域的离散导频同步装置的电路方块图。基于时域的离散导频同步装置60配置于DVB系统的接收端,其本身接收快速傅立叶转换(Fast Fourier Transform,简称为FFT)电路所输出的符号。基于时域的离散导频同步装置60包括了第一通道脉冲响应估计电路61、第二通道脉冲响应估计电路62、第三通道脉冲响应估计电路63、第四通道脉冲响应估计电路64、位置决定电路65与符号寄存器66。其中,符号寄存器66耦接于DVB系统的接收端的FFT电路与第一至第四通道脉冲响应估计电路61~64,而第一至第四通道脉冲响应估计电路61~64则耦接于位置决定电路65。
符号寄存器66用以接收FFT电路的所输出的符号,此符号位频域上的符号。接着,第一通道脉冲响应估计电路61会粹取符号中第一位置集合的多个子载波信号,并根据这些子载波信号获得第一通道脉冲响应。其中,第一位置集合指的是频域上第12、24、...、12×pmax+12个位置。换言之,第一通道脉冲响应估计电路61利用符号的频域上第12、24、...、12×pmax+12个子载波信号来计算第一通道脉冲响应。
第二通道脉冲响应估计电路62会粹取符号中第二位置集合的多个子载波信号,并根据这些子载波信号获得第二通道脉冲响应。其中,第二位置集合指的是频域上第3、15、...、12×pmax+3个位置。换言之,第二通道脉冲响应估计电路62利用符号的频域上第3、15、...、12×pmax+3个子载波信号来计算第二通道脉冲响应。
第三通道脉冲响应估计电路63会粹取符号中第三位置集合的多个子载波信号,并根据这些子载波信号获得第三通道脉冲响应。其中,第三位置集合指的是频域上第6、18、...、12×pmax+6个位置。换言之,第三通道脉冲响应估计电路63利用符号的频域上第6、18、...、12×pmax+6个子载波信号来计算第三通道脉冲响应。
第四通道脉冲响应估计电路64会粹取符号中第四位置集合的多个子载波信号,并根据这些子载波信号获得第四通道脉冲响应。其中,第四位置集合指的是频域上第9、21、...、12×pmax+9个位置。换言之,第四通道脉冲响应估计电路64利用符号的频域上第9、21、...、12×pmax+9个子载波信号来计算第四通道脉冲响应。
接着,位置决定电路65根据第一至第四通道脉冲响应决定符号中的离散导频的正确位置。如同前面所述,因为利用正确位置的离散导频所计算出的通道脉冲响应会具有多径信息,因此,第一至第四通道脉冲响应中会有一个通道脉冲响应具有多径信息,而其他的通道脉冲响应则像是白色噪声。
第一通道脉冲响应估计电路61包括了第一萃取器611、通道频率响应计算电路612与IFFT电路613。其中,第一萃取器611耦接于符号寄存器66与通道频率响应计算电路612,而IFFT电路613则耦接于通道频率响应计算电路612与位置决定电路65。第一萃取器611用以萃取符号中第一位置集合的多个子载波信号,频率响应计算电路612则根据第一萃取器611所输出的多个子载波信号计算第一通道频率响应。接着,IFFT电路613对第一通道频率响应进行补零的动作,并对补零后的第一通道频率响应作IFFT,以产生第一通道脉冲响应。
第二通道脉冲响应估计电路62包括了第二萃取器621、通道频率响应计算电路622与IFFT电路623。其中,第二萃取器621耦接于符号寄存器66与通道频率响应计算电路622,而IFFT电路623则耦接于通道频率响应计算电路622与位置决定电路65。第二萃取器621用以萃取符号中第二位置集合的多个子载波信号,频率响应计算电路622则根据第二萃取器621所输出的多个子载波信号计算第二通道频率响应。接着,IFFT电路623对第二通道频率响应进行补零的动作,并对补零后的第二通道频率响应作IFFT,以产生第二通道脉冲响应。
第三通道脉冲响应估计电路63包括了第三萃取器631、通道频率响应计算电路632与IFFT电路633。其中,第三萃取器631耦接于符号寄存器66与通道频率响应计算电路632,而IFFT电路633则耦接于通道频率响应计算电路632与位置决定电路65。第三萃取器631用以萃取符号中第三位置集合的多个子载波信号,频率响应计算电路632则根据第三萃取器631所输出的多个子载波信号计算第三通道频率响应。接着,IFFT电路633对第三通道频率响应进行补零的动作,并对补零后的第三通道频率响应作IFFT,以产生第三通道脉冲响应。
第四通道脉冲响应估计电路64包括了第四萃取器641、通道频率响应计算电路642与IFFT电路643。其中,第四萃取器641耦接于符号寄存器66与通道频率响应计算电路642,而IFFT电路643则耦接于通道频率响应计算电路642与位置决定电路65。第一萃取器641用以萃取符号中第四位置集合的多个子载波信号,频率响应计算电路642则根据第四萃取器641所输出的多个子载波信号计算第四通道频率响应。接着,IFFT电路643对第四通道频率响应进行补零的动作,并对补零后的第四通道频率响应作IFFT,以产生第四通道脉冲响应。
位置决定电路65包括数据处理电路651~654与比较电路656。其中,数据处理电路651耦接于IFFT电路613与比较电路656,数据处理电路652耦接于IFFT电路623与比较电路656,数据处理电路653耦接于IFFT电路633与比较电路656,数据处理电路654耦接于IFFT电路643与比较电路656。
数据处理电路651~654分别用以处理第一至第四通道脉冲响应,并分别产生第一至第四时域数据。其中,第一至第四时域数据可以是根据第一至第四通道脉冲响应所计算出来的第一至第四方差、第一至第四PAPR、第一至第四最大径的值、大于门槛值的第一至第四径数或减去最大径后的第一至第四残余能量等。其中,第一至第四方差分别是第一至第四通道脉冲响应所对应的方差,第一至第四PAPR是第一至第四通道脉冲响应所对应的PAPR,第一至第四最大径的值分别是第一至第四通道脉冲响应中的最大径的值,大于门槛值的第一至第四径数分别是第一至第四通道脉冲响应中大于某一门槛值的径数,减去最大径后的第一至第四残余能量分别是第一至第四通道脉冲响应中减去最大径后所剩余的能量。
接着,比较电路656比较第一至第四时域数据,并根据第一至第四时域数据决定离散导频在符号中的正确位置。举例来说,若数据处理电路651~654用以根据第一至第四通道脉冲响应来计算其对应的第一至第四方差时,则比较电路656会自第一至第四方差找出最大方差,并决定最大方差所对应的位置为离散导频的位置。
值得一提的是,当第一至第四时域数据是第一至第四方差、第一至第四PAPR、第一至第四最大径的值或大于门槛值的第一至第四径数时,比较电路656是找出第一至第四时域数据中最大者;然而,当第一至第四时域数据是减去最大径后的第一至第四残余能量时,比较电路656是找出第一至第四时域数据中最小者。理论上来说,第一至第四通道脉冲响应的能量大约相等,而具有多径信息的那一个通道脉冲响应的能量会集中在其最大径,因此,具有多径信息的通道脉冲响应在减去最大径后的残余能量最小。当然,多径信息有可能包含多个径,因此比较减去最大径后的第一至第四残余能量也可以修改成比较减去超过某一门槛值的径后的第一至第四残留能量。
另外,上述的实施方式虽然是以并行地方是来计算第一至第四通道脉冲响应,但事实上,基于时域的离散导频装置60亦可以仅有一个通道脉冲响应计算电路,此通道脉冲响应计算电路会分别萃取第一至第四位置集合的多个子载波,并计算出第一至第四通道脉冲响应。同理,基于时域的离散导频装置60亦可以仅有一个数据处理电路,此数据处理电路会分别根据第一至第四通道脉冲响应计算出第一至第四时域数据。
接着,请参照图7,图7是本发明的示范实施例所提供的基于时域的离散导频同步方法的流程图。首先,在步骤S70,利用符号寄存器接收来自于FFT电路的符号,此符号为频域上的符号。在步骤S71中,利用第一通道脉冲响应估计电路萃取符号的频域上的第一位置集合的多个子载波信号,并分别根据第一位置集合的多个子载波信号获得第一通道脉冲响应。步骤S71获得第一通道脉冲响应的详细方式是先根据第一位置集合的多个子载波信号计算出第一通道频率响应,接着,再将第一通道频率响应补零,并将补零后的第一通道频率响应作IFFT,以获得第一通道脉冲响应。
在步骤S72中,利用第二通道脉冲响应估计电路萃取符号的频域上的第二位置集合的多个子载波信号,并分别根据第二位置集合的多个子载波信号获得第二通道脉冲响应。步骤S72获得第二通道脉冲响应的详细方式是先根据第二位置集合的多个子载波信号计算出第二通道频率响应,接着,再将第二通道频率响应补零,并将补零后的第二通道频率响应作IFFT,以获得第二通道脉冲响应。
在步骤S73中,利用第三通道脉冲响应估计电路萃取符号的频域上的第三位置集合的多个子载波信号,并分别根据第三位置集合的多个子载波信号获得第三通道脉冲响应。步骤S73获得第三通道脉冲响应的详细方式是先根据第三位置集合的多个子载波信号计算出第三通道频率响应,接着,再将第三通道频率响应补零,并将补零后的第三通道频率响应作IFFT,以获得第三通道脉冲响应。
在步骤S74中,利用第四通道脉冲响应估计电路萃取符号的频域上的第四位置集合的多个子载波信号,并分别根据第四位置集合的多个子载波信号获得第四通道脉冲响应。步骤S74获得第四通道脉冲响应的详细方式是先根据第四位置集合的多个子载波信号计算出第四通道频率响应,接着,再将第四通道频率响应补零,并将补零后的第四通道频率响应作IFFT,以获得第四通道脉冲响应。
接着,在步骤S75中,使用位置决定电路根据第一至第四通道脉冲响应决定符号的离散导频的正确位置。步骤S75决定符号的离散导频的正确位置的详细方式是先根据第一至第四通道脉冲响应分别计算出第一至第四时域数据,接着,再比较第一至第四时域数据,以决定符号的离散导频的正确位置。其中,第一至第四时域数据可以是第一至第四方差、第一至第四PAPR、第一至第四最大径的值、大于门槛值的第一至第四径数、减去最大径后的第一至第四残余能量或减去大于门槛值的径后的第一至第四残留能量。值得一提的是,当第一至第四时域数据是第一至第四方差、第一至第四PAPR、第一至第四最大径的值或大于门槛值的第一至第四径数时,步骤S75是找出第一至第四时域数据中最大者,并决定最大者所对应的位置为符号的离散导频的正确位置;然而,当第一至第四时域数据是减去最大径后的第一至第四残余能量或减去大于门槛值的径后的第一至第四残留能量时,步骤S75是找出第一至第四时域数据中最小者,并决定最小者所对应的位置为符号的离散导频的正确位置。
值得一提的是,虽然上述的基于时域的离散导频同步方法是以数字广播电视系统的例子来说明。然而,示范实施例所提供的基于时域的离散导频同步方法亦并不受限于数字广播电视系统,换言之,本发明的示范实施例所提供的基于时域的离散导频同步方法、装置可以用于任何需要进行离散导频同步的通信系统中。
请参照图8,图8是中国移动多媒体广播系统(China Mobile MultimediaBroadcasting,CMMB)中的符号的离散导频在频域上的位置的示意图。在此CMMB系统中第2×i+1个符号的离散导频的位置为8×p1+1与8×p2+3,而第2×i+2个符号的离散导频的位置为8×p1+5与8×p2+7,其中,i为大于等于0的正整数。在2MHz的带宽模式下,Nv等于627,p1为0至38的整数,p2为39至77的整数;在8MHz的带宽模式下,Nv等于3075,p1为0至191的整数,p2为192至383的整数。
由图8及上述的叙述可知,在CMMB系统中不管是哪种带宽模式,离散导频都有两种可能的位置,修正上述示范实施例所提供的基于时域的离散导频同步方法与装置,使其在一个符号内分别提取这两种可能位置上的离散导频,来求得两个通道频率响应。接着,将此两个通道频率响应通过IFFT计算得到两个通道脉冲响应,再通过计算两个通道脉冲响应的一些时域特性,便能判断得到符号内正确的离散导频位置。简单地说,只要通信系统中的符号内的离散导频分布是具有规则性的,则本发明的示范实施例所提供的基于时域的离散导频同步方法与装置皆可以适用。
接着,请参照图9,图9是基于能量的离散导频同步方法、基于相关性的离散导频同步方法、基于时频相关性的离散导频同步方法与基于时域的离散导频同步方法在AWGN通道下的效能示意图。基于能量的离散导频同步方法的保护比定义为正确位置所对应的总能量与错误位置所对应的总能量的最大者的比值,基于相关性的离散导频同步方法的保护比定义正确位置所对应的相关性与错误位置所对应的相关性的最大者的比值,基于时频相关性的离散导频同步方法的保护比定义正确位置所对应的时频相关性与错误位置所对应的时频相关性的最大者的比值。载波噪声比的定义为载波的功率与噪声的功率比值。
在第一至第四时域数据是第一至第四方差、第一至第四PAPR、第一至第四最大径的值或大于门槛值的第一至第四径数时,基于时域的离散导频同步方法的保护比定义为正确位置所对应的时域数据与错误位置所对应的时域数据的最大者的比值;当第一至第四时域数据是减去最大径后的第一至第四残余能量或减去大于门槛值的径后的第一至第四残留能量时,基于时域的离散导频同步方法的保护比定义为错误位置所对应的时域数据的最小者与正确位置所对应的时域数据的比值。
在图9中,第一至第四时域数据是第一至第四方差,因此,基于时域的离散导频同步方法的保护比定义为正确位置所对应的时域数据与错误位置所对应的时域数据的最大者的比值。曲线801表示基于时域的离散导频同步方法的保护比与载波噪声比的对应曲线,曲线802表示基于能量的离散导频同步方法的保护比与载波噪声比的对应曲线,曲线803表示基于相关性的离散导频同步方法的保护比与载波噪声比的对应曲线,曲线804表示基于时频相关性的离散导频同步方法的保护比与载波噪声比的对应曲线。由图9可知,在载波噪声比很低的情况下(小于0dB),基于能量的离散导频同步方法的保护比会低于1,而无法正确地找出离散导频的位置。
接着,请参照图10,图10是基于能量的离散导频同步方法、基于相关性的离散导频同步方法、基于时频相关性的离散导频同步方法与基于时域的离散导频同步方法在时变通道下的效能示意图。时变通道是TU6通道,共定义了六条路径,第一至第六条路径的时间延迟分别是0、0.2、0.5、1.6、2.3与5.0微秒,第一至第六条路径的通道衰减值分别是3、0、2、6、8与10dB,而且第一至第六条路径的多普勒(Doppler)频谱皆采用瑞利(Rayleigh)分布模型。各离散频导方法的保护比的定义与图9相同,曲线901表示基于时域的离散导频同步方法的保护比与多普勒频率的对应曲线,曲线902表示基于能量的离散导频同步方法的保护比与多普勒频率的对应曲线,曲线903表示基于相关性的离散导频同步方法的保护比与多普勒频率的对应曲线,曲线904表示基于时频相关性的离散导频同步方法的保护比与多普勒频率对应曲线。由图10可知,在多普勒频率较高的情况下(大于400Hz),基于相关性的离散导频同步方法的保护比会低于1,而无法正确地找出离散导频的位置。
接着,请参照图11,图11是基于能量的离散导频同步方法、基于相关性的离散导频同步方法、基于时频相关性的离散导频同步方法与基于时域的离散导频同步方法在多径通道下的效能示意图。多径通道是静态的双径通道,第一与第二条路径的时间延迟分别是0与35微秒,第一与第二条路径的通道衰减值皆为0dB。各离散频导方法的保护比的定义与图9相同,曲线1001表示基于时域的离散导频同步方法的保护比与载波噪声比的对应曲线,曲线1002表示基于能量的离散导频同步方法的保护比与载波噪声比的对应曲线,曲线1003表示基于相关性的离散导频同步方法的保护比与载波噪声比的对应曲线,曲线1004表示基于时频相关性的离散导频同步方法的保护比与载波噪声比的对应曲线。由图11可知,在载波噪声比很低的情况下(小于5dB),基于时频相关性的离散导频同步方法的保护比会低于1,而无法正确地找出离散导频的位置。
综上所述,本发明的示范实施例所提供的基于时域的离散导频同步方法、装置及其DVB系统在各种不同的通道环境皆能够具有良好的效能,并能够准确地检测出离散导频的正确位置。
虽然本发明已以实施例揭露如上,然其并非用以限定本发明,任何所属技术领域的技术人员,在不脱离本发明的精神和范围内,当可作些许的更动与润饰,故本发明的保护范围当以权利要求所界定范围为准。

Claims (15)

1.一种离散导频同步方法,其特征在于,所述离散导频同步方法用于一通信系统,所述方法包括:
接收一频域符号;
分别萃取所述频域符号的第一至第x位置集合的多个子载波,其中,x为大于1的整数;
获得第一至第x通道脉冲响应,其中,第i通道脉冲响应是根据第i位置的子载波所获得,i为1至x的整数;以及
根据所述第一至第x通道脉冲响应决定所述频域符号的多个离散导频的位置。
2.如权利要求1所述的离散导频同步方法,其特征在于,所述第一至第x通道脉冲响应中具有多径信息的通道脉冲响应所对应的多个位置为所述频域符号的所述这些离散导频的位置。
3.如权利要求1所述的离散导频同步方法,其特征在于,获得所述第一至第x通道脉冲响应的步骤更包括:
获得第一至第x通道频率响应,其中,所述第i通道频率响应是根据所述第i位置的子载波计算而获得;
分别对所述第一至第x通道频率响应进行补零的动作;
分别将补零后的所述第一至第x通道频率响应进行反快速傅立叶转换,以获得所述第一至第x通道脉冲响应。
4.如权利要求1所述的离散导频同步方法,其特征在于,决定所述频域符号的所述这些离散导频的位置的步骤更包括:
分别对所述第一至第x通道脉冲响应进行处理,以获得第一至第x时域数据;以及
比较所述第一至所述第x时域数据,以决定所述频域符号的所述这些离散导频的正确位置。
5.如权利要求4所述的离散导频同步方法,其特征在于,所述第一至第x时域数据是第一至第x方差、第一至第x PAPR、第一至第x最大径的值、大于一门槛值的第一至第x径数、减去一最大径后的第一至第x残余能量或减去大于所述门槛值的径后的第一至第x残留能量。
6.一种离散导频同步装置,其特征在于,所述离散导频同步装置用于一通信系统系统,所述装置包括:
一符号寄存器,用以接收一频域符号;
至少一通道脉冲响应计算电路,耦接于所述符号寄存器,用以萃取所述频域符号的第一至第x位置集合的多个子载波,并用以获得第一至第x通道脉冲响应,其中,x为大于1的整数,第i通道脉冲响应是根据第i位置的子载波计算得来,i为1至x的整数;以及
一位置决定电路,耦接于所述通道脉冲响应计算电路,根据所述第一至第x通道脉冲响应决定所述频域符号的多个离散导频的位置。
7.如权利要求6所述的离散导频同步装置,其特征在于,位置决定电路将所述第一至第x通道脉冲响应中具有多径信息的通道脉冲响应所对应的多个位置决定为所述频域符号的所述这些离散导频的位置。
8.如权利要求6所述的离散导频同步装置,其特征在于,所述通道脉冲响应计算电路包括:
一第i萃取器,耦接于所述符号寄存器,用以萃取所述频域符号的所述第i位置集合的所述这些子载波;
一通道频率响应计算电路,耦接于所述第i萃取器,用以根据所述第i位置的子载波计算一第i通道频率响应;
一IFFT电路,用以对所述第i通道频率响应进行补零的动作,并将补零后的所述第i通道频率响应进行IFFT,以获得所述第i通道脉冲响应。
9.如权利要求6所述的离散导频同步装置,其特征在于,所述位置决定电路包括:
至少一数据处理电路,每一所述至少一数据处理电路耦接于其所对应的通道脉冲响应计算电路,且每一所述至少一数据处理电路用以分别对所对应的所述第一至第x通道脉冲响应中的一者进行处理,从而所述至少一数据处理电路获得第一至第x时域数据;以及
一比较电路,耦接于所述数据处理电路,用以比较所述第一至所述第x时域数据,以决定所述频域符号的所述这些离散导频的正确位置。
10.如权利要求9所述的离散导频同步装置,其特征在于,所述第一至第x时域数据是第一至第x方差、第一至第x PAPR、第一至第x最大径的值、大于一门槛值的第一至第x径数、减去一最大径后的第一至第x残余能量或减去大于所述门槛值的径后的第一至第x残留能量。
11.一种通信系统,其特征在于,所述通信系统包括一接收端以及一发射端,所述接收端包括一离散导频同步装置,所述离散导频同步装置包括:
一符号寄存器,用以接收一频域符号;
至少一通道脉冲响应计算电路,耦接于所述符号寄存器,用以萃取所述频域符号的第一至第x位置集合的多个子载波,并用以获得第一至第x通道脉冲响应,其中,x为大于1的整数,第i通道脉冲响应是根据第i位置的子载波计算得来,i为1至x的整数;以及
一位置决定电路,耦接于所述通道脉冲响应计算电路,根据所述第一至第x通道脉冲响应决定所述频域符号的多个离散导频的位置。
12.如权利要求11所述的通信系统,其特征在于,位置决定电路将所述第一至第x通道脉冲响应中具有多径信息的通道脉冲响应所对应的多个位置决定为所述频域符号的所述这些离散导频的位置。
13.如权利要求11所述的通信系统,其特征在于,所述通道脉冲响应计算电路包括:
一第i萃取器,耦接于所述符号寄存器,用以萃取所述频域符号的第i位置的子载波;
一通道频率响应计算电路,耦接于所述第i萃取器,用以根据所述第i位置的子载波计算一第i通道频率响应;
一IFFT电路,用以对所述第i通道频率响应进行补零的动作,并将补零后的所述第i通道频率响应进行IFFT,以获得所述第i通道脉冲响应。
14.如权利要求11所述的通信系统,其特征在于,所述位置决定电路包括:
至少一数据处理电路,每一所述至少一数据处理电路耦接于其所对应的通道脉冲响应计算电路,且每一所述至少一数据处理电路用以分别对所对应的所述第一至第x通道脉冲响应中的一者进行处理,从而所述至少一数据处理电路获得第一至第x时域数据;以及
一比较电路,耦接于所述数据处理电路,用以比较所述第一至所述第x时域数据,以决定所述频域符号的所述这些离散导频的正确位置。
15.如权利要求14所述的通信系统,其特征在于,所述第一至第x时域数据是第一至第x方差、第一至第x PAPR、第一至第x最大径的值、大于一门槛值的第一至第x径数、减去一最大径后的第一至第x残余能量或减去大于所述门槛值的径后的第一至第x残留能量。
CN 200910006341 2009-02-10 2009-02-10 离散导频同步方法与装置及其通信系统 Expired - Fee Related CN101800719B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200910006341 CN101800719B (zh) 2009-02-10 2009-02-10 离散导频同步方法与装置及其通信系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200910006341 CN101800719B (zh) 2009-02-10 2009-02-10 离散导频同步方法与装置及其通信系统

Publications (2)

Publication Number Publication Date
CN101800719A CN101800719A (zh) 2010-08-11
CN101800719B true CN101800719B (zh) 2013-03-20

Family

ID=42596218

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200910006341 Expired - Fee Related CN101800719B (zh) 2009-02-10 2009-02-10 离散导频同步方法与装置及其通信系统

Country Status (1)

Country Link
CN (1) CN101800719B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103078819B (zh) * 2011-10-26 2015-09-30 扬智电子科技(上海)有限公司 精符号同步方法及其装置
CN103501284B (zh) * 2013-10-22 2016-09-28 湖南国科微电子股份有限公司 一种dvb-t2系统中精细定时同步的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060088133A1 (en) * 2004-10-22 2006-04-27 Industrial Technology Research Institute Time-frequency correlation-based synchronization for coherent OFDM receiver
CN1905546A (zh) * 2005-07-25 2007-01-31 上海原动力通信科技有限公司 时分双工ofdm系统确定导频符号位置的方法
WO2008045845A1 (en) * 2006-10-09 2008-04-17 Intel Corporation Link adaptation for retransmission error-control technique transmissions

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060088133A1 (en) * 2004-10-22 2006-04-27 Industrial Technology Research Institute Time-frequency correlation-based synchronization for coherent OFDM receiver
CN1905546A (zh) * 2005-07-25 2007-01-31 上海原动力通信科技有限公司 时分双工ofdm系统确定导频符号位置的方法
WO2008045845A1 (en) * 2006-10-09 2008-04-17 Intel Corporation Link adaptation for retransmission error-control technique transmissions

Also Published As

Publication number Publication date
CN101800719A (zh) 2010-08-11

Similar Documents

Publication Publication Date Title
CN101447970B (zh) 利用训练序列进行lofdm系统定时和载波同步的方法
US7961816B2 (en) Device for and method of signal synchronization in a communication system
CN101425999B (zh) 正交频分复用接收机的载频偏差同步的方法及装置
CN109660478A (zh) 一种基于改进的Park频域训练序列的定时频率同步方法
CN101809930A (zh) Dmb-t系统的时间频率同步和帧号检测
CN109743278A (zh) 物理帧中前导符号的接收处理方法
CN101588327B (zh) 用于ofdm系统的数据指导信道估计方法及其实现装置
CN101465833A (zh) 一种正交频分复用信号定时同步方法和装置
CN102130871A (zh) 信道估计方法及装置
CN102461102A (zh) Lte通信网络中根据pss信号的信道估计过程及其接收器
CN101605120A (zh) 用于正交频分复用系统的迭代内码译码和信道估计方法及装置
CN105187352A (zh) 一种基于ofdm前导的整数频偏估计方法
CN101606367B (zh) 使多载波系统同步的方法和多载波系统
CN104022996A (zh) 一种基于信道估计的ofdm系统定时同步方法
CN101552758B (zh) 正交频分复用系统中精符号时序同步方法
CN101257471B (zh) 移动多媒体数字广播ofdm系统信道估计方法
CN101800719B (zh) 离散导频同步方法与装置及其通信系统
CN101374129B (zh) 基于正交频分复用的同步序列生成方法、同步方法及系统
CN102265573B (zh) 用于估计正交频分复用传输系统中相位噪声的方法和装置
CN103236993B (zh) 一种基于多径延时分布的信道估计方法
CN102377701B (zh) 基于循环相关的信道估计方法及系统
CN103078819B (zh) 精符号同步方法及其装置
CN102065035B (zh) 多带正交频分复用超宽带系统的信道估计方法
CN101447969A (zh) 一种多带正交频分复用超宽带系统的信道估计方法
CN104735012B (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: 20130320

Termination date: 20190210

CF01 Termination of patent right due to non-payment of annual fee