2024深圳杯东三省A题完稿 多个火箭残骸的准确定位

摘 要

对于问题一：本文着眼于精确确定火箭残骸发生音爆时的位置和时间，采用了多边定位法建立定位模型。该模型利用了多个监测设备的数据，并通过优化算法求解，得到了引爆位置的精确坐标。这种方法能够高效地解决单个残骸定位问题，且模型的准确度较高。

对于问题二与问题三：考虑到可能出现多个残骸同时发生音爆的情况，使用了差分进化算法建立了复杂的多残骸定位模型。在这个模型中，通过对各残骸到各监测设备的预期到达时间与实际接收到的震动波到达时间进行比较，确定了各残骸的位置和时间。这种方法有效应对了多残骸定位问题的复杂性，提高了定位的准确性。

对于问题四：考虑了监测设备记录时间存在的随机误差，引入了粒子群优化算法来进一步提高定位精度。该算法能够有效处理随机误差带来的影响，使得定位结果更加可靠。

最后，本文对使用的模型进行了优缺点分析：多边定位法在精度和适用范围方面表现出色，能够通过多个监测设备的数据计算目标物体的位置，具有更高的定位精度，并且适用于不同环境和条件下。与此同时，它不要求监测设备的时间同步，从而降低了系统的复杂性。此外，差分进化算法在全局搜索能力方面表现突出，能够在解空间中搜索到较好的解，适用于复杂的多峰问题和高维问题。

关键词：粒子群算法；差分进化算法；多边定位法；优化模型

一、问题重述

1.1 问题背景

在现代航天领域，火箭扮演着实现空间探索的关键角色。然而，由于火箭发射的高成本和复杂性，对火箭残骸的回收和重复利用变得至关重要。多级火箭发射是常见的方法之一，它涉及到在任务完成后分离并使各个阶段的助推器返回地球。这些坠落的部件不仅具有潜在的再利用价值，而且可能对地面设施和安全造成威胁。因此，对这些火箭残骸的准确定位和迅速回收至关重要。通常情况下，多级火箭的每个阶段在完成其推进任务后会自动分离，通过预设的分离装置进行。这些分离的组件包括助推器和发动机等，在预定的轨道上返回地球。在返回过程中，由于高速穿越大气层，残骸可能会达到超音速，导致产生跨音速音爆。这种音爆是由于残骸的运动速度超过声音传播速度而产生的强烈震动波。

1.2 问题要求

基于上述背景以及所获取的数据需要建立数学模型解决以下问题：

（1）建立数学模型以研究如何精确确定空中单个火箭残骸发生音爆时的位置坐标（包括经度、纬度和高程）及时间是至关重要的。假设在预计的落点附近布置了7台监测设备，每台设备都记录了音爆抵达的时间。现在的问题是如何利用这些数据来选取适当的信息，以计算音爆发生的具体位置和时间。

（2）当多架火箭（例如一级残骸和多个助推器）同时发生音爆时，监测设备可能会收到多组不同的震动波数据。为了确定每组数据对应的残骸，需要建立数学模型。此外，需要讨论至少需要多少台监测设备才能准确确定所有残骸。

（3）建立数学模型，对附件三所给出的非连续时间段的电磁辐射与声发射信号数据进行分析，对每个时间段最后出现前兆特征的数据进行预测。

二、问题分析

针对问题一：分析在精准确定空中单个残骸发生音爆时的位置坐标和时间时，至少需要布置几台监测设备。考虑到需要解出三个未知数（经度、纬度、高程），因此至少需要三个方程，即至少需要三个监测设备。在提供的情景中，火箭一级残骸分离后，在落点附近布置了7台监测设备，因此满足了至少三个监测设备的条件。因此，利用这些监测设备的数据可以精准确定空中单个残骸发生音爆时的位置坐标和时间。

针对问题二：分析如何确定监测设备接收到的震动波是来自哪一个残骸，以及在空中多个残骸发生音爆时至少需要布置多少台监测设备。可以通过计算每个残骸到每个监测设备的预期到达时间，并与实际接收到的震动波到达时间进行比较，从而确定来自哪个残骸的震动波。在提供的情景中，若空中有四个残骸产生音爆，每个设备按时间先后顺序收到四组震动波，则至少需要布置四个监测设备，以确保有足够的方程来解出这些未知数。

针对问题三：通过问题二的数学模型，可以利用这些数据确定四个残骸在空中发生音爆时的位置和时间。算法将比对实际接收到的震动波到达时间与预期到达时间，以确定每个残骸的位置和时间。

针对问题四：要求修正模型以较精确地确定残骸在空中发生音爆时的位置和时间，考虑设备记录时间存在0.5秒的随机误差。可以通过引入时间误差修正来提高精度，一种方法是使用最小二乘法拟合监测设备接收到的震动波到达时间和残骸预期到达时间之间的关系，从而减小误差。另外，可以使用更准确的定位设备和更精细的时间同步系统来提高定位精度。在模拟时，可以考虑添加随机误差，然后利用修正后的模型来验证定位结果的精确度。

四、模型的建立与求解

5.1 数据处理

5.1.1 数据可视化

为了刻画附件中的地理位置信息，本文首先对附件进行可视化，结合经度、纬度与高度进行可视化绘图，如图1、图2。结果显示，各个设备处于不同的高度上，并不规律，同时对各设备所覆盖度区域进行对应标记，为后续的模型建立做准备。

图1-图2 数据可视化

5.2问题一模型的建立与求解

5.2.1模型的准备

为使得后续的计算更高效有序，在此建立笛卡尔坐标系，使用近似法将经纬度映射到此坐标系之中，具体的转换公式如式（1）-（3）。

为该点的经度，为了更方便计算，减小问题的计算难度，在此进一步进行简化计算，具体转化公式如式（4）-（6）。

由于篇幅限制，如需获取成品论文，联系作者，具体成品全览如下图