CN115078847A - 深空探测微波载荷的干扰源确定方法、装置、设备及介质 - Google Patents

Abstract

本公开提供一种深空探测微波载荷的干扰源确定方法、装置、设备及介质，方法包括：计算不同极化天线在同一时刻接收的微波信号之间的时域相关系数；计算不同极化天线接收的同一频点的微波信号之间的频域相关系数；计算同一极化天线在不同时刻接收的微波信号之间的交叉相关系数；根据时域相关系数和/或频域相关系数和/或交叉相关系数确定干扰信号的来源和频率。该方法能够实现对微波干扰来源直观、可靠的分析，解决了由于干扰情况复杂难以确定干扰来源的技术问题。

Description

深空探测微波载荷的干扰源确定方法、装置、设备及介质

技术领域

本公开涉及微波信号处理技术领域，尤其涉及一种深空探测微波载荷的干扰源确定方法、装置、设备及介质。

背景技术

为探测目标天体表面、次表层特性及空间电磁环境，深空探测器上通常会搭载一些能进行微波探测的仪器，如雷达、低频射电频谱仪等，这些仪器被广泛应用于月球与火星探测中。但实际探测经验表明，此类仪器对微波极为敏感，容易受到电磁干扰的影响而降低探测数据质量，而探测器本身以及其他载荷等原因使得星上电磁干扰情况极其复杂，难以确定电磁干扰的来源，对微波载荷数据的干扰定位与消除都造成一定困难。

发明内容

鉴于上述技术问题，本公开第一方面提供一种深空探测微波载荷的干扰源确定方法，包括：计算不同极化天线在同一时刻接收的微波信号之间的时域相关系数；计算不同极化天线接收的同一频点的微波信号之间的频域相关系数；计算同一极化天线在不同时刻接收的微波信号之间的交叉相关系数；根据时域相关系数和/或频域相关系数和/或交叉相关系数确定干扰信号的来源和频率。

根据本公开的实施例，计算不同极化天线在同一时刻接收的微波信号之间的时域相关系数，具体包括：根据

计算时域相关系数r(X，Y)，其中，X，Y分别表示不同极化天线在同一时刻接收的微波信号，Cov(X，Y)为X，Y之间的协方差，Var[X]为X的方差，Var[Y]为Y的方差。

根据本公开的实施例，计算不同极化天线接收的同一频点的微波信号之间的频域相关系数，具体包括：计算每一微波数据的频谱；将同一极化天线接收到的微波数据的频谱按照时间排序，得到每个频点对应的频谱数据；将不同极化天线同一频点的频谱数据或复数频谱数据进行相关性计算，得到频域相关系数。

根据本公开的实施例，计算同一极化天线在不同时刻接收的微波信号之间的交叉相关系数，具体包括：从微波信号中确定一组具有代表性的微波信号；对具有代表性的微波信号进行循环圆周移动，每移动一个时间采样点产生一个微波信号序列；计算每个微波信号序列与相同极化天线接收到的微波信号的相关系数，得到交叉相关系数。

根据本公开的实施例，根据时域相关系数和/或频域相关系数和/或交叉相关系数确定干扰信号的来源和频率，具体包括：响应于时域相关系数小于第一阈值，确定得到干扰信号来源于天线内部；响应于时域相关系数大于第一阈值且小于第二阈值，确定得到干扰信号来源于天线内部和天线外部；响应于时域相关系数大于第二阈值，确定得到干扰信号来源于天线外部。

根据本公开的实施例，根据时域相关系数和/或频域相关系数和/或交叉相关系数确定干扰信号的来源和频率，具体包括：根据频域相关系数确定干扰信号的频率；或根据交叉相关系数估计微波信号的重复周期，并根据重复周期估计干扰信号的频率。

根据本公开的实施例，计算每一微波数据的频谱，具体包括：对微波数据进行傅里叶变换，然后计算傅里叶变换结果的幅度，得到微波数据的频谱；或直接采用傅里叶变换的复数结果作为所述微波数据的频谱。

