CN103188190B - 一种测量载波信号干扰的方法和系统 - Google Patents

Abstract

本发明公开了一种测量载波信号干扰的方法和系统。该方法包括：接收载波信号的多个有效样本；基于所述载波信号的所述多个有效样本，通过迭代计算获得所述干扰信号的常数部分的量值和零均值部分的均方差；计算干扰信号的总功率值。采用这种测量方法，考虑了干扰信号中常数部分，适用于载波干扰为非零均值的情况，所得到的载波信号干扰的测量结果相比现有方法更为准确。

Description

一种测量载波信号干扰的方法和系统

技术领域

本发明涉及移动通信领域，特别涉及测量载波信号干扰的方法和系统。

背景技术

多载波系统是通过把数据流分解为若干个子比特流，每个子数据流具有很低的比特速率，用低速率多状态符号去调制相应的子载波，构成多路低速率符号串并行发送。多载波系统中具有代表性的调制方式是正交频分复用(OFDM)技术，其子载波在频段上均匀分布并相互正交。子载波间可以重叠，从而提高频谱利用率，已经成为第四代移动通信系统的物理层候选技术。

在实际系统中，OFDM系统的调制和解调可以用离散傅里叶逆变换和离散傅里叶变换实现。在发送端，通过N点IDFT运算(N为OFDM子载波数)，将待发送的信号变换为时间域信号，经过上变频调制后，发送到信道。在接收端，对接收信号进行相干解调，得到基带信号后进行N点DFT。各资源块(或子载波)，这样可以把多径频率选择性信道转化为N个独立的平坦衰落信道，每个子载波对应一个子信道。

OFDMA系统中，载波增益和载波干扰水平是衡量各子载波信道质量的关键因素。而现有载波干扰测量方法将载波干扰看成是零均值的随机过程。

申请人对现有的载波干扰测量方法进行了研究，发现现有的载波干扰测量方法存在一定的测量误差。这是因为，现有载波干扰测量方法将载波干扰看成是零均值的随机过程。而在实际应用中，由于邻区干扰等因素，载波干扰不一定满足零均值的假设。因此，采用现有的载波干扰测量方法所得到的结果并不准确，有必要研究针对非零均值载波干扰的测量方法。

发明内容

本发明的发明人发现现有的载波干扰测量方法存在一定的测量误差。本发明的一个目的是提供一种新的技术方案，考虑载波信号的非零均值情况，能够准确测量载波干扰信号的总功率值。

根据本发明的第一方面，提供测量载波信号干扰的方法。该方法包括：接收载波信号的多个有效样本，该载波信号包括发送信号部分和干扰信号部分，该干扰信号部分包括常数部分和零均值部分；基于该载波信号的多个有效样本，通过迭代计算获得干扰信号的常数部分的量值和零均值部分的均方差；由干扰信号的常数部分的量值和零均值部分的均方差得到干扰信号的总功率值。

优选地，干扰信号的常数部分的迭代初始值I₀是采用如下公式计算的：其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值；干扰信号零均值部分的均方差的迭代初始值是采用如下公式计算的：其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值；进行迭代计算，以得出I_k+1和σ_k+1 ²。

对于发送信号采用BPSK方式调制的情况，优选地，在所述迭代计算中，

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\tanh \left(\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\right),$

${\mathrm{\σ}}_{k+1}^2=1+\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}({(x_n-I_k)}^2-\tanh \left(\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\right)2(x_n-I_k)),$ 其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为干扰信号的常数部分的第k次迭代值，为干扰信号零均值部分的均方差的第k次迭代值。

对于发送信号采用QAM方式进行调制的情况，优选地，在所述迭代计算中，

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\frac{{\mathrm{\Σ}}_{m=1}^M(2m-1)S_{m,n}(I_k,{\mathrm{\σ}}_k^2)}{{\mathrm{\Σ}}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)},$

${\mathrm{\σ}}_{k+1}^2={\mathrm{\σ}}_0^2-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\left(\frac{{\mathrm{\Σ}}_{m=1}^M(4m(m-1)C_{m,n}(I_k,{\mathrm{\σ}}_k^2)-2(2m-1)(x_n-I_k)S_{m,n}(I_k,{\mathrm{\σ}}_k^2))}{{\mathrm{\Σ}}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)}\right)$ 其中，

M表示QAM实部(虚部)的星座图阶数，

