CN1339215A - 利用查表代替乘法运算的最大似然序列估计 - Google Patents

利用查表代替乘法运算的最大似然序列估计 Download PDF

Info

Publication number
CN1339215A
CN1339215A CN99816504A CN99816504A CN1339215A CN 1339215 A CN1339215 A CN 1339215A CN 99816504 A CN99816504 A CN 99816504A CN 99816504 A CN99816504 A CN 99816504A CN 1339215 A CN1339215 A CN 1339215A
Authority
CN
China
Prior art keywords
branch metric
value
precomputation
equalizer
filter
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
Application number
CN99816504A
Other languages
English (en)
Other versions
CN1145330C (zh
Inventor
S·陈纳克舒
K·赞吉
G·E·博顿利
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.)
Ericsson Inc
Original Assignee
Ericsson 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 Ericsson Inc filed Critical Ericsson Inc
Publication of CN1339215A publication Critical patent/CN1339215A/zh
Application granted granted Critical
Publication of CN1145330C publication Critical patent/CN1145330C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • 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
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03184Details concerning the metric
    • H04L25/03197Details concerning the metric methods of calculation involving metrics
    • 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
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03248Arrangements for operating in conjunction with other apparatus
    • H04L25/03299Arrangements for operating in conjunction with other apparatus with noise-whitening circuitry

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Abstract

本发明包括存在符号干扰失真时对以M元(M-ary)调制所调制的无线电信号进行解调。本发明提出一种方法,用于减少最大似然序列估计(MLSE)均衡器对以M元调制所调制的信号进行解调时所需乘法运算的次数。在本发明的示例性实施例中,通过预计算确定分支量度所需的某些值并把这些值存储在乘积表中来减少乘法运算的次数。需要进行分支量度计算时,无论是Euclidean分支量度计算还是Ungerboeck分支量度计算,某些乘法运算都可以用简单的查表操作来代替,因而增加了系统的功效和速度。任何对以M元调制所调制的信号进行解调的接收机都可以使用本发明。由此得到的解调器的复杂度低于现有解调器的复杂度。

Description

