一种卫星频谱检测方法、装置、设备和存储介质与流程

本申请涉及卫星通信领域，特别是涉及一种卫星频谱检测方法、装置、设备和存储介质。

背景技术：

作为地面通信方式的有效补充，卫星移动通信已经在人们日常生活中广泛应用。而卫星移动通信系统所集中应用的频段范围恰是地面通信应用的重点区域，加之多用户多卫星同处于一个频段，因此，为了提高频谱利用率，对卫星频谱进行频谱检测是本领域技术人员需要重点研究的技术问题之一。然而，传统的频谱检测手段主要应用于地面通信网络，且在频谱检测过程中的计算复杂度较大，计算时间较长，难以满足计算资源受限的卫星以及卫星对频谱检测的实时性要求。

技术实现要素：

本申请实施例提供一种卫星频谱检测方法、装置、设备和存储介质。

第一方面，本申请实施例提供一种卫星频谱检测方法，包括：

步骤a：分别获取对当前待检测频段按照对应的采样间隔采样后得到的各个频点的能量值以及所述当前待检测频段对应的第一门限值和第二门限值，其中，所述第一门限值大于所述第二门限值；

步骤b：针对每个频点，根据所述频点的能量值、所述第一门限值和所述第二门限值的比较结果，确定下一待检测频段，并将所述下一待检测频段作为所述当前待检测频段；

步骤c：当所述第一门限值或者所述第二门限值不满足对应的预设收敛门限值时，缩小所述采样间隔，并按照缩小后的采样间隔继续执行步骤a，直至所述第一门限值和所述第二门限值均满足对应的预设收敛门限值为止；

步骤d：在所述下一待检测频段的循环平稳特征检测结果为无峰值时，确定所述下一待检测频段为空闲频段。

第二方面，本申请实施例提供一种卫星频谱检测装置，包括：

获取模块，用于分别获取对当前待检测频段按照对应的采样间隔采样后得到的各个频点的能量值以及所述当前待检测频段对应的第一门限值和第二门限值，其中，所述第一门限值大于所述第二门限值；

第一确定模块，用于针对每个频点，根据所述频点的能量值、所述第一门限值和所述第二门限值的比较结果，确定下一待检测频段，并将所述下一待检测频段作为所述当前待检测频段；

处理模块，用于当所述第一门限值或者所述第二门限值不满足对应的收敛门限值时，缩小所述采样间隔，并按照缩小后的采样间隔继续执行所述获取当前待检测频段按照对应的采样间隔采样后得到的各个频点的能量值以及所述当前待检测频段对应的第一门限值和第二门限值，直至所述第一门限值和所述第二门限值均满足对应的收敛门限值为止；

第二确定模块，用于在所述下一待检测频段的循环平稳特征检测结果为无峰值时，确定所述下一待检测频段为空闲频段。

第三方面，本申请实施例提供一种检测设备，安装于卫星中，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现本申请实施例第一方面提供的一种卫星频谱检测方法的步骤。

第四方面，本申请实施例提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现本申请实施例第一方面提供的一种卫星频谱检测方法的步骤。

本申请实施例提供的卫星频谱检测方法、装置、设备和存储介质，检测设备分别获取对当前待检测频段按照对应的采样间隔采样后得到的各个频点的能量值以及所述当前待检测频段对应的第一门限值和第二门限值，针对每个频点，根据所述频点的能量值、所述第一门限值和所述第二门限值的比较结果，从所述当前待检测频段中确定下一待检测频段，当所述第一门限值或者所述第二门限值不满足对应的预设收敛门限值时，缩小所述采样间隔，并按照缩小后的采样间隔继续对下一待检测频段进行能量检测。在第一门限值和所述第二门限值均满足对应的预设收敛门限值时，根据下一待检测频段的循环平稳特征检测结果来确定下一待检测频段的干扰情况。也就是说，在频谱检测过程中，检测设备能够基于对应的能量检测门限对当前待检测频段进行能量检测，并基于检测结果逐步缩小频段检测范围，同时，也能够根据能量检测门限是否满足预设收敛门限的收敛情况逐步缩小采样间隔，使得频谱检测适用于卫星上，既节省了卫星的计算资源，又能够快速准确判断出频段是否被干扰，具有较高的可靠性和鲁棒性。

附图说明

图1为本申请实施例提供的卫星频谱检测方法的一种流程示意图；

图2为本申请实施例提供的卫星频谱检测方法的另一种流程示意图；

