一种适用于宽带抗干扰系统的信号捕获方法与流程

1.本发明涉及卫星通信领域中的一种适用于宽带跳频系统的信号捕获方法。


背景技术：

2.宽带抗干扰系统是保障军事通信安全可靠的重要手段，目前常用的宽带抗干扰系统，通常采用宽带跳频体制。跳频卫星通信具有通信距离远、覆盖面积大、机动灵活、抗干扰能力强等优点，因此在应急通信、海事通信、机动通信等方面具有广泛的作用，是卫星通信的重要发展方向。
3.地面终端接收卫星跳频信号时，首先要进行天线对星。天线对星，是使天线波束中心对准目标卫星的过程，它是地球站进行入网验证和日常业务传输的基础。天线只有准确地对准目标卫星，才能保证地球站稳定可靠地工作。但是，在跳频系统中，由于对星时跳频图案尚未同步，因此收到的信号时有时无，这就导致接收的信号呈现出非平稳状态，很难为自适应对星提供稳定的功率参考，使得地面终端自适应对星难度增加。
4.因此，捕获到跳频信号，并进行功率估计，然后输出稳定的功率值，为天线对星提供可靠地功率反馈，是跳频卫星通信的首要任务。跳频信号是一个非平稳、非线性、随时间变化的信号，一般的fft分析方法达不到分析非平稳信号的要求。如典型的线性时频变化：gabor展开，短时傅立叶变换(stft)，小波变换(wt)等，满足不了捕获跳频信号的需求。


技术实现要素：

