一种基于KT-LVT的抗符号跳变直接扩频信号动态聚焦方法技术

本发明专利技术公开了一种基于KT‑LVT的抗符号跳变直接扩频信号动态聚焦方法,本发明专利技术采用基于延迟相关抗符号跳变的KT‑LVT动态聚焦方法，在消除伪码相位距离徙动和多普勒频率徙动的同时，消除了符号跳变的影响，实现了低信噪比、高动态场景下直扩信号动态聚焦。与非搜索类动态聚焦方法相比，本方法中非线性运算造成的信噪比损失小，处理增益高，具有良好的抗噪声性能，更加适用于弱信号的捕获；同时，由于本方法不需要实时遍历补偿本地伪码速率和多普勒频率，也无需对符号跳变沿进行搜索，所以相对于非迭代非搜索类动态聚焦方法本方法运算量低、存储量小、实时性高。

【技术实现步骤摘要】
一种基于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)]2；(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，然后对本文档来自技高网...

【技术保护点】
一种基于KT‑LVT的抗符号跳变直接扩频信号动态聚焦方法，其特征在于，包括如下步骤：步骤一、将接收端的复基带信号表示为sr(Ts,Tk)，本地伪码序列表示为sl(Ts,Tk)；Ts＝sTc为段内延迟时间，Tk＝kTc为段间延迟时间，Tc为部分相关时间，s＝0,1,2,…,2N‑1,k＝0,1,2,…,K K为部分相关段数；N为部分相关点数；步骤二、将复基带信号sr(Ts,Tk)与本地伪码序列sl(Ts,Tk)进行相关处理，得到部分相关函数矩阵xcor(τ,Tk),τ为本地伪码序列的初始时间延迟；步骤三、将部分相关函数矩阵xcor(τ,Tk)进行平方处理得到

【技术特征摘要】
1.一种基于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变换后的矩阵进行二维鉴频，获得动态聚焦结果。2.如权利要求1所述的一种基于KT-LVT的抗符号跳变直接扩频信号动态聚焦方法，其特征在于，步骤一至步骤五由以下步骤具体实现：(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)...

【专利技术属性】
技术研发人员：申宇瑶，汪文韬，王永庆，吴嗣亮，于秀丽，马淑芬，陈岩，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：北京,11