利用查表代替乘法运算的最大似然序列估计
发明背景:
数字无线通信系统被用于在多个位置间传输各种信息。由于数字通信,信息为通信目的而转换为数字或二进制形式,称之为比特。发射机把这类比特流映射为调制符号流,数字接收机可以检测到这些符号流并把它们映射回比特和信息。
在数字通信中,无线环境中有很多阻碍成功通信的困难。一个困难是由于信号的多径传输而导致的信号幅度衰落。结果,到达接收机天线的信号影像的相位并不相同。这种类型的衰落称为Rayleigh衰落或快衰落。当信号发生衰落时,信噪比变低,导致通信链路的质量下降。
当多条信号路径长度差别很大时会产生另一个问题。这种情况下会发生时间扩散,其中多个衰落的信号影像在不同的时间到达接收机天线,因而造成信号回波或辐射。这就引起符号间干扰(ISI),其中一个符号的回波会干扰后面的符号。
在接收机中,因为相干解调提供最好的性能,所以希望得到相干解调。这需要知道多径信道。在许多无线应用中,由于发射机移动、接收机移动和/或散射体移动,信道是随时间变化的。因而,有必要对随时间变化的多径信道进行跟踪。
为了提供多径信号的相干解调,可以使用最大似然序列估计(MLSE)均衡器。这种均衡器考虑传输符号序列的各种假设,并通过色散信道模型,确定最适合接收数据的假设。这可以利用Viterbi算法来有效地实现。这种均衡技术为本领域的技术人员所熟知,而且在J.C.Proakis 1989年的“数字通信”中可以查阅到该技术。
可以用简单的例子来解释传统的MLSE均衡器。假设发射机发射符号流s(n),符号取值为“+B”或“-B”,分别对应于比特值“0”或“1”。这种比特流通过二进制位移键(BPSK)调制。在接收机中,接收信号经过滤波,放大,并使用I和Q载波混频,然后每个符号周期(T)抽样一次,从而产生接收信号流r(n)。在此示例中,干扰信道包括两种射线组成,即一种主射波和一种反射波,反射波比主射波晚T秒到达,其中T是符号周期。随后,接收信号可以模拟为
r(n)=c0s(n)+c1s(n-1)+η(n)           (1)
其中c0和c1是复信道的抽头值,η(n)是加性噪声或干扰。
在MLSE均衡器中,第n次迭代时会有两个不同的先前“状态”0和1,分别对应于先前符号的两个可能值:
Figure A9981650400051
Figure A9981650400052
每个先前状态都有一个累加量度与之相关,是从以前迭代中累加得到的,分别产生先前状态0和1的累加量度A0(n-1)和A1(n-1)。
可能还有两个当前状态对应于s(n)的两个可能值。先前状态和当前状态的每个可能的配对对应于假设序列{sh(n-1),sh(n)}。对每个这样的假设,在n时刻都有一个相应的假设接收信号值:
rh(n)=c0sh(n)+c1sh(n-1)               (2)
而且,对每个这样的假设,有一个相应的“分支”量度,表示为
Mh(n)=|r(n)-rh(n)|2                   (3)
在示例中,有四个可能的假设序列,这四个序列可以表示为h∈(00,01,10,11),如表1所示:
先前状态 当前状态 h
0(B) 0(B) 00
1(-B) 0(B) 10
0(B) 1(-B) 01
1(-B) 1(-B) 11
表1:状态命名惯例
每个可能的当前状态的候选量度是相应分支量度和sh(n-1)相关的先前累加量度的和。每个当前状态都有两个可能的先前状态。对于每个当前状态,选择能得到最小候选量度的先前状态作为前驱状态,而候选量度成为那个当前状态的累加量度。
因而,对于当前状态0(即s(n)=B),有两个假设序列{sh(n-1)=B,sh(n)=B)或(sh(n-1)=-B,sh(n)=B),表示为假设h=00和h=10。这产生了当前状态0的两个候选量度:
C00(n)=A0(n-1)+M00(n)    (4)
C10(n)=A1(n-1)+M10(n)    (5)
这两个候选量度中较小的一个作为当前状态0的累加量度。相应的先前状态成为n时刻状态0的前驱状态。
1.如果C10(n)<C00(n),A0(n)=C10(n)
2.如果C00(n)<C10(n),A0(n)=C00(n)
当前状态1(即s(n)=-B)使用类似的过程,在本例中,两个假设是h=01和h=11。
对于每个n,必须根据等式(2)计算四个假设接收信号的值。计算每个假设值rh(n)包括了两次乘法运算。每次乘法运算包括一个信道抽头与一个假设符号值的乘积。
如果设B为+1,乘法就成了简单的符号变化。但是,使用幅度调制时就不可能了。还有,假设符号值归一化为
Figure A9981650400061
(即8PSK调制)。在这种情况下,不能把必要的乘法当作简单的符号变化。
增加数字信道吞吐量的一种方法就是允许每个发射符号的取值多于两个。一般来说,每个符号选自M个可能值的集合{B1,B2,...,BM),从而增加了每个符号的比特数。MLSE均衡器传统实现的复杂度随所允许符号数(即M)呈指数级增长。特别地,如果有L个信道抽头,均衡器的传统实现所需乘法的次数与ML-1成正比(比例常数大于1)。
即使M和L的值(例如分别为8和5)适中,在低功耗通用DSP中实现均衡器也是不现实的。而功率敏感型应用的趋势是从DSP向ASIC转移。ASIC可以最佳用于以低功耗和小尺寸(即少的门数)执行一个任务。但用ASIC实现乘法器比实现加法或减法器更为昂贵(根据功耗和门数)。因而,非常希望避免与MLSE传统实现相关的大量乘法运算。
因为所需乘法的次数有增加的趋势,而技术正从乘法操作容易的DSP向乘法相对昂贵的ASIC转移,因而希望减少MLSE实际执行的乘法次数。
发明概述:
本公开涉及有符号间干扰失真的M元(M-ary)调制的无线电信号的解调。本发明提出一种方法,它减少实现用于M元调制所调制的信号的MLSE均衡器所需乘法运算的次数。
在本发明的示例性实施例中,通过预计算分支长度分支量度判决所需的某些值并把这些预计算的值存储在乘积表中来减少乘法运算的次数。计算分支量度时,无论是Euclidean分支量度计算还是Ungerboeck分支量度计算,某些乘法运算都可以用简单的查表操作来代替,从而减小了均衡器的功耗和尺寸。
任何对M元调制所调制的信号进行解调的接收机都可以使用本发明来实现。由此得到的解调器的复杂度要低于现有解调器的复杂度。
附图简述:
通过参考附图对最佳实施例的下述描述,本发明的上述目的和特征将更加清楚,附图中:
图1示出可以实现按照本发明的技术的无线通信系统;
图2示出图1中基带处理器的功能方框图;
图3示出图1中传输功能的功能方框图;
图4示出用于计算Euclidean分支量度的系统,其中可以实现按照本发明的技术;
图5示出用于根据现有技术计算Euclidean量度的滤波器的功能方框图;
图6示出根据本发明一个实施例的Euclidean量度的乘积表;
图7示出根据本发明一个实施例的用于计算Euclidean量度的滤波器的功能方框图;
图8示出用于计算Ungerboeck分支量度的系统,其中可以实现按照本发明的技术;
图9示出根据现有技术计算Ungerboeck量度的传统滤波器的功能方框图;
图10示出根据本发明一个实施例的Ungerboeck量度的乘积表;和
图11示出根据本发明一个实施例计算Ungerboeck量度的滤波器的功能方框图。
详细说明:
图1示出无线通信系统,其中实现本发明。在图1中,提供无线通信系统的无线电发射机和无线电接收机。无线通信系统可以工作于FDMA、TDMA、或CDMA方式,或工作于它们的组合方式。发射机有数字符号发生器102,它可以接收携带信号101的信息并产生相应的数字符号序列S。根据已知技术,数字符号S经过数模(D/A)转换、调制、脉冲成形滤波和放大后,作为模拟信号Y由数字发射机103发射。
信号Y经过无线信道传输,由接收机的天线104截取。天线104还截取了热噪声n。
无线电单元105根据已知方法对接收信号进行放大、下变频和滤波以产生模拟输出。该模拟输出耦合到A/D转换器106,它把模拟信号转换为接收信号抽样流r(kTs),其中Ts是抽样周期,k是整数计数器。抽样周期Ts可以小于符号周期T。处理器107收集接收和处理接收信号抽样流并产生发射数字符号流S的估计。在后面的说明中,传输功能109被用来表示整个通过数字发射机103、无线电传输信道105和A/D106的信号路径。
传输功能109产生接收信号抽样流r(kTs)并把它送入处理单元107,在处理单元107中根据本发明对其进行处理。
图2示出基带处理单元107的功能方框图。接收信号抽样流r(kTs)耦合到信号预处理器或同步单元206,其中接收信号抽样流根据已知技术与已知定时/同步序列相关联。对于符号间隔解调,如果抽样周期Ts小于符号周期T,信号预处理器206会对接收信号抽样流r(kTs)进行抽选,以便每个符号产生一个抽样,指定为r(n)。对于分隔解调(fractionally spaced demodulation),每个符号产生一个以上的抽样。
估计电路202根据已知技术产生用于模仿无线传输信道的信道抽头估计c(τ)。这可能包括一个初始信道估计周期,然后是一个跟踪周期。信道抽头估计c(τ)耦合到分支量度处理器203的输入端。分支量度处理器203耦合到提供数字符号流估计S的序列估计处理器204。
图3更为详细地示出传输功能109,其中为了简便,接收天线的数目限制为一。本领域的技术人员将知道本发明还可以用于两个或两个以上天线的情况。在图3中,符号序列S输入到发射模拟信号Y的数字发射机103。模拟信号Y通过无线信道传到无线电单元105。无线信道301会引入衰落和时间扩散。同时接收普遍存在的热噪声。无线电单元105根据已知技术对接收信号进行放大、下变频和滤波以产生模拟信号。该模拟信号耦合到A/D 106,并由A/D 106把该模拟信号转变为接收信号抽样r(kTs)。
在MLSE均衡器中,要有效考虑所有可能的传输符号序列S。在一种实现中,由信道抽头估计c(τ)对假设符号值sh(n)例进行过滤以产生假设接收抽样rh(n)。假设的rh(n)和实际的r(n)接收信号抽样流之间的差别,称为假设误差e(n),可以表示一个特别假设的好坏程度。假设误差的平方可以作为评价某个特别假设的量度。对不同的假设进行量度累加,用于确定使用序列估计算法时哪个假设较好。这可以利用Viterbi算法有效地实现。在G.Forney,Jr.的“维特比算法”(Proceedings of the IEEE,vol.61,no.3,March 1973,pp.267-278)一文中可以查阅到对Viterbi算法的描述。本领域的技术人员将明白,也可以使用其它序列估计算法。
在MLSE均衡器中,存在与不同的发射序列假设sh(n)相关的状态。在给定的迭代中,存在先前状态,并且每一个先前状态都与一个累加量度相关。每一对先前状态和当前状态产生一个分支量度Mh(n)。此分支量度Mh(n)和先前累加量度相加得到当前状态的候选量度。对于每个当前状态,选择得到最小候选量度的先前状态作为前驱状态,而最小的候选量度成为当前状态的累加量度。分支量度可以表示为:
Mh(n)=|r(n)-rh(n)|2              (6)其中 r h ( n ) = Σ k = 0 N 1 - 1 c ( k ) s h ( n - k ) - - - ( 7 )
由c(τ)指定信道抽头估计,其中τ是延迟(即,τ=0是主射波,τ=1是第一反射波,等等)。Nl是估计的信道抽头数。对于每个n,必须根据等式(7)计算假设接收信号rh(n)。大多数情况下,每个假设接收值rh(n)的计算包括Nl次复数乘法(即4Nl次实数乘法)。每次复数乘法包括一个估计信道抽头和假设发射信号的乘积。
图4示出计算Euclidean分支量度Mh(n)的系统,其中可以实现按照本发明的技术。在图4中,符号序列发生器410产生的符号假设序列sh(n)耦合到滤波器400以产生假设接收抽样rh(n)。假设接收抽样rh(n)和实际接收信号抽样流r(n)之间的差是假设误差e(n)。单元403执行假设误差的平方运算以产生分支量度Mh(n)。
图5示出用于计算Euclidean分支量度的传统滤波器的功能方框图。在图5中,计算每个rh(n)需要执行Nl次复数乘法。如果发射符号限制在M个可能值的集合{B1,B2,...,BM}中,本发明提供一种实现避免任何乘法运算的滤波器400的方法。
对于Euclidean量度,假定每个假设发射符号sh(n)都在集合{B1,B2...,BM}中。等式(7)中加法的第一项在集合{Bjc(0)}j=0 M中。该集合中的所有项均被预计算并存储在M×Nl表(见图6)的第一列中。
此表的第j列,对应于估计信道抽头c(j-1),存储了(sh(·)c(j-1))的所有可能值,即{B1,c(j-1),B2c(j-1),...,BMc(j-1)}。通过所述乘积表(见图7)中合适项的简单相加便可以计算每一个假设接收值rh(n)。
作为一个具体的例子,考虑两抽头信道(c0,c1)和8PSK调制,即。这种情况下,对应于信道抽头c0,图6所示表的第一列有8项: [ e j 2 π 8 0 c 0 , e j 2 π 8 c 0 , e j 2 π 8 2 c 0 . . . , e j 2 π 8 7 c 0 ] - - - ( 8 ) 类似地,对应于信道抽头c1,图6中的表的第二列有8项: [ e j 2 π 8 0 c 1 , e j 2 π 8 c 1 , e j 2 π 8 2 c 1 , · · · , e j 2 π 8 7 c 1 ] - - - ( 9 )
可以使用对称来减少存储要求。首先,使用负对称可以使存储的项减半。基本上, e j 2 π 8 l + π c 0 = - e j 2 π 8 l c 0 - - - ( 10 )
因而,只要存储l=0,1,2,3的值。此外, e j 2 π 8 2 c 0 = je j 2 π 8 0 c 0 - - - ( 11 )
以及 e j 2 π 8 3 c 0 = je j 2 π 8 ( 1 ) c 0 - - - ( 12 )
因而,通过对实部和虚部进行交换以及对新实部取反,只需存储l=0和1的值。对于l=0,不必进行乘法运算,因为 e j 2 π 8 0 c 0 = c 0 - - - ( 13 )
因而,可以只存储c0
Figure A9981650400124
并利用取反和对实部和虚部进行
交换的逻辑来导出其它值。这些节省来自于8PSK信号星座图中的象限对称。
由等式(3)的Euclidean量度Mh(n),经过两步可以得到Ungerboeck量度。第一步是展开Mh(n)以得到
Mh(n)=A(n)+B(n)+C(n)+D(n)              (14)其中
A(n)=|r(n)|2                           (15) B ( n ) = - 2 Re { r ( n ) c 0 * s h * ( n ) } - 2 Re { r ( n ) c 1 * s h * ( n - 1 ) } - - - ( 16 )
C(n)=|c0|2|sh(n)|2+|c1|2|sh(n-1)|2    (17) D ( n ) = 2 Re { c 0 c 1 * s h ( n ) s h * ( n - 1 ) } . - - - ( 18 ) Ungerboeck方法省去了对所有假设都相同的项A(n)。第二步是组合不同次迭代中与sh *(n)成比例的项。在第(n+1)次迭代中,该项变为: B ( n + 1 ) = - 2 Re { r ( n + 1 ) c 0 * s h * ( n + 1 ) } - 2 Re { r ( n + 1 ) c 1 * s h * ( n ) } - - - ( 19 )
  C(n+1)=|c0|2|sh(n+1)|2+|c1|2|sh(n)|2            (20) D ( n + 1 ) = 2 Re { c 0 c 1 * s h ( n + 1 ) s h * ( n ) } . - - - ( 21 )
