CN105471774A - 一种相位噪声的估计方法、估计装置、接收机及通信设备 - Google Patents

Abstract

本发明提供一种相位噪声的估计方法、估计装置、接收机及通信设备，属于无线通信领域。所述相位噪声的估计方法包括如下步骤：获取接收数据的均衡后符号r(k)对应的发送符号其中，k＝0,...,K-1，K为一帧接收数据包括的符号数目；根据r(k)和获取公共相位噪声；根据r(k)、和所述公共相位噪声获取残余相位噪声；根据所述公共相位噪声、所述残余相位噪声和离散余弦变换的基函数矩阵获取相位噪声。本发明能够提高相位噪声估计的准确度，从而提高接收机的接收性能。

Description

一种相位噪声的估计方法、估计装置、接收机及通信设备

技术领域

本发明涉及无线通信领域，特别涉及一种相位噪声的估计方法、估计装置、接收机及通信设备。

背景技术

在无线通信中，常常由于发射机和接收机本振的不稳定性，使得载波信号不是单一频率的正弦波，这种不稳定性，通常用相位噪声来衡量。由于这个原因，在无线通信系统的载波频率比较高时，以及，在无线通信系统使用高阶的调制方式时，相位噪声的影响通常是不可忽略的。因此，接收机必须对整个链路引入的相位噪声进行估计，对受相位噪声影响的接收信号进行相噪补偿，以达到更加准确的检测发射数据符号的目的。

目前进行相位噪声估计的方法，主要有以下几种：

方法1、利用锁相环对相位噪声进行跟踪，目前一般采用二阶锁相环。

方法2、将接收信号的时间分为多个时间片，将每个时间片内的相位噪声近似为常数，通过适当的统计算法来获得在特定时间片的相位平均值。

方法3、基于插值的算法：构造特殊的帧结构，比如，每一帧的首尾都有导频，先对每一帧数据的导频位置的相位噪声进行估计，帧内数据符号部分的相位噪声通过插值来实现。

方法4、基于信道和相位噪声的联合检测估计算法，此方法的复杂度非常高，信号处理的时延很大。

上述的方法1，锁相环需要一定的捕获时间才能锁定相位，造成相位噪声的补偿不及时。

上述的方法2～4，需要知道相位噪声的统计特性，相位噪声估计的准确度不高。

发明内容

本发明所要解决的技术问题是提供一种相位噪声的估计方法、估计装置、接收机及通信设备，以提高相位噪声估计的准确度，从而提高接收机的接收性能。

为解决上述技术问题，本发明提供技术方案如下：

一种相位噪声的估计方法，包括如下步骤：

获取接收数据的均衡后符号r(k)对应的发送符号其中，k＝0,...,K-1，K为一帧接收数据包括的符号数目；

根据r(k)和获取公共相位噪声；

根据r(k)、和所述公共相位噪声获取残余相位噪声；以及

根据所述公共相位噪声、所述残余相位噪声和离散余弦变换的基函数矩阵获取相位噪声。

优选地，所述获取接收数据的均衡后符号r(k)对应的发送符号包括：

对r(k)进行硬判决，得到或者

对r(k)进行软判决，并对软判决结果进行重构得到

优选地，按照如下公式获取公共相位噪声：

${\widetilde{\mathrm{\θ}}}_{avg}=\arg \left(\underset{k=0}{\overset{K-1}{\Σ}}r\right(k\left){\tilde{a}}^{*}\right(k\left)\right)$

其中，为公共相位噪声，()^*表示求共轭运算，arg()表示求复数的辐角运算。

优选地，按照如下公式获取残余相位噪声：

$\widetilde{\mathrm{\φ}}\left(k\right)=\arg \left(r\right(k\left){\tilde{a}}^{*}\right(k\left)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}\right)$

其中，为第k个符号的残余相位噪声。

优选地，按照如下公式获取相位噪声：

$\widetilde{\mathrm{\θ}}\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{k^{\′}=0}{\overset{K-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\mathrm{\Ψ}}_K^T\right)}_{k,k^{\′}}\widetilde{\mathrm{\φ}}\left(k^{\′}\right)$

其中，为相位噪声，Ψ_K为离散余弦变换的基函数矩阵，为Ψ_K的转置，(Ψ_K)_k,n＝ψ_n(k)，ψ_n(_k)为离散余弦变换的第n个基函数，n＝0,...,N-1，N为所述基函数的数目。

