罗 威,王良斯,林贤斌
(空军预警学院六系,湖北 武汉 430019)
半主动激光制导武器具有制导精度高、抗干扰能力强、效费比高等优点,在现代战争中大量用于对重要目标的空地精确打击。激光角度欺骗干扰是对付激光制导武器的有效手段[1],其前提条件是对激光制导信号的侦察告警和解码识别。目前,对单一激光制导信号的解码技术和激光告警设备发展较为成熟[2]。然而,现代战争空对地打击行动中,如攻击机场、桥梁、坦克集群时,就涉及发射多枚激光制导导弹同时打击多个目标的情况,这时将产生在时域、空域、频域交叠重合的混合激光制导信号。若激光告警设备不能对其进行正确分选和解码,就无法实施有效的激光角度欺骗干扰以保护重要目标免遭打击。
在混合激光制导信号分选与解码技术研究方面,基于信号脉冲间隔、脉宽等时间域特征,先后有文献提出序列差值直方图和变步长序列搜索法[3]、加权归一自相关方法[4]、码元周期搜索法[5]、基于统计方法[6]、基于频率方法[7]、基于深度自编码网络方法[8]等系列方法。但是,这些方法一般都是先识别相对简单或具备先验知识的单路信号,而后再对混合信号进行分选。本文基于骨架周期分析方法,编写了相应MATLAB程序,可实现对无先验知识、相对复杂混合激光制导信号的自动分选和解码。
激光编码技术,是激光半主动制导武器避免自扰和对抗转发式激光干扰的关键技术。激光脉冲编码类型包括:脉宽编码、精确频率码、二间隔码、脉冲调制码、有限位随机周期码、伪随机编码等。
目前激光制导武器中主要采用的是精确频率码和有限位随机周期码。精确频率码即脉冲重复频率(Pulsed Repetition Frequency,PRF)激光编码,是指激光脉冲的频率重复不变,相邻激光脉冲间隔为固定值。有限位随机周期码即脉冲时间间隔调制(Pulse Interval Modulate,PIM)激光编码,是指各激光脉冲时间间隔是随机无关联的,但其位数有限,具有可重复性。例如设置一个4位PIM激光编码,对应有4个脉冲间隔{t1,t2,t3,t4},其重复周期也可称为骨架周期,为同一周期内各脉冲间隔之和,即T=t1+t2+t3+t4。
本文主要针对PRF信号和PIM编码信号所组成的混合激光制导信号开展研究。混合信号包括PRF与PRF混合信号、PRF与PIM混合信号、PIM与PIM混合信号等三种类型。混合信号时序分布分别如图1、图2、图3所示。信号1重复周期为T1,信号2重复周期为T2,假设信号1先到达激光告警器,两组信号的达到时间差为Δt。
图1 PRF与PRF混合信号示意图Fig.1 PRF mixed PRF signals
图2 PRF与PIM混合信号示意图Fig.2 PRF mixed PIM signals
图3 PIM与PIM混合信号示意图Fig.3 PIM mixed PIM signals
由于激光制导武器所使用的激光目标指示器性能相同,因此两组激光制导信号的频率、脉宽相同,照射方位也基本一致,故两组信号会形成一个混合信号,且很难直接进行区分。混合信号会形成一个新的骨架周期Tm,Tm为T1和T2的最小公倍数。为避免自扰影响,T1和T2会避免设置为简单倍数关系,甚至互为质数,故大部分情况下,Tm为一较大值,即两组信号混合后的激光脉冲间隔近似为伪随机序列,脉冲间隔值在很长时间内都不重复。
另外,当两组信号的达到时间间隔较大且满足|Δt|>2T1时,激光告警器在信号2到达之前就可以完成对信号1的解码,对后续混合信号的分选与解码难度也会相应减小。本文中主要考虑的是|Δt|<2T1这种更复杂的情况,上述情况也可兼容处理。
当混合信号中包含PRF信号时,即对于PRF与PRF混合信号、PRF与PIM混合信号,其分选和解码相对简单,本文采用时差矩阵统计方法进行处理。对于PIM与PIM混合信号,其分选和解码相对复杂,本文针对性提出骨架周期分析方法,较好地解决了问题。
两个PRF信号组成混合信号,其重复周期分别为T1和T2,如图1所示。假设共有L个脉冲,每个脉冲信号的到达时间XTOA已知,设ΔTmn为第m个脉冲和第n个脉冲到达时间XTOAm、XTOAn之差,则有:
ΔTmn=XTOAm-XTOAn(m>n)
(1)
从L个脉冲中选取相邻的k个脉冲,可以构建一个(k-1)阶的时差矩阵ΔT:
如果第m个脉冲和第n个脉冲是同一个重频信号相邻脉冲,则此时的ΔTmn=T1或T2。如果不是,则ΔTmn为非相关值,且基本不重复。对该时差矩阵进行统计分析,在统计结果中出现次数最多的值即为T1或T2。
k的取值会直接影响信号分选准确度和分选时间。当k取值较小时,时差中相关值可能偏少,影响信号分选准确度,一般取k>5。k取值也不宜过大,否则运算量会大大增加,分选时间也相应增大。在时差矩阵中还会出现很多T1和T2的倍频值,在时差矩阵基础上,还需要对时差矩阵的统计结果进一步处理,将统计结果作为第二个矩阵,称为统计矩阵,对统计矩阵进行相关运算,可以更加准确提取T1或T2。
对于PRF与PIM混合信号,同样可利用上述矩阵统计方法,建立时差矩阵对混合信号进行分选,从中分选出重频信号,剩下的信号就是编码信号,在此不再详述。也可用下文即将介绍的动态自相关算法对混合信号进行分选,从中分选出编码信号,剩下的信号就是重频信号。
假设PIM与PIM混合信号如图3所示,其中一组为4位编码信号,脉冲间隔为t1、t2、t3、t4,骨架周期为T1=t1+t2+t3+t4,另一组同样为4位编码信号,脉冲间隔为t5、t6、t7、t8,骨架周期为T2=t5+t6+t7+t8。如前文所述,混合信号形成的新的骨架周期很长,相邻激光脉冲间隔近似为伪随机序列,即脉冲间隔值基本不重复。若同样建立时差矩阵并进行统计分析,相比包含PRF信号的混合信号,分选出t1、t2、t3、t4这4个值要困难很多,原因一是要分析的脉冲数量大大增加,二是统计结果差距不显著。
因此,本文将骨架周期作为PIM与PIM混合信号分选的主要参数。其流程描述如下:(1)确定分选脉冲数动态范围(x,y)。以4位编码信号为例,其骨架周期最少需要5个脉冲序列构成,即最小值x取值5即可。根据实际情况,激光指示器无论是PRF信号还是PIM信号,其脉冲间隔大小一般不会有数量级的差异,因此y取值为x的2倍或略大即可,就可以保证(x,y)内存在一个完整的骨架周期。(2)建立时差矩阵。从第n个脉冲开始,在(x,y)的脉冲长度范围内,构建时差序列,注意此时时差并非是相邻脉冲的时差,而是相距x个以上、y个以下脉冲的时差,再从(n+1)个脉冲开始,重复上述流程,形成时差矩阵。(3)统计分析获取骨架周期。对时差矩阵进行统计分析,可以得到编码信号的骨架周期T1和T2。(4)根据骨架周期进行信号分选。从第1个脉冲序列开始分析,其到达时间加上T1,如果在混合信号脉冲序列中存在相应脉冲,则说明这2个脉冲均属于信号1的脉冲;如果不存在,则该脉冲属于信号2。遍历完整混合信号脉冲序列,就可以完成混合信号的分选。(5)信号解码。对分选后得到的信号进行解码,得出其位数和脉冲间隔。
根据上述混合激光制导信号分选与解码方法,基于MATLAB软件平台编写了混合激光制导信号分选与解码软件,软件界面如图4所示。其功能包括:模拟产生PRF与PRF、PRF与PIM、PIM与PIM等三种类型混合信号,对混合信号进行自动分选和解码,分选与解码结果的图文可视化显示。
图4 混合激光制导信号分选与解码软件界面图Fig.4 Interface of the sorting and decoding software
以相对最复杂的PIM与PIM混合信号为例介绍混合信号产生操作,其信号设置和产生的混合信号如图5所示。其中PIM信号1的骨架周期是T1=69 ms+85 ms+67 ms,PIM信号2的骨架周期是T2=42 ms+65 ms+81 ms。
图5 PIM与PIM混合信号生成及可视化Fig.5 PIM mixed PIM signals generation and visualization
模拟产生PRF与PRF、PRF与PIM、PIM与PIM等三种类型混合信号,每种类型各100组,并利用混合激光制导信号分选与解码软件分别对其进行分选与解码,其中信号1和信号2达到时间间隔随机生成。对100组PRF与PRF混合信号的典型分选与解码结果如图6所示,平均分选用时0.453 s。对100组PRF与PIM混合信号的典型分选与解码结果如图7所示,其中PIM信号为3位编码信号,平均分选用时1.190 s。对100组PIM与PIM混合信号的典型分选与解码结果如图8所示,PIM信号均为3位编码信号,平均分选用时2.206 s。
图6 PRF与PRF混合信号分选与解码结果Fig.6 Sorting and decoding result of PRF mixed PRF signals
图7 PRF与PIM混合信号分选与解码结果Fig.7 Sorting and decoding result of PRF mixed PIM signals
图8 PIM与PIM混合信号分选与解码结果Fig.8 Sorting and decoding result of PIM mixed PIM signals
本文主要研究混合激光制导信号分选和解码方法,并编写了混合激光制导信号分选与解码软件,对PRF与PRF、PRF与PIM、PIM与PIM等三类混合信号进行了分选和解码。其中,PIM信号由于包含编码位数和脉冲间隔值等多重未知参数,分选难度较大,本文利用骨架周期分析方法,能够正确分选包含PIM信号的混合信号,且分选和解码总用时较短,能够较好满足实际需求。研究工作为下一步在实际环境下对混合激光制导信号进行分选与解码打下了良好基础。