本公开第二方面提供一种深空探测微波载荷的干扰源确定装置，包括：第一计算模块，用于计算不同极化天线在同一时刻接收的微波信号之间的时域相关系数；第二计算模块，用于计算不同极化天线接收的同一频点的微波信号之间的频域相关系数；第三计算模块，用于计算同一极化天线在不同时刻接收的微波信号之间的交叉相关系数；确定模块，用于根据时域相关系数和/或频域相关系数和/或交叉相关系数确定干扰信号的来源和频率。

本公开第三方面还提供一种电子设备，包括：一个或多个处理器；存储器，用于存储一个或多个程序，其中，当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现上述所述的方法。

本公开第四方面还提供一种计算机可读存储介质，其特征在于，其上存储有可执行指令，该指令被处理器执行时使处理器实现上述所述的方法。

根据本公开实施例提供的深空探测微波载荷的干扰源确定方法、装置、设备及介质，至少能够实现以下技术效果：

通过对微波数据进行时域相关性分析和频域相关性分析来确定干扰信号的来源和频率，相比起直接对比微波载荷不同极化的时域数据以及频谱来评估干扰，使用相关方法来进行干扰评估更为直观、可靠，解决了高频电磁波的波形复杂，难以进行比较的问题。

进一步地，采用频域相关性分析可以确定哪些频率的频谱幅度随时间变化的趋势相关性较强，除了有助于分析干扰来源，这可能为确定噪声的频率位置提供更可靠的参考。

更进一步地，交叉相关系数的求取与应用可以确定信号中是否存在某些分量重复出现，或者具有较强的相关性，有助于定位噪声波形以及重复周期，辅助进行噪声源头确定与噪声压制。

附图说明

通过以下参照附图对本公开实施例的描述，本公开的上述以及其他目的、特征和优点将更为清楚，在附图中：

图1示意性示出了根据本公开实施例的深空探测微波载荷的干扰源确定方法的处理方法图。

图2示意性示出了根据本公开一实施例的频域相关系数计算方法流程图。

图3示意性示出了根据本公开一实施例的交叉相关系数计算方法流程图。

图4示意性示出了根据本公开一实施例的干扰信号的来源确定方法流程图。

图5示意性示出了根据本公开实施例的低频射电频谱仪干扰信号的处理装置的框图。

图6示意性示出了根据本发明实施例的适于实现上文描述的方法的电子设备的框图。

具体实施方式

为使本公开的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本公开进一步详细说明。显然，所描述的实施例是本公开一部分实施例，而不是全部的实施例。基于本公开中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本公开保护的范围。

在此使用的术语仅仅是为了描述具体实施例，而并非意在限制本公开。在此使用的术语“包括”、“包含”等表明了所述特征、步骤、操作和/或部件的存在，但是并不排除存在或添加一个或多个其他特征、步骤、操作或部件。

在本公开中，除非另有明确的规定和限定，术语“安装”、“相连”“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接或可以互相通讯；可以是直接连接，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本公开中的具体含义。

在本公开的描述中，需要理解的是，术语“纵向”、“长度”、“周向”、“前”、“后”、“左”、“右”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本公开和简化描述，而不是指示或暗示所指的子系统或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本公开的限制。

贯穿附图，相同的元素由相同或相近的附图标记来表示。可能导致本公开的理解造成混淆时，将省略常规结构或构造。并且图中各部件的形状、尺寸、位置关系不反映真实大小、比例和实际位置关系。另外，在本公开中，不应将位于括号之间的任何参考符号构造成对本公开的限制。

类似地，为了精简本公开并帮助理解各个公开方面中的一个或多个，在上面对本公开示例性实施例的描述中，本公开的各个特征有时被一起分到单个实施例、图或者对其描述中。参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本公开的至少一个实施例或示例中。本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或者多个实施例或示例中以合适的方式结合。

此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。因此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本公开的描述中，“多个”的含义是至少两个，例如两个、三个等，除非另有明确具体的限定。

针对现有技术的不足，本公开实施例提供一种深空探测微波载荷的干扰源确定方法，借助微波载荷上不同极化方式、指向的天线同一时刻接收到的信号以及不同极化天线接收的同一频点的信号，同一极化天线接收的不同时刻的信号，基于相关方法实现对微波干扰来源直观、可靠的分析。下面结合具体的实施例进行详细介绍。