图3为本申请实施例提供的卫星频谱检测方法的又一种流程示意图；

图4为本申请实施例提供的卫星频谱检测装置的结构示意图；

图5为本申请实施例提供的一种检测设备的结构示意图。

具体实施方式

为了使本申请的目的、技术方案及优点更加清楚明白，通过下述实施例并结合附图，对本申请实施例中的技术方案进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本申请，并不用于限定本申请。

需要说明的是，下述方法实施例的执行主体可以是卫星频谱检测装置，该装置可以通过软件、硬件或者软硬件结合的方式实现成为卫星中安装的检测设备的部分或者全部。下述方法实施例以执行主体是检测设备为例进行说明。

图1为本申请实施例提供的卫星频谱检测方法的一种流程示意图。本实施例涉及的是检测设备如何对卫星频谱进行干扰检测的具体过程。如图1所示，该方法可以包括：

s101、分别获取对当前待检测频段按照对应的采样间隔采样后得到的各个频点的能量值以及所述当前待检测频段对应的第一门限值和第二门限值。

其中，所述第一门限值大于所述第二门限值。当前待检测频段为本次需要进行干扰检测的待检测频段，其可以为最初始的待检测频段，也可以为最初始的待检测频段中的子频段。也就是说，当首次进行卫星频谱检测时，当前待检测频段为检测设备获取的最初始的待检测频段，反之，则当前待检测频段为最初始的待检测频段中的子频段。

同时，针对当前待检测频段可以设置对应的采样间隔，在当前待检测频段为最初始的待检测频段中的子频段时，所对应的第一采样间隔小于在当前待检测频段为最初始的待检测频段时的第二采样间隔，即在当前待检测频段逐步缩小时，所对应的采样间隔也逐步缩小，从而提高了频谱分辨率，进而基于更高分辨率的频谱进行干扰检测，可以提高频谱检测的准确性。可选的，采样间隔大于卫星的最小信道带宽。

进一步的，能量检测门限可以包括第一门限值和第二门限值，第一门限值大于第二门限值，且能量检测门限与采样后得到的各个频点的能量值相关。同时，能量检测门限也与当前待检测频段相对应，在当前待检测频段发生变化时，所对应的能量检测门限也随之发生变化，即本申请实施例所提供的方案，检测设备能够在干扰信号不确定的情况下自适应改变能量检测门限，使其不断逼近预设收敛门限。

在实际应用中，检测设备按照当前待检测频段对应的采样间隔对当前待检测频段进行均匀采样，得到n′个离散的时域信号。其中，n′等于当前待检测频段的带宽与采样间隔之比。接着，检测设备对得到的时域信号进行傅里叶变换，得到n′个频域信号，即得到n′个频点。可选的，为了提高检测结果的准确性，可以对当前待检测频段进行多次采样，并累加计算上述得到的n′个频点中的每个频点的能量值。其中，每个频点的能量值的获取过程可以为：假设对当前待检测频段进行第i次傅里叶变换后的结果为xi＝(xi,1,xi,2,…,xi,n)，则n′个频点中的第k个频点的能量值ek可以通过下述公式1或者公式1的变型计算得到：

公式1：

其中，nt为采样次数，xj,k为第k个频点在第j次采样时的能量值。

接着，在得到各个频点的能量值之后，检测设备便可以基于各个频点的能量值确定当前待检测频段对应的能量检测门限，即基于各个频点的能量值确定第一门限值和第二门限值。

s102、针对每个频点，根据所述频点的能量值、所述第一门限值和所述第二门限值的比较结果，从所述当前待检测频段中确定下一待检测频段，并将所述下一待检测频段作为所述当前待检测频段。

其中，下一待检测频段为当前待检测频段的子频段。在得到各个频点的能量值以及当前待检测频段对应的能量检测门限之后，针对每个频点，检测设备将频点的能量值分别于第一门限值和第二门限值进行比较，并基于比较结果从当前待检测频段中确定下一待检测频段。具体的比较过程可以为以下几种情况：

(1)若频点的能量值大于第一门限值，则将该频点和该频点的下一频点之间的频段确定为被干扰频段；

(2)若频点的能量值小于第二门限值，则将该频点和该频点的下一频点之间的频段确定为空闲频段；

(3)若频点的能量值大于第二门限值且小于第一门限值，则确定该频点和该频点的下一频点之间的频段为下一待检测频段。

s103、当所述第一门限值或者所述第二门限值不满足对应的预设收敛门限值时，缩小所述采样间隔，并按照缩小后的采样间隔继续执行步骤s101，直至所述第一门限值和所述第二门限值均满足对应的预设收敛门限值为止。