优选地，按照如下公式获取残余相位噪声：

${\widetilde{\mathrm{\φ}}}_B\left(m\right)=\frac{1}{B}\arg \left(\underset{b=0}{\overset{B-1}{\Σ}}r\left(b+mB\right)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}{\tilde{a}}^{*}\left(b+mB\right)\right)$

其中，一帧接收数据的K个符号被分成M段，每段包括B个符号，为第m段符号的残余相位噪声，m＝0,...,M-1。

优选地，按照如下公式获取相位噪声：

${\widetilde{\mathrm{\θ}}}_B\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{m=0}{\overset{M-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\left({\mathrm{\Ψ}}_{avg}^T{\mathrm{\Ψ}}_{avg}\right)}^{-1}{\mathrm{\Ψ}}_{avg}\right)}_{k,m}{\widetilde{\mathrm{\φ}}}_B\left(m\right)$

其中，为相位噪声，Ψ_K为离散余弦变换的基函数矩阵，(Ψ_K)_k,n＝ψ_n(k)，ψ_n(k)为离散余弦变换的第n个基函数，n＝0,...,N-1，N为所述基函数的数目，为Ψ_avg的转置， ${\left({\mathrm{\Ψ}}_{avg}\right)}_{m,n}=\frac{1}{B}\underset{b=0}{\overset{B-1}{\Σ}}{\mathrm{\ψ}}_n(b+mB).$

一种相位噪声的估计装置，包括：

发送符号估计单元，用于获取接收数据的均衡后符号r(k)对应的发送符号其中，k＝0,...,K-1，K为一帧接收数据包括的符号数目；

公共相位噪声估计单元，用于根据r(k)和获取公共相位噪声；

残余相位噪声估计单元，用于根据r(k)、和所述公共相位噪声获取残余相位噪声；以及

相位噪声估计单元，用于根据所述公共相位噪声、所述残余相位噪声和离散余弦变换的基函数矩阵获取相位噪声。

优选地，所述发送符号估计单元进一步用于：

对r(k)进行硬判决，得到或者

对r(k)进行软判决，并对软判决结果进行重构得到

优选地，所述公共相位噪声估计单元进一步用于，按照如下公式获取公共相位噪声：

${\widetilde{\mathrm{\θ}}}_{avg}=\arg \left(\underset{k=0}{\overset{K-1}{\Σ}}r\right(k\left){\tilde{a}}^{*}\right(k\left)\right)$

其中，为公共相位噪声，()^*表示求共轭运算，arg()表示求复数的辐角运算。

优选地，所述残余相位噪声估计单元进一步用于，按照如下公式获取残余相位噪声：

$\widetilde{\mathrm{\φ}}\left(k\right)=\arg \left(r\right(k\left){\tilde{a}}^{*}\right(k\left)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}\right)$

其中，为第k个符号的残余相位噪声。

优选地，所述相位噪声估计单元进一步用于，按照如下公式获取相位噪声：

$\widetilde{\mathrm{\θ}}\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{k^{\′}=0}{\overset{K-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\mathrm{\Ψ}}_K^T\right)}_{k,k^{\′}}\widetilde{\mathrm{\φ}}\left(k^{\′}\right)$

其中，为相位噪声，Ψ_K为离散余弦变换的基函数矩阵，为Ψ_K的转置，(Ψ_K)_k,n＝ψ_n(k)，ψ_n(k)为离散余弦变换的第n个基函数，n＝0,...,N-1，N为所述基函数的数目。

优选地，所述残余相位噪声估计单元进一步用于，按照如下公式获取残余相位噪声：：

${\widetilde{\mathrm{\φ}}}_B\left(m\right)=\frac{1}{B}\arg \left(\underset{b=0}{\overset{B-1}{\Σ}}r\left(b+mB\right)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}{\tilde{a}}^{*}\left(b+mB\right)\right)$

其中，一帧接收数据的K个符号被分成M段，每段包括B个符号，为第m段符号的残余相位噪声，m＝0,...,M-1。

优选地，所述相位噪声估计单元进一步用于，按照如下公式获取相位噪声：

${\widetilde{\mathrm{\θ}}}_B\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{m=0}{\overset{M-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\left({\mathrm{\Ψ}}_{avg}^T{\mathrm{\Ψ}}_{avg}\right)}^{-1}{\mathrm{\Ψ}}_{avg}\right)}_{k,m}{\widetilde{\mathrm{\φ}}}_B\left(m\right)$

