CN101228704B - 用于参量通用rake接收机中比例因数估计的方法和装置 - Google Patents
用于参量通用rake接收机中比例因数估计的方法和装置 Download PDFInfo
- Publication number
- CN101228704B CN101228704B CN2005800512418A CN200580051241A CN101228704B CN 101228704 B CN101228704 B CN 101228704B CN 2005800512418 A CN2005800512418 A CN 2005800512418A CN 200580051241 A CN200580051241 A CN 200580051241A CN 101228704 B CN101228704 B CN 101228704B
- Authority
- CN
- China
- Prior art keywords
- proportionality factor
- ratio
- rake
- transmission
- proportionality
- 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
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B2201/00—Indexing scheme relating to details of transmission systems not covered by a single group of H04B3/00 - H04B13/00
- H04B2201/69—Orthogonal indexing scheme relating to spread spectrum techniques in general
- H04B2201/707—Orthogonal indexing scheme relating to spread spectrum techniques in general relating to direct sequence modulation
- H04B2201/7097—Direct sequence modulation interference
- H04B2201/709727—GRAKE type RAKE receivers
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
一种无线通信接收机(诸如包含在无线通信网络的基站或移动站中实施的无线通信收发机中的接收机)包括一种参量G-RAKE接收机电路以及一种在每个传输间隔基础上计算参量比例因数的方法。在一个实施例中,为单独的传输时隙获得测量的损伤相关并且将其用于估计比例因数的瞬时值。然后根据一个或多个所定义的限制来约束那些瞬时值的其中一个或两个。在其他实施例中,使用多个传输时隙来增加可用于估计比例因数的测量的数目,可选地应用参数约束。其他实施例使用迭代方法和/或求解一个参数,并使用结果获得其他一个(或多个)参数。通过使用纠错/检测信息可以改进这些实施例的其中一个或多个。
Description
技术领域
本发明涉及无线通信接收机,并且具体涉及通用RAKE接收机中的比例因数(scaling parameter)估计。
背景技术
在多径传播的情况下,接收机以对应于不同传播路径的不同延迟接收所发送的信号的多个图像。“标准”RAKE接收机通过时间上将一个或多个RAKE耙指(finger)的每一个与相应的其中一个传播路径延迟对准来提高多径环境中的接收性能。每个RAKE耙指从对应于那个传播路径延迟的信号图像输出解扩通信量值(despread traffic value),并且合并电路在根据一组组合权值对耙指输出进行组合的基础上形成组合的信号。分配给每个耙指的组合权值通常被计算为那个耙指的延迟路径所计算的信道抽头的复共轭。
如果接收的信号损伤在RAKE耙指上是不相关的,则标准RAKE接收机解决方案是最佳的。然而,接收信号中的相关损伤的存在降低了标准RAKE接收机解决方案的性能,并且性能退化可能相当大。
由于各种原因,诸如对增加的数据速率使用更高的传输带宽,打算在发展的无线通信网络中使用的接收机在它们的接收信号中更可能“看到”相当显著的有色噪声电平和其他相关损伤。照此,标准RAKE接收机对于打算在这种环境中工作的接收机来说通常并不是合适的候选体系结构。
经常被称为“通用”RAKE(G-RAKE)接收机的一种更新的方法通过提供干扰抵消和信道均衡的组合来改进色散信道环境中的接收性能。为了在标准RAKE接收机体系结构上获得这些改进,G-RAKE接收机以许多重要方式修改其信号处理操作。首先,将其耙指的一个或多个置于分离路径(off-path),即偏离于为接收信号所测量的路径延迟。其次,其组合权值并不是严格的信道抽头共轭。相反,组合权值是至少部分基于损伤相关估计,其允许组合过程减小组合信号中的相关干扰。
一种类型的G-RAKE认识到损伤相关可以分解成组成元件,并且可以对这些组成元件的底层“结构”进行建模。如这里所使用的,术语“损伤”具有广泛的定义,包括但不限于以下项的一个或多个:相邻系统干扰、自身和多用户干扰、以及噪声。因为在它们的损伤相关估计中使用模型化的项,所以这些类型的G-RAKE接收机被统称为“参量”G-RAKE接收机。作为示例,为接收信号测量的损伤相关可以表示为由相应的模型拟合参数进行比例缩放的模型化的干扰相关与由相应的模型拟合参数进行比例缩放的模型化的噪声相关之和。由于建模项的结构是已知的,并且可以从一组解扩导频(pilot)值中测量短期损伤相关,例如,损伤相关估计任务被减为确定合适的模型拟合参数,所述模型拟合参数也称为比例因数。
然而,在参量G-RAKE处理中还存在问题。更具体地,某些条件使得难以适当地保持更新的比例因数。特别地,比例因数依赖于信号的每码片周期的总能量Ec以及白噪声或仅仅噪声功率(热噪声加其他干扰)N0。在某些环境下这些值迅速变化,使得参量G-RAKE接收机难以适当地保持更新的比例因数。示例包括变化的传输条件,例如与高速共享数据信道上的预定的和非预定的用户相关的大范围变化的传输功率,所述高速共享数据信道诸如在宽带CDMA(WCDMA)标准的版本5中引入的用于高速下行链路分组接入(HSDPA)的高速下行链路共享信道(HS-DSCH)。另外,在WCDMA标准的版本6中引入的增强型上行链路技术(EUL)(其也包括严格的用户调度)预示着G-RAKE接收机应用中的参数估计的复杂化。
发明内容
在如这里所教导的参量通用RAKE(G-RAKE)接收机电路的至少一个实施例中,基于每个传输间隔,估计用于对模型化的损伤相关进行比例缩放的参量比例因数,模型化的损伤相关又用于产生损伤相关估计,以用于组合权值计算。照此,涉及比例因数计算的信号条件中的间隔到间隔的潜在急剧变化并不会干扰比例因数估计过程。
因此,在一个实施例中,包括参量G-RAKE接收机电路的一个或多个处理电路被配置为基于每个传输间隔来估计一个(或多个)比例因数。这样做防止在接收的通信信号的一个传输间隔中估计的比例因数受到在接收的通信信号的另一个传输间隔中估计的比例因数的影响。
在另一个实施例中,一种在参量通用RAKE(G-RAKE)接收机中估计用于对模型化的损伤相关进行比例缩放的比例因数的方法包括,对接收的通信信号的一个或多个传输间隔的每一个,在所有或部分传输间隔上获取测量的损伤相关。该方法继续通过将测量的损伤相关表示为由比例因数的未知瞬时值进行比例缩放的模型化的损伤相关的函数并且求解该瞬时值来估计比例因数。
例如,测量的损伤相关可以表示为由第一比例因数进行比例缩放的模型化的干扰相关和由第二比例因数进行比例缩放的模型化的噪声相关的函数。可以通过执行最小二乘方(LS)拟合过程来确定这两个比例因数的瞬时值。在一个变型中,一个参数用另一个参数来表示,并且LS拟合过程简化为求解一个比例因数。在其他变型中,可以使用迭代方法来求解一个(或多个)参数。
此外,接收的通信信号中的纠错/检测信息,诸如循环冗余检验(CRC)值,可以被用来从一组候选值中选择候选比例因数。例如,候选比例因数值可以用于产生组合权值,并且可以检验从那些组合权值产生的接收符号估计中获取的解码信息中的CRC。
当然,本发明并不限于上述特征和优势。实际上,在阅读了以下详细说明以及附图之后,本领域的技术人员将会认识到其他特征和优势。
附图说明
图1是无线通信收发机的框图,它们的至少其中一个实施如这里所教导的参量G-RAKE比例因数估计。
图2是参量G-RAKE接收机电路的一个实施例的框图。
图3是图2的G-RAKE接收机电路的其他细节的框图。
图4是比例因数估计的一个实施例的逻辑流程图。
图5是用于接收的通信信号的典型时隙和间隔结构的框图。
图6是示出图4中示例的处理逻辑的更详细实施方式的逻辑流程图。
图7是比例因数估计的另一个实施例的逻辑流程图。
图8是比例因数估计的又一个实施例的逻辑流程图。
具体实施方式
图1是无线通信收发机10的框图,其经由一个或多个天线12向另一个无线通信收发机14发送通信信号,该另一个无线通信收发机14经由一个或多个天线16接收所发送的信号加上其他信号和噪声。尽管收发机10和14两者可能配置相似,但是更详细地示出了收发机14,以突出如这里所教导的参量G-RAKE处理的一个或多个实施例。
以简化的形式,收发机10包括一个或多个发送机电路18、一个或多个接收机电路20以及如所需的其他控制/处理电路22。收发机14(虽然为清楚起见进行了简化)示出了更多的细节,作为讨论的基础,并且包括一个或多个发送机电路24以及如所需的其他处理/控制电路26,并且更特别地,包括一个或多个接收机电路30,接收机电路30包括无线电处理器电路32、参量G-RAKE接收机电路34和如所需的其他接收机处理电路36,诸如被配置为从通信量符号估计中提取所发送的信息比特的解调/解码电路。
应当理解的是,例如,收发机10可以包括无线通信网络中的无线电基站,并且收发机14可以包括移动站-例如蜂窝无线电话、寻呼机、掌上型计算机、PDA或其他移动通信设备。更特别地,收发机10可以包括根据可用的WCDMA标准配置的基站,并且收发机14可以包括根据可用的WCDMA标准配置的移动站。更一般地,收发机14基本上可以包括任何类型的无线通信设备。
在WCDMA环境中,所接收的感兴趣的通信信号可以包括HS-DSCH信号(用在HSDPA中)。HS-DSCH是用于基于调度哪些HS-DSCH的传输间隔被用来向特定用户发送数据以时间多路复用的方式服务多个用户的共享分组数据信道。包括三个传输时隙的每个传输间隔因此可被用来服务不同的预定用户。当然,依赖于向给定用户发送的数据量,可以对相同的用户调度多个传输间隔。对于共享分组信道的其他示例,传输间隔可以包括不同数目的传输时隙。
在其他实施例中,收发机14可以包括无线电基站,例如WCDMA无线电基站。这里所教导的参量G-RAKE接收机电路34到基站的应用提供了许多优势,尤其在某些反向链路条件下。特别地,当影响G-RAKE比例因数的物理量在反向链路传输之间急剧变化时,接收机电路34是有利的。增强型上行链路(EUL)到WCDMA标准的扩展表示了其中影响比例因数的估计的信号条件在基站接收的反向链路信号的传输间隔之间可以较大变化的一个示例。
广泛地,那么图1应当被理解为说明性而非限制性的,并且本领域的技术人员将认识到这里所教导的用于参量G-RAKE处理的比例因数估计可以在各种接收机应用以及各种网络类型中得以实施。更特别地,这里所教导的方法提供了比例因数的每个传输间隔估计,以用于对模型化的损伤相关进行比例缩放,其是G-RAKE处理的一部分。
作为理解参量G-RAKE接收机电路34的操作的起点,应当理解为其包括被配置为在每个传输间隔的基础上估计比例因数的一个或多个处理电路,以防止在所接收的通信信号的一个传输间隔中估计的比例因数受到所接收的通信信号的另一个传输间隔中估计的比例因数所影响。
术语“处理电路”应当被广泛理解为表示功能性的结构,而不必是物理结构。例如,参量G-RAKE接收机电路34可以包括被配置为执行这里所教导的信号处理的一个或多个集成电路设备。照此,包括参量G-RAKE接收机电路34的一个(或多个)处理电路可以包括一个或多个微处理器或者DSP,以及相应的存储的计算机程序指令,该计算机程序指令配置所述一个(或多个)处理电路执行想要的信号处理。可选地,或附加地,可以使用专用集成电路(ASIC)、现场可编程门阵列(FPGA)、复杂可编程逻辑器件(CPLD)等实施参量G-RAKE接收机电路34。
转到与处理相关的细节,注意在用于解扩所接收的通信信号的通信量信道分量的RAKE耙指上得到的总损伤相关可以表示为损伤相关矩阵R,其中总损伤等于干扰加噪声。损伤相关矩阵R可以分析地表示为:
R=ECRISI+N0Rn Eq.1
其中Ec是信号的每码片周期的总能量,N0是噪声功率,RISI以及Rn分别是模型化的干扰和噪声相关矩阵。干扰矩阵RISI依赖于脉冲形状函数、码片周期、信道系数和延迟。噪声矩阵Rn依赖于脉冲形状函数和延迟。注意假设零均值,损伤相关矩阵是协方差矩阵。
因为可以从已知量或者可以(例如从接收信号的导频信道分量)了解的量获得模型化的损伤相关矩阵RISI和Rn的矩阵元素,R的估计可以认为是估计Ec和N0的问题。作为实际问题,参量G-RAKE接收机电路34可以使用导频信道对总损伤相关R进行短期估计。作为短期“快照”的这些测量的损伤相关可以表示为并且可以使用参量比例因子α和β重写方程1,以说明由于导频信道估计所需要的任何其他比例缩放。用参量比例因子,方程1可以重写为:
R=αRISI+βRn Eq.2
其中α是用于对在RISI中获取的模型化的干扰相关进行比例缩放的第一参量G-RAKE比例因数,β是用于对在Rn中获取的模型化的噪声相关进行比例缩放的第二参量G-RAKE比例因数。
在方程2中,RISI由下式给出:
并且Rn由下式给出:
其中波浪字符表示可以对应于与d1不同的接收天线16的一个,并且其中 表示复(导频)信道介质模型,RP(τ)是(滤波器)脉冲形状自相关函数,TC是接收信号的码片周期,dk是第k个G-RAKE耙指的延迟,a(d)表示d对应的是哪个接收天线16(例如接收天线1、2等)。而且,注意g值是对应于导频信道的信道系数(即由相应路径延迟的导频信道强度进行比例缩放的信道系数)。
图2和3提供了电路实施细节的非限制性示例,以支持上述处理。在图2中,可以看到参量G-RAKE接收机电路34的运转功能实施方式包括通信量相关电路40、组合器电路42、耙指放置电路44和G-RAKE处理电路46。
通信量相关电路40包括若干通信量信道“耙指”(解扩相关器)。每个通信量信道耙指从由无线电处理器电路32输出的所接收的通信信号r输出解扩通信量值y-例如,对应于复合的、天线接收的信号的时域样值的n比特数字值。因此,对有关的每个信道化代码,通信量相关电路40输出解扩通信量值的矢量y,该矢量的长度由分配给那个信道化代码的耙指的数量设定。
对应于相同信道化代码的相同码元周期的解扩通信量值由组合器电路42使用减小干扰的组合权值进行组合。更具体地,组合器电路42使用组合权值,该组合权值是至少部分基于参数地确定的接收的信号损伤相关形成的。组合权值由w=R-1h给出,其中h是对应于G-RAKE耙指延迟的净信道系数的矢量。这些组合权值用来组合解扩通信量值矢量元素,以形成通信量码元估计z,该通信量码元估计被提供到其他接收机电路36,例如用于解调/解码(z=WHy)。
G-RAKE处理电路46基于来自耙指放置电路44的耙指放置信息和信道估计信息以及基于访问接收的通信信号r的样值来确定合适的组合权值。图3示出了G-RAKE处理电路46的一个实施例的功能电路细节,其中它包括导频相关电路50、信道跟踪器52、测量的损伤相关计算器54、结构化的元素计算器56、比例因数计算器58、模型化的损伤相关计算器60、组合权值计算器62和信号质量估计器64,所述信号质量估计器64可以是信干比(SIR)估计器。
在操作中,接收的通信信号的样值被提供给导频相关电路50,其将接收的样值与导频或其他参考信号扩频序列进行相关,并且去除码元调制以生成(导频)相关值。在耙指放置电路44提供的耙指延迟处进行相关。信道跟踪器电路52接收这些导频相关并且使用它们估计或另外跟踪接收的信号的信道系数。将这些系数提供给测量的损伤相关计算器54,损伤相关计算器54也接收导频相关。计算器54被配置为在每个传输时隙的基础上获得用于接收的通信信号的损伤相关的短期测量-即“快照”损伤相关测量例如-通过从导频相关中减去信道估计以获得损伤样值并且将损伤样值彼此相关和自相关。
此外,结构化的元素计算器56接收信道估计以及来自耙指放置电路44的耙指延迟,并且使用这些来构造对应于模型化的损伤相关矩阵RISI和Rn的元素。短期损伤相关测量和模型化的损伤相关RISI和Rn被提供给比例因数计算电路58,比例因数计算电路58使用它们计算需要的参量比例因数,例如对RISI矩阵进行比例缩放的α以及对Rn进行比例缩放的β。
模型化的损伤计算器60又被配置为基于比例因数和模型化的损伤相关为接收的通信信号产生损伤相关的模型化估计,例如R=αRISI+βRn。从那个计算中,权值计算器62产生组合权值,用于组合通信量相关电路40输出的解扩通信量值的矢量。还要注意,信道估计和损伤相关的模型化估计也可以提供给SIR计算器64,除了其他应用之外,SIR计算器64生成SIR估计,用于功率控制的目的。
谨记上述框架,图4概括地示出了根据定义的时隙和/或间隔对发送的或另外处理的接收的通信信号进行如这里所教导的比例因数估计的实施例。例如,图5示出了对于有关的接收的通信信号的给定时间窗,并且示出了将信号逻辑划分为固定长度的重复的传输时隙,所述时隙被一起组成为重复的传输间隔。
在图示中,一个传输间隔包括三个传输时隙,并且这个示出的配置对应于为HSDPA传输采用每个传输时间间隔(TTI)三个时隙的WCDMA标准。更特别地,HSDPA传输包括2ms宽的三个时隙TTI。每个TTI可以用于向不同的预定用户发送数据。此外,预定的传输可以利用两个或更多个连续的TTI,依赖于将要发送给特定用户的数据量。相似的预定传输可以在WCDMA反向链路上至少与WCDMA标准的版本6的EUL规定一起使用。
在这种环境中-例如短TTI和急剧变化的传输功率-进行用来在传输间隔上计算比例因数的测量和/或计算会有问题。因此,这里所教导的比例因数计算提供了一种以确保对短TTI和(潜在的)急剧发送机功率变化有良好接收机性能的方式计算参数α和β的方法。
更一般地,这里所教导的方法提供了在每个传输间隔的基础上计算鲁棒的(可靠的)比例因数的基础,以防止在一个传输间隔计算的比例因数被另一个传输间隔上的它们的计算所影响。概括来讲,接收的信号的传输间隔/时隙结构可以由传输方法决定或者可以由参量G-RAKE接收机电路34逻辑地决定。
无论如何,图4具体示出了比例因数计算包括估计一个或多个比例因数的瞬时值(步骤100)然后约束(bounding)或另外限制该一个(或多个)瞬时值(步骤102)的实施例。在这种方法的一个具体实施例中,每个传输间隔包括一个或多个传输时隙,并且在每个传输间隔的基础上估计比例因数包括在每个时隙基础上估计比例因数。
在每个时隙基础上估计比例因数包括对传输时隙获得测量的损伤相关,基于测量的损伤相关为比例因数产生瞬时值,并且根据一个或多个定义的限制约束比例因数。更特别地,图1中的参量G-RAKE接收机电路34(例如图3中的比例因数计算器58)被配置为将测量的损伤相关表示为由第一比例因数(α)进行比例缩放的模型化的干扰相关(RISI)和由第二比例因数(β)进行比例缩放的模型化的噪声相关(Rn)的近似函数,并且求解该函数得到第一和第二比例因数。也就是说,假定接收机能够知道或确定RISI和Rn,其包括信道的了解或确定(系数g1和延迟τ1),以及G-RAKE耙指延迟d,下式中的每个量都已知,除了α和β:
其中等于测量的损伤相关矩阵的第i列,rISI,i等于模型化的干扰相关矩阵RISI的第i列,并且rn,i等于模型化的噪声相关矩阵Rn的第i列。可以经由最小二乘方(LS)拟合过程求解该方程系统。LS拟合过程可以被表示为:
其中
注意,这个最小二乘方拟合不需要应用到测量的损伤相关矩阵中的所有元素。例如,只使用对角线元素和第一非对角线元素来进行拟合操作。概括来说,是Hermitian对称的,因此只使用对角元素和上三角或下三角的其中一个。自始至终都可以使用这个性质,因此只需要计算和存储唯一的矩阵元素。
上述LS拟合过程有几个方面值得详细描述。首先,损伤协方差矩阵R是N×N维的,其中N是由参量G-RAKE接收机电路34所使用的探测耙指的数目。参量G-RAKE处理的模拟工作表明只要N≥8,α和β的瞬时估计提供了良好的解调性能。使用八个或更多探测耙指,参量G-RAKE接收机电路34可以在短期的基础上(例如每个时隙、每隔一个时隙等)采用方程7和方程8描述的拟合过程。
然而,方程7和方程8所述的过程的两个方面还需要另外解释。首先,给出方程仅仅是为了描述简单。实际上,参量G-RAKE接收机电路34只使用对应于测量的和模型化的损伤相关矩阵的上(下)三角部分的矩阵元素。一般应该使用这些矩阵的所有上(下)三角元素,以确保良好的解调性能。
值得注意的第二点是参量G-RAKE接收机电路34具有一些先验知识,它可以利用这些先验知识简化计算,即它可以利用比例因数α和β是实数值的知识,以简化LS拟合问题。为了说明这点,可以重写方程5如下:
Re{R}+jIm{R}=α(Re{RISI}+jIm{RISI})+β(Re{Rn}+jIm{Rn}) Eq.9
其中Re{}是返回复数值的实部的操作符,Im{}是返回复数值的虚部的操作符,并且j是-1的平方根。
通过合并实部和虚部项,重写方程6如下:
因为方程10的所有元素都是实数,α和β的解变为:
其中
注意RISI和Rn的对角元素产生纯实数方程。
用上述的处理步骤,如图4中总结的那样,参量G-RAKE接收机电路34允许在每个时隙的基础上进行测量和计算,该测量和计算支持比例因数的瞬时值的基于LS的估计。每个时隙计算又使得在一个传输间隔上计算的比例因数与在另一个传输间隔计算的那些比例因数之间的独立成为可能。
然而,为了确保瞬时值是鲁棒的,该方法包括示出的约束步骤。也就是说,示出的方法根据一个或多个定义的限制通过限制或另外约束瞬时值的一个或两个来约束比例因数。例如,可以根据第一比例因数α的值来约束为第二比例因数β计算的瞬时值。
对收发机14的单天线实施来说,对比例因数的一个或两个定阈值(限制)的步骤是可选的。在这种情况中,产生α和β的瞬时估计的LS拟合步骤可以提供足够的解调性能。然而,对使用多个接收天线16的实施例的模拟暗示当没有约束地使用α和/或β的瞬时值时,块出错率(BLER)曲线可能具有不希望的高错误平层(error f1oor)。错误平层的上升涉及估计的比例因数的比值上的变化(例如φ=α/β)。
好像最小二乘方估计过程会导致β变得相当小(相对于α)或者甚至为负值。这种现象在EUL和HS-DSCH情景中更为突出,因为在那些情况中干扰项是主要损伤。在任何情况中,参数比值的变化可以导致收发机14中解码器电路所使用的软信息的较大变化。这种变化因此导致潜在较差的解码器性能,导致高解码错误平层。
参数约束的一个示例是基于监视α和β的相对值,并且确定该比值是否超过某个(定义的)阈值ρ。该阈值可以根据经验确定,例如通过模拟来确定。利用这种方法进行参数约束,LS拟合过程对给定时隙或间隔产生的α的瞬时值保持不变(只要其为正),但是β的瞬时值被设为β=α(1/ρ)。这里,应当注意,在实际实施中,参量G-RAKE接收机电路34在α和β的估计中使用导频信道,并且这个实际细节意味着α不再表示码片能量Ec而是码片能量与导频能量的比值,由Ec/Ep给出。
由于用来约束或另外限制β的一个(或多个)阈值在Ec/Ep上预测,因此参量G-RAKE接收机电路34必须估计Ep。估计Ep的一个简单方法是在时间上对hHh进行平均。这个平均值可以利用具有指数滤波器的形式来计算:
Ep(n)=γEp(n-1)+(1-γ)hHh Eq.13
其中γ是0到1之间的常数。
将Ep的估计引入这个实施例得到关于β的以下关系:
在方程14中,函数max(A,B)返回A和B的最大(带符号的)值。
图6示出了上述比例因数约束的方法的示例,其中在每个时隙或每个间隔的基础上计算的α和/或β的瞬时估计被约束,这样α和β的估计足够鲁棒,用以产生组合权值和随后的码元估计和解码。
处理开始于根据上述LS拟合过程例如在每个时隙的基础上产生α和β的瞬时估计(步骤110),并且继续如上所述估计导频能量Ep(步骤112)。该处理继续评估表达式:
其中ρ是希望的(固定的或可调的)阈值常数(步骤114)。如果比例因数β大于评估量,那么将用于随后计算的β的工作值设为β的瞬时值(步骤116)。如果不是,那么将工作值β设为约束值β=αEp/ρ(步骤118)。
当然,图6的处理逻辑代表非限制性示例,并且应当理解可以按照需要或希望采用不同的方法生成比例因数的鲁棒瞬时值。例如,另一个实施例监视从最小二乘方过程获得的β值。如果β是负的或零,可将β设为标称值ε。可选地,可将β设为α的一小部分。例如,如果β<ε,则设β为ε。
在另一个实施例中,参量G-RAKE接收机电路34监视比值φ和β,以及阈值β。例如,如果比值φ超过某个阈值ρ,或者如果β<ε,则将β设为β=α/ρ。
在另一个实施例中,参量G-RAKE接收机电路34监视比值φ。如果比值φ超过阈值ρ,那么用先前时隙中计算的瞬时值替换对当前的传输时隙根据LS拟合过程计算的α和β的瞬时值。当然,这个实施例通常不用于给定传输间隔的第一时隙,以防止比例因数值在传输间隔上重复使用。然而,有此警告,很显然从相同传输间隔的先前传输时隙计算的约束之内的比例因数值可以被重复使用以代替在随后的传输时隙中计算的位于约束之外的比例因数值。
在另一个实施例中,参量G-RAKE接收机电路34为给定传输间隔(例如给定的TT I)的第一时隙计算比值φ。如果该比值小于阈值ρ,那么α和β的值被冻结(freeze)并将其用于传输间隔中的所有时隙。这个实施例因此可以节省计算周期,因为α和β的LS估计每个传输间隔只需要进行一次。然而,通常不冻结先前传输间隔的α和β的值以转移到随后的传输间隔。这种情况的一个例外是参量G-RAKE接收机电路34在连续的传输间隔上处理预定的传输(指向它的)。例如,如果收发机14是WCDMA网络中共享HS-DSCH的几个预定用户的其中之一时,会发生这种情况。而且,作为通用的建议希望不考虑阈值而希望对整个传输间隔冻结α和β和/或限制φ的变化。
在另一个实施例中,参量G-RAKE接收机电路34监视比值φ。如果该比值超过阈值ρ,那么重新应用LS拟合过程来只估计参数α。对于α的重新估计,方程2变为:
类似于之前对α和β的估计所描述的LS过程,α的LS解由下式给出:
其中
其中k=1,...,K。注意现在A是只有一列的矩阵,并且没有涉及矩阵求逆,因为要求逆的量是标量。
在另一个实施例中,其是上述方法的变型,参量G-RAKE接收机电路34利用这样的条件:根据一些定义的关系α与β成比例。在该基础上,通过将测量的损伤相关表示为只是α的函数,然后使用方程16-方程19求解α来修改LS拟合过程。一旦确定了α,就使用定义的关系计算β。
在上述实施例中,例如通过生成α和β的瞬时值并且对瞬时值的一个或两个应用约束而在每个时隙的基础上获得比例因数的鲁棒估计。这种约束例如防止比例因数β与α相比变得太小,并且因此保证对于每个传输时隙中的损伤相关估计和随后的组合权值生成有α和β的“良好”的工作值可用。
在一个或多个其他实施例中,通过在给定传输间隔内的两个或更多传输时隙上累积数据/测量来获得比例因数的鲁棒估计。也就是说,不是使用单个时隙估计α和β,而是可以将参量G-RAKE接收机电路34配置为使用来自多个时隙的信息来执行采用的LS拟合步骤以在给定的传输间隔内确定比例因数α和β的估计。
图7示出了在给定传输间隔的多于一个传输时隙上的数据/测量的累积(步骤104),这样用于比例因数估计的样值集的大小变得更大。作为示例,对于三个时隙传输间隔而言,参量G-RAKE接收机电路34可被配置为从所有三个时隙上接收的解扩导频(或通信量)值获得测量的损伤相关,以此提供了更大数量的样值用于产生由矩阵表示的快照损伤相关测量。
参量G-RAKE接收机电路34使用累积的测量估计比例因数α和β的瞬时值(步骤106),并且可选地,可以对瞬时值的一个或两个应用约束/限制(步骤108)。然而,使用更大样值集估计α和β的瞬时值(即使用多于一个时隙的测量)倾向于降低比例因数的估计“噪声”。随着估计噪声的减小,用于获得比例因数的LS拟合过程通常生成内在更鲁棒的结果。从而不需要约束。
数学上,参量G-RAKE接收机电路34求解由下式给出的多时隙估计问题:
其中
并且Aj和Pj代表给定传输间隔的第j个时隙的A矩阵和P矢量。
作为该方法的变型,参量G-RAKE接收机电路34可以被配置为对传输间隔中的每个时隙获得α和β的LS估计(例如B=Aj,x=Pj),然后一起平均这些估计。在数学形式上,这个变量由下式给出:
以及
其中k代表当前时隙索引(index),并且j从传输间隔(例如TTI)的一开始向当前时隙索引进行索引。同时注意在多个时隙上的平均仍代表每个传输间隔的解,因为时隙到时隙的平均通常不在传输间隔上执行。
在其他实施例中,参量G-RAKE接收机电路34可以有利地配置为与多个接收天线16一起使用。例如,考虑两个接收天线16,其中方程2变为:
相关矩阵Rn的非对角块的元素对应于不同天线的耙指之间的噪声相关,并且等于0。所有先前的实施例可以在多个接收天线的情况下使用。另一个实施例考虑将测量的损伤相关矩阵R的非对角元素设为零,如下所示:
然后使用前述的LS拟合过程估计比例因数α和β。
应当注意这个方法将参数α的估计值向更低的值偏离,但是它性能良好,因为它减小了比例因数的比值的变化(相比R矩阵的所有元素都用在LS估计中的情况)。如果LS拟合过程对比例因数产生负值,参量G-RAKE接收机电路34可以被配置为简单地将它们设为零。
还要注意,上述方法代表获得比例因数α和β的鲁棒估计的一步过程。然而,将理解的是约束步骤可被执行为可选的第二步骤,其中对α和/或β的值进行约束或另外加以限制以确保良好的接收机性能。
在另一个实施例中,参量G-RAKE接收机电路34使用迭代过程估计比例因数。考虑收发机14使用两个接收天线16的情况。相关矩阵Rn的非对角块的元素对应于不同天线的耙指之间的噪声相关。这些值等于零。因此,对于两个天线的情况,方程5变为方程24,在此重写如下:
从上面的表达式可以看出β的值不依赖于非对角块。因此,可以只使用对角块估计β。采用两个步骤来估计α和β。首先,只使用对角块估计这两个参数,将这些初始估计标记为α1和β1。对此,一个适当的表达式由下给出:
然后,将β设为估计值β=β1,并且使用方程26重新估计参数α,记为α2。然后通过设定α=α2并且使用方程27重新估计比例因数β以确定β2。在固定β和使用方程26重新估计α,然后固定α,使用方程27重新估计β之间的迭代,可以以预定的次数重复或者重复直到满足定义的停止标准。通常来说,在每个传输间隔基础上估计比例因数可以基于初始估计第一和第二比例因数,然后迭代地使用第一和第二比例因数的其中之一的估计修改对第一和第二比例因数的另一个的估计。
图8示出了迭代确定比例因数的另一个实施例,其是上述方法的变型。处理开始于如上所述从对角块估计α1和β1(步骤120)。处理继续计算比值φ1=α1/β1(步骤122)。然后使用β2=β1的值经由LS拟合过程估计α2(步骤124),并且计算新比值φ2=α2/β2(步骤126)。两个比值φ1和φ2与定义的阈值ρ进行比较测试(步骤128)。如果两个比值都超过ρ值,那么使用对应于两个比值中最小一个的比例因数α和β,作为α的函数根据定义的阈值(约束)限制β(步骤130)。
可选地,如果只有一个比值没有超过阈值ρ,那么使用对应于那个比值的比例因数α和β的值。因此,如果φ1超过阈值而φ2低于阈值(由步骤128和132的组合检查),那么使用对应于φ2的比例因数(步骤134)。相反,如果φ2超过阈值而φ1低于阈值(由步骤128和132的组合检查),那么使用对应于φ1的比例因数(步骤136)。
在另一个迭代求解实施例中,参量G-RAKE接收机电路34被配置为通过采用其可用的额外信息改进比例因数α和β的估计。例如,接收的通信信号通常包括纠错码和/或检错码,其可用于改进比例因数估计。
例如,假设对每个(接收信号)编码块提供了循环冗余检验(CRC)值,参量G-RAKE接收机电路34可以通过尝试比例因数α和β的不同候选对来改进它的比例因数估计。对于每个候选对,解调并解码所接收的信号,并且检查CRC。对应于没有被检查的CRC的所有参数候选对都被抛弃。对所有候选对继续该过程直到CRC检查或者处理了所有候选对。注意一个编码块可以对应于多个传输时隙,并且不同的传输时隙可以使用不同的比例因数。在这些情况中,候选对指的是一个块中的所有时隙的参数对的组。
有几种方法决定一个好的候选组。例如,假设参量G-RAKE接收机电路获得α和β的初始估计。用这些初始估计,它然后可以通过修改α和β的标称值来产生一组候选对,例如通过α=α+ε1以及β=β+ε2,其中ε1和ε2可以分别选为α和β的一小部分。另一种方法是尝试不同的α和β比值。当然,可以使用其他方法获得比例因数对的候选组,并且本领域的技术人员将意识到上述示例是非限制性的。
比例因数估计的又一个实施例使用两步骤比例因数确定过程,其对实值接收滤波器有效。(例如,对大多数WCDMA接收机实施方式来说实值接收机滤波器是合理的假设。)回忆方程2给出的损伤相关表达式,如果接收滤波器是实值的,那么Rn为实值。由于β必须为实,损伤相关矩阵R的虚部项只是α的函数。概念上,参量G-RAKE接收机电路34可以首先使用虚部项求解α,然后使用实部项加上对α计算的值求解β。按照详细过程,这个实施例由以下步骤组成:(1)给定超定方程系统 对α求解最小二乘方问题;以及(2)给定超定方程系统 使用步骤(1)的α,对β求解最小二乘方问题。
不管实施上述实施例的哪个或其组合,本发明教导了基于比例因数的每个间隔估计使用参量G-RAKE处理的无线通信设备或系统,例如收发机14。每个间隔估计改进了在涉及比例因数估计的条件在传输间隔之间潜在地急剧变化的情况中作为G-RAKE处理的一部分执行的损伤相关以及相应的组合权值计算。照此,在每个间隔基础上的比例因数估计通常防止在接收的通信信号的一个传输间隔中估计的比例因数受到在接收的通信信号的另一个传输间隔中估计的比例因数的影响。
在那点上,这里所使用的术语“在每个间隔基础上”应当给予宽泛的解释。
例如,这里所教导的比例因数估计的一些实施例确实从一个时隙到另一个时隙携带比例因数,即在一个时隙估计的一个(或多个)比例因数被重用在一个或多个后续时隙中。然而,就在给定传输间隔内计算比例因数的新值来说,通过在所有或部分所述给定传输间隔期间主要的信号状态来驱动那些计算。
谨记该宽泛的理解,本领域的技术人员将理解本发明并不限于上述特征和优势。相反,本发明仅由随附的权利要求和它们的法律等同物限定。
Claims (26)
1.一种在参量通用RAKE(G-RAKE)接收机中估计用于分别对模型化的干扰相关和模型化的噪声相关进行比例缩放的第一比例因数和第二比例因数的方法,所述方法特征在于:
在每个传输间隔基础上估计第一比例因数和第二比例因数;
至少由第一比例因数形成比值;以及
基于该比值选择性地修正第一比例因数和第二比例因数的至少其中一个。
2.如权利要求1所述的方法,其中至少由第一比例因数形成比值包括形成第一比例因数和第二比例因数的比值。
3.如权利要求2所述的方法,其中基于该比值选择性地修正第一比例因数和第二比例因数的至少其中一个包括基于将该比值与定义的阈值相比较而选择性地修正第二比例因数,使得第二比例因数与第一比例因数成比例。
4.如权利要求3所述的方法,其中基于将该比值与所定义的阈值相比较而选择性地修正第二比例因数包括如果该比值超过所定义的阈值,则设定第二比例因数与第一比例因数成比例。
5.如权利要求2所述的方法,其中每个传输间隔包括多个传输时隙,并且其中在每个传输间隔基础上估计第一比例因数和第二比例因数包括在每个传输间隔的一个或多个传输时隙中估计第一比例因数和第二比例因数。
6.如权利要求5所述的方法,其中基于该比值选择性地修正第一比例因数和第二比例因数的至少其中一个包括如果在传输时隙的第一时隙中计算的比值小于所定义的阈值,则在传输间隔的剩余传输时隙上冻结第一比例因数和第二比例因数。
7.如权利要求5所述的方法,其中基于该比值选择性地修正第一比例因数和第二比例因数的至少其中一个包括如果在传输间隔的当前传输时隙中计算的比值超过所定义的阈值,则使用在传输间隔的先前传输时隙中估计的第一比例因数和第二比例因数。
8.如权利要求2所述的方法,其中基于该比值选择性地修正第一比例因数和第二比例因数的至少其中一个包括如果该比值超过所定义的阈值,则重新计算第一比例因数的最小二乘方估计。
9.如权利要求1所述的方法,其中至少根据第一比例因数形成比值包括通过用第一比例因数和导频能量的估计的乘积除以所定义的阈值来形成比值,并且其中基于该比值选择性地修正第一比例因数和第二比例因数的至少其中一个包括基于将第二比例因数与该比值相比较而选择性地修正第二比例因数,使得如果第二比例因数大于该比值,则将第二比例因数设定为第二比例因数的瞬时值,如果第二比例因数不超过该比值,则将第二比例因数设定为等于该比值。
10.如权利要求1所述的方法,进一步包括在一个或多个随后的传输间隔中重新使用在先前传输间隔中估计的比例因数。
11.如权利要求10所述的方法,其中在一个或多个随后的传输间隔中重新使用在先前传输间隔中估计的比例因数包括在第一预定的传输间隔中估计比例因数,并且对一个或多个连续的预定传输间隔重新使用在那个第一传输间隔中估计的比例因数。
12.如权利要求1所述的方法,其中每个传输间隔包括一个或多个传输时隙,并且其中在每个传输间隔的基础上估计第一比例因数和第二比例因数包括在传输间隔的多于一个传输时隙上测量损伤相关,将测量的损伤相关表示为由第一比例因数进行比例缩放的模型化的干扰相关和由第二比例因数进行比例缩放的模型化的噪声相关的函数,以及为第一比例因数和第二比例因数求解该函数。
13.如权利要求1所述的方法,其中每个传输间隔包括一个或多个传输时隙,并且其中在每个传输间隔的基础上估计第一比例因数和第二比例因数包括在每个传输间隔内在每个传输时隙基础上估计第一比例因数和第二比例因数。
14.一种参量通用RAKE(G-RAKE)接收机电路(34),包括一个或多个处理电路(46),其被配置为在参量通用RAKE(G-RAKE)接收机电路(34)中估计用于分别对模型化的干扰相关和模型化的噪声相关进行比例缩放的第一比例因数和第二比例因数,其特征在于该一个或多个处理电路(46)被配置为:
在每个传输间隔基础上估计第一比例因数和第二比例因数;
至少由第一比例因数形成比值;以及
基于该比值选择性地修正第一比例因数和第二比例因数的至少其中一个。
15.如权利要求14所述的参量通用RAKE(G-RAKE)接收机电路(34),其中该一个或多个处理电路(46)被配置为形成所述比值作为第一比例因数与第二比例因数的比值。
16.如权利要求15所述的参量通用RAKE(G-RAKE)接收机电路(34),其中该一个或多个处理电路(46)被配置为基于将该比值与所定义的阈值相比较而选择性地修正第二比例因数,使得第二比例因数与第一比例因数成比例。
17.如权利要求16所述的参量通用RAKE(G-RAKE)接收机电路(34),其中该一个或多个处理电路(46)被配置为如果该比值超过所定义的阈值,通过设定第二比例因数与第一比例因数成比例,而选择性地修正第二比例因数。
18.如权利要求15所述的参量通用RAKE(G-RAKE)接收机电路(34),其中每个传输间隔包括多个传输时隙,并且其中该一个或多个处理电路(46)被配置为在每个传输间隔的一个或多个传输时隙中估计第一比例因数和第二比例因数。
19.如权利要求18所述的参量通用RAKE(G-RAKE)接收机电路(34),其中该一个或多个处理电路(46)被配置为如果在传输时隙的第一时隙中计算的比值小于所定义的阈值,则通过在传输间隔的剩余传输时隙上冻结第一比例因数和第二比例因数,而基于该比值选择性地修正第一比例因数和第二比例因数的至少其中一个。
20.如权利要求18所述的参量通用RAKE(G-RAKE)接收机电路(34),其中该一个或多个处理电路(46)被配置为如果在传输间隔的当前传输时隙中计算的比值超过所定义的阈值,则通过使用在传输间隔的先前传输时隙中估计的第一比例因数和第二比例因数,而基于该比值选择性地修正第一比例因数和第二比例因数的至少其中一个。
21.如权利要求15所述的参量通用RAKE(G-RAKE)接收机电路(34),其中该一个或多个处理电路(46)被配置为如果该比值超过所定义的阈值,则通过重新计算第一比例因数的最小二乘方估计,而基于该比值选择性地修正第一比例因数和第二比例因数的至少其中一个。
22.如权利要求14所述的参量通用RAKE(G-RAKE)接收机电路(34),其中该一个或多个处理电路(46)被配置为通过用第一比例因数和导频能量的估计的乘积除以所定义的阈值来形成比值,并且通过基于将第二比例因数与该比值相比较而选择性地修正第二比例因数,来基于该比值选择性地修正第一比例因数和第二比例因数的至少其中一个,使得如果第二比例因数大于该比值,则将第二比例因数设定为第二比例因数的瞬时值,如果第二比例因数不超过该比值,则将第二比例因数设定为等于该比值。
23.如权利要求14所述的参量通用RAKE(G-RAKE)接收机电路(34),其中该一个或多个处理电路(46)被配置为在一个或多个随后的传输间隔中重新使用在先前传输间隔中估计的第一比例因数和第二比例因数。
24.如权利要求14所述的参量通用RAKE(G-RAKE)接收机电路(34),其中每个传输间隔包括一个或多个传输时隙,并且其中该一个或多个处理电路(46)被配置为在每个传输间隔的基础上估计第一比例因数和第二比例因数,所述估计是通过以下进行的:在传输间隔的多于一个传输时隙上测量损伤相关,将测量的损伤相关表示为由第一比例因数进行比例缩放的模型化的干扰相关和由第二比例因数进行比例缩放的模型化的噪声相关的函数,以及为第一比例因数和第二比例因数求解该函数。
25.如权利要求14所述的参量通用RAKE(G-RAKE)接收机电路(34),其中每个传输间隔包括一个或多个传输时隙,并且其中该一个或多个处理电路(46)被配置为在每个传输间隔内在每个传输时隙基础上估计第一比例因数和第二比例因数。
26.如权利要求14所述的参量通用RAKE(G-RAKE)接收机电路(34),其中该一个或多个处理电路(46)包括:
导频相关电路(50),配置成接收所接收的通信信号的样值,将所接收的样值与导频或其他参考信号扩频序列进行相关并且去除码元调制以生成导频相关值;
信道跟踪器电路(52),配置成接收导频相关值并且使用导频相关值估计或另外跟踪所接收的通信信号的信道系数;
测量的损伤相关计算器(54),被配置来为所接收的通信信号获取短期损伤相关测量;
结构化的元素计算器(56),被配置来获取模型化的干扰相关矩阵和模型化的噪声相关矩阵的元素;
比例因数计算器(58),被配置来基于短期损伤相关测量和模型化的干扰和噪声相关矩阵计算第一比例因数和第二比例因数;以及
模型化的损伤相关计算器(60),被配置来基于第一比例因数和第二比例因数以及模型化的干扰相关矩阵和模型化的噪声相关矩阵为所接收的通信信号产生损伤相关的模型化估计;
权值计算器(62),配置成产生组合权值用于组合通信量相关电路(40)输出的解扩通信量值的矢量;
信噪比SIR计算器(64),配置成接收信道估计和损伤相关的模型化估计并生成信噪比SIR估计用于功率控制目的。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US68582505P | 2005-05-31 | 2005-05-31 | |
US60/685,825 | 2005-05-31 | ||
US11/219,626 US7848389B2 (en) | 2004-03-12 | 2005-09-02 | Method and apparatus for scaling parameter estimation in parametric generalized rake receivers |
US11/219,626 | 2005-09-02 | ||
PCT/EP2005/009512 WO2006128490A1 (en) | 2005-05-31 | 2005-09-05 | A method and apparatus for scaling parameter estimation in parametric generalized rake receivers |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101228704A CN101228704A (zh) | 2008-07-23 |
CN101228704B true CN101228704B (zh) | 2011-05-18 |
Family
ID=39859690
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2005800512418A Expired - Fee Related CN101228704B (zh) | 2005-05-31 | 2005-09-05 | 用于参量通用rake接收机中比例因数估计的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101228704B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1411635A (zh) * | 1999-10-19 | 2003-04-16 | 艾利森公司 | 用于选择瑞克接收机中相关时序的设备和方法 |
-
2005
- 2005-09-05 CN CN2005800512418A patent/CN101228704B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1411635A (zh) * | 1999-10-19 | 2003-04-16 | 艾利森公司 | 用于选择瑞克接收机中相关时序的设备和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101228704A (zh) | 2008-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101048949B (zh) | 在通用rake接收机中用于参数估计的方法和设备 | |
KR101110899B1 (ko) | 파라미터 일반화된 rake 수신기에서의 파라미터 추정을스케일링하는 방법 및 장치 | |
CN101278489A (zh) | 用于处理接收到的通信信号的方法和设备 | |
EP1206843B1 (en) | Adaptive channel estimation in a wireless communication system | |
US7590167B2 (en) | Method and apparatus for QAM demodulation in a generalized rake receiver | |
CN101617476B (zh) | 基于具偏差去除的信道化码功率估计抑制干扰的方法和设备 | |
EP1210803B1 (en) | Determination of data rate, based on power spectral density estimates | |
CN101601193B (zh) | 通信接收机中的mmse信道估计 | |
EP2162995B1 (en) | Method and apparatus for removing pilot channel amplitude dependencies from rake receiver output | |
EP2181508B1 (en) | Selective signal demodulation method and apparatus | |
CN100385816C (zh) | 使用相互关联后软门槛预处理已接收cdma信号的功率量测 | |
KR20070060116A (ko) | 확산 스펙트럼 수신기에서의 지연 선택 장치, 방법 및컴퓨터 프로그램 제품 | |
EP2514107B1 (en) | Data to pilot power ratio estimation | |
US8295417B2 (en) | Method and apparatus for efficient estimation of interference in a wireless receiver | |
CN101213762B (zh) | 无线通信接收器中的减损相关估算的方法和装置 | |
US7577190B2 (en) | Method for prediction of a channel coefficient | |
JP4448847B2 (ja) | 複雑さを低減させたスライディングウィンドウ方式による等化器 | |
CN1890891A (zh) | 一种用于无线通信体系的噪声方差估算方法及装置 | |
US7929591B2 (en) | Method and apparatus for multiple scrambling code impairment compensation | |
CN101228704B (zh) | 用于参量通用rake接收机中比例因数估计的方法和装置 | |
CN1943255B (zh) | Cdma类型的通信系统中优化调度的方法 | |
CN102460986B (zh) | 用于通信接收器的方法和装置 | |
EP1428338A2 (en) | Method and apparatus for interference signal code power and noise variance estimation | |
RU2533159C2 (ru) | Способ и система для улучшенного подавления помех посредством выбора пути | |
JP2006514452A (ja) | 上りリンクsinrの推定 |
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: 20110518 Termination date: 20190905 |
|
CF01 | Termination of patent right due to non-payment of annual fee |