一种基于KT‑LVT的抗符号跳变直接扩频信号动态聚焦方法与流程

本发明属于航天测控通信技术领域，具体涉及一种基于kt-lvt的抗符号跳变直接扩频信号动态聚焦方法。

背景技术：

直接序列扩频(dsss)信号具有抗截获和抗干扰能力强的优点，可以实现码分多址，广泛应用于航天测控通信、卫星导航和移动通信等领域。微弱直扩信号的捕获需要长时间有效的能量积累。而增长信号积累时间长度，则需要考虑符号跳变和收发双方相对运动引入的接收信号动态的影响。其中，积累时间内的符号跳变会造成对应不同符号的信号能量相互抵消，出现谱峰分裂现象；接收信号动态会造成信号能量在伪码相位维和频率维发生徙动和扩散。此时常规积累方法无法实现能量的有效积累。

高动态直扩信号捕获的关键在于通过能量聚焦获得良好的信号能量积累性能。由于高阶动态模型下非搜索类动态聚焦方法采用相关运算获得高阶核函数，其抗噪声性能恶化严重，不适用于弱信号的捕获。典型的非迭代非搜索类动态聚焦方法包括：基于一阶相邻互相关函数的动态聚焦方法和基于keystone变换-lv(吕氏)变换(kt-lvt)的动态聚焦方法，这两种方法均可用于对二阶动态模型信号的进行聚焦。其中，基于一阶相邻互相关函数的动态聚焦方法直接对接收信号进行延迟相关运算，该非线性运算在输入信噪比低时会造成较大的信噪比损失，不适用于弱信号的捕获。而kt-lvt方法在延迟和瞬时时间二维积累，可以减小非线性运算带来的损失，具有较高的信噪比增益，适合弱信号捕获。基于此，在2015年ietradar,sonar&navigation第9卷第5期第600页至607页由jintian等人发表的“parameterestimationofmanoeuvringtargetsbasedonsegmentintegrationandlv’stransform”一文中，提出了kt-lvt动态聚焦方法，对动态场景下雷达信号进行动态聚焦，实现了雷达的信号检测和参数估计。该方法应用在雷达信号处理领域，不需要考虑符号调制对动态聚焦方法的影响。但在直扩信号模型下，由于长时间积累过程中存在符号跳变的现象，常规动态聚焦方法会出现能量抵消和谱峰分裂的问题，因此要kt-lvt动态聚焦方法用于直扩信号，还需要考虑符号调制的影响，对kt-lvt算法进行改进，设计基于kt-lvt的抗符号跳变的动态聚焦方法。

技术实现要素：

有鉴于此，本发明的目的是提供一种基于平方运算抗符号跳变的kt-lvt的动态聚焦方法,利用dbzp消除部分相关时间段内的数据跳变的影响，并在kt-lvt算法的基础上进行改进，利用平方运算消除符号跳变造成的部分相关时间段间能量抵消问题，具有较高的信噪比增益，适合微弱直扩信号的捕获。

一种基于kt-lvt的抗符号跳变直接扩频信号动态聚焦方法，包括如下步骤：

步骤一、将接收端的复基带信号表示为sr(ts,tk)，本地伪码序列表示为sl(ts,tk)；ts＝stc为段内延迟时间，tk＝ktc为段间延迟时间，tc为部分相关时间，s＝0,1,2,…,2n-1,k＝0,1,2,…,kk为部分相关段数；n为部分相关点数；

步骤二、将复基带信号sr(ts,tk)与本地伪码序列sl(ts,tk)进行相关处理，得到部分相关函数矩阵xcor(τ,tk),τ为本地伪码序列的初始时间延迟；

步骤三、将部分相关函数矩阵xcor(τ,tk)进行平方处理得到