图1示意性示出了根据本公开实施例的深空探测微波载荷的干扰源确定方法的处理方法图。

如图1所示，该深空探测微波载荷的干扰源确定方法例如可以包括操作S101～操作S104。

在操作S101，计算不同极化天线在同一时刻接收的微波信号之间的时域相关系数。

一般情况下，深空探测微波载荷会安装多个不同极化方式的天线，以接收电磁波中不同极化方向的分量，而微波干扰如果来自探测器外部及空间环境，则不同天线同时接收到的电磁波数据应显示出较高的相关性。因此，在本实施例中，可以计算不同极化天线在同一时刻接收的微波信号之间的时域相关系数来确定干扰信号的来源。

在本实施例中，可以通过计算不同极化天线接收到的微波信号的皮尔逊相关系数得到时域相关系数。具体地，可以根据

计算时域相关系数r(X，Y)，其中，X，Y分别表示不同极化天线在同一时刻接收的微波信号，Cov(X，Y)为X，Y之间的协方差，Var[X]为X的方差，Var[Y]为Y的方差。

如果接收到的微波信号下传后为复数，则可以直接将复数的微波信号输入上述公式，输出的r(X，Y)也为复数，取实部为最终求得的时域相关系数。

在操作S102，计算不同极化天线接收的同一频点的微波信号之间的频域相关系数。

图2示意性示出了根据本公开一实施例的频域相关系数计算方法流程图。

如图2所示，频域相关系数的计算过程例如可以包括操作S201～操作S203。

在操作S201，计算每一微波数据的频谱。

在本公开一实施例中，可以通过对微波数据进行傅里叶变换，并计算傅里叶变换结果的幅度，得到微波数据的频谱。或直接采用傅里叶变换的复数结果作为微波数据的频谱，则相关系数计算后需取实部得到最终结果。如果微波信号在下传前经过IQ解调去除载频，则可以进行频谱移动，将原来频段为[0，Fs]的频谱转换为频段为[-Fs/2，Fs/2]的频谱，将中心频率置于0。

在操作S202，将同一极化天线接收到的微波数据的频谱按照时间排序，得到每个频点对应的频谱数据。

在本公开一实施例中，可以将同一天线接收到的数据的频谱按照接收时间依次排列，则每个频率点都可以获得一条长度等于道数的频谱数据，代表了该频点能量随时间变化的趋势。

在操作S203，将不同极化天线同一频点的频谱数据进行相关性计算，得到频域相关系数。

在本公开一实施例中，可以将不同天线接收，同一频点的频谱数据进行频域相关系数计算。如使用复数频谱数据进行相关性计算，则计算结果需取实部得到频域相关系数。该参数代表不同极化方向下该频率能量随时间变化趋势的相关性，最终每个频率采样点都能计算出一个频域相关系数。

在操作S103，计算同一极化天线在不同时刻接收的微波信号之间的交叉相关系数。

图3示意性示出了根据本公开一实施例的交叉相关系数计算方法流程图。

如图3所示，交叉相关系数的计算过程例如可以包括操作S301～操作S303。

在操作S301，从微波信号中确定一组具有代表性的微波信号。

在本公开一实施例中，代表性可以指幅度较大，频谱结构与其他信号较为相似等，具体可以根据实际应用需求，本公开不做限制。

在操作S302，对具有代表性的微波信号进行循环圆周移动，每移动一个时间采样点产生一个微波信号序列。

在本公开一实施例中，可以对具有代表性的微波信号进行循环圆周移动每移动一个时间采样点，原始微波信号采样点发生变化，从而产生一个新的信号序列。

在操作S303，计算每个微波信号序列与相同极化天线接收到的微波信号的相关系数，得到交叉相关系数。

在本公开一实施例中，计算每个信号序列与其他相同天线接收到的信号的相关系数，可以取相关系数的实部，这样就能产生一条时间采样点移动数与相关系数实部的对应曲线，每道其他信号都能产生这样一条曲线。

在操作S104，根据时域相关系数和/或频域相关系数和/或交叉相关系数确定干扰信号的来源和频率。

图4示意性示出了根据本公开一实施例的干扰信号的来源确定方法流程图。

如图4所示，干扰信号的来源确定方法例如可以包括操作S401～操作S404。

