一种适用于无人机通信链路的干扰方法

1.本发明涉及通信干扰技术，特别涉及无人机通信链路的干扰技术。
技术背景
2.随着无人机相关技术的发展，无人机在现代化社会中发挥着越来越重要的作用，尤其是在城市管理、农业、地质等行业。但在给人们带来便利的同时，无人机的不规范使用也暴露出一系列公共安全隐患。基于无人机飞行管控的要求，国家空管委、中国民航局等部门出台了一系列无人机监管政策，电磁干扰是反无人机系统中不可或缺的一部分。根据无人机系统特点，干扰方式分为阻塞干扰、跟踪干扰以及多频连续波干扰。阻塞干扰是一种预先不需要知道跳变频率，将干扰分布在整个或者部分信号频段的定频干扰方式。跟踪干扰是一种预先需要知道跳变频率，在其通信信道实施窄带干扰的干扰方式。多频连续波干扰是一种预先不需要知道跳变频率，将干扰分布在部分特定频率上的干扰方式。本专利包括信号帧同步、中心频率估计、发送干扰信号这三个模块，不需要预先知道跳变规律，即可在无人机跳频通信信号出现的时间和频点处实施干扰，既能提高干扰效率、节省功率，也能减少对其他通信设备的影响。
3.在现有的无人机信号检测方法如下：
4.1)带宽阻塞干扰
5.这种方法只需要知道跳频信号的频段信息，不需要知道跳频信号的跳变规律，一般采用随机噪声波形，使用大功率对整个或部分信号频段进行压制。该方法的干扰效果完全由干扰能量大小来决定。若跳频系统采用bfsk调制方式，则系统的误码率为
[0006][0007]
其中，p为系统误码率，eb为每个码元的能量，n0为背景热噪声功率谱密度，nj为干扰机产生的高斯白噪声功率谱密度。
[0008]
可见，只要能够达到一定的干信比，宽带阻塞干扰能使跳频系统达到一定的误码率。
[0009]
2)部分频带阻塞干扰
[0010]
这种方法只需要知道跳频信号的频段信息，不需要知道跳频信号的跳变规律，与宽带阻塞干扰相比，该方法降低了频带覆盖率，提升了部分频带干扰信号的功率密度，通过提升覆盖频段内的干扰效果，达到干扰目标信号的目的。该干扰的干扰效果与干扰机带宽与总跳频带宽的比例有关。
[0011]
干扰机驻留的频谱中的误码率为
[0012]
[0013]
其中，ec为码元的能量，n0为背景热噪声功率谱密度，是提高了的干扰机噪声密度。
[0014]
可见，在相同的干扰功率下，选择合适的干扰频带覆盖比例r，部分频带阻塞干扰能使跳频系统达到一定的误码率。
[0015]
3)跟踪干扰
[0016]
这种方法通过对信号进行快速的侦察分析，在较短的时间内对准跳频频点并实施窄带干扰。该方法需要知道跳频信号的跳变规律，在跳频通信过程中对驻留的一个频点进行截获分析并进行窄带干扰，当跳频系统频率跳变后，干扰系统也随之改变干扰频率。当跳频系统属于慢速跳频系统，干扰机频率对准时间小于跳频驻留时间，跟踪式干扰机能使跳频系统达到一定的误码率。
[0017]
4)多频连续波干扰
[0018]
这种方法只需要知道跳频信号的频点信息，不需要知道跳频信号的跳变规律，一般采用窄带噪声波形的组合，对准跳频频点进行干扰。就采用正交mfsk调制的跳频通信系统而言，干扰频点数为
[0019][0020]
其中，j为干扰机发射总功率，p为收到的信号功率。为向下取整。
[0021]
当多频连续波对一个跳频信号进行干扰时，误符号率为
[0022][0023]
其中，w为总跳频带宽，wd为跳频间隔。
[0024]
可见，当干扰机发射总功率小于信号收到的功率时，不会产生传输错误；当干扰机发射总功率大于信号收到的功率，且干扰噪声在跳频带宽内时，会产生传输错误。
[0025]
在现有的无人机信号检测方法中，带宽阻塞干扰的干扰功率被平均分配到整个通信频带上，跳频系统的跳频带宽越宽，单位频带上的干扰功率越小，跳频系统受到的干扰也就越小，而且如果跳频通信信号功率较高且频段较宽，那么需要消耗较大的干扰能量。部分窄带阻塞干扰中，如果想要达到较好的干扰效果，干扰覆盖的频段需要达到足够的比例，若被干扰目标使用了自适应技术，那么干扰方式将失效。现有的跟踪式干扰中，需要预先知道跳频信号的跳变规律，而且在快速跳频系统，如果干扰机频率对准时间大于跳频驻留时间，那么跟踪式干扰将失效。


技术实现要素：