因而,两次迭代中都有与sh *(n)成比例的项。通过定义一个新的量度可以组合这些项,Mh 1(n),其表示如下: M h 1 ( n ) = Re { s h * ( n ) [ - 2 z ( n ) + q 0 s h ( n ) + q 1 s h ( n - 1 ) ] } - - - ( 22 )
其中 2 ( n ) = c 0 * r ( n ) + c 1 * r ( n + 1 ) - - - ( 23 )
  q0=|c0|2+|c1|2                                  (24) q 1 = 2 c 0 c 1 * . - - - ( 25 )
这里,q0和q1称为s参数。
因而对于Ungerboeck量度,分支量度Mh 1(n)定义为: M h l ( n ) = Re { s h * ( h ) [ - 2 t ( n ) + t h ( n ) ] } , - - - ( 26 )
其中 t h ( n ) = Σ k = 0 N 1 - 1 q ( k ) s h ( n - k ) . - - - ( 27 )
图8示出计算Ungerboeck分支量度的系统,其中可以实现按照本发明的技术。在图8中,假设符号sh(n)与脉冲响应q(τ)一起耦合到滤波器600以产生th(n)。q(τ)由估计信道抽头c(τ)’确定;因而,q(τ)’不取决于假设符号sh(n)。如图9所示,用于实现图8的滤波器的技术的当前状态需要Nl个复数乘法。
假定每个假设发射符号sh(n)都在集合{B1,B2……,BM}中。等式(27)中加法的第一项在集合{biq(0)}i=0 M中。所述集合中的所有项都可以预计算并存储在M×Nl表(见图10)的第一列中。
对应于q(j-1),所述表的第j列存储了sh(·)q(j-1)所有可能的值,即{B1q(j-1),B2q(j-1),...,BMq(j-1)}。通过所述乘积表(见图11)中合适项的简单相加便可以计算每一个假设th(n)。
如前所述,可以利用象限对称把存储在每一列中项的数目减少到2。
部分Ungerboeck量度也包括象前面部分提出的qi的项;因此,本发明可以用来来减少用于实现使用部分Ungerboeck量度的解调器所需的乘法次数。在给Bottomly的美国专利5499272中公开了部分Ungerboeck量度。
对于使用Euclidean量度的多信道MLSE,通过几个信道接收发射符号流。如果有两个信道,就有分别对应于信道a和b的两个接收信号ra(n)和rb(n):
       ra(n)=ca0s(n)+ca1s(n-1)+ηa(n)      (28)rb(n)=cb0s(n)+cb1s(n-1)+ηb(n)。    (29)
