导弹与航天运载技术（中英文）

设备	经度/（°）	纬度/（°）	高程/m	音爆抵达时间/s
A	110.241	27.204	824	100.767
B	110.780	27.456	727	112.220
C	110.712	27.785	742	188.020
D	110.251	27.825	850	258.985
E	110.524	27.617	786	118.443
F	110.467	27.921	678	266.871
G	110.047	27.121	575	163.024

设备	经度/（°）	纬度/（°）	高程/m	音爆抵达时间/s
A	110.241	27.204	824	100.767
B	110.780	27.456	727	112.220
C	110.712	27.785	742	188.020
D	110.251	27.825	850	258.985
E	110.524	27.617	786	118.443
F	110.467	27.921	678	266.871
G	110.047	27.121	575	163.024

设备	经度 /（°）	纬度 /（°）	高程/m	音爆抵达时间/s
A	110.241	27.204	824	100.767	164.229	214.850	270.065
B	110.783	27.456	727	92.453	112.220	169.36	196.583
C	110.762	27.785	742	75.560	110.696	156.936	188.020
D	110.251	28.025	850	94.653	141.409	196.517	258.985
E	110.524	27.617	786	78.600	86.216	118.443	126.669
F	110.467	28.081	678	67.274	166.270	175.482	266.871
G	110.047	27.521	575	103.738	163.024	206.789	210.306

设备	经度 /（°）	纬度 /（°）	高程/m	音爆抵达时间/s
A	110.241	27.204	824	100.767	164.229	214.850	270.065
B	110.783	27.456	727	92.453	112.220	169.36	196.583
C	110.762	27.785	742	75.560	110.696	156.936	188.020
D	110.251	28.025	850	94.653	141.409	196.517	258.985
E	110.524	27.617	786	78.600	86.216	118.443	126.669
F	110.467	28.081	678	67.274	166.270	175.482	266.871
G	110.047	27.521	575	103.738	163.024	206.789	210.306

设备	经度/（°）	纬度/（°）	高程/m	音爆抵达时间/s
H	109.850	19.620	102	85.323
I	110.200	19.910	88	92.152
J	110.550	20.150	75	105.671
K	110.002	20.430	120	116.841

设备	经度/（°）	纬度/（°）	高程/m	音爆抵达时间/s
H	109.850	19.620	102	85.323
I	110.200	19.910	88	92.152
J	110.550	20.150	75	105.671
K	110.002	20.430	120	116.841

设备	经度/（°）	纬度/（°）	高程/m	残骸1音爆抵达时间/s	残骸2音爆抵达时间/s
L	112.400	25.650	95	76.214	82.451
M	112.732	25.900	85	81.341	88.920
N	113.000	26.150	105	85.672	92.112
O	112.200	25.450	110	90.120	96.342
P	112.500	25.750	90	94.563	98.781

设备	经度/（°）	纬度/（°）	高程/m	残骸1音爆抵达时间/s	残骸2音爆抵达时间/s
L	112.400	25.650	95	76.214	82.451
M	112.732	25.900	85	81.341	88.920
N	113.000	26.150	105	85.672	92.112
O	112.200	25.450	110	90.120	96.342
P	112.500	25.750	90	94.563	98.781

残骸	x/km	y/km	h/km
长征三号	17.358	81.265	0.910
长征五号B遥三	44.290	55.743	0.105
长征七号甲遥五残骸1	25.016	25.924	0.090
长征七号甲遥五残骸2	46.205	9.346	0.095

残骸	x/km	y/km	h/km
长征三号	17.358	81.265	0.910
长征五号B遥三	44.290	55.743	0.105
长征七号甲遥五残骸1	25.016	25.924	0.090
长征七号甲遥五残骸2	46.205	9.346	0.095

设备	x/km	y/km	h/km	音爆抵达时间/s
A	18.884	9.235	0.824	100.767
B	71.350	37.273	0.727	112.220
C	64.731	73.879	0.742	188.020
D	19.857	78.329	0.850	258.985
E	46.431	55.186	0.786	118.443
F	40.883	89.010	0.678	266.871
G	0	0	0.575	163.024
H	0	0	0.102	85.323
I	34.069	32.266	0.088	92.152
J	68.138	58.969	0.075	105.671
K	14.796	90.123	0.120	116.841
设备	x/km	y/km	h/km	残骸1音爆抵达时间/s	残骸2音爆抵达时间/s
L	19.468	22.253	0.095	76.214	82.451
M	51.785	50.068	0.085	81.341	88.920
N	77.872	77.884	0.105	85.672	92.112
O	0	0	0.110	90.120	96.342
P	29.202	33.379	0.090	94.563	98.781

设备	x/km	y/km	h/km	音爆抵达时间/s
A	18.884	9.235	0.824	100.767
B	71.350	37.273	0.727	112.220
C	64.731	73.879	0.742	188.020
D	19.857	78.329	0.850	258.985
E	46.431	55.186	0.786	118.443
F	40.883	89.010	0.678	266.871
G	0	0	0.575	163.024
H	0	0	0.102	85.323
I	34.069	32.266	0.088	92.152
J	68.138	58.969	0.075	105.671
K	14.796	90.123	0.120	116.841
设备	x/km	y/km	h/km	残骸1音爆抵达时间/s	残骸2音爆抵达时间/s
L	19.468	22.253	0.095	76.214	82.451
M	51.785	50.068	0.085	81.341	88.920
N	77.872	77.884	0.105	85.672	92.112
O	0	0	0.110	90.120	96.342
P	29.202	33.379	0.090	94.563	98.781

残骸类别	修正前	修正后	距离误差/km
1	（37.121，60.844，0.778 3）	81.041	（37.473，61.752，0.779 6）	80.854	0.34
2	（29.445，50.015，0.757 1）	75.274	（29.521，49.314，0.765 9）	77.243	0.41
3	（17.856，43.978，0.701 1）	75.703	（17.735，44.833，0.687 4）	81.572	0.27
4	（27.337，49.626，0.353 1）	77.016	（27.128，48.992，0.385 7）	82.052	0.65

残骸类别	修正前	修正后	距离误差/km
1	（37.121，60.844，0.778 3）	81.041	（37.473，61.752，0.779 6）	80.854	0.34
2	（29.445，50.015，0.757 1）	75.274	（29.521，49.314，0.765 9）	77.243	0.41
3	（17.856，43.978，0.701 1）	75.703	（17.735，44.833，0.687 4）	81.572	0.27
4	（27.337，49.626，0.353 1）	77.016	（27.128，48.992，0.385 7）	82.052	0.65

设备	经度/（°）	纬度/（°）	高程/km	音爆抵达时间/s
Q	99.450	27.890	1.500	145.323
R	99.780	28.120	0.980	152.670
S	99.692	28.350	2.200	138.451
T	99.600	28.250	1.800	142.890
实际残骸位置	99.875	28.212	1.200	—

