一种多站雷达异步多帧联合检测方法技术

本发明专利技术公开了一种多站雷达异步多帧联合检测方法，属于雷达目标检测领域。在每个时刻，本发明专利技术利用雷达组网中多个雷达站点提供的多帧异步回波数据计算联合广义对数似然比，然后通过求解最优化问题找到一组使联合似然比最大的状态序列作为目标航迹的估计。另外，本发明专利技术结合了维特比算法，能够快速高效地实现该最优化求解问题，在保证算法检测实时性的同时，提高系统微弱目标的检测概率，改善雷达系统的探测性能。

【技术实现步骤摘要】

本专利技术属于雷达目标检测领域。
技术介绍
单基地雷达系统只利用单部雷达对监视区域进行扫描，获得的采样数据单一，在强杂波、目标起伏等情况下，对目标的检测概率低。而多站协同探测联合处理多个雷达节点提供的信息，由于利用了更高的数据维度，该系统能够捕获目标RCS多样性、频移多样性，进而改善雷达系统的探测性能，提供更可靠的检测结果。但随着电子技术的不断发展，越来越多的隐身目标逐渐出现。这些目标的信噪比低且回波弱，导致协同探测系统的探测性能和威力范围大大下降。利用单帧数据检测的协同探测方法已经不能满足系统对这类目标的探测需求。另外，由于雷达组网中不同雷达节点的起始扫描时刻、扫描周期、扫描范围等都不相同，融合中心接收到的是一系列异步数据，传统的多帧检测方法并不适用。因此，寻找一种能够适用于多站雷达协同探测的多帧数据联合处理方法已经成为一个亟待解决的问题。
技术实现思路
本专利技术所要解决的技术问题是，提供一种适用于多站雷达协同探测的多帧数据联合处理方法。本专利技术为解决上述技术问题所采用的技术方案是，一种多站雷达异步多帧联合检测方法，包括以下步骤：步骤1:初始化系统参数：初始化系统参数，主要包括：雷达节点数M，雷达m的发射信号sm(t)为简单脉冲串，载波频率fc，回波信号采样间隔Ts，采样次数N，各雷达节点向融合中心传输的数据帧数K，系统检测门限VT；初始化系统变量：i＝1；步骤2：将各雷达节点在K个扫描时刻接收到的回波信号传送到融合中心；步骤3：时间校准：在融合中心，将所有雷达节点不同时刻的回波信号按时间顺序重新排序；重新排序后的信号序列表示为[r1(t),…,rMK(t)]，并记录重新排序后第i(i＝1,…,MK)个信号的所属雷达标号为βi；步骤4：计算对数似然比：步骤4.1：计算第i次扫描时，量化状态对应的回波时延τ和多普勒频移fd，即：其中，(xi,yi)代表目标位置，代表目标速度，(xm,ym)为第βi个雷达的坐标，为分辨单元(xi,yi)与雷达βi的夹角，c表示光速；步骤4.2：如果第i次扫描，分辨单元(xi,yi)在雷达βi的扫描范围内，则确定目标状态为θi时的反射信号d(t,θi)，执行步骤4.3；否则令状态θi的对数似然比等于0，执行步骤4.5；步骤4.3：对信号采样：按照采样间隔Ts对反射信号d(t,θi)和回波信号ri(t)，在[τ,…,τ+(N-1)Ts]处进行N次采样，得到采样后信号d(θi)和ri；其中d(θi)＝[d(τ,θi),…,d(τ+(N-1)Ts,θi)]，ri＝[ri(τ),…,ri(τ+(N-1)Ts)]；步骤4.4：计算目标状态为θi时的对数似然比，即：其中，C为噪声协方差矩阵；步骤4.5：采用步骤4.1到步骤4.4的方法计算所有雷达节点不同时刻的回波信号的对数似然比；步骤5：计算值函数：从第一帧回拨信号开始，i＝1，用各个量化状态对应的对数似然比初始化各状态的值函数，即：I(θ1)＝lnΛ(r1|θ1)；如果i>1，首先确定i时刻与i-1时刻的扫描间隔，然后根据这个扫描间隔确定i-1时刻能够转移到θi的有效状态转移范围δ(θi)，并利用该范围内状态的值函数最大值更新θi的值函数，即：其中，I(θi)表示状态θi的值函数；步骤6：如果i≠MK，则令i＝i+1，返回步骤5；否则，执行步骤7；步骤7：如果值函数的最大值超过检测门限，即则宣布检测到目标，否则宣布没有检测到目标。进一步的，所述步骤7中检测门限的确定方法为:设定虚警概率值为pfa，联合处理l次扫描数据的门限的仿真生成方式如下：1)生成l帧噪声数据；2)利用l帧数据计算状态空间中每个状态的值函数，并记录值函数最大值；3)重复执行步骤1)，2)N次，将N次的值函数最大值按降序排列存入集合F中；4)取F中第Npfa个元素，作为门限。本专利技术的优点是能够联合处理多个雷达节点提供的多帧异步回波数据，提供更可靠的检测结果，同时改善系统对微弱目标的检测性能。另外，本专利技术结合了维特比算法，可以快速高效的实现算法中涉及的最优化求解问题，保证系统检测的实时性。附图说明图1为本专利技术的流程图。图2为单个扫描时刻的检测统计量平面与本专利技术计算的检测统计量平面对比图。(a)为单个扫描时刻的检测统计量平面；(b)为本专利技术计算的检测统计量平面。图3为本专利技术与单帧检测的多站雷达系统的检测概率对比图。具体实施方式本专利技术主要采用计算机仿真的方法进行验证，所有步骤、结论都在MATLAB-R2012b上验证确认。具体实施步骤如下：步骤1)初始化系统参数：在本实例中：雷达节点数M＝3，雷达发射波形为简单脉冲串，即载波频率fc＝109，回波信号采样间隔等于脉冲重复频率，即Ts＝Tp＝10-3，采样数N＝10，每个雷达节点向融合中心传输的数据帧数K＝3，对应的检测门限VT由蒙特卡洛仿真得到为82.92。初始化系统变量：i＝1。步骤2)将各雷达节点在K个扫描时刻接收到的回波信号传送到融合中心。步骤3)时间校准：在融合中心，将不同雷达节点的不同扫描时刻的回波信号按时间顺序重新排序。重新排序后的信号序列表示为[r1(t),…,rMK(t)]，并记录重新排序后第i(i＝1,……,MK)个信号的所属雷达标号为βi。步骤4)计算对数似然比：4.1)计算第i次扫描时，量化状态对应的回波时延τ和多普勒频移fd，即：其中，(xi,yi)代表目标位置，代表目标速度，为第βi个雷达的坐标，为分辨单元(xi,yi)与雷达βi的夹角，c表示光速。4.2)如果第i次扫描，分辨单元(xi,yi)在雷达βi的扫描范围内，则确定目标状态为θi时的反射信号d(t,θi)，执行步骤4.3)；否则令状态θi的对数似然比等于0，执行步骤4.5)。4.3)对信号采样：按照采样间隔Ts对反射信号d(t,θi)和回波信号ri(t)，在[τ,…,τ+(N-1)Ts]处进行N次采样。得到采样后信号d(θi)和ri。其中d(θi)＝[d(τ,θi),…,d(τ+(N-1)Ts,θi)]，ri＝[ri(τ),…,ri(τ+(N-1)Ts)]。4.4)计算目标状态为θi时的对数似然比，即：其中，C为噪声协方差矩阵。4.5)如果所有量化状态θi的对数似然比都已计算完，执行步骤4.6)；否则，返回步骤4.1)计算下一个量化状态θi的对数似然比。4.6)如果i<MK，则令i＝i+1，返回步骤4.1)；否则，令i＝0，i＝i+1，执行步骤5)。步骤5)计算值函数：如果i＝1，用各个量化状态对应的对数似然比初始化各状态的值函数，即：I(θ1)＝lnΛ(r1|θ1)。如果i>1，首先确定i时刻与i-1时刻的扫描间隔，然后根据这个扫描间隔确定i-1时刻能够转移到θi的有效状态转移范围δ(θi)，并利用该范围内状态的值函数最大值更新θi的值函数，即：其中，I(θi)表示状态θi的值函数。步骤6)如果i≠MK，则令i＝i+1，返回步骤5)；否则，执行步骤7)。步骤7)如果值函数的最大值超过检测门限，即则宣布检测到目标，否则宣布没有检测到目标。图2为单个扫描时刻的检测统计量平面与本专利技术计算的检测统计量平面对比图。从图2中可以看出，利用单个扫描时刻的数据，目标状态的检测统计量淹没在噪声中，但利用本专利技术联合处理多个雷达节点提本文档来自技高网...