其中，检测设备可以根据预设的虚警率pf和检测率pd，确定频谱检测过程中所需要的收敛门限值。在计算预设收敛门限值时，设二元假设模型表示为：

公式2：

其中，n为检测样本数，w(n)为噪声信号，x(n)为通信信号，r(n)为卫星上的接收机的接收信号，h0为接收信号中不存在通信信号，h1为接收信号中存在通信信号，r(n)的能量值y近似符合高斯分布，则在当前待检测频段存在干扰时，虚警率pf和检测率pd的表达式为：

公式3：

公式4：

其中，为高斯白噪声的方差，为信号的平均功率，当当前待检测频段实际存在干扰时，虚警率pf为未检测出干扰的概率，检测率pd为检测出干扰的概率。这样，检测设备便可以根据设定的pf和pd，通过上述公式3和公式4，或者公式3和公式4的变型确定出第一门限值对应的预设收敛门限ηh，以及第二门限值对应的预设收敛门限ηl。

接着，检测设备比较第一门限值和第一门限值对应的预设收敛门限，以及比较第二门限值和第二门限值对应的预设收敛门限，若第一门限值大于第一门限值对应的预设收敛门限，或者第二门限值小于第二门限值对应的预设收敛门限，则缩小上述采样间隔，并按照缩小后的采样间隔对下一待检测频段进行更高分辨率的均匀采样，并继续执行上述s101-s103的过程，直至第一门限值和第二门限值均满足对应的预设收敛门限值为止，即直至第一门限值和第二门限值无限逼近对应的预设收敛门限值为止。

s104、在所述下一待检测频段的循环平稳特征检测结果为无峰值时，确定所述下一待检测频段为空闲频段。

其中，在第一门限值和第二门限值均满足对应的预设收敛门限值时，检测设备对下一待检测频段进行循环平稳特征检测，在下一待检测频段的循环平稳特征检测结果为无峰值时，确定所述下一待检测频段为空闲频段。可选的，在下一待检测频段的循环平稳特征检测结果为有峰值时，确定所述下一待检测频段为被干扰频段。

具体的循环平稳特征检测过程可以为：检测设备计算下一待检测频段中的接收信号的循环功率谱密度，并判断接收信号的循环功率谱密度在非零循环频率处是否有峰值，若在非零循环频率处有峰值，则确定下一待检测频段为被干扰频段；若在非零循环频率处无峰值，则确定下一待检测频段为空闲频段。

本申请实施例提供的卫星频谱检测方法，检测设备分别获取当前待检测频段按照对应的采样间隔采样后得到的各个频点的能量值以及所述当前待检测频段对应的第一门限值和第二门限值，针对每个频点，根据所述频点的能量值、所述第一门限值和所述第二门限值的比较结果，从所述当前待检测频段中确定下一待检测频段，当所述第一门限值或者所述第二门限值不满足对应的预设收敛门限值时，缩小所述采样间隔，并按照缩小后的采样间隔继续对下一待检测频段进行能量检测。在第一门限值和所述第二门限值均满足对应的预设收敛门限值时，根据下一待检测频段的循环平稳特征检测结果来确定下一待检测频段的干扰情况。也就是说，在频谱检测过程中，检测设备能够基于对应的能量检测门限对当前待检测频段进行能量检测，并基于检测结果逐步缩小频段检测范围，同时，也能够根据能量检测门限是否满足预设收敛门限的收敛情况逐步缩小采样间隔，使得频谱检测适用于卫星上，既节省了卫星的计算资源，又能够快速准确判断出频段是否被干扰，具有较高的可靠性和鲁棒性。

在一个实施例中，还提供了当前待检测频段对应的第一门限值和第二门限值的获取过程。在上述实施例的基础上，可选的，如图2所示，上述检测设备获取所述当前待检测频段对应的第一门限值和第二门限值的过程可以为：

s201、按照从大到小的顺序，对各个频点的能量值进行排序。

其中，为了进一步提高检测结果的准确性，可以从当前待检测频段中选择更多的采样样本进行分析，为此，可选的，在上述s201之前，检测设备还可以对所有频点中的相邻两个频点进行l倍的插值估计，其中，所述l为自然数。

以对所有频点中的第i个频点和第i+1个频点进行l倍的插值估计为例，l倍插值是指在第i个频点和第i+1个频点之间插入l个频点，所插入的每个频点的能量值e插值频点按照下述公式5或者公式5的变型计算得到：

