CN101116298B - 用于信道均衡的方法和系统 - Google Patents

用于信道均衡的方法和系统 Download PDF

Info

Publication number
CN101116298B
CN101116298B CN2006800039607A CN200680003960A CN101116298B CN 101116298 B CN101116298 B CN 101116298B CN 2006800039607 A CN2006800039607 A CN 2006800039607A CN 200680003960 A CN200680003960 A CN 200680003960A CN 101116298 B CN101116298 B CN 101116298B
Authority
CN
China
Prior art keywords
signal
channel
characteristic
value
equipment
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.)
Active
Application number
CN2006800039607A
Other languages
English (en)
Other versions
CN101116298A (zh
Inventor
D·P·马拉迪
B-h·金
魏永斌
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.)
Qualcomm Inc
Original Assignee
Qualcomm Inc
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 Qualcomm Inc filed Critical Qualcomm Inc
Publication of CN101116298A publication Critical patent/CN101116298A/zh
Application granted granted Critical
Publication of CN101116298B publication Critical patent/CN101116298B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00Arrangements for detecting or preventing errors in the information received
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L12/00Data switching networks
    • H04L12/28Data switching networks characterised by path configuration, e.g. LAN [Local Area Networks] or WAN [Wide Area Networks]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L2025/0335Arrangements for removing intersymbol interference characterised by the type of transmission
    • H04L2025/03375Passband transmission
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L2025/03777Arrangements for removing intersymbol interference characterised by the signalling
    • H04L2025/03783Details of reference signals
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/022Channel estimation of frequency response
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0224Channel estimation using sounding signals
    • H04L25/0228Channel estimation using sounding signals with direct estimation from sounding signals
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

一个实施例包括一种接收发送信号的方法。该方法包括接收在一信道上发送的信号。该信号包含一已知信号和一信息信号。该方法还包括至少部分基于该部分的已知信号来确定信道特性的至少一个指标。该方法进一步包括至少部分基于信道特性的该至少一个指标来生成指示该信息信号的第一值。该第一值包含一误差信号。该方法进一步包括至少部分基于该部分的已知信号从该信号的第一估计消除该误差信号。其它实施例包括用于执行该方法的系统以及制造此类系统的方法。

Description