$C_{m,n}(I,{\mathrm{\σ}}^2)=\exp (-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\cosh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

$S_{m,n}(I,{\mathrm{\σ}}^2)=\mathrm{exo}(-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\sinh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为干扰信号的常数部分的第k次迭代值，为所述干扰信号零均值部分的均方差的第k次迭代值。

优选地，比较|I_k+1-I_k|与第一门限值，比较|σ_k+1 ²-σ_k ²|与第二门限值；当|I_k+1-I_k|小于第一门限值，且|σ_k+1 ²-σ_k ²|小于第二门限值时，结束该迭代计算，I_k+1和σ_k+1 ²即为通过迭代计算获得的所述干扰信号的常数部分的量值和零均值部分的均方差。

优选地，干扰信号的总功率是通过如下公示计算的：干扰信号的总功率值＝I_k+1 ²+σ_k+1 ²。

根据本发明的第二方面，提供一种测量载波信号干扰的系统。该系统包括：载波信号接收单元，用于接收载波信号的多个有效样本，该载波信号包括发送信号部分和干扰信号部分，该干扰信号部分包括常数部分和零均值部分；迭代计算单元，基于载波信号的多个有效样本，通过迭代计算获得干扰信号的常数部分的量值和零均值部分的均方差；总功率计算单元，由干扰信号的常数部分的量和零均值部分的均方差得到干扰信号的总功率值。

优选地，迭代计算单元采用如下公式计算干扰信号的常数部分的迭代初始值I₀：其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值；迭代计算单元采用如下公式计算干扰信号零均值部分的均方差的迭代初始值σ₀ ²：

其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值；迭代计算单元进行迭代计算，以得出I_k+1和σ_k+1 ²。

可选地，对于采用BPSK方式调制的发送信号，在迭代计算中，

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\tanh \left(\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\right),$

${\mathrm{\σ}}_{k+1}^2=1+\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}({(x_n-I_k)}^2-\tanh \left(\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\right)2(x_n-I_k)),$ 其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为所述干扰信号的常数部分的第k次迭代值，为干扰信号零均值部分的均方差的第k次迭代值；

对于采用QAM方式调制的发送信号，在迭代计算中，

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\frac{{\mathrm{\Σ}}_{m=1}^M(2m-1)S_{m,n}(I_k,{\mathrm{\σ}}_k^2)}{{\mathrm{\Σ}}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)},$

${\mathrm{\σ}}_{k+1}^2={\mathrm{\σ}}_0^2-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\left(\frac{{\mathrm{\Σ}}_{m=1}^M(4m(m-1)C_{m,n}(I_k,{\mathrm{\σ}}_k^2)-2(2m-1)(x_n-I_k))}{{\mathrm{\Σ}}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)}\right),$ 其中，

M表示QAM实部(虚部)的星座图阶数，