公式5：

其中，ei为第i个频点的能量值，ei+1为第i+1个频点的能量值。

接着，检测设备按照从大到小的顺序，对各个频点的能量值进行排序，得到的排序结果为{e1,e2…,en}。

s202、根据排序结果中的最大预设数量个频点的能量值，确定所述当前待检测频段对应的第一门限值。

其中，该预设数量可以根据实际应用需求进行相应的设置。可选的，可以将该预设数量设置为n′*m％，m为小于100的自然数。以该预设数量为n′*m％为例，检测设备可以根据下述公式6或者公式6的变型，确定当前待检测频段对应的第一门限值λh：

公式6：

s203、根据所述排序结果中的最小预设数量个频点的能量值，确定所述当前待检测频段对应的第二门限值。

其中，继续以该预设数量为n′*m％为例，检测设备可以根据下述公式7或者公式7的变型，确定当前待检测频段对应的第二门限值λl：

公式7：

综上，可以看出，在当前待检测频段发生变化时，当前待检测频段对应的能量检测门限(即第一门限值和第二门限值)也会发生相应的变化，并随着当前待检测频段以及采样间隔得不断缩小，能量检测门限不断逼近对应的预设收敛门限，从而提高了检测结果的准确性。

在本实施例中，检测设备可以根据对当前待检测频段进行采样后得到的各个频点的能量值的排序结果，来确定当前待检测频段对应的能量检测门限，使得能量检测门限与当前待检测频段相对应，从而在不断缩小当前待检测频段的范围时，也能自适应调整对应的能量检测门限，使能量检测门限不断逼近对应的预设收敛门限，这样，基于无限逼近的能量检测门限得到的检测结果更加准确，进一步提高了检测结果的准确性。

在上述实施例的基础上，可选的，在对下一待检测频段进行循环平稳特征检测之前，检测设备还可以对下一待检测频段进行扫描式或者串行式频谱感知，以获得更高分辨率的频谱。

对应地，上述s104可以包括：对所述下一待检测频段进行扫描式或者串行式频谱感知，得到新的下一待检测频段；在所述新的下一待检测频段的循环平稳特征检测结果为无峰值时，确定所述新的下一待检测频段为空闲频段。其中，所述新的下一待检测频段的分辨率高于所述下一待检测频段，且新的下一待检测频段与下一待检测频段的频段范围相同。

在本实施例中，检测设备对下一待检测频段进行扫描式或者串行式频谱感知，以获得更高分辨率的频谱，进而对更高分辨率的频谱进行循环平稳特征检测，可以得到更加准确的干扰检测结果，进一步提高了检测结果的准确性。

为了便于本领域技术人员的理解，以下结合图3具体介绍检测设备进行卫星频谱检测的过程，具体的，如图3所示，该方法可以包括：

s301、检测设备按照对应的采样间隔对当前待检测频段进行均匀采样。

其中，采样间隔fc与卫星的最小信道带宽fmin之间的关系如下：fc＝α·fmin，α为扩展系数，α的取值可以为(1，2]。经过采样后，检测设备得到n′个离散的时域信号。其中，n′等于当前待检测频段的带宽与采样间隔之比。另外，可以对当前待检测频段进行多次采样。

s302、对采样得到的时域信号进行傅里叶变换，得到多个(n′个)频点，并基于多次采样结果分别计算每个频点的能量值。

s303、对多个(n′个)频点中的每相邻两个频点进行l倍的插值估计。

s304、按照从大到小的顺序，对各个频点的能量值进行排序，并基于排序结果确定当前待检测频段对应的第一门限值和第二门限值。

s305、确定频谱检测过程中所需要的收敛门限值。

s306、针对每个频点，根据频点的能量值、第一门限值和第二门限值的比较结果，从当前待检测频段中确定下一待检测频段。

s307、判断第一门限值和第二门限值是否满足对应的预设收敛门限值。

其中，若第一门限值或者第二门限值不满足对应的预设收敛门限值，则执行s308，若第一门限值和第二门限值均满足对应的预设收敛门限值，则执行s309-s310。

s308、缩小采样间隔，并将下一待检测频段作为当前待检测频段，以及按照缩小后的采样间隔继续执行步骤s301，直至第一门限值和第二门限值均满足对应的预设收敛门限值为止。

s309、对下一待检测频段进行扫描式或者串行式频谱感知，得到新的下一待检测频段。