用于信道均衡的方法和系统
在U.S.C.§119下的优先权要求
本发明申请要求2005年2月4日提交、且被转让给本发明受让人并因此被明确地援引包含在本文中的题为“FREQUENCY DOMAIN EQUALIZATION IN THEPRESENCE OF HIGH DOPPLER(存在高多普勒效应情况下的频域均衡)”的临时申请No.60/650,042号的优先权。
背景
领域
本发明涉及无线通信方法和系统,尤其涉及频域均衡方法和系统。
背景
在具有较高数据率——例如,在50-200Mbps范围里的数据率——的通信系统中,接收信号中的误差可能导致阻碍可用带宽的充分利用的重发延迟。减少重发延迟的一种方式是减少接收信号中的比特差错率(BER)。
更具体地,在无线信道上通信期间,信道行为随时间变化,由此影响在信道上发送的信号(例如,增加其中的差错)。需要预测或表征随时间变化的信道行为以在接收发送信号时补偿信道特性的此类变动。在某些系统中,在时域中实现信道均衡以更准确地估计发送信号。但是,时域中的均衡要求接收机中有强大的计算能力,并使得电路复杂化。此外,诸如信道上的白噪声等噪声的存在进一步使得时域中的校正复杂化。因此,需要简化计算信道行为和校正在信道上发送的信号中的计算复杂度。
概述
本发明的一个方面包括一种处理信号的方法。该方法包括接收一在信道上发送的信号。该信号包含一信息信号以及至少一部分的已知信号。该方法进一步包括至少部分基于该部分已知信号来确定信道特性的至少一个指标。该方法进一步包括至少部分基于信道特性的该至少一个指标来生成指示该信息信号的第一值。该第一值包含一误差信号。该方法进一步包括至少部分基于该部分已知信号来从该信号的第一估计消除该误差信号。
本发明的另一个方面包括一种配置接收在信道上发送的信号的设备。该信号包含一信息信号以及至少一部分的已知信号。该设备包括第一电路,该第一电路配置成至少部分基于该部分的已知信号来确定信道特性的至少一个指标。该设备进一步包括一第二电路,该第二电路配置成至少部分基于信道特性的该至少一个指标来生成指示该信息信号的第一值。该第一值包含一误差信号。该设备进一步包括一第三电路,该第三电路被配置成至少部分基于该部分已知信号来从该信号的第一估计消除该误差信号。
本发明的另一个方面包括一种配置成接收在信道上发送的信号的设备。该信号包含一信息信号以及至少一部分的已知信号。该设备包括至少部分基于该部分已知信号来确定信道特性的至少一个指标的装置。该设备进一步包括至少部分基于信道特性的该至少一个指标来生成指示该信息信号的第一值的装置。该第一值包含一误差信号。该设备进一步包括至少部分基于该部分已知信号从该信号的第一估计消除该误差信号的装置。
附图简述
图1示出一种示例性无线通信系统的概观。
图2是示出用于在图1的系统的一个实施例中传送的一种示例性数据帧结构的框图。
图3是更详细地示出图1的示例性无线通信系统的框图。
图4是示出在图3所示的示例性系统中传送数据的方法的一个实施例的流程图。
图5是示出在诸如图4所示的方法的一部分中的一种估计信道特性的示例性方法的流程图。
图6是示出在诸如图4所示的方法的一部分中的一种执行频域均衡的示例性方法的流程图。
图7是示出在诸如图4所示的方法的一部分中的一种消除残余码元间干扰(ISI)的示例性方法的流程图。
具体说明
以下具体说明针对本发明的某些具体实施例。但是,本发明可用如所附权利要求定义和涵盖的众多不同方式实施。在此说明中对附图进行参考,贯穿所有附图,相似的部分被标以相同的标号。
在一个实施例中,一接收机接收在射频(RF)信道上发送的信号。这些信号包含例如导频信号的一已知信号、以及例如数据信号的一未知信号。该接收机至少部分基于该已知信号的至少一接收部分来对该通信信道的特性作出估计。该接收机至少部分基于所估计的信道特性来生成该未知信号的第一估计。该未知信号的第一估计包含例如残余码元间干扰(ISI)的一误差信号。该接收机至少部分基于该已知信号来确定该误差信号的估计。该接收机至少部分基于该误差信号估计来从该信号的第一估计消除该误差信号。
图1示出示例性无线通信系统100的概观。在该示例性实施例中,通信系统100包括一个或多个基站102以及一个或多个用户终端104。在该示例性实施例中,该通信系统被配置成作为一蜂窝无线电网络来工作。蜂窝无线电网络包含一个或多个基站102。每个基站102向不同区域(“小区”)(它们可能重叠)提供通信,以提供其区域大于一个小区的区域的无线电覆盖。用户终端104的位置可以是固定的或移动的。可以使用各种切换技术来允许移动着的用户终端104在诸如进入或穿过各小区时与不同基站102通信。在其它实施例中,通信系统100可包括用户终端102之间的点对点通信或者是单向通信链路。另外,本文中某些实施例是参考使用射频(RF)载波的无线通信来讨论的。但是,在其它实施例中,该通信网络可包括其它通信介质,诸如光信号或是有线连接上的通信等。
系统100的各种实施例可在一个或多个RF频带——例如,一般称为800MHz、850MHz、900MHz、1800MHz、1900MHz或2000MHz的频带——中的一个或多个信道上通信。通信系统100的实施例包括决定该通信系统如何操作基站102与用户终端104之间的无线电链路的空中接口。例如,通信系统100可使用基于码分多址(CDMA)的空中接口或是时分多址(TDMA)空中接口。在一个示例性实施例中,通信系统100包括使用1900MHz频带中的一5MHz信道的宽带CDMA(W-CDMA)空中接口。通常,此5MHz频带中仅3.84MHz可供使用。
某些类型的数据传输可能对重发延迟很敏感,例如语音和诸如TCP(传输控制协议)等常用因特网协议。例如,如果传输时间(包括在较低通信层中重发丢失的数据的时间)太大,则TCP连接一般没有充分利用可用信道带宽。例如,在一个实施例中,可实现在50Mbps范围里的数据率,其TCP SER(分段差错率)在2×10-5到5×10-6范围里。在一个实施例中,通过对接收信号执行信道均衡,通信系统100的比特差错率(BER)得以降低,由此最终诸如TCP等较高层协议的分段差错率也得以降低。
信道是指可在其上传输信号的通信介质。一般而言,信道不是理想的,例如信道一般具有影响在该信道上传送的信号的时间和/或频率相关特性。在数学上,信道可由将输入——例如发送——到该信道的信号与从该信道输出——例如,接收——的信号相关的信道冲激响应h(t)表示或表征。信道均衡一般是指响应于在其上发送信号的信道的动态特性来调整接收信号的过程。
均衡可在时域或频域中执行。但是,时域均衡器可能计算复杂度高。另外,有限长度时域均衡器在某些信道状况下易有性能损失。在一个实施例中,本发明提供一种包含计算效率高且能补偿多普勒效应的频域均衡器的设备。
在以下将进一步详细阐述的一个实施例中,系统100通过将一已知信号随一数据信号传送给接收机来执行信道均衡。该数据信号可包含一种或多种形式的数据,诸如语音或其它数据。例如,该数据信号可包括至少部分在IP(网际协议)网络上传送的TCP数据。该接收机标识该已知信号,并利用其来估计信道的特性。此信道估计被用于在频域中执行信道均衡。
码元间干扰(ISI)一般是指在一信号的脉冲(对应于数据中的不同码元)之间的干扰,在相邻脉冲在时间上弥散从而相互重叠的情况下会发生这种干扰。当此重叠变得太大时,接收机可能不再能够正确地区别或标识每个脉冲。在一个实施例中,经均衡信号可能包含尚未被均衡消除的码元间干扰(ISI)。在一个此类实施例中,通过使用该已知信号与该已知信号的接收版本之差估计此残余ISI来消除此残余ISI。
图2是示出用于在通信系统100的一个实施例——例如,高速下行链路分组接入(HSDPA)系统——中传送的一种示例性数据帧结构的框图。HSDPA是宽带码分多址(W-CDMA)系统中的一种基于分组的数据服务,其中数据在时隙110期间被传送给接收机。在一个实施例中,该数据包括表示诸如语音、图像、视频或其它数据的用户数据的一个或多个比特的“码片”。例如,可使用前向纠错码将用户数据比特编码成若干码元。通过应用“扩展”码,这些码元中的每一个可被进一步编码成更大数目的码片。在诸如CDMA等空中接口中使用扩展码来创建扩频信号,例如在显著宽于用户数据的频率内容的频带中传输的信号。
在图2的示例性结构中,时隙110包括各有1280个码片的两个子隙。每个子隙包括三个数据块112,两个为448个码片,一个为128个码片。示例性时隙110还包括分隔这些数据块112的64比特前缀114和64比特后缀116。在前缀114和后缀116中传送一已知或导频信号xp(t),从而发射机302在每一时隙内传送这些已知信号xp(t)和数据信号xd(t)。在一个实施例中,该已知信号在每一时隙中是相同的。例如,该已知信号可被选为跨整个带宽具有恒定频率响应。这样一个序列可包括多相序列, x p [ n ] = e - | πn 2 / P ∀ 0 ≤ n ≤ P - 1 , 其中P是表示该已知或导频序列的长度的有限整数,例如在图2中所示的实施例中为64。在一个实施例中,该导频序列包括无频谱零点的信号,例如,Xp(f)(xp(t)在频率空间中的形式)被选为异于0。
在另一个实施例中,已知信号xp(t)是在一不同码信道上传送的信号,例如CDMA系统的导频信道信号。就在数据信号之前、之后或与之并发接收的导频信号的一部分可被用作该已知信号。在一个实施例中,CDMA系统的导频信道被以伪随机数(PN)序列编码。由此,系统100的这样一个实施例中的该已知信号xp(t)是时变的,并且依赖于在该PN序列内的位置。在其它实施例中,导频信号可包括任何已知的信号。
图3是更详细地示出示例性无线通信系统100的框图。在各种实施例中,图3中所示的功能块可由执行软件指令的处理器实现,被实现为数字电路,被实现为模拟电路,或由其组合实现。此框图示出系统100中的发射机302和接收机304两者的一部分。特别地,在一个实施例中,发射机302可包括块122-124,而接收机304可包括块130-146。在一个实施例中,接收机304包括使用半导体制造技术形成的一个或多个集成电路。这些集成电路可包括通用处理器、专用集成电路(ASIC)、数字信号处理器(DSP)、其它合适硬件或其组合当中的至少一个。在一个实施例中,块130-146中的两个或多个被形成在单块集成电路上。在一个实施例中,发射机302包括基站102,而接收机304包括用户终端104。在一个实施例中,系统100中的通信是非对称的,例如,仅基站102包括图3的发射机302部分,并且仅用户终端104包括接收机304部分。在其它实施例中,系统100中的通信是对称的,基站102和用户终端104两者皆包括图3的发射和接收机功能块两者。
在一个实施例中,发射机302包括将已知信号xp(t)与数据信号xd(t)组合以形成信号x(t)的组合器122。在一个实施例中,组合器122将已知信号xp(t)与数据信号xd(t)组合以形成如图2的时隙110中所示的数据流。在另一CDMA系统实施例中,组合器122通过依从CDMA系统将已知信号xp(t)置于例如导频信道的一个码信道中,并将数据信号xd(t)置于例如数据信道的另一个码信道中来将已知信号xp(t)与数据信号xd(t)组合。
发射机302进一步包括脉冲整形器124,脉冲整形器124通过有效地将脉冲整形信号p(t)与信号x(t)卷积来生成带宽受限的信号。该带宽受限信号在由信道冲激响应h(t)表示、并且在图3中被示为块126的信道上被传送。信道的作用可被表示为信道冲激响应h(t)与经整形的信号x(t)*p(t)卷积,例如,x(t)*p(t)*h(t)。
在一个实施例中,在接收机上,对该信号应用匹配滤波器130(与脉冲整形块124匹配)来接收发送的信号。匹配滤波器130对接收自该信道的信号x(t)*p(t)*h(t)应用该脉冲整形信号的复共轭p*(-t)来产生接收信号y(t)。在数学上,从匹配滤波器130输出的接收信号y(t)可被表示为:
y(t)=x(t)*p(t)*h(t)*p*(-t)。(式1)
接收信号y(t)包含对应于发送数据信号xd(t)和发送已知信号xp(t)的接收数据信号yd(t)和接收已知信号yp(t)。
接收机304进一步包括快速傅立叶变换(FFT)模块132,该模块132将接收信号y(t)从时域变换成频域信号Y(f)。通过将脉冲整形和匹配滤波视为构成一有效传输信道,频域中的接收信号可被表示为发送信号X(f)与信道的传递函数H(f)之积:
Y(f)=X(f)·H(f),其中“·”表示相乘。(式2)
在一个实施例中,是基于信道特性的估计
Figure G06803960720070807D000061
来获得发送信号的估计
Figure G06803960720070807D000062
特别地,接收机304包括基于分别与发送已知信号xP(t)和XP(F)对应的这部分接收已知信号yP(t)和YP(F)来确定信道特性的估计
Figure G06803960720070807D000063
的信道估计器134。例如,在一个此类实施例中,已知信号xP(t)和该已知信号的接收版本yP(t)可在离散频率空间中分别被表示为XP[k]和YP[k]。类似地,信道的传递函数H(f)在离散频率空间中可被表示为H[k]。假定系统中存在表示为W[k]的加性白噪声,并且假定该已知信号XP[k]包括具有整数长度的P个码片或脉冲的信号,则接收已知信号YP[k]可用离散形式被表达为:
Y P [ k ] = P · H [ k ] · X P [ k ] + W [ k ] . (式3)
可使用该已知发送信号(在离散频率空间中)XP[k]与已知接收信号(在离散频率空间中)YP[k]之间的关系来求解此式以确定信道特性
Figure G06803960720070807D000065
的估计。求解该式的一种方法包括应用称为最小均方估计(MMSE)技术的统计解。在一个实施例中,可在离散频率空间中得到信道估计
Figure G06803960720070807D000071
的MMSE解如下:
H ^ ( k ) = R cc [ k ] · P · X P * [ k ] P · R cc [ k ] · | X P [ k ] | 2 + R ww [ k ] · Y [ k ] ,
其中平均信道统计表示为Rcc=E{|H[k]|2}(式4)
而平均噪声方差表示为Rww[k]=E{|W[k]|2}。
在一个实施例中,信道和噪声方差统计Rcc和Rww可从先前接收的数据推导。例如,Rcc可从
Figure G06803960720070807D000073
的长期平均——诸如在若干帧上的平均来计算。Rww可基于导频样本中的误差来估计。在另一个实施例中,可得到迫零(ZF)(其中噪声方差被取为其零极限)MMSE离散解如下:
H ^ [ k ] = 1 P 1 X P [ k ] Y P [ k ] (式5)
使用信道估计器134的输出,接收机304进一步包括均衡器136,均衡器136使用该信道估计
Figure G06803960720070807D000075
来对接收信号Y(f)进行滤波以得到发送数据信号的估计
Figure G06803960720070807D000076
如以上提及的,接收信号Y(f)等于发送信号X(f)与信道传递函数H(F)之积。由此,通过取接收信号Y(f)与信道传递函数的复共轭H*(f)之积可得到发送信号x(t)的频域估计
Figure G06803960720070807D000077
考虑到加性白噪声W(f)的存在,估计的发送信号可被表示为:
X ^ ( f ) = X ( f ) · H ( f ) · H * ( f ) + W ( f ) = Y ( f ) · H * ( f ) + W ( f ) . (式6)
接收信号的这一估计的最小均方估计(MMSE)均衡器解
Figure G06803960720070807D000079
可被表示为:
X ^ ( f ) = Y ( f ) · H ^ * ( f ) | H ^ ( f ) | 2 + σ 2 , (式7)
其中σ2表示与该噪声项W(f)相关的噪声方差项。
MMSE解中的噪声方差项σ2起到偏项的作用,以消去可能会在此式中引入不连续的零点。在一个实施例中,接收信号
Figure G06803960720070807D0000711
的以上表示的离散频率空间MMSE解可被表达如下:
X ^ ( k ) = R bb [ k ] · P · H ^ * ( k ) P · R bb [ k ] · | H ^ [ k ] | 2 + R ww [ k ] · Y [ k ] ,
其中Rbb[k]=E{|X[k]|2}并且Rww[k]=E{|W[k]|2}并且(式8)
其中W表示数据信号中的噪声。
在这一实施例中,数据统计Rbb和噪声统计Rww可从估计接收数据的信号能量或从先前接收的数据来推导。Rbb可从接收数据信号星座能量或从话务或数据能量与导频能量之比来推导。Rww可基于导频样本中的误差来估计。
如以上所提及的,估计的信号包含执行快速傅立叶逆变换之后的残余ISI。在一个实施例中,此残余ISI被校正以得到发送数据信号xd(t)的估计
Figure G06803960720070807D000081
特别地,接收机304进一步包括快速傅立叶逆变换(IFFT)模块140,该模块140将经均衡信号
Figure G06803960720070807D000082
转换成时域信号
Figure G06803960720070807D000083
接收机304进一步包括ISI校正模块142,该模块142对残余ISI进行校正以推导数据信号xd(t)的估计
Figure G06803960720070807D000084
在参考图7进一步详细讨论的一个实施例中,ISI校正模块142基于已知信号与该已知信号的接收版本——例如xP(t)与
Figure G06803960720070807D000085
——之差对残余ISI进行校正。接收机304进一步包括对估计的数据信号
Figure G06803960720070807D000086
作进一步处理的解调器144。此处理可包括例如使用诸如QPSK、16-QAM、64-QAM或任何其它合适的方案的解调方案来从调制信号恢复该信号。
图4是示出在图3的示例性系统100中传送数据的方法200的一个实施例的流程图。方法200始于框202,在此例如基站102的发射机标识要进行传输的数据信号xD(t)。在一个实施例中,该数据信号至少被部分地调制,例如,该数据信号xD(t)包括已用一个或多个纠错码编码的诸如声音、视频、语音或其它数据的数据。移到框204,发射机302生成一已知或导频信号xP(t)。在一个实施例中,该导频信号xP(t)包括作为前缀(例如,在前)或后缀(例如,在后)插入到要在一时隙期间发送的数据块的数据脉冲、码片。在一个实施例中,该导频信号的长度为64个码片。
前进至框206,发射机302在信道上将该数据信号xD(t)和导频信号xP(t)传送给例如用户终端104的接收机作为接收信号y(t)=yD(t)+yp(t)。接下来在框210,接收机304执行快速傅立叶逆变换以将该接收信号y(t)从时域转换到频域中成为Y(f)。
移到框212,接收机304基于至少一部分的该已知信号XP(f)和相应的Yp(f)来推导频域中的信道特性估计
Figure G06803960720070807D000087
此估计过程将参考图5更详细地说明。接下来在框214,接收机304基于该信道估计
Figure G06803960720070807D000088
生成接收数据信号的包含残余ISI的估计
Figure G06803960720070807D000089
此过程将参考图6更详细地说明。
移到框216,接收机304执行快速傅立叶逆变换以将该频域接收信号Y(t)转换到时域信号y(t)。接下来在框220,接收机304校正残余ISI以产生数据信号xd(t)的估计信号
Figure G06803960720070807D0000810
在一个实施例中,接收机304基于至少一部分该已知信号xP(t)与对应于该已知信号的接收信号yP(t)之差来校正残余ISI。此过程将参考图7来更详细地说明。前进至框222,接收机304通过对纠错或其它编码进行解码来解调接收数据信号
Figure G06803960720070807D0000811
以提供原始声音、视频、语音或其它数据。
图5是示出在诸如图4所示的方法200的一部分中的一种估计信道特性
Figure G06803960720070807D0000812
的示例性方法212的流程图。在框250,接收机304标识导频信号yp(t)的P个样本,并执行FFT以将那些样本转换到频域,例如Yp(f)。接下来在框252,接收机304如以上参考图3的信道估计器134所述地基于导频信号Yp(f)来推导信道估计
Figure G06803960720070807D000091
在一个实施例中,可使用式4或5来推导
Figure G06803960720070807D000092
的离散形式
一般而言,数据块要比导频样本的数目P长得多。例如,在图2的示例性实施例中,数据块包含D=448或128个码片,而前缀或后缀块仅为P=64个码片长。由此,信道估计
Figure G06803960720070807D000094
的离散数据点或称抽头(tap)的数目与接收信号Y(f)的离散数据点或称抽头的数目之间存在着失配。由此,在将参考框254-258进一步说明的一个实施例中,P-抽头(例如,有P个离散值)的频域信道估计被转换到D-抽头(例如,具有D个离散值)的频域信道估计以便于在图4的均衡块214中使用。特别地,在图5的框254,接收机304执行IFFT以将信道估计
Figure G06803960720070807D000096
转换到时域成为
Figure G06803960720070807D000097
接下来在框256,接收机304使用合适的阈值来执行噪声抑制(迫零)。在一个实施例中,落在阈值——例如,数据中的最大值的10-30%——之下的数据点被迫零或设为零,或从信道估计中被移除。移到框258,接收机304执行FFT以将信道估计
Figure G06803960720070807D000098
转换到D-抽头(例如,在一个实施例中有D值=448)频域信道估计该信道估计
Figure G06803960720070807D0000910
由此可供进一步处理接收信号使用。
图6是示出在诸如图4所示的方法200的一部分中的一种执行频域均衡的示例性方法214的流程图。在框270,接收机304标识接收信号X(f)的一部分。在一个实施例中,该估计的信号X(f)包含一部分的导频信号XP(f)。所包括的这一部分导频信号由此如参考图6所讨论地被均衡以便于在后面的处理步骤中使用。接下来在框272,接收机304基于信道估计
Figure G06803960720070807D0000911
推导发送数据信号的包含残余ISI的估计
Figure G06803960720070807D0000912
例如
Figure G06803960720070807D0000913
在一个实施例中,该估计信号
Figure G06803960720070807D0000914
包含数据信号部分
Figure G06803960720070807D0000915
和导频信号部分
Figure G06803960720070807D0000916
在一个实施例中,如以上参考图3的估计器136所讨论的,接收机304使用下式的MMSE或ZF-MMSE离散解来推导该估计信号:
X ^ ( f ) = Y ( f ) · H ^ * ( f ) | H ^ ( f ) | 2 + σ 2 , (式7)
由此提供经均衡信号
Figure G06803960720070807D0000918
以便于进行进一步处理。
图7是示出在诸如图4所示的方法200的一部分中的一种消除残余码元间干扰(ISI)的示例性方法220的流程图。在框280,接收机304标识已均衡的一部分导频信号
Figure G06803960720070807D0000919
Figure G06803960720070807D0000920
在框282到286,接收机304使用标识出的导频信号
Figure G06803960720070807D0000921
来推导滤波器抽头(数据样本)和系数以便于执行从经均衡信号
Figure G06803960720070807D0000922
滤除残余ISI的时域均衡。特别地,在框282,接收机304至少部分基于信道估计
Figure G06803960720070807D000101
该已知导频信号xP(t)(已转换到频域成为XP(f))以及该估计的导频信号
Figure G06803960720070807D000102
来标识滤波器抽头位置(滤波器中使用的数据点)。在一个实施例中,在时域滤波器中使用的滤波器抽头可用离散形式表示为:
在频域中的 E [ k ] = 1 P · 1 H ^ [ k ] · 1 X P [ k ] · Y P [ k ] , (式9);以及
在时域中的 e [ n ] = 1 P · Σ k = - ( P / 2 ) ( P / 2 ) - 1 E [ k ] · e - j · 2 · π · k · n / P , (式10)。
在一个实施例中,对这些数据点应用一合适的阈值以使得仅较显著的数据点被选择在滤波器中使用。在一个这样的实施例中,该阈值为一被选择成使得只有那些落在主抽头(最强的,例如最大值)的约6db或10%范围内的数据点被使用的预定阈值。在另一个实施例中,该阈值被选择成使得只有那些落在主抽头的1/10与1/4之间的阈值范围内的抽头被使用。在一个实施例中,仅使用主抽头。
接下来在框284,接收机304基于该已知导频信号xp(t)与该估计的导频信号
Figure G06803960720070807D000105
之差来确定滤波器系数。假定在框282中选择了表示为τk的整数k个不同的抽头位置或数据点,则可使用以下最优化问题的MMSE或ZF-MMSE解找到使指示该导频信号xp(t)与该估计的导频信号
Figure G06803960720070807D000106
之差的值最小化的滤波器系数值a[k]的最优化值,例如如下式所示的滤波器系数aopt[k],从而来推导离散形式的滤波器系数a[k]:
a opt [ k ] = arg min a Σ n = 0 P - 1 | y ^ P ( n ) - Σ k a [ k ] · y ^ P [ n - τ k ] - x P [ n ] | 2 (式11)
移到框286,这些抽头和系数被用于从接收信号滤去残余ISI。同样可使用下式来推导移除了残余ISI的数据信号的离散形式:
x ^ D [ n ] = y ^ D [ n ] - Σ k a opt [ k ] y ^ D [ n - τ k ] . (式12)
经均衡和滤波的信号由此被提供给解调器以进行进一步的处理。
应认识到,根据本实施例,本文中所描述的这些方法的某些动作或事件可按不同顺序执行,可被添加、合并或一并省去(例如,所描述的动作或事件并非都是实施各个方法所必须的)而不会偏离本发明的范围。另外,在某些实施例中,动作或事件可被并发地执行,例如通过多线程处理、中断处理、或多处理器来并发地处理,而不是顺序处理。
本领域技术人员将可理解,信息和信号可使用各种不同技术和方法中的任意一种来表示。例如,贯穿以上说明可能被引述的数据、指令、命令、信息、信号、比特、码元以及码片可由电压、电流、电磁波、磁场或磁粒子、光场或光粒子、或其任意组合来表示。
本领域技术人员将可进一步认识到,结合本文中公开的实施例描述的各种说明性逻辑块、模块、电路以及算法步骤可被实现为电子硬件、计算机软件、或两者的组合。为清楚地说明硬件与软件的这一可互换性,以上各种说明性组件、框、模块、电路以及步骤以其功能的方式作一般化的说明。此类功能是被实现为硬件还是软件取决于具体的应用以及强加在整个系统上的的设计约束。本领域技术人员可针对每种具体应用以不同方式来实现所说明的功能,但是此类实现决策不应被解释为致使偏离本发明的范围。
结合本文中公开的实施例描述的各种说明性逻辑块、模块以及电路可用通用处理器、数字信号处理器(DSP)、专用集成电路(ASIC)、现场可编程门阵列(FPGA)或其它可编程逻辑器件,分立门或晶体管逻辑、分立硬件组件、或其设计成执行本文中所描述的这些功能的任意组合来实现或执行。通用处理器可以是微处理器,但替换地,该处理器可以是任何常规处理器、控制器、微控制器或者状态机。处理器也可被实现为计算设备的组合,例如,DSP与微处理器的组合、多个微处理器、与DSP核心协作的一个或多个微处理器、或任何其它此类配置。
结合本文中公开的实施例描述的方法或算法的步骤可直接在硬件中、在由处理器执行的软件模块中、或在这两者的组合中实施。软件模块可驻留在RAM存储器、快闪存储器、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动盘、CD-ROM、或本领域中已知的任何其它形式的存储介质中。示例性的存储介质被耦合到处理器,从而使得该处理器能从该存储介质读取信息以及向其写入信息。在替换方案中,存储介质可被集成到处理器。处理器和存储介质可驻留在ASIC中。ASIC可驻留在用户终端中。在替换方案中,处理器和存储介质可作为分立组件驻留在用户终端中。
提供所公开的实施例的以上说明是为使本领域任何普通技术人员皆可制作或使用本发明。对这些实施例的各种修改易为本领域技术人员显而易见,并且本文中所定义的一般化原理可被应用于其它实施例而不会偏离本发明的精神或范围。由此,本发明并不旨在被限定于本文中所示出的这些实施例,而是应符合与本文中公开的原理和新颖特征一致的最广义的范围。