$C_{m,n}(I,{\mathrm{\σ}}^2)=\exp (-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\cosh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

$S_{m,n}(I,{\mathrm{\σ}}^2)=\exp (-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\sinh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为干扰信号的常数部分的第k次迭代值，为所述干扰信号零均值部分的均方差的第k次迭代值。

优选地，迭代计算单元还用于：比较|I_k+1-I_k|与第一门限值，比较|σ_k+1 ²-σ_k ²|与第二门限值；当|I_k+1-I_k|小于第一门限值，且|σ_k+1 ²-σ_k ²|小于第二门限值时，结束迭代计算，I_k+1和σ_k+1 ²即为通过迭代计算获得的干扰信号的常数部分的量值和零均值部分的均方差，则干扰信号的总功率值＝I_k+1 ²+σ_k+1 ²。

本发明的一个优点在于，所计算的干扰信号的总功率值为干扰信号的常数部分的量值和零均值部分的均方差。这种测量方法是针对载波干扰信号非零均值的情况，考虑了干扰信号中常数部分。因此，所得到的载波信号干扰的测量结果更为准确。

通过以下参照附图对本发明的示例性实施例的详细描述，本发明的其它特征及其优点将会变得清楚。

附图说明

构成说明书的一部分的附图描述了本发明的实施例，并且连同描述一起用于解释本发明的原理。

参照附图，根据下面的详细描述，可以更加清楚地理解本发明，其中：

图1是示出根据本发明的一个实施例测量载波信号干扰的方法的流程图；

图2A示出了采用本发明的方法和对比方法的估计误差，图2B是采用本发明的方法的迭代收敛测试图。

图3是示出根据本发明的一个实施例的测量载波信号干扰的系统的结构图。

具体实施方式

现在将参照附图来详细描述本发明的示例性实施例。应注意到：除非另外具体说明，否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值不限制本发明的范围。

同时，应当明白，为了便于描述，附图中所示出的各个部分的尺寸并不是按照实际的比例关系绘制的。

以下对至少一个示例性实施例的描述实际上仅仅是说明性的，决不作为对本发明及其应用或使用的任何限制。

对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论，但在适当情况下，所述技术、方法和设备应当被视为说明书的一部分。

在这里示出和讨论的所有示例中，任何具体值应被解释为仅仅是示例性的，而不是作为限制。因此，示例性实施例的其它示例可以具有不同的值。

应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步讨论。

图1示出了根据本发明的测量载波信号干扰的方法的一个实施例的流程图。

在步骤S101中，接收载波信号的多个有效样本。

所接收的载波信号可包括发送信号部分和干扰信号部分。其中，干扰信号部分可包括常数部分和零均值部分。

有效样本的数目可根据实际需要进行确定，在此不做限定。例如，如何如何。

在步骤S102中，计算干扰信号的常数部分的量值和零均值部分的均方差。

基于所接收的载波信号的多个有效样本，可以通过迭代计算获得干扰信号的常数部分的量值和零均值部分的均方差。

对于干扰信号的常数部分的迭代初始值I₀可以采用如下公式进行计算：

其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值。

对于干扰信号零均值部分的均方差的迭代初始值可采用如下公式计算：

其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值。

然后进行迭代计算，以得出I_k+1和σ_k+1 ²。

既可以设定迭代次数，当达到迭代次数时，即可完成迭代计算。也可以对迭代结果进行判断。当迭代结果达到一定条件时，即可完成迭代计算。本实施例以后一方式进行详细说明。

在迭代计算中，对迭代结果进行第一判断和第二判断。具体地，可设有第一门限值和第二门限值，比较|I_k+1-I_k|与第一门限值的大小，以及比较|σ_k+1 ²-σ_k ²|与第二门限值的大小。

当|I_k+1-I_k|小于第一门限值，且|σ_k+1 ²-σ_k ²|小于第二门限值时，结束迭代计算。此时，I_k+1和σ_k+1 ²即为通过迭代计算获得的所述干扰信号的常数部分的量值和零均值部分的均方差。

可以根据计算条件以及误差要求合理设定迭代计算中的门限值。例如，门限值的取值范围可以是10^-3～10^-5。在本实施例中，所采用的门限值为10^-3。

需要说明的是，发送信号可以采用不同方式进行调制。调制方式可包括，但不限于，BPSK方式和QAM方式。因此，对于采用不同方式调制的发送信号，迭代计算中的常数部分的k次迭代值I_k和零均值部分的均方差的第k次迭代值的表达式是不同的。下面进行详细说明。

对于采用BPSK方式进行调制的发送信号，在迭代计算中，可以采用

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\tanh \left(\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\right),$

${\mathrm{\σ}}_{k+1}^2=1+\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}({(x_n-I_k)}^2-\tanh \left(\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\right)2(x_n-I_k)),$ 其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为所述干扰信号的常数部分的第k次迭代值，为所述干扰信号零均值部分的均方差的第k次迭代值。

而对于QAM方式进行调制的发送信号，在迭代计算中，可以采用：

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\frac{{\mathrm{\Σ}}_{m=1}^M(2m-1)S_{m,n}(I_k,{\mathrm{\σ}}_k^2)}{{\mathrm{\Σ}}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)},$

${\mathrm{\σ}}_{k+1}^2={\mathrm{\σ}}_0^2-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\left(\frac{{\mathrm{\Σ}}_{m=1}^M(4m(m-1)C_{m,n}(I_k,{\mathrm{\σ}}_k^2)-2(2m-1)(x_n-I_k)S_{m,n}(I_k,{\mathrm{\σ}}_k^2))}{{\mathrm{\Σ}}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)}\right)$ 其中，

M表示QAM实部(虚部)的星座图阶数，

