基于ILSP-CMA的同步DS-CDMA信号伪码序列和信息序列联合盲估计制造技术

本发明专利技术请求保护一种基于ILSP‑CMA的同步DS‑CDMA信号伪码序列和信息序列联合盲估计，属于信号处理技术领域。利用已知的用户数，伪码列周期和信息码序列周期构造接收的同步DS‑CDMA信号模型，然后对信号的协方差矩阵进行特征值分解，得到对应的特征向量即为信号子空间，并将接收信号投影到信号子空间，建立对数似然函数模型，最后利用迭代最小二乘投影恒模算法求解信息矩阵和伪码矩阵，最终实现伪码序列和信息序列的联合盲估计。本方法可以在较低信噪比下精确地估计出同步DS‑CDMA信号信号的伪码序列和信息序列，且该方法易于操作，具有重要的工程实际意义。

Blind estimation of pseudo code sequence and information sequence of synchronous DS-CDMA signal based on ilsp-cma

【技术实现步骤摘要】
基于ILSP-CMA的同步DS-CDMA信号伪码序列和信息序列联合盲估计
本专利技术涉及直接序列扩频通信信号处理，具体为一种基于迭代最小二乘投影恒模(ILSP-CMA)算法的同步DS-CDMA信号伪码序列和信息序列联合盲估计方法。
技术介绍
直接序列扩频(DirectSequenceSpreadSpectrum,DSSS)通信信号是利用伪码序列(或称为PN码)技术将原始窄带数字信号扩展到一个非常宽的频谱范围后再进行传输，使得信号频谱被展宽并降低信号功率谱密度，可工作于负信噪比环境即信号淹没于噪声之中，具有强的抗干扰、抗多径、低截获概率、多址复用等优点，在民用和军事通信领域中得到了广泛的应用。如现在的第三代通信中的CDMA系统就采用了直接序列扩频技术，而在美国最新的联合战术分配系统(JTIDS)中，也利用了直接序列扩频和跳频的混合方式进行通信。针对同步直接序列码分多址(DS-CDMA)信号，不同用户使用不同的伪码序列，且各用户之间的伪码序列相互正交，使其具有多址(Multi-Access,MA)特性。在非合作通信中，要实现对接收的DS-CDMA信号的分析，需要对信号中的伪码序列和信息序列进行获取。因此开展DS-CDMA信号伪码序列和信息序列的研究具有重要意义。DSSS信号根据用户数不同可分为单用户直接序列扩频信号和多用户直接序列扩频信号。针对单用户信号的伪码序列和信息序列盲估计问题已经得到了广泛而深入的研究，主要方法有特征分解法、三阶相关法、信息子空间估计法和神经网络算法等。针对多用户信号的研究，文献(YaoY,PoorV.EavesdroppinginthesynchronousCDMAchannel:anEM-basedapproach[J].IEEETransactionsonSignalProcessing,2015,49(8):1748-1756.)提出了基于EM算法的DS-CDMA信号盲估计算法，这也是目前针对同步DS-CDMA信号较为有效的盲估计算法，但是该算法需要利用估计出来的伪码序列对信号进行解扩得到信息序列估计并且该算法的复杂度会随着用户数增加呈指数增长。文献(AvitzourD.DetectionofasynchronousCDMAwithunknownuserwaveforms[J].IEEESignalProcessingLetters,2004,11(2):209-211.)提出利用迭代最小二乘投影(IterativeLeastSquaresProjective,ILSP)算法，通过循环迭代对扩频矩阵和信息矩阵进行估计，但在扩频序列未知情况下，信息序列估计性能受限。鉴于此，本专利技术将ILSP-CMA算法应用在DS-CDMA信号估计上，通过迭代计算可直接估计得到伪码序列和信息序列。
技术实现思路
本专利技术所要解决的技术问题，是针对当前DS-CDMA信号在低信噪比下，难以对伪码序列和信息序列进行联合盲估计的问题，提出一种基于ILSP-CMA算法的估计方法。该方法不需要利用估计的伪码序列对信号进行解扩来得到信息序列，并且该方法利用信息序列的恒模特性，最小化输出信号的模量变化，能够快速收敛到门限，解决了在低信噪比条件下，此类信号伪码序列和信息序列联合盲估计的难题。本专利技术解决上述技术问题而提出的技术方案的实施办法为：首先对接收信号的协方差矩阵进行特征分解，并根据用户数得到信号子空间，然后将接收信号投影到信号子空间，并建立关于信息矩阵和伪码矩阵的对数似然函数，最后利用迭代最小二乘法对函数求解，最终得到信息序列和伪码序列。接收的DS-CDMA信号经过同步和基带处理后，可表示为：其中，K表示用户个数，M表示信息码个数，也是伪码序列周期数，Ak表示第k个用户的信号幅度，v(n)表示均值为零、方差为σ2的高斯白噪声，L和TC分别表示伪码序列的长度和码片宽度。{bk(m)∈±1,m∈Z+}为码元宽度为Ts的第k个用户的信息序列且服从等概率分布，ck(t)表示第k个用户的伪码波形，N为信号样本长度，N＝ML。(1)式写为矩阵形式为Y＝CAB+V(2)其中C＝[c1,c2,…,cK]∈RL×K为伪码矩阵，A＝diag(A1,A2,…,AK)为信号幅度对角矩阵，B＝[b1,b2,…,bK]T∈RK×M为信息码矩阵，V代表方差为σ2的高斯白噪声矩阵。由式(2)可构建Y的协方差矩阵为RY＝E[YmYHm]＝CA2CH+σ2I(3)对于(3)式可由样本协方差近似得到。然后对进行特征值分解，可得到前K个最大特征值对应的特征向量所表示的信号子空间US的估计值K表示用户数。将观测矩阵Y投影到其信号子空间对Y进行降维处理得到其中，X∈RK×M,仍然为高斯白噪声矩阵。X的对数似然函数可由下式给出当(5)式最后一项取得最小值时对应的H和B为该似然函数所对应的最大似然估计，此时等价于求解下面的最小值问题：其中，||·||F表示矩阵的Frobenius范数。对于式(6)的求解，总体思路是利用最小二乘法进行求解。令首先给定H的初始值利用最小二乘法求出的最小值对应的并将每一个元素bk(m)投影到单位圆，以保证bk(m)为单位长度。固定得到的再次利用最小二乘法求出的最小值对应的H。重复该过程，直至收敛或者达到最大迭代次数。附图说明图1本专利技术ILSP-CMA算法流程示意图；图2本专利技术不同用户数下伪码序列盲估计随信噪比变化性能图；图3本专利技术不同用户数下信息序列盲估计随信噪比变化性能图；图4本专利技术算法的收敛速度曲线图；具体实施方式以下结合附图和具体实例，对本专利技术的实施作进一步的描述。图1所示为本专利技术ILSP-CMA算法流程示意图，具体步骤如下：(1)初始化H0∈RK×K，p＝0，收敛门限为threshold；(2)p＝p+1，更新B矩阵：并将Bp投影：(3)更新H矩阵：(4)计算(5)判断当th小于threshold(或者达到最大迭代次数)时停止迭代，否则重复步骤(2)、(3)、(4)。通过ILSP-CMA算法最终可得到信息码序列B的估计H的估计并由其得到第k个用户信号幅度和伪码序列的估计：其中，下标p表示第p次迭代的结果，符号表示矩阵的伪逆，sgn{·}表示求符号运算。图2和图3为不同用户数下的伪码序列和信息序列估计性能图。二者都是在信息码个数为200，用户数分别为4、5、6，信噪比为-10dB—0dB时，通过300次蒙特卡洛仿真得到的误码率曲线图。由图2可知，本专利技术所提算法在信噪比为-6dB和-2dB下能实现对4路用户和7路用户的伪码序列的正确估计。由图3可知，本专利技术所提算法在信噪比为-4dB和-2dB下能实现对4路用户和7路用户的信息序列的正确估计。同时，也可以看出在同一信噪比下，用户数越少，误码率越低，算法性能增强。图4为不同用户数下的算法收敛速度曲本文档来自技高网...

【技术保护点】
1.一种基于ILSP-CMA的同步DS-CDMA信号伪码序列和信息序列联合盲估计，具体步骤：首先对接收信号的协方差矩阵进行特征分解，并根据用户数得到信号子空间，然后将接收信号投影到信号子空间，并建立关于信息矩阵和伪码矩阵的对数似然函数，最后利用迭代最小二乘法对函数求解，最终得到信息序列和伪码序列。/n

【技术特征摘要】
1.一种基于ILSP-CMA的同步DS-CDMA信号伪码序列和信息序列联合盲估计，具体步骤：首先对接收信号的协方差矩阵进行特征分解，并根据用户数得到信号子空间，然后将接收信号投影到信号子空间，并建立关于信息矩阵和伪码矩阵的对数似然函数，最后利用迭代最小二乘法对函数求解，最终得到信息序列和伪码序列。


2.根据...

【专利技术属性】
技术研发人员：张天骐，喻盛琪，赵健根，张刚，李鑫凯，
申请(专利权)人：重庆邮电大学，
类型：发明
国别省市：重庆;50