一种无直达波条件下分布式无源雷达目标检测的方法
【专利摘要】本发明公开了一种无直达波条件下分布式无源雷达目标检测的方法，该方法属于分布式无源雷达目标检测技术领域。现有的无源雷达目标检测方法是基于经典匹配滤波理论，前提是其可以实时地获取信噪比比较高的直达波，可利用参考通道接收到的直达波信号高质量地估计出发射信号，才能获得近似最优的检测性能。为解决实际中接收不到直达波信号时的目标检测问题，本发明提出了一种无直达波条件下分布式无源雷达目标检测的方法，构造了无直达波条件下的集中式目标检测器，可用于多发多收几何架构中直达波信号信噪比比较低或根本接收不到直达波信号时开展目标检测，同时还间接避免了分布式无源雷达不同接收站间很难解决的目标匹配问题。
【专利说明】
一种无直达波条件下分布式无源雷达目标检测的方法
技术领域
[0001] 本发明属于分布式无源雷达目标检测技术领域，特别涉及一种无直达波条件下分 布式无源雷达目标检测的方法。
【背景技术】
[0002] 国内外许多实验已证明，基于第三方非合作辐射源的无源雷达具有很好的探测潜 力，能够实现对空中小目标和低空飞行目标的探测和定位。为进一步提高目标探测稳定性 和跟踪连续性，欧洲和美国许多研究机构已开展了基于FM、DAB、DVB-T、GPS信号的多发多收 体制的分布式无源雷达技术研究。
[0003] 目前，分布式无源雷达的目标检测方法是基于经典的匹配滤波理论，即利用每个 发射机-目标-接收机双基地对，无源雷达接收机将直达波和目标回波分开接收，分别对应 参考通道和目标监视通道，然后通过使天线定向指向发射机和预计目标出现的区域来获取 目标回波信号。在监视通道利用自适应滤波技术来抑制直达波和其他杂波的干扰，而在参 考通道利用均衡技术进一步分离直达波信号，然后计算参考通道和监视通道接收信号间的 互模糊函数或广义互相关。若互模糊函数图中有信号过门限则认为检测到目标，而利用对 应的目标双基地距离、双基地多普勒和到达角信息实现对目标的定位与跟踪。由于目标检 测是在每个双基地发射接收对内分别进行的，因此是分布式目标检测，检测结果在后续进 行融合处理。由于这与有源雷达采用的匹配滤波处理过程基本类似，因此这种处理方法在 实际无源雷达目标检测中常用。计算参考通道和监视通道接收信号间的互模糊函数的过程 与有源雷达中进行的匹配滤波处理近似，参考通道的信号就是对先验未知的发射信号的一 种估计。
[0004] 事实上，有源雷达中的匹配滤波是一种在奈曼-皮尔逊准则下，对发射信号准确已 知时的最优滤波器，而分布式无源雷达系统采用经典的匹配滤波理论进行目标检测的前提 是其能够实时地获取直达波，且直达波信噪比比较高，能够利用参考通道接收到的直达波 信号高质量地估计出发射信号，从而获得近似最优的检测性能。
[0005] 然而，在分布式无源雷达中，如果直达波信号信噪比比较低，无法高质量地估计出 发射信号，则导致与真实的发射信号失配，检测性能将下降，甚至无法检测到目标。虽然通 过通道均衡处理在某种程度上减轻失配影响，但不可能完全消除失配。此外，在实际中由于 复杂地形环境的影响，如果在非合作辐射源和接收机间的视距传播路径中还存在物理遮 挡，无源雷达将接收不到来自非合作辐射源的直达波信号，采用经典的匹配滤波理论无法 完成目标检测。因此，通过计算参考通道和监视通道接收信号间的互模糊函数或广义互相 关的目标检测方法不适合于直达波信噪比很低或根本接收不到直达波信号的情况。
【发明内容】
[0006] 本发明的目的是提供一种无直达波条件下分布式无源雷达目标检测的方法，能够 用于解决多发多收几何架构下，分布式无源雷达接收不到直达波信号情况下的目标检测问 题，其中要解决的技术问题包括：
[0007] (1)给出多发多收几何架构下，建立无直达波条件下分布式无源雷达目标回波信 号分量的实施步骤；
[0008] (2)在多发多收几何架构下，构建无直达波条件下分布式无源雷达集中式目标检 测统计量和目标检测的实施步骤。
[0009] 本发明所述的一种无直达波条件下分布式无源雷达目标检测的方法，包括以下步 骤：
[0010] (1)建立多发多收几何架构下目标回波信号分量的实施步骤，具体包括如下子步 骤：
[0011] A1.构建目标运动时，分布式无源雷达系统中第i个辐射源发射的信号经目标反射 后被第j个接收机阵列的第n阵元接收到并经基带处理后的目标回波信号分量
[0012] A2.对第ij个双基地对中，第j个接收机的第n个阵元接收到的由目标回波和接收 机噪声组成的信号进行量化采样，采用时延多普勒算子表示其离散形式，并给出波束 形成后的目标回波信号
[0013] A3.利用分布式无源雷达系统中Nr个接收机阵列接收到的与Nt个非合作雷达辐射 源对应的所有目标回波信号采样，构造回波信号矩阵s;
[0014] (2)构建无直达波条件下分布式无源雷达集中式目标检测统计量，开展目标检测， 具体包括如下子步骤；
[0015] B1.利用回波信号矩阵s构建二元备择假设检验；
[0016] B2.利用已构建的假设检验，推导其广义对数似然比，得到无直达波条件下分布式 无源雷达目标检测的集中式检测统计量Iss ;
[0017] B3.根据分布式无源雷达系统中所有目标回波信号的实际采样，计算出目标检测 统计量U，通过比较检测统计量U与门限k的大小，然后判定目标是否存在，完成目标检 测 。
[0018] 优选的，步骤(2)中构造无直达波条件下分布式无源雷达目标检测统计量U寸，弓丨 入目标位置矢量和速度矢量，作为开展二元备择假设检验时的目标检测单元；
[0019]优选的，步骤(2)中无直达波条件包括根本接收不到直达波信号的情况和直达波 信噪比非常低无法估计出非合作辐射源发射信号波形的情况。
[0020] 与目前分布式无源雷达中基于经典匹配滤波理论的分布式目标检测方法相比，本 发明提供的分布式无源雷达目标检测的方法能够在无直达波或者直达波信噪比非常低的 情况下实现目标检测。此外，在无直达波条件下构造集中式目标检测器时，引入了目标位置 矢量和速度矢量，使得在实现目标检测的同时，间接解决了目标的定位问题，避免了多发多 收几何架构中存在的不同接收机阵列间很难解决的定位模糊问题和目标匹配问题。
【附图说明】
[0021] 附图1是本发明的无直达波条件下分布式无源雷达的拓扑结构示意图。
[0022] 附图2是本发明的分布式无源雷达第ij对非合作辐射源-目标-接收机阵列的几何 关系图。
[0023] 附图3是本发明的第ij个双基地对对应目标回波信号的波束形成示意图。
[0024] 附图4是本发明的分布式无源雷达目标检测方法实施流程图。
[0025] 附图5是本发明实施例的计算机仿真结果示意图。
【具体实施方式】
[0026] 以下结合说明书附图对本发明作进一步详细描述。
[0027]如图1所示，分布式无源雷达系统包括Nt个发射机，发射机在无源雷达领域中也称 为非合作辐射源，Nr个接收机阵列，1个目标，其中Nt彡2，Nr>2。
[0028] 如图2所示，第ij个双基地对的几何关系和信号环境，第i个发射机的位置和速度 分别记为d1和七，i = 1，…，N t，第j个接收机阵列的位置和速度分别记为r 和F，j = 1，…， Nr，而目标的位置和速度分别记为t和i，其中^、^、^、^都是时间的函数广般情况 下，发射机和接收机、目标都是运动的。第i个发射机到第j个接收机的距离为 切〇 = |卜;|。类似地，< (/) = |t -心||和祀⑴=|r^卜 离和目标到第j个接收机的距离。第j个接收天线是有W个阵元的阵列，」V/k2，l<j$Nr， 第n个阵元的位置为匕=r/+W，其中是参考阵元的位置，S,(是第n个阵元 相对参考阵元的指向偏移矢量，而5丨=0 ;接收天线阵列的阵元数都相同，g卩=；Vp，j = 1，…，Nr。第j个接收机的第n个阵元到位置x的单位指向矢量为kh；)，即
。在 远场，对于给定的X，，即从阵列阵元到远场某一位置的单位指向矢量 近似相等。
[0029] 如图3所示，多发多收几何架构下分布式无源雷达系统中，所有无源接收机均采用 阵列天线，通过波束形成的方法形成目标监视通道，从而实现目标回波信号的接收。
[0030] 第i个辐射源发射信号表示为
[0031] ，/'(/) = "'(/),'./e[o,r] (1)
[0032] 其中，表示载频，T表示信号持续时间，表示第i个发射机对应的信号，uKt) 为对应的复包络，频域表示为1^( ? )，带宽为，且当| ? | SjtB1时，UU ? )~0， {/?(/):/_ = 1.…在频域没有重叠。
[0033] 信号沿着目标路径通道传播到第j个接收机，第j个接收机阵列的第n个阵元接收 到的信号g 是来自所有在接收机带内目标回波与接收机噪声之和，即
[0034] .V,； (/) = (/)?'(? - r；；" (〇) + ?! (〇 (2) /=1
[0035] 其中，< ⑷为目标路径通道的幅度系数，P为与第i j个双基地对对应目标的复 双基地反射系数，⑷为目标路径通道的传播时延，圮⑷是功率谱密度为 疔（《) = 乂，|?- < | 的广义平稳高斯白噪声，带宽为Bj，载频为<。通道系数⑷ 考虑了发射、传播以及目标路径通道的影响，定义为
，其中^ 为第i个发射机指向X的有效辐射功率，f =r/./；:为第i个发射机发射信号的波长，C为光速， /2冗。在[0，T ]内，被⑷、i?;'⑷和祀(〇都不会有显著变化，则af全< ⑴匕。信号 <(〖；)经下变频和频域信道化处理后，提取每个发射信号的复基带信号。记第i个通道的复基 带信号为 <，利用式(1)和(2)，第j个接收机的第n个阵元接收到的信号表示为
[0036] s；； (?) = ?,/a；,e，(e，"〇；r；；"(，t)M , (/ - r；；" (〇) + ?;； (?) (3)
[0037] 其中，03为第j个接收机下变频处理时本振的未知相位，表明接收机阵列间是非相 参的。
[0038] 如图4所示，本发明提供一种无直达波条件下实现分布式无源雷达目标检测方法 的实施流程图，【具体实施方式】包括以下子步骤：
[0039] A1.构建多发多收几何架构下目标运动时，分布式无源雷达系统中第i个辐射源发 射的信号经目标反射后被第j个接收机阵列的第n阵元接收到并经基带处理后的目标回波 信号分量#,&；)，具体过程如下：
[0040] 4 (/) = (t - rf,, (?)) (4)
[0041] 进一步表示为
(5)
[0043]式中（a)为幅度尺度因子，（b)为未知本振相位，（c)为第n个阵元接收目标回波信 号相对参考阵元的相位差，（d)参考载频相位因子，（e)时延后复基带信号，（f)多普勒频率 因子，而 < 为从第i个发射机到目标到第j个接收机的双基地时间延迟，即
(6)
[0045] 劣〇：}为第n个阵元接收目标回波信号相对参考阵元的相位差，即
[0046] 巧(，)兰-?'(〇 二 (7)：
[0047] 为目标的双基地多普勒频移
(8)
[0049] 因此，第i个辐射源发射的信号经目标反射后被第j个接收机阵列的第n阵元接收 到，并经基带处理后的目标回波信号分量构建为
[0050] sl；n(t) = r：，e'3：Ut}u,(t-T,/)e ，H'1 (9)
[0051 ] 其中，rf是第ij个双基地对对应的目标路径通道系数：
[0052] A2.对第ij个双基地对中，第j个接收机的第n个阵元接收到的由目标回波和接收 机噪声的和组成的信号进行量化采样，采用时延多普勒算子表示其离散形式，并给出 波束形成后的目标回波信号<，具体过程如下：
[0053] 以/) =1/7； =F的采样频率进行量化采样，则二彳，得离散信号形式为
[0054] < [/]二（,;[/] + ":〈[/],/ = 0,…，1 (10)
[0055] 其中，￡ =^/)」为总采样点数，目标回波感的量化形式为
[0056] [/] = [l - ei ]eJV?，J = L-\ (11)
[0057] 其中，.vf =r?f///分别为每个样本的归一化多普勒频率，单位为弧度， 《卜4乂1 = </」，分别为每个样本的归一化时延。记[*〖1=€[/]，[€>4[/1，1=0，-丄 1-1。噪声采样序列<[丨]~以(0,(12)，〇 2 =關1为平均噪声功率，五|<卞](<[/])1=〇^_,^,1，^^ 为Kronecker符号。i € C￡'xl 为发射波形，第1个元素为卜'1, = 〃'（0 ?= 〃' J =()，…，￡ -1。
[0058] 定义Df? C1气Dl(x) =diag( [ej(°U1)x，…，ej(L-1)x])，其中diag(x)中对角 线元素xeCf，是LXL的方阵，因此[diag(x)] n,n=[x]n。最后，令WieC^为酉离散傅里叶 变换矩阵，其第(m，n)个元素为
(12)
[0060]其中，m = 0,…，L_1，n = 0,…，L-1 ?，则
[0061] C = rfem)Dr (tf )(wfDf (-2^f/Z/)Wr )uJ (13)
[0062] 引入时延多普勒算子，定义为= D￡,㈦，由于 +{七），M (Q今=!)(-4-U)，因此时延多普勒算子及(4# 为酉算子，即Z^(〇)Z)(〇)=Ifi，其中&是1^1/的单位矩阵。
[0063] 因此，目标回波信号分量的离散形式为
[0064] 4? =rfej3l!it)D((Iluf)ui (14)
[0065] 第j个接收机阵列第n阵元接收到信号的离散形式为
[0066] 别)D( 0)u;. + nf (15)
[0067] 其中，< ~ GV(〇z, ?夂\是长度为L1的零矢量。通过波束形成而获得的目标回 波信号<为
[0068] sf =幺[bf]:sf (16): n=l
[0069] bf为目标监视通道的波束形成器。
[0070] A3.利用分布式无源雷达系统中Nr个接收机阵列接收到的与Nt个非合作雷达辐射 源对应的所有目标回波信号采样，构造回波信号矩阵s，具体过程如下：
[0071] 记af =a"'(t)eCA〃xl表示在t方向的空间指向矢量a"⑴=Z⑷Z⑴乂及⑴ 7， 则第j个接收机阵列所有个阵元接收到的与第i个辐射源对应的信号矢量 sg =⑷1'…⑷7 r化啊为
[0072] f .j = \,.'Nnj = l.、Nr (17)
[0073] 其中，，nf - 0f,)，方差为〇2，矩阵Mf = <(a:f ?巧卜 _?表示Kronecker积，时延多普勒算子为W=Dt",'，p,v)
[0074] 因此，所有Nr个接收机阵列收到的与第i个辐射源对应的采样s 1为 s;= (sn)r,---,(sw')r e€5z><1，而与所有Nt个非合作辐射源和Nr个接收机阵列对应的所有采 样组成的矩阵为
[0075] s= (8丫，."，(8~广 reCWxl (18)
[0076] 即s为与所有Nt个非合作辐射源对应的s1组成的矩阵。
[0077] B1.引入待检测目标的位置矢量和速度矢量，作为目标检测单元，利用回波信号矩 阵s构建二元备择假设检验，具体过程如下：
[0078] 令待检测目标的位置和速度对应在(P，P)单元，即检测单元，其中p，|j分别表示目 标的位置和速度。构建二元备择假设检验，BP
[0079] A : S9 = My + n"; K : s" = n3 (19)
[0080] 其中，1 = 1，"，队，」=1，"，心，1^=<(3〇1)；：)￡15￡'3卜￥?(1： ；￥，1表 示在P方向的空间指向矢量= e胃(R)…而#为目标位置为P时对应目 标路径通道的系数，/)|为目标状态为对应的时延多普勒算子。
[0081 ] B2.利用已构建的假设检验，推导其广义对数似然比，得到分布式无源雷达目标检 测的集中式检测统计量U，具体过程如下；
[0082] 由于接收机噪声与发射机通道无关，所以在邱假设下的条件概率密度Pl(s| YP， u)为 Nr
[0083] A (su) = npi (s' \fP^u') (20)
[0084] 其中，
，而在和假设下的条件概 率密度PQ(S)为
(」）
[0086]其中，发射信号U和通道系数yP是确定性的未知参量。因此，M是复合假设，因为 pi(s| yP，u)是以发射信号U和通道系数参数的。因此，将在似然比检验中的未知参数 用其最大似然估计替换即可得到广义似然比。令li( y p，u| s) = logpi(s I y p，U)，lo(s) = log P〇(s)，则广义对数似然比为
(22)
[0088]其中，k由分布式无源雷达系统的虚警概率决定。分别推导得到|s)和l〇 (s)的解析式后，得到集中式目标检测统计量|%为
(2^)
[0090] 其中，M ?)为矩阵参量的最大特征值，eCA^为Gram矩阵，（?）H 表示厄米特转置，而奶=[^，…,0>(0￡>^ :，§丨是时延多普勒补偿后的目标回波信号sf， 即时延和多普勒频移补偿后的监视通道目标回波信号。由于矩阵的元素 = 是由监视通道的目标回波信号的内积构成，所以|ss能够用于接收不到直 达波或者直达波参考信噪比非常低的情况下的目标检测。
[0091] B3.根据实际分布式无源雷达系统中所有目标回波信号的采样，计算出目标检测 统计量U，通过比较检测统计量U与门限k的大小，然后判定目标是否存在，完成目标检 测，具体过程如下：
[0092]利用分布式无源雷达系统中Nr个接收机阵列对应Nt个非合作雷达辐射源的所有采 样构造回波信号矩阵，计算目标检测统计量Iss，当Iss多K时，则判断目标存在，且目标位置 和速度状态为而当lss<K时，则判断目标不存在。
[0093]如图5所示，本发明实施例的计算机仿真结果示意图，图中BRng等值线表示双基地 距离等值线，BDop等值线表示双基地多普勒等值线，目标出现在真实目标位置附近，同时还 实现了对目标的定位。在仿真环境背景中，两个发射机的位置分别为diitO.SJlkm和d 2 = [-0.5，-4]1〇11，三个接收机的位置分别为1'1 = [-4,2]，1'2 = [-4,0.5]和1'3 = [-4，-2.5]1〇11，目 标在七=[4,0]1〇11，目标运动速度为丨=[/,1，/1.] = [_375.375]111/8。发射机信号载频分别为 8.0GHz和8.1GHz，全向辐射功率均为= 50 W 均匀线性阵列由6个阵元组成，所有接收机 天线都是阵元间隔为1.875cm的均匀线性阵列，指向+px方向，各阵元的波瓣图为G/G；) = l。 复基带信号<采样率乜=5001^，积累时间为1 = 21]18,〇^=2、0019\10-14'^，111 = 61口{_]_91}， 0; eRixl是相互独立的随机相位矢量，在[0,2jt]间服从均勾分布，L = fsT=1000,目标截面 积为1 OdBsm。接收机间没有相位同步，目标回波的平均信噪比为SNRavg = -15dB。
【主权项】
1. 一种无直达波条件下分布式无源雷达目标检测的方法，其特征在于，包括以下步骤： (1) 建立多发多收几何架构下目标回波信号分量的实施步骤，具体包括如下子步骤： Al.构建目标运动时，分布式无源雷达系统中第i个辐射源发射的信号经目标反射后被 第j个接收机阵列的第η阵元接收到并经基带处理后的目标回波信号分量 Α2.对第ij个双基地对中，第j个接收机的第η个阵元接收到的由目标回波和接收机噪 声组成的信号私?进行量化采样，采用时延多普勒算子表示其离散形式，并给出波束形成 后的目标回波信号sf ; A3.利用分布式无源雷达系统中Nr个接收机阵列接收到的与Nt个非合作雷达辐射源对 应的所有目标回波信号采样，构造回波信号矩阵s; (2) 构建无直达波条件下分布式无源雷达集中式目标检测统计量并开展目标检测的实 施步骤，具体包括如下子步骤： Bl.利用目标回波信号矩阵s构建二元备择假设检验； Β2.利用已构建的假设检验，推导其广义对数似然比，构造无直达波条件下分布式无源 雷达目标检测统计量Iss ; Β3.根据分布式无源雷达系统中所有目标回波信号的实际采样，计算出目标检测统计 量Iss，通过比较检测统计量Iss与门限IC的大小，然后判定目标是否存在，完成目标检测。2. 如权利要求1所述的一种无直达波条件下分布式无源雷达目标检测的方法，其特征 在于，所述步骤(2)中构造无直达波条件下分布式无源雷达目标检测统计量I ss时，引入目标 位置矢量和速度矢量，作为开展二元备择假设检验时的目标检测单元。3. 如权利要求1所述的一种无直达波条件下分布式无源雷达目标检测的方法，其特征 在于，所述步骤(2)中无直达波条件包括根本接收不到直达波信号的情况和直达波信噪比 非常低无法估计出非合作辐射源发射信号波形的情况。
【文档编号】G01S13/56GK105891817SQ201610410562
【公开日】2016年8月24日
【申请日】2016年6月8日
【发明人】张财生, 唐小明, 朱洪伟, 刘明春, 李国君, 张涛, 宋洪良
【申请人】中国人民解放军海军航空工程学院