CN101022433B - 高速数字接收机并行自适应盲均衡方法 - Google Patents
高速数字接收机并行自适应盲均衡方法 Download PDFInfo
- Publication number
- CN101022433B CN101022433B CN200710064140A CN200710064140A CN101022433B CN 101022433 B CN101022433 B CN 101022433B CN 200710064140 A CN200710064140 A CN 200710064140A CN 200710064140 A CN200710064140 A CN 200710064140A CN 101022433 B CN101022433 B CN 101022433B
- Authority
- CN
- China
- Prior art keywords
- data
- parallel
- vector
- equalizer
- alpha
- 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
- 238000000034 method Methods 0.000 title claims abstract description 56
- 230000003044 adaptive effect Effects 0.000 title claims description 24
- 239000013598 vector Substances 0.000 claims abstract description 78
- 238000004891 communication Methods 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 17
- 230000004044 response Effects 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 6
- BGFTWECWAICPDG-UHFFFAOYSA-N 2-[bis(4-chlorophenyl)methyl]-4-n-[3-[bis(4-chlorophenyl)methyl]-4-(dimethylamino)phenyl]-1-n,1-n-dimethylbenzene-1,4-diamine Chemical compound C1=C(C(C=2C=CC(Cl)=CC=2)C=2C=CC(Cl)=CC=2)C(N(C)C)=CC=C1NC(C=1)=CC=C(N(C)C)C=1C(C=1C=CC(Cl)=CC=1)C1=CC=C(Cl)C=C1 BGFTWECWAICPDG-UHFFFAOYSA-N 0.000 claims description 2
- 102000006942 B-Cell Maturation Antigen Human genes 0.000 claims description 2
- 108010008014 B-Cell Maturation Antigen Proteins 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 5
- 230000006870 function Effects 0.000 description 20
- 238000010586 diagram Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 9
- 238000012545 processing Methods 0.000 description 9
- 238000012549 training Methods 0.000 description 8
- 238000006243 chemical reaction Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 6
- 230000008676 import Effects 0.000 description 6
- 230000006978 adaptation Effects 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000012546 transfer Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Landscapes
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
本发明属于数字通信技术领域,其特征在于利用第k次迭代计算的误差向量ek,即利用[y1 (k)(R2-|y1 (k)|2),…,yL (k)(R2-|yL (k)|2)]T来更新均衡器抽头权系数,并根据最小均方误差准则对输入串行数据进行L路并行自适应盲均衡处理。其中L为并行路数;y1 (k),…,yL (k)分别为第k个均衡器输出向量的第1,…,L个分量;R2为一个正实数,定义为:R2=E(|an|4)/E(|an|2),其中E(·)为数学期望函数,an为通信系统发射端发送的第n个数据。该方法可应用于在高速数字通信中多径、群时延等衰落信道条件下的并行盲均衡,并可在一定程度上克服硬件实现的速度限制,适用于不需要训练序列的保密通信、单点对多点的广播通信等场合。
Description
技术领域
本发明涉及一种高速数字接收机的并行自适应盲均衡方法,属于数字通信技术领域。
背景技术
在数字通信领域中衰落信道条件下接收信号需要考虑信道畸变造成的码间干扰问题,而解决这一问题需要采用均衡,即从所收到的信号中可靠地恢复符号流的方法。在这种方法中,均衡处理一般是自适应的,即需要均衡的符号流输入到这样的系统。这些系统自动调整均衡器的传递函数使其适配于信道条件,来恢复该符号流的星座点。在很多环境下,自适应均衡处理至少以两种阶段或操作模式进行:(a)初始捕获阶段和(b)跟踪阶段。在均衡器启动时执行的初始捕获阶段,这里的算法可以是基于训练序列的(使用训练序列作为基准),也可以是盲的(不使用训练序列)。在初始阶段,均衡器试图可靠地初始收敛其输出符号流于该符号星座点的某个邻域内。在符号初始收敛之后,均衡处理进入跟踪阶段,以便将符号继续保持在符号星座点的某个邻域内。
目前,公知的数字接收机的自适应均衡方法主要有基于最小均方误差的有数据辅助自适应均衡方法和利用信号自身统计特性的无数据辅助自适应盲均衡方法。同时,自适应均衡的实现还受制于FPGA硬件速度的限制,在高速情况下需要采用并行实现方法。其中,采用数据辅助的自适应均衡方法已经有并行的实用方法,如块最小均方(BLMS,Block Least Mean Square)方法,参见Gregory A.Clark,Sanjit K.Mitra,Sydney R.Parker,“Block implementation of adaptive digital filters”,IEEE Trans.on Acoustic,Speech and SignalProcessing,vol.ASSP-29,June 1981:744-752。但采用数据辅助的自适应均衡方法具有以下弱点:使用训练序列易被截获不利于保密通信,且训练序列占用了一部分传输带宽,不适用于单点对多点通信的广播通信等情况。所以也有人提出自适应盲均衡方法,如常模方法(CMA,Constant Modulus Algorithm),参见C.Richard Johnson et al,“Blind equalizationusing the constant modulus criterion:a review”,Proc.IEEE,vol.86,Oct.1998:1927-1950。但自适应盲均衡方法由于其算法的复杂性,目前尚未见实用的并行实现方法。
发明内容
本发明针对上述问题提供一种并行的自适应盲均衡方法。它能在无训练序列辅助的情况下,利用信号自身的统计特性,自适应地调整均衡器的抽头系数,实现对输入的信号的并行盲均衡处理。通过均衡后,系统可以输出残余码间干扰非常小的信号,提高了接收信号的误码性能。而且并行的方式还可以克服因硬件速度限制而不能实现高速均衡的问题。
本发明的特征在于,它含有:
一个并行有限冲击响应(FIR,Finite-Impulse Response)滤波器,它主要是对输入数据进行并行滤波的滤波器。如设L为自然数,为并行FIR滤波器每次并行滤波处理的信号数量,则可以设该并行FIR滤波器的L个滤波器系数,即均衡器的权系数,为{w1,w2,…,wL}。在滤波过程中,它先将均衡器的输入中同相(I),正交(Q)两路信号序列分别作为实部序列与虚部序列,构成一个复信号数据序列{…,x(k-1)L+1,x(k-1)L+2,…,xkL+L,…}。再通过以下所述规则,利用串并变换及单位延迟操作将串行输入数据转换成一个2L-1路并行输入数据序列。相应的数据输入规则如下:2L-1路输入数据采用单位延迟的方式输入,即第k次并行输入L路共L个数据{xkL-L+1,xkL-L+2,…,xkL},并由单位延迟处理,等待下一次,即第k+1次新输入的L个新数据{x(k+1)L-L+1,x(k+1)L-L+2,…,x(k+1)L},以共同形成并行的2L-1路输入数据{xkL-L+1,xkL-L+2,…,xkL+L-1}。其中k为自然数,表示为第k次均衡处理;该连续的数据序列{xkL-L+1,xkL-L+2,…,xkL+L-1}分别表示为均衡器的第K*L-L+1,K*L-L+2,…,k*L+L-1个对应下标序号的输入数据。在构成一组2L-1路并行输入数据之后,可将该2L-1路数据输入到并行FIR滤波器,使其与滤波器系数进行相乘计算,得到滤波器输出结果:{y1 (k),y2 (k),…,yL (k)}。其中k为自然数,表示为第k次均衡处理;数据序列{y1 (k),y2 (k),…,yL (k)}又等同于{ykL-L+1,ykL-L+2,…,ykL},即该两个向量对应的分量是相同的。而{ykL-L+1,ykL-L+2,…,ykL}为均衡器的第k*L-L+1,k*L-L+2,…,k*L个对应下标序号的输出数据。由这种方式可产生并行处理的2*L-1个数据,以构造该并行FIR滤波器所需的输入数据块矩阵。
一个并行均衡器权向量自适应调整器,它是一个均衡权向量更新的实现方法:它通过计算由并行FIR滤波器的输出数据和一个正实数之间关系的代价函数J,得到由最小均方误差准则导出的梯度下降算法来实现权向量J的更新来实现权向量的更新。最后在此基础上对该权向量的更新进行并行化扩展,实现并行权向量的自适应调整功能。该代价函数J定义为:其中n为自然数,主要表示通信系统中信号序列的序号;yn为并行FIR滤波器处理后的第n个输出数据,也作为均衡器的第n个输出数据;E(·)为数学期望函数;|·|为求绝对值函数,对于复数则是求模函数;R2为一个正实数,定义为:R2=E(|an|4)/E(|an|2)。其中an为通信系统发射端发送的第n个数据。
该方法主要是通过处理并行FIR滤波器处理均衡器的输入信号,结合均衡器的权向量,计算得到均衡器的输出信号。其中,均衡器权向量主要通过以上权向量自适应调整器得到:通过计算代价函数J相对于权向量的梯度,并按梯度下降方式更新权向量,而使J达到最小就能得到合适的均衡器权向量。
本发明的意义在于:数字通信系统在恶劣多径信道下未经过均衡的信号出现严重的码间干扰,使接收信号的误码性能急剧恶化。而经过该均衡器均衡后的信号在很大程度上纠正了码间干扰,使接收信号的误码性能得到很大的改善。例如,在数字通信系统采用正交相移键控(QPSK)调制的情况下,由图9所示的仿真结果很好的体现了这一点。同时相对于其他传统的串行均衡方式,它所采用的并行方式可在一定程度上克服硬件实现的速度限制,适合于数字通信系统需要高速均衡的场合。
附图说明
本发明最好通过参考以下详细描述和图解本发明的附图来理解,在这些附图中;
图1图解了均衡器应用的框图;
图2图解了整个并行自适应盲均衡方法实现的框图;
图3图解了L路并行数据通过单位延时构造输入矩阵所需的2L-1个输入数据的框图;
图4图解了并行FIR滤波器的框图;
图5图解了更新均衡器权系数的框图;
图6图解了计算梯度估计向量的框图;
图8图解了由均衡器输出向量Yk的第i个元素和一个正实数R2计算得到误差向量的第i个分量的框图;
具体实施方式
高速数字通信在信道快速时变、难以准确跟踪并且高速传输会带来较大的信号群时延畸变的情况下,可以让均衡器的抽头系数的初始调整与训练序列无关,即,自恢复或“盲”恢复。参见J.G.Proakis,Digital Communications,4th ed.,McGraw-Hill Inc.N.Y.,2001。另外,它对每个数据符号进行操作,而且训练过程是不中断的,因此盲算法将更好的适应这种信道时变及群时延畸变较大的情况。盲均衡最常用的算法为CMA,是由D.N.Godard在D.N.Godard,”Self-Recovering Equalization and Carrier Tracking in Two-Dimensional Data CommunicationSystems”,IEEE Trans.on Communications,vol.COM-28,Nov.1980:1867-1875中提出的。该CMA方法主要是最小化所定义的非凸代价函数,并根据这一准则进行相应的迭代处理让均方误差最小来自适应的更新权系数。本并行盲均衡方法主要是对该方法进行并行化扩展,以实现并行盲均衡。
对于包括具有L路并行均衡滤波的自适应盲均衡器,它有L个单位延迟处理单元,用于产生每次运算所需的2*L-1个数据。在本发明中,每次新输入L个数据,并通过单位延迟处理等待下一次并行输入的L个数据,并取该L个数据的前L-1个数据。两者共同构成2*L-1个数据用于处理。该2*L-1个数据可排列成以下矩阵形式χk:
其中矩阵χk的数据分量xkL-L+1,xkL-L+2,…,xkL+L-1分别为均衡器的第k*L-L+1,k*L-L+2,…,k*L+L-1个对应下标序号的输入数据。例如第1次和第2次并行输入的数据分别为{x1,…,xL}和{xL+1,…,x2L},则构成的输入数据块矩阵如下:
由以上χk的定义,则均衡器输出向量Yk表示为:
Yk=[y1 (k),…,yL (k)]=χk·Wk *=χk·[w1,…,wL]H (2)
其中k为自然数,在这里表示第k次均衡处理;y1 (k),…,yL (k)分别为第k次均衡处理时均衡器输出向量Yk的第1,…,L个分量;同时通过均衡器的输出数据的顺序关系,可以得到y1 (k),y2 (k),…,yL (k)为ykL-L+1,ykL-L+2,…,ykL,即yi (k)为ykL-L+i,i为自然数,i=1,…,L,而ykL-L+1,ykL-L+2,…,ykL为均衡器的第k*L-L+1,k*L-L+2,…,k*L个对应下标序号的输出数据。;Wk为第k次均衡处理时的均衡器抽头权向量系数(tap weight vectorcoefficient):[w1,w2,…,wL]T,简称权向量系数;(·)*为共轭函数;(·)H为共轭转置函数。
均衡器抽头权向量系数Wk的更新按以下方式实现:
按照本发明的原理,自适应盲均衡方法借助误差信号向量更新其抽头权向量系数,该误差信号向量是:
其中ek为第k次均衡处理时的误差向量;e1 (k),…,eL (k)分别为ek的第1,…,L个分量,且为标量;y1 (k),…,yL (k)分别为第k次均衡处理后,均衡器的第1,…,L个分量输出;R2为一个正实数,也可以称为常模统计量,参见Dominique N.Godard,”Self-RecoveringEqualization and Carrier Tracking in Two-Dimensional Data Communication Systems”,IEEE Trans.on Communications,vol.COM-28,Nov.1980:1867-1875,其定义如下:
R2=E(|an|4)/E(|an|2) (4)
E(·)为数学期望函数;|·|为求绝对值函数;对于复数则是求模函数;an为通信系统发射端发送的第n个数据。根据以上的误差信号向量ek得到抽头权向量系数更新的公式为:
其中,Wk、Wk+1为第k次和第k+1次均衡处理时的均衡器权向量系数,其初始值W1设置为:[1,0,…,0]T;μBCMA为一个非常小的正实数,作为权系数更新的步长,一般取值的数量级为10-6;则由以下给出:
如将输入的2L-1个数据根据下标位置排列,则χk也可以表示为:
其中矩阵χk的数据分量xkL-L+1,xkL-L+2,…,xkL+L-1分别为均衡器的第k*L-L+1,k*L-L+2,…,k*L+L-1个对应下标序号的输入数据。为了简化公式形式,可将αi,j (k)表示为第k次均衡处理时矩阵χk的第i行,第j列的元素。这样可以根据随机梯度代替准确梯度,用以下的方法得到权向量的梯度估计:根据公式(3)中的误差εk,运用最小均方误差准则将其平方后对权向量进行求梯度计算,得到权向量的梯度估计为:
以下结合附图说明,对该L路并行自适应盲均衡方法进行说明。
图1是该方法所使用的应用场合。其中n为自然数,作为数据序列的下标;an为通信系统发射端发送的第n个数据;h为信道滤波器的单位冲击响应系数:如另设m和M分别为自然数,m为信道滤波器的单位冲击响应系数的下标,M为信道滤波器的单位冲击响应系数的长度,则信道滤波器的单位冲击响应系数h可以表示为(hm,m=1,…,M);vn为发送数据序列中第n个数据经过信道后被附加的高斯白噪声。则该通信系统发射端发送的数据序列经过信道滤波器并附加高斯白噪声后得到均衡器的第n个输入数据xn可表示为:其中i,n,M为自然数;∑i=1 M(·)表示对括号中的变量,按下标i=1,…,M进行求和的函数。yn为均衡器输入序列经过均衡器处理后得到的均衡器第n个输出数据。
图2是整个并行自适应盲均衡方法的总体实现。输入的信号数据流首先串并变换成L路并行数据。然后通过图3所示的单位延时方式,产生2L-1个公式(1)所示并行FIR滤波器输入矩阵数据所需的元素,并将其输入到并行FIR滤波器。再通过图4所示的方法,将该2L-1个输入数据与均衡器权向量进行相应的处理,实现公式(2)所示的矩阵与权向量相乘的功能,并得到均衡后的输出数据。同时通过以上实施方式中所采用的权系数更新方式(公式(5),其图解说明见图5),将前一次均衡后的输出数据、当前均衡器的输入数据、系统外部输入的正实数R2和权系数更新的步长μBCMA进行处理,实现权向量的自适应更新。另外图2中的关键部件也可由以下图所描述:图2中的“单位延迟处理:L路数据变换为2L-1路数据”框图可以由图3描述;图2中的“L路并行FIR滤波器”框图可以由图4描述;图2中的“L路并行均衡器权向量自适应调整器”框图可以由图5描述。其他,图2中的“1:L串并变换”框图和“L:1并串变换”框图分别表示串行1路数据转换为L路并行数据的串并变换及L并行路数据转换为串行1路数据的并串变换。
图3是在该方法中,将连续2次的L路并行输入数据组成公式(1)所示的并行FIR滤波器输入矩阵2*L-1个数据元素,其中L为自然数,是并行路数;k为自然数,表示为第k次均衡处理。图中均衡器的输入数据序列{xkL-L+1,xkL-L+2,…,xkL+L}分别为均衡器的第k*L-L+1,k*L-L+2,…,k*L+L个连续输入数据。该模块每次并行输入L路共L个数据{x(k-1)L+1,,xkL},并结合单位延迟单元D延迟等待下一次新输入的L个新数据{xkL+1,,xkL+L},由这种方式可构造公式(1)所需要的输入数据块矩阵。在图3中“1:L路串并变换”框图描述了1路数据变换到L路数据的串并转换过程;D为单位延迟处理单元:对于经过D的串行数据序列,则在下一次均衡处理时,该数据仍保持不变输出,相当于该数据序列延迟了一个处理单位时间。
图4是该方法在第k次均衡处理时,将2*L-1个输入数据xkL-L+1,xkL-L+2,…,xkL+L-1与均衡器权向量Wk的各个分量w1,w2,…,wL以图4的方式进行相应的处理,得到输出向量Yk的各个分量y1 (k),…,yL (k)。即实现公式(2)中的输入数据矩阵χk与均衡器权向量Wk相乘,得到输出向量Yk的功能。
图5是该方法实现公式(5)中的权向量系数更新。该图描述了在第k次均衡处理时,将其中已计算得到的权向量Wk减去权系数更新步长μBCMA与梯度估计向量的乘积并得到新的权向量Wk+1,以实现公式(2)中权系数的更新功能。而且本次权向量计算输出结果Wk+1是作为下一次均衡处理时的输入权向量,在图中体现为本次计算所得的权向量通过并行延迟处理单元PD后,可以作为下一次均衡处理时所用到的权向量。其中,Wk、Wk+1为第k次和第k+1次均衡处理时的均衡器权向量系数;PD为并行L路单位延迟处理单元:对于经过PD的L路并行数据序列,则在下一次均衡处理时,该L路数据仍保持不变输出,相当于该L路数据序列延迟了一个处理单位时间。另外因为第1次权向量W1之前没有其他权向量的计算,则需要设置权向量的初值,即设置W1=[1,0,…,0]T。其他,该图所描述的系统在第k次均衡处理时向图4所示的并行FIR滤波器输出权向量系数Wk,并在该并行FIR滤波器计算得y1 (k),…,yL (k),作为权向量系数更新所需的输入数据输入该图所示的系统;第k次均衡处理时权向量更新所需的输入数据序列xkL-L+1,xkL-L+2,…,xkL+L-1,权系数更新步长μBCMA和一个正实数R2也由外部输入;梯度估计向量的计算则由图5所描述的梯度估计向量的计算模块实现,而该模块的具体实现可由图6来具体描述。
图6是该方法计算第k次均衡处理时的梯度估计向量该图主要描述了在第k次均衡处理输入数据序列xkL-L+1,xkL-L+2,…,xkL+L-1和y1 (k),…,yL (k)时,利用梯度估计向量的分量计算模块,实现公式(7)的功能。其中梯度估计向量的分量计算模块可由图7来具体描述。
图7是该方法实现公式(7)中计算第k次均衡处理时梯度估计向量的第j个分量δj (k)的功能。其中,j为自然数,取值为1,…,L;conj(·)是求共轭函数;|·|为求绝对值函数,对于复数是求模函数。根据公式(7)可知,且由公式(6)可得α1,j (k),α2,j (k),…,αL,j (k)即为xkL-L+1,xkL-L+i+1,…,xkL+L-1。则δj (k)可以由图7中所描述的方法实现。其中,(ei (k))*为ei (k)的共轭,而ei (k)的计算可由图8来具体描述。
图8是由均衡器输出向量Yk的一个分量yi (k)和正实数R2根据公式(3)计算得到误差向量的分量ei (k)。其中图中i是自然数,其取值范围为1,…,L;正实数R2在公式(4)中定义,但为了简化实现,一般都是根据调制方式已经在该方法实现前通过公式(4)计算完成的,在本方法中可以根据该计算结果把R2当作一个固定的正常实数。通过L个ei (k)(i=1,…,L)的计算,可以得到公式(3)中第k次均衡处理时的误差向量ek。
图9是为了评价所建议的并行自适应盲均衡方法,以MATLAB编程语言进行仿真。通用的仿真系统图解如图1所示,而图9显示了该方法的仿真性能结果。其中曲线表示在QPSK调制方式下信号在经过恶劣多径信道及高斯白噪声(AWGN)影响下经过该均衡方法处理后的误码曲线;曲线“——”表示理论上在QPSK调制方式下信号经过只存在AWGN影响时的误码曲线;曲线表示在QPSK调制方式下信号在经过恶劣多径信道及高斯白噪声(AWGN)影响下未经过该均衡方法处理的误码曲线。在该通信系统中,发送端所发送的为QPSK信号,采用随机数生成的方式产生QPSK符号,恶劣的多径信道滤波器为有限脉冲响应(FIR,Finite-Impulse Response)型滤波器,其FIR冲击响应为{0.227,0.460,0.688,0.460,0.227},参见J.G. Proakis,Digital Communications,4th ed.,McGraw-Hill Inc.N.Y.,2001(P.630,图.10.2-5c)。仿真采用独立50次的蒙特卡罗仿真,采用4倍分数间隔上采样,在均衡后再1/4抽取的方式实现数据均衡。仿真条件如下:16000个数据样本点;32路并行;{0.227,0.460,0.688,0.460,0.227}的多径信道滤波器;附加AWGN高斯白噪声;权系数初值为32维向量[1,0,…,0]T;μBCMA取为2*10-7。根据以上条件得到图8所示的该方法的QPSK误码性能结果。仿真结果显示,在低信噪比(Eb/N0<7dB)情况下,该方法得到的误码曲线相距理论值不超过3dB。同时它相对于未经过该均衡方法处理后的误码曲线可获得更好的误码性能,以实现较好的抗ISI功能。
Claims (1)
1.高速数字接收机并行自适应盲均衡方法,其特征在于,该方法是在数字集成电路上依次按以下步骤实现的:
步骤(1).从均衡器第k次输入以供均衡处理的长度为2*L的数据序列{xkL-L+1,xkL-L+2,…,xkL,xkL+1,xkL+2,…,xkL+L-1,xkL+L}中提取前2*L-1个数据,去除第k*L+L个数据xkL+L的数据序列,其步骤如下:
步骤(1.1).把该数据序列{xkL-L+1,xkL-L+2,…,xkL,xkL+1,xkL+2,…,xkL+L-1,xkL+L}输入到一个1∶L串/并变换电路,其中,L是一个自然数,为并行的路数;
步骤(1.2).该1∶L串/并变换电路对其中前L个数据{xkL-L+1,xkL-L+2,…,xkL}分别用单位延迟开关进行延迟,并和后面输出的L-1个数据{xkL+1,xkL+2,…,xkL+L-1}共同形成路数为2*L-1的并行数据序列;
步骤(1.3).把步骤(1.2)中的2L-1个并行的数据排列成以下矩阵形式:
其中矩阵χk的数据分量xkL-L+1,xkL-L+2,…,xkL+L-1分别为均衡器的第k*L-L+1,k*L-L+2,…,k*L+L-1个对应下标序号的输入数据;
步骤(2).把步骤(1.3)中得到的输入矩阵2*L-1个数据元素输入到一个L路并行有限冲击响应滤波器,使该2*L-1个数据与L个滤波器系数即均衡器的权系数向量Wk的各个分量,即w1,w2,…,wL相乘,得到由所述1∶L串/并变换电路和L路并行有限冲击响应滤波器构成的均衡器输出向量Yk=[y1 (k),…,yL (k)]=χk·Wk *=χk·[w1,…,wL]H,其中Wk *为权系数向量Wk的共轭函数,(·)H为共轭转置函数;
步骤(3).把步骤(2)得到的输出向量Yk送入一个L∶1并/串变换电路,得到均衡器的串行输出数据序列;
步骤(4).更新步骤(2)所述的权系数向量Wk=[w1,w2,…,wL]T,以便进行第K+1次均衡处理,其步骤如下:
步骤(4.1).把步骤(1.3)所述的矩阵χk改写成以下形式:
即把αi,j (k)表示为第k次均衡处理时矩阵χk的第i行,第j列的元素;
其中,R2为常模统计量,是一个正实数,其按下式计算:
R2=E(|an|4)/E(|an|2)
an为通信系统发射端发送的第n个数据,
E(·)为数学期望函数,
|·|为求绝对值函数,对复数是求模函数,
ek为第k次均衡时的误差向量,e1 (k),…,eL (k)分别为ek的第1,...,L个分量,ek由下式表示:
步骤(4.4).按下式计算均衡器第k+1次作均衡处理时的权向量更新系数Wk+1:
在k=1时,W1的取值为L维向量[1,0,…,0]T;
步骤(4.5).在设定的时间延迟后,按步骤(1)~(4.4)做第k+1次的均衡处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710064140A CN101022433B (zh) | 2007-03-02 | 2007-03-02 | 高速数字接收机并行自适应盲均衡方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710064140A CN101022433B (zh) | 2007-03-02 | 2007-03-02 | 高速数字接收机并行自适应盲均衡方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101022433A CN101022433A (zh) | 2007-08-22 |
CN101022433B true CN101022433B (zh) | 2010-05-19 |
Family
ID=38710088
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200710064140A Expired - Fee Related CN101022433B (zh) | 2007-03-02 | 2007-03-02 | 高速数字接收机并行自适应盲均衡方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101022433B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101567863B (zh) * | 2008-04-24 | 2011-07-20 | 魏昕 | 浅海水声通信系统的间接自适应均衡方法 |
CN101478509B (zh) * | 2009-01-20 | 2011-05-18 | 南京信息工程大学 | 正交小波变换和时间分集技术相融合的盲均衡方法 |
CN101599930B (zh) * | 2009-07-23 | 2011-09-21 | 西安空间无线电技术研究所 | 高速并行均衡器及均衡方法 |
CN101616376B (zh) * | 2009-07-30 | 2011-12-28 | 中兴通讯股份有限公司 | 一种码分多址通信系统中的数据接收方法及装置 |
CN101808412B (zh) * | 2009-10-27 | 2012-09-05 | 清华大学 | 多基站多用户网络的调度方法 |
CN102209055B (zh) | 2010-03-30 | 2014-10-01 | 富士通株式会社 | 用于相移键控信号的自适应盲均衡方法、均衡器及接收机 |
CN102082749B (zh) * | 2010-12-24 | 2014-01-15 | 西安空间无线电技术研究所 | 一种均衡器的权系数更新装置及方法 |
DE102014201233B4 (de) * | 2014-01-23 | 2023-12-14 | Rohde & Schwarz GmbH & Co. Kommanditgesellschaft | Entzerrungsfilter |
CN104144139B (zh) * | 2014-08-18 | 2017-09-01 | 西安烽火电子科技有限责任公司 | 一种改进的时域自适应盲均衡方法 |
CN107786475B (zh) * | 2016-08-26 | 2020-04-10 | 深圳市中兴微电子技术有限公司 | 盲均衡误差计算方法和装置 |
CN107809400B (zh) * | 2017-09-19 | 2020-07-28 | 北京理工大学 | 一种基于fpga的高速低复杂度并行盲均衡实现方法 |
CN108616466B (zh) * | 2018-05-09 | 2020-12-01 | 清华大学 | 信道并行均衡方法及装置 |
CN110224953B (zh) * | 2019-05-31 | 2021-10-22 | 西安理工大学 | 用于高速串行接口的均衡器及其实现盲均衡自适应的方法 |
CN110247642B (zh) * | 2019-06-13 | 2023-07-11 | 江苏卓胜微电子股份有限公司 | 一种fir滤波方法及滤波器 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6049574A (en) * | 1998-04-17 | 2000-04-11 | Trustees Of Tufts College | Blind adaptive equalization using cost function that measures dissimilarity between the probability distributions of source and equalized signals |
US6490318B1 (en) * | 1999-06-24 | 2002-12-03 | Agere Systems Inc. | Phase-compensating constant modulus algorithm |
CN1574812A (zh) * | 2003-06-09 | 2005-02-02 | 三星电子株式会社 | 多载波接收装置的盲均衡器及其方法 |
-
2007
- 2007-03-02 CN CN200710064140A patent/CN101022433B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6049574A (en) * | 1998-04-17 | 2000-04-11 | Trustees Of Tufts College | Blind adaptive equalization using cost function that measures dissimilarity between the probability distributions of source and equalized signals |
US6490318B1 (en) * | 1999-06-24 | 2002-12-03 | Agere Systems Inc. | Phase-compensating constant modulus algorithm |
CN1574812A (zh) * | 2003-06-09 | 2005-02-02 | 三星电子株式会社 | 多载波接收装置的盲均衡器及其方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101022433A (zh) | 2007-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101022433B (zh) | 高速数字接收机并行自适应盲均衡方法 | |
US7027504B2 (en) | Fast computation of decision feedback equalizer coefficients | |
CN111082870B (zh) | 相干光通信中时钟恢复和自适应均衡器联合装置及方法 | |
NO317331B1 (no) | Fremgangsmate og anordning for digitalt radiosignalmottak | |
JP2770839B2 (ja) | 適応受信機 | |
US7830956B2 (en) | Method and system for processing a sampled signal | |
CN113904898B (zh) | 一种基于等效时间采样的均衡与载波同步方法 | |
CN112468419B (zh) | 一种自适应的双模式盲均衡方法及系统 | |
WO2015103804A1 (zh) | 基于lms的信道均衡和频偏估计联合并行方法 | |
US20040151269A1 (en) | Receiver sampling in an ultra-wideband communications system | |
Papadias et al. | A space-time constant modulus algorithm for SDMA systems | |
CN104410593B (zh) | 基于判决反馈模型的数字符号非线性误差修正均衡方法 | |
CN103841064B (zh) | 一种基于采样率数据的siso均衡方法及其均衡器 | |
CA2601383A1 (en) | Channel estimation enhanced lms equalizer | |
US7289555B2 (en) | Method and system for signal processing using vector output from scalar data | |
CN101599930B (zh) | 高速并行均衡器及均衡方法 | |
WO2000044141A1 (en) | Mlse using look up tables for multiplication | |
CN111817992B (zh) | 数字化通道群延时均衡器及其实现方法、装置 | |
EP1520384B1 (en) | Method and system for concurrent estimation of frequency offset and modulation index | |
CN110213184B (zh) | 一种基于修正代价函数的自适应信道盲均衡方法 | |
Yang et al. | A vector constant modulus algorithm for shaped constellation equalization | |
CN106254286B (zh) | 高速并行处理的自适应盲均衡方法及装置 | |
CN109818894A (zh) | 多径信道下gmsk信号检测方法及检测装置 | |
US20190182082A1 (en) | Partially disjoint equalization and carrier recovery | |
JP2885618B2 (ja) | 適応受信機 |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100519 Termination date: 20110302 |