设备	经度/（°）	纬度/（°）	高程/km	音爆抵达时间/s
Q	99.450	27.890	1.500	145.323
R	99.780	28.120	0.980	152.670
S	99.692	28.350	2.200	138.451
T	99.600	28.250	1.800	142.890
实际残骸位置	99.875	28.212	1.200	—

设备	x/km	y/km	h/km	音爆抵达时间/s
Q	0	0	1.500	145.323
R	32.122	25.590	0.980	152.670
S	23.56	51.181	2.200	138.451
T	14.601	40.055	1.800	142.890
实际残骸位置	41.369	35.827	1.200	—

设备	x/km	y/km	h/km	音爆抵达时间/s
Q	0	0	1.500	145.323
R	32.122	25.590	0.980	152.670
S	23.56	51.181	2.200	138.451
T	14.601	40.055	1.800	142.890
实际残骸位置	41.369	35.827	1.200	—

设备	经度/（°）	纬度/（°）	高程/ km	音爆抵达时间/s	修正后时间/s
U	115.230	30.550	0.080	88.760	90.212
V	115.622	31.850	0.095	92.342	93.790
W	116.000	31.105	0.105	96.451	97.903
X	115.400	30.700	0.090	89.120	90.571
实际残骸位置	115.805	30.938	0.093	—	—

设备	经度/（°）	纬度/（°）	高程/ km	音爆抵达时间/s	修正后时间/s
U	115.230	30.550	0.080	88.760	90.212
V	115.622	31.850	0.095	92.342	93.790
W	116.000	31.105	0.105	96.451	97.903
X	115.400	30.700	0.090	89.120	90.571
实际残骸位置	115.805	30.938	0.093	—	—

设备	x/km	y/km	h/km	音爆抵达时间/s	修正后时间/s
U	0	0	0.080	88.760	90.212
V	38.157	144.642	0.095	92.342	93.790
W	74.952	61.751	0.105	96.451	97.903
X	16.548	16.689	0.090	89.120	90.571
实际残骸位置	55.970	43.170	0.090	—	—

设备	x/km	y/km	h/km	音爆抵达时间/s	修正后时间/s
U	0	0	0.080	88.760	90.212
V	38.157	144.642	0.095	92.342	93.790
W	74.952	61.751	0.105	96.451	97.903
X	16.548	16.689	0.090	89.120	90.571
实际残骸位置	55.970	43.170	0.090	—	—

火箭残骸精准定位算法优化

PDF下载

赵铭 , 腾云

导弹与航天运载技术（中英文） | 运载器及导弹总体技术 2025,48(5): 13-23

收起

导弹与航天运载技术（中英文） | 运载器及导弹总体技术 2025, 48(5): 13-23

火箭残骸精准定位算法优化

全屏

赵铭, 腾云

作者信息

东北农业大学工程学院，哈尔滨，150030

赵铭（2004—），男，本科，主要研究方向为系统工程。

腾云（1980—），女，博士研究生，教授，主要研究方向为系统工程。

Optimization of Precise Location Algorithm for Rocket Debris

Ming ZHAO, Yun TENG

Affiliations

College of Engineering, Northeast Agricultural University, Harbin, 150030

出版时间: 2025-10-25 doi: 10.7654/j.issn.2097-1974.20250508

文章导航

摘要

收起

针对多级火箭残骸定位精度不足的问题，提出一种融合粒子群算法与牛顿迭代-梯度下降法的时间差定位方法（Time Difference of Arrival，TDOA）联合优化模型，通过全局搜索与局部优化的协同机制提升定位性能，并构建多残骸信号分离约束模型与环境干扰补偿模型。试验表明，该模型在火箭残骸回收任务中，定位误差由传统单级优化算法的1~10 km降低至0.5 km以内，多残骸信号分离率达96.2%，山地及强风干扰下仍保持亚千米级精度。结合Chan-Taylor算法与最小二乘法的验证表明，其抗干扰性与定位可靠性显著优于现有方法。本算法可拓展至移动通信、无人驾驶等领域，兼具理论创新与工程应用价值。

关键词

火箭残骸 / 精准定位 / TDOA技术 / 联合算法优化 / Chan-Taylor算法

Abstract

收起

A TDOA joint optimization model (TDOA-PSO-NI-GD) that integrates Particle Swarm Optimization (PSO) and Newton's Iteration-Gradient Descent (NI-GD) method is proposed in response to the problem of insufficient positioning accuracy of multi-stage rocket debris, which improves the localization performance through the synergistic mechanism of global search and local optimization and constructs the multi-debris signal separation constraint model and the environmental interference compensation model. Experiments show that the model reduces the positioning error from 1~10 km to less than 0.5 km in the traditional single-stage optimization algorithm in the rocket debris recovery mission, and the multi-debris signal separation rate reaches 96.2%, and maintains the sub-kilometer accuracy under the mountainous terrain and strong wind disturbances; the validation of the Chan-Taylor algorithm combined with the least-squares method shows that its anti-jamming and positioning reliability are significantly better than that of the existing methods. The algorithm can be extended to mobile communications, unmanned vehicles and other fields, with both theoretical innovation and engineering application value.

Key words

rocket debris / precise localization / TDOA technology / joint algorithm optimization / Chan-Taylor algorithm

引用本文

赵铭, 腾云. 火箭残骸精准定位算法优化. 导弹与航天运载技术（中英文）, 2025 , 48 (5) : 13 -23 . DOI: 10.7654/j.issn.2097-1974.20250508

Ming ZHAO, Yun TENG. Optimization of Precise Location Algorithm for Rocket Debris[J]. Missiles and Space Vehicles, 2025 , 48 (5) : 13 -23 . DOI: 10.7654/j.issn.2097-1974.20250508

正文

收起

0 引言

收起

在国家战略需求和科技进步的双重驱动下，中国航天事业正迎来从“航天大国”迈向“航天强国”的关键转型期。随着“一带一路”倡议的深入推进以及中国空间站的持续建设，中国航天任务的规模不断扩大，频率不断提升，涵盖了通信卫星、导航卫星、气象卫星、探月工程及深空探测等多类型发射任务。航天火箭发射次数的增加带来了更多残骸，不仅造成环境污染，而且可能会对地面设施和人员安全构成潜在威胁^［1-2］。在多级发射任务中，火箭及助推器各级结构按照设计要求，在完成各阶段的既定任务后会与主火箭分离，并依据发射计划落入预定区域。这些残骸虽然在发射前通过落区划分得到了初步规划，但在高空高速下落的过程中，会受到复杂环境因素的干扰，偏离预设区域。因此，火箭残骸精准定位的研究日益紧迫。