其中，为相位噪声，Ψ_K为离散余弦变换的基函数矩阵，(Ψ_K)_k,n＝ψ_n(k)，ψ_n(k)为离散余弦变换的第n个基函数，n＝0,...,N-1，N为所述基函数的数目，为Ψ_avg的转置， ${\left({\mathrm{\Ψ}}_{avg}\right)}_{m,n}=\frac{1}{B}\underset{b=0}{\overset{B-1}{\Σ}}{\mathrm{\ψ}}_n(b+mB).$

一种包括上述的估计装置的接收机。

一种包括上述的接收机的通信设备。

与现有技术相比，本发明的有益效果是：

本发明的技术方案不需要捕获时间，相噪补偿可以快速启动，也不需要知道相位噪声的统计特性，相位噪声估计的准确度较高。利用本发明的技术方案，能够有效的对相位噪声进行估计，从而对受到相位噪声影响的接收信号进行补偿，达到准确检测发射符号的目的，进而能够提高系统的性能。

附图说明

图1为根据本发明实施例的相位噪声的估计方法的流程图；

图2为根据本发明实施例的相位噪声的估计装置的结构图；

图3为根据本发明实施例的接收机的结构图；

图4为根据本发明另一实施例的接收机的结构图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图及具体实施例对本发明进行详细描述。

首先对本发明实施例的相位噪声估计的原理进行介绍。

由于带宽远低于采样率的离散随机过程，可以用基函数进行展开，而且展开项的数目通常不会很多，一般无线通信系统的相位噪声在数字基带可以视为一个这样的离散随机过程。因此，本发明实施例将相位噪声过程，用离散余弦变换(DCT)函数进行正交展开。

不失一般性，设一帧基带数据包括K个符号。令k时刻(对应第k个符号)，接收端设备的接收机中经过均衡器输出的均衡后符号为r(k)，r(k)对应的发送端设备发送符号为a(k)，遭受到的高斯白噪声为n(k)，相位噪声为θ(k)，则均衡后符号r(k)可以表示为：

r(k)＝a(k)e^jθ(k)+n(k)，k＝0，...，K-1(1)

本发明的目的在于通过一种合理的算法，来对相位噪声进行估计得到第k个符号的相位噪声估计值在获得后，对r(k)进行补偿，得到进行相位噪声补偿以后的信号z(k)，如下：

$z\left(k\right)=r\left(k\right)e^{-j\widetilde{\mathrm{\θ}}\left(k\right)}---\left(2\right)$

后续的均衡解调等在z(k)上进行，相位噪声被消除了。

下面给出求的算法。

定义公共相位噪声

${\widetilde{\mathrm{\θ}}}_{avg}=\arg \left(\underset{k=0}{\overset{K-1}{\Σ}}r\left(k\right){\tilde{a}}^{*}\left(k\right)\right)---\left(3\right)$

其中，()^*表示求共轭运算，arg()表示求复数的辐角运算，表示对发送符号a(k)的估计。对于导频部分，可以直接采用相应的导频符号，对于数据部分，需要进行均衡判决得到。此判决可以是硬判决，见图3所示，也可以是软判决加重构得到，见图4所示。其中软判决相对性能要更好一些，因为软判决是有编码增益的，但是基于软判决的方法复杂度高，时延也相对大一些。另外，公共相位噪声表征的是一帧数据所包括的所有符号对应的相位噪声的统计平均值，为了获取各符号分别对应的相位噪声，还需要获得各自的残余相位噪声φ(k)。

第k个符号对应的残余相位噪声φ(k)为一个低通随机过程，表示如下：

$\mathrm{\φ}\left(k\right)=\mathrm{\θ}\left(k\right)-{\widetilde{\mathrm{\θ}}}_{avg}---\left(4\right)$

φ(k)可以用DCT函数基进行展开如下：

$\mathrm{\φ}\left(k\right)\≈\underset{n=0}{\overset{N-1}{\Σ}}{x}_n{\mathrm{\ψ}}_n\left(k\right)---\left(5\right)$

其中，ψ_n(k)为DCT基函数，N为所述基函数的数目，x_n为第n个基函数的系数，ψ_n(k)的取值如下：