其中，所述新的下一待检测频段的分辨率高于所述下一待检测频段，且新的下一待检测频段与下一待检测频段的频段范围相同。

s310、对新的下一待检测频段进行循环平稳特征检测。

其中，在新的下一待检测频段的循环平稳特征检测结果为无峰值时，确定该新的下一待检测频段为空闲频段；在新的下一待检测频段的循环平稳特征检测结果为有峰值时，确定该新的下一待检测频段为被干扰频段。

图4为本申请实施例提供的卫星频谱检测装置的结构示意图。如图4所示，该装置可以包括：获取模块10、第一确定模块11、处理模块12和第二确定模块13。

具体的，获取模块10用于分别获取对当前待检测频段按照对应的采样间隔采样后得到的各个频点的能量值以及所述当前待检测频段对应的第一门限值和第二门限值，其中，所述第一门限值大于所述第二门限值；

第一确定模块11用于针对每个频点，根据所述频点的能量值、所述第一门限值和所述第二门限值的比较结果，从所述当前待检测频段中确定下一待检测频段，并将所述下一待检测频段作为所述当前待检测频段；

处理模块12用于当所述第一门限值或者所述第二门限值不满足对应的收敛门限值时，缩小所述采样间隔，并按照缩小后的采样间隔继续执行所述获取对当前待检测频段按照对应的采样间隔采样后得到的各个频点的能量值以及所述当前待检测频段对应的第一门限值和第二门限值，直至所述第一门限值和所述第二门限值均满足对应的收敛门限值为止；

第二确定模块13用于在所述下一待检测频段的循环平稳特征检测结果为无峰值时，确定所述下一待检测频段为空闲频段。

本申请实施例提供的卫星频谱检测装置，检测设备分别获取当前待检测频段按照对应的采样间隔采样后得到的各个频点的能量值以及所述当前待检测频段对应的第一门限值和第二门限值，针对每个频点，根据所述频点的能量值、所述第一门限值和所述第二门限值的比较结果，从所述当前待检测频段中确定下一待检测频段，当所述第一门限值或者所述第二门限值不满足对应的预设收敛门限值时，缩小所述采样间隔，并按照缩小后的采样间隔继续对下一待检测频段进行能量检测。在第一门限值和所述第二门限值均满足对应的预设收敛门限值时，根据下一待检测频段的循环平稳特征检测结果来确定下一待检测频段的干扰情况。也就是说，在频谱检测过程中，检测设备能够基于对应的能量检测门限对当前待检测频段进行能量检测，并基于检测结果逐步缩小频段检测范围，同时，也能够根据能量检测门限是否满足预设收敛门限的收敛情况逐步缩小采样间隔，使得频谱检测适用于卫星上，既节省了卫星的计算资源，又能够快速准确判断出频段是否被干扰，具有较高的可靠性和鲁棒性。

在上述实施例的基础上，可选的，获取模块10包括：排序单元、第一确定单元和第二确定单元；

具体的，排序单元用于按照从大到小的顺序，对各个频点的能量值进行排序；

第一确定单元用于根据排序结果中的最大预设数量个频点的能量值，确定所述当前待检测频段对应的第一门限值；

第二确定单元用于根据所述排序结果中的最小预设数量个频点的能量值，确定所述当前待检测频段对应的第二门限值。

在上述实施例的基础上，可选的，获取模块10还包括：插值单元；

具体的，插值单元用于在排序单元按照从大到小的顺序，对各个频点的能量值进行排序之前，对所有频点中的相邻两个频点进行l倍的插值估计，其中，所述l为自然数。

在上述实施例的基础上，可选的，第一确定模块11具体用于若所述频点的能量值大于所述第二门限值且小于所述第一门限值，则确定所述频点和所述频点的下一频点之间的频段为下一待检测频段。

在上述实施例的基础上，可选的，第二确定模块13具体用于对所述下一待检测频段进行扫描式或者串行式频谱感知，得到新的下一待检测频段；在所述新的下一待检测频段的循环平稳特征检测结果为无峰值时，确定所述新的下一待检测频段为空闲频段，其中，所述新的下一待检测频段的分辨率高于所述下一待检测频段。

在上述实施例的基础上，可选的，第二确定模块13还用于在所述下一待检测频段的循环平稳特征检测结果为有峰值时，确定所述下一待检测频段为被干扰频段。

可选的，所述采样间隔大于卫星的最小信道带宽。