目前，火箭残骸定位回收主要依赖于视觉、雷达等传统监测手段，但这些方法在实际操作中存在诸多局限性。由于火箭残骸通常坠落在地形复杂、偏远的山区或海域，视觉和雷达监测易受到地形遮挡、天气干扰等外部因素的影响，难以实现精确定位。此外，火箭残骸坠落速度极快，音爆信号传播时间极短，传统监测技术在短时精确捕捉、处理信号等技术性问题上遭遇瓶颈。为克服这些难点，基于音爆震动波到达时间差（Time Difference of Arrival，TDOA）的定位技术逐渐成为研究热点^［3-6］。TDOA技术通过布设多台监测设备记录音爆信号到达时间差，利用信号传播速度推算音爆源位置，其核心优势在于不需要依赖光学或雷达设备，可有效适应复杂地形和恶劣天气。该技术已在无线传感器网络定位^［7］、灾害源追踪^［8］和移动通信基站优化^［9］等领域取得显著成果。例如，Chen等^［10］通过联合TDOA与频差技术实现了对移动目标厘米级定位，Qu等^［11］提出改进的加权最小二乘法，显著降低了TDOA定位误差，公式化问题全局最优解达到CRLB精度，20次迭代内收敛百分比超过96%。然而，现有研究多聚焦于单目标场景，在火箭残骸多源音爆交织、信号非线性叠加等复杂情境下，传统TDOA算法易因雅可比矩阵奇异或初值敏感导致定位发散^［12-13］。

近年来，学者们尝试将智能优化算法与传统定位模型结合以突破上述限制。粒子群算法（Particle Swarm Optimization，PSO）因其全局搜索能力强、参数设置简单等特性，在危险气体泄漏源定位^［14］和无人机集群路径规划^［15］中表现出色；牛顿迭代法（Newton Iteration Method，NIM）与梯度下降法（Gradient Descent，GD）则凭借局部收敛速度快的特点，被应用于光伏功率预测^［16］和非线性方程组求解^［17］。然而，单一算法在火箭残骸定位场景中仍存在明显短板，PSO局部最优但收敛精度不足^［18］，NI-GD对初始值依赖性强，在噪声干扰下稳定性较差^［19］。现有研究尚未系统探讨多算法协同优化在TDOA定位中的可行性，特别是针对火箭多残骸场景的信号分离与抗干扰机制仍缺乏理论支撑。

综上，本文基于TDOA定位技术提出一种融合PSO全局搜索与NI-GD局部优化的混合定位模型TDOA-PSO-NI-GD，以长征系列火箭残骸回收为研究对象，重点解决3个关键问题：a）设计两阶段优化框架，通过PSO快速锁定潜在解域，再经NI-GD精细化修正，突破单一算法精度瓶颈；b）构建多残骸信号分离约束模型，引入聚类分析与时间序列匹配机制，实现交织音爆信号的精准识别；c）建立环境干扰补偿模型，通过引入高程权重因子与风速漂移修正项，提升复杂地形下算法鲁棒性。本研究不仅为火箭残骸高效回收提供了新方法，其优化框架还可拓展至地质灾害监测、海上搜救等多元定位场景，具有一定的工程应用价值。

1 材料与方法

收起

1.1 研究区概况及残骸位点总览

1.1.1 西昌卫星发射中心

研究区1为西昌卫星发射中心，位于四川省凉山彝族自治州，是中国重要的航天发射场之一，海拔约为1 500 m。该中心地处大凉山腹地，其区域地形以山地为主。

西昌卫星发射中心的周边区域自然条件复杂，包括山地、河流、森林和部分农业区。虽然山区的地形使得发射活动产生的影响相对较小，但对回收火箭残骸可能带来一定的挑战。火箭一级和二级残骸可能会落在周边的四川、云南、贵州等地。因此，可以建立西昌周围的山脉模型，并分离出西昌、云南和贵州等区域，标注出残骸大致区域。运用Python编程软件生成西昌卫星发射中心布局图及其火箭残骸分布区域，如图1所示。

1.1.2 文昌航天发射中心

研究区2为文昌航天发射中心，位于海南省文昌市龙楼镇，坐落于海南岛东北部沿海地区。

相比其他内陆发射场，文昌航天发射中心海拔高度较低，约10 m。其周边多为平原和台地，便于大型设备的运输和安装，为发射场的建设及后续航天活动的开展提供了良好的基础条件。同时，靠近海洋的地理位置使其在火箭发射后，残骸可坠落在海洋区域，大大降低了残骸对地面人员和设施的安全风险。

近些年来，文昌航天中心火箭发射任务激增，分别于2022年7月24日14时22分和2022年9月13日21时18分发射了长征五号B遥三运载火箭以及长征七号甲遥五运载火箭，应用Python中Folium库将其火箭残骸的坠落路线和残骸落点的大致区域可视化，如图2和图3所示。

1.2 数据采集

a）机器设备采集。

为保证采集数据的广泛性和准确性，监测设备需要安装在预定区域内的多个地理位置，每个设备点选址均基于火箭残骸可能坠落区域的地形、飞行轨迹预测，以及音爆在空气中的传播模型。应在既定区域内选取多个监测点，分别记录每个监测点的经度、纬度和高程，形成网格化分布，确保涵盖不同地形和高度的区域。每个监测点安装声学传感器、GPS定位设备、气象监测仪等。声学传感器能够以毫秒级分辨率检测火箭残骸音爆到达的时间，GPS能精确记录监测设备的地理坐标，气象监测仪能记录温度、气压及风速等数据。由于音爆传播会受到地形和天气条件的影响，不同设备记录的抵达时间可能有所不同，但均在误差范围之内。根据Badawi等^［20］的方法，模拟机械设备监测数据如表1—4所示。

b）现场勘探。

现场勘探通常在火箭发射后，通过实地勘查残骸坠落区域，结合地理和物理测量技术，获取火箭残骸的位置信息以及相关环境数据。勘探区域主要集中在音爆抵达时间较短、地形适中且人员易抵达的地点。现场勘探残骸实际位点数据见表5^［21］。

1.3 数据处理

在建模前需要对数据进行转化，将每台设备的经度、纬度、高程转化到基准设备的距离坐标。假设A的经纬度为

(x A, y A)

，其余设备经纬度为

(x i, y i) (i = B, C, ⋯, G)

。根据纬度间每度距离值近似为111.263 km，经度间每度距离值近似为97.304 km，转化公式为

97.34 × x i - x A 111.263 × y i - y A

（1）

转化后的设备坐标值如表6所示。

1.4 研究方法

1.4.1 时差定位

三维空间内的TDOA定位技术主要基于双曲线定位原理和信号传播的时间与距离关系，从而实现复杂条件下既定目标的准确定位^［22］。

三维空间内定位需要4个以上的基站。设基站坐标为A

x 1, y 1, z 1

、B

x 2, y 2, z 2

、C

x 3, y 3, z 3

、D

x 4, y 4, z 4

等，移动台坐标为M

x, y, z

。

移动台M到基站A、B的距离分别为

d 1 = x - x 1 2 + y - y 1 2 + z - z 1 2

（2）

d 2 = x - x 2 2 + y - y 2 2 + z - z 2 2

（3）

信号到达基站A、B的时间差为

Δ t A B

，则有：

d 2 - d 1 = v ⋅ Δ t A B

（4）

同理，对于其他基站对也可以得到类似的结果，通过联立多个上述原理的方程，即可求解出三维空间中移动台的坐标

x, y, z

。

1.4.2 粒子群算法

PSO是一种基于群体智能的优化算法，灵感来自鸟群或鱼群等社会性生物的群体行为，它通过模拟粒子在解空间中的移动和交流来搜索最优解^［23-24］。

基本PSO算法：

a）D维空间中，有m个粒子，则：

粒子

i

位置：

x i = x i 1, x i 2, …, x i D

；

粒子

i

速度：

v i = v i 1, v i 2, …, v i D ， 1 ≤ i ≤ m, 1 ≤ d ≤ D

；

粒子

i

经历过的历史最好位置：

p i = p i 1, p i 2, …, p i D

；

群体内（或领域内）所有粒子所经历过的最好位置：

p g = p g 1, p g 2, …, p g D

。

一般来说，粒子的位置和速度都是在连续的实数内进行取值。

b）基本PSO公式^［18］：

v i D k + 1 = v i D k + c 1 r 1 p i D k - x i D k + c 2 r 2 p g D k - x i D k

（5）

x i D k + 1 = x i D k + v i D k + 1

（6）

式中

c 1

，

c 2

为学习因子或加速系数，一般为正常数，通常等于2；

r 1

，

r 2

的取值范围为

0 ， 1

，该区间内是均匀分布的伪随机数。

c）基本PSO流程图。

粒子

i

第

d

维速度更新公式^［25］：

v i d k + 1 = ω * v i d k + c 1 r 1 p i d k - x i d k + c 2 r 2 p g d k - x i d k

（7）

粒子

i

第

d

维位置更新公式^［25］：

x i d k + 1 = x i d k + v i d k + 1

（8）

式中

v i d k

为第

k

次迭代粒子

i

飞行速度矢量的第

d

维分量；

x i d k

为第

k

次迭代粒子

i

位置矢量的第

d

维分量；

c 1

，

c 2

为学习因子；

r 1

，

r 2

为两个随机函数，取值范围是

0,1

，以增加搜索随机性；

ω

为惯性权重，调节对解空间的搜索能力。PSO具体流程如图4所示。

1.4.3 梯度下降法

梯度下降法是一种用于最小化函数的数值方法，通过迭代方式，沿着函数梯度的反方向逐步调整参数，直至找到函数的局部最小值或全局最小值^［26］。

其基本原理及相关公式如下：

a）初始化：随机或给定一个初始参数向量 θ。

b）计算梯度：计算目标函数关于参数向量 θ 的梯度，即目标函数的偏导数，表示为

∇ J θ

，其中

J θ

是目标函数。

c）更新参数：根据梯度的反方向调整参数向量 θ，以减小目标函数的值，其定义形式如式（9）所示^［26］：

θ = θ - α ⋅ ∇ J θ

（9）

式中 α为学习率（也称为步长），控制参数更新步伐大小。

d）重复迭代：重复执行步骤b和步骤c，直至满足停止条件，例如达到最大迭代次数或梯度小于某个阈值。

1.4.4 牛顿迭代法

NIM是一种用于求解方程的迭代数值方法，尤其适用于求解非线性方程或优化问题^［27］。其基本原理及相关公式如下：

a）初始化：给定初始近似值

x 0

。

b）计算下一个近似值：基于泰勒级数展开，用函数在当前近似值处的一阶导数（梯度）和二阶导数（Hessian矩阵）来逼近函数。下一个近似值通过式（10）计算：

x 1 = x 0 - f' x 0 / f'' x 0

（10）

式中

f' x 0

为函数f在点x₀处的一阶导数（梯度）；

f ″ x 0

为函数f在点x₀处的二阶导数（Hessian矩阵）。

c）迭代更新：重复步骤b，使用新计算的近似值x₁作为下一个迭代的初始值，直至满足停止条件（如达到指定的精度或迭代次数）。

综上所述，传统TDOA技术在定位场景下面临两大挑战：信号来源难以区分以及非线性方程求解效率低。传统TDOA技术、PSO与NI-GD结合的核心是协同发挥两种算法的互补优势，即同时具备PSO的全局搜索能力与NI-GD的局部收敛效率。因此，本研究流程优化为：

a）第一阶段，借助PSO群体智能特性，快速搜索解空间，初步确定残骸位置的潜在区域。通过设置较大搜索范围与较高粒子多样性，防止陷入局部最优。

b）第二阶段，将PSO输出最优解作为NI-GD初始值，利用牛顿法二阶导数信息和梯度下降法的稳定性，进行高精度局部优化，修正残骸坐标的细微偏差。

2 结果与分析

收起

2.1 单残骸定位研究

若有n个监测设备，第i个监测设备的位置坐标为

A i x i, y i, h i

，其中x_i 为第i个监测设备经度在水平面的投影横坐标，y_i 为第i个监测设备纬度在水平面的投影纵坐标，h_i 为第i个监测设备的高程。B为音爆位置

x ， y ， h

，由TDOA定位技术推导如下：

d i = x i - x 2 + y i - y 2 + h i - h 2, i = 1, 2, …, n

（11）

音爆时各站点接收信号的时间t与距离存在关系为

d i = v t i

（12）

式中 v为声速；t_i 为从音爆开始至各站点接收信号的时间。

根据式（12）以及TDOA定位原理可得：

x 1 - x 2 + y 1 - y 2 + h 1 - h 2 = v 2 t 12 x 2 - x 2 + y 2 - y 2 + h 2 - h 2 = v 2 t 22 ⋮ x n - x 2 + y n - y 2 + h n - h 2 = v 2 t n 2

（13）

d i - d i' = v t i, (i ∈ 1,2, ⋯, n; i' ∈ 1,2, …, n)

（14）

分别采用上述粒子群算法及NI-GD对式（12）求解。由于三维空间中定位残骸至少需要4台监测设备，因此为寻找音爆位置，将前4、5、6、7个监测设备的数据带入式（13），结果如图5所示。

图5显示，在4种不同情形下，TDOA-PSO定位误差大致为0.92 km、13.45 km、0.782 5 km和2.872 km；TDOA-NI-GD的优化误差约为1.25 km、10.85 km、1.85 km和0.982 km。由此可见，不同情况下，TDOA分别与PSO和NI-GD结合的定位结果偏差较大，稳定性欠佳。为提升定位准确性，将TDOA定位技术与PSO、NI-GD优化技术相结合。

结合TDOA定位技术的两阶段优化模型：

阶段1：应用PSO全局搜索最小化残差平方和：

f x, y, h, T = ∑ i = 1 n x i - x 2 + y i - y 2 + h i - h 2 - v t i - T 2

（15）

其更新公式如下：

v i d k + 1 = w ⋅ v i d k + c 1 r 1 P i d - x i d k + c 2 r 2 g d - x i d k

（16）

x i d k + 1 = x i d k + v i d k + 1

（17）

阶段2：以PSO最优解

x p s o, y p s o, h p s o, T p s o

为初始值进行牛顿迭代和梯度下降：

x k + 1 = x k - H - 1 x k ∇ f x k; x k + 1 = x k - α ∇ f x k

（18）

式中 H 为Hessian矩阵。

以长征五号B遥三运载火箭一级残骸为研究对象，应用TDOA-PSO、TDOA-NI-GD、TDOA-PSO-NI-GD三种不同的优化算法结果可视化呈现于图6。

图6表明，在单残骸定位中，TDOA-PSO与TDOA-NI-GD定位结果仍不稳定。相比之下，TDOA定位技术结合粒子群与NI-GD综合优化后，定位误差降至0.428 5 km，远低于平均值。

2.2 多残骸定位研究

除了一级残骸，火箭残骸还有2个或4个助推器^［28］。在多个火箭残骸发生音爆时，监测设备在监测范围内可能会采集到几组音爆数据。

2.2.1 确定最小监测设备数

若设备监测到4组信号，则需要对各监测设备的数据信号进行分类，且各数据只能选择一次。即假设设备接收到残骸音爆抵达时间为 T，其中：

T = t 11 t 12 t 13 t 14 t 21 t 22 t 23 t 24 ⋮ ⋮ ⋮ ⋮ t n 1 t n 2 t n 3 t n 4

（19）

首先，定义变量a_ij 是0~1变量，即：

a i j = 0 ， 第 i 个音 爆抵 达时 间不 来自 第 j 个残 骸 1 ， 第 i 个音 爆抵 达时 间来 自第 j 个残 骸

（20）

由发生爆炸时间+音爆传播时间=设备接收到的时间，可建立如下模型：

m i n ∑ i = 1 4 ∑ j = 1 4 ∑ k = 1 7 t i j k - t k - d j k / 340 2 ⋅ a i j

（21）

∑ j = 1 4 x i j = 1 ∑ i = 1 4 x i j = 1

（22）

式中

∑ j = 1 4 x i j = 1

表示每个音爆只能来自一个残骸；

∑ i = 1 4 x i j = 1

表示每个残骸只能产生一个音爆。假设M点为音爆位置，N与P两点为设备位置（M，N，P位置在残骸可能坠落区域内随机选取），可建立三角形模型。

由三角形法则可得：

M N - M P < M N < M P + N P

（23）

即：

v (t p - t M) - d N P < v (t N - t M) < v (t p - t M) + d N P

（24）

整理可得：

t k ″ + t j - t k' k ″ / v < t k' - t j < t k ″ - t j + d k' k ″ / v

（25）

对于每台设备k，假设已知音爆的抵达时间顺序为p₁，p₂，p₃，p₄（其中p_i 为第i个音爆编号）。需要确保音爆p_i 对应的残骸编号在音爆p_j 对应的残骸编号之前或相同（但由于残骸编号唯一，所以实际上不会相同），音爆抵达时间顺序约束的不等式为

∑ l = 1 j - 1 x p i l ≥ ∑ l = 1 j x p i l, ∀ j = 2, 3, 4; ∀ k

（26）

综上，此模型表述为

m i n ∑ i = 1 4 ∑ j = 1 4 ∑ k = 1 7 t i j k - t k - d j k / 340 2 ⋅ a i j ∑ j = 1 4 x i j = 1, ∑ i = 1 4 x i j = 1 t k ″ + t j - t k' k ″ / v < t k' - t j < t k ″ - t j + d k' k ″ / v ∑ l = 1 j - 1 x p i l ≥ ∑ l = 1 j x p i l, ∀ j = 2, 3, 4; ∀ k

（27）

计算最小监测设备数：根据上述模型可知，未知量共有4n个，等式约束共有2n个，不等式约束共有

C n + 1 2

个，建立方程如下：

q = 4 n - 2 n - C n 2 - 1

（28）

根据线性方程组的求解原理可得，q>0时，整个模型可能存在多个解，q<0时，整个模型可能存在唯一解。当n的值在4到5之间接近于5时，q=0，唯一解不存在，当n大于这个数时，q<0，可能存在唯一解，由于n取整数，因此最小设备数为5。

2.2.2 TDOA-PSO-NI-GD的多残骸定位优化

在多残骸定位场景下，不同残骸发射的信号可能相互交织、干扰，若不能有效分离，会导致定位信息的混淆与错误，进而严重影响定位的准确性与可靠性。因此，为实现精准的多残骸定位，信号分离是关键的考量因素。

在原有模型的基础上，需要加入信号分离约束：

a）输入矩阵

T = t i j N × M

，其中t_ij 表示第i个设备接收到第j个残骸音爆的时间。

b）将 T 聚类为M组信号，每组对应一个残骸。

每个 音爆 仅属 于一 个残 骸： ∑ k = 1 M a j k = 1 ∀ j 每个 残骸 至少 分配 一个 信号 ： ∑ j = 1 M a j k ≥ 1 ∀ k

（29）

即最终多残骸定位的TDOA-PSO-NI-GD的综合优化方程为

阶段 1 P S O ： m i n {a j k} ∑ k = 1 M ∑ i = 1 N ∑ j = 1 M a j k t i j - t k ¯ 2 + λ ∑ l ≠ k t k ¯ - t l ¯ - 1 阶段 2 N I - G D ： F o r e a c h k ∈ 1,2, ⋯, M : S o l v e x k - x i 2 + y k - y i 2 + h k - h i 2 - x k - x 1 2 + y k - y 1 2 + h k - h 1 2 = v t i k - t 1 k ∀ i v i a : θ n + 1 = θ n - α ∇ J θ n i f ∇ J > ε θ n - H - 1 ∇ J θ n o t h e r w i s e

（30）

以长征七号甲遥五运载火箭多级残骸为研究对象，应用TDOA-PSO、TDOA-NI-GD、TDOA-PSO-NI-GD三种不同的优化算法结果可视化呈现于图7。

如图8所示，经综合优化后，相较于单一算法优化，误差大幅下降，从接近2 km降至0.5 km以内，同时信号分离率显著提升，由68.5%升至96.2%。

2.3 整体模型检验

2.3.1 模型精度校验

由于实际情况存在随机误差，为修正误差所带来的影响，建立修正系数

τ

，其中

τ

是均值为0、标准差为0.5的正态分布。分析机器设备采集及现场勘探得到的数据，并检查监测设备，由于设备记录同时存在0.5 s的固定误差和随机误差，当所需最少设备中每个设备未精确识别出4个随机残骸在空中发生音爆时的位置和时间时，考虑在原有设备基础上增加一个设备，建立模型如式（31）所示。

m i n ∑ i = 1 4 ∑ j = 1 6 ∑ k = 1 7 t i j k - t k - d j k 340) 2 ⋅ a i j ∑ j = 1 6 x i j = 1, ∑ i = 1 4 x i j = 1 t i - t j ≤ 0.5 r 1 = (x 1 - x) 2 + (y 1 - y) 2 + (z 1 - z) 2 r i = (x i - x) 2 + (y i - y) 2 + (z i - z) 2 Δ r i = r i - r 1 = c Δ t i