Claims (19)

1.一种处理信号以简化计算信道行为和校正在信道上发送的信号中的计算复杂度的方法,所述方法包括:
a)接收在信道上传送的信号,所述信号包含一信息信号和已知信号的至少一部分;
b)至少部分基于所述已知信号的至少一部分来确定信道特性的至少一个指标;
c)至少部分基于信道特性的所述至少一个指标来生成指示所述信息信号的第一值,其中所述第一值包含误差信号;以及
d)至少部分基于所述已知信号的至少一部分从所述信号的第一值消除所述误差信号,
其中确定信道特性的至少一个指标包括,在频域中确定所述信道特性的估计,
其中消除所述误差信号包括,通过基于所述信道估计特性的估计、所述已知信号以及所接收的所述已知信号的至少一部分来标识滤波器抽头位置并且至少部分基于所接收的所述已知信号的至少一部分和所述已知信号之差而为所述滤波器抽头位置中每一个确定滤波器系数,来确定指示所述误差信号的值,并且
其中生成所述第一值包括,至少部分基于信道特性的所述至少一个指标在频域中均衡所接收的信号。
2.如权利要求1所述的方法,其特征在于,所述误差信号至少部分归因于码元间干扰。
3.如权利要求1所述的方法,其特征在于,所述已知信号由频域中的非零信号表征。
4.如权利要求1所述的方法,其特征在于,所述已知信号包括导频信号。
5.如权利要求1所述的方法,其特征在于,确定信道特性的所述至少一个指标进一步包括,将所述信道特性的估计从包括第一多个离散值的信号变换成包括第二多个离散值的信号,其中所述第二多个离散值的数目对应于所述信息信号的离散值的数目。
6.如权利要求1所述的方法,其特征在于,消除所述误差信号包括,对所述信息信号的所述第一值进行滤波。
7.一种配置成接收在信道上传送的信号以简化计算信道行为和校正在信道上发送的信号中的计算复杂度的设备,所述信号包含一信息信号和已知信号的至少一部分,所述设备包括:
第一电路,配置成至少部分基于所述已知信号的至少一部分来确定信道特性的至少一个指标;
第二电路,配置成至少部分基于信道特性的所述至少一个指标来生成指示所述信息信号的第一值,其中所述第一值包含误差信号;以及
第三电路,配置成至少部分基于所述已知信号的至少一部分来从所述信号的第一值消除所述误差信号,
其中所述第一电路还配置成在频域中确定所述信道特性的估计作为所述信道特性的至少一个指标,
其中所述第三电路还配置成通过基于所述信道估计特性的估计、所述已知信号以及所接收的所述已知信号的至少一部分来标识滤波器抽头位置并且至少部分基于所接收的所述已知信号的至少一部分和所述已知信号之差而为所述滤波器抽头位置中每一个确定滤波器系数,来确定指示所述误差信号的值,并且
其中所述第二电路被配置成至少部分基于信道特性的所述至少一个指标在频域中均衡所接收的信号。
8.如权利要求7所述的设备,其特征在于,所述误差信号至少部分归因于码元间干扰。
9.如权利要求7所述的设备,其特征在于,所述已知信号由频域中的非零信号表征。
10.如权利要求7所述的设备,其特征在于,所述第一电路被配置成将所述信道特性的估计从包括第一多个离散值的信号变换成包括第二多个离散值的信号,其中所述第二多个离散值的数目对应于所述信息信号的离散值的数目。
11.如权利要求7所述的设备,其特征在于,所述第三电路被配置成对所述信息信号的所述第一值进行滤波。
12.如权利要求7所述的设备,其特征在于,所述第一、第二和第三电路的至少一部分被形成在一集成电路上。
13.如权利要求7所述的设备,其特征在于,所述第一、第二和第三电路的至少一部分可被实现为一处理器。
14.一种配置成接收在信道上传送的信号以简化计算信道行为和校正在信道上发送的信号中的计算复杂度的设备,所述信号包含一信息信号和已知信号的至少一部分,所述设备包括:
至少部分基于所述已知信号的至少一部分来确定信道特性的至少一个指标的装置;
至少部分基于信道特性的所述至少一个指标来生成指示所述信息信号的第一值的装置,其中所述第一值包含误差信号;以及
至少部分基于所述已知信号的至少一部分来从所述信号的第一值消除所述误差信号的装置,
其中确定信道特性的至少一个指标的装置包括在频域中确定所述信道特性的估计的装置,
其中消除所述误差信号的装置包括通过基于所述信道估计特性的估计、所述已知信号以及所接收的所述已知信号的至少一部分来标识滤波器抽头位置并且至少部分基于所接收的所述已知信号的至少一部分和所述已知信号之差而为所述滤波器抽头位置中每一个确定滤波器系数,来确定指示所述误差信号的值的装置,并且
其中所述用于生成的装置包括至少部分基于信道特性的所述至少一个指标在频域中均衡所接收的信号的装置。
15.如权利要求14所述的设备,其特征在于,所述误差信号至少部分归因于码元间干扰。
16.如权利要求14所述的设备,其特征在于,所述已知信号包括导频信号。
17.如权利要求14所述的设备,其特征在于,所述已知信号由频域中的非零信号表征。
18.如权利要求14所述的设备,其特征在于,所述用于生成的装置包括将所述信道特性的估计从包括第一多个离散值的信号变换成包括第二多个离散值的信号的装置,其中所述第二多个离散值的数目对应于所述信息信号的离散值的数目。
19.如权利要求14所述的设备,其特征在于,所述用于消除的装置包括对所述信息信号的所述第一值进行滤波的装置。
CN2006800039607A 2005-02-04 2006-02-06 用于信道均衡的方法和系统 Active CN101116298B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US65004205P 2005-02-04 2005-02-04
US60/650,042 2005-02-04
US11/293,527 2005-12-02
US11/293,527 US8019032B2 (en) 2005-02-04 2005-12-02 Method and system for channel equalization
PCT/US2006/004248 WO2006084267A1 (en) 2005-02-04 2006-02-06 Method and system for channel equalization