$C_{m,n}(I,{\mathrm{\σ}}^2)=\exp (-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\cosh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

$S_{m,n}(I,{\mathrm{\σ}}^2)=\exp (-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\sinh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为干扰信号的常数部分的第k次迭代值，为干扰信号零均值部分的均方差的第k次迭代值。

为了清楚描述本发明的技术方案和效果，下面结合图2A至图2B对本发明的测量载波信号干扰的方法和效果进行详细描述。

在该例子中，给出计算机仿真的测量结果。所用到的仿真参数如下：发送符号是二进制随机信号，确定性干扰I＝0.5，噪声的方差是0.2512，对应的信噪比为6dB。

在本例子中，进行了200次仿真。每次的发送信号和噪声都随机产生，然后求参数估计的平均误差。

图2A中纵坐标所示的估计误差是用均方误差进行衡量的，即估计误差＝E{(估计结果-实际参数)²}。其中，“E{}”表示对多次重复试验求平均值。

图中的对比方法只考虑随机噪声误差，应用基于统计量的方法进行计算。具体地，不考虑干扰信号的常数部分，不进行迭代计算，将零均值部分的均方差的初始值作为估计结果。

在图2A中，图线A为采用现有方法的估计误差，图线B为采用本发明的测量方法的估计误差。从图2A中可以看出，采用对比方法所得到的估计误差远大于采用本发明进行测量的估计误差。

在图2B中，图线C为干扰信号常数部分的量值估值；图线D为干扰信号零均值部分的均方差的估值。由图2B可以看出，本发明所采用的迭代计算在10次迭代后已经收敛。

在步骤S103中，由干扰信号的常数部分的量值和零均值部分的均方差得到干扰信号的总功率值。

具体地，干扰信号的总功率值＝I_k+1 ²+σ_k+1 ²。

采用本发明的载波干扰测量方法，所计算的干扰信号的总功率值为干扰信号的常数部分的量值和零均值部分的均方差。这种测量方法是针对载波干扰信号非零均值的情况，考虑了干扰信号中常数部分。因此，所得到的载波信号干扰的测量结果更为准确。

本发明还公开了一种测量载波信号干扰的系统。图3示出了根据本发明的测量载波信号干扰的系统的一个实施例的结构示意图。

该系统可包括载波信号接收单元11、迭代计算单元12和总功率计算单元13。其中，

载波信号接收单元11接收载波信号的多个有效样本。所接收的载波信号包括发送信号部分和干扰信号部分。该干扰信号部分可包括常数部分和零均值部分。

迭代计算单元12基于载波信号接收单元11所接收的多个有效样本，通过迭代计算获得干扰信号的常数部分的量值和零均值部分的均方差。

迭代计算单元12采用如下公式计算干扰信号的常数部分的迭代初始值I₀：

其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值；

迭代计算单元12采用如下公式计算干扰信号零均值部分的均方差的迭代初始值

其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值；

迭代计算单元12通过迭代计算得出I_k+1和σ_k+1 ²。

对于不同调制方式的发送信号，迭代计算中所用到的常数部分的第k次迭代值和干扰信号零均值部分的均方差的第k次迭代值是不同的。

具体地，对于采用BPSK方式调制的发送信号，在迭代计算中，

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\tanh \left(\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\right),$

${\mathrm{\σ}}_{k+1}^2=1+\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}({(x_n-I_k)}^2-\tanh \left(\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\right)2(x_n-I_k)),$ 其中，N为所接收的载

波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为干扰信号的常数部分的第k次迭代值，为干扰信号零均值部分的均方差的第k次迭代值。

对于采用QAM方式调制的发送信号，在所述迭代计算中，

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\frac{{\mathrm{\Σ}}_{m=1}^M(2m-1)S_{m,n}(I_k,{\mathrm{\σ}}_k^2)}{{\mathrm{\Σ}}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)},$

${\mathrm{\σ}}_{k+1}^2={\mathrm{\σ}}_0^2-\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}\left(\frac{{\mathrm{\Σ}}_{m=1}^M(4m(m-1)C_{m,n}(I_k,{\mathrm{\σ}}_k^2)-2(2m-1)(x_n-I_k)S_{m,n}(I_k,{\mathrm{\σ}}_k^2))}{{\mathrm{\Σ}}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)}\right),$ 其中，

M表示QAM实部(虚部)的星座图阶数，