5.本发明的目的在于，针对上述背景技术中的缺点，提供一种跳频同步前信号捕获方法。本发明提出了跳频信号快速捕获和功率估计算法，解决了地面终端在跳频图案同步前对星难的问题。本发明采用时域分析、频域分析、带宽检测、功率检测等多维域分析方法，对跳频信号进行了准确的捕获，并通过多跳平均进行功率估计，为对星提供功率参考，解决了跳频同步前由于信号时有时无导致的自动对星难题。该方法具有稳定可靠、检测准确度高、功率估计精度高的优良性能。
6.本发明的目的是这样实现的：
7.一种适用于宽带抗干扰系统的信号捕获方法，用于捕获带宽为b1和带宽为b2的目标信号，包括以下步骤：
8.(1)接收机接收中频信号，选择中心频点f0，将信号变频为零中频信号；
9.(2)对零中频信号进行采样，采样率fs不小于检测带宽b
detec
的2倍，采样后的序列为y(n)；
10.(3)对采样后的序列y(n)，进行n点fft，得到y(k)；其中，n为每跳数据分为m段后，每段数据的长度，m取不小于4的整数；
11.(4)将频谱信息y(k)进行平方运算，得到y(k)2；
12.(5)将每个fft，与前m-1个fft，按对应频点进行相加，输出累加后的结果r1(k)，并存储；
13.(6)对r1(k)进行长度为max(l1，l2)的滑动累加，并输出累加结果r2(k)；其中，l1＝b1/(fs/n)，表示带宽为b1的目标信号的谱线数量，l2＝b2/(fs/n)，表示带宽为b2的目标信号的谱线数量；
14.(7)在步骤(6)输出的累加结果中，找到最大值p
max
，判断p
max
是否大于门限th，如果p
max
》th，则认为找到信号，并记录步骤(5)中存储结果的相应频点位置；
15.(8)通过大于门限的频谱根数nm，判断信号带宽，识别出所找到的信号是哪一个目标信号；
16.(9)将大于门限th的p
max
作为所找到信号的有效功率值并进行存储；
17.(10)将存储的p
max
值累加求平均，输出给天线跟踪系统；
18.完成宽带跳频信号的捕获。
19.进一步的，步骤(7)中，门限th由底噪乘以加权系数来确定。
20.进一步的，步骤(7)中，若有多个相邻p
max
大于门限th，则利用时域特性对相邻的p
max
进行处理，找到相邻p
max
的最大值，去除多余值。
21.进一步的，步骤(7)中，若时域上有连续多个p
max
大于门限，则认为是干扰；若所选择的采样带宽长期存在干扰，则返回步骤(1)改变中心频点，重新执行步骤(1)～(7)。
22.本发明相比背景技术具有如下优点：
23.1、本发明的设计非常适合跳频系统的信号检测，不需要提前进行跳频同步。相对于其他信号检测方法要求工作于固定频点，该方法不受跳频频点的变化的影响。
24.2、本发明能够适应多种速率混合信号的检测，且能够工作于低信噪比的情况。
25.3、本发明采用的fft方法和功率统计方法，适合fpga实现，便于工程应用。
附图说明
26.图1是本发明方法的流程图。
具体实施方式
27.下面对本发明做进一步的详细说明。
28.一种适用于宽带抗干扰系统的信号捕获方法，用于捕获带宽为b1和带宽为b2的目标信号，包括以下步骤：
29.(1)接收机接收中频信号，选择中心频点f0，将信号变频为零中频信号；
30.(2)对零中频信号进行采样，采样率fs不小于检测带宽b
detec
的2倍，采样后的序列为y(n)；
31.(3)对采样后的序列y(n)，进行n点fft，得到y(k)；其中，n为每跳数据分为m段后，每段数据的长度，m取不小于4的整数；
32.(4)将频谱信息y(k)进行平方运算，得到y(k)2；
33.(5)将每个fft，与前m-1个fft，按对应频点进行相加，输出累加后的结果r1(k)，并存储；
34.(6)对r1(k)进行长度为max(l1，l2)的滑动累加，并输出累加结果r2(k)；其中，l1＝b1/(fs/n)，表示带宽为b1的目标信号的谱线数量，l2＝b2/(fs/n)，表示带宽为b2的目标信号的谱线数量；
35.(7)在步骤(6)输出的累加结果中，找到最大值p
max
，判断p
max
是否大于门限th，如果p
max
》th，则认为找到信号，并记录步骤(5)中存储结果的相应频点位置；
36.(8)通过大于门限的频谱根数nm，判断信号带宽，识别出所找到的信号是哪一个目标信号；
37.(9)将大于门限th的p
max
作为所找到信号的有效功率值并进行存储；
38.(10)将存储的p
max
值累加求平均，输出给天线跟踪系统；
39.完成宽带跳频信号的捕获。
40.进一步的，步骤(7)中，门限th由底噪乘以加权系数来确定。
41.进一步的，步骤(7)中，若有多个相邻p
max
大于门限th，则利用时域特性对相邻的p
max
进行处理，找到相邻p
max
的最大值，去除多余值。
42.进一步的，步骤(7)中，若时域上有连续多个p
max
大于门限，则认为是干扰；若所选择的采样带宽长期存在干扰，则返回步骤(1)改变中心频点，重新执行步骤(1)～(7)。
43.下面以捕获带宽为500khz和带宽为1mhz的目标信号为例进行说明：
44.(1)接收机接收中频信号并变频。接收机接收中频信号，按中心频点f0，将信号变频为零中频信号。其中，f0为跳频带宽内的某一频点，可根据需要调整；
45.(2)ad采样。对零中频“信号”(也可能不含信号，只有噪声。此处统称“信号”)采样，采样率fs不小于检测带宽b
detec
的2倍，采样后的序列用y(n)表示；其中，检测带宽b
detec
为跳频带宽内的某一段带宽；
46.(3)fft。对采样后的序列y(n)，进行n点fft，得到y(k)；其中，n为每跳数据分为m段后，每段数据的长度，m取不小于4的整数。
47.(4)计算功率。将fft后的频谱信息y(k)，进行平方运算，得到y(k)2；
48.(5)通过时域滑动相加方式，综合整跳的功率信息。采用将m个fft相加的方式，将每个fft，与前m-1个fft，按对应频点进行相加，输出累加后的结果r1(k)，并存储；
[0049][0050]
(6)频域滑动相加，对r1(k)按一定长度滑动相加，即计算目标信号的功率。目标信号为出现概率较为稳定的某种跳频信号。假设有2种信号存在：带宽分别为500khz、1mhz。对于500khz信号，可认为有l1＝500/(fs/n)根谱线；对于1msps信号，可认为有l2＝1000/(fs/n)根谱线。对r1(k)进行长度为max(l1，l2)的滑动累加，并输出累加结果r2(k)；
[0051]
(7)信号捕获：找功率最大值p
max
，并与门限比较、判决。在步骤(6)输出的(n-l1)个结果中，找到最大值p
max
，然后判断p
max
是否大于门限th(门限根据底噪乘以加权系数来定)。如果p
max
》th，则认为找到信号，并记录步骤(5)存储结果的相应位置；
[0052]
(8)去除p
max
冗余。由于时域上的滑动平均会产生有多个相邻p
max
大于th，因此需要利用时域特性对相邻的p
max
进行处理，找到相邻p
max
的最大值，去除多余值。
[0053]
(9)干扰识别。利用p
max
的时域特点进行干扰判断，时域上连续多个大于门限则认为是干扰。如果所选择的某段采样带宽，长时间存在干扰，则可以通过改变步骤(1)中的中心频点，避开该干扰所在的频点位置。
[0054]
(10)不同带宽信号识别，以及宽带信号剔除。识别两种不同带宽信号。在步骤(6)
的存储数据中，通过大于门限的频谱根数nm，判断信号带宽。
[0055]
(11)存储有效功率值p
max
。
[0056]
(12)将多个有效p
max
值累加求平均，输出给天线跟踪系统。
[0057]
完成宽带跳频信号捕获方法。
[0058]
以下为一个更具体的例子：
[0059]
参照图1，一种适用于宽带抗干扰系统的信号捕获方法，主要步骤如下：
[0060]
(1)接收机接收中频信号并变频。以中心频点f0变频。
[0061]
(2)ad采样。假设采样带宽60mhz，采样时钟128mhz，跳频带宽1.2ghz，跳速12500跳/秒，在60mhz带宽提取跳频信号功率(跳频频点在1.2ghz带宽内随机均匀分布，所以频点落入60mhz带宽概率约为60mhz/1.2ghz＝0.05)；
[0062]
(3)fft。对采样后的信号，每16us进行一次fft运算，即1/5跳做一次fft。进行2048fft点数为2048；每根谱线代表128mhz/2048＝62.5khz；
[0063]
(4)计算各频点功率。将fft后的2048个频谱信息，进行取模、平方运算；
[0064]
(5)综合整跳的功率信息。采用5个fft相加的方法，将每个fft，与前4个fft，按对应频点进行相加，输出累加后的结果。结果为2048的数据，并存储；
[0065]
(6)频域滑动相加，计算信号功率。假设只有2种信号存在：500khz、1mhz。对于500khz信号，有8根谱线；对于1mhz信号，有16根谱线。在2048点上进行长度为16的滑动累加，并输出累加结果；
[0066]
(7)信号捕获。通过比较功率大小进行信号识别。在步骤(6)输出的(2048-16)个结果中，找最大值p
max
，然后判断p
max
是否大于门限th。如果p
max
》th，则认为找到信号，并返回到步骤(5)存储结果的相应位置；
[0067]
(8)由于时域上的滑动平均会产生有多个相邻p
max
大于th，因此需要利用时域特性对相邻的p
max
进行处理，找到相邻p
max
的最大值，去除多余值。
[0068]
(9)干扰识别。利用p
max
的时域特点进行干扰判断，时域上连续20个值大于门限则认为时干扰。去除频点，返回步骤(6)重新计算。
[0069]
(10)不同带宽信号识别，以及宽带信号剔除。识别两种不同带宽信号，在步骤(5)的存储数据中，通过大于门限的频谱根数nth，判断信号带宽。如果nm大于16，则判定为1mhz信号，并剔除。
[0070]
(11)存储有效功率值pmax。
[0071]
(12)将多个pmax值累加求平均，输出给天线跟踪系统。
[0072]
总之，本发明提出了跳频信号快速捕获和功率估计算法，解决了地面终端在跳频图案同步前对星难的问题。本发明采用时域分析、频域分析、带宽检测、功率检测等多维域分析方法，对跳频信号进行了准确的捕获，并通过多跳平均进行功率估计，为对星提供功率参考，解决了跳频同步前由于信号时有时无导致的自动对星难题。该方法具有稳定可靠、检测准确度高、功率估计精度高的优良性能，特别适用于宽带抗干扰系统的信号捕获。
[0073]
本发明书中未进行详细描述的内容为本领域内的公知内容。
[0074]
以上所述仅为本发明的一种具体实施方式，但本发明的保护范围并不局限于此。任何熟悉本技术领域的技术人员根据本发明的技术方案及其发明构思作出的等同替换或改变，都应涵盖在本发明的保护范围之内。