CN103326743A - Method and device for acquiring merging weight - Google Patents
Method and device for acquiring merging weight Download PDFInfo
- Publication number
- CN103326743A CN103326743A CN2012100729192A CN201210072919A CN103326743A CN 103326743 A CN103326743 A CN 103326743A CN 2012100729192 A CN2012100729192 A CN 2012100729192A CN 201210072919 A CN201210072919 A CN 201210072919A CN 103326743 A CN103326743 A CN 103326743A
- Authority
- CN
- China
- Prior art keywords
- receiver
- channel estimate
- channel
- interpolation operation
- weight
- 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 38
- 230000003044 adaptive effect Effects 0.000 claims description 32
- 238000012804 iterative process Methods 0.000 claims description 26
- 239000013598 vector Substances 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000009191 jumping Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 6
- 238000005562 fading Methods 0.000 abstract description 5
- 238000001914 filtration Methods 0.000 description 13
- 230000008859 change Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 11
- 230000004044 response Effects 0.000 description 5
- 238000013213 extrapolation Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000012549 training Methods 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005352 clarification Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011176 pooling Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
Landscapes
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
本发明提供了一种合并权重的获取方法及装置,其中,上述方法包括;选取用于计算信道的信道估值的时间段;计算所述信道在所述时间段的信道估值;对所述信道估值进行插值运算,并根据进行插值运算后得到的信道估值获得接收机的合并权重。采用本发明提供的上述技术方案,解决了相关技术中,存在的合并权重更新速度慢,无法准确反映信道变化等技术问题,提高了合并权重更新速度,以及使合并权重的粒度更小,从而达到了能够抵抗快速时变的衰落信道与一定程度的残留频偏,降低了接收机的误码率。
The present invention provides a method and device for acquiring combined weights, wherein the method includes: selecting a time period for calculating channel estimation of the channel; calculating the channel estimation of the channel in the time period; An interpolation operation is performed on the channel estimate, and a combination weight of the receiver is obtained according to the channel estimate obtained after the interpolation operation. Adopting the above-mentioned technical solution provided by the present invention solves the technical problems in related technologies such as the slow updating speed of combining weights and the inability to accurately reflect channel changes, improves the updating speed of combining weights, and makes the granularity of combining weights smaller, thereby achieving In order to resist the fast time-varying fading channel and a certain degree of residual frequency offset, the bit error rate of the receiver is reduced.
Description
技术领域 technical field
本发明涉及通信领域,具体而言,涉及一种合并权重的获取方法及装置。The present invention relates to the field of communication, in particular, to a method and device for acquiring combining weights.
背景技术 Background technique
在码分多址系统中,瑞克(Rake)接收机被广泛应用。Rake接收机对搜索到的多径进行解扰解扩,然后以一定的权重进行最大比合并。通常的合并权重即某一条多径的信道估值的共轭。随着高级接收机技术的发展,合并多径不再局限于搜索到的能量径,还有一定数量的噪声径。此时,合并权重不再是信道估值,而是基于一定算法获得的,例如最小均方误差(MMSE)或高级Rake接收机(Arake)。Arake接收机是Rake接收机的扩展,它合并的多径数多于真实多径的条数,多出的部分称为噪声径。当然,合并权重仍然是以信道估值为基础的。无论是Rake接收机还是高级接收机,合并权重的计算都至关重要,决定着最终的误码率。In CDMA systems, Rake receivers are widely used. The Rake receiver performs descrambling and despreading on the searched multipath, and then performs maximum ratio combination with a certain weight. The usual combination weight is the conjugate of the channel estimate of a certain multipath. With the development of advanced receiver technology, combined multipath is no longer limited to the searched energy path, but also has a certain number of noise paths. At this time, the combined weight is no longer channel estimation, but obtained based on a certain algorithm, such as minimum mean square error (MMSE) or advanced Rake receiver (Arake). The Arake receiver is an extension of the Rake receiver. The number of multipaths it combines is more than the number of real multipaths, and the extra part is called the noise path. Of course, the pooling weights are still based on channel estimates. Whether it is a Rake receiver or an advanced receiver, the calculation of combining weights is very important, which determines the final bit error rate.
为了对抗噪声,信道估值一般都是一段时间内平均的结果,平均后的信道估值会应用到该时间段内。当信道是快速时变的,或有较大残留频偏时,平均时间段内的信道响应也是时变的。将一个不变的信道估值应用于该时间段是不合适的,无法准确反映信道变化。信道快速时变时误码率会提高,例如车载移动台。当存在较大残留频偏时则造成星座图的旋转。In order to combat noise, the channel estimation is generally the result of averaging over a period of time, and the averaged channel estimation will be applied to this period of time. When the channel is rapidly time-varying, or has a large residual frequency offset, the channel response in the average time period is also time-varying. Applying a constant channel estimate to this time period is inappropriate and does not accurately reflect channel variation. The bit error rate will increase when the channel is rapidly time-varying, such as in-vehicle mobile stations. When there is a large residual frequency offset, the rotation of the constellation diagram is caused.
针对相关技术中的上述问题,目前尚未提出有效的解决方案。For the above-mentioned problems in related technologies, no effective solution has been proposed yet.
发明内容 Contents of the invention
针对相关技术中存在的合并权重更新速度慢,无法准确反映信道变化等技术问题,本发明提供一种合并权重的获取方法及装置,以至少解决上述技术问题。Aiming at the technical problems existing in related technologies such as slow update speed of combining weights and inability to accurately reflect channel changes, the present invention provides a method and device for acquiring combining weights to at least solve the above technical problems.
根据本发明的一个方面,提供了一种合并权重的获取方法,包括;选取用于计算信道的信道估值的时间段;计算所述信道在所述时间段的信道估值;对所述信道估值进行插值运算,并根据进行插值运算后得到的信道估值获得接收机的合并权重。According to one aspect of the present invention, there is provided a method for acquiring combined weights, including: selecting a time period for calculating the channel estimate of the channel; calculating the channel estimate of the channel in the time period; An interpolation operation is performed on the estimate, and a combination weight of the receiver is obtained according to the channel estimate obtained after the interpolation operation.
上述接收机为瑞克Rake接收机;根据进行插值运算后得到的信道估值获得接收机的合并权重,包括:对进行插值运算后得到的信道估值进行共轭运算,得到所述Rake接收机的合并权重。The above-mentioned receiver is a Rake receiver; the combination weight of the receiver is obtained according to the channel estimate obtained after the interpolation operation, including: performing a conjugate operation on the channel estimate obtained after the interpolation operation to obtain the Rake receiver The combined weight of .
上述接收机为高级瑞克Arake接收机;在对所述信道估值进行插值运算之前,上述方法还包括:获取所述Arake接收机的初始合并权重;The above-mentioned receiver is an advanced Rake Arake receiver; before performing an interpolation operation on the channel estimate, the above-mentioned method further includes: obtaining an initial combination weight of the Arake receiver;
根据进行插值运算后得到的信道估值获得接收机的合并权重,包括:以所述初始合并权重为初始值,以进行插值运算后得到的信道估值为自适应滤波器的输入,对所述初始合并权重进行调整,获取插值运算后的所述Arake接收机的合并权重。Obtaining the combination weight of the receiver according to the channel estimate obtained after the interpolation operation includes: using the initial combination weight as an initial value, and using the channel estimate obtained after the interpolation operation as the input of the adaptive filter, and performing the above-mentioned The initial combination weight is adjusted to obtain the combination weight of the Arake receiver after the interpolation operation.
通过以下迭代过程获取插值运算后的所述Arake接收机的合并权重,包括:Obtain the combined weight of the Arake receiver after the interpolation operation through the following iterative process, including:
w(0)w(0)
e(n)=d(n)-wH(n)u(n)e(n)=d(n)-w H (n)u(n)
w(n+1)=w(n)+k(n)e*(n),其中,w(0)为所述初始合并权重,w(n)为第n次迭代后的合并权重,wH(n)表示w(n)的共轭转置,u(n)为进行插值运算后得到的信道估值,d(n)为所述自适应滤波器输出的期望值,e(n)为所述自适应滤波器输出与期望值之间的差值即误差,e(n)表示e(n)的共轭,k(n)为权重更新向量。w(n+1)=w(n)+k(n)e * (n), where w(0) is the initial combined weight, w(n) is the combined weight after the nth iteration, and w H (n) represents the conjugate transpose of w(n), u(n) is the channel estimate obtained after the interpolation operation, d(n) is the expected value of the adaptive filter output, and e(n) is The difference between the adaptive filter output and the expected value is an error, e(n) represents the conjugate of e(n), and k(n) is a weight update vector.
在所述迭代过程中,所述方法还包括:在当前迭代过程中的所述误差值大于前一次迭代的误差值时,重新产生所述权重更新向量以使得迭代过程趋于收敛。During the iterative process, the method further includes: when the error value in the current iterative process is greater than the error value in the previous iterative process, regenerating the weight update vector so that the iterative process tends to converge.
所述迭代过程在满足以下之一条件时,跳出所述迭代过程:所述迭代过程的迭代次数到达预设门限值;所述误差值满足跳出迭代过程的预定条件。When the iterative process satisfies one of the following conditions, exit the iterative process: the number of iterations of the iterative process reaches a preset threshold value; the error value satisfies a predetermined condition for exiting the iterative process.
上述信道在选取的所述时间段内呈线性变化。The above channel changes linearly within the selected time period.
对所述信道估值进行插值运算,包括:采用线性样条插值算法对所述信道估值进行插值运算。Performing an interpolation operation on the channel estimate includes: performing an interpolation operation on the channel estimate by using a linear spline interpolation algorithm.
根据本发明的另一个方面,提供了一种合并权重的获取装置,位于接收机中,包括;选取模块,用于选取用于计算信道的信道估值的时间段;计算模块,用于计算所述信道在所述时间段的信道估值;插值模块,用于对所述信道估值进行插值运算;获取模块,用于根据进行插值运算后得到的信道估值获得接收机的合并权重According to another aspect of the present invention, there is provided a combination weight acquisition device, which is located in the receiver and includes: a selection module, used to select a time period for calculating the channel estimate of the channel; a calculation module, used to calculate the The channel estimate of the channel in the time period; the interpolation module is used to interpolate the channel estimate; the acquisition module is used to obtain the combination weight of the receiver according to the channel estimate obtained after the interpolation operation
上述获取模块,还用于对进行插值运算后得到的信道估值进行共轭运算,得到瑞克Rake接收机的合并权重。The acquisition module is further configured to perform a conjugate operation on the channel estimate obtained after the interpolation operation to obtain the combination weight of the Rake receiver.
上述获取模块包括:第一获取单元,用于对所述信道估值进行插值运算之前,获取高级瑞克Arake接收机的初始合并权重;第二获取单元,用于以所述初始合并权重为初始值,以进行插值运算后得到的信道估值为自适应滤波器的输入,对所述初始合并权重进行调整,获取插值运算后的所述Arake接收机的合并权重。The above acquisition module includes: a first acquisition unit, used to acquire the initial combination weight of the advanced Rake Arake receiver before performing interpolation operation on the channel estimate; a second acquisition unit, used to use the initial combination weight as the initial value, the channel estimate obtained after the interpolation operation is used as the input of the adaptive filter, and the initial combination weight is adjusted to obtain the combination weight of the Arake receiver after the interpolation operation.
通过本发明,采用基于插值后的信道估值获取合并权重的技术方案,解决了相关技术中,存在的合并权重更新速度慢,无法准确反映信道变化等技术问题,提高了合并权重更新速度,以及使合并权重的粒度更小,从而达到了能够抵抗快速时变的衰落信道与一定程度的残留频偏,降低了接收机的误码率。Through the present invention, the technical solution of obtaining combining weights based on channel estimation after interpolation is adopted, which solves the technical problems of slow updating speed of combining weights and inability to accurately reflect channel changes in related technologies, and improves the updating speed of combining weights, and The granularity of combining weights is made smaller, so as to be able to resist fast time-varying fading channels and a certain degree of residual frequency offset, and reduce the bit error rate of the receiver.
附图说明 Description of drawings
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:The accompanying drawings described here are used to provide a further understanding of the present invention and constitute a part of the application. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations to the present invention. In the attached picture:
图1为根据本发明实施例的合并权重的获取方法的流程图;FIG. 1 is a flowchart of a method for obtaining a combination weight according to an embodiment of the present invention;
图2为根据本发明优选实施例的对快速时变信道划分平均区间的示意图;FIG. 2 is a schematic diagram of dividing an average interval for a fast time-varying channel according to a preferred embodiment of the present invention;
图3为根据本发明优选实施例的获得平均区间内平均信道估值的示意图;Fig. 3 is a schematic diagram of obtaining an average channel estimate in an average interval according to a preferred embodiment of the present invention;
图4为根据本发明优选实施例的自适应滤波的原理示意图;FIG. 4 is a schematic diagram of the principle of adaptive filtering according to a preferred embodiment of the present invention;
图5为根据本发明优选实施例的采用线性样条插值算法进行插值的示意图;5 is a schematic diagram of interpolation using a linear spline interpolation algorithm according to a preferred embodiment of the present invention;
图6为基于图4所示实施例的优化后的自适应滤波的原理示意图;FIG. 6 is a schematic diagram of the principle of optimized adaptive filtering based on the embodiment shown in FIG. 4;
图7为根据本发明实施例的合并权重的获取装置的结构框图;FIG. 7 is a structural block diagram of a device for obtaining a combination weight according to an embodiment of the present invention;
图8为基于图7所示实施例的获取模块的结构框图;Fig. 8 is a structural block diagram based on the acquisition module of the embodiment shown in Fig. 7;
图9为根据本发明优选实施例的最大比合并权重的装置的结构示意图;FIG. 9 is a schematic structural diagram of a device for combining weights with a maximum ratio according to a preferred embodiment of the present invention;
图10为根据本发明优选实施例的采用LMS实现自适应滤波的原理示意图。Fig. 10 is a schematic diagram of the principle of implementing adaptive filtering by using LMS according to a preferred embodiment of the present invention.
具体实施方式 Detailed ways
下文中将参考附图并结合实施例来详细说明本发明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。Hereinafter, the present invention will be described in detail with reference to the drawings and examples. It should be noted that, in the case of no conflict, the embodiments in the present application and the features in the embodiments can be combined with each other.
图1为根据本发明实施例的合并权重的获取方法的流程图。如图1所示,该方法包括:FIG. 1 is a flow chart of a method for acquiring a combination weight according to an embodiment of the present invention. As shown in Figure 1, the method includes:
步骤S102,选取用于计算信道的信道估值的时间段;Step S102, selecting a time period for calculating the channel estimate of the channel;
步骤S104,计算所述信道在所述时间段的信道估值;Step S104, calculating the channel estimate of the channel in the time period;
步骤S106,对所述信道估值进行插值运算,并根据进行插值运算后得到的信道估值获得接收机的合并权重。Step S106, performing an interpolation operation on the channel estimate, and obtaining a combination weight of the receiver according to the channel estimate obtained after the interpolation operation.
无论信道是否为快速时变的或者有较大残留频偏,均可以采用上述处理过程获取接收机的合并权重。尤其在信道为快速时变或有较大残留频偏时,可以解决由于合并权重更新速度慢无法准确反映信道变化的技术问题,使合并权重(例如最大比合并权重)更新速度更快,粒度更小,能抵抗快速时变的衰落信道与一定程度的残留频偏,降低了接收机的误码率。Regardless of whether the channel is fast time-varying or has a large residual frequency offset, the above processing procedure can be used to obtain the combination weight of the receiver. Especially when the channel is fast time-varying or has a large residual frequency offset, it can solve the technical problem that the combination weight cannot accurately reflect the channel change due to the slow update speed of the combination weight, so that the combination weight (such as the maximum ratio combination weight) can be updated faster and more granular. Small, can resist fast time-varying fading channels and a certain degree of residual frequency offset, which reduces the bit error rate of the receiver.
在步骤S102和步骤S104中,可以采用以下处理过程实现:当信道是时变的,如图2所示,选取一定的平均区间(即时间段),使得该区间内信道变化是近似线性的。此时平均后的信道估值反映的是平均时隙中心的信道响应。基于近似线性的假设,对信道估值进行插值(例如采用线性样条插值算法来进行插值)来减小信道估值的粒度。In step S102 and step S104, the following process can be adopted: when the channel is time-varying, as shown in FIG. 2, a certain average interval (ie time period) is selected so that the channel variation in this interval is approximately linear. At this time, the averaged channel estimate reflects the channel response at the center of the average time slot. Based on the assumption of approximate linearity, interpolation is performed on the channel estimate (for example, interpolation is performed by using a linear spline interpolation algorithm) to reduce the granularity of the channel estimate.
在步骤S106中,可以采用以下处理过程对信道估值进行插值:In step S106, the channel estimate may be interpolated using the following process:
如图3所示,假设信道估值的平均窗口时间为T,那么平均的结果反映的是窗口中心位置的信道响应,令窗口中心位置对应索引n1,n2,n3,相应的信道估值为wn1,wn2,wn3。采用线性样条算法进行插值,令插值位置索引为n。As shown in Figure 3, assuming that the average window time of the channel estimation is T, then the average result reflects the channel response at the center of the window, so that the center of the window corresponds to the indices n1, n2, n3, and the corresponding channel estimation is w n1 , w n2 , w n3 . Use the linear spline algorithm for interpolation, and let the interpolation position index be n.
如果n1,n2对应的窗口分别是第一个和最后一个平均窗口,则对于n<n1或n>n3来说,插值属于外插,对于n1≤n≤n3来说,插值属于内插。外插属于预测,因此外插方法必须符合已知的趋势,否则误差会比较大。已知趋势是:频偏或是衰落会引起信道估值的旋转,但是在平均窗口内旋转不会太大,可以用一次曲线来近似。If the windows corresponding to n1 and n2 are the first and last averaging windows respectively, then for n<n1 or n>n3, the interpolation belongs to extrapolation, and for n1≤n≤n3, the interpolation belongs to interpolation. Extrapolation belongs to prediction, so the extrapolation method must conform to the known trend, otherwise the error will be relatively large. The known trend is: frequency offset or fading will cause the rotation of channel estimation, but the rotation will not be too large in the average window, and can be approximated by a linear curve.
在步骤S106中,合并权重的获取方案对Rake接收机(或称为普通Rake接收机)和高级Rake接收机而言是不同的。具体而言,In step S106, the acquisition schemes of combining weights are different for the Rake receiver (or called the common Rake receiver) and the advanced Rake receiver. in particular,
当接收机为瑞克Rake接收机时,上述合并权重可以通过以下处理过程获取:对进行插值运算后得到的信道估值进行共轭运算,得到所述Rake接收机的合并权重。并且,由于信道估值是(近似)线性变化的,因此,Rake接收机的合并权重也满足线性关系,这样,对于插值后的信道估值可以直接应用于最大比合并。When the receiver is a Rake receiver, the combination weight may be obtained through the following process: performing a conjugate operation on the channel estimate obtained after the interpolation operation to obtain the combination weight of the Rake receiver. Moreover, since the channel estimate is (approximately) linearly changing, the combining weights of the Rake receiver also satisfy a linear relationship, so that the channel estimate after interpolation can be directly applied to the maximum ratio combining.
当接收机为高级瑞克Arake接收机时,上述合并权重可以通过以下处理过程获取:在对所述信道估值进行插值运算之前,获取所述Arake接收机的初始合并权重;以所述初始合并权重为初始值,以进行插值运算后得到的信道估值为自适应滤波器的输入,对所述初始合并权重进行调整,获取插值运算后的所述Arake接收机的合并权重。When the receiver is an advanced Rake Arake receiver, the above-mentioned combined weights can be obtained through the following process: before interpolation operation is performed on the channel estimate, the initial combined weights of the Arake receiver are obtained; The weight is an initial value, and the channel estimate obtained after the interpolation operation is used as the input of the adaptive filter, and the initial combination weight is adjusted to obtain the combination weight of the Arake receiver after the interpolation operation.
其中,可以通过以下迭代过程获取插值运算后的所述Arake接收机的合并权重,包括:Wherein, the combination weight of the Arake receiver after the interpolation operation can be obtained through the following iterative process, including:
w(0)w(0)
e(n)=d(n)-wH(n)u(n)e(n)=d(n)-w H (n)u(n)
w(n+1)=w(n)+k(n)e*(n)w(n+1)=w(n)+k(n)e * (n)
其中,w(0)为所述初始合并权重,w(n)为第n次迭代后的合并权重,wH(n)表示w(n)的共轭转置,u(n)为进行插值运算后得到的信道估值,d(n)为所述自适应滤波器输出的期望值,e(n)为所述自适应滤波器输出与期望值之间的差值即误差,e(n)表示e(n)的共轭,k(n)为权重更新向量。不同的自适应算法有不同的表现形式。Among them, w(0) is the initial combined weight, w(n) is the combined weight after the nth iteration, w H (n) represents the conjugate transpose of w(n), and u(n) is the interpolation The channel estimate obtained after the operation, d(n) is the expected value of the adaptive filter output, and e(n) is the difference between the adaptive filter output and the expected value, which is the error, and e(n) represents The conjugate of e(n), k(n) is the weight update vector. Different adaptive algorithms have different manifestations.
一般而言,如图4所示,自适应滤波器的输入是接收到的训练序列,因为训练序列的期望值是已知的。经过一定长度的训练序列的滑动,横向滤波器的系数会趋于稳定。在上述迭代过程中,自适应滤波器的输入是插值后的信道估值,其实就是信道冲击响应,即训练序列相当于一个单位脉冲函数。经过信道衰落后,接收到的信号u(即进行插值运算后得到的信道估值)不再是冲击函数,但是经过横向滤波后期望能还原为单位脉冲函数。所以,如果功率归一的话wHu的输出就应该等于1,功率不归一的情况下至少输出值的虚部为零。所以,滤波器输出的期望值即wHu的实部,误差值为:In general, as shown in Figure 4, the input to the adaptive filter is the received training sequence, since the expected value of the training sequence is known. After a certain length of training sequence sliding, the coefficients of the transversal filter tend to be stable. In the above iterative process, the input of the adaptive filter is the channel estimate after interpolation, which is actually the channel impulse response, that is, the training sequence is equivalent to a unit impulse function. After channel fading, the received signal u (that is, the channel estimate obtained after interpolation) is no longer an impulse function, but it is expected to be restored to a unit impulse function after transversal filtering. Therefore, if the power is normalized, the output of w Hu u should be equal to 1, and if the power is not normalized, at least the imaginary part of the output value is zero. Therefore, the expected value of the filter output is the real part of w H u, and the error value is:
e(n)=-i*imag(wHu)e(n)=-i*imag(w H u)
其中imag表示虚部,即复数中的虚部。“i”为虚数单位,它是“-1”的一个平方根。Where imag represents the imaginary part, that is, the imaginary part in a complex number. "i" is the imaginary unit, which is a square root of "-1".
在上述迭代过程中,由于权重更新向量决定了上述迭代过程的收敛过程,为了保证收敛速度,因此,需要考虑权重更新向量的产生,权重更新向量的产生需要考虑误差值的收敛过程。初始权重更新向量的产生以最大化收敛速度为原则,一旦发现某一次迭代过程中(或当前迭代过程中)误差值大于前一次迭代的误差值,即不收敛,则需要重新产生权重更新向量以使得迭代过程趋于收敛。In the above iterative process, since the weight update vector determines the convergence process of the above iterative process, in order to ensure the convergence speed, the generation of the weight update vector needs to be considered, and the generation of the weight update vector needs to consider the convergence process of the error value. The generation of the initial weight update vector is based on the principle of maximizing the convergence speed. Once the error value in a certain iteration process (or in the current iteration process) is found to be greater than the error value of the previous iteration, that is, it does not converge, it is necessary to regenerate the weight update vector to make the iterative process tend to converge.
为了控制计算复杂度,当自适应滤波器的输出结果(插值后的信道估值)满足要求时需要结束迭代过程,上述迭代过程在满足以下之一条件时,可以跳出所述迭代过程:所述迭代过程的迭代次数到达预设门限值;所述误差值满足跳出迭代过程的预定条件。In order to control the computational complexity, the iterative process needs to be ended when the output result of the adaptive filter (the channel estimate after interpolation) meets the requirements, and the above iterative process can jump out of the iterative process when one of the following conditions is met: the The number of iterations of the iterative process reaches a preset threshold; the error value satisfies a predetermined condition for exiting the iterative process.
在步骤S102中,为了使合并权重也呈现线性变化,在选取的上述时间段内,信道特性可以呈线性变化,此处的线性变化可以是在一定允许偏移范围内的线性变化,并不限于绝对的线性变化。In step S102, in order to make the combination weight also show a linear change, the channel characteristics can change linearly within the selected time period, and the linear change here can be a linear change within a certain allowable offset range, and is not limited to Absolutely linear change.
对上述信道估值进行插值运算可以采用多种插值算法,在本发明的一个优选实施方式中,优选采用线性样条插值算法对所述信道估值进行插值运算。A variety of interpolation algorithms may be used to perform the interpolation operation on the channel estimate. In a preferred embodiment of the present invention, a linear spline interpolation algorithm is preferably used to perform the interpolation operation on the channel estimate.
如图5所示,线性样条插值与线性拟合不同的是,前者是将采样点连接成线段获得内插点,把线段向外延伸获得外插点。该插值算法为:As shown in Figure 5, the difference between linear spline interpolation and linear fitting is that the former connects the sampling points into line segments to obtain interpolation points, and extends the line segments outward to obtain extrapolation points. The interpolation algorithm is:
为了更好地理解上述实施例中的自适应滤波过程(迭代过程),以下简要说明该过程的原理。In order to better understand the adaptive filtering process (iterative process) in the above embodiment, the principle of the process is briefly described below.
如图6所示,优化后的自适应滤波原理如下,在该优化方案中,增加了自适应迭代退出机制与权重更新向量产生机制。具体如下:As shown in Figure 6, the optimized adaptive filtering principle is as follows. In this optimization scheme, an adaptive iterative exit mechanism and a weight update vector generation mechanism are added. details as follows:
步骤(1):对信道估值进行插值,每一组插值后的信道估值对应一个自适应滤波器;Step (1): Interpolating the channel estimates, each group of interpolated channel estimates corresponds to an adaptive filter;
步骤(2):利用插值后的信道估值来计算初始权重更新向量;Step (2): Use the interpolated channel estimate to calculate the initial weight update vector;
步骤(3):对某一个时刻的所有径的信道估值u进行横向滤波,滤波系数即当前合并权重w_now,即sum(u.*conj(w_now));其中sum为求和,conj为共轭运算,“.*”为点乘运算,即两个向量对应元素相乘。Step (3): Transverse filtering is performed on channel estimates u of all paths at a certain moment, and the filter coefficient is the current combined weight w_now, that is, sum(u.*conj(w_now)); where sum is the sum and conj is the total The yoke operation, ".*" is the dot multiplication operation, that is, the corresponding elements of two vectors are multiplied.
步骤(4):取横向滤波的结果的虚部为误差,满足e(n)=-i*imag(wHu);Step (4): get the imaginary part of the result of transverse filtering as error, satisfy e(n)=-i*imag(w H u);
步骤(5):如果满足最大迭代次数或跳出条件(即误差满足条件)则跳出迭代;Step (5): If the maximum number of iterations or the jumping condition is met (that is, the error satisfies the condition), then jump out of the iteration;
步骤(6):跳出迭代时如果是满足跳出条件的,则最终的权重值为最后一次迭代的权重,否则为上一次迭代的权重值;Step (6): If the exit condition is satisfied when exiting the iteration, the final weight value is the weight of the last iteration, otherwise it is the weight value of the previous iteration;
步骤(7):如果继续迭代则进而判断本次迭代是否发散,如果发散,则重新产生权重更新向量,忽略本次迭代;否则更新权重值;Step (7): If you continue to iterate, then judge whether this iteration diverges. If it diverges, regenerate the weight update vector and ignore this iteration; otherwise, update the weight value;
步骤(8):当前合并权重值的更新,满足以下公式:Step (8): The update of the current combined weight value satisfies the following formula:
w(0)w(0)
e(n)=d(n)-wH(n)u(n)e(n)=d(n)-w H (n)u(n)
w(n+1)=w(n)+k(n)e*(n);w(n+1)=w(n)+k(n)e * (n);
步骤(9):如果继续迭代则重复步骤(3)~步骤(8)。Step (9): Repeat steps (3) to (8) if the iteration continues.
在本实施例中还提供了一种合并权重的获取装置,该装置位于接收机中,用于实现上述实施例及优选实施方式,已经进行过说明的不再赘述,下面对该装置中涉及到模块进行说明。如以下所使用的,术语“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。图7为根据本发明实施例的合并权重的获取装置的结构框图。如图7所示,该装置包括:In this embodiment, there is also provided a combination weight acquisition device, which is located in the receiver, and is used to implement the above-mentioned embodiments and preferred implementation modes. What has been explained will not be described in detail. to the module for clarification. As used below, the term "module" may be a combination of software and/or hardware that realizes a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, implementations in hardware, or a combination of software and hardware are also possible and contemplated. Fig. 7 is a structural block diagram of an apparatus for acquiring combination weights according to an embodiment of the present invention. As shown in Figure 7, the device includes:
选取模块70,连接至计算模块,用于选取用于计算信道的信道估值的时间段;A
计算模块72,连接至插值模块74,用于计算所述信道在所述时间段的信道估值;
插值模块74,连接至获取模块76,用于对所述信道估值进行插值运算;An
获取模块76,用于根据进行插值运算后得到的信道估值获得接收机的合并权重。The obtaining
正如上述方法实施例中所述的,合并权重的获取方案对Rake接收机(或称为普通Rake接收机)和高级Rake接收机而言是不同的。对于Rake接收机,可以采用以下方案实现:获取模块76,还用于对进行插值运算后得到的信道估值进行共轭运算,得到Rake接收机的合并权重。As described in the foregoing method embodiments, the acquisition schemes of combining weights are different for a Rake receiver (or called a common Rake receiver) and an advanced Rake receiver. For the Rake receiver, the following solution can be adopted: the
为了使合并权重也呈现线性变化,在选取模块70选取的上述时间段内,信道特性可以呈线性变化,此处的线性变化可以是在一定允许偏移范围内的线性变化,并不限于绝对的线性变化。In order to make the combination weight also show a linear change, within the above-mentioned time period selected by the
对于高级瑞克(Arake)接收机,可以采用以下方案实现:如图8所示,上述获取模块76包括:For advanced Rick (Arake) receiver, can adopt following scheme to realize: As shown in Figure 8, above-mentioned
第一获取单元760,用于对所述信道估值进行插值运算之前,获取高级瑞克Arake接收机的初始合并权重;The first obtaining
第二获取单元762,用于以所述初始合并权重为初始值,以进行插值运算后得到的信道估值为自适应滤波器的输入,对所述初始合并权重进行调整,获取插值运算后的所述Arake接收机的合并权重The second obtaining
为了更好地理解上述实施例以下结合优选实施例和相关附图详细说明。需要说明的是,在以下实施例中,主要涉及宽带码分多址接入(Wideband Code Division Multiple Access,简称为WCDMA)基带接收机技术,与Rake接收机中的信道估计与均衡尤其相关。以下实施例详细描述了一种最大比合并权重的插值方案。In order to better understand the above embodiments, detailed descriptions will be given below in conjunction with preferred embodiments and related drawings. It should be noted that, in the following embodiments, it mainly relates to wideband code division multiple access (WCDMA for short) baseband receiver technology, and is particularly related to channel estimation and equalization in a Rake receiver. The following embodiment describes in detail an interpolation scheme of maximum ratio combining weights.
如图9所示,最大比合并权重的装置包括以下模块:As shown in Figure 9, the device for maximally combining weights includes the following modules:
控制信道解扰解扩模块90:用来对控制信道做解扰解扩处理,每一条多径一个解扰解扩资源。Control channel descrambling and despreading module 90: used to perform descrambling and despreading processing on the control channel, and each multipath has one descrambling and despreading resource.
数据信道解扰解扩模块92:用来对数据信道做解扰解扩处理,每一条多径一个解扰解扩资源,包括能量径与噪声径。Data channel descrambling and despreading module 92: used to perform descrambling and despreading processing on the data channel, each multipath has one descrambling and despreading resource, including energy path and noise path.
信道估计模块94:用于估计每一条多径的信道响应。Channel estimation module 94: for estimating the channel response of each multipath.
高级接收机合并权重计算模块96:在Arake中用于计算能量径与噪声径的合并权重,在MMSE接收机中则用于计算横向滤波器的所有抽头系数。这里泛指所有与Rake接收机不同的合并权重计算,计算出的合并权重不等于信道估值的共轭。Advanced receiver combination weight calculation module 96: used to calculate the combination weight of energy path and noise path in Arake, and used to calculate all tap coefficients of transversal filter in MMSE receiver. This generally refers to all combined weight calculations different from the Rake receiver, and the calculated combined weight is not equal to the conjugate of the channel estimate.
最大比合并模块98:对需要合并的多径的符号级数据进行最大比合并。Maximum ratio combining module 98: performing maximum ratio combining on the multipath symbol-level data that needs to be combined.
解调译码模块100:对最大比合并结果进行解调,译码,获得发送端的信息比特。Demodulation and decoding module 100: demodulate and decode the maximum ratio combination result to obtain information bits at the sending end.
信道估值的线性插值模块102:对平均后的信道估计进行插值,获得粒度更小的信道估值。Channel estimation linear interpolation module 102: Interpolate the averaged channel estimation to obtain channel estimation with smaller granularity.
自适应滤波插值模块104:以高级接收机合并权重计算的结果作为初始值,以插值后的信道估值为自适应滤波器的输入,通过自适应迭代对合并权重进行调整,获得插值后的合并权重。Adaptive filtering interpolation module 104: take the result of the combined weight calculation of the advanced receiver as the initial value, take the interpolated channel estimate as the input of the adaptive filter, adjust the combined weight through adaptive iteration, and obtain the combined Weights.
控制、数据信道的解扰解扩模块,信道估计与高级接收机合并权重计算模块、最大比合并模块、解调译码模块构成了整个信道估计与均衡的必要模块,而信道估值的插值与自适应滤波插值则是合并权重插值的核心模块,用于提高系统性能,抵抗快速时变信道与残留频偏。The descrambling and despreading module of the control and data channels, the channel estimation and advanced receiver combination weight calculation module, the maximum ratio combination module, and the demodulation and decoding module constitute the necessary modules for the entire channel estimation and equalization, and the channel estimation interpolation and Adaptive filter interpolation is the core module of combined weight interpolation, which is used to improve system performance and resist fast time-varying channels and residual frequency offset.
需要说明的是,本实施例中的信道估值的线性插值模块102相当于上述实施例中的选取模块70、计算模块72和插值模块74;本实施例中的自适应滤波插值模块104相当于上述实施例中的获取模块76。It should be noted that the channel estimation linear interpolation module 102 in this embodiment is equivalent to the
以下基于图9所示实施例中的原理来说明最大比合并权重的处理流程:The following is based on the principle in the embodiment shown in FIG. 9 to illustrate the processing flow of the maximum ratio combining weight:
第一步:对所有多径进行控制信道、数据信道的解扰解扩,获得符号级数据。Step 1: Descramble and despread the control channel and data channel for all multipaths to obtain symbol-level data.
第二步:对控制信道的符号数据进行处理,获得每条径的信道估值,并在平均窗口内进行平均。Step 2: Process the symbol data of the control channel to obtain channel estimates for each path and average them within the averaging window.
第三步:利用第二步获得的信道估值等信道统计特性计算高级接收机的合并权重。Step 3: Calculate the combining weights of the advanced receiver by using the channel statistics obtained in the second step, such as the channel estimate.
第四步:对第二步获得的信道估值进行线性样条插值,获得插值后的Rake接收机合并权重,满足以下公式:Step 4: Perform linear spline interpolation on the channel estimate obtained in the second step to obtain the combined weight of the Rake receiver after interpolation, which satisfies the following formula:
第五步:以第四步获得的插值后的信道估值为输入,以第三步获得的高级接收机的合并权重为初始值,进行自适应滤波插值,获得插值后的高级接收机合并权重。Step 5: Take the interpolated channel estimate obtained in the fourth step as input, and use the combined weight of the advanced receiver obtained in the third step as the initial value, perform adaptive filtering interpolation, and obtain the combined weight of the advanced receiver after interpolation .
第六步:对需要合并多径的数据信道的符号数据进行最大比合并,合并权重来自第四步或第五步。Step 6: Perform maximum ratio combining on the symbol data of the data channel that needs to combine multipaths, and the combining weight comes from the fourth or fifth step.
第七步:对合并后的数据进行解调与译码,获得信息比特。Step 7: Demodulate and decode the combined data to obtain information bits.
将第五步进一步拆分为如下步骤,如图10所示,自适应滤波采用LMS实现,LMS自适应滤波的迭代公式为:The fifth step is further divided into the following steps, as shown in Figure 10, the adaptive filtering is implemented by LMS, and the iterative formula of the LMS adaptive filtering is:
w(0)w(0)
e(n)=d(n)-wH(n)u(n)e(n)=d(n)-w H (n)u(n)
w(n+1)=w(n)+μu(n)e*(n)w(n+1)=w(n)+μu(n)e * (n)
μ是迭代步长,μu(n)构成了权重更新向量k(n)。μ is the iteration step size, and μu(n) constitutes the weight update vector k(n).
步骤1:将信道估值进行插值,将高级接收机合并权重复制到插值后的粒度,每一组插值后的信道估值构成一个LMS自适应滤波。Step 1: Interpolate the channel estimates, copy the advanced receiver combining weights to the interpolated granularity, and each group of interpolated channel estimates constitutes an LMS adaptive filter.
步骤2:利用插值后的信道估值u来计算初始步长,满足Step 2: Use the interpolated channel estimate u to calculate the initial step size, satisfying
其中S为输入u的功率谱密度的最大值,N为滤波器阶数,mean表示取均值,用mean{|u|2}*3来近似S。Among them, S is the maximum value of the power spectral density of the input u, N is the order of the filter, mean means to take the mean value, and use mean{|u| 2 }*3 to approximate S.
步骤3:对某一个时刻的所有径的信道估值u进行横向滤波,滤波系数即当前合并权重Step 3: Transversely filter the channel estimates u of all paths at a certain moment, and the filter coefficient is the current combined weight
步骤4:取横向滤波的结果的虚部为误差,满足e(n)=-i*imag(wHu)Step 4: Take the imaginary part of the result of transversal filtering as the error, satisfying e(n)=-i*imag(w H u)
步骤5:如果满足最大迭代次数或跳出条件则跳出LMS迭代,跳出条件为Step 5: If the maximum number of iterations or the exit condition is met, exit the LMS iteration, and the exit condition is
其中th为阈值,e为误差值。其中real为复数中的实部,imag为虚部。Where th is the threshold and e is the error value. Where real is the real part of the complex number, and imag is the imaginary part.
步骤6:跳出迭代时如果是满足跳出条件的,则最终的权重值为最后一次迭代的权重,否则为上一次迭代的权重值Step 6: If the exit condition is met when exiting the iteration, the final weight value is the weight of the last iteration, otherwise it is the weight value of the previous iteration
步骤7:如果继续迭代则进而判断本次迭代是否发散,即满足Step 7: If you continue to iterate, then judge whether this iteration diverges, that is, satisfy
|e(n)|>|e(n-1)||e(n)|>|e(n-1)|
如果发散,步长除2,重新产生权重更新向量;否则更新前一次权重值与前一次迭代误差,即令前一次权重等于当前权重,前一次迭代误差等于当前误差。If it diverges, the step size is divided by 2, and the weight update vector is regenerated; otherwise, the previous weight value and the previous iteration error are updated, that is, the previous weight is equal to the current weight, and the previous iteration error is equal to the current error.
步骤8:当前合并权重值的更新,满足以下公式:Step 8: The update of the current combined weight value satisfies the following formula:
w(0)w(0)
e(n)=d(n)-wH(n)u(n)e(n)=d(n)-w H (n)u(n)
w(n+1)=w(n)+μu(n)e*(n)。w(n+1)=w(n)+μu(n)e * (n).
步骤9:如果继续迭代则重复步骤3~步骤8。Step 9: Repeat steps 3 to 8 if iteration continues.
在另外一个实施例中,还提供了一种软件,该软件用于执行上述实施例及优选实施方式中描述的技术方案。In another embodiment, software is also provided, and the software is used to implement the technical solutions described in the above embodiments and preferred implementation manners.
在另外一个实施例中,还提供了一种存储介质,该存储介质中存储有上述软件,该存储介质包括但不限于:光盘、软盘、硬盘、可擦写存储器等。In another embodiment, there is also provided a storage medium, in which the software is stored, the storage medium includes but not limited to: optical discs, floppy disks, hard disks, rewritable memories, and the like.
显然,本领域的技术人员应该明白,上述的本发明的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明不限制于任何特定的硬件和软件结合。Obviously, those skilled in the art should understand that each module or each step of the above-mentioned present invention can be realized by a general-purpose computing device, and they can be concentrated on a single computing device, or distributed in a network formed by multiple computing devices Alternatively, they may be implemented in program code executable by a computing device so that they may be stored in a storage device to be executed by a computing device, and in some cases in an order different from that shown here The steps shown or described are carried out, or they are separately fabricated into individual integrated circuit modules, or multiple modules or steps among them are fabricated into a single integrated circuit module for implementation. As such, the present invention is not limited to any specific combination of hardware and software.
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and changes. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included within the protection scope of the present invention.
Claims (11)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210072919.2A CN103326743B (en) | 2012-03-19 | 2012-03-19 | Method and device for acquiring merging weight |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210072919.2A CN103326743B (en) | 2012-03-19 | 2012-03-19 | Method and device for acquiring merging weight |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103326743A true CN103326743A (en) | 2013-09-25 |
CN103326743B CN103326743B (en) | 2017-02-08 |
Family
ID=49195305
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210072919.2A Expired - Fee Related CN103326743B (en) | 2012-03-19 | 2012-03-19 | Method and device for acquiring merging weight |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103326743B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1524351A (en) * | 2001-05-17 | 2004-08-25 | �����ɷ� | System and method for adjusting combiner weights using an adaptive algorithm in a wireless communications system |
WO2008055946A1 (en) * | 2006-11-10 | 2008-05-15 | Telefonaktiebolaget Lm Ericsson (Publ) | Method and apparatus for efficient signal interpolation |
CN101375517A (en) * | 2006-01-26 | 2009-02-25 | 三星电子株式会社 | Apparatus and method for controlling dynamic range of weight vectors according to combining methods in a mobile station equipped with multiple antennas in high rate packet data system using code divis |
CN102006114A (en) * | 2004-12-17 | 2011-04-06 | 华为技术有限公司 | Intelligent antenna system and spread-spectrum signal receiver |
-
2012
- 2012-03-19 CN CN201210072919.2A patent/CN103326743B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1524351A (en) * | 2001-05-17 | 2004-08-25 | �����ɷ� | System and method for adjusting combiner weights using an adaptive algorithm in a wireless communications system |
CN102006114A (en) * | 2004-12-17 | 2011-04-06 | 华为技术有限公司 | Intelligent antenna system and spread-spectrum signal receiver |
CN101375517A (en) * | 2006-01-26 | 2009-02-25 | 三星电子株式会社 | Apparatus and method for controlling dynamic range of weight vectors according to combining methods in a mobile station equipped with multiple antennas in high rate packet data system using code divis |
WO2008055946A1 (en) * | 2006-11-10 | 2008-05-15 | Telefonaktiebolaget Lm Ericsson (Publ) | Method and apparatus for efficient signal interpolation |
Also Published As
Publication number | Publication date |
---|---|
CN103326743B (en) | 2017-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8654913B2 (en) | Equalisation processing | |
US8761323B2 (en) | Impairment covariance and combining weight updates during iterative turbo interference cancellation reception | |
US7801210B2 (en) | Apparatus and methods for implementing a split equalizer filter for sparse channels | |
TWI232662B (en) | Method and apparatus for varying the length of an adaptive equalizer based on Doppler frequency | |
JP3145295B2 (en) | Data receiving device | |
CN101359926B (en) | Receiver, channel estimation method and apparatus | |
JP2011507443A (en) | Time reversal pre-equalization method | |
CN102106093A (en) | Method and apparatus for communication signal processing based on mixed parametric and non-parametric estimation of impairment correlations | |
CN104519001B (en) | A kind of channel equalization method and balanced device based on RLS and LMS unified algorithms | |
WO2011092255A1 (en) | Signal processing in wireless communication receivers | |
US8565359B2 (en) | Refinement of channel estimation with a bank of filters | |
US8098717B1 (en) | Hybrid equalizer | |
US8170089B2 (en) | Apparatus for equalizing channel in frequency domain and method therefor | |
CN101222460A (en) | A Novel Frequency Domain Equalization Joint Partial Parallel Interference Cancellation Reception Method | |
CN103368605B (en) | Suppress the system and method for the interference in community | |
CN104967584B (en) | A kind of symbol level LMS adaptive equilibrium methods for short spreading code communication system | |
CN100362756C (en) | Equalization technique and associated detection technique combined receiver and receiving method thereof | |
CN103326743B (en) | Method and device for acquiring merging weight | |
CN102882816B (en) | A kind of method of multiple signals equilibrium and device | |
CN101188436A (en) | receiver | |
JP2004312145A (en) | Adaptive equalizer | |
KR101758107B1 (en) | Apparatus and method for equalizing a signal using canceling interference in a wiless communication system | |
CN101741411B (en) | Environment-adaptable equalizing receiver | |
Torabi et al. | Pre-equalization for pre-rake MISO DS-UWB systems | |
Liu et al. | A cross-layer Kalman-PDA approach to soft-decision equalization for FIR MIMO channels |
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: 20170208 |
|
CF01 | Termination of patent right due to non-payment of annual fee |