Publications (2)

Publication Number Publication Date
CN101116298A CN101116298A (zh) 2008-01-30
CN101116298B true CN101116298B (zh) 2012-07-04

Family

ID=36407918

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2006800039607A Active CN101116298B (zh) 2005-02-04 2006-02-06 用于信道均衡的方法和系统

Country Status (6)

Country Link
US (2) US8019032B2 (zh)
EP (2) EP3664393B1 (zh)
JP (2) JP2008530864A (zh)
KR (1) KR101220717B1 (zh)
CN (1) CN101116298B (zh)
WO (1) WO2006084267A1 (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8019032B2 (en) * 2005-02-04 2011-09-13 Qualcomm Incorporated Method and system for channel equalization
US20080219343A1 (en) * 2007-03-09 2008-09-11 Sharp Laboratories Of America, Inc. Systems and methods for processing a signal within a communications system with a superimposed reference signal
US20100111229A1 (en) * 2008-08-08 2010-05-06 Assaf Kasher Method and apparatus of generating packet preamble
CN101478510B (zh) * 2009-02-17 2013-06-19 上海高清数字科技产业有限公司 一种自适应均衡器及使用该均衡器的接收机系统
WO2010150356A1 (ja) * 2009-06-23 2010-12-29 三菱電機株式会社 光アクセスシステム、局側終端装置および加入者側終端装置
US8385387B2 (en) * 2010-05-20 2013-02-26 Harris Corporation Time dependent equalization of frequency domain spread orthogonal frequency division multiplexing using decision feedback equalization
US9065895B2 (en) * 2012-02-22 2015-06-23 Broadcom Corporation Non-linear echo cancellation
KR102355979B1 (ko) * 2015-04-06 2022-01-26 삼성전자주식회사 통신 디바이스 및 그 제어 방법
EP3776897A4 (en) * 2018-04-13 2022-03-02 Nokia Technologies Oy TAP LOCATION INFORMATION COMPRESSION FOR EXPLICIT CHANNEL STATE INFORMATION FEEDBACK IN TIME DOMAIN IN A NEW RADIO

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6434084B1 (en) * 1999-11-22 2002-08-13 Halliburton Energy Services, Inc. Adaptive acoustic channel equalizer & tuning method
CN1435037A (zh) * 1999-12-22 2003-08-06 汤姆森特许公司 在正交频分多路复用系统中采样频率偏移的纠正

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4709374A (en) * 1984-07-05 1987-11-24 American Telephone And Telegraph Company Technique for decision-directed equalizer train/retrain
JPH08181638A (ja) 1994-12-27 1996-07-12 Matsushita Electric Ind Co Ltd 等化器及びその性能評価方法
US6314147B1 (en) * 1997-11-04 2001-11-06 The Board Of Trustees Of The Leland Stanford Junior University Two-stage CCI/ISI reduction with space-time processing in TDMA cellular networks
JP2000224148A (ja) 1999-01-29 2000-08-11 Matsushita Electric Ind Co Ltd ピンポン伝送装置
US6396886B1 (en) 1999-02-12 2002-05-28 Nec Usa, Inc. DMT time-domain equalizer algorithm
JP3084368B1 (ja) 1999-03-30 2000-09-04 株式会社次世代デジタルテレビジョン放送システム研究所 Ofdm用受信装置
EP2151962A1 (en) * 1999-04-02 2010-02-10 NTT DoCoMo, Inc. Channel estimation device and method and demodulation device
US6985432B1 (en) * 2000-01-28 2006-01-10 Zion Hadad OFDM communication channel
JP2001196978A (ja) 2000-01-11 2001-07-19 Hitachi Kokusai Electric Inc 適応等化方式及びダイバーシチ受信方式並びに適応等化器
FI108265B (fi) * 2000-02-21 2001-12-14 Tellabs Oy Menetelmä ja laitteisto adaptiivisen kanavakorjauksen opetusvaiheen toteuttamiseksi digitaalisella tietoliikenneyhteydellä
US7042937B2 (en) 2001-04-23 2006-05-09 Koninklijke Philips Electronics N.V. Hybrid frequency-time domain equalizer
US7310304B2 (en) * 2001-04-24 2007-12-18 Bae Systems Information And Electronic Systems Integration Inc. Estimating channel parameters in multi-input, multi-output (MIMO) systems
US7023928B2 (en) * 2001-08-06 2006-04-04 Lucent Technologies Inc. Synchronization of a pilot assisted channel estimation orthogonal frequency division multiplexing system
FR2829326A1 (fr) * 2001-09-06 2003-03-07 France Telecom Procede et systeme de reception iterative sous optimale pour systeme de transmission haut debit cdma
JP4298320B2 (ja) * 2002-11-08 2009-07-15 富士通株式会社 Ofdm伝送方式における受信装置
US7308063B2 (en) 2002-12-20 2007-12-11 Nokia Corporation Apparatus, and associated method, for effectuating post-FFT correction of fine frequency offset
US20040203812A1 (en) 2003-02-18 2004-10-14 Malladi Durga Prasad Communication receiver with an adaptive equalizer that uses channel estimation
JP4365125B2 (ja) 2003-03-27 2009-11-18 株式会社日立国際電気 適応等化方式
US6873596B2 (en) 2003-05-13 2005-03-29 Nokia Corporation Fourier-transform based linear equalization for CDMA downlink
KR101018381B1 (ko) * 2003-07-15 2011-03-02 삼성전자주식회사 수신 성능이 향상된 디지털 방송 송/수신 시스템 및 그의신호처리방법
US7412005B2 (en) * 2004-02-19 2008-08-12 Nokia Corporation Method and apparatus providing time domain interpolated channel estimation with advanced noise suppression for multicarrier transmissions
US8019032B2 (en) * 2005-02-04 2011-09-13 Qualcomm Incorporated Method and system for channel equalization

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6434084B1 (en) * 1999-11-22 2002-08-13 Halliburton Energy Services, Inc. Adaptive acoustic channel equalizer & tuning method
CN1435037A (zh) * 1999-12-22 2003-08-06 汤姆森特许公司 在正交频分多路复用系统中采样频率偏移的纠正

Also Published As

Publication number Publication date
KR20070110325A (ko) 2007-11-16
EP1849276B1 (en) 2020-03-18
EP3664393A1 (en) 2020-06-10
JP5174141B2 (ja) 2013-04-03
WO2006084267A1 (en) 2006-08-10
US8019032B2 (en) 2011-09-13
EP3664393B1 (en) 2023-11-08
US20060176987A1 (en) 2006-08-10
US20100310027A1 (en) 2010-12-09
JP2011125030A (ja) 2011-06-23
KR101220717B1 (ko) 2013-01-09
CN101116298A (zh) 2008-01-30
JP2008530864A (ja) 2008-08-07
EP1849276A1 (en) 2007-10-31
US8594260B2 (en) 2013-11-26

Similar Documents

Publication Publication Date Title
CN101116298B (zh) 用于信道均衡的方法和系统
JP5662528B2 (ja) ソフトパイロットシンボルを用いる無線信号を処理するための受信機および方法
US7260055B2 (en) Method for reducing channel estimation error in an OFDM system
US8483338B2 (en) Channel estimation method and system for inter-carrier interference-limited wireless communication network
EP1416688A1 (en) Iterative channel estimation in multicarrier receivers
US7065371B1 (en) Channel order selection and channel estimation in wireless communication system
US7949327B2 (en) Method and apparatus for reception of long range signals in bluetooth
US20090285086A1 (en) OFDM Receiver and OFDM Receiving Method
JP4470377B2 (ja) 移動通信システムにおける伝搬路推定方法
US7969947B2 (en) Mobile TV broadcast systems and methods based on TD-SCDMA network
JP2008206045A (ja) 無線通信システム及び無線機
US20020181430A1 (en) Space-time transmit diversity scheme for time-dispersive propagation media
JP5246771B2 (ja) 位相雑音補償受信機
US7920662B2 (en) Equalizer-based receiver to eliminate training sequence interference in TDD systems
CN102281224B (zh) Ofdm接收机及其实现方法
Faran et al. Comparison of Different Channel Estimation Techniques in OFDM Systems
CN102594742B (zh) 单载波系统中基于导频的信道岭估计方法
EP2479946A1 (en) Low complexity adaptive channel interpolation
CN104601290A (zh) 一种适合硬件实现的csi辅助判决译码方法
MX2008005275A (es) Corrección de frecuencia interactiva utilizando secuencia de entrenamiento y bits de datos

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
REG Reference to a national code

Ref country code: HK

Ref legal event code: DE

Ref document number: 1113872

Country of ref document: HK

C14 Grant of patent or utility model
GR01 Patent grant
REG Reference to a national code

Ref country code: HK

Ref legal event code: WD

Ref document number: 1113872

Country of ref document: HK