简介:这是发表在《电子与信息学报》2017年第1期的学术论文PDF,面向雷达信号处理、电子对抗相关专业的研究者与工程师。内容围绕分布式多站雷达体制下转发式欺骗干扰的有效性展开,从空间分辨单元物理含义出发,推导了多站雷达在特定空间分辨率下可被有效干扰的条件,并区分干扰机位于远场与近场两种情形,给出最大欺骗距离范围与控制转发延迟量的具体结论。论文同时为干扰策略优化和雷达抗干扰算法改进提供了理论支撑,适合作为科研选题、论文写作或课程设计的参考文献。压缩包内共1个文件,为PDF格式,包含中英文摘要、公式推导及详细参考文献,大小仅414KB,便于阅读与打印。该资源已有134人学习下载,对正在研究分布式雷达干扰与抗干扰技术的读者具有直接参考价值。
1. 分布式多站雷达转发式欺骗干扰:一块看似铁板的区域,实际有缝
分布式多站雷达转发式欺骗干扰,这个组合词在电子对抗圈里经常被默认成“不好使”——多站雷达做空间配准时,会把不在同一个空间分辨单元里的假目标自动剔除,转发式干扰辛辛苦苦造出来的距离假目标瞬间变废。但这篇2017年的期刊论文,恰恰在这块铁板上敲出一条缝:只要干扰机和各接收站的几何关系满足特定条件,转发式干扰不仅能骗,而且在远场条件下任意转发时延都能骗成功。它把“骗不骗得了”这个问题,从功率对抗拉回到纯粹的几何层面。适合做多站雷达抗干扰仿真的人、写DRFM干扰机参数表的雷达工程师,以及刚入门电子对抗、想找一套完整推导当参考的从业者。
2. 空间分辨单元与多站雷达工作流程:假目标被“空间配准”剔除的几何原因
2.1 从发射站到接收站:正交波形与MN个信道
分布式多站雷达的典型配置是M个发射站、N个接收站,各发射站发射相互正交的信号。对第m个发射站的发射信号来说,它与其他发射信号的内积在时间对齐后趋向于零,这样每个接收站收到混合回波后,只需用一组由各发射信号构造的匹配滤波器做分离,就能得到M×N个发射-接收通道对,论文里称之为信道。每个信道本质上是一部双基雷达:发射站m负责照射,接收站n负责接收。
这个结构带来的直接好处是探测视角增加:同一个目标同时被多个接收站从不同角度观察,空间分集让检测概率和参数估计精度都优于单站,这也是分布式多站雷达抗干扰能力的来源。但代价也随之而来——所有接收站的数据必须拼到同一个空间坐标系里才有意义,这个环节就是空间配准。
常见的做法是:各接收站先完成脉冲压缩和匹配滤波,把时域回波变换成距离-方位图,再送到融合中心。融合中心拿到的是多个“视角下的距离像”,而不是一个统一的点迹列表。一个容易忽略的工程点是:信道分离的前提是正交波形设计,这个条件一旦被破坏,比如多普勒频移把正交性拉歪,通道间会产生串扰,空间配准的第一步就埋下隐患。我在做仿真时,通常会先把各波形的互相关系数打印出来看一眼,峰值旁瓣差个三五分贝,后面判别假目标时就会处处被动。
为了讲清楚后面的几何推导,可以把信道和雷达体制对应起来,方便核对参数传递关系:
| 信道编号 | 发射站 | 接收站 | 等效体制 |
|---|---|---|---|
| 信道m-n | 发射站m | 接收站n | 双基雷达 |
| 信道1-1 | 发射站1 | 接收站1 | 双基雷达1 |
| 信道1-2 | 发射站1 | 接收站2 | 双基雷达2 |
2.2 SRC的物理含义:扁长分辨单元决定了哪些假目标能留下来
空间分辨单元(Space Resolution Cell,SRC)是整篇论文的核心,定义为各信道距离分辨单元的重叠区域。以1个发射站、2个接收站组成的多站雷达为例,两个信道各对应一部双基雷达,每部双基雷达的距离分辨单元是相邻两条等距离环之间的带状区域,两个带状区域交叉出来的那块阴影,就是该多站雷达的SRC。
从几何上看,SRC是扁长的:两条等距离环的切线方向决定了沿法线方向的分辨能力,交叉区域的径向尺寸通常小于横向尺寸,也就是说多站雷达的径向分辨率优于横向分辨率,而且两个接收站相对目标的空间探测夹角越小,横向分辨率越差,SRC被拉得更长。这一点是后面近场、远场结论的物理基础。可以理解成:多站雷达对“同一目标”的认定是一个空间格子,假目标只要能落进这个格子,空间配准阶段就不会被优先剔除。
判别两个假目标是否落在同一个SRC里,论文给了一组非常朴素的不等式:目标A和目标B到接收站1的距离差绝对值不超过该站距离分辨率,同时到接收站2的距离差绝对值也不超过该站距离分辨率。两个条件同时成立,才算有效欺骗。这个判据看起来简单,实际用起来有个隐含前提:距离分辨率的定义要一致。两站如果带宽不同,分辨率r1、r2数值不同,判据也要分别写,不能拿一个值代两个条件。
距离分辨率的工程估算一般是c/(2B),B是接收站信号带宽。真实系统里还要考虑加窗后的主瓣展宽,通常在理论值上乘一个1.2~1.5的系数。论文的理论分析用的是理想矩形谱,仿真时建议按实际窗函数修正后再代入不等式,否则边界上容易误判。
2.3 空间配准的四步处理流程:网格划分、信号向量组装、联合检测与参数估计
多站雷达的联合处理流程可以拆成四步,这四步在代码里对应四个功能模块。第一步,按SRC把各信道覆盖的空间区域划分成网格。网格大小直接决定运算量和配准精度,网格划得比SRC小对性能没有额外收益,反而成倍增加计算量;网格划得太大,两个不同目标会被揉进同一个单元,空间分辨率优势就废了。我一般取网格边长等于SRC尺寸,必要时在目标区域加密到SRC的一半。
第二步,对每一个SRC网格,把所有信道在该单元的信号幅度组装成回波信号向量r = [r11, r12, ..., rMN],这个向量是后续联合检测的输入。组装顺序要和信道编号一致,否则融合中心拿到的向量跟发射-接收关系对不上,定位解算时会出现系统性偏差。实际工程中,这一步往往还要做幅度归一化,因为不同接收站的噪声基底、天线增益不一样,直接拼原始幅度会让强站主导整个检测统计量。
第三步,利用该向量做联合检测,判断这个SRC里有没有目标。论文提到常用非相参积累检测器,理由是回波相位关系复杂,双基几何下不同信道的相位差未必符合单站模型,非相参积累实现简单、对相位误差不敏感。第四步,若检测到目标,用联合参数估计方法做精确定位,甚至可以估算目标速度矢量。
我踩过的一个坑是:一开始把“时间对齐”和“空间配准”看成一步。实际上时间对齐只需要按统一时间基准搬齐脉冲,空间配准是按SRC网格对齐空间位置,两者顺序反了,回波信号向量组装出来就是错位的,后续检测性能掉得莫名其妙。先把时间上同一时刻的脉冲拿出来,再做SRC网格映射,这是标准做法。
注意:假目标被剔除的过程发生在第二步到第三步之间。多站雷达不会人为判断“这个目标是不是假的”,它只是按SRC对信号幅度做联合积累,目标分散到不同单元后,每个单元里的幅度都不足以过检测门限,自然就丢了。理解这一点,后面推导欺骗条件时思路会顺很多。
3. 转发式干扰模型与远场结论:假目标点位怎么算,为什么远场“任意时延”必骗
3.1 DRFM存储转发:从延迟τ到欺骗距离Δd的换算
转发式干扰一般用数字射频存储器(DRFM)实现。DRFM把截获的雷达信号快速数字化存储,需要时经调制再转发出去,延迟量τ对应一个欺骗距离Δd = c·τ/2,c是光速。单站场景里,这个欺骗距离直接叠加在真实目标距离上,雷达看到的假目标方向与真实目标一致,只是距离被平移了一段。
多站场景比这复杂得多。干扰机转发一次信号,每个接收站各自从自己的角度测量距离,由于各接收站空间位置不同,同一转发延迟在不同接收站产生的假目标并不是同一个点,而是分别落在“接收站-干扰机”连线的延长线上,且假目标到干扰机的距离近似等于Δd/2。这个“近似等于Δd/2”的成立条件是各接收站相距不远,这也是论文所有推导的前提。若接收站间距拉到几十公里,双程几何关系改变,这个近似就要重新推导,不能照搬结论。
这里有个容易绕进去的点:为什么不是Δd而是Δd/2?单站雷达里转发延迟τ对应双程距离差Δd,而多站模型里各接收站测的是干扰机到接收站的单程路径加上DRFM延迟带来的额外单程路径差,推导到最后落在假目标与干扰机的间距上时,自然就剩Δd/2。我建议不要死记这个系数,直接看后面的坐标公式,用向量算一遍就明白了。
3.2 两个不等式判据:把“同处一个SRC”变成坐标几何问题
以1发2收为例建立坐标系:接收站1在原点,接收站2在(L, 0),干扰机在(x0, y0),欺骗距离为Δd。接收站1看到的假目标A在干扰机与接收站1连线的延长线上;接收站2看到的假目标B在干扰机与接收站2连线的延长线上。设干扰机到接收站1的距离为ρ1,到接收站2的距离为ρ2,则:
- A点到接收站1的距离 = ρ1 + Δd/2;
- A点到接收站2的距离,需要根据三角形边角关系由余弦定理推出;
- B点到接收站2的距离 = ρ2 + Δd/2;
- B点到接收站1的距离,同样由余弦定理推出。
有效欺骗条件就是上一章的两个不等式:|AR1 - BR1| ≤ r1 且 |AR2 - BR2| ≤ r2,其中r1、r2是两个接收站各自的距离分辨率。这个判据的物理意义很直白:同一个假目标,在两个接收站眼里必须出现在“差不多”的位置,差得超过一个距离分辨单元,空间配准时就会被当成两个不同目标。
实际工程里,我不太建议每次都手推余弦定理公式。更稳的做法是先在坐标系里写出A、B两点坐标,再用距离函数直接算。A点的计算方法是:从接收站1指向干扰机的单位向量,乘上(ρ1 + Δd/2)后叠加到接收站1坐标上,就是A点坐标;B点同理,从接收站2指向干扰机的单位向量乘上(ρ2 + Δd/2),再叠加到接收站2坐标上。这种向量写法比三角函数更直观,也更容易改写成任意接收站数量的通用代码,后面的仿真脚本就是按这个思路写的。
3.3 远场极限值γ1、γ2:距离差上限只由布站和干扰机位置决定
论文最关键的一步是证明:随着Δd增大,|AR1 - BR1|和|AR2 - BR2|不会无限增大,而是各自收敛到极限值γ1、γ2。极限值由两个因素决定:一是两接收站之间的距离L,二是干扰机位置坐标(x0, y0),与转发延迟量本身无关。
这一点相当反直觉:转发时延明明是干扰方自己控制的,增大转发延迟可以让假目标越推越远,但两个接收站看到的假目标彼此之间的“距离差”却是有上限的。几何上的解释是:假目标被推远的同时,两个接收站视角下的假目标方向逐渐趋近于平行,最后两站测得的距离差稳定在一个固定值。论文里把这个极限值写成了两项:两站到干扰机距离差,叠加一项与基线在视线方向投影有关的修正量,所以极限值不是简单的ρ2 - ρ1,而是经过投影修正后的结果。
原文用L = 1 km、干扰机位于(10, 30) km做了数值验证,得到两个极限值分别约为14.5092 m和14.3707 m。如果两站距离分辨率都大于这两个极限值,那就意味着无论转发时延取多大,假目标都稳稳落在同一个SRC内。论文的仿真曲线显示,两个距离差随Δd增大时单调上升并逐渐逼近极限,上升的“速度”在Δd较小时最快,后面越来越平缓。
据此可以画出等高线:在干扰机可能出现的区域网格上逐点计算γ1、γ2,再画出与距离分辨率相等的等高线。L = 1 km时,大部分区域两个极限值都小于10 m;L = 1.5 km时,大于10 m的区域明显扩大。这说明一个硬结论:接收站间距越大,远场“任意时延都能骗”的干扰机可布局区域越收缩;距离分辨率越高,这个区域同样收缩。论文还指出,当接收站不止两个时,应以相距最远的两站作为判据基准,因为它们的基线最长、极限值通常也最大,这对最苛刻的组合都成立,其余组合自然成立。
物理解释也不难:SRC是扁长的,径向分辨率远好于横向分辨率,而转发式干扰产生的假目标分布在“干扰机-接收站”延长线的横向方向,不同接收站看到的假目标横向分散程度天然小于径向分散程度,所以反而容易落进同一个扁长单元里。这解释了为什么单站雷达里效果显著的转发式干扰,到了多站雷达这里反而要重新讨论。
4. 近场区域的欺骗距离上限:公式化简与数值仿真验证
4.1 远场/近场的边界:用γ与距离分辨率的比较画等高线
远场结论成立的条件是γ1、γ2都小于等于对应距离分辨率。如果干扰机处在γ等高线以内,比如γ1或γ2中任何一个大于对应r,说明在不加约束的情况下,时延取到无穷大时假目标终会分开到不同SRC。但这不代表一定骗不了——距离差随Δd是连续增大的,只要把Δd限制在某个上限内,距离差仍然压得住。
所以近场策略的本质是:用限制最大转发时延来换取欺骗成功。论文给出的操作路径是:根据两个不等式分别反解出各自允许的Δd上限,然后取两者中的较小值。近场区域里,干扰机离各接收站的距离不同,两个不等式一般不会同时达到临界,取交集的操作保证了最苛刻的那个站也不会把假目标剔除。
从仿真角度看,判断远场还是近场只需要一步:在干扰机坐标处计算γ1、γ2,并与r1、r2比较,任何一个超限就进入近场流程。这个判断不需要跑复杂的搜索,几行代码就能完成。我自己习惯把这一步放在所有转发干扰仿真最前面,输出一个“远场/近场”标志位和一个提示信息,避免后续参数设置时忘了当前所处的区域。
4.2 最大欺骗距离:两个不等式各自解出上限后取交集
论文在近场条件下对两个不等式做了化简,式(14)、式(15)分别给出从接收站1和接收站2角度算出的Δd上限。完整解析式里包含ρ1、ρ2、L、x0以及对应的距离分辨率,形式不算复杂,但每一项的符号要非常小心——方向差一点点,算出来的上限可能差一个量级。我不建议手算,把公式抄进脚本里让程序代劳更稳妥。
从工程角度,我更关心这个上限的两个性质。第一,Δd上限对距离分辨率是线性敏感的,分辨率收紧一半,允许的欺骗距离大约也收紧一半;这意味着抗干扰方提高距离分辨率是直接压缩欺骗空间的硬手段。第二,Δd上限随干扰机到多站雷达距离增大而增大,论文的仿真图用指数坐标展示了近场区域的最大欺骗距离分布,结论一致:干扰机离得越远,近场区域内的最大欺骗距离越大,欺骗越容易成功。
最终判据是Δd < min(Δd1, Δd2)。只要转发时延产生的欺骗距离小于这个最小值,两个不等式同时满足,假目标就能留在同一个SRC内。注意这里说的是“欺骗距离”,不是转发延迟τ,两者相差因子2,换算时别搞混。DRFM设备上设置的延迟参数要除以2再和这个阈值比较。
4.3 一个可复现的数值仿真:二分法求最大转发时延
公式看着繁琐,我会直接用几何脚本验证,顺便把近场上限算出来。下面这个脚本基于向量法求A、B点坐标,再用二分法找出近场条件下允许的最大Δd。输入只有4组参数:接收站间距L、干扰机位置(xj, yj)、两个接收站的距离分辨率r1、r2。
import math def unit_vec(px, py, qx, qy): # 从点p指向点q的单位向量 dx, dy = qx - px, qy - py n = math.hypot(dx, dy) return dx / n, dy / n def target_point(rx, ry, jx, jy, dd): # 接收站视角的假目标:在接收站-干扰机延长线上,距干扰机 dd/2 ux, uy = unit_vec(rx, ry, jx, jy) rj = math.hypot(jx - rx, jy - ry) return jx + ux * dd / 2.0, jy + uy * dd / 2.0 def dist(ax, ay, bx, by): return math.hypot(ax - bx, ay - by) def diff_at_dd(L, jx, jy, dd): # 接收站1在原点,接收站2在(L,0) ax, ay = target_point(0.0, 0.0, jx, jy, dd) bx, by = target_point(L, 0.0, jx, jy, dd) d1 = abs(dist(ax, ay, 0.0, 0.0) - dist(bx, by, 0.0, 0.0)) d2 = abs(dist(ax, ay, L, 0.0) - dist(bx, by, L, 0.0)) return d1, d2 def max_deception_range(L, jx, jy, r1, r2, hi=5000.0, tol=0.01): lo = 0.0 # 距离差随dd单调递增,满足条件的dd区间是[0, max_dd] while hi - lo > tol: mid = (lo + hi) / 2.0 d1, d2 = diff_at_dd(L, jx, jy, mid) if d1 <= r1 and d2 <= r2: lo = mid else: hi = mid return (lo + hi) / 2.0 # 参数:L=1km,干扰机在(0.3,0.8)km,两站距离分辨率均为10m L = 1000.0 jx, jy = 300.0, 800.0 r1 = r2 = 10.0 dd_max = max_deception_range(L, jx, jy, r1, r2) print("max dd =", dd_max, "m")脚本逻辑分三层:先定义“假目标在干扰机与接收站连线延长线上、与干扰机距离为Δd/2”的向量计算函数;再对一个候选Δd,分别算出A、B两点,计算两个接收站视角下的距离差;由于距离差随Δd单调增加,满足条件的区域一定是从0开始的一段连续区间,二分法可以稳定找到边界。
参数说明:hi = 5000.0 是搜索上界,表示只关心5000 m以内的欺骗距离;tol = 0.01 是二分精度,实际工程中1 m精度足够,不需要更细。输出结果会打印一个具体数值,以我给出的坐标跑出来大概是几十米量级;但你换成自己的布站坐标后,结果会完全不同。干扰机稍微挪动几十米,Δd上限可能差出两三倍,这也是我不建议用某个固定经验值的原因。把脚本里的干扰机坐标换成真实场景,跑一遍再定延迟量,才是稳妥做法。
5. 常见问题与踩坑:多站欺骗仿真中容易翻车的五个点
5.1 把“远场任意时延必骗”当成全局结论
现象:看到远场证明后,做转发式干扰仿真一律不控延迟,近场场景也直接上大延迟,结果假目标在空间配准阶段被剔光,还以为是多站雷达抗干扰能力太强。
原因:远场必胜的前提是γ1、γ2都小于等于各自距离分辨率。干扰机一旦落入等高线以内区域,极限值超过分辨率,任意时延就不再成立,这是边界条件被忽略导致的问题。
解决:先算γ1、γ2,再决定策略。把第一步写成脚本:输入L和干扰机坐标,输出γ1、γ2以及它们与r1、r2的比值,比值全小于1才允许放开延迟,否则走近场最大欺骗距离流程。这一步在4.1节的脚本框架里顺手就能实现,没必要省。
5.2 忽略两站距离分辨率不同
现象:一号站距离分辨率10 m,二号站距离分辨率30 m,按10 m的标准去推最大欺骗距离,结果二号站那边很快失效,联合检测时假目标被剔除。
原因:判据是两个不等式取交集,必须分别满足。只盯分辨率更高的一站,另一站的分辨率更宽松或更严格,都会造成判据偏差,不能用一个r值代替所有站。
解决:两个站分别用各自的r1、r2计算Δd1和Δd2,再取min(Δd1, Δd2)。4.3节脚本里,r1和r2要作为独立参数传入,而不是写死成同一个值。这是最常见的参数化错误,尤其在接收站数量超过两个时,每个站的分辨率都要单独维护。
5.3 用单站DRFM延迟量直接代入多站模型
现象:单站雷达里延迟τ对应200 m欺骗距离,拿到多站模型里原样使用,假目标位置对不上,仿真里连假目标都找不到。
原因:单站场景里假目标方向就是视线方向,延迟直接换算距离;多站场景里每个接收站看到的假目标在各自延长线上,且距离增量近似Δd/2,不是Δd。方向一变,距离差分布完全不同。
解决:不要用“延迟→距离”的直觉,先做向量几何换算,再用4.3节脚本验证。转发延迟与欺骗距离之间因子2的换算,在单站里是常识,在多站里却要重新过一遍坐标系,把它当成两个独立问题处理最稳妥。
5.4 多站场景选错基准站
现象:对几个接收站单独仿真,每个站都显示假目标正常,但整系统联合处理时假目标被剔除,问题表现得很玄学。
原因:欺骗成功与否不能靠局部视角替代。论文明确指出,多接收站时要选择相距最远的两部接收站作为基准;如果最远两站都能被欺骗,整个多站系统基本都能被欺骗。很多仿真只验证相邻站,漏了最远站对。
解决:把接收站两两组合后,选距离最大的那对代入不等式。接收站数量较多时,建议直接遍历所有接收站对,找到距离最大的一对再套脚本,不要凭感觉选站。
5.5 把近场最大欺骗距离当成常数
现象:上一次仿真算出Δd上限是几十米,换了个布站方案后还继续用同一个值,结果假目标时灵时不灵,参数调半天找不到规律。
原因:Δd上限是干扰机坐标、站距、距离分辨率共同决定的曲面,不是常数。L从1 km变到1.5 km,10 m等高线内区域明显扩大,最大欺骗距离也跟着变,固定值当然失效。
解决:每次改布站参数,重跑一遍扫描。把干扰机位置网格化,将Δd_max存成二维表,仿真时直接查表,避免每次重复二分。这也是把论文结论落到工程里的推荐方式。
6. 从论文结论到工程判断:等高线安全区与布站设计技巧
6.1 干扰方视角:先画等高线,再谈时延
把论文的结论整理成一张判据表,贴在工作位上比反复翻PDF好用得多:
| 场景 | 成立条件 | 转发时延策略 |
|---|---|---|
| 远场 | γ1 ≤ r1 且 γ2 ≤ r2 | 任意时延,放开转发 |
| 近场 | γ1 > r1 或 γ2 > r2 | Δd < min(Δd1, Δd2),限制时延 |
| 多站 | 接收站多于2部 | 以相距最远的两部接收站为基准 |
实际操作顺序是固定的:第一步,在干扰机可能出现的区域网格上计算γ1、γ2;第二步,用距离分辨率画等高线,确定干扰机落在等值线内还是外;第三步,线外直接放开延迟,线内用二分脚本求Δd上限。这套流程我后来写成了固定模板,每次项目复现都先跑一遍,不再凭经验定转发延迟。
6.2 抗干扰方视角:基线、距离分辨率与假目标筛查
抗干扰方读这篇论文的收获更直接。结论反向翻译过来就是:要压缩转发式干扰的可用空间,把接收站基线拉大,或者提高距离分辨率。拉大基线会让γ等高线向外扩,干扰机可布局的“安全区”变小;提高距离分辨率等价于把10 m等高线向外推,能让更多原本放松的区域变成受限区。
另一个可落地的技巧是:既然最大欺骗距离受限于近场判据,抗干扰系统可以在时间维上检测异常长的转发延迟。如果假目标对应的距离超过了当前布站几何条件下允许的Δd上限,直接判为欺骗。这个逻辑不需要复杂的特征提取,只需要一张布站几何表和一次距离测量值,工程实现成本很低。我在做多站雷达抗干扰仿真时,把这个判据加进了融合中心的预筛选模块,密集假目标场景下的虚警率明显下降。
从那以后,我每次做多站雷达抗干扰仿真,都强制走一遍“画γ等高线→标距离分辨率→算Δd上限”的流程,再也没出现过延迟失控导致整轮仿真空跑的情况。这篇论文的推导和仿真条件都列得比较完整,值得下载下来把式(14)到式(16)对着抄一遍,用起来会更踏实。希望帮到你。
本文还有配套的精品资源,点击获取