对应于每个接收信号和每个假设符号序列,有一个分支错误信号:
       eah(n)=ra(n)-rah(n)                   (30)
       ebh(n)=rb(n)-rbh(n),                 (31)其中rah(n)和rbh(n)是信道a和b的假设接收信号(两者都对应于相同的假设序列{sh(n-1),sh(n-1)}):
   rah(h)=ca0sh(n)+Ca1sh(n-1)      (32)
   rbh(n)=cb0sh(n)+Cb1sh(n-1)。    (33)注意信道a的一个分支量度是eah(n)eah *(n),信道b的一个分支量度是ebh(n)ebh *(n)。
两个信道输出的联合均衡的不同分支量度可以通过以不同的方式组合eah(n)和ebh(n)得到。联合均衡问题中几种引起关注的分支量度可以表示为如下形式:
       Mh(n)=Eh(n)HQEh(n),        (34)其中 E h ( n ) = e ah ( n ) e bh ( n ) - - - ( 35 ) Q是2×2的加权矩阵。
利用两个预计算表可以有效的计算等式(34)给出的所述形式的任何分支量度。可以通过利用两个预计算表计算每个假设接收信号rah(n)或rbh(n)而不需执行乘法操作。
本领域的技术人员可以理解接收多于两个信道时如何使用本发明。类似地,本领域的技术人员可以理解符号取M个可能值之一时如何使用本发明。
在量度结合时,加权矩阵Q是对角矩阵。该矩阵对角线上的项是每个信道的加权系数。
在干扰抑制结合(IRC)中,加权矩阵Q是减损的相关矩阵的逆矩阵。特别地,减损矢量i(n)可以定义为: i ( n ) = η a ( n ) η b ( n ) , - - - ( 36 ) 令Rηη=E{i(n)i(n))H}。在IRC中,加权矩阵Q简单地为Q=Rηη -1。美国专利5680419中更详细地介绍了IRC,通过引用将其结合在此。
对于使用Ungerboeck量度的多信道MLSE,考虑具有两个信道的情况;因此,对于每个n都有由等式(28)至(29)给出的两个接收抽样值ra(n)和rb(n)。
在这种情况下,Ungerboeck量度由下式给出: M h 1 ( n ) = Re { s h * ( n ) [ - 2 z ( n ) + q 0 s h ( n ) + q 1 s h ( n - 1 ) ] } - - - ( 37 ) 其中z(n)是从两个接收信号和两个信道抽头得到的,(q0,q1)是从两个信道抽头和减损相关矩阵的逆矩阵得到的(见美国专利5680419)。值得注意的是(q0,q1)固定值,并不取决于假设符号。
从等式(37)看出,多信道情况下计算Mh 1(n)显然需要计算:
                th(n)=q0sh(n)+q1sh(n-1)。     (38)计算th(n)通常需要两次乘法运算。通过预计算可能的符号值s的值{sq0,sq1}并将其存储在存储器的表中可以避免这些乘法运算。通过此预计算表中恰当项的相加可以计算出任何假设th