$C_{m,n}(I,{\mathrm{\σ}}^2)=\exp (-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\cosh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

$S_{m,n}(I,{\mathrm{\σ}}^2)=\exp (-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\sinh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为所述干扰信号的常数部分的第k次迭代值，为所述干扰信号零均值部分的均方差的第k次迭代值。

需要说明的是，对发送信号的调制方式不限于以上两种，相应地，在迭代计算中，干扰信号的常数部分的量值的表达式以及零均值部分的均方差的表达式可进行相应调整。

对于迭代计算中的结束，既可以设定一定迭代次数，当迭代计算达到所设定的计算次数时，即可完成迭代计算。也可以对迭代计算结果进行判断。当迭代计算结果满足一定条件时，完成迭代计算。

在本实施例中，迭代计算单元12对|I_k+1-I_k|与第一门限值的大小以及|σ_k+1 ²-σ_k ²|与第二门限值的打下分别进行比较。当|I_k+1-I_k|小于第一门限值，且|σ_k+1 ²-σ_k ²|小于第二门限值时，结束迭代计算，I_k+1和σ_k+1 ²即为通过迭代计算获得的所述干扰信号的常数部分的量值和零均值部分的均方差。

总功率计算单元13根据干扰信号的常数部分的量值和零均值部分的均方差计算干扰信号的总功率值。

具体地，干扰信号的总功率值＝I_k+1 ²+σ_k+1 ²。

至此，已经详细描述了根据本发明的测量载波信号干扰的方法和系统。为了避免遮蔽本发明的构思，没有描述本领域所公知的一些细节。本领域技术人员根据上面的描述，完全可以明白如何实施这里公开的技术方案。

虽然已经通过示例对本发明的一些特定实施例进行了详细说明，但是本领域的技术人员应该理解，以上示例仅是为了进行说明，而不是为了限制本发明的范围。本领域的技术人员应该理解，可在不脱离本发明的范围和精神的情况下，对以上实施例进行修改。本发明的范围由所附权利要求来限定。

Claims (10)

1.一种测量载波信号干扰的方法，其特征在于，包括：

接收载波信号的多个有效样本，所述载波信号包括发送信号部分和干扰信号部分，所述干扰信号部分包括常数部分和零均值部分；

基于所述载波信号的所述多个有效样本，通过迭代计算获得所述干扰信号的常数部分的量值和零均值部分的均方差；

由所述干扰信号的常数部分的量值和零均值部分的均方差得到所述干扰信号的总功率值。

2.如权利要求1所述的方法，其特征在于，

所述干扰信号的常数部分的迭代初始值I₀是采用如下公式计算的：

其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值；

所述干扰信号零均值部分的均方差的迭代初始值是采用如下公式计算的：

其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值；

进行迭代计算，以得出I_k+1和σ_k+1 ²。

3.如权利要求2所述的方法，其特征在于，所述发送信号是采用二进制相移键控BPSK方式调制的，在所述迭代计算中，

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}\tanh \left(\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\right),$

${\mathrm{\σ}}_{k+1}^2=1+\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}({\left(x_n-I_k\right)}^2-\tanh (\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\left)2\right(x_n-I_k\left)\right),$ 其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为所述干扰信号的常数部分的第k次迭代值，为所述干扰信号零均值部分的均方差的第k次迭代值。

4.如权利要求2所述的方法，其特征在于，所述发送信号是采用正交振幅调制QAM方式调制的，在所述迭代计算中，

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}\frac{{\Σ}_{m=1}^M(2m-1)S_{m,n}(I_k,{\mathrm{\σ}}_k^2)}{{\Σ}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)},$

${\mathrm{\σ}}_{k+1}^2={\mathrm{\σ}}_0^2-\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}\left(\frac{{\mathrm{\Σ}}_{m=1}^M\left(4m\right(m-1\left)C_{m,n}\right(I_k,{\mathrm{\σ}}_k^2)-2(2m-1\left)\right(x_n-I_k\left)S_{m,n}\right(I_k,{\mathrm{\σ}}_k^2\left)\right)}{{\mathrm{\Σ}}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)}\right)$ 其中，

M表示QAM实部和虚部的星座图阶数，

