[发明专利]一种有色噪声的干扰消除方法、装置、电子设备及存储介质

权利要求书

1.一种有色噪声的干扰消除方法，其特征在于，包括：

利用有色噪声协方差矩阵，初始化干扰消除相关参数；

所述初始化干扰消除相关参数，是计算复合矩阵G＝HAP，对有色噪声的协方差矩阵Σ做Cholesky分解，得到下三角矩阵Δ，满足Σ＝ΔΔ^H，其中，(·)^-H表示矩阵的共轭转置；利用Cholesky分解，有色噪声n能够表示为n＝Δw，其中，向量w＝Δ^-1n是对有色噪声白化变换后得到的高斯白噪声，其服从复高斯分布，均值为协方差矩阵为其中，表示大小为O_sN的单位方阵；

初始化迭代次数I＝0，设定最大迭代次数I_max，设定GAMP最大迭代次数T_max；发送符号向量β的第m个星座符号β_m的先验分布为p_a(β_m＝χ_q)＝1/Q，m＝0,...,M-1，q＝1,...,Q，其中，Q是调制阶数，χ_q是星座集合中的第q个星座点；

根据所述干扰消除相关参数，使用GAMP算法对接收信号进行干扰消除，得到编码比特的外部对数似然比；

对编码比特的外部对数似然比使用BCJR译码，得到编码比特的先验对数似然比；

根据编码比特的先验对数似然比，计算星座符号的先验概率；

将星座符号的先验概率回代入所述使用GAMP算法对接收信号进行干扰消除的步骤，进行迭代循环，并在达到设定的最大迭代次数后，输出判决信息比特，完成干扰消除。

2.根据权利要求1所述的有色噪声的干扰消除方法，其特征在于，所述干扰消除方法用于线性调制系统，该系统中发送星座符号所构成的向量为β，大小为N×1，其中，N是对向量β预编码后得到向量的大小；接收端的基带接收信号为：y≈HAPβ+n，其中，H是大小为O_sN×O_sN的多径信道和低通滤波器组合得到的矩阵，A是大小为O_sN×N的上采样矩阵，O_s是过采样倍数，P是大小为N×M的预编码矩阵，向量n表示有色噪声，其服从复高斯分布，均值为协方差矩阵为Σ，二者在接收端均已知，是大小为O_s×1的全零向量。

3.根据权利要求1所述的有色噪声的干扰消除方法，其特征在于，所述GAMP干扰消除算法，具体包括：

初始化一阶导数构成的向量其中，[·]内的值表示GAMP的迭代次数；初始化符号向量β的外部均值外部方差所构成的向量μ^e,β[0]＝∞1_M×1，其中，1_M×1表示大小为M×1的全1矩阵，∞表示无穷大的数值；根据概率分布q＝1,...,Q，计算发送符号向量的后验均值和方差m＝0,...,M-1，其中，和分别表示和μ^e,β[0]的第m个元素；初始化白噪声w的后验均值外部方差所构成的向量初始化GAMP迭代次数t＝0；

计算线性调制发送信号的先验均值向量和先验方差向量μ^p[t]：其中，|·|²表示对括号内的矩阵逐元素取模的平方，表示矩阵间逐元素相乘，向量向量

计算一阶和二阶导数构成的向量和μ^s[t]：其中，./表示矩阵间逐元素相除；

计算星座符号向量β的外部均值和方差μ^e,β[t]：

根据概率分布q＝1,...,Q，m＝0,...,M-1，计算发送符号向量的后验均值和方差μ^β[t]；

计算白噪声向量w的外部均值和方差μ^e,w[t]：

计算白噪声向量w的后验均值和方差μ^d,w[t]：