本领域的技术人员可以理解接收多于两个信道时如何使用本发明。类似地,本领域的技术人员可以理解符号可以取M个可能值之一时如何使用本发明。
对于分隔MLSE均衡的情况,(T/2)分隔MLSE均衡器在长度为T的每个抽样间隔中要接收两个抽样。用ra(n)表示接收信号的偶数抽样,用rb(n)表示接收信号的奇数抽样。注意ra(n)和rb(n)都是符号间隔的。在分隔MLSE均衡器的具体实现中,ra(n)和rb(n)被视为两个独立的接收信号(由相同的发射符号流产生)。根据此示例性实施例,可以用上一段所述双输入MLSE均衡器检测发射符号流。上面公开的减少MLSE均衡器分集复杂度的技术明显可用于分隔MLSE均衡器的这种实现。
本领域的技术人员可以理解本发明如何用于其它形式的分隔均衡。例如,在分隔均衡器的一个特定实现中,对每个子抽样序列使用预白化滤波器(见Hamied和Stuber,“分隔MLSE接收机”,IEEE1995)。这个特定的表述非常类似于量度结合。在分隔均衡器的另一个表述中,使用Ungerboeck量度;因此,也可以用本发明减小这种均衡器的复杂度。本领域的技术人员可以理解每个符号接收两个以上抽样时如何使用本发明。
对于使用时变信道的MLSE,每次迭代都考虑不同的信道抽头集合。对于Euclideau量度,假设接收信号可以表示为:
rh(n)=c0(n)sh(n)+c1(n)sh(n-1)              (39)
在这种情况下,在第n次迭代中,该表可以被更新并用于计算本次迭代中所有不同假设接收值。第“n+1”次迭代计算一个新表,等等。
本领域的技术人员可以理解时变信道中使用Ungerboeck时如何使用本发明。每次qi更新时表也要更新。
对于解调利用偏置M元调制(offset M-ary modulation)所调制的信号,可以修改本发明来进一步减少解调器的复杂度。在偏置M元调制中,偶数“n”的发射符号是纯实数,奇数“n”时的发射符号是纯虚数:
    s(n)=f*Bi  i∈{1,...,M},                     (40)其中每个Bi是纯实数。因此,sh(n)可以取2M个值{Bi,jBi}i=1 M中的一个。
