基于快速傅立叶变换的频域迭代均衡方法
专利汇可以提供基于快速傅立叶变换的频域迭代均衡方法专利检索，专利查询，专利分析的服务。并且基于快速 傅立叶变换 的频域 迭代 均衡方法可以有效实现 频率 选择性信道下的迭代均衡。该方法为：首先，对于接收序列，取出K+L个样值，将前L个样值 叠加 到最后的L个样值上，其中K表示每次均衡的发送符号 块 的长度，L为信道阶数；然后，初次均衡采用多项式展开的频域 软输出 均衡；解交织后，根据均衡输出的软信息进行 软输入 软输出解码；解码输出的软信息进行交织后，计算符号的统计量；利用反馈的符号均值进行块间干扰抵消，然后采用频域软输入软输出均衡进入下一轮迭代。该方法可以在没有块间保护的系统进行频域迭代均衡，可以采用 快速傅立叶变换 来实现，在提高系统性能的同时，可以大大降低迭代接收机的复杂度。，下面是基于快速傅立叶变换的频域迭代均衡方法专利的具体信息内容。
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.)进行软解调，解调输出经解交织后传递给软输入软输出译码器；
6.)译码器输出的软信息经过交织后，计算发送信号的统计量，包括均值 和方差，并对方差求平均，表示为v； 迭代均衡：
7.)对K+L个样值y，根据译码器获得的发送符号的信息，抵消前一个块 和后一个块对当前块的干扰，然后将前L个样值叠加到最后的L个样值上，表 示为r2；
8.)采用频域实现的线性MMSE软输入软输出均衡，并对均衡输出做高斯 近似，计算等效幅度以及等效噪声方差；
9.)进行软解调，解调输出经解交织后传递给软输入软输出译码器；如果 没有达到给定的迭代次数，则译码器输出的软信息经过交织后，计算发送信号 的统计量，包括均值和方差，并对方差求平均，并跳转到步骤7)。
3.根据权利要求2所述的快速傅立叶变换的频域迭代均衡方法，其特征 在于第4)步所述的多项式展开软输出频域均衡，它包括以下几个步骤：
41.)计算频域均衡所用到的循环矩阵∑，该循环矩阵的第一列为：
的K点非归一化的IFFT变换，其中σ2为噪声方差；
42.)计算多项式展开时的归一化因子
43.)计算：(2-γ)r1-γR∑r1；
44.)将步骤43)的结果进行K点FFT变换，并采用如下系数的单点均衡： 其中k＝1，2，...，K；
45.)将步骤44)的结果进行K点IFFT变换，并计算均衡器输出的有偏系 数以及等效噪声方差，最后进行软解调。
4.根据权利要求2所述的快速傅立叶变换的频域迭代均衡方法，其特征 在于第8)步所述的频域实现的线性MMSE软输入软输出均衡，它包括以下几 个步骤：
81.)对接收矢量r2进行K点FFT变换；
82.)对反馈的均值信号进行K点FFT变换并对其输出的每一个点乘以λk
其中k＝1，2，...，K；
83.)将步骤81)的结果减去步骤82)的结果，并对结果进行使用系数的 单点均衡：其中k＝1，2，...，K；
84.)均衡之后的结果进行K点IFFT变换，并加上，其中
85.)计算均衡器输出的有偏系数以及等效噪声方差，最后进行软解调。
技术领域
本发明是一种应用于无线通信系统的联合均衡和译码的方法，属于移动通 信中的均衡技术领域。
背景技术
发明内容
技术方案：本发明的基于快速傅立叶变换的频域迭代均衡方法的初次均衡 采用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.4.)将步骤4.3)的结果进行K点FFT变换，并采用如下系数的单点均 衡：其中k＝1，2,...，K；
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]
其中，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：estimation theory，Prentice Hall，1993.]，[公式4]的线性MMSE均衡可以表示为
[公式5]
上式可以进一步表示为：
[公式6]
定义
由于Hc是循环矩阵，那么∑也是循环矩阵，循环矩阵和矢量相乘可以用两次 FFT和一次IFFT来实现，那么[公式6]中一个关键步骤就是求(IK+R∑)-1。根据 矩阵求逆的多项式逼近，我们采用下式来近似它，
(IK+R∑)-1≈(2-γ)IK-γR∑
根据文献[2，G.M.A.Sessler and F.K.Jondral，“Rapidly converging polynomial expansion multiuser detector with low complexity for CDMA systems,”Electronics Letters，vol.38，pp.997-998，2002.]，γ可以表示为，
‖·‖∞表示矩阵的无穷范数，即矩阵行的和的最大值。
[公式7]
上式可以通过两次循环卷积来实现，也可以用两次FFT，两次IFFT求出。
上面均衡器是有偏估计，假设均衡器的第k个输出值可以表示为有偏的幅 度ρk以信号sk再加上一个均值为0方差为的高斯噪声：
[公式8]
那么根据[公式7]可以得到有偏的幅度为，
[公式9]
噪声方差为，
[公式10]
根据[公式8]可以进行软解调，从而得到发送比特的对数似然比。
把解调得到的比特似然比进行解交织，传递给软输入软输出译码器，译码 输出的软信息经过交织后进行下一次的软输入软输出检测。根据该软信息，可 以计算出发送符号的统计量，包括均值和方差。由于经过译码后，我们可以较 准确地得到前一个块的干扰s′，以及下一个块的干扰s″，因此，我们可以将其 抵消掉。这里我们假设理想抵消，这样在第二次迭代时，接收信号可以重新表 示为，
y＝Hs+n
将y的前L个样值叠加到后L个样值，可得，
[公式11]
虽然的协方差矩阵是对角阵，但是对角线的前K-L个元素为σ2最后L个元素 2σ2。为降低复杂度，我们假设的协方差矩阵是对角线元素相等的对角阵，且 对角线元素等于σ2。这样，在已知发送符号s的一阶和二阶统计量的情况下， 我们可以针对[公式11]进行软输入软输出检测。
发送信号的软输入软输出线性MMSE可以表示为，
[公式12]
上式中，s是由根据解码反馈的软信息得到的发送信号的均值，v是发送信号 方差的平均值，而
由于Hc是循环阵，[公式12]可以用两次FFT，一次IFFT求出。
检测器输出仍然做高斯近似，如[公式8]，有偏幅度为ρ，等效噪声方差为 ρ(1-vρ)。这样可以计算出比特软信息供译码器译码，从而完成第二次迭代。 以后的迭代可以重复第二次迭代步骤。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 可变长快速傅立叶变换电路及实现方法 | 2020-05-13 | 204 |
| 快速傅立叶变换电路 | 2020-05-11 | 679 |
| 快速傅立叶变换/快速傅立叶反变换处理器及其处理方法 | 2020-05-11 | 219 |
| 基于Turbo结构的快速傅立叶DFT变换系统 | 2020-05-12 | 940 |
| 一种并行快速傅立叶变换处理方法 | 2020-05-12 | 314 |
| 快速傅立叶变换运算的设备 | 2020-05-12 | 842 |
| 可变长度的快速傅立叶变换装置 | 2020-05-11 | 885 |
| 可变长度的快速傅立叶变换装置 | 2020-05-11 | 834 |
| 补偿输出位信号的快速傅立叶变换设备及方法 | 2020-05-13 | 952 |
| 一种基于快速傅立叶变换的打盹行为检测方法 | 2020-05-13 | 38 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
