TOA的定位,建模与解算的步骤、公式推导

在这里插入图片描述

TOA(到达时间)定位的核心是通过测量信号从目标到多个基站的传播时间,将其转换为距离信息,并利用几何关系解算目标位置。本文给出具体的建模与解算步骤及公式推导

通用模型建立

假设条件

  • 已知 N N N个基站坐标 ( x i , y i ) (x_i, y_i) (xi,yi) i = 1 , 2 , . . . , N i=1,2,...,N i=1,2,...,N)。
  • 目标位置为未知参数 ( x , y ) (x, y) (x,y)
  • 信号传播速度为光速 c c c(空中场景通常取 c = 3 × 1 0 8   m/s c=3 \times 10^8 \, \text{m/s} c=3×108m/s)。

距离模型
信号到达时间 t i t_i ti对应的距离为:
r i = c ⋅ t i = ( x − x i ) 2 + ( y − y i ) 2 + n i r_i = c \cdot t_i = \sqrt{(x - x_i)^2 + (y - y_i)^2} + n_i ri=cti=(xxi)2+(yyi)2 +ni
其中 n i n_i ni为测量噪声(通常假设为高斯白噪声),实际距离为 d i = ( x − x i ) 2 + ( y − y i ) 2 d_i = \sqrt{(x - x_i)^2 + (y - y_i)^2} di=(xxi)2+(yyi)2

非线性方程组构建

每个基站对应一个圆方程:
( x − x i ) 2 + ( y − y i ) 2 = r i 2 ( i = 1 , 2 , . . . , N ) (x - x_i)^2 + (y - y_i)^2 = r_i^2 \quad (i=1,2,...,N) (xxi)2+(yyi)2=ri2(i=1,2,...,N)
由于噪声存在,方程需通过优化方法求解目标位置 ( x , y ) (x, y) (x,y)

线性化处理(最小二乘法)

将非线性方程转化为线性方程组以减少计算复杂度:

  1. 展开方程
    对第$i)个基站:
    [
    x^2 - 2x_i x + x_i^2 + y^2 - 2y_i y + y_i^2 = r_i^2
    ]

  2. 消去二次项
    选取第1个基站为参考,其余方程减去参考方程:
    − 2 ( x i − x 1 ) x − 2 ( y i − y 1 ) y = r i 2 − r 1 2 − ( x i 2 + y i 2 − x 1 2 − y 1 2 ) -2(x_i - x_1)x - 2(y_i - y_1)y = r_i^2 - r_1^2 - (x_i^2 + y_i^2 - x_1^2 - y_1^2) 2(xix1)x2(yiy1)y=ri2r12(xi2+y

在脉冲星导航(XNAV)中,模型预测的TOA(Time of Arrival)是导航系统中用于计位置的关键参数。为了获得模型预测的TOA,通常需要结合脉冲星的射电或X射线辐射特性、轨道动力学模型以及观测数据进行建模和预测。以下是主要的计方法和相关软件工具: ### 脉冲星导航中模型预测TOA的计方法 1. **脉冲星计时模型** 模型预测TOA的核心在于建立高精度的脉冲星计时模型。该模型通常包括脉冲星自转参数(如自转频率、频率导数)、轨道参数(如脉冲星太阳系质心之间的相对运动)以及星际介质延迟(如色散延迟)。通过这些参数,可以预测脉冲星信号到达太阳系质心的理论TOA。 在实际导航应用中,还需考虑航天器相对于太阳系质心的运动,利用相对论效应(如引力时间延迟、多普勒频移)修正TOA预测值。 2. **色散延迟修正** 脉冲星信号在传播过程中会受到星际介质的影响,导致不同频率的信号以不同速度传播,产生色散延迟。通过测量不同频段的TOA差异,结合色散量(DM, Dispersion Measure)可以修正这一延迟。 3. **相对论时间延迟模型** 航天器太阳系质心之间的相对运动会导致相对论时间延迟,包括夏皮罗延迟(Shapiro delay)和引力红移效应。这些效应需要通过广义相对论模型进行计,并对TOA进行相应修正。 4. **TOA残差分析** 模型预测TOA实际观测TOA之间的差值称为TOA残差。通过最小化TOA残差,可以优化脉冲星参数模型,并进一步提升导航精度。通常使用最小二乘法或最大似然估计方法进行残差拟合。 ### 常用软件工具 1. **TEMPO2** TEMPO2 是一个广泛用于脉冲星计时分析的软件包,能够处理高精度TOA数据并建立脉冲星计时模型。它支持多种脉冲星参数模型(如DE430、DE436等行星历表)以及相对论修正模型,适用于XNAV中的TOA预测任务。 2. **PINT (PINT is Not TEMPO)** PINT 是一个基于Python的现代脉冲星计时分析工具,旨在提供更灵活和可扩展的TOA建模能力。它支持TEMPO2兼容的模型文件,并能够处理复杂的相对论效应和星际介质延迟修正。 3. **XSPEC/XRONOS(用于X射线数据分析)** 在X射线脉冲星导航中,常使用HEASOFT软件包中的XSPEC和XRONOS工具进行光谱和时间序列分析。XRONOS可用于提取X射线脉冲星的TOA,而XSPEC用于分析光谱特性以辅助TOA建模。 4. **NavXNAV** NavXNAV 是专为X射线脉冲星导航设计的仿真分析工具,集成了TOA预测、轨道动力学建模和导航滤波法。该工具可用于生成模型预测TOA并进行导航性能评估。 ### 示例代码(使用PINT计模型预测TOA) ```python from pint.models import get_model from pint.toa import get_TOAs from pint.fitter import Fitter # 加载脉冲星计时模型文件(如B1855+09.par) model = get_model("B1855+09.par") # 加载观测TOA数据文件(如B1855+09.tim) toas = get_TOAs("B1855+09.tim") # 使用模型预测TOA predicted_toas = model.predict_toas(toas) # 拟合模型参数并优化TOA预测 fitter = Fitter(toas, model) fitter.fit() # 输出优化后的模型预测TOA print(predicted_toas) ``` ###
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

MATLAB卡尔曼

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值