但是,通过对接收信号r(n)预旋转j-n,我们可以减少解调纯实符号的解调问题,其中每个实符号都在集合{Bi}i=1 M中。使用预旋转,我们只需要存储信道抽头和每个{Bi}i=1 M的乘积。
本领域的普通技术人员可以理解,在脱离这些特定细节的其它实施例中可以实践本发明。在其它例子中,省去对众所周知的方法、设备、电路的详细说明,以免不必要的细节使本发明的说明不明显。
在不脱离后附权利要求书所述的本发明精神或范围的条件下,本领域的技术人员可以设计上述技术的多种变体和组合。

Claims (13)

1.一种用于确定中分支量度的一种方法,其中所述最大似然估计均衡器接收至少一个以M元(M-ary)调制所调制的天线信号,所述方法包括以下步骤:
预计算等于复数和假设符号值的乘积的值;
把所述预计算值存储在乘积表中;
把来自所述乘积表的选定预计算值相加,以产生一个结果;以及
利用所述结果确定所述分支量度。
2.权利要求1的方法,其特征在于所述复数对应信道系数。
3.权利要求1的方法,其特征在于所述复数对应s参数。
4.一种最大似然序列估计均衡器中的滤波器,用于产生用于确定分支量度的假设接收信号抽样,其中所述最大似然序列估计均衡器对至少一个以M元调制所调制的接收无线电信号进行解调,所述滤波器包括:
存储器,用于存储乘积表,所述乘积表具有等于不同次迭代的信道抽头估计与假设符号值乘积的预计算值;和
加法器,用于把来自所述乘积表的选定项进行相加以产生假设接收信号抽样。
5.一种最大似然序列估计均衡器中用于M元调制的滤波器,所述滤波器包括:
用于预计算用于确定分支量度的多个可能值的装置;
用于存储所述多个预计算的可能值的存储器;和
用于组合来自所述存储器的选定预计算值的装置。
6.根据权利要求5的滤波器,其特征在于所述分支量度是Ungerboeck分支量度。
7.根据权利要求5的滤波器,其特征在于所述分支量度是Euclidean分支量度。
8.根据权利要求5的滤波器,其特征在于所述分支量度是部分Ungerboeck分支量度。
9.一种用于在对M元调制信号进行解调的最大似然序列估计均衡器中计算分支量度的方法,所述方法包括如下步骤:
预计算用于所述分支量度计算的多个可能值;
把所述多个预计算的可能值存储在存储器中;
把来自所述存储器的选定预计算值相加;和
利用所述相加的选定预计算值计算所述分支量度。
10.根据权利要求9的方法,其特征在于所述分支量度是Ungerboeck分支量度。
11.根据权利要求9的方法,其特征在于所述分支量度是部分Ungerboeck分支量度。
12.根据权利要求9的方法,其特征在于所述分支量度是Euclidean分支量度。
13.一种用于在多信道最大似然序列估计(MLSE)均衡器中计算分支量度的方法,其中所述多信道最大似然序列估计均衡器对M元调制信号进行解调,所述方法包括如下步骤:
预计算所述多信道MLSE中用于分支量度计算的各个信道的多个可能值;
把所述各个信道的多个可能值存储在分开的乘积表中;
将来自所述分开的乘积表的选定值相加;以及
利用所述相加的选定值计算所述分支量度。
CNB998165042A 1999-01-26 1999-12-29 利用查表代替乘法运算的最大似然序列估计 Expired - Fee Related CN1145330C (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US09/237,356 1999-01-26
US09/237,356 US7099410B1 (en) 1999-01-26 1999-01-26 Reduced complexity MLSE equalizer for M-ary modulated signals
US09/237356 1999-01-26

Publications (2)

Publication Number Publication Date
CN1339215A true CN1339215A (zh) 2002-03-06
CN1145330C CN1145330C (zh) 2004-04-07

Family

ID=22893386

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB998165042A Expired - Fee Related CN1145330C (zh) 1999-01-26 1999-12-29 利用查表代替乘法运算的最大似然序列估计

Country Status (8)

Country Link
US (1) US7099410B1 (zh)
EP (1) EP1147641A1 (zh)
JP (1) JP2003518340A (zh)
CN (1) CN1145330C (zh)
AR (1) AR025511A1 (zh)
AU (1) AU2596000A (zh)
TR (1) TR200102153T2 (zh)
WO (1) WO2000044141A1 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101213762A (zh) * 2005-06-10 2008-07-02 艾利森电话股份有限公司 无线通信接收器中的减损相关估算的方法和装置
CN1819690B (zh) * 2006-03-31 2011-03-23 北京邮电大学 一种应用于fdd/ofdm系统的下行干扰抑制方法
CN104503711A (zh) * 2014-11-17 2015-04-08 杭州先临三维科技股份有限公司 一种3d打印的自适应分层方法
CN105573712A (zh) * 2014-10-31 2016-05-11 Arm有限公司 用于计算重复迭代和的结果的装置、方法和程序
CN115242584A (zh) * 2022-07-21 2022-10-25 苏州大学 一种优化基于查找表的mlse算法复杂度的方法和装置

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6115406A (en) 1999-09-10 2000-09-05 Interdigital Technology Corporation Transmission using an antenna array in a CDMA communication system
US6278726B1 (en) 1999-09-10 2001-08-21 Interdigital Technology Corporation Interference cancellation in a spread spectrum communication system
DE50203032D1 (de) * 2001-11-09 2005-06-09 Infineon Technologies Ag Viterbi-entzerrung unter verwendung eines tabellenspeichers zur bereithaltung von rekonstruierten signalwerten für die berechnung von übergangsmetriken
WO2003075528A1 (en) * 2002-03-07 2003-09-12 Telefonaktiebolaget L M Ericsson (Publ) Soft value calculation for multilevel signals
GB2392066B (en) 2002-08-16 2005-11-09 Toshiba Res Europ Ltd Equaliser apparatus and methods
US7590175B2 (en) * 2003-05-20 2009-09-15 Rambus Inc. DFE margin test methods and circuits that decouple sample and feedback timing
US7627029B2 (en) 2003-05-20 2009-12-01 Rambus Inc. Margin test methods and circuits
US7292650B2 (en) * 2003-07-01 2007-11-06 Silicon Integrated Systems Corp. OFDM receiver and metric generator thereof
GB0420183D0 (en) 2004-09-10 2004-10-13 Ttp Communications Ltd Metric calculation utilising pre-stored values
US20070127608A1 (en) * 2005-12-06 2007-06-07 Jacob Scheim Blind interference mitigation in a digital receiver
CN101296350B (zh) * 2007-04-25 2011-07-27 普诚科技股份有限公司 解交错器的编译码方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2215567B (en) 1988-03-05 1992-11-18 Plessey Co Plc Improvements in or relating to equalisers
JP2960436B2 (ja) 1989-06-26 1999-10-06 株式会社日立製作所 非線形データ伝送システム用受信器
JP3258067B2 (ja) 1992-03-25 2002-02-18 株式会社東芝 最尤系列推定方式
JP3560991B2 (ja) * 1993-09-20 2004-09-02 株式会社東芝 適応型最尤系列推定装置
US5488635A (en) * 1993-10-29 1996-01-30 General Electric Company Low complexity adaptive equalizer radio receiver employing reduced complexity branch metric calculation
US5499272A (en) 1994-05-31 1996-03-12 Ericsson Ge Mobile Communications Inc. Diversity receiver for signals with multipath time dispersion
US5481572A (en) 1994-08-02 1996-01-02 Ericsson Inc. Method of and apparatus for reducing the complexitiy of a diversity combining and sequence estimation receiver
US5680419A (en) 1994-08-02 1997-10-21 Ericsson Inc. Method of and apparatus for interference rejection combining in multi-antenna digital cellular communications systems
US5949796A (en) * 1996-06-19 1999-09-07 Kumar; Derek D. In-band on-channel digital broadcasting method and system
US6108517A (en) * 1997-07-28 2000-08-22 Ericsson Inc. Methods and apparatus for joint demodulation of adjacent channel signals in digital communications systems

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101213762A (zh) * 2005-06-10 2008-07-02 艾利森电话股份有限公司 无线通信接收器中的减损相关估算的方法和装置
CN101213762B (zh) * 2005-06-10 2013-02-06 艾利森电话股份有限公司 无线通信接收器中的减损相关估算的方法和装置
CN1819690B (zh) * 2006-03-31 2011-03-23 北京邮电大学 一种应用于fdd/ofdm系统的下行干扰抑制方法
CN105573712A (zh) * 2014-10-31 2016-05-11 Arm有限公司 用于计算重复迭代和的结果的装置、方法和程序
CN105573712B (zh) * 2014-10-31 2020-09-04 Arm 有限公司 用于计算重复迭代和的结果的装置、方法和存储介质
CN104503711A (zh) * 2014-11-17 2015-04-08 杭州先临三维科技股份有限公司 一种3d打印的自适应分层方法
CN104503711B (zh) * 2014-11-17 2017-10-20 杭州先临三维科技股份有限公司 一种3d打印的自适应分层方法
CN115242584A (zh) * 2022-07-21 2022-10-25 苏州大学 一种优化基于查找表的mlse算法复杂度的方法和装置
CN115242584B (zh) * 2022-07-21 2023-12-01 苏州大学 一种优化基于查找表的mlse算法复杂度的方法和装置

Also Published As

Publication number Publication date
EP1147641A1 (en) 2001-10-24
WO2000044141A1 (en) 2000-07-27
JP2003518340A (ja) 2003-06-03
AR025511A1 (es) 2002-12-04
TR200102153T2 (tr) 2001-12-21
US7099410B1 (en) 2006-08-29
AU2596000A (en) 2000-08-07
CN1145330C (zh) 2004-04-07

Similar Documents

Publication Publication Date Title
CN1078410C (zh) 用于多径时间离散信号的分集接收机
CN1113482C (zh) 利用分集信号矩阵处理降低码元间干扰的方法和设备
CN1145330C (zh) 利用查表代替乘法运算的最大似然序列估计
CN1193563C (zh) 多载波接收机和含此类接收机的发射系统
CN1122394C (zh) 联合解调码分多址访问信号的系统和方法
CN1248515C (zh) 自适应均衡设备和方法
CN1063890C (zh) 用于解调下行链路码分多址信号的方法和系统
CN1284317C (zh) Ofdm发送和接收设备
CN1245607A (zh) 通过使用传输媒体响应估值进行数字符号检测的方法和设备
CN1689292A (zh) 具有虚拟并联均衡器的通信接收机
CN1502175A (zh) 使用频道关联矩阵快速傅利叶转换的低复杂性数据检测
CN1820424A (zh) 带有混合均衡器和rake接收机的接收设备以及相应的接收方法
CN1280995C (zh) 使用最大似然序列估计于无线接收的均衡器装置与方法
CN1348642A (zh) 使用发射分集技术的方法和系统
CN101317411A (zh) 相位噪声消除ofdm接收器
CN1898925A (zh) 用于格雷-映射qam的快速软值计算方法
CN1172498C (zh) 自适应最大似然序列检测的方法
CN1063600C (zh) 用于码分多址通信系统接收机的检波器
CN1207879C (zh) 实现维特比平衡的度量增量预先计算法
CN1239033C (zh) 无线传输中基于离散付氏变换的空时联合检测装置和方法
CN1294805A (zh) 数字均衡方法和一台实施这种方法的无线电通信接收机
US7386035B2 (en) Method of compensating for energy loss and eliminating inter-symbol interference and inter-chip interference and rake receiver for WLAN adopting the same
CN1599364A (zh) 用于突发通信的数据均衡方法
CN1355978A (zh) 分数间隔衰落无线信道的符号间隔估计和/或跟踪方法
US7760797B2 (en) Method of reducing a number of computations in an equalization process and apparatus for performing the same

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C06 Publication
PB01 Publication
C14 Grant of patent or utility model
GR01 Patent grant
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee