本发明属于雷达
技术领域：
，特别涉及脉冲多普勒雷达的动目标检测技术。
背景技术：
：空时自适应处理(space-timeadaptiveprocessing,简称stap)是应用于机载地面动目标指示(groundmovingtargetindication，简称gmti)一项关键技术。其根据机载雷达的地物杂波呈现的空时二维耦合谱特性，对空间不同位置采集信号(空间采样信号)的处理就是利用波达方向信息进行区分的方向滤波，而同时对时域和空域采样信号进行处理，以期同时利用多普勒谱和波达方向信息来区分运动目标和静止的地杂波。stap技术可以应用于预警飞机，其中包含stap技术的机载预警雷达对动目标的检测能力在战争中意义重大，即使在和平时期，其对国家的安全也发挥着重要作用。stap技术虽然有卓越的性能，但是由于传统stap算法会产生较大的数据量，在发展之初，由于stap技术面临很大的计算压力而很少应用于实际雷达系统中。众所周知，stap技术自适应权值的计算是动目标检测的必要前提。但是在实际应用中的系统却很难实现，一是因为系统的输出信杂噪比(signaltointerferenceplusnoiseratio，sinr)是由估计杂波协方差矩阵的近距杂波环的数目决定，即估计杂波协方差矩阵的近距杂波环的数目必须遵循reed，mallett，brennan所提出准则，其要求的距离门数之多难满足，详见参考文献l.s.reed，j.d.mallett，andl.e.brennan,“rapidconvergencerateinadaptivearrays,”ieeetransactionsonaerospaceandelectronicsystems,vol.47,no.1pp.569-585,2011；二是自适应权值的计算需要对高维的杂波协方差矩阵进行估计和求逆，其运算量之大难处理。为了降低运算量，减少辅助距离门数据数量，提高stap在实际系统中的实用性，研究者对传统stap方法进行了改进，提出了多种改进stap方法。近年来，基于学习的stap方法成为研究者们研究的热点，代表性文献可见文献1：elkhatiba,assalehkandmirh,“space-timeadaptiveprocessingusingpatternclassification,”ieeetransactionsonsignalprocessing,vol63,no.3pp.766-779,2015该方法利用二项式分类器进行动目标检测，能够在距离门较少的情况下相对于传统stap有较大的性能提升。然而，该方法要求回波有较高的信杂比，当回波信号信杂比较低时，其输出性能较差。2016年，研究者提出一种基于支持向量回归机(supportvectorregression，svr)的方法(中国发明专利“一种基于支持向量机的空时自适应处理方法”，申请号201610814854.2)，该方法可以克服基于二项式stap方法的缺陷，在信杂比较低的情况下，仍可获得较好的检测性能。然而，该方法权值不稀疏，在检测动目标时存在运算量大的缺点。技术实现要素：本发明为了进一步提升stap性能，克服基于二项式分类器stap、基于svr的stap方法的不足，提出了一种基于稀疏贝叶斯学习的stap方法，从stap问题模型出发，将原始的动目标检测的滤波问题转化模式识别的问题，然后采用稀疏贝叶斯学习(sparsebayesianlearning,sbl)的方法实现动目标检测，在距离门较少的情况下，不仅有良好的输出性能，而且在检测时，运算量较小。本发明采用的技术方案是：一种基于稀疏贝叶斯学习的空时自适应处理方法，包括：s1、初始化pd雷达动目标检测系统的参数，包括：阵列天线的阵元个数，记做n；一个相干处理时间内发射脉冲个数，记做m；雷达发射信号波长，记做λ；接收信号的杂噪比，记做cnr；接收信号的信噪比，记做snr；接收信号的干噪比，记做jnr；雷达系统脉冲重复频率，记做prf；雷达平台的速度，记做v；阵列天线的阵元间隔，记做d；训练数据中每一个多普勒的训练次数，记做nc；待检测的距离门数，记做nt；待检测的回波数据矩阵为l，为nt行mn列的矩阵；s2、根据步骤s1中初始化的参数构建训练数据矩阵；s3、根据步骤s2的训练数据矩阵构建目标矩阵；s4、根据步骤s2得到的训练数据矩阵以及步骤s3得到的目标矩阵，通过稀疏贝叶斯学习的方法，计算得到权系数矩阵；s5、根据步骤s4得到的权系数矩阵，以及待检测的回波数据矩阵，检测各距离门的动目标。进一步地，步骤s2所述训练数据矩阵为：f＝[re(h)im(h)]；其中，f为(m+1)nc行2mn列的矩阵，h为(m+1)nc行mn列的矩阵，且h＝s+[t0]*；则f中共(m+1)nc个训练样本，上标*表示矩阵转置，s为mnc行mn列的矩阵，t为mnc行mn列的矩阵。更进一步地，所述矩阵s由下式确定：其中，sβ＝xβ+yβ+zβ，β＝1,…,(m+1)nc；sβ、xβ、zβ均为mn×1的列向量，sβ表示矩阵s中的元素，xβ表示杂波，yβ表示干扰，zβ表示噪声。更进一步地，所述矩阵t由下式确定：其中，为mn×1的列向量，表示第i个多普勒位置的第j个动目标的空时快拍，i＝1,…,m，j＝1,…nc。进一步地，步骤s3所述目标矩阵为：p＝[p1…pm…pm+1]；其中，pm为1×(m+1)nc列的向量，m＝1,…,m+1，pm的第ξ个元素为pm(ξ)，ξ＝1,…,(m+1)nc。更进一步地，所述pm(ξ)的值确定过程为：m＝1,…m时，目标矩阵p的第ξ行第m列元素pm(ξ)的值根据步骤s2所构建的训练数据矩阵f中第ξ行训练数据的动目标速度来确定；具体为：当该行训练数据中包含有动目标，且动目标的速度为时，pm(ξ)＝1；否则pm(ξ)＝0，为第m个多普勒类的归一化多普勒频率；m＝m+1时，目标矩阵p的第ξ行第m列元素pm(ξ)的值根据步骤s2所构建的训练数据矩阵f中第ξ行训练数据是否含有动目标确定；具体为：当该行训练数据不含有动目标，且vt＝0时，pm(ξ)＝1；否则pm(ξ)＝0。进一步地，所述步骤s4具体包括：s41、由下式计算权系数向量；pm(ξ)＝hξwm；其中，wm为2mn×1的列向量，m＝1,…m+1，hξ为1×2mn的行向量，ξ＝1,…,(m+1)nc，hξ为训练数据矩阵f的第ξ行向量，表示第ξ个训练样本，pm(ξ)是目标矩阵p的第m列的第ξ个元素；s42、根据步骤s41计算得到的权系数向量构建权系数矩阵：w＝[w1…wm…wm+1]；其中，w是2mn行m+1列的矩阵。进一步地，所述步骤s5具体为：矩阵q＝lw；当矩阵q中的元素且l＝1,…m时，矩阵q中的元素则判断第个距离门含有动目标，否则判断第个距离门无动目标；其中，l则为待检测距离门的回波数据构成的nt行2mn列的矩阵，w是2mn行m+1列的矩阵，q为nt行m+1列的矩阵，…nt，δ为门限值。本发明的有益效果：一种基于稀疏贝叶斯学习的空时自适应处理方法，充分利用空时自适应处理回波解调采样后空时快拍数据的特点，产生动目标空时快拍加入干扰(interference)环境，形成训练数据。将归一化多普勒频率进行平均划分，在目标矩阵中用1表示有目标，0表示无动目标的方式，巧妙的将空时自适应处理中的杂波抑制的问题进行了转化，不仅避免了传统stap方法中使用较多的辅助数据对协方差矩阵的估计，而且检测运算量小。此外，本发明检测权值一旦计算得到，适合该参数下其他同一相干处理时间内任何一个距离门的检测，因此，本发明结构简洁适用于实际中空时自适应处理中动目标检测等应用；本发明填补了现阶段在少量的回波数据的以及回波信杂比较小的情况下下无法正确检测检测动目标的空白。附图说明图1为本发明的方案流程图。具体实施方式为便于本领域技术人员理解本发明的技术内容，下面结合附图对本
发明内容进一步阐释。为方便描述本发明的内容，首先作一下术语定义：定义1、稀疏贝叶斯学习(sparsebayesianlearning,sbl)稀疏贝叶斯的数学模型为：其中，为样本的输入向量，xn为ns×1的列向量且ns为一个样本中特征向量的维数，n为大于或等于1的正整数，代表样本的个数。w权重为1×ns行向量，φ(x)为基底函数。令tn＝y(xn；w)+εn，其中为样本xn对应的训练值。当权重w中不为0的个数小于ns时，则称权重稀疏。稀疏贝叶斯学习(sbl)方法给出一种权重w的方法，使得权重稀疏。具体可参考论文michaele.tipping“sparsebayesianlearningandtherelevancevectormachine”journalofmachinelearningresearch,1(jun),pp211-244,2011。如图1所示为本发明的方案流程图，本发明的技术方案为：一种基于稀疏贝叶斯学习的空时自适应处理方法，包括：s1、初始化参数初始化pd雷达动目标检测系统的参数包括：阵列天线的阵元个数n＝10；一个相干处理时间内发射脉冲个数m＝10；雷达发射信号波长λ＝0.6m；接收信号的杂噪比cnr＝30db；接收信号的信噪比snr＝20db；接收信号的干噪比jnr＝10db；雷达系统脉冲重复频率prf＝1000hz；雷达平台的速度v＝150m/s；阵列天线的阵元间隔d＝0.3m；训练数据中每一个多普勒的训练次数nc＝10；待检测的距离门数nt＝50；待检测的回波为l，为50行200列的矩阵。s2、构建训练数据矩阵根据步骤s1中的初始化参数，构建训练数据矩阵，本实施例中训练数据矩阵f为110行200列的矩阵，根据下式构建：f＝[re(h)im(h)]；其中，h为110行100列的矩阵，且h＝s+[t0]*；则f中共110个训练样本，上标*表示矩阵转置，s为110行100列的矩阵，称为辅助数据(secondarydata)s＝[s1s2…s110]*；其中，sβ＝xβ+yβ+zβ，β＝1,…,(m+1)nc；sβ、xβ、zβ均为mn×1的列向量，sβ表示矩阵s中的元素，xβ表示杂波，yβ表示干扰，zβ表示噪声。xβ是一个100×1的列向量，可以通过下式得到其中，k是杂波单元总数；是第β个距离门内第k个杂波单元的幅度；表示空间频率为归一化多普勒频率为时杂波的空时快拍；其中，分别表示产生归一化多普勒频率为的时间快拍和空间频率为的空间快拍，θβ为第β个距离门时的俯仰角，φk为第k个杂波单元的方位角，d为阵列天线的阵元间隔，λ为雷达发射信号的波长，prf为雷达系统的脉冲重复频率，表示kronecker积，上标*表示矩阵的转置。yβ为一个100×1的列向量，可以由下式得到其中，bβ是10×1列向量，表示干扰幅度；aβ是10×1的列向量，表示干扰的空间快拍；其中，d＝0.3m为阵列天线的阵元间隔，λ＝0.6m为雷达发射信号的波长，μβ、分别为第β距离门干扰的俯仰角和方位角。zβ为一个100×1列向量zβ＝κβ；其中，κβ是一个100×1随机向量，其服从的分布类型通常为复高斯分布。详情参见参考文献j.ward，space-timeadaptiveprocessingforairborneradar.lincolnlab-oratory，mti，cam-bridge，ma，usa，tech.rep.1015，1994.17-24t为100行100列的矩阵其中，为100×1的列向量，表示第i个多普勒位置的第j个动目标的空时快拍；其中，ai为动目标的幅度，为动目标的第i类的归一化多普勒频率，vt为动目标的速度，λ＝0.6m为雷达发射信号的波长，prf＝1000hz为雷达系统的脉冲重复频率。为第i类的动目标的空间频率，ρi、σi分别为第i类动目标的俯仰角和方位角，d＝0.3m为阵列天线的阵元间隔，λ＝0.6m为雷达发射信号的波长。详情参见参考文献j.ward，space-timeadaptiveprocessingforairborneradar.lincolnlaboratory，mti，cambridge，ma，usa，tech.rep.1015，1994.13-17。s3、构建目标矩阵目标矩阵p表达式为：p＝[p1…pm…p11]；其中，pm为1×110列的向量，m＝1,…,11，其第ξ个元素为pm(ξ)，ξ＝1,…,110。所述pm(ξ)的值确定过程为：m＝1,…10时，目标矩阵p的第ξ行第m列元素pm(ξ)的值根据步骤s2所构建的训练数据矩阵f中第ξ行训练数据的动目标速度来确定；具体为：当该行训练数据中包含有动目标，且动目标的速度为时，pm(ξ)＝1；否则pm(ξ)＝0，为第m个多普勒类的归一化多普勒频率；m＝11时，目标矩阵p的第ξ行第m列元素pm(ξ)的值根据步骤s2所构建的训练数据矩阵f中第ξ行训练数据是否含有动目标确定；具体为：当该行训练数据不含有动目标，且vt＝0时，pm(ξ)＝1；否则pm(ξ)＝0。pm(ξ)的值确定详见参考文献alaaelkhatib，khaledassaleh,spacetimeadaptiveprocessingusingpatternclassification，ieeetransactionsonsignalprocessing,vol63,no.3pp.766-779,2015。s4、计算权系数利用由步骤s2得到的训练数据矩阵以及由步骤s3得到的目标矩阵，通过稀疏贝叶斯学习的方法，由下式计算动目标检测模型中的权系数向量wm；pm(ξ)＝hξwm；其中，wm为200×1的列向量，m＝1,…,11，hξ为1×200的行向量，ξ＝1,…,110，hξ为训练数据矩阵f的第ξ行向量，表示第ξ个训练样本，pm(ξ)是目标矩阵p的第m列的第ξ个元素；根据得到的各个wm构建权系数矩阵：w＝[w1…wm…w11]；其中，w是200行11列的矩阵。s5、检测动目标利用步骤s4得到的权系数矩阵w，以及待检测的回波数据矩阵l，计算下式q＝lw其中，l则为待检测距离门的回波数据构成的50行200列的矩阵，w是200行11列的矩阵，q为50行11列的矩阵，q中的某一个元素表示为δ是门限值，根据目标矩阵的设置，δ一般略小于1，本实施例中δ＝0.9。当矩阵q中的元素且l＝1,…10时，矩阵q中的元素则判断第个距离门含有动目标，否则判断第个距离门无动目标。如表1所示的仿真结果，比较了传统stap、基于二项式的stap、基于svr的stap以及本申请方法的输出信号干扰噪声比如表1所示的仿真结果，在使用100个辅助距离门(远小于rmb准则要求的197个)的情况下，传统stap方法无法正确检测动目标，基于二项式的stap方法输出性能较低，本论文所提出的方法比基于svr的stap方法输出性能略大，验证了本申请所提出方法有效性同时，也验证了本申请方法具有较好的输出性能。表1仿真结果通过表2的非0权值个数的统计，在使用步骤s5检测动目标时，相对于基于svr的stap方法，本发明可以少做84次乘法运算，所以本发明所提出的方法检测时间更短。表2权值统计表权值基于svr的stap本发明非0权值(个)200116本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。当前第1页12