在一个实施例中，提供了一种检测设备，其结构示意图可以如图5所示。该检测设备包括通过系统总线连接的处理器、存储器。其中，该检测设备的处理器用于提供计算和控制能力。该检测设备的存储器用于存储计算机程序。该计算机程序被处理器执行时以实现一种卫星频谱检测方法。

本领域技术人员可以理解，图5中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的检测设备的限定，具体的检测设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。

在一个实施例中，提供了一种检测设备，安装于卫星中，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：

步骤a：分别获取对当前待检测频段按照对应的采样间隔采样后得到的各个频点的能量值以及所述当前待检测频段对应的第一门限值和第二门限值，其中，所述第一门限值大于所述第二门限值；

步骤b：针对每个频点，根据所述频点的能量值、所述第一门限值和所述第二门限值的比较结果，从所述当前待检测频段中确定下一待检测频段，并将所述下一待检测频段作为所述当前待检测频段；

步骤c：当所述第一门限值或者所述第二门限值不满足对应的预设收敛门限值时，缩小所述采样间隔，并按照缩小后的采样间隔继续执行步骤a，直至所述第一门限值和所述第二门限值均满足对应的预设收敛门限值为止；

步骤d：在所述下一待检测频段的循环平稳特征检测结果为无峰值时，确定所述下一待检测频段为空闲频段。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：按照从大到小的顺序，对各个频点的能量值进行排序；根据排序结果中的最大预设数量个频点的能量值，确定所述当前待检测频段对应的第一门限值；根据所述排序结果中的最小预设数量个频点的能量值，确定所述当前待检测频段对应的第二门限值。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：对所有频点中的相邻两个频点进行l倍的插值估计，其中，所述l为自然数。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：若所述频点的能量值大于所述第二门限值且小于所述第一门限值，则确定所述频点和所述频点的下一频点之间的频段为下一待检测频段。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：对所述下一待检测频段进行扫描式或者串行式频谱感知，得到新的下一待检测频段；在所述新的下一待检测频段的循环平稳特征检测结果为无峰值时，确定所述新的下一待检测频段为空闲频段，其中，所述新的下一待检测频段的分辨率高于所述下一待检测频段。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：在所述下一待检测频段的循环平稳特征检测结果为有峰值时，确定所述下一待检测频段为被干扰频段。

可选的，所述采样间隔大于卫星的最小信道带宽。

在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：

步骤a：分别获取对当前待检测频段按照对应的采样间隔采样后得到的各个频点的能量值以及所述当前待检测频段对应的第一门限值和第二门限值，其中，所述第一门限值大于所述第二门限值；

步骤b：针对每个频点，根据所述频点的能量值、所述第一门限值和所述第二门限值的比较结果，从所述当前待检测频段中确定下一待检测频段，并将所述下一待检测频段作为所述当前待检测频段；

步骤c：当所述第一门限值或者所述第二门限值不满足对应的预设收敛门限值时，缩小所述采样间隔，并按照缩小后的采样间隔继续执行步骤a，直至所述第一门限值和所述第二门限值均满足对应的预设收敛门限值为止；

步骤d：在所述下一待检测频段的循环平稳特征检测结果为无峰值时，确定所述下一待检测频段为空闲频段。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：按照从大到小的顺序，对各个频点的能量值进行排序；根据排序结果中的最大预设数量个频点的能量值，确定所述当前待检测频段对应的第一门限值；根据所述排序结果中的最小预设数量个频点的能量值，确定所述当前待检测频段对应的第二门限值。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：对所有频点中的相邻两个频点进行l倍的插值估计，其中，所述l为自然数。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：若所述频点的能量值大于所述第二门限值且小于所述第一门限值，则确定所述频点和所述频点的下一频点之间的频段为下一待检测频段。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：对所述下一待检测频段进行扫描式或者串行式频谱感知，得到新的下一待检测频段；在所述新的下一待检测频段的循环平稳特征检测结果为无峰值时，确定所述新的下一待检测频段为空闲频段，其中，所述新的下一待检测频段的分辨率高于所述下一待检测频段。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：在所述下一待检测频段的循环平稳特征检测结果为有峰值时，确定所述下一待检测频段为被干扰频段。

可选的，所述采样间隔大于卫星的最小信道带宽。

上述实施例中提供的卫星频谱检测装置、设备以及存储介质可执行本申请任意实施例所提供的卫星频谱检测方法，具备执行该方法相应的功能模块和有益效果。未在上述实施例中详尽描述的技术细节，可参见本申请任意实施例所提供的卫星频谱检测方法。

以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本申请专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。