[发明专利]一种侦察卫星定位估计方法

摘要

本发明提出了一种侦察卫星定位估计方法，所述方法利用北斗卫星获取目标侦察卫星在零时刻的初始位置和初始速度，建立目标侦察卫星的简化运动方程，并将目标侦察卫星的简化运动方程进行分解，建立能确定目标侦查卫星估计位置的非线性微分方程模型；对所建立的非线性微分方程模型进行求解，得出需观测的侦察卫星在各个时间点的估计位置，进行定位估计。本发明能对侦察卫星进行更精确、更可靠、更自主地定位估计，能为判断与掌握目标侦察卫星的飞行意图、飞行状态及飞行趋势提供更好的信息基础，该方法简单易行，实用性极强。

主权项

1.一种侦察卫星定位估计方法,其特征在于，包括如下步骤：步骤A，利用北斗卫星获取目标侦察卫星零时刻在基础坐标系下的初始位置和初始速度：x_c(0)=x₀,y_c(0)=y₀,z_c(0)=z₀x·c(0)=x00,y·c(0)=y00,z·c(0)=z00]]>其中，x_c(0)、y_c(0)、z_c(0)表示零时刻目标侦察卫星在基础坐标系下的初始位置；表示零时刻侦察卫星在基础坐标系下的初始速度；x₀、y₀、z₀为利用北斗卫星获取的目标侦察卫星在零时刻在基础坐标系下的初始位置值，x₀₀、y₀₀、z₀₀为利用北斗卫星获取的目标侦察卫星在零时刻在基础坐标系下的初始速度值；步骤B，建立目标侦察卫星在基础坐标系下的简化运动方程：r→··c(t)=F→e=-Gm|r→c(t)|3r→c(t)rc(t)=xc2(t)+yc2(t)+zc2(t)Gm=3.986005×1014]]>其中，为飞行器所受的外力加速度之和，G_m为地球引力常数，为在基础坐标系下的位置矢量，为对时间t的二阶导数，即加速度，r_c(t)为的绝对值，即x_c(t)、y_c(t)、z_c(t)为t时刻目标侦察卫星在基础坐标系下的位置；步骤C，将目标侦察卫星的简化运动方程进行分解，分解为以下方程组：x··c(t)=-Gmrc(t)3xc(t)y··c(t)=-Gmrc(t)3yc(t)z··c(t)=-Gmrc(t)3zc(t)]]>其中，为t时刻目标侦察卫星在基础坐标系下的加速度；步骤D，将目标侦察卫星初始位置和初始速度的方程结合分解的简化运动方程，建立确定目标侦察卫星的估计位置的非线性微分方程模型：x··c(t)=-Gmrc(t)3xc(t)y··c(t)=-Gmrc(t)3yc(t)z··c(t)=-Gmrc(t)3zc(t)Gm=3.986005×1014rc(t)=xc2(t)+yc2(t)+zc2(t)xc(0)=x0,yc(0)=y0,zc(0)=z0x·c(0)=x00,y·c(0)=y00,z·c(0)=z00]]>步骤E，对所建立的非线性微分方程模型进行求解，得出目标侦察卫星在各个时间点在基础坐标系下的估计位置，从而实现定位估计。