CN101136896A - 基于快速傅立叶变换的频域迭代均衡方法 - Google Patents
基于快速傅立叶变换的频域迭代均衡方法 Download PDFInfo
- Publication number
- CN101136896A CN101136896A CNA2007101326910A CN200710132691A CN101136896A CN 101136896 A CN101136896 A CN 101136896A CN A2007101326910 A CNA2007101326910 A CN A2007101326910A CN 200710132691 A CN200710132691 A CN 200710132691A CN 101136896 A CN101136896 A CN 101136896A
- Authority
- CN
- China
- Prior art keywords
- soft
- output
- frequency domain
- equilibrium
- sample value
- 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
Images
Landscapes
- Noise Elimination (AREA)
Abstract
基于快速傅立叶变换的频域迭代均衡方法可以有效实现频率选择性信道下的迭代均衡。该方法为:首先,对于接收序列,取出K+L个样值,将前L个样值叠加到最后的L个样值上,其中K表示每次均衡的发送符号块的长度,L为信道阶数;然后,初次均衡采用多项式展开的频域软输出均衡;解交织后,根据均衡输出的软信息进行软输入软输出解码;解码输出的软信息进行交织后,计算符号的统计量;利用反馈的符号均值进行块间干扰抵消,然后采用频域软输入软输出均衡进入下一轮迭代。该方法可以在没有块间保护的系统进行频域迭代均衡,可以采用快速傅立叶变换来实现,在提高系统性能的同时,可以大大降低迭代接收机的复杂度。
Description
技术领域
本发明是一种应用于无线通信系统的联合均衡和译码的方法,属于移动通信中的均衡技术领域。
背景技术
为适应未来发展,移动通信系统必须能够支持高数据速率。这不可避免地会遇到频率选择性摔落信道。在传统的单载波系统中,针对频率选择性摔落信道,首先进行均衡,通常采用时域的线性最小均方误差均衡器(Minimum MeanSquared Error,MMSE),或软输出维特比均衡器,然后再进行译码。这种传统的做法并不是最优方法,但是最优的联合均衡和译码的实现极其复杂。Turbo迭代接收机可以很好地逼近联合最优接收机。但是,目前,针对频率选择性信道,Turbo接收机中,均衡器的实现通常在时域实现,例如采用软输入软输出线性MMSE均衡器,或软输入软输出维特比均衡器。前者通常需要较大的矩阵求逆运算,后者的复杂度通常随着信道阶数、调制阶数的增加而指数增长。在初次均衡,我们做适当的近似,在频域实现线性MMSE均衡器,之后的迭代均衡中,我们首先进行块间干扰抵消,然后采用频域软输入软输出MMSE均衡,可以采用FFT(Fast Fourier Transform Algorithm)快速实现,这样均衡器复杂度可以大大降低。
发明内容
技术问题:本发明的目的是提供一种基于快速傅立叶变换的频域迭代均衡方法,即是一种单载波系统在频率选择性摔落信道下的低复杂度的迭代接收机的实现方法。
技术方案:本发明的基于快速傅立叶变换的频域迭代均衡方法的初次均衡采用FFT/IFFT(Fast Fourier Transform/Inverse Fast Fourier Transform)快速实现的多项式展开的线性MMSE均衡,之后的每次迭代采用FFT/IFFT快速实现的软输入软输出线性MMSE均衡。具体有以下几个步骤:
初次均衡:
1.)取出K+L个样值,表示为y,先抵消前一个块的干扰,然后将前L个样值叠加到最后的L个样值上,表示为r1;
2.)对信道的时域响应进行K点FFT变换,表示为λ=[λ1λ2…λk]T,并保存下来,其中T表示矢量的转置;
3.)计算下一个块的干扰加上叠加引入的噪声的相关矩阵R;
4.)采用频域实现的多项式展开软输出频域均衡,并对均衡输出做高斯近似,计算等效幅度以及等效噪声方差;
5.)进行软解调,解调输出经解交织后传递给软输入软输出译码器;
6.)译码器输出的软信息经过交织后,计算发送信号的统计量,包括均值s和方差,并对方差求平均,表示为v;
迭代均衡:
7.)对K+L个样值y,根据译码器获得的发送符号的信息,抵消前一个块和后一个块对当前块的干扰,然后将前L个样值叠加到最后的L个样值上,表示为r2;
8.)采用频域实现的线性MMSE软输入软输出均衡,并对均衡输出做高斯近似,计算等效幅度以及等效噪声方差;
9.)进行软解调,解调输出经解交织后传递给软输入软输出译码器;如果没有达到给定的迭代次数,则译码器输出的软信息经过交织后,计算发送信号的统计量,包括均值和方差,并对方差求平均,并跳转到步骤7)。
具体地,第4)步所述的多项式展开软输出频域均衡,它包括以下几个步骤:
4.1.)计算频域均衡所用到的循环矩阵∑,该循环矩阵的第一列为:
的K点非归一化的IFFT变换,其中σ2为噪声方差;
4.2.)计算多项式展开时的归一化因子γ:
4.3.)计算:(2-γ)r1-γR∑r1;
4.5.)将步骤4.4)的结果进行K点IFFT变换,并计算均衡器输出的有偏系数以及等效噪声方差,最后进行软解调。
第8)步所述的频域实现的线性MMSE软输入软输出均衡,它包括以下几个步骤:
8.1.)对接收矢量r2进行K点FFT变换;
8.2.)对反馈的均值信号s进行K点FFT变换并对其输出的每一个点乘以λk,其中k=1,2,...,K;
8.3.)将步骤8.1)的结果减去步骤8.2)的结果,并对结果进行使用系数的单点均衡:其中k=1,2,...,K;
8.4.)均衡之后的结果进行K点IFFT变换,并加上其中
8.5.)计算均衡器输出的有偏系数以及等效噪声方差,最后进行软解调。
有益效果:本发明的主要优点在于可以在没有加保护间隔的系统中实现频域迭代接收机,可以用快速傅立叶变换实现,从而降低迭代接收机的复杂度。同时该方法可以逼近复杂度较高的时域迭代均衡接收机,可以大大提高系统性能。
本发明提出的Turbo迭代均衡方法能用于单载波系统在频率选择性衰落信道下的均衡,这些系统包括目前广泛应用的全球移动通信系统(GSM)。由于本发明针对一般的频率选择性信道,该方法也适用于其它系统的均衡,例如码分多址接入系统在扩频比较低时,存在符号间干扰。
附图说明
图1是本发明迭代均衡器的初次均衡的过程。它包括6个子图:图1(a)表示计算r1及其相应的FFT变换,它包括块间干扰抵消器和叠加器;图1(b)计算信道的频域响应,完成时域信道参数向频域的变换;图1(c)计算相关矩阵R;图1(d)计算矩阵Σ;图1(e)计算多项式展的归一化因子γ;图1(f)是软输出均衡,它包括计算多项式展开后的矩阵与发送信号的相乘运算器、FFT变换器、频域单点均衡器、IFFT变换器和软解调器。经过上述步骤,完成初次软输出均衡。
图2是本发明第二次以及之后的迭代均衡中所用到的频域软输入软输出线性MMSE均衡器框图。它包括FFT和IFFT变换装置、频域单点软输入软输出检测装置和软解调装置。
具体实施方式
下面结合图1、图2和图3对本发明的各个组成部分做进一步的详细说明。
假设在发送端,信息比特经过纠错编码、交织和调制后,通过频率选择性衰落信道。接收端接收信号可以表示为:
其中,hl是信道第1径的衰落系数,L为信道的阶数,Sk-1为第k-l时刻的发送符号,nk是均值为0方差为σ2的加性高斯白噪声。
假设接收端采用块处理,每次待均衡的发送信号的块长度为K,表示为s,相应的每次处理的接收信号样值个数为K+L,表示为y。例如,假设发送序列总共有N块,总长度为N×K,接收信号总长度为N×K+L,第一次块处理取出y1,y2,...,yK+L,第二次块处理取出yK,yK+1,...,y2K+L。这样,第n次处理时,接收信号y可以写为如下矩阵—矢量形式:
y=Hs+H′s′+H″s″+n [公式2]
其中:n是长度为K+L的噪声矢量,s′表示前一个块,s″表示下一个块,而H′s′表示前一个块对当前块的干扰,H″s″表示下一个块队当前块的干扰。上式中的信道矩阵如下所示:
我们假设前一帧的干扰H′s′可以被完全消除(例如前面符号为导频符号,或前一个块已经被均衡)。那么接收信号可以表示成:
y=Hs+H″s″+n [公式3]
将接收矢量y的前L个元素叠加在v的最后L个元素上。此时,[公式3]可以写为:
r1=Hcs+z [公式4]
其中,Hc为K×K的循环阵,它的的第一列为
它的非归一化的FFT为
λ=[λ1 λ2 … λK]T
噪声z均值为E(z)=0,协方差矩阵有如下形式:
定义:
z的协方差矩阵可以重新写为:σ2IK+R。
根据文献[1,S.Kay,Fundamental of statistical signal processing:estimationtheory,Prentice Hall,1993.],[公式4]的线性MMSE均衡可以表示为
上式可以进一步表示为:
定义
由于Hc是循环矩阵,那么∑也是循环矩阵,循环矩阵和矢量相乘可以用两次FFT和一次IFFT来实现,那么[公式6]中一个关键步骤就是求(IK+R∑)-1。根据矩阵求逆的多项式逼近,我们采用下式来近似它,
(IK+R∑)-1≈(2-γ)IK-γR∑
根据文献[2,G.M.A.Sessler and F.K.Jondral,“Rapidly converging polynomialexpansion multiuser detector with low complexity for CDMA systems,”ElectronicsLetters,vol.38,pp.997-998,2002.],γ可以表示为,
‖·‖∞表示矩阵的无穷范数,即矩阵行的和的最大值。
上式可以通过两次循环卷积来实现,也可以用两次FFT,两次IFFT求出。
那么根据[公式7]可以得到有偏的幅度为,
根据[公式8]可以进行软解调,从而得到发送比特的对数似然比。
把解调得到的比特似然比进行解交织,传递给软输入软输出译码器,译码输出的软信息经过交织后进行下一次的软输入软输出检测。根据该软信息,可以计算出发送符号的统计量,包括均值和方差。由于经过译码后,我们可以较准确地得到前一个块的干扰s′,以及下一个块的干扰s″,因此,我们可以将其抵消掉。这里我们假设理想抵消,这样在第二次迭代时,接收信号可以重新表示为,
y=Hs+n
将y的前L个样值叠加到后L个样值,可得,
虽然的协方差矩阵是对角阵,但是对角线的前K-L个元素为σ2最后L个元素2σ2。为降低复杂度,我们假设的协方差矩阵是对角线元素相等的对角阵,且对角线元素等于σ2。这样,在已知发送符号s的一阶和二阶统计量的情况下,我们可以针对[公式11]进行软输入软输出检测。
发送信号的软输入软输出线性MMSE可以表示为,
上式中,s是由根据解码反馈的软信息得到的发送信号的均值,v是发送信号方差的平均值,而
由于Hc是循环阵,[公式12]可以用两次FFT,一次IFFT求出。
检测器输出仍然做高斯近似,如[公式8],有偏幅度为ρ,等效噪声方差为ρ(1-vρ)。这样可以计算出比特软信息供译码器译码,从而完成第二次迭代。以后的迭代可以重复第二次迭代步骤。
Claims (4)
1.基于快速傅立叶变换的频域迭代均衡方法,其特征在于:初次均衡采用多项式展开的频域软输出均衡,随后的迭代中,抵消块间干扰,采用快速傅立叶变换实现的频域软输入软输出线性最小均方误差均衡;在初次均衡中,对于接收序列,取出K+L个样值,先抵消前一个块的干扰,然后将前L个样值叠加到最后的L个样值上,其中K表示每次均衡的发送符号块的长度,L为信道阶数;然后,采用多项式展开近似线性最小均方误差均衡中的矩阵求逆,并且利用信道的循环卷积特性,采用快速傅立叶变换实现;之后进行逐符号的软解调,得到发送比特的对数似然比,并输入到软输入软输出译码器译码,根据译码输出的软信息计算发送符号的统计量,包括均值和方差,进行第二次迭代均衡;根据反馈的信息,对接收序列的K+L个样值,抵消前一个块和后一个块对当前块的干扰,然后将前L个样值叠加到最后的L个样值上,接着进行频域线性MMSE软输入软输出均衡,输出软信息提供给译码器译码。
2.根据权利要求1所述的快速傅立叶变换的频域迭代均衡方法,其特征在于该迭代均衡方法包括以下几个步骤:
初次均衡:
1.)取出K+L个样值,表示为y,先抵消前一个块的干扰,然后将前L个样值叠加到最后的L个样值上,表示为r1;
2.)对信道的时域响应进行K点FFT变换,表示为λ=[λ1 λ2…λK]T,并保存下来,其中T表示矢量的转置;
3.)计算下一个块的干扰加上叠加引入的噪声的相关矩阵R;
4.)采用频域实现的多项式展开均衡,并对均衡输出做高斯近似,计算等效幅度以及等效噪声方差;
5.)进行软解调,解调输出经解交织后传递给软输入软输出译码器;
7.)对K+L个样值y,根据译码器获得的发送符号的信息,抵消前一个块和后一个块对当前块的干扰,然后将前L个样值叠加到最后的L个样值上,表示为r2;
8.)采用频域实现的线性MMSE软输入软输出均衡,并对均衡输出做高斯近似,计算等效幅度以及等效噪声方差;
9.)进行软解调,解调输出经解交织后传递给软输入软输出译码器;如果没有达到给定的迭代次数,则译码器输出的软信息经过交织后,计算发送信号的统计量,包括均值和方差,并对方差求平均,并跳转到步骤7)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101326910A CN101136896B (zh) | 2007-09-18 | 2007-09-18 | 基于快速傅立叶变换的频域迭代均衡方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101326910A CN101136896B (zh) | 2007-09-18 | 2007-09-18 | 基于快速傅立叶变换的频域迭代均衡方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101136896A true CN101136896A (zh) | 2008-03-05 |
CN101136896B CN101136896B (zh) | 2011-06-29 |
Family
ID=39160730
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007101326910A Expired - Fee Related CN101136896B (zh) | 2007-09-18 | 2007-09-18 | 基于快速傅立叶变换的频域迭代均衡方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101136896B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010015104A1 (en) * | 2008-08-04 | 2010-02-11 | Nxp B.V. | Iterative channel estimation method and apparatus for ici cancellation in multi-carrier systems |
CN102163976A (zh) * | 2011-01-20 | 2011-08-24 | 东南大学 | 一种基于gdft-ii变换的快速解码方法 |
CN105306396A (zh) * | 2015-11-03 | 2016-02-03 | 西安烽火电子科技有限责任公司 | 一种无线宽带通信信道迭代均衡的优化方法 |
CN106302282A (zh) * | 2016-08-31 | 2017-01-04 | 高拓讯达(北京)科技有限公司 | 判决反馈均衡处理方法及均衡器 |
WO2017050089A1 (zh) * | 2015-09-25 | 2017-03-30 | 电信科学技术研究院 | 参数处理方法、音频信号播放方法、装置及音频均衡器 |
CN110471018A (zh) * | 2019-09-19 | 2019-11-19 | 中南大学 | 一种频谱校正方法 |
CN111901262A (zh) * | 2020-07-30 | 2020-11-06 | 西安烽火电子科技有限责任公司 | 一种适用于短波通信的高阶调制Turbo时域均衡算法 |
CN112988732A (zh) * | 2021-04-14 | 2021-06-18 | 湖南工程学院 | 一种观测数据中异常值的处理方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100461116C (zh) * | 2003-04-02 | 2009-02-11 | 高通股份有限公司 | 用于在块相干通信系统中的交织的方法和设备 |
ATE503328T1 (de) * | 2004-08-13 | 2011-04-15 | Agency Science Tech & Res | Verfahren zum entzerren eines digitalsignals und entzerrer |
US8102907B2 (en) * | 2005-04-28 | 2012-01-24 | Qualcomm Incorporated | Space-frequency equalization for oversampled received signals |
-
2007
- 2007-09-18 CN CN2007101326910A patent/CN101136896B/zh not_active Expired - Fee Related
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010015104A1 (en) * | 2008-08-04 | 2010-02-11 | Nxp B.V. | Iterative channel estimation method and apparatus for ici cancellation in multi-carrier systems |
CN102113286A (zh) * | 2008-08-04 | 2011-06-29 | Nxp股份有限公司 | 用于多载波系统中的ici消除的迭代信道估计方法和设备 |
US8374266B2 (en) | 2008-08-04 | 2013-02-12 | Nxp B.V. | Iterative channel estimation method and apparatus for ICI cancellation in multi-carrier |
CN102113286B (zh) * | 2008-08-04 | 2014-06-18 | Nxp股份有限公司 | 用于多载波系统中的ici消除的迭代信道估计方法和设备 |
CN102163976A (zh) * | 2011-01-20 | 2011-08-24 | 东南大学 | 一种基于gdft-ii变换的快速解码方法 |
CN102163976B (zh) * | 2011-01-20 | 2013-02-27 | 东南大学 | 一种基于gdft-ii变换的快速解码方法 |
WO2017050089A1 (zh) * | 2015-09-25 | 2017-03-30 | 电信科学技术研究院 | 参数处理方法、音频信号播放方法、装置及音频均衡器 |
US10326422B2 (en) | 2015-09-25 | 2019-06-18 | China Academy Of Telecommunications Technology | Parameter processing method, audio signal playing method and device, and audio equalizer |
CN105306396A (zh) * | 2015-11-03 | 2016-02-03 | 西安烽火电子科技有限责任公司 | 一种无线宽带通信信道迭代均衡的优化方法 |
CN105306396B (zh) * | 2015-11-03 | 2018-11-09 | 西安烽火电子科技有限责任公司 | 一种无线宽带通信信道迭代均衡的优化方法 |
CN106302282A (zh) * | 2016-08-31 | 2017-01-04 | 高拓讯达(北京)科技有限公司 | 判决反馈均衡处理方法及均衡器 |
CN106302282B (zh) * | 2016-08-31 | 2019-09-13 | 高拓讯达(北京)科技有限公司 | 判决反馈均衡处理方法及均衡器 |
CN110471018A (zh) * | 2019-09-19 | 2019-11-19 | 中南大学 | 一种频谱校正方法 |
CN111901262A (zh) * | 2020-07-30 | 2020-11-06 | 西安烽火电子科技有限责任公司 | 一种适用于短波通信的高阶调制Turbo时域均衡算法 |
CN112988732A (zh) * | 2021-04-14 | 2021-06-18 | 湖南工程学院 | 一种观测数据中异常值的处理方法 |
CN112988732B (zh) * | 2021-04-14 | 2023-10-20 | 湖南工程学院 | 一种观测数据中异常值的处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101136896B (zh) | 2011-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101136896B (zh) | 基于快速傅立叶变换的频域迭代均衡方法 | |
US8325588B2 (en) | Iterative reception method and iterative receiver | |
CN102624652B (zh) | Ldpc解码方法和装置及接收终端 | |
CN101621327A (zh) | 单载波发送系统中的无线电通信方法和装置 | |
CN1989750A (zh) | 用于mimo cdma下行链路的基于傅立叶变换的线形均衡 | |
JP4827695B2 (ja) | 無線受信装置 | |
WO2007139145A1 (ja) | 無線受信機、無線通信システムおよび無線受信方法 | |
EP3229429B1 (en) | Methods and devices for symbols detection in multi antenna systems | |
KR20080059014A (ko) | 다중 입출력 무선통신 시스템에서 수신 장치 및 방법 | |
CN106982106B (zh) | 递归子块解码 | |
US8139669B2 (en) | Space domain filter detecting method in a multi-antenna wireless communication system | |
CN101005299B (zh) | 信号处理的方法和系统 | |
US8107546B2 (en) | Detection method of space domain maximum posteriori probability in a wireless communication system | |
CN102868434B (zh) | 一种mimo检测方法及装置 | |
CN104184505A (zh) | 发射信号的多输入多输出mimo检测方法、装置及系统 | |
US7974334B2 (en) | Systems and methods for hybrid-MIMO equalization | |
US20190253090A1 (en) | Stochastic interference cancellation | |
CN101807939A (zh) | 多用户检测方法和系统 | |
Kim et al. | MMSE-based lattice-reduction-aided fixed-complexity sphere decoder for low-complexity near-ML MIMO detection | |
US8018985B2 (en) | CDMA receivers for the golden code | |
CN101005298A (zh) | 信号处理的方法和系统 | |
Chang et al. | Cancellation of ICI by Doppler effect in OFDM systems | |
JP2020141173A (ja) | 無線通信システム、無線通信方法、送信局装置および受信局装置 | |
CN102664852B (zh) | 多输入多输出正交频分复用系统中的软输入软输出检测方法 | |
JP2011139294A (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: 20110629 Termination date: 20120918 |