【技术保护点】
一种多站雷达异步多帧联合检测方法，包括以下步骤：步骤1:初始化系统参数：初始化系统参数，主要包括：雷达节点数M，雷达m的发射信号sm(t)为简单脉冲串，载波频率fc，回波信号采样间隔Ts，采样次数N，各雷达节点向融合中心传输的数据帧数K，系统检测门限VT；初始化系统变量：i＝1；步骤2：将各雷达节点在K个扫描时刻接收到的回波信号传送到融合中心；步骤3：时间校准：在融合中心，将所有雷达节点不同时刻的回波信号按时间顺序重新排序；重新排序后的信号序列表示为[r1(t),…,rMK(t)]，并记录重新排序后第i(i＝1,…,MK)个信号的所属雷达标号为βi；步骤4：计算对数似然比：步骤4.1：计算第i次扫描时，量化状态对应的回波时延τ和多普勒频移fd，即：其中，(xi,yi)代表目标位置，代表目标速度，(xm,ym)为第βi个雷达的坐标，为分辨单元(xi,yi)与雷达βi的夹角，c表示光速；步骤4.2：如果第i次扫描，分辨单元(xi,yi)在雷达βi的扫描范围内，则确定目标状态为θi时的反射信号d(t,θi)，执行步骤4.3；否则令状态θi的对数似然比等于0，执行步骤4.5；步骤4.3：对信号采样：按照采样间隔Ts对反射信号d(t,θi)和回波信号ri(t)，在[τ,…,τ+(N‑1)Ts]处进行N次采样，得到采样后信号d(θi)和ri；其中d(θi)＝[d(τ,θi),…,d(τ+(N‑1)Ts,θi)]，ri＝[ri(τ),…,ri(τ+(N‑1)Ts)]；步骤4.4：计算目标状态为θi时的对数似然比，即：其中，C为噪声协方差矩阵；步骤4.5：采用步骤4.1到步骤4.4的方法计算所有雷达节点不同时刻的回波信号的对数似然比；步骤5：计算值函数：从第一帧回拨信号开始，i＝1，用各个量化状态对应的对数似然比初始化各状态的值函数，即：I(θ1)＝lnΛ(r1|θ1)；如果i>1，首先确定i时刻与i‑1时刻的扫描间隔，然后根据这个扫描间隔确定i‑1时刻能够转移到θi的有效状态转移范围δ(θi)，并利用该范围内状态的值函数最大值更新θi的值函数，即：其中，I(θi)表示状态θi的值函数；步骤6：如果i≠MK，则令i＝i+1，返回步骤5；否则，执行步骤7；步骤7：如果值函数的最大值超过检测门限，即则宣布检测到目标，否则宣布没有检测到目标。...