${\mathrm{\ψ}}_n\left(k\right)=\left\{\begin{array}{c}\sqrt{\frac{1}{K}},n=0\\ \sqrt{\frac{2}{K}}cos\left(\frac{\mathrm{\π}n}{K}\right(k+\frac{1}{2}\left)\right),n>0\end{array}\right.---\left(6\right)$

实际上，可以按照如下公式得到φ(k)的最大似然估计

$\begin{array}{c}\widetilde{\mathrm{\φ}}\left(k\right)=\arg \left(r\right(k\left){\tilde{a}}^{*}\right(k\left)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}\right)\\ \≈\mathrm{\θ}\left(k\right)-{\widetilde{\mathrm{\θ}}}_{avg}+n\left(k\right)\end{array}---\left(7\right)$

另外，还可以表示为：

$\widetilde{\mathrm{\φ}}\left(k\right)=\underset{n=0}{\overset{N-1}{\Σ}}{\tilde{x}}_n{\mathrm{\ψ}}_n\left(k\right)---\left(8\right)$

其中，为对进行DFT展开时的第n个基函数的系数。

构造一个K×N阶基函数矩阵Ψ_K，Ψ_K的第k行第n列的元素为：

(Ψ_K)_k,n＝ψ_n(k),k＝0,...,K-1；n＝0,...,N-1(9)

由公式(4)得到：

$\widetilde{\mathrm{\θ}}\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\widetilde{\mathrm{\φ}}\left(k\right)---\left(10\right)$

将公式(8)代入公式(10)得到相位噪声的估计值如下：

$\widetilde{\mathrm{\θ}}\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{k^{\′}=0}{\overset{K-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\mathrm{\Ψ}}_K^T\right)}_{k,k^{\′}}\widetilde{\mathrm{\φ}}\left(k^{\′}\right)---\left(11\right)$

其中，为Ψ_K的转置。

需要说明的是，公式(7)在信噪比相对比较高的时候是近似成立的，在信噪比比较低时不够准确。实际上，在信噪比比较低时，可以通过分段来降低噪声的影响。具体可以这样来处理，将每一帧数据的K个符号分成M段，每段包括B个符号，即满足K＝BM。

然后，令

$\begin{array}{c}{\widetilde{\mathrm{\φ}}}_B\left(m\right)=\frac{1}{B}\arg \left(\underset{b=0}{\overset{B-1}{\Σ}}r\left(b+mB\right)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}{\tilde{a}}^{*}\left(b+mB\right)\right)\\ \≈\frac{1}{B}\underset{b=0}{\overset{B-1}{\Σ}}\mathrm{\φ}(b+mB)+u\left(m\right)\\ \≈\underset{n=0}{\overset{N-1}{\Σ}}{x}_n{\mathrm{\ψ}}_n^{avg}\left(m\right)+u\left(m\right)\end{array}$

${\mathrm{\ψ}}_n^{avg}\left(m\right)=\frac{1}{B}\underset{b=0}{\overset{B-1}{\Σ}}{\mathrm{\ψ}}_n(b+mB)$

其中，表示m段符号的残余相位噪声的最大似然估计，u(m)为第m段符号遭受到的高斯白噪声，m＝0,...,M-1。

经过推导可以得到第k个符号的相位噪声的估计值为：

${\widetilde{\mathrm{\θ}}}_B\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{m=0}{\overset{M-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\left({\mathrm{\Ψ}}_{avg}^T{\mathrm{\Ψ}}_{avg}\right)}^{-1}{\mathrm{\Ψ}}_{avg}\right)}_{k,m}{\widetilde{\mathrm{\φ}}}_B\left(m\right)$

其中，()^T为求矩阵的转置运算。

需要说明的是，以上公式推导考虑的是高斯白噪声(AWGN)信道环境下，加入相位噪声的处理方法。对于多径信道环境，可以通过合理的转化将其转化为AWGN信道环境下的模型来处理，具体如下：

对于多径信道环境的情况，信道的多径为{h_l},l＝0,...,L-1，L为符号数，可以做如下的转化，转化为AWGN模型。

y(k)＝v(k)e^jθ(k)+n(k)

v(k)＝a(k)*h(k)

其中，h(k)为信道冲激响应，v(k)为理想接收信号，y(k)为实际受相位噪声影响的接收信号，*表示卷积。

另外需要特别指出的是，接收机中均衡器输出的符号在做功率归一化处理之后，一般可以近似认为是AWGN信道。这里的均衡器可以是时域均衡器，也可以是频域均衡器，可以是利用信道估计的均衡器，也可以是不使用信道估计的自适应均衡器，比如最小均方(LMS)均衡器，递归式最小均方(RLS)均衡器。

基于以上相位噪声的估计原理，本发明实施例提供一种相位噪声的估计方法。

参照图1，本发明实施例的相位噪声的估计方法，可以包括如下步骤：

步骤101：获取接收数据的均衡后符号r(k)对应的发送符号其中，k＝0,...,K-1，K为一帧接收数据包括的符号数目；

步骤102：根据r(k)和获取公共相位噪声；

步骤103：根据r(k)、和所述公共相位噪声获取残余相位噪声；以及

步骤104：根据所述公共相位噪声、所述残余相位噪声和离散余弦变换的基函数矩阵获取相位噪声。

在步骤101中，接收机首先对接收数据进行均衡处理，得到均衡后符号r(k)，然后，可以对r(k)进行硬判决得到对应的发送符号估计也可以对r(k)进行软判决，并对软判决结果进行重构得到对应的发送符号估计

可选地，对于多径信道环境，可以对均衡后符号r(k)进行功率归一化处理，以进一步提高相位噪声估计的准确性。

在得到一帧数据的K个符号分别对应的发送符号估计后，在步骤102中，可以按照如下公式获取这一帧数据对应的公共相位噪声：

${\widetilde{\mathrm{\θ}}}_{avg}=\arg \left(\underset{k=0}{\overset{K-1}{\Σ}}r\right(k\left){\tilde{a}}^{*}\right(k\left)\right)$

其中，为公共相位噪声，()^*表示求共轭运算，arg()表示求复数的辐角运算。

在步骤103中，可以按照如下公式获取第k个符号对应的残余相位噪声估计：

$\widetilde{\mathrm{\φ}}\left(k\right)=\arg \left(r\right(k\left){\tilde{a}}^{*}\right(k\left)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}\right)$

在步骤104中，可以按照如下公式获取第k个符号对应的相位噪声估计：

$\widetilde{\mathrm{\θ}}\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{k^{\′}=0}{\overset{K-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\mathrm{\Ψ}}_K^T\right)}_{k,k^{\′}}\widetilde{\mathrm{\φ}}\left(k^{\′}\right)$

其中，Ψ_K为离散余弦变换的基函数矩阵，为Ψ_K的转置，(Ψ_K)_k,n＝ψ_n(k)，ψ_n(k)为离散余弦变换的第n个基函数，n＝0,...,N-1，N为所述基函数的数目。

另外，在步骤103中，还可以按照如下公式获取残余相位噪声估计：

${\widetilde{\mathrm{\φ}}}_B\left(m\right)=\frac{1}{B}\arg \left(\underset{b=0}{\overset{B-1}{\Σ}}r\left(b+mB\right)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}{\tilde{a}}^{*}\left(b+mB\right)\right)$

其中，一帧接收数据的K个符号被分成M段，每段包括B个符号，为第m段符号的残余相位噪声估计，m＝0,...,M-1。

相应地，在步骤104中，按照如下公式获取相位噪声估计：

${\widetilde{\mathrm{\θ}}}_B\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{m=0}{\overset{M-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\left({\mathrm{\Ψ}}_{avg}^T{\mathrm{\Ψ}}_{avg}\right)}^{-1}{\mathrm{\Ψ}}_{avg}\right)}_{k,m}{\widetilde{\mathrm{\φ}}}_B\left(m\right)$

其中，为第k个符号的相位噪声估计，Ψ_K为离散余弦变换的基函数矩阵，(Ψ_K)_k,n＝ψ_n(k)，ψ_n(k)为离散余弦变换的第n个基函数，n＝0,...,N-1，N为所述基函数的数目，为Ψ_avg的转置， ${\left({\mathrm{\Ψ}}_{avg}\right)}_{m,n}=\frac{1}{B}\underset{b=0}{\overset{B-1}{\Σ}}{\mathrm{\ψ}}_n(b+mB).$

上述步骤中，ψ_n(k)的取值如下：

${\mathrm{\ψ}}_n\left(k\right)=\left\{\begin{array}{c}\sqrt{\frac{1}{K}},n=0\\ \sqrt{\frac{2}{K}}cos\left(\frac{\mathrm{\π}n}{K}\left(k+\frac{1}{2}\right)\right),n>0\end{array}\right.$

对应于上述相位噪声的估计方法，本发明实施例还提供一种相位噪声的估计装置。

参照图2，本发明实施例的相位噪声的估计装置，可以包括：

发送符号估计单元10，用于获取接收数据的均衡后符号r(k)对应的发送符号其中，k＝0,...,K-1，K为一帧接收数据包括的符号数目；

公共相位噪声估计单元20，与发送符号估计单元10耦合，用于根据r(k)和获取公共相位噪声；

残余相位噪声估计单元30，与发送符号估计单元10和公共相位噪声估计单元20耦合，用于根据r(k)、和所述公共相位噪声获取残余相位噪声；以及

相位噪声估计单元40，与公共相位噪声估计单元20和残余相位噪声估计单元30耦合，用于根据所述公共相位噪声、所述残余相位噪声和离散余弦变换的基函数矩阵获取相位噪声。

该估计装置能够执行方法实施例中所描述的多个处理，其具体的工作过程以及工作原理在方法部分已经进行了详细描述，在此不再赘述，参照方法中相应部分的描述即可。

以下给出上述相位噪声估计方法及估计装置在无线通信设备接收机中的具体应用，所述无线通信设备可以是数字微波系统中的通信设备；也可以是其他移动通信系统中的通信设备，例如手机、PDA等。

参照图3，本发明实施例的接收机可以包括：相位噪声补偿模块31、FFT模块32、信道估计模块33、频域均衡器34、IFFT模块35、硬判决模块36和相位噪声估计模块37。在本实施例中，相位噪声估计基于硬判决。

该接收机的工作过程如下：

(1)首先，接收基带信号z(k)，第一次均衡处理时，相位噪信息未知，此时，z(k)＝y(k)，y(k)为未经相噪补偿的基带信号，然后，信道估计模块33进行信道估计，得到信道估计的频域表示形式H(k)，其中，k＝0,...,K-1，K为一帧接收数据包括的符号数目。

(2)FFT模块32对时域的基带信号z(k)作FFT变换(快速傅里叶变换)，得到基带信号z(k)的频域表示形式Z(k)。

(3)频域均衡器34使用迫零或者MMSE(最小均方误差)准则等对Z(k)进行频域均衡，得到频域的均衡输出符号R(k)。

(4)IFFT模块35对频域均衡器34输出的频域符号R(k)作IFFT变换(快速傅里叶逆变换)，得到均衡后符号的时域表示形式r(k)。

(5)硬判决模块36对r(k)进行硬判决得到时域符号即为发送符号的估计值。

(6)相位噪声估计模块37根据r(k)和进行相位噪声的估计，具体如下：

首先，利用公式(3)求得公共相位噪声

其次，利用公式(6)计算DCT基函数ψ_n(k)；

然后，利用公式(7)计算利用公式(9)计算Ψ_K；

最后，利用公式(11)计算相位噪声

(7)在获得相位噪声以后，由相位噪声补偿模块31对基带信号按照公式(2)进行相位噪声补偿，得到补偿后的基带信号z(k)。

(8)频域均衡器34对补偿后的基带信号z(k)进行均衡处理，得到新的均衡输出r(k)。

上述过程可以进行迭代多次，以更好的消除相位噪声。

参照图4，本发明另一实施例的接收机可以包括：相位噪声补偿模块41、FFT模块42、信道估计模块43、频域均衡器44、IFFT模块45、软判决模块46、相位噪声估计模块47、解交织模块48、信道译码模块49、信道编码模块50、交织模块51和调制模块52。在本实施例中，相位噪声估计基于软判决。

该接收机的工作过程如下：

(1)首先，接收基带信号z(k)，第一次均衡处理时，相位噪信息未知，此时，z(k)＝y(k)，y(k)为未经相噪补偿的基带信号，然后，信道估计模块43进行信道估计，得到信道估计的频域表示形式H(k)，其中，k＝0,...,K-1，K为一帧接收数据包括的符号数目。

(2)FFT模块42对时域的基带信号z(k)作FFT变换(快速傅里叶变换)，得到基带信号z(k)的频域表示形式Z(k)。

(3)频域均衡器44使用迫零或者MMSE(最小均方误差)准则等对Z(k)进行频域均衡，得到频域的均衡输出符号R(k)。

(4)IFFT模块45对频域均衡器44输出的频域符号R(k)作IFFT变换(快速傅里叶逆变换)，得到均衡后符号的时域表示形式r(k)。

(5)软判决模块46对r(k)进行软判决得到比特的软信息

(6)由解交织模块48、信道译码模块49、信道编码模块50、交织模块51和调制模块52对进行解交织、信道译码、信道编码、交织、调制，得到时域符号，即为发送符号的估计值。

(7)相位噪声估计模块47根据r(k)和进行相位噪声的估计，具体如下：

首先，利用公式(3)求得公共相位噪声

其次，利用公式(6)计算DCT基函数ψ_n(k)；

然后，利用公式(7)计算利用公式(9)计算Ψ_K；

最后，利用公式(11)计算相位噪声

(8)在获得相位噪声以后，由相位噪声补偿模块41对基带信号按照公式(2)进行相位噪声补偿，得到补偿后的基带信号z(k)。

(9)频域均衡器44对补偿后的基带信号z(k)进行均衡处理，得到新的均衡输出r(k)。

上述过程可以进行迭代多次，以更好的消除相位噪声。

最后应当说明的是，以上实施例仅用以说明本发明的技术方案而非限制，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神范围，其均应涵盖在本发明的权利要求范围当中。

Claims (16)

1.一种相位噪声的估计方法，其特征在于，包括如下步骤：

获取接收数据的均衡后符号r(k)对应的发送符号其中，k＝0,...,K-1，K为一帧接收数据包括的符号数目；

根据r(k)和获取公共相位噪声；

根据r(k)、和所述公共相位噪声获取残余相位噪声；以及

根据所述公共相位噪声、所述残余相位噪声和离散余弦变换的基函数矩阵获取相位噪声。

2.如权利要求1所述的估计方法，其特征在于，所述获取接收数据的均衡后符号r(k)对应的发送符号包括：

对r(k)进行硬判决，得到或者

对r(k)进行软判决，并对软判决结果进行重构得到

3.如权利要求1所述的估计方法，其特征在于，按照如下公式获取公共相位噪声：

${\widetilde{\mathrm{\θ}}}_{avg}=\arg \left(\underset{k=0}{\overset{K-1}{\Σ}}r\right(k\left){\tilde{a}}^{*}\right(k\left)\right)$

其中，为公共相位噪声，()^*表示求共轭运算，arg()表示求复数的辐角运算。

4.如权利要求3所述的估计方法，其特征在于，按照如下公式获取残余相位噪声：

$\widetilde{\mathrm{\φ}}\left(k\right)=\arg \left(r\right(k\left){\tilde{a}}^{*}\right(k\left)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}\right)$

