非周期长码直扩信号伪码估计方法技术

本发明专利技术涉及一种非周期长码直扩信号伪码估计方法。现有直扩信号伪码估计技术无法应用于接收信号长度不足一个长码周期的非周期长码直扩信号。本发明专利技术方法是：计算非周期长码直扩信号的三阶相关函数；然后构建三阶相关峰检测的二元假设检验模型和概率分布函数，通过循环去均值处理和拟合优度检验实现三阶相关峰粗提取和精检测；最后利用提取到的三阶相关峰估计信号长伪码m序列的本原多项式。基于m序列三阶相关峰对应多项式的最大公约式为m序列本原多项式的基本原理，本发明专利技术通过循环去均值处理和拟合优度检验更为精确地提取三阶相关峰，实现对非周期长码直扩信号的伪码盲估计。

【技术实现步骤摘要】

本专利技术属于通信对抗中直接序列扩频信号的盲参数估计领域，具体涉及非合作通信条件下非周期长码直扩信号伪码估计方法。
技术介绍
直接序列扩频(DirectSequenceSpreadSpectrum，DSSS)是应用较为广泛的一种扩展频谱系统，它具有抗干扰、抗多径、保密和码分多址等优点，在军事、导航和民用通信中得到广泛应用。非周期长码直扩信号的一个长码周期调制多个信息符号，结构复杂，且接收样本数不足一个长码周期，导致传统的扩频码估计方法很难估计非周期长码直扩信号的伪随机码(简称伪码)，给非合作通信中信息截获带来极大的挑战。研究表明高阶统计量对信息缺失不敏感，在信号检测与估计领域有着良好的应用前景。非周期长码直扩信号的扩频伪码一般采用m序列，因此可以利用m序列特有的三阶相关函数(TripleCorrelationFunction，TCF)特性实现非周期长码直扩信号的伪码盲估计。
技术实现思路
本专利技术的目的是针对非合作通信中无法估计非周期长码直扩信号的伪码的问题，提出一种基于循环去均值处理、拟合优度检验和三阶相关的非周期长码直扩信号伪码估计方法，从而解决了非周期长码直扩信号的伪码盲估计问题。本专利技术中非周期长码直扩信号伪码盲估计方法的步骤是：1、将非周期长码直扩信号以扩频码码片速率采样转化为基带信号，消除信息码影响后计算得到信号的三阶相关函数。2、将信号的三阶相关函数矩阵进行预处理和循环去均值处理，实现信号三阶相关峰的粗提取。3、对于粗提取得到三阶相关峰，利用信号三阶相关函数值的概率分布特性，通过拟合优度检验实现三阶相关峰的精检测。4、将检测到的三阶相关峰坐标点按共轭系分组，属于不同共轭系的峰值点坐标表示为多项式形式，并两两计算最大公约式，排除掉不合理的公因式后，出现次数最多的公因式即为长伪码m序列本原多项式。本专利技术通过循环去均值处理和拟合优度检验更为精确地提取三阶相关峰，再利用m序列三阶相关峰对应多项式的最大公约式为m序列本原多项式的基本原理，实现对非周期长码直扩信号的伪码盲估计。本专利技术利用非周期长码直扩信号中信息符号宽度远大于扩频码码片长度的特点，结合m序列移位相加特性，通过将信号延迟一个扩频码码片后与自身相乘，可以有效消除信息码对非周期长码直扩信号三阶相关函数矩阵估计的影响。本专利技术利用m序列三阶相关峰特性对估计得到的非周期长码直扩信号三阶相关函数矩阵进行预处理和循环去均值处理，粗提取得到信号三阶相关峰，降低了直接设置阈值提取峰值点的误差。本专利技术利用信号峰值点和非峰值点的三阶相关函数值概率分布的差异，构建三阶相关峰检测的二元假设检验模型，从粗提取到的三阶相关峰中排除伪峰。通过拟合优度检验，将极大提高峰值点检测的准确性，提高伪码估计性能。本专利技术将精确检测到的三阶相关峰按共轭系分组，对属于不同共轭系的峰值点多项式两两求最大公因式，可以提高伪码本原多项式估计的效率和准确率。本专利技术的有益效果是：1、利用三阶统计量对信息缺失不敏感的特性，通过三阶相关法解决了接收信号不足一个长码周期而无法估计伪码的问题。2、通过对信号三阶相关函数矩阵进行循环去均值处理，实现了三阶相关峰的初步提取，避免了传统门限值提取可能出现的漏检等问题。3、构建了非周期长码直扩信号三阶相关峰检测的二元假设检验模型，通过拟合优度KS检验实现三阶相关峰的准确检测。具体实施方式下面进一步详细说明本专利技术的实施步骤。步骤1，本专利技术非周期长码直扩信号伪码估计方法在使用时，首先将接收到的非周期长码直扩信号以扩频码码片速率采样后，转化为基带信号y(l)：y(l)＝Ad(l)c(l)＋v(l)(1)其中，l为采样时刻，l＝0,1,…,L-1；A、d(l)和c(l)分别表示接收信号的幅度、信息码、扩频长码的采样值；v(l)为零均值高斯白噪声，方差为σ2；扩频增益为G(等于信息符号宽度)，长扩频码周期为N，N＞＞G；接收信号长度为L，L＜N。将基带信号延迟一个扩频码码片后与原基带信号相乘，得到y1(l)：y1(l)＝y(l)y(l＋1)(2)在L×L范围内计算的三阶相关函数，得到信号的三阶相关函数：C^y1(p,q)=1LΣl=1Ly1(l)y1(l+p)y1(l+q)---(3)]]>其中，p和q为延迟量。步骤2，是关于p,q的二元函数，在p,q＝1,2,…,L处的L×L个取值可构成L×L维的三阶相关函数矩阵C：其中，cpq表示三阶相关函数矩阵函数C在p行q列处的取值，cp是矩阵C的第p行的行向量。三阶相关函数矩阵C中提取峰值点坐标可以描述为如下问题：给定起点p(p＝1,2,…,L)，在行向量cp中搜索是否存在相应的q，使得cpq对应三阶相关函数峰值。首先利用m序列三阶相关函数特性对三阶相关函数矩阵做预处理，将矩阵C中满足p＝q或p＝L或q＝L的元素置零。然后再对三阶相关函数矩阵循环去均值处理，具体方法为：①令α＝0，统计中非零元素个数②令α＝α＋1，对做去均值处理：cp(α+1)=cp(α)-1NΣq=1Lcpq(α)---(5)]]>若中的元素值小于零，则将其置为零，同时统计中非零元素的个数③若则循环结束，否则重复执行②。三阶相关函数矩阵C循环去均值处理后记为C1，检测C1中的I个非零元素，其行列号即为峰值点坐标，将其记为集合Ψ1，步骤3，令则非周期长码直扩信号伪码三阶相关峰检测可表示为如下的二元假设检验：H0:ω~N(-A6/N,σ22)H1:ω~N(A6,σ22)---(6)]]>其中，H0表示(p,q)是非峰值点的假设；H1表示(p,q)是峰值点的假设。由坐标点(p,q)可得L×L范围内的K个共轭系坐标点序列各点对应的三阶相关函数值序列为其中当(p,q)是非峰值点时，其所对应的三阶相关函数值序列样本相互独立，且服从正态分布，由于N＞＞A，所以，分布函数近似为：F0(t)=12πσ2∫-∞texp(-z2/2σ22)dz---(7)]]>因此判断(p,q)是否为峰值点可以等效为如下二元假设检验模型：H0:F(t)=F0(t)H1:F(t)≠F0(t)---(8)]]>给定的虚警概率Pf，计算判决门限为γPf=Δ-1(Pf)/(K+0.12+0.11K)---(9)]]>对粗估计得到本文档来自技高网...

