[发明专利]一种MLT区域流星尾迹径向漂移速度的解算方法

说明书

本发明提出了一种MLT区域流星尾迹径向漂移速度的解算方法。本发明根据流星雷达所部署的所有接收机天线接收到的同一流星事件的回波复信号的复数乘积，得到平均多普勒频移复信号；根据得到的平均多普勒频移复信号通过傅里叶变换将时域信号转换为频域信号，得到该平均多普勒频移复信号的频率谱；根据平均多普勒频移复信号的频率谱，从中找到特征频率，再将该特征频率除以接收机天线数，定义为流星尾迹径向漂移速度所对应的多普勒频移；根据得到的多普勒频移和它与该流星事件对应的流星尾迹径向漂移速度之间的线性关系，解算得到该流星事件中流星尾迹的径向漂移速度。本发明较好地提高了流星尾迹径向漂移速度的计算精度，降低了计算所需的时间。

技术领域

本发明涉及空间探测技术领域，具体涉及一种MLT区域流星尾迹径向漂移速度的解算方法。

背景技术

在流星体高速飞行穿入大气并剧烈升温电离之后，流星尾迹中等离子体的双极性扩散过程成为最为重要的物理过程。伴随着双极性扩散的，是高温等离子体中的大量无规则运动的电子整体伴随着背景大气风场的漂移，这个反映了背景大气风场信息的漂移是可以被部署在地面的流星雷达所探测到的，对雷达接收的复信号解算得到的等离子体中的电子径向漂移速度和流星事件位置信息构成了基于流星雷达测风的流程中所需的基本参量。考虑到流星尾迹所在位置与雷达天线部署位置的距离往往在100km的量级，而雷达探测波和流星回波两者既近似平行又都垂直于流星的飞行方向，因此，基于高等数学上的近似，垂直于雷达观测视线上的分量对回波信号相位的影响太过微小，忽略不计，因而回波信号的相位变化集中反映了与流星尾迹径向漂移速度有关的信息，且与背景风场在垂直于径向上的分量无关。之前有大量学者的研究表明，回波信号相位的斜率与流星尾迹径向漂移速度之间存在近似的线性关系。

国际上主流的求解流星尾迹径向漂移速度的方法是通过计算不同接收机复信号之间的互相关函数的Poly-pulse pair method来求解的，这种方法有效地提高了弱信噪比条件下径向漂移速度的估计质量。然而，即便是采用这种方法，仍然存在如下问题：噪声的存在依然会给速度的求解带来较大的偶然误差；2π模糊的存在使得在对速度求解的过程中，所选取点的位置和数量必须是动态调整的，这对于需要对大量数据进行实时运算的雷达软件系统而言无疑是沉重的负担，但是简单的静态选取又往往会导致斜率的拟合出现错误；当同时存在有除流星之外的其他空间事件的回波时，也会对速度的解算形成干扰，但是因为其持续时间不长或信号能量不强，所以单纯因此而排除该流星事件又十分可惜；因为流星飞行速度的影响，相位的线性变化会伴随着振幅逐渐减小的振荡，但该振荡并不能在互相关计算中被全部抵消；大批量的流星事件数据的互相关计算消耗了大量的计算资源，因而耗时相对较长。为了更好地在测风雷达系统中求解MLT区域的大气动力学过程，需要对流星尾迹的径向漂移速度进行精度更高，耗时更短的新的解算方法。但是目前尚未有成熟的、精度更高的和耗时更短的流星尾迹径向漂移速度解算方法。

发明内容

本发明的目的是针对现有技术的技术缺陷,提供一种MLT区域流星尾迹径向漂移速度的解算方法,以对流星尾迹的径向漂移速度进行精度更高，耗时更短的解算。

为实现以上技术目的，本发明的技术方案为一种MLT区域流星尾迹径向漂移速度的解算方法，具体步骤如下：

步骤1，根据流星雷达所部署的所有接收机天线接收到的同一流星事件的回波复信号的复数乘积得到平均多普勒频移复信号；

步骤2，对平均多普勒频移复信号通过傅里叶变换将时域信号转换为频域信号得到平均多普勒频移复信号的频率谱；

步骤3，从平均多普勒频移复信号的频率谱中搜索特征频率，将特征频率除以接收机天线数，定义为流星尾迹径向漂移速度所对应的多普勒频移；

步骤4，根据流星尾迹径向漂移速度所对应的多普勒频移以及流星事件对应的流星尾迹径向漂移速度之间的线性关系，解算得到流星事件中流星尾迹的径向漂移速度。