基于PARAFAC模型的新型DS-CDMA盲接收机

摘要：根据平行因子(PARAFAC)模型，研究DS-CDMA盲多用户检测算法。将直接三线性分解算法(DTLD)与三线性交替最小二乘(TALS)算法结合，提出一种新的DTALS-PARAFAC盲接收机，解决了三线性交替最小二乘(TALS)算法中因为初始值估计不当引起的收敛速度差的问题。仿真结果表明，与TALS-PARAFAC接收机相比，DTALS-PARAFAC接收机改善了误码率性能，并且具有更快的收敛速度。
关键词：盲信号处理；平行因子分析；直接三线性分解；三线性交替最小二乘；直接序列码分多址

    在信号处理和通信领域，盲信号处理已经逐渐成为研究热点。肓信号处理算法能够在未知信道状态和系统参数的情况下实现对信号的检测和参数的辨识。PARAFAC(parallelfactor，平行因子)方法最早是用于心理计量学和应用化学领域的一种三维数据分析方法。21世纪初，Nicholas D．Sidiro-poulos将其应用到盲信号处理领域。PARAFAC方法能够充分利用信号的代数性质和分集特性对接收的信号进行处理，通过二三维数据的拟合得到信号处理中需要的各种参数。因其良好的性能，近几年来相继被应用到信号处理的各种领域，如阵列信号处理、正交频分复用(orthogonal frequency divisionmultiplexng，OFDM)系统、多输入多输出(multi-input multi-output，MIMO)系统。
    文献中提出一种性能收敛于非盲MMSE的基于PARAFAC模型的盲接收机。这种盲接收机综合利用了不同用户信号在时间、空间、扩频码方面的分集特性，在不需要已知扩频码、多径传输、DOA校准信息的情况下将用户有用信息重建出来。本文将对这种基于TALS算法的盲PARAFAC接收机进行改进，将TALS算法与DTLD(direct triIinear decomposition)算法结合起来，提出一种新的DTALS算法，改进了TALS算法中用随机矩阵进行初始估计引起的不稳定性和收敛速度差的问题。仿真结果表明，使用DTALS算法具有更好的拟合精度和收敛速度，在一定程度上降低了算法的复杂度，提高了算法性能。
    文中用到的符号与算子说明如下：diag(．)是对角化算子，表示括号中向量为对角元素形成的对角阵；为Kronecker积；为Khatri-Rao积，它是Kronecker积的列形式；AB=[a1b1，…，aRbR]；‖·‖F为矩阵的Frobenius范数，简称F范数；(·)+表示矩阵Moore-Penro se逆矩阵；(·)T表示矩阵的转置。

1 数据模型
    Nicholas D．Sidiropoulos将PARAFAC模型应用到DS-CDMA系统中，将接收信号构造成一个在时间、扩频、天线3个方向分集的三线性立体模型，如图1所示。

    假设发送端有F个用户发送数据，系统采用基带信号进行直接传输，信号经过一个加性高斯白信道进入接收端，接收端基站(BS)有I个接收天线，用码片速率对基带天线输出进行采样，总共收集到J个符号周期的数据。此时，接收端接收到的信号为：

    根据文献可知，在一定条件下，PARAFAC模型能够唯一分解。

2 DTALS-PARAFAC接收机
    文献中提出的TALS-PARAFAC盲接收机性能上接近于最小均方误差接收机。TALS-PARAFAC算法中使用随机矩阵来初始化矩阵A、B，基于式子(3)、(4)、(5)，根据最小二乘原理迭代实现PARAFAC模型的三线性分解。然而当初始矩阵A、B估计不当时，迭代过程容易陷入局部最优的计算“沼泽”，收敛相当缓慢，并可能产生错误的解。
    DTLD(direct trilinear decomposition)是一种直接三线性分解方法，也是一种基于PARAFAC模型的三维分解方法，它是非迭代的，具有直接快速的优点。但是没有明确的最优化界限，在信噪比不高的情况下，若三维数据不严格服从三线性模型或随机误差较大，有可能出现无意义的虚数解，可靠性差。
    在应用化学领域，已有学者提出并验证，在基于PARAFAC模型的分析化学算法中，使用适当的初始值对矩阵进行初始估计(如奇异值分解矩阵、特征分析)可以有效地提高算法的收敛速度和拟合精度。因此，在基于PARAFAC模型的TALS算法巾融入DTLD算法，能够解决TALS算法中初始值选取不当导致的收敛缓慢并且可能出现错误的解的问题。仿真结果表明，使用DTLD对TALS算法进行初始化，不仅能够提高原本算法的收敛速度，而且能够在一定程度上提高算法的精度，使分解出的数据更加接近于真实数据。

3 基于DTLD的交替最小二乘算法

    其中ε为一个很小的值(通常取1e-6)。
    算法收敛后，得到估计矩阵A、B、C。DTALS算法利用步骤1、2、3得到的矩阵进行迭代初始化，使初始值更加精确，能够减少迭代的次数，降低运算的复杂度。仿真结果表明，DTALS算法在拟合精度上比原有的TALS算法有所提高，并且具有更快的收敛速度。

4 仿真结果与分析
    通过仿真，分析DTALS-PARAFAC接收机的误码率(biterror rate，BER)性能和收敛性能，将仿真精度与仿真速度等因素综合考虑，假设M=4个用户的DS-CDMA系统，一次处理的用户信息样本长度为I=50，扩频码采用Hadamard码，长度为K=7，天线数J=4，噪声为加性高斯白噪声，用户与天线的衰落因子服从高斯随机分布。Monte Carlo仿真1 000次取其平均误码率。
4．1 DTALS-PARAFAC、TALS-PARAFAC、MMSE接收机的BER性能比较
    将DTALS-PARAFAC盲接收机与文献中TALS-PARAFAC接收机的误比特率性能比较。同时对非盲的线性MMSE接收机进行仿真以作参考。MMSE接收机需要已知信道衰落矩阵A和扩频码矩阵C。仿真结果如图2所示。

    由图2可以看出，在未知信道衰落矩阵A和扩频码矩阵C的基础上，DTALS-PARAFAC接收机的性能要优于TALS-PARAFAC接收机，且与非盲MMSE接收机性能相近。说明基于PARAFAC的盲接收机性能优越，两种算法的误码性能反映了其迭代算法的拟合精度，即接收机的误码率越低，拟合精度越高，反之亦然。因此，在多用户DS-CDMA系统中，DTALS算法拟合精度更高。
4．2 DTALS-PARAFAC接收机和TALS-PARAFAC接收机的收敛速度比较
    图3给出了在不同信噪比条件下，两种接收机的平均迭代次数。可以看出2种接收机的迭代次数均随着信噪比的增加而下降，在相同信噪比的情况下，DTALS-PARAFAC的迭代次数远小于TALS-PARAFAC接收机。因此，在保证误码率的前提下，DTALS算法具有更优越的收敛速度，性能更加稳定。

5 结论
    文中将FARAFAC模型与多用户多天线DS-CDMA系统相结合，在原有的TALS算法的基础上结合一种DTLD算法，提出了一种新的DTALS-PARAFAC盲接收机，通过DTLD算法将更加符合模型分解的初始值引入到TALS迭代过程中，优化了算法的性能。结果表明，DTALS-PARAFAC接收机在处理DS-CDMA信号的过程中具有更好的拟合精度和收敛速度，更适用于DS-CDMA系统。