【技术保护点】
非周期长码直扩信号伪码估计方法，其特征在于该方法包括以下步骤：(1)将非周期长码直扩信号以扩频码码片速率采样转化为基带信号，消除信息码影响后计算得到信号的三阶相关函数；(2)将信号的三阶相关函数矩阵进行预处理和循环去均值处理，实现信号三阶相关峰的粗提取；(3)对于粗提取得到三阶相关峰，利用信号三阶相关函数值的概率分布特性，通过拟合优度检验实现三阶相关峰的精检测；(4)将检测到的三阶相关峰坐标点按共轭系分组，属于不同共轭系的峰值点坐标表示为多项式形式，并两两计算最大公约式，排除掉不合理的公因式后，出现次数最多的公因式即为长伪码m序列本原多项式，该本原多项式唯一标识非周期长码直扩信号的伪码。

【技术特征摘要】
1.非周期长码直扩信号伪码估计方法，其特征在于该方法包括以下步
骤：
(1)将非周期长码直扩信号以扩频码码片速率采样转化为基带信号，
消除信息码影响后计算得到信号的三阶相关函数；
(2)将信号的三阶相关函数矩阵进行预处理和循环去均值处理，实现
信号三阶相关峰的粗提取；
(3)对于粗提取得到三阶相关峰，利用信号三阶相关函数值的概率分
布特性，通过拟合优度检验实现三阶相关峰的精检测；
(4)将检测到的三阶相关峰坐标点按共轭系分组，属于不同共轭系的
峰值点坐标表示为多项式形式，并两两计算最大公约式，排除掉不合理的公
因式后，出现次数最多的公因式即为长伪码m序列本原多项式，该本原多项
式唯一标识非周期长码直扩信号的伪码。
2.根据权利要求1所述的非周期长码直扩信号伪码估计方法，其特征
在于：步骤(1)将接收到的非周期长码直扩信号以扩频码码片速率采样后，
转化为基带信号y(l)：
y(l)＝Ad(l)c(l)+v(l)(1)
其中，l为采样时刻，l＝0,1,…,L-1；A、d(l)和c(l)分别表示接收信号
的幅度、信息码、扩频长码的采样值；v(l)为零均值高斯白噪声，方差为σ2；
扩频增益为G，长扩频码周期为N，N＞＞G；接收信号长度为L，L＜N；
将基带信号y(l)延迟一个扩频码码片后与原基带信号相乘，得到y1(l)：
y1(l)＝y(l)y(l+1)(2)
在L×L范围内计算的三阶相关函数，得到信号的三阶相关函数：
C^y1(p,q)=1LΣl=1Ly1(l)y1(l+p)y1(l+q)---(3)]]>其中，p和q为延迟量。
3.根据权利要求2所述的非周期长码直扩信号伪码估计方法，其特征
在于：步骤(2)具体是：
步骤(1)得到的信号三阶相关函数是关于p,q的二元函数，在
p,q＝1,2,…,L处的L×L个取值可构成L×L维的三阶相关函数矩阵C：
其中，cpq表示三阶相关函数矩阵函数C在p行q列处的取值，cp是矩阵
C的第p行的行向量；
首先利用m序列三阶相关函数特性对三阶相关函数矩阵C做预处理，将
矩阵C中满足p＝q或p＝L或q＝L的元素置零；
然后再对三阶相关函数矩阵循环去均值处理，具体方法为：
①令初始值α＝0，统计中非...

【专利技术属性】
技术研发人员：赵知劲，李淼，尚俊娜，徐春云，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江;33