其中，为第k个符号的残余相位噪声。

5.如权利要求4所述的估计方法，其特征在于，按照如下公式获取相位噪声：

$\widetilde{\mathrm{\θ}}\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{k^{\′}=0}{\overset{K-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\mathrm{\Ψ}}_K^T\right)}_{k,k^{\′}}\widetilde{\mathrm{\φ}}\left(k^{\′}\right)$

其中，为相位噪声，Ψ_K为离散余弦变换的基函数矩阵，为Ψ_K的转置，(Ψ_K)_k,n＝ψ_n(k)，ψ_n(k)为离散余弦变换的第n个基函数，n＝0,…,N-1，N为所述基函数的数目。

6.如权利要求3所述的估计方法，其特征在于，按照如下公式获取残余相位噪声：

${\widetilde{\mathrm{\φ}}}_B\left(m\right)=\frac{1}{B}\arg \left(\underset{b=0}{\overset{B-1}{\Σ}}r\left(b+mB\right)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}{\tilde{a}}^{*}\left(b+mB\right)\right)$

其中，一帧接收数据的K个符号被分成M段，每段包括B个符号，为第m段符号的残余相位噪声，m＝0,…,M-1。

7.如权利要求4所述的估计方法，其特征在于，按照如下公式获取相位噪声：

${\widetilde{\mathrm{\θ}}}_B\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{m=0}{\overset{M-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\left({\mathrm{\Ψ}}_{avg}^T{\mathrm{\Ψ}}_{avg}\right)}^{-1}{\mathrm{\Ψ}}_{avg}\right)}_{k,m}{\widetilde{\mathrm{\φ}}}_B\left(m\right)$

其中，为相位噪声，Ψ_K为离散余弦变换的基函数矩阵，(Ψ_K)_k,n＝ψ_n(k)，ψ_n(k)为离散余弦变换的第n个基函数，n＝0,…,N-1，N为所述基函数的数目，为Ψ_avg的转置， ${\left({\mathrm{\Ψ}}_{avg}\right)}_{m,n}=\frac{1}{B}\underset{b=0}{\overset{B-1}{\Σ}}{\mathrm{\ψ}}_n(b+mB).$