$C_{m,n}(I,{\mathrm{\σ}}^2)=\exp (-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\cosh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

$S_{m,n}(I,{\mathrm{\σ}}^2)=\exp (-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\sinh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为所述干扰信号的常数部分的第k次迭代值，为所述干扰信号零均值部分的均方差的第k次迭代值。

5.如权利要求3或4所述的方法，其特征在于，还包括：

比较|I_k+1-I_k|与第一门限值，比较|σ_k+1 ²-σ_k ²|与第二门限值；

当|I_k+1-I_k|小于第一门限值，且|σ_k+1 ²-σ_k ²|小于第二门限值时，结束所述迭代计算，I_k+1和σ_k+1 ²即为通过第k+1次迭代计算获得的所述干扰信号的常数部分的量值和零均值部分的均方差。

6.如权利要求5所述的方法，其特征在于，所述干扰信号的总功率是通过如下公示计算的：

干扰信号的总功率值＝I_k+1 ²+σ_k+1 ²。

7.一种测量载波信号干扰的系统，其特征在于，包括：

载波信号接收单元，用于接收载波信号的多个有效样本，所述载波信号包括发送信号部分和干扰信号部分，所述干扰信号部分包括常数部分和零均值部分；

迭代计算单元，基于所述载波信号的所述多个有效样本，通过迭代计算获得所述干扰信号的常数部分的量值和零均值部分的均方差；

总功率计算单元，由所述干扰信号的常数部分的量值和零均值部分的均方差得到所述干扰信号的总功率值。

8.如权利要求7所述的系统，其特征在于，

所述迭代计算单元采用如下公式计算所述干扰信号的常数部分的迭代初始值I₀：

其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值；

所述迭代计算单元采用如下公式计算所述干扰信号零均值部分的均方差的迭代初始值σ₀ ²：

其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值；

所述迭代计算单元进行迭代计算，以得出I_k+1和σ_k+1 ²。

9.如权利要求8所述的系统，其特征在于，

对于采用二进制相移键控BPSK方式调制的发送信号，在所述迭代计算中，

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}\tanh \left(\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\right),$

${\mathrm{\σ}}_{k+1}^2=1+\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}({\left(x_n-I_k\right)}^2-\tanh (\frac{(x_n-I_k)}{{\mathrm{\σ}}_k^2}\left)2\right(x_n-I_k\left)\right),$ 其中，N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为所述干扰信号的常数部分的第k次迭代值，为所述干扰信号零均值部分的均方差的第k次迭代值；

对于采用正交振幅调制QAM方式调制的发送信号，在所述迭代计算中，

$I_{k+1}=I_0-\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}\frac{{\mathrm{\Σ}}_{m=1}^M(2m-1)S_{m,n}(I_k,{\mathrm{\σ}}_k^2)}{{\mathrm{\Σ}}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)},$

${\mathrm{\σ}}_{k+1}^2={\mathrm{\σ}}_0^2-\frac{1}{N}\underset{n=1}{\overset{N}{\Σ}}\left(\frac{{\mathrm{\Σ}}_{m=1}^M\left(4m\right(m-1\left)C_{m,n}\right(I_k,{\mathrm{\σ}}_k^2)-2(2m-1\left)\right(x_n-I_k\left)S_{m,n}\right(I_k,{\mathrm{\σ}}_k^2\left)\right)}{{\mathrm{\Σ}}_{m=1}^M{C}_{m,n}(I_k,{\mathrm{\σ}}_k^2)}\right),$ 其中，

M表示QAM实部和虚部的星座图阶数，

$C_{m,n}(I,{\mathrm{\σ}}^2)=\exp (-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\cosh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

$S_{m,n}(I,{\mathrm{\σ}}^2)=\exp (-\frac{2m(m-1)}{{\mathrm{\σ}}^2})\sinh \left(\frac{(2m-1)(x_n-I)}{{\mathrm{\σ}}^2}\right),$

N为所接收的载波信号的有效样本数目，x_n为第n个时刻所接收的样本值，I_k为所述干扰信号的常数部分的第k次迭代值，为所述干扰信号零均值部分的均方差的第k次迭代值。

10.如权利要求9所述的系统，其特征在于，所述迭代计算单元还用于：

比较|I_k+1-I_k|与第一门限值，比较|σ_k+1 ²-σ_k ²|与第二门限值；

当|I_k+1-I_k|小于第一门限值，且|σ_k+1 ²-σ_k ²|小于第二门限值时，结束所述迭代计算，I_k+1和σ_k+1 ²即为通过迭代计算获得的所述干扰信号的常数部分的量值和零均值部分的均方差，

则干扰信号的总功率值＝I_k+1 ²+σ_k+1 ²。