在操作S401，判断时域相关系数的值与预设阈值的大小。

在操作S402，响应于时域相关系数小于第一阈值，确定得到干扰信号来源于天线内部。

在操作S403，响应于时域相关系数大于第一阈值且小于第二阈值，确定得到干扰信号来源于天线内部和天线外部。

在操作S404，响应于时域相关系数大于第二阈值，确定得到干扰信号来源于天线外部。

示例性地，第一阈值可以为0.2，第二阈值可以为0.6。如果相关系数小于等于0.2，则说明X，Y数据间基本不具有相关性，干扰来源为外界而被天线接收的可能性很低；如果相关系数大于0.2小于0.6，则显示弱相关性，可能存在内外多种干扰来源且强度差异不大；如果相关系数大于等于0.6，则说明相关性较强，干扰可能来自于探测器外部，统一被天线所接收。

应当理解，第一阈值和第二阈值的大小可以根据实际应用需求确定，本公开不做限制。

在本公开另一实施例中，干扰信号的频率的确定方法例如可以包括：根据频域相关系数确定干扰信号的频率。具体地，根据不同频点的频域相关系数，可以初步确定出干扰信号所在的频率，便于后续的干扰抑制。

在本公开另一实施例中，干扰信号的频率的确定方法例如还可以包括：根据交叉相关系数估计微波信号的重复周期，并根据重复周期估计干扰信号的频率。具体的，在操作S303中获取了时间采样点移动数与相关系数实部的对应曲线，此时可以寻找曲线中较为突出的尖峰或梳状结构，以此为依据估计信号的重复周期，并推测对应干扰的频率。

综上所述，本公开实施例提供的深空探测微波载荷的干扰源确定方法，通过对微波数据进行时域相关性分析和频域相关性分析来确定干扰信号的来源和频率，相比起直接对比微波载荷不同极化的时域数据以及频谱来评估干扰，使用相关方法来进行干扰评估更为直观、可靠，解决了高频电磁波的波形复杂，难以进行比较的问题。采用频域相关性分析可以确定哪些频率的频谱幅度随时间变化的趋势相关性较强，除了有助于分析干扰来源，这可能为确定噪声的频率位置提供更可靠的参考。交叉相关系数的求取与应用可以确定信号中是否存在某些分量重复出现，或者具有较强的相关性，有助于定位噪声波形以及重复周期，辅助进行噪声源头确定与噪声压制。

图5示意性示出了根据本公开实施例的低频射电频谱仪干扰信号的处理装置的框图。

如图5所示，深空探测微波载荷的干扰源确定装置500可以包括第一计算模块510、第二计算模块520、第三计算模块530以及确定模块540。

第一计算模块510，用于计算不同极化天线在同一时刻接收的微波信号之间的时域相关系数。

第二计算模块520，用于计算不同极化天线接收的同一频点的微波信号之间的频域相关系数。

第三计算模块530，用于计算同一极化天线在不同时刻接收的微波信号之间的交叉相关系数。

确定模块540，用于根据时域相关系数和/或频域相关系数和/或交叉相关系数确定干扰信号的来源和频率。

根据本发明的实施例的模块、子模块、单元、子单元中的任意多个、或其中任意多个的至少部分功能可以在一个模块中实现。根据本发明实施例的模块、子模块、单元、子单元中的任意一个或多个可以被拆分成多个模块来实现。根据本发明实施例的模块、子模块、单元、子单元中的任意一个或多个可以至少被部分地实现为硬件电路，例如现场可编程门阵列(FPGA)、可编程逻辑阵列(PLA)、片上系统、基板上的系统、封装上的系统、专用集成电路(ASIC)，或可以通过对电路进行集成或封装的任何其他的合理方式的硬件或固件来实现，或以软件、硬件以及固件三种实现方式中任意一种或以其中任意几种的适当组合来实现。或者，根据本发明实施例的模块、子模块、单元、子单元中的一个或多个可以至少被部分地实现为计算机程序模块，当该计算机程序模块被运行时，可以执行相应的功能。