（31）

式中 r₁为目标残骸到第1个设备的距离；r_i 为目标到其余各设备的距离；

Δ r i

为目标残骸到第i个设备和到第1个设备之间距离差；

Δ t i

为目标残骸到第

i

个设备与到第1个设备时间差测量值；

c ⋅ Δ t i

为第i个设备的距离差测量值，c为传播速度。

在考虑修正误差所带来的影响下，基于以上算法在原有的设备数上再增加一个设备，应用Chan-Taylor算法结合最小二乘法修正此误差。假设有k个监测设备，每个设备的位置已知为

x i, y i, z i

，其中i=1，2，3，…，k。对于每个设备，记录到音爆的时间t_i，但由于存在误差，实际时间应为

t i + τ

，从火箭残骸位置

x, y, z

到第i个监测设备位置

x i, y i, z i

的直线距离为r_i，根据音爆所到达时间来计算残骸到设备的距离，并更新位置和时间的估计值，以最小化产生的残差值，Chan-Taylor算法求解修正误差模型如式（32）所示：

r i = (x i - x) 2 + (y i - y) 2 + (z i - z) 2 r i = c * T + τ - t i e i = r i - c * T + τ - t i

（32）

式中 T为音爆产生的实际时间；t_i 是第i个监测设备记录到的音爆到达时间；

τ i

为该设备的时间误差；c为声速；

e i

为残差。

生成1个-0.5~0.5 s之间的随机偏移量，累加到原始的音爆抵达时间上。运用调整后的模型并借助Chan-Taylor算法进行求解。在7组样本数据中任意抽取6组，根据设备A、C、D、E、F、G模拟所需要的监测设备位置和音爆抵达时间数据的求解，结果如表7所示。通过优化的模型检验结果可知，标定为样本的4个残骸精准定位误差均小于1 km。

2.3.2 多示例验证

考虑到实际中多受山地地形与强风环境干扰，模拟相应环境以验证模型的完备性与准确性。

a）山地地形情景验证。山地地形考虑到高程维度变化，具体数据如表8所示，转换后坐标见表9。

山地地形下，应在欧氏距离中引入高程权重因子，其中α取平均值0.8：

d i = x i - x 2 + y i - y 2 + α h i - h 2

（33）

图9为山地地形下经TDOA-PSO-NI-GD多级优化后的误差可视化结果。

受山地地形影响，TDOA-PSO-NI-GD优化结果误差为0.611 7 km，略大于0.5 km，但仍显著小于TDOA-PSO（1.113 2 km）与TDOA-NI-GD（1.941 9 km）。

b）强风干扰情景验证。强风干扰环境下考虑到时间修正维度，数据详见表10。

转换后坐标见表11。

在强风干扰环境下，应在时间计算中引入水平漂移量

Δ x = v w i n d ⋅ t c o s θ

，修正后时间为

t i' = t i + Δ x v w i n d

。

强风干扰下TDOA-PSO-NI-GD多级优化后的误差可视化如图10所示。

显然，TDOA-PSO-NI-GD综合算法的定位优化误差（0.742 km）仍显著低于单级优化算法。以上实例有力验证了本模型的完备性与精确性。

3 结论

收起

本研究通过优化TDOA定位技术，提出了一种结合粒子群算法、NI-GD的多残骸定位方法，并以长征系列火箭残骸回收为案例验证了其有效性。得到以下结论：

a）对于单残骸定位，引入PSO和NI-GD的TDOA定位技术使单残骸定位误差被有效缩减至1 km范围内。此外，单残骸定位至少需要布置4台监测设备，以保证定位结果的稳定性和精确性，合理配置的设备数量可以实现监测效率提升和资源再利用的双重功能。

b）对于多残骸定位，相较于单残骸定位的TDOA基础模型，进一步引入信号分离技术和多音爆信号识别方法，显著提升了定位精度与信号分离率。此外，多残骸定位至少需要布置5台监测设备，以便有效应对信号重叠问题，保证模型在复杂环境下的定位精度。

c）本文模型能够通过灵敏度调节修正误差以确保定位结果的稳定性，对风速、地形变化等外部因素的抗干扰能力较强，具备良好的适应性。

综上所述，本文方法平均误差显著低于传统的结合单级优化算法的TDOA定位技术，在复杂地形和多变天气中表现出良好的定位精度及环境适应性。该模型不仅适用于火箭残骸的精准定位，还可推广到地质灾害监测、海上搜救、物流及供应链管理等场景，为复杂环境中的定位需求提供全面支持，展现出较强的实用性和应用前景。

基金

收起

国家社会科学基金项目(22BJY240)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

庞之浩. 回收火箭的重要意义与关键技术[J]. 科技导报, 2016, 34(1): 15-19.

PANG

Zhihao

. Importance and key technologies for recovering rockets[J]. Science & Technology Review, 2016, 34(1): 15-19.

[2]

张骁, 刘丙杰. 基于蒙特卡洛方法的火箭残骸落点范围预测[J]. 计算机测量与控制, 2023, 31(3): 102-109.

ZHANG

Xiao

, LIU

Bingjie

. Landing point range prediction for rocket wreckage based on monte carlo method[J]. Computer Measurement & Control, 2023, 31(3): 102-109.

[3]

Dexiu

, HUANG

, ZHANG

, et al. Joint TDOA, FDOA and differential Doppler rate estimation: method and its performance analysis[J]. Chinese Journal of Aeronautics, 2018, 31(1): 137-147.

[4]

Wanchun

, CHEN

Ruibin

, GUO

Yuning

, et al. Closed form algorithm of double-satellite TDOA+AOA localization based on WGS-84 model[J]. Chinese Journal of Aeronautics, 2019, 32(10): 2354-2367.

[5]

DELCOURT

, BOUDEC

J Y L

. TDOA source-localization technique robust to time-synchronization attacks[J]. IEEE Transac-tions on Information Forensics and Security, 2021(16): 4249-4264.

[6]

KIRMAZ

, ŞAHIN

, MICHALOPOULOS

D S

, et al. TOA and TDOA estimation using artificial neural networks for high-accuracy ranging[J]. IEEE Journal on Selected Areas in Communications, 2023, 41(12): 3816-3830.

[7]

, LIU

G P

. A robust high-accuracy ultrasound indoor positioning system based on a wireless sensor network[J]. Sensors, 2017, 17(11): 2554.

[8]

PATE

B L

. Identifying and tracking disaster victims: state-of-the-art technology review[J]. Family & Community Health, 2008, 31(1): 23-34.

[9]

LYU J, ZENG

, ZHANG

, et al. Placement optimization of UAV-mounted mobile base stations[J]. IEEE Communications Letters, 2016, 21(3): 604-607.

[10]

CHEN

, CHEN

, CAO

, et al. Acoustic indoor localization system integrating TDMA+FDMA transmission scheme and positioning correction technique[J]. Sensors, 2019, 19(10): 2353.