【技术特征摘要】
1.一种多站雷达异步多帧联合检测方法，包括以下步骤：步骤1:初始化系统参数：初始化系统参数，主要包括：雷达节点数M，雷达m的发射信号sm(t)为简单脉冲串，载波频率fc，回波信号采样间隔Ts，采样次数N，各雷达节点向融合中心传输的数据帧数K，系统检测门限VT；初始化系统变量：i＝1；步骤2：将各雷达节点在K个扫描时刻接收到的回波信号传送到融合中心；步骤3：时间校准：在融合中心，将所有雷达节点不同时刻的回波信号按时间顺序重新排序；重新排序后的信号序列表示为[r1(t),…,rMK(t)]，并记录重新排序后第i(i＝1,…,MK)个信号的所属雷达标号为βi；步骤4：计算对数似然比：步骤4.1：计算第i次扫描时，量化状态对应的回波时延τ和多普勒频移fd，即：其中，(xi,yi)代表目标位置，代表目标速度，(xm,ym)为第βi个雷达的坐标，为分辨单元(xi,yi)与雷达βi的夹角，c表示光速；步骤4.2：如果第i次扫描，分辨单元(xi,yi)在雷达βi的扫描范围内，则确定目标状态为θi时的反射信号d(t,θi)，执行步骤4.3；否则令状态θi的对数似然比等于0，执行步骤4.5；步骤4.3：对信号采样：按照采样间隔Ts对反射信号d(t,θi)和回波信号ri(t)，在[τ,…,τ+(N-1)Ts]处进行N次采样，得到采样后信号d(θi)和ri；其中d(θi)＝...

【专利技术属性】
技术研发人员：易伟，王经鹤，李溯琪，姜萌，方梓成，孔令讲，崔国龙，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51