步骤四、对xsq(τ,k)进行keystone变换，消除一阶动态造成的伪码相位徙动，得到xsq(τ,k')；

步骤五、忽略二阶动态距离维徙动，将部分相关结果xsq(τ,k')近似为：

f0为射频频率，rc(τ)为伪码互相关函数，tk'＝tk(f0+fτ)/f0，a0为雷达和目标之间的初始距离，a1为一阶目标动态参数，a2为二阶目标动态参数，a'为keystone变换后信号幅值，xsq(τ,k')的对称延迟相关函数为：

rs(τ,k',m)＝xsq(τ,k'+m+b)[xsq(τ,k'-m-b)]*

其中，fτ表示变量τ经傅里叶变换后的频域变量；b为固定延迟，m为延迟变量，m＝0,1,2,…,m-1；m-1为最大延迟变量；

步骤六、对延迟相关函数进行keystone变换，再对keystone变换后的矩阵进行二维鉴频，获得动态聚焦结果。

步骤一至步骤五由以下步骤具体实现：

(1)设置行搜索编号l的初始值为l＝0，设置列搜索编号n的初始值为n＝0；

(2)对本地伪码序列进行采样，并将采样序列分段，每段n个采样点，共k+2m+2b段，并在每段段末补n个0，将数据段按行存入2n×(k+2m+2b)维的矩阵作为本地时域矩阵，对本地时域矩阵的行向量进行2n点傅里叶变换，取共轭，得到本地频域共轭矩阵

(3)对复基带信号进行采样，并对采样序列分段处理，每段数据长n点，共2(k+m+b)段，然后在每段段尾补n点数据，每段段尾数据为下一段数据的前n点，将数据段按行存入2n×2(k+m+b)维的矩阵作为接收拓展矩阵，然后对接收拓展矩阵的行向量进行2n点的傅里叶变换，得到复基带信号的频域值矩阵；

(4)截取所述复基带信号的频域值矩阵的第l+1行至l+(k+2m+2b)行作为接收信号频域值矩阵sr(fτ,k)；

(5)将接收信号频域值矩阵sr(fτ,k)与本地频域共轭矩阵对应位置相乘，得到频域值矩阵；

(6)对复基带信号的频域值矩阵xcor(fτ,k)中fτ维进行傅里叶逆变换，截取部分相关矩阵的前n列，得到部分相关函数时域值矩阵xcor(τ,k)，其中，

τ＝0,1,2,…,n-1；

(7)对部分相关函数时域值矩阵xcor(τ,k)进行平方处理，得到平方矩阵xsq(τ,k)＝[x(τ,k)]²；

(8)对平方矩阵xsq(τ,k)中τ维进行傅里叶变换，得到频域值矩阵xsq(fτ,k)；

(9)对矩阵xsq(fτ,k)按列进行keystone变换，得到keystone变换后频域值矩阵xsq(fτ,k')；

(10)对频域值矩阵xsq(fτ,k')中fτ维进行逆傅里叶变换，得到运动补偿后的部分相关矩阵xsq(τ,k')；

(11取出相关矩阵xsq(τ,k')的第n列的数据xsq(n,k')；

(12)根据xsq(n,k')计算对称延迟相关函数矩阵：

其中，1≤m≤m，0≤k'≤k-1；

(13)将延迟相关函数矩阵rs(n,k',m)按列进行keystone变换，得到时域值矩阵rs(n,k″,m)；

(14)将rs(n,k″,m)的列向量做m点傅里叶变换，得到矩阵rs(n,k″,fm)；

(15)对rs(n,k″,fm)矩阵的行向量做k点傅里叶变换，得到动态聚焦矩阵rs(n,fk″，fm)；

(16)对动态聚焦矩阵rs(n,fk″,fm)取模，找出|rs(n,fk″,fm)|峰值所在的位置(ia,iv)，并计算相应的一阶动态参数和二阶动态参数；然后计算矩阵|rs(n,fk″,fm)|的峰值和均值之比作为检测判决量vdec，将其与设定门限vt进行比较。若检测判决量vdec大于门限vt，则输出输出动态参数和动态聚焦完成标志，n×l+n为码相位搜索结果，hfsiv/(2k)为多普勒估计结果；反之，调整列搜索序号n＝n+1，后对其进行判断，若n小于n，则跳转至步骤(11)；若n等于n，则调整行搜索编号l＝l+1，然后对l进行判断，若l小于k，则跳转至步骤(4)；若l等于k，则输出动态聚焦失败标志。

本发明具有如下有益效果：

本发明采用基于延迟相关抗符号跳变的kt-lvt动态聚焦方法，在消除伪码相位距离徙动和多普勒频率徙动的同时，消除了符号跳变的影响，实现了低信噪比、高动态场景下直扩信号动态聚焦。与非搜索类动态聚焦方法相比，本方法中非线性运算造成的信噪比损失小，处理增益高，具有良好的抗噪声性能，更加适用于弱信号的捕获。

同时，由于本方法不需要实时遍历补偿本地伪码速率和多普勒频率，也无需对符号跳变沿进行搜索，所以相对于非迭代非搜索类动态聚焦方法本方法运算量低、存储量小、实时性高。

附图说明

图1为本发明的方法流程图。

具体实施方式

下面结合附图并举实施例，对本发明进行详细描述。

基于平方运算的kt-lvt方法原理如下：

设n为部分相关点数，tc为部分相关时间，k为部分相关段数，且ktc大于等于一个伪码周期，在无噪声干扰的情况下，接收端复基带信号可表示为

其中，ts＝stc为段内延迟时间，tk＝ktc为段间延迟时间，c为光速，a为信号幅值，r0(tk)为发射机与接收机之间的距离

其中，a0为雷达和目标之间的初始距离，a1为一阶目标动态参数，a2为二阶目标动态参数。

本地信号可表示为

经过相关处理后得到部分相关函数矩阵

其中，f0为射频频率，τ为本地信号的初始时间延迟，rc(τ)为伪码互相关函数，rd(τ)为符号互相关函数。由于按照公式(4)每次仅能计算矩阵中的某一个元素，为了节省运算量，可用基于ft方法实现，二者的处理结果相同。由于复基带信号的符号位和伪码是相参关系，则当接收信号和本地信号伪码对齐时，在部分相关时间内不存在符号位跳变，由此消除了瞬时时间维符号跳变的影响。

对xcor(τ,tk)进行平方处理得到

对xsq(τ,k)进行kt，消除一阶动态造成的伪码相位徙动，即

其中，tk'＝tk(f0+fτ)/f0，a'为kt后信号幅值。

在假定二阶动态距离维徙动可忽略的情况下，xsq(τ,tk')可看做是一个关于tk'的定常幅线性调频信号，进而用lvt进行多普勒频率维二阶动态联合聚焦。部分相关结果经kt可近似为

其对称延迟相关函数为

rs(τ,k',m)＝xsq(τ,k'+m+b)[xsq(τ,k'-m-b)]^*(8)

其中，为b固定延迟，m为延迟变量，m＝0,1,2,…,m-1。

对延迟相关函数进行kt和二维鉴频，获得动态聚焦结果

其中，tk″＝[h(m+b)tn]tk'，tm＝mtc，h为尺度因子，且h越大，lvt方法中kt的性能越好，但动态参数a2的估计精度差，需要折衷选取h的值。

通过公式(9)可以看到信号能量集中在τ＝a0/c,fm＝4a1f0,fk″＝4a1f0/h的位置，进一步的可以求得估计值a0＝cτ,fm＝4a1f0,fk″＝4a1f0/ha1＝fm/(4f0),a2＝fk″h/(4f0)。

由于本方法是在某个固定的延迟时间τ下计算的延迟相关函数时(参考公式(8))，因此需要对τ进行遍历搜索，当且仅当τ＝a0/c时，才能获得能量聚焦结果。

可见，kt-lvt方法可以消除符号跳变和高阶动态的影响，信号能量集中在目标单元格内。

本发明结构框图如图1所示，包括数据预处理模块、本地频域值矩阵计算模块、行矩阵ift1模块、平方模块、行矩阵ft1模块、行矩阵ift2模块、频域kt、相关结果存储模块、延迟相关函数计算模块，时域kt模块，矩阵列ft模块和矩阵行ft2模块、检测判决模块和逻辑控制模块。具体步骤如下：

(1)预置逻辑控制模块内行搜索编号l＝0输出给数据预处理模块，列搜索编号为n＝0输出给相关结果存储模块。

(2)本地伪码频域值矩阵计算模块对本地伪码序列按采样率fs进行采样，将采样序列分段，每段n个采样点，共k+2m+2b段，并在每段段末补n个0，将数据段按行存入2n×(k+2m+2b)的矩阵作为本地时域矩阵，对本地时域矩阵的行向量进行2n点ft计算，取共轭，得到本地频域共轭矩阵sl^*(fτ,k)输出给乘积模块。

(3)数据预处理模块按采样率fs对复基带信号进行采样、分段处理，每段数据长n点，共2(k+m+b)段，然后在每段段尾补n点数据，每段段尾数据为下一段数据的前n点，将数据段按行存入2n×2(k+m+b)的矩阵作为接收拓展矩阵，然后对拓展矩阵的行向量进行2n点的ft计算，得到复基带信号的频域值矩阵。

(4)数据预处理模块根据逻辑控制模块输出的行搜索编号l截取复基带信号的频域值矩阵的第l+1行～l+(k+2m+2b)行做为接收信号频域值矩阵sr(fτ,k)输出给乘积模块。

(5)乘积模块接收矩阵sr(fτ,k)并将其与本地频域共轭矩阵对应位置相乘，得到频域值矩阵xcor(fτ,k)，输出给矩阵行ift1模块。

(6)矩阵行ift1模块接收复基带信号的频域值矩阵xcor(fτ,k)，对fτ进行ift，截取部分相关矩阵的前n列，得到部分相关函数时域值矩阵xcor(τ,k)(即令

τ＝0,1,2,…,n-1)输出给平方模块。

(7)平方模块接收部分相关函数时域值矩阵xcor(τ,k)，计算平方矩阵

xsq(τ,k)＝[x(τ,k)]²输出给矩阵行ft1模块。

(8)矩阵行ft1模块接收平方矩阵xsq(τ,k)，对τ进行ft，得到频域值矩阵xsq(fτ,k)输出给频域kt模块。

(9)频域kt模块接收矩阵xsq(fτ,k)，按列进行kt，输出kt后频域值矩阵xsq(fτ,k')。

(10)矩阵行ift模块接收频域值矩阵xsq(fτ,k')，对fτ进行ift计算，得到运动补偿后的部分相关矩阵xsq(τ,k')，输出给相关结果存储模块。

(11)相关结果存储模块接收xsq(τ,k')，按逻辑控制模块输出的列搜索编号n，取出相关矩阵第n列的数据xsq(n,k')给延迟相关函数计算模块。

(12)延迟相关函数计算模块接收信号xsq(n,k')，计算延迟相关函数矩阵

d(m,k')＝xsq(n,k'+2m+2b)×ml^*(n,k')(6)

其中，1≤m≤m，0≤k'≤k-1。步骤五矩阵rs(n,k',m)是三维矩阵，由于步骤(11)中变量n被设为定值，在步骤(12)中该矩阵减少了一个维度，步骤(12)中的rs(n,k',m)是步骤五中rs(n,k',m)的一部分，每一次迭代对应不同的n值，迭代完成后即可得到完整的三维矩阵。

(13)时域kt模块接收矩阵d(m,k')，按列进行kt，得到时域值矩阵d(m,k″)，输出给矩阵列ft模块。

(14)矩阵列ft模块接收到d(m,k″)后对该矩阵的列向量做m点ft，得到矩阵d(fm,k″)，并将该矩阵输出给矩阵行ft模块。

(15)矩阵行ft模块进一步对d(fm,k″)矩阵的行向量做k点行ft，得到矩阵d(fm,fk″)，将其输出给检测判决模块。

(16)检测判决模块接收步骤(11)输出的动态聚焦矩阵d(fm,fk″)，取模，找出|d(fm,fk″)|峰值所在的位置(ia,iv)，并计算相应的一阶和二阶动态参数，输出给逻辑控制模块。计算矩阵的峰值和均值之比作为检测判决量vdec，将其与门限vt进行比较。若检测判决量vdec大于门限vt则输出检测结果标志flag＝1；否则输出flag＝0。

(17)逻辑控制模块接收检测结果标志flag，对flag的值进行判断，若flag＝1，则输出动态参数和动态聚焦完成标志，n×l+n为码相位搜索结果，hfsiv/(2k)为多普勒估计结果；若flag＝0，调整列搜索序号n＝n+1，后对其进行判断，若n小于n，则跳转至步骤(11)；若n等于n，则调整行搜索编号l＝l+1，然后对l进行判断，若l小于k，则跳转至步骤(4)；若l等于k，则输出动态聚焦失败标志。

综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。