[0026]
本发明所要解决的技术问题是，针对目前无人机干扰技术在跳频通信信号功率较高且频段较宽或者预先对跳变规律未知时干扰效果较差且需要消耗较大的干扰能量的情况，提供一种适用于无人机反制的跟踪干扰方法。
[0027]
本发明为解决上述技术问题所采用的技术方案是，一种适用于无人机通信链路的
干扰方法，包括步骤：
[0028]
1)信号帧同步步骤：
[0029]
1-1)从信号检测系统中获取信号帧结束时刻t
frame_end
计算同步时刻t
syn
；
[0030]
1-2)以同步时刻t
syn
为基准采集信号，求得采集信号对应的功率谱从而得到同步时刻前功率谱变化值p
differ_pre
、同步时刻处功率谱变化值p
differ_mid
以及同步时刻后功率谱变化值p
differ_pos
；
[0031]
判断p
differ_mid
是否大于p
differ_th
，如是，认为在同步时刻存在跳变信号帧，则进入步骤1-3)；否则更新t
syn
＝t
syn
period
fh
t
fft
，返回步骤1-2)；
[0032]
1-3)计算同步时刻的调整量δt
syn
，并判断当前是否同步成功，如是，则进入步骤2)，否则更新t
syn
后返回步骤1-2)，t
syn
＝t
syn
δt
syn
；
[0033]
2)中心频率估计步骤：
[0034]
分别在时刻t
syn
period
fh-t
frame
、t
syn
period
fh-t
frame
2t
fft
处采集信号并求得功率谱spec6(f)、spec7(f)；再计算功率谱差分spec
differ_pilot
(f)＝abs(spec6(f)-spec7(f))；period
fh
为配置的跳频周期，t
fft
为配置的一次fft运算的信号时长；abs表示取绝对值；
[0035]
将满足功率谱差分spec
differ_pilot
(f)最大值时的频率作为导频频率f
pilot
；导频频率f
pilot
为发射干扰的中心频率；
[0036]
3)发送干扰信号步骤：在中心频率处发送时长小于信号帧时长t
frame
的干扰信号。
[0037]
本发明通过构建信号帧同步、中心频率估计、发送干扰信号这三个步骤，利用同步时刻功率谱变化值和功率谱变化门限的关系来判定跳变信号帧是否存在，不需要知道跳频信号的跳变规律，使用导频估计后得到的导频频率生成特定的干扰信号使得干扰有效，可以在较短的时间内对准跳频频点并实施窄带干扰，既能提高干扰效率、节省功率，也能减少对其他通信设备的影响。
[0038]
本发明的有益效果是，能节省干扰功率，提高干扰精度；计算复杂度低，能够快速响应跳频信号；干扰系统中参数可调，适用于不同型号的无人机。
附图说明
[0039]
图1为本发明示意图。
具体实施方式
[0040]
实施步骤如图1所示：
[0041]
1)信号帧同步步骤：
[0042]
估计信号帧结束时刻，若同步到精确信号帧结束时刻则进入步骤2)，步骤1)具体过程如下：
[0043]
1-1)信号帧结束时刻确定：加载干扰信号基带波形，配置跟踪参数，其中flag
track
置零。跟踪参数包括以下几个内容
[0044]
参数名称参数说明f
start
信号起始频点f
end
信号终止频点period
fh
跳频周期
p
differ_th
功率变化门限flag
track
跟踪标记位t
frame
信号帧时长t
fft
一次fft运算的信号时长th
track
跟踪指标的门限
[0045]
从信号检测系统中获取最新一次信号帧结束时刻t
frame_end
，令同步时刻
[0046]
1-2)跳变信号帧判定：在t
syn
、t
syn
t
fft
、处采集信号，分别求信号的功率谱得到spec0(f),spec1(f),spec2(f),spec3(f),spec4(f),spec5(f)。本方法最终会将信号帧结束时刻同步在spec1(f)结束、spec3(f)开始的时刻。对功率谱进行谱运算，计算表达式如下
[0047][0048][0049][0050]
其中，p
differ_pre
表示同步时刻前功率谱变化值，p
differ_mid
表示同步时刻处功率谱变化值，p
differ_pos
表示同步时刻后功率谱变化值。判断p
differ_mid
是否大于p
differ_th
。若同步时刻功率谱变化值大于功率谱变化门限，则认为在同步时刻存在跳变信号帧，跳转到步骤1-3)；否则认为在该时刻并没有无人机帧信号发生跳变，此时令t
syn
＝t
syn
period
fh
t
fft
，跳转到步骤1-2)进行下一次同步。
[0051]
1-3)同步：计算跟踪指标值indicate
track_0
，indicate
track_1
表达式如下
[0052][0053][0054]
根据跟踪指标值indicate
track_1
更改下一次的t
syn
值，t
syn
＝t
syn
δt
syn
。
[0055]
δt
syn
具体表达式为
[0056][0057]
判断是否满足indicate
track_0
》th
track
，如是，则flag
track
＝flag
track
1。反之则flag
track
＝0。当flag
track
》5时，认为同步成功，在每一次跟踪过程后进入步骤2)。
[0058]
3)中心频率估计步骤：
[0059]
在t
syn
period
fh-t
frame
、t
syn
period
fh-t
frame
2t
fft
采集两次信号，分别求信号的功率谱得到spec6(f)、spec7(f)。对功率谱进行谱运算spec
differ_pilot
(f)，计算表达式如下
[0060]
spec
differ_pilot
(f)＝abs(spec6(f)-spec7(f))
[0061]
abs表示取绝对值；
[0062]
通过功率谱差分求导频频率，计算表达式为：
[0063]fpilot
＝argmax(spec
differ_pilot
(f))
[0064]
得到导频频率f
pilot
后，将导频频率f
pilot
作为发射干扰的中心频率；
[0065]
3)发送干扰信号步骤：在该中心频率处发送时长小于t
frame
的干扰信号实施干扰，之后跳转到步骤1-2)，继续跟踪步骤。