例如，第一计算模块510、第二计算模块520、第三计算模块530以及确定模块540中的任意多个可以合并在一个模块/单元/子单元中实现，或者其中的任意一个模块/单元/子单元可以被拆分成多个模块/单元/子单元。或者，这些模块/单元/子单元中的一个或多个模块/单元/子单元的至少部分功能可以与其他模块/单元/子单元的至少部分功能相结合，并在一个模块/单元/子单元中实现。根据本发明的实施例第一计算模块510、第二计算模块520、第三计算模块530以及确定模块540中的至少一个可以至少被部分地实现为硬件电路，例如现场可编程门阵列(FPGA)、可编程逻辑阵列(PLA)、片上系统、基板上的系统、封装上的系统、专用集成电路(ASIC)，或可以通过对电路进行集成或封装的任何其他的合理方式等硬件或固件来实现，或以软件、硬件以及固件三种实现方式中任意一种或以其中任意几种的适当组合来实现。或者，第一计算模块510、第二计算模块520、第三计算模块530以及确定模块540中的至少一个可以至少被部分地实现为计算机程序模块，当该计算机程序模块被运行时，可以执行相应的功能。

需要说明的是，本发明的实施例的深空探测微波载荷的干扰源确定装置部分与本发明的实施例的深空探测微波载荷的干扰源确定方法部分是相对应的，其具体实施细节及带来的技术效果也是相同的，在此不再赘述。

图6示意性示出了根据本发明实施例的适于实现上文描述的方法的电子设备的框图。图6示出的电子设备仅仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。

如图6所示，根据本发明实施例的电子设备600包括处理器601，其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM)603中的程序而执行各种适当的动作和处理。处理器601例如可以包括通用微处理器(例如CPU)、指令集处理器和/或相关芯片组和/或专用微处理器(例如，专用集成电路(ASIC))，等等。处理器601还可以包括用于缓存用途的板载存储器。处理器601可以包括用于执行根据本发明实施例的方法流程的不同动作的单一处理单元或者是多个处理单元。

在RAM 603中，存储有电子设备600操作所需的各种程序和数据。处理器601、ROM602以及RAM603通过总线604彼此相连。处理器601通过执行ROM 602和/或RAM603中的程序来执行根据本发明实施例的方法流程的各种操作。需要注意，所述程序也可以存储在除ROM602和RAM 603以外的一个或多个存储器中。处理器601也可以通过执行存储在所述一个或多个存储器中的程序来执行根据本发明实施例的方法流程的各种操作。

根据本发明的实施例，电子设备600还可以包括输入/输出(I/O)接口605，输入/输出(I/O)接口605也连接至总线604。电子设备600还可以包括连接至I/O接口605的以下部件中的一项或多项：包括键盘、鼠标等的输入部分606；包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分607；包括硬盘等的存储部分608；以及包括诸如LAN卡、调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口605。可拆卸介质611，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器610上，以便于从其上读出的计算机程序根据需要被安装入存储部分608。

根据本发明的实施例，根据本发明实施例的方法流程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机程序产品，其包括承载在计算机可读存储介质上的计算机程序，该计算机程序包含用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过通信部分609从网络上被下载和安装，和/或从可拆卸介质611被安装。在该计算机程序被处理器601执行时，执行本发明实施例的系统中限定的上述功能。根据本发明的实施例，上文描述的系统、设备、装置、模块、单元等可以通过计算机程序模块来实现。

本发明还提供了一种计算机可读存储介质，该计算机可读存储介质可以是上述实施例中描述的设备/装置/系统中所包含的；也可以是单独存在，而未装配入该设备/装置/系统中。上述计算机可读存储介质承载有一个或者多个程序，当上述一个或者多个程序被执行时，实现根据本发明实施例的方法。

根据本发明的实施例，计算机可读存储介质可以是非易失性的计算机可读存储介质。例如可以包括但不限于：便携式计算机磁盘、硬盘、随机访问存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本发明中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。

例如，根据本发明的实施例，计算机可读存储介质可以包括上文描述的ROM602和/或RAM 603和/或ROM 602和RAM 603以外的一个或多个存储器。