8.一种相位噪声的估计装置，其特征在于，包括：

发送符号估计单元，用于获取接收数据的均衡后符号r(k)对应的发送符号其中，k＝0,...,K-1，K为一帧接收数据包括的符号数目；

公共相位噪声估计单元，用于根据r(k)和获取公共相位噪声；

残余相位噪声估计单元，用于根据r(k)、和所述公共相位噪声获取残余相位噪声；以及

相位噪声估计单元，用于根据所述公共相位噪声、所述残余相位噪声和离散余弦变换的基函数矩阵获取相位噪声。

9.如权利要求8所述的估计装置，其特征在于，所述发送符号估计单元进一步用于：

对r(k)进行硬判决，得到或者

对r(k)进行软判决，并对软判决结果进行重构得到

10.如权利要求8所述的估计装置，其特征在于，所述公共相位噪声估计单元进一步用于，按照如下公式获取公共相位噪声：

${\widetilde{\mathrm{\θ}}}_{avg}=\arg \left(\underset{k=0}{\overset{K-1}{\Σ}}r\right(k\left){\tilde{a}}^{*}\right(k\left)\right)$

其中，为公共相位噪声，()^*表示求共轭运算，arg()表示求复数的辐角运算。

11.如权利要求10所述的估计装置，其特征在于，所述残余相位噪声估计单元进一步用于，按照如下公式获取残余相位噪声：

$\widetilde{\mathrm{\φ}}\left(k\right)=\arg \left(r\right(k\left){\tilde{a}}^{*}\right(k\left)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}\right)$

其中，为第k个符号的残余相位噪声。

12.如权利要求11所述的估计装置，其特征在于，所述相位噪声估计单元进一步用于，按照如下公式获取相位噪声：

$\widetilde{\mathrm{\θ}}\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{k^{\′}=0}{\overset{K-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\mathrm{\Ψ}}_K^T\right)}_{k,k^{\′}}\widetilde{\mathrm{\φ}}\left(k^{\′}\right)$

其中，为相位噪声，Ψ_K为离散余弦变换的基函数矩阵，为Ψ_K的转置，(Ψ_K)_k,n＝ψ_n(k)，ψ_n(k)为离散余弦变换的第n个基函数，n＝0,…,N-1，N为所述基函数的数目。

13.如权利要求10所述的估计装置，其特征在于，所述残余相位噪声估计单元进一步用于，按照如下公式获取残余相位噪声：

${\widetilde{\mathrm{\φ}}}_B\left(m\right)=\frac{1}{B}\arg \left(\underset{b=0}{\overset{B-1}{\Σ}}r\left(b+mB\right)e^{-j{\widetilde{\mathrm{\θ}}}_{avg}}{\tilde{a}}^{*}\left(b+mB\right)\right)$

其中，一帧接收数据的K个符号被分成M段，每段包括B个符号，为第m段符号的残余相位噪声，m＝0,…,M-1。

14.如权利要求13所述的估计装置，其特征在于，所述相位噪声估计单元进一步用于，按照如下公式获取相位噪声：

${\widetilde{\mathrm{\θ}}}_B\left(k\right)={\widetilde{\mathrm{\θ}}}_{avg}+\underset{m=0}{\overset{M-1}{\Σ}}{\left({\mathrm{\Ψ}}_K{\left({\mathrm{\Ψ}}_{avg}^T{\mathrm{\Ψ}}_{avg}\right)}^{-1}{\mathrm{\Ψ}}_{avg}\right)}_{k,m}{\widetilde{\mathrm{\φ}}}_B\left(m\right)$

其中，为相位噪声，Ψ_K为离散余弦变换的基函数矩阵，(Ψ_K)_k,n＝ψ_n(k)，ψ_n(k)为离散余弦变换的第n个基函数，n＝0,…,N-1，N为所述基函数的数目，为Ψ_avg的转置， ${\left({\mathrm{\Ψ}}_{avg}\right)}_{m,n}=\frac{1}{B}\underset{b=0}{\overset{B-1}{\Σ}}{\mathrm{\ψ}}_n(b+mB).$

15.一种包括如权利要求8～14任一项所述的估计装置的接收机。

16.一种包括如权利要求15所述的接收机的通信设备。