[11]

, XIE

, TAN

. Iterative constrained weighted least squares source localization using TDOA and FDOA measurements[J]. IEEE Transactions on Signal Processing, 2017, 65(15): 3990-4003.

[12]

DELCOURT

, BOUDEC

J Y L

. TDOA source-localization technique robust to time-synchronization attacks[J]. IEEE Transac-tions on Information Forensics and Security, 2020(16): 4249-4264.

[13]

田昕, 魏国亮, 王甘楠. 无线传感器网络定位综述[J]. 信息与控制, 2022, 51(1): 69-87.

TIAN

Xin

, WEI

Guoliang

, WANG

Gannan

. Review of wireless sensor network localization[J]. Information and Control, 2022, 51(1): 69-87.

[14]

, LIU

, XIAO

. Location of natural gas leakage sources on offshore platform by a multi-robot system using particle swarm optimization algorithm[J]. Journal of Natural Gas Science and Engineering, 2020(84): 103636.

[15]

PHUNG

M D

, HA

Q P

. Safety-enhanced UAV path planning with spherical vector-based particle swarm optimization[J]. Applied Soft Computing, 2021(107): 107376.

[16]

HERBAZI

, AMECHNOUE

, KHOUYA

, et al. Investigation of photovoltaic output characteristics with iterative methods[J]. Proce-dia Manufacturing, 2019(32): 794-801.

[17]

刘兰冬, 刘铭. 一类二次矩阵方程的牛顿迭代法及其收敛性[J]. 工程数学学报, 2024, 41(3): 587-594.

LIU

Landong

, LIU

Ming

. Newton iterative methods for a class of quadratic matrix equations and its convergence[J]. Chinese Journal of Engineering Mathematics, 2024, 41(3): 587-594.

[18]

WANG

, TAN

, LIU

. Particle swarm optimization algorithm: an overview[J]. Soft computing, 2018, 22(2): 387-408.

[19]

CHUN

. Iterative methods improving Newton’s method by the decomposition method[J]. Computers & Mathematics with Applica-tions, 2005, 50(10-12): 1559-1568.

[20]

BADAWI

, BASSI

, OZEV

, et al. Deep-learning-based gas leak source localization from sparse sensor data[J]. IEEE Sensors Journal, 2022, 22(21): 20999-21008.

[21]

中国国家航天局/CNSA. 中国国家航天局[CNSA]-Organizations-Roberts Space Industries. (n.d.)[EB/OL]. (2024-12-27)[2025-02-13] https://robertsspaceindustries.com/orgs/CNSA. https://robertsspaceindustries.com/orgs/CNSA

[22]

WANG

, CAI

, LI

, et al. A bias-reduced nonlinear WLS method for TDOA/FDOA-based source localization[J]. IEEE Transactions on Vehicular Technology, 2015, 65(10): 8603-8615.

[23]

ZHANG

, LU

, CHE

, et al. Moving-distance-minimized PSO for mobile robot swarm[J]. IEEE Transactions on Cybernetics, 2021, 52(9): 9871-9881.

[24]

万邦银, 蒯念生, 何雄元, 等. 免疫粒子群算法在修正高斯模型下的源强反演[J]. 中国安全科学学报, 2024, 34(7): 132-138.

WAN

Bangyin

, KUAI

Niansheng

, HE

Xiongyuan

, et al. Source strength inversion of PSO-IA under modified Gaussian models[J]. China Safety Science Journal, 2024, 34(7): 132-138.

[25]

MISTRY

, ZHANG

, NEOH

S C

, et al. A micro-GA embedded PSO feature selection approach to intelligent facial emotion reco-gnition[J]. IEEE Transactions on Cybernetics, 2016, 47(6): 1496-1509.

[26]

, ZHAO

, CUI

, et al. DC cascaded energy storage system based on DC collector with gradient descent method[J]. IEEE Transactions on Industrial Electronics, 2023, 71(2): 1594-1605.

[27]

MISHCHENKO

. Regularized Newton method with global convergence[J]. SIAM Journal on Optimization, 2023, 33(3): 1440-1462.

[28]

路荣, 张高鹏, 曹剑中, 等. 基于单目视觉的火箭回收高度测量技术研究[J]. 光学精密工程, 2024, 32(14): 2166-2188.

Rong

, ZHANG

Gaopeng

, CAO

Jianzhong

, et al. Research on measurement technology of rocket recovery height based on monocular vision[J]. Optics and Precision Engineering, 2024, 32(14): 2166-2188.

2025年第48卷第5期

PDF下载

590

288

引用本文

BibTeX

文章信息

doi: 10.7654/j.issn.2097-1974.20250508

接收时间：2024-12-27
首发时间：2025-11-27
出版时间：2025-10-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2024-12-27
修回日期：2025-02-13

基金

国家社会科学基金项目(22BJY240)

作者信息

东北农业大学工程学院，哈尔滨，150030

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/ddyht/CN/10.7654/j.issn.2097-1974.20250508

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

设备	经度/（°）	纬度/（°）	高程/m	音爆抵达时间/s
A	110.241	27.204	824	100.767
B	110.780	27.456	727	112.220
C	110.712	27.785	742	188.020
D	110.251	27.825	850	258.985
E	110.524	27.617	786	118.443
F	110.467	27.921	678	266.871
G	110.047	27.121	575	163.024

设备

经度/（°）

纬度/（°）

高程/m

音爆抵达时间/s

110.241

27.204

824

100.767

110.780

27.456

727

112.220

110.712

27.785

742

188.020

110.251

27.825

850

258.985

110.524

27.617

786

118.443

110.467

27.921

678

266.871

110.047

27.121

575

163.024

设备	经度 /（°）	纬度 /（°）	高程/m	音爆抵达时间/s
A	110.241	27.204	824	100.767	164.229	214.850	270.065
B	110.783	27.456	727	92.453	112.220	169.36	196.583
C	110.762	27.785	742	75.560	110.696	156.936	188.020
D	110.251	28.025	850	94.653	141.409	196.517	258.985
E	110.524	27.617	786	78.600	86.216	118.443	126.669
F	110.467	28.081	678	67.274	166.270	175.482	266.871
G	110.047	27.521	575	103.738	163.024	206.789	210.306

设备

经度

/（°）

纬度

/（°）

高程/m

音爆抵达时间/s

110.241

27.204

824

100.767

164.229

214.850

270.065

110.783

27.456

727

92.453

112.220

169.36

196.583

110.762

27.785

742

75.560

110.696

156.936

188.020

110.251

28.025

850

94.653

141.409

196.517

258.985

110.524

27.617

786

78.600

86.216

118.443

126.669

110.467

28.081

678

67.274

166.270

175.482

266.871

110.047

27.521

575

103.738

163.024

206.789

210.306

设备	经度/（°）	纬度/（°）	高程/m	音爆抵达时间/s
H	109.850	19.620	102	85.323
I	110.200	19.910	88	92.152
J	110.550	20.150	75	105.671
K	110.002	20.430	120	116.841