附图中的流程图和框图，图示了按照本发明各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分，上述模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个接连地表示的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图或流程图中的每个方框、以及框图或流程图中的方框的组合，可以用执行规定的功能或操作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。本领域技术人员可以理解，本发明的各个实施例中记载的特征可以进行多种组合和/或结合，即使这样的组合或结合没有明确记载于本发明中。特别地，在不脱离本发明精神和教导的情况下，本发明的各个实施例记载的特征可以进行多种组合和/或结合。所有这些组合和/或结合均落入本发明的范围。

Claims (10)

1.一种深空探测微波载荷的干扰源确定方法，其特征在于，包括：

计算不同极化天线在同一时刻接收的微波信号之间的时域相关系数；

计算不同极化天线接收的同一频点的微波信号之间的频域相关系数；

计算同一极化天线在不同时刻接收的微波信号之间的交叉相关系数；

根据所述时域相关系数和/或所述频域相关系数和/或所述交叉相关系数确定干扰信号的来源和频率。

2.根据权利要求1所述的干扰源确定方法，其特征在于，所述计算不同极化天线在同一时刻接收的微波信号之间的时域相关系数，具体包括：

根据

计算所述时域相关系数r(X，Y)，其中，X，Y分别表示不同极化天线在同一时刻接收的微波信号，Cov(X，Y)为X，Y之间的协方差，Var[X]为X的方差，Var[Y]为Y的方差。

3.根据权利要求1所述的干扰源确定方法，其特征在于，所述计算不同极化天线接收的同一频点的微波信号之间的频域相关系数，具体包括：

计算每一微波数据的频谱；

将同一极化天线接收到的微波数据的频谱按照时间排序，得到每个频点对应的频谱数据；

将不同极化天线同一频点的频谱数据进行相关性计算，得到所述频域相关系数。

4.根据权利要求1所述的干扰源确定方法，其特征在于，所述计算同一极化天线在不同时刻接收的微波信号之间的交叉相关系数，具体包括：

从所述微波信号中确定一组具有代表性的微波信号；

对所述具有代表性的微波信号进行循环圆周移动，每移动一个时间采样点产生一个微波信号序列；

计算每个微波信号序列与相同极化天线接收到的微波信号的相关系数，得到所述交叉相关系数。

5.根据权利要求1所述的干扰源确定方法，其特征在于，所述根据所述时域相关系数和/或频域相关系数和/或所述交叉相关系数确定干扰信号的来源和频率，具体包括：

响应于所述时域相关系数小于等于第一阈值，确定得到干扰信号来源于天线内部；

响应于所述时域相关系数大于第一阈值且小于第二阈值，确定得到干扰信号来源于天线内部和天线外部；

响应于所述时域相关系数大于等于第二阈值，确定得到干扰信号来源于天线外部。

6.根据权利要求1所述的干扰源确定方法，其特征在于，所述根据所述时域相关系数和/或频域相关系数和/或所述交叉相关系数确定干扰信号的来源和频率，具体包括：

根据所述频域相关系数确定所述干扰信号的频率；或

根据所述交叉相关系数估计所述微波信号的重复周期，并根据所述重复周期估计所述干扰信号的频率。

7.根据权利要求3所述的干扰源确定方法，其特征在于，所述计算每一微波数据的频谱，具体包括：

对所述微波数据进行傅里叶变换，并计算傅里叶变换结果的幅度，得到所述微波数据的频谱；

或直接采用傅里叶变换的复数结果作为所述微波数据的频谱。

8.一种深空探测微波载荷的干扰源确定装置，其特征在于，包括：

第一计算模块，用于计算不同极化天线在同一时刻接收的微波信号之间的时域相关系数；

第二计算模块，用于计算不同极化天线接收的同一频点的微波信号之间的频域相关系数；

第三计算模块，用于计算同一极化天线在不同时刻接收的微波信号之间的交叉相关系数；

确定模块，用于根据所述时域相关系数和/或所述频域相关系数和/或所述交叉相关系数确定干扰信号的来源和频率。

9.一种电子设备，其特征在于，包括：

一个或多个处理器；

存储器，用于存储一个或多个程序，

其中，当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现权利要求1至7中任一项所述的方法。

10.一种计算机可读存储介质，其特征在于，其上存储有可执行指令，该指令被处理器执行时使处理器实现权利要求1至7中任一项所述的方法。

Images