设备

经度/（°）

纬度/（°）

高程/m

音爆抵达时间/s

109.850

19.620

102

85.323

110.200

19.910

92.152

110.550

20.150

105.671

110.002

20.430

120

116.841

设备	经度/（°）	纬度/（°）	高程/m	残骸1音爆抵达时间/s	残骸2音爆抵达时间/s
L	112.400	25.650	95	76.214	82.451
M	112.732	25.900	85	81.341	88.920
N	113.000	26.150	105	85.672	92.112
O	112.200	25.450	110	90.120	96.342
P	112.500	25.750	90	94.563	98.781

设备

经度/（°）

纬度/（°）

高程/m

残骸1音爆抵达时间/s

残骸2音爆抵达时间/s

112.400

25.650

76.214

82.451

112.732

25.900

81.341

88.920

113.000

26.150

105

85.672

92.112

112.200

25.450

110

90.120

96.342

112.500

25.750

94.563

98.781

残骸	x/km	y/km	h/km
长征三号	17.358	81.265	0.910
长征五号B遥三	44.290	55.743	0.105
长征七号甲遥五残骸1	25.016	25.924	0.090
长征七号甲遥五残骸2	46.205	9.346	0.095

残骸

x/km

y/km

h/km

长征三号

17.358

81.265

0.910

长征五号B遥三

44.290

55.743

0.105

长征七号甲遥五残骸1

25.016

25.924

0.090

长征七号甲遥五残骸2

46.205

9.346

0.095

设备	x/km	y/km	h/km	音爆抵达时间/s
A	18.884	9.235	0.824	100.767
B	71.350	37.273	0.727	112.220
C	64.731	73.879	0.742	188.020
D	19.857	78.329	0.850	258.985
E	46.431	55.186	0.786	118.443
F	40.883	89.010	0.678	266.871
G	0	0	0.575	163.024
H	0	0	0.102	85.323
I	34.069	32.266	0.088	92.152
J	68.138	58.969	0.075	105.671
K	14.796	90.123	0.120	116.841
设备	x/km	y/km	h/km	残骸1音爆抵达时间/s	残骸2音爆抵达时间/s
L	19.468	22.253	0.095	76.214	82.451
M	51.785	50.068	0.085	81.341	88.920
N	77.872	77.884	0.105	85.672	92.112
O	0	0	0.110	90.120	96.342
P	29.202	33.379	0.090	94.563	98.781

设备

x/km

y/km

h/km

音爆抵达时间/s

18.884

9.235

0.824

100.767

71.350

37.273

0.727

112.220

64.731

73.879

0.742

188.020

19.857

78.329

0.850

258.985

46.431

55.186

0.786

118.443

40.883

89.010

0.678

266.871

0.575

163.024

0.102

85.323

34.069

32.266

0.088

92.152

68.138

58.969

0.075

105.671

14.796

90.123

0.120

116.841

设备

x/km

y/km

h/km

残骸1音爆抵达时间/s

残骸2音爆抵达时间/s

19.468

22.253

0.095

76.214

82.451

51.785

50.068

0.085

81.341

88.920

77.872

77.884

0.105

85.672

92.112

0.110

90.120

96.342

29.202

33.379

0.090

94.563

98.781

残骸类别	修正前	修正后	距离误差/km
1	（37.121，60.844，0.778 3）	81.041	（37.473，61.752，0.779 6）	80.854	0.34
2	（29.445，50.015，0.757 1）	75.274	（29.521，49.314，0.765 9）	77.243	0.41
3	（17.856，43.978，0.701 1）	75.703	（17.735，44.833，0.687 4）	81.572	0.27
4	（27.337，49.626，0.353 1）	77.016	（27.128，48.992，0.385 7）	82.052	0.65

残骸类别

修正前

修正后

距离误差/km

音爆发生位置/km

音爆发生时间/s

音爆发生位置/km

音爆发生时间/s

（37.121，60.844，0.778 3）

81.041

（37.473，61.752，0.779 6）

80.854

0.34

（29.445，50.015，0.757 1）

75.274

（29.521，49.314，0.765 9）

77.243

0.41

（17.856，43.978，0.701 1）

75.703

（17.735，44.833，0.687 4）

81.572

0.27

（27.337，49.626，0.353 1）

77.016

（27.128，48.992，0.385 7）

82.052

0.65

设备	经度/（°）	纬度/（°）	高程/km	音爆抵达时间/s
Q	99.450	27.890	1.500	145.323
R	99.780	28.120	0.980	152.670
S	99.692	28.350	2.200	138.451
T	99.600	28.250	1.800	142.890
实际残骸位置	99.875	28.212	1.200	—

设备

经度/（°）

纬度/（°）

高程/km

音爆抵达时间/s

99.450

27.890

1.500

145.323

99.780

28.120

0.980

152.670

99.692

28.350

2.200

138.451

99.600

28.250

1.800

142.890

实际残骸位置

99.875

28.212

1.200

—

设备	x/km	y/km	h/km	音爆抵达时间/s
Q	0	0	1.500	145.323
R	32.122	25.590	0.980	152.670
S	23.56	51.181	2.200	138.451
T	14.601	40.055	1.800	142.890
实际残骸位置	41.369	35.827	1.200	—

设备

x/km

y/km

h/km

音爆抵达时间/s

1.500

145.323

32.122

25.590

0.980

152.670

23.56

51.181

2.200

138.451

14.601

40.055

1.800

142.890

实际残骸位置

41.369

35.827

1.200

—

设备	经度/（°）	纬度/（°）	高程/ km	音爆抵达时间/s	修正后时间/s
U	115.230	30.550	0.080	88.760	90.212
V	115.622	31.850	0.095	92.342	93.790
W	116.000	31.105	0.105	96.451	97.903
X	115.400	30.700	0.090	89.120	90.571
实际残骸位置	115.805	30.938	0.093	—	—

设备

经度/（°）

纬度/（°）

高程/

音爆抵达

时间/s

修正后

时间/s

115.230

30.550

0.080

88.760

90.212

115.622

31.850

0.095

92.342

93.790

116.000

31.105

0.105

96.451

97.903

115.400

30.700

0.090

89.120

90.571

实际残骸位置

115.805

30.938

0.093

—

设备	x/km	y/km	h/km	音爆抵达时间/s	修正后时间/s
U	0	0	0.080	88.760	90.212
V	38.157	144.642	0.095	92.342	93.790
W	74.952	61.751	0.105	96.451	97.903
X	16.548	16.689	0.090	89.120	90.571
实际残骸位置	55.970	43.170	0.090	—	—

设备

x/km

y/km

h/km

音爆抵达时间/s

修正后时间/s

0.080

88.760

90.212

38.157

144.642

0.095

92.342

93.790

74.952

61.751

0.105

96.451

97.903

16.548

16.689

0.090

89.120

90.571

实际残骸位置

55.970

43.170

0.090

—