首页> 中国专利> 一种编队飞行卫星绝对和相对轨道确定方法

一种编队飞行卫星绝对和相对轨道确定方法

页面导航

摘要
著录项
法律信息
说明书
相似文献

摘要

一种编队飞行卫星绝对和相对轨道确定方法，针对星间存在相对距离测量设备的编队飞行任务，通过对星间距离的相对测量值进行解算，获得编队双星的绝对轨道根数和相对轨道根数。由于解算过程仅需要相对距离的测量值，且无需第三方测量信号源，因此星间距离的接收方无需地面测控站以及GPS等辅助设备支持，即可单独确定此类编队飞行任务的绝对和相对轨道根数。

著录项

公开/公告号CN102322862A

专利类型发明专利
公开/公告日2012-01-18

原文格式PDF
申请/专利权人航天东方红卫星有限公司;
展开▼

申请/专利号CN201110182467.9
发明设计人谭田;徐明;李延东;陶成华;蒙薇;
展开▼

申请日2011-06-29
分类号G01C21/24(20060101);
代理机构11009 中国航天科技专利中心;
代理人安丽
地址 100094 北京市5616信箱
入库时间 2023-12-18 04:21:34

法律信息

法律状态公告日

法律状态信息

法律状态
2013-07-24

授权

授权
2012-03-14

实质审查的生效 IPC(主分类):G01C21/24 申请日:20110629

实质审查的生效
2012-01-18

公开

公开

说明书

技术领域

本发明涉及一种编队飞行卫星的轨道确定方法。

背景技术

编队飞行是未来小卫星领域的重要发展方向，多颗卫星以一定的组织形式完成公里至百公里级的集群飞行，彼此之间通过信息交换来共同实现对地观测以及目标定位等任务。

为了构成稳定有效的测量基线，各组成卫星需要基于相对状态来完成相对状态的控制，而相对状态的解算精度将是决定编队飞行任务成败的重要因素。

现有编队飞行卫星的绝对轨道确定，往往需要借助地面测站以及以GPS 和GLONASS为代表的全球导航卫星系统等第三方信号源实现，而相对轨道或相对位置确定，大多采用测距或测角雷达直接获取相对状态或采用GPS(或 GLONASS)进行差分解算。上述研究方法对测量设备要求较高：或者需要借助第三方测量设备，或者要求相对测量设备提供测距和测角等功能；而对于仅依靠相对距离测量进行相对以及绝对轨道解算，尚无相关研究。

发明内容

本发明的技术解决问题是：克服现有技术的不足，提供了一种基于星间距离且不引入第三方测量信息的编队飞行卫星相对和绝对轨道确定方法。

本发明的技术解决方案是：一种编队飞行卫星绝对和相对轨道确定方法，步骤如下：

(1)获取待解算的两颗卫星在初始时刻t₀的瞬时轨道根数估计值，记为σ_a0和σ_b0，a和b分别代表两颗卫星；

(2)根据初始时刻的瞬时轨道根数估计值σ_a0和σ_b0，分别计算两颗卫星在指定时刻t_i的瞬时轨道根数，记为σ_ai和σ_bi，i＝1，2，3，...，n；

(3)根据两颗卫星在指定时刻t_i的瞬时轨道根数σ_ai和σ_bi，分别计算惯性坐标系下两颗卫星的绝对位置的笛卡儿坐标分量[x_ai y_ai z_ai]^T和[x_bi y_bi z_bi]^T；

(4)根据两颗卫星在惯性坐标系下绝对位置的笛卡儿坐标分量，得到在各指定时刻两颗卫星的星间距离计算值，以及该星间距离计算值与测量真值的偏差δL_i；

(5)计算卫星a在初始时刻t₀的瞬时轨道根数对星间距离计算值的状态转移矩阵；

(6)计算卫星b在初始时刻t₀的瞬时轨道根数对星间距离计算值的状态转移矩阵；

(7)根据星间相对距离的计算值与测量真值的偏差δL_i，利用步骤(5)和步骤(6)得到的两个状态转移矩阵计算卫星a和卫星b在初始时刻t₀瞬时轨道根数估计值的修正量δσ_a0和δσ_b0：

(8)将σ_a0-δσ_a0和σ_b0-δσ_b0分别作为卫星a和卫星b在初始时刻t₀瞬时轨道根数的新的估计值，重复步骤(1)～(7)的计算过程进行迭代，直至两颗卫星的半长轴修正量δa_a0和δa_b0满足定位精度要求；

(9)利用满足定位精度要求的卫星a和卫星b在初始时刻t₀瞬时轨道根数的迭代结果σ_a0和σ_b0，按照第(2)步的计算方法得到两颗卫星在任意时刻的瞬时轨道根数σ_a和σ_b，从而得到相对瞬时轨道根数Δσ＝σ_b-σ_a。

本发明与现有技术相比的优点在于：本发明方法针对星间存在相对距离测量设备的编队飞行任务，通过对星间距离的相对测量值进行解算，可以获得编队双星的绝对和相对轨道根数。由于解算过程仅需要相对距离的测量，且无需第三方测量信号源，因此星间距离的接收方无需地面测控站以及GPS等辅助设备支持，即可以单独确定此类编队飞行任务的绝对和相对轨道根数。方法简便、易于实现。

附图说明

图1为本发明方法的流程框图；

图2为本发明实施例的星间距离真值变化图。

具体实施方式

如图1所示，为本发明方法的流程框图，主要步骤如下：

(1)输入待解算的两颗卫星(卫星a和卫星b)在初始时刻(t＝t₀)的瞬时轨道根数的估计值，记为σ_a0和σ_b0。本发明方法中涉及到的轨道根数，均采用适用于近圆轨道的无奇点轨道根数，即半长轴a、偏心率矢量e_x和e_y、倾角 i、升交点赤经Ω、平纬度幅角λ；

(2)根据初始时刻的瞬时轨道根数估计值σ_a0和σ_b0，分别计算指定时刻 t_i的瞬时轨道根数，记为σ_ai和σ_bi，i＝1，2，3，...，n；计算时根据考虑J₂、J₃、J₄和J₅等非球型引力项的摄动微分方程进行数值求解：

$(\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) = (\begin{matrix} 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ \sqrt{\frac{μ}{a^{3}}} \end{matrix}) + (\begin{matrix} \frac{{2 a}^{2}}{\sqrt{μp}} (e_{x} \sin u - e_{y} \cos u) & \frac{{2 a}^{2}}{\sqrt{μp}} (1 + e_{x} \cos u + e_{y} \sin u) & 0 \\ 0 & 0 & \frac{r \sin u}{\sqrt{μp} \sin i} \\ 0 & 0 & \frac{r \cos u}{\sqrt{μp}} \\ \sqrt{\frac{p}{μ}} \sin u & \sqrt{\frac{p}{μ}} [(1 + \frac{r}{p}) \cos u + \frac{e_{x} r}{p}] & \sqrt{\frac{p}{μ}} \frac{e_{y} r}{p} \frac{\sin u}{\tan i} \\ - \sqrt{\frac{p}{μ}} \cos u & \sqrt{\frac{p}{μ}} [(1 + \frac{r}{p}) \sin u + \frac{e_{y} r}{p}] & - \sqrt{\frac{p}{μ}} \frac{e_{x} r}{p} \frac{\sin u}{\tan i} \\ \sqrt{\frac{p}{μ}} [- \frac{1}{2} (e_{x} \cos u + e_{y} \sin u) - \frac{2 r}{p} \sqrt{1 - e_{x}^{2} - e_{y}^{2}}] & \frac{1}{2} \sqrt{\frac{p}{μ}} (e_{x} \sin u - e_{y} \cos u) & - \sqrt{\frac{p}{μ}} \frac{r}{p} \frac{\sin u}{\tan i} \end{matrix}) (\begin{matrix} f_{r} \\ f_{u} \\ f_{h} \end{matrix})$

其中，μ为引力系数(＝3.986005×10¹⁴m³/s²)，其中u为纬度幅角，p和r为分别为半通径和地心距，即

u＝λ+2e_xsinλ-2e_ycosλ

$p = a (1 - e_{x}^{2} - e_{y}^{2})$

$r = \frac{a (1 - e_{x}^{2} - e_{y}^{2})}{1 + e_{x} \cos λ + e_{y} \sin λ}$

摄动加速度在径向、切向和法向的分量为

f_r＝Δg_r

f_u＝Δg_mcosζ+Δg_psinζ

f_h＝Δg_msinζ-Δg_pcosζ

其中，ζ可由球面三角学得到，即

sinφ＝sinλsin i

cosζ＝ctanλtanφ

其中，φ为地心纬度，ctan为余切函数。

摄动加速度在的分量为在径向、子午向和纬线向分量为

$Δ g_{r} = \frac{\partial ΔU}{\partial r}$

$Δ g_{m} = \frac{1}{r} \frac{\partial ΔU}{\partial φ}$

$Δ g_{p} = \frac{1}{r \cos φ} \frac{\partial ΔU}{\partial λ}$

前5阶带形谐函数在径向、子午向和纬线向分量为：

$\frac{\partial ΔU}{\partial r} = \frac{μ}{r^{2}} [\frac{3}{2} J_{2} {(\frac{R_{E}}{r})}^{2} (3 \sin^{2} φ - 1)$

$+ 2 J_{3} {(\frac{R_{E}}{r})}^{3} (5 \sin^{3} φ - 3 \sin φ)$

$+ \frac{5}{2} J_{4} {(\frac{R_{E}}{r})}^{4} (35 \sin^{4} φ - 30 \sin^{2} φ + 3)$

$+ 3 J_{5} {(\frac{R_{E}}{r})}^{5} (63 \sin^{5} φ - 70 \sin^{3} φ + 15 \sin φ)]$

$\frac{\partial ΔU}{\partial r} = - \frac{μ}{r} \cos φ [\frac{1}{2} J_{2} {(\frac{R_{E}}{r})}^{2} (6 \sin φ)$

$+ \frac{1}{2} J_{3} {(\frac{R_{E}}{r})}^{3} (15 \sin^{2} φ - 3)$

$+ \frac{1}{8} J_{4} {(\frac{R_{E}}{r})}^{4} (140 \sin^{3} φ - 60 \sin φ)$

$+ \frac{1}{8} J_{5} {(\frac{R_{E}}{r})}^{5} (325 \sin^{4} φ - 210 \sin^{2} φ + 15)]]$

$\frac{\partial ΔU}{\partial λ} = 0$

其中J₂、J₃、J₄、J₅带形谐系数，Re为地球赤道半径。

(3)根据上步计算得到的指定时刻t_i的瞬时轨道根数σ_ai和σ_bi，分别计算惯性坐标系下卫星a和卫星b的绝对位置的笛卡儿坐标分量[x_ai y_ai z_ai]^T和 [x_bi y_bi z_bi]^T，计算方法如下：

$(\begin{matrix} x \\ y \\ z \end{matrix}) = (\begin{matrix} \cos u & \sin u & 0 \\ - \sin u & \cos u & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} 1 & 0 & 0 \\ 0 & \cos i & \sin i \\ 0 & - \sin i & \cos i \end{matrix}) (\begin{matrix} \cos Ω & \sin Ω & 0 \\ - \sin Ω & \cos Ω & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} r \\ 0 \\ 0 \end{matrix})$

其中，惯性坐标系的坐标原点在地心，x轴由地心指向平春分点，z轴由地心指向平赤道法线方向，y轴由右手坐标法则确定。

(4)根据卫星a和卫星b在惯性坐标系下的位置坐标分量，得到在各指定时刻的星间距离计算值 $L_{i} = \sqrt{{(x_{bi} - x_{ai})}^{2} + {(y_{bi} - y_{ai})}^{2} + {(z_{bi} - z_{ai})}^{2}};$ 该计算值与测量真值的偏差记为

(5)计算卫星a在初始时刻t₀瞬时轨道根数对星间距离计算值的状态转移矩阵由于星间距离L_i是x_ai、y_ai、z_ai的函数，而x_ai、y_ai、z_ai分别是t_i时刻瞬时轨道根数σ_ai的函数，σ_ai是t₀时刻瞬时轨道根数σ_a0函数，则根据上述复合函数传递规律，可得

$\frac{\partial L_{i}}{\partial σ_{a 0}} = \frac{\partial L_{i}}{\partial x_{ai}} \frac{\partial x_{ai}}{\partial σ_{ai}} \frac{\partial σ_{ai}}{\partial σ_{a 0}} + \frac{\partial L_{i}}{\partial y_{ai}} \frac{\partial y_{ai}}{\partial σ_{ai}} \frac{\partial σ_{ai}}{\partial σ_{a 0}} + \frac{\partial L_{i}}{\partial z_{ai}} \frac{\partial z_{ai}}{\partial σ_{ai}} \frac{\partial σ_{ai}}{\partial σ_{a 0}}$

其中，相对距离对各坐标分量的偏导数如下：

$\frac{\partial L_{i}}{\partial x_{ai}} = \frac{x_{ai} - x_{bi}}{L_{i}}$

$\frac{\partial L_{i}}{\partial y_{ai}} = \frac{y_{ai} - y_{bi}}{L_{i}}$

$\frac{\partial L_{i}}{\partial z_{ai}} = \frac{z_{ai} - z_{bi}}{L_{i}}$

卫星在惯性坐标系下的位置坐标分量(x_ai、y_ai、z_ai)、指定时刻的瞬时轨道根数 σ_ai均可视为初始时刻瞬时轨道根数σ_a0关于时间的函数。初始时刻瞬时轨道根数的微小改变量引起星间距离的变化量可由得到。因此，通过设置合适的σ_a0改正量，可将L_i改正-δL_i，即达到期望值

坐标分量对各轨道根数的偏导数采用下面的方法求取：

a.坐标分量对a的偏导数

$\frac{\partial}{\partial a} (\begin{matrix} x \\ y \\ z \end{matrix}) = (\begin{matrix} \cos u & \sin u & 0 \\ - \sin u & \cos u & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} 1 & 0 & 0 \\ 0 & \cos i & \sin i \\ 0 & - \sin i & \cos i \end{matrix}) (\begin{matrix} \cos Ω & \sin Ω & 0 \\ - \sin Ω & \cos Ω & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} \frac{1 - e_{x}^{2} - e_{y}^{2}}{1 + e_{x} \cos >+e_{y}\sin>00} \end{matrix})$

b.坐标分量对e_x的偏导数

$\frac{\partial}{\partial e_{x}} (\begin{matrix} x \\ y \\ z \end{matrix}) = (\begin{matrix} - 2 \sin u \sin λ & 2 \cos u \sin λ & 0 \\ - 2 \cos u \sin λ & - 2 \sin u \sin λ & 0 \\ 0 & 0 & 0 \end{matrix}) (\begin{matrix} 1 & 0 & 0 \\ 0 & \cos i & \sin i \\ 0 & - \sin i & \cos i \end{matrix}) (\begin{matrix} \cos Ω & \sin Ω & 0 \\ - \sin Ω & \cos Ω & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} r \\ 0 \\ 0 \end{matrix})$

$+ (\begin{matrix} \cos u & \sin u & 0 \\ - \sin u & \cos u & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} 1 & 0 & 0 \\ 0 & \cos i & \sin i \\ 0 & - \sin i & \cos i \end{matrix}) (\begin{matrix} \cos Ω & \sin Ω & 0 \\ - \sin Ω & \cos Ω & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} \frac{- 2 a e_{x} - r (\cos u - 2 e_{x} \sin u \sin λ + 2 e_{y} \cos u \sin λ)}{1 + e_{x} \cos u + e_{y} \sin u} \\ 0 \\ 0 \end{matrix})$

c.坐标分量对e_y的偏导数

$\frac{\partial}{\partial e_{y}} (\begin{matrix} x \\ y \\ z \end{matrix}) = (\begin{matrix} 2 \sin u \cos λ & - 2 \cos u \cos λ & 0 \\ 2 \cos u \cos λ & 2 \sin u \cos λ & 0 \\ 0 & 0 & 0 \end{matrix}) (\begin{matrix} 1 & 0 & 0 \\ 0 & \cos i & \sin i \\ 0 & - \sin i & \cos i \end{matrix}) (\begin{matrix} \cos Ω & \sin Ω & 0 \\ - \sin Ω & \cos Ω & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} r \\ 0 \\ 0 \end{matrix})$

$+ (\begin{matrix} \cos u & \sin u & 0 \\ - \sin u & \cos u & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} 1 & 0 & 0 \\ 0 & \cos i & \sin i \\ 0 & - \sin i & \cos i \end{matrix}) (\begin{matrix} \cos Ω & \sin Ω & 0 \\ - \sin Ω & \cos Ω & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} \frac{- 2 a e_{y} - r (\sin u - 2 e_{y} \cos u \cos λ + 2 e_{x} \sin u \cos λ)}{1 + e_{x} \cos u + e_{y} \sin u} \\ 0 \\ 0 \end{matrix})$

d.坐标分量对i的偏导数

$\frac{\partial}{\partial i} (\begin{matrix} x \\ y \\ z \end{matrix}) = (\begin{matrix} \cos u & \sin u & 0 \\ - \sin u & \cos u & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} 1 & 0 & 0 \\ 0 & - \sin i & \cos i \\ 0 & - \cos i & - \sin i \end{matrix}) (\begin{matrix} \cos Ω & \sin Ω & 0 \\ - \sin Ω & \cos Ω & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} r \\ 0 \\ 0 \end{matrix})$

e.坐标分量对Ω的偏导数

$\frac{\partial}{\partial Ω} (\begin{matrix} x \\ y \\ z \end{matrix}) = (\begin{matrix} \cos u & \sin u & 0 \\ - \sin u & \cos u & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} 1 & 0 & 0 \\ 0 & \cos i & \sin i \\ 0 & - \sin i & \cos i \end{matrix}) (\begin{matrix} - \sin Ω & \cos Ω & 0 \\ - \cos Ω & - \sin Ω & 0 \\ 0 & 0 & 0 \end{matrix}) (\begin{matrix} r \\ 0 \\ 0 \end{matrix})$

f.坐标分量对λ的偏导数

$\frac{\partial}{\partial λ} (\begin{matrix} x \\ y \\ z \end{matrix}) = (1 + 2 e_{x} \cos λ + 2 e_{y} \sin λ) (\begin{matrix} - \sin u & \cos u & 0 \\ - \cos u & - \sin u & 0 \\ 0 & 0 & 0 \end{matrix}) (\begin{matrix} 1 & 0 & 0 \\ 0 & \cos i & \sin i \\ 0 & - \sin i & \cos i \end{matrix}) (\begin{matrix} \cos Ω & \sin Ω & 0 \\ - \sin Ω & \cos Ω & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} r \\ 0 \\ 0 \end{matrix})$

$+ (\begin{matrix} \cos u & \sin u & 0 \\ - \sin u & \cos u & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} 1 & 0 & 0 \\ 0 & \cos i & \sin i \\ 0 & - \sin i & \cos i \end{matrix}) (\begin{matrix} \cos Ω & \sin Ω & 0 \\ - \sin Ω & \cos Ω & 0 \\ 0 & 0 & 1 \end{matrix}) (\begin{matrix} \frac{- r (1 + 2 e_{x} \cos λ + 2 e_{y} \sin λ) (e_{y} \cos u - e_{x} \sin u)}{1 + e_{x} \cos u + e_{y} \sin u} \\ 0 \\ 0 \end{matrix})$

将瞬时轨道根数σ_ai＝[a_i，e_xi，e_yi，i_i，Ω_i，λ_i]^T对初始轨道根数σ_a0＝[a₀，e_x0，e_y0，i₀，Ω₀，λ₀]^T的偏导数写成矩阵形式为：

$X (t, σ_{0}) = \frac{\partial σ_{i}}{\partial σ_{0}} = (\begin{matrix} \frac{\partial a_{i}}{\partial a_{0}} & \frac{\partial a_{i}}{\partial e_{x 0}} & \frac{\partial a_{i}}{\partial e_{y 0}} & \frac{\partial a_{i}}{\partial i_{0}} & \frac{\partial a_{i}}{\partial Ω_{0}} & \frac{\partial a_{i}}{\partial λ_{0}} \\ \frac{\partial e_{xi}}{\partial a_{0}} & \frac{\partial e_{xi}}{\partial e_{x 0}} & \frac{\partial e_{xi}}{\partial e_{y 0}} & \frac{\partial e_{xi}}{\partial i_{0}} & \frac{\partial e_{xi}}{\partial Ω_{0}} & \frac{\partial e_{xi}}{\partial λ_{0}} \\ \frac{\partial e_{yi}}{\partial a_{0}} & \frac{\partial e_{yi}}{\partial e_{x 0}} & \frac{\partial e_{yi}}{\partial e_{y 0}} & \frac{\partial e_{yi}}{\partial i_{0}} & \frac{\partial e_{yi}}{\partial Ω_{0}} & \frac{\partial e_{yi}}{\partial λ_{0}} \\ \frac{\partial i_{i}}{\partial a_{0}} & \frac{\partial i_{i}}{\partial e_{x 0}} & \frac{\partial i_{i}}{\partial e_{y 0}} & \frac{\partial i_{i}}{\partial i_{0}} & \frac{\partial i_{i}}{\partial Ω_{0}} & \frac{\partial i_{i}}{\partial λ_{0}} \\ \frac{\partial Ω_{i}}{\partial a_{0}} & \frac{\partial Ω_{i}}{\partial e_{x 0}} & \frac{\partial Ω_{i}}{\partial e_{y 0}} & \frac{\partial Ω_{i}}{\partial i_{0}} & \frac{\partial Ω_{i}}{\partial Ω_{0}} & \frac{\partial Ω_{i}}{\partial λ_{0}} \\ \frac{\partial λ_{i}}{\partial a_{0}} & \frac{\partial λ_{i}}{\partial e_{x 0}} & \frac{\partial λ_{i}}{\partial e_{y 0}} & \frac{\partial λ_{i}}{\partial i_{0}} & \frac{\partial λ_{i}}{\partial Ω_{0}} & \frac{\partial λ_{i}}{\partial λ_{0}} \end{matrix})$

X(t，σ₀)采用数值算法求解矩阵微分方程，即

$(\begin{matrix} \overset{\cdot}{X} (t, σ_{0}) = F (t, σ_{0}) X (t, σ_{0}) \\ X (t_{0}, σ_{0}) = I_{6 \times 6} \end{matrix})$

其中F(t，σ₀)的表达式为

$F (t, σ_{0}) = (\begin{matrix} \frac{\partial}{\partial a} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) & \frac{\partial}{\partial e_{x}} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) & \frac{\partial}{\partial e_{y}} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) & \frac{\partial}{\partial i} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) & \frac{\partial}{\partial Ω} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) & \frac{\partial}{\partial λ} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) \end{matrix})$

其中I_6×6为6维单位方阵。

a.瞬时轨道根数对a的偏导数

$\frac{\partial}{\partial a} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) = (\begin{matrix} 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ - \frac{3}{2} \sqrt{\frac{μ}{a^{5}}} \end{matrix})$

$+ (\begin{matrix} \frac{3 a}{\sqrt{μp}} (e_{x} \sin u - e_{y} \cos u) & \frac{3 a}{\sqrt{μp}} (1 + e_{x} \cos u + e_{y} \sin u) & 0 \\ 0 & 0 & \frac{r \sin u}{2 a \sqrt{μp} \sin i} \\ 0 & 0 & \frac{r \cos u}{\sqrt{μp}} \\ \frac{1}{2 a} \sqrt{\frac{p}{μ}} \sin u & \frac{1}{2 a} \sqrt{\frac{p}{μ}} [(1 + \frac{r}{p}) \cos u + \frac{e_{x} r}{p}] & \frac{1}{2 a} \sqrt{\frac{p}{μ}} \frac{e_{y} r}{p} \frac{\sin u}{\tan i} \\ - \frac{1}{2 a} \sqrt{\frac{p}{μ}} \cos u & \frac{1}{2 a} \sqrt{\frac{p}{μ}} [(1 + \frac{r}{p}) \sin u + \frac{e_{y} r}{p}] & - \frac{1}{2 a} \sqrt{\frac{p}{μ}} \frac{e_{x} r}{p} \frac{\sin u}{\tan i} \\ \frac{1}{2 a} \sqrt{\frac{p}{μ}} [- \frac{1}{2} (e_{x} \cos u + e_{y} \sin u) - \frac{2 r}{p} \sqrt{1 - e_{x}^{2} - e_{y}^{2}}] & \frac{1}{4 a} \sqrt{\frac{p}{μ}} (e_{x} \sin u - e_{y} \cos u) & - \frac{1}{2 a} \sqrt{\frac{p}{μ}} \frac{r}{p} \frac{\sin u}{\tan i} \end{matrix}) (\begin{matrix} f_{r} \\ f_{u} \\ f_{h} \end{matrix})$

$- \frac{4}{a} (\begin{matrix} \frac{{2 a}^{2}}{\sqrt{μp}} (e_{x} \sin u - e_{y} \cos u) & \frac{{2 a}^{2}}{\sqrt{μp}} (1 + e_{x} \cos u + e_{y} \sin u) & 0 \\ 0 & 0 & \frac{r \sin u}{\sqrt{μp} \sin i} \\ 0 & 0 & \frac{r \cos u}{\sqrt{μp}} \\ \sqrt{\frac{p}{μ}} \sin u & \sqrt{\frac{p}{μ}} [(1 + \frac{r}{p}) \cos u + \frac{e_{x} r}{p}] & \sqrt{\frac{p}{μ}} \frac{e_{y} r}{p} \frac{\sin u}{\tan i} \\ - \sqrt{\frac{p}{μ}} \cos u & \sqrt{\frac{p}{μ}} [(1 + \frac{r}{p}) \sin u + \frac{e_{y} r}{p}] & - \sqrt{\frac{p}{μ}} \frac{e_{x} r}{p} \frac{\sin u}{\tan i} \\ \sqrt{\frac{p}{μ}} [- \frac{1}{2} (e_{x} \cos u + e_{y} \sin u) - \frac{2 r}{p} \sqrt{1 - e_{x}^{2} - e_{y}^{2}}] & \frac{1}{2} \sqrt{\frac{p}{μ}} (e_{x} \sin u - e_{y} \cos u) & - \sqrt{\frac{p}{μ}} \frac{r}{p} \frac{\sin u}{\tan i} \end{matrix}) (\begin{matrix} f_{r} \\ f_{u} \\ f_{h} \end{matrix})$

b.瞬时轨道根数对e_x的偏导数

$\frac{\partial}{\partial e_{x}} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) = (\begin{matrix} P_{11} & P_{12} & 0 \\ 0 & 0 & P_{23} \\ 0 & 0 & P_{33} \\ P_{41} & P_{42} & P_{43} \\ P_{51} & P_{52} & P_{53} \\ P_{61} & P_{62} & P_{63} \end{matrix}) (\begin{matrix} f_{r} \\ f_{u} \\ f_{h} \end{matrix}) + Φ (\begin{matrix} \frac{3 J_{2} {Re}^{2} μ}{{2 r}^{4}} [6 \sin^{2} i \sin 2 u \sin λ + (4 \frac{P_{2}}{P_{3}} + \frac{8 e_{x}}{1 - e_{x}^{2} - e_{y}^{2}}) (3 \sin^{2} i \sin^{2} u - 1)] \\ - \frac{3 J_{2} {Re}^{2} μ}{{2 r}^{4}} [4 \sin^{2} i \cos 2 u \sin λ + (4 \frac{P_{2}}{P_{3}} + \frac{8 e_{x}}{1 - e_{x}^{2} - e_{y}^{2}}) \sin^{2} i \sin 2 u] \\ - \frac{3 J_{2} {Re}^{2} μ}{{2 r}^{4}} [2 \sin 2 i \cos u \sin λ + (4 \frac{P_{2}}{P_{3}} + \frac{8 e_{x}}{1 - e_{x}^{2} - e_{y}^{2}}) \sin 2 i \sin u] \end{matrix})$

其中

P₁＝sin u+2e_xcos u sinλ+2e_ysin u sinλ

P₂＝cos u-2e_xsin u sinλ+2e_ycos u sinλ

P₃＝1+e_xcos u+e_ysin u

$P_{11} = \frac{{2 a}^{2}}{\sqrt{μp}} [P_{1} + \frac{e_{x} (e_{x} \sin u - e_{y} \cos u)}{1 - e_{x}^{2} - e_{y}^{2}}]$

$P_{41} = 2 \sqrt{\frac{p}{μ}} \cos u \sin λ - \sqrt{\frac{p}{μ}} \frac{e_{x}}{1 - e_{x}^{2} - e_{y}^{2}} \sin u$

$P_{51} = 2 \sqrt{\frac{p}{μ}} \sin u \sin λ + \sqrt{\frac{p}{μ}} \frac{e_{x}}{1 - e_{x}^{2} - e_{y}^{2}} \cos u$

$P_{61} = \sqrt{\frac{p}{μ}} [- \frac{1}{2} P_{2} - \frac{{2 e}_{x} P_{3} + 2 (1 - e_{x}^{2} - e_{y}^{2}) P_{2}}{\sqrt{1 - e_{x}^{2} - e_{y}^{2}} P_{3}^{2}}] - \sqrt{\frac{p}{μ}} \frac{e_{x}}{1 - e_{x}^{2} - e_{y}^{2}} [- \frac{1}{2} (e_{x} \cos u + e_{y} \sin u) - \frac{2 r}{p} \sqrt{1 - e_{x}^{2} - e_{y}^{2}}]$

$P_{12} = \frac{{2 a}^{2}}{\sqrt{μp}} [P_{2} + \frac{e_{x} P_{3}}{1 - e_{x}^{2} - e_{y}^{2}}]$

$P_{42} = \sqrt{\frac{p}{μ}} [- 2 \sin u \sin λ + \frac{(1 - 2 \sin u \sin λ) P_{3} - (\cos u + e_{x}) P_{2}}{{P_{3}}^{2}}] - \sqrt{\frac{p}{μ}} \frac{e_{x}}{1 - e_{x}^{2} - e_{y}^{2}} [(1 + \frac{r}{p}) \cos u + \frac{e_{x} r}{p}]$

$P_{52} = \sqrt{\frac{p}{μ}} [2 \cos u \sin λ + \frac{2 \cos u \sin λ P_{3} - (\sin u + e_{y}) P_{2}}{{P_{3}}^{2}}] - \sqrt{\frac{p}{μ}} \frac{e_{x}}{1 - e_{x}^{2} - e_{y}^{2}} [(1 + \frac{r}{p}) \sin u + \frac{e_{y} r}{p}]$

$P_{62} = \frac{1}{2} \sqrt{\frac{p}{μ}} [P_{1} - \frac{e_{x} (e_{x} \sin u - e_{y} \cos u)}{1 - e_{x}^{2} - e_{y}^{2}}]$

$P_{23} = \sqrt{\frac{p}{μ}} \frac{1}{\sin i} \frac{2 \cos u \sin λ P_{3} - \sin u P_{2}}{{P_{3}}^{2}} - \sqrt{\frac{p}{μ}} \frac{1}{\sin i} \frac{e_{x}}{1 - e_{x}^{2} - e_{y}^{2}} \frac{\sin u}{P_{3}}$

$P_{33} = \sqrt{\frac{p}{μ}} \frac{- 2 \sin u \sin λ P_{3} - \cos u P_{2}}{{P_{3}}^{2}} - \sqrt{\frac{p}{μ}} \frac{e_{x}}{1 - e_{x}^{2} - e_{y}^{2}} \frac{\cos u}{P_{3}}$

$P_{43} = e_{y} [\sqrt{\frac{p}{μ}} \frac{1}{\tan i} \frac{2 \cos u \sin λ P_{3} - \sin u P_{2}}{{P_{3}}^{2}} - \sqrt{\frac{p}{μ}} \frac{1}{\tan i} \frac{e_{x}}{1 - e_{x}^{2} - e_{y}^{2}} \frac{\sin u}{P_{3}}]$

$P_{53} = \sqrt{\frac{p}{μ}} \frac{1}{\tan i} \frac{(\sin u + 2 e_{x} \cos u \sin λ) P_{3} - e_{x} \sin u P_{2}}{{P_{3}}^{2}} + \sqrt{\frac{p}{μ}} \frac{1}{\tan i} \frac{e_{x}}{1 - e_{x}^{2} - e_{y}^{2}} \frac{e_{x} \sin u}{P_{3}}$

$P_{63} = - \sqrt{\frac{p}{μ}} \frac{1}{\tan i} \frac{2 \cos u \sin λ P_{3} - \sin u P_{2}}{{P_{3}}^{2}} + \sqrt{\frac{p}{μ}} \frac{1}{\tan i} \frac{e_{x}}{1 - e_{x}^{2} - e_{y}^{2}} \frac{\sin u}{P_{3}}$

$Φ = (\begin{matrix} \frac{{2 a}^{2}}{\sqrt{μp}} (e_{x} \sin u - e_{y} \cos u) & \frac{{2 a}^{2}}{\sqrt{μp}} (1 + e_{x} \cos u + e_{y} \sin u) & 0 \\ 0 & 0 & \frac{r \sin u}{\sqrt{μp} \sin i} \\ 0 & 0 & \frac{r \cos u}{\sqrt{μp}} \\ \sqrt{\frac{p}{μ}} \sin u & \sqrt{\frac{p}{μ}} [(1 + \frac{r}{p}) \cos u + \frac{e_{x} r}{p}] & \sqrt{\frac{p}{μ}} \frac{e_{y} r}{p} \frac{\sin u}{\tan i} \\ - \sqrt{\frac{p}{μ}} \cos u & \sqrt{\frac{p}{μ}} [(1 + \frac{r}{p}) \sin u + \frac{e_{y} r}{p}] & - \sqrt{\frac{p}{μ}} \frac{e_{x} r}{p} \frac{\sin u}{\tan i} \\ \sqrt{\frac{p}{μ}} [- \frac{1}{2} (e_{x} \cos u + e_{y} \sin u) - \frac{2 r}{p} \sqrt{1 - e_{x}^{2} - e_{y}^{2}}] & \frac{1}{2} \sqrt{\frac{p}{μ}} (e_{x} \sin u - e_{y} \cos u) & - \sqrt{\frac{p}{μ}} \frac{r}{p} \frac{\sin u}{\tan i} \end{matrix})$

c.瞬时轨道根数对e_y的偏导数

$\frac{\partial}{\partial e_{y}} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) = (\begin{matrix} Q_{11} & Q_{12} & 0 \\ 0 & 0 & Q_{23} \\ 0 & 0 & Q_{33} \\ Q_{41} & Q_{42} & Q_{43} \\ Q_{51} & Q_{52} & Q_{53} \\ Q_{61} & Q_{62} & Q_{63} \end{matrix}) (\begin{matrix} f_{r} \\ f_{u} \\ f_{h} \end{matrix}) + Φ (\begin{matrix} \frac{3 J_{2} {Re}^{2} μ}{{2 r}^{4}} [- 6 \sin^{2} i \sin 2 u \cos λ + (4 \frac{Q_{2}}{P_{3}} + \frac{8 e_{y}}{1 - e_{x}^{2} - e_{y}^{2}}) (3 \sin^{2} i \sin^{2} u - 1)] \\ - \frac{3 J_{2} {Re}^{2} μ}{{2 r}^{4}} [- 4 \sin^{2} i \cos 2 u \cos λ + (4 \frac{Q_{2}}{P_{3}} + \frac{8 e_{y}}{1 - e_{x}^{2} - e_{y}^{2}}) \sin^{2} i \sin 2 u] \\ - \frac{3 J_{2} {Re}^{2} μ}{{2 r}^{4}} [2 \sin 2 i \cos u \cos λ + (4 \frac{Q_{2}}{P_{3}} + \frac{8 e_{y}}{1 - e_{x}^{2} - e_{y}^{2}}) \sin 2 i \sin u] \end{matrix})$

Q₁＝-cos u-2e_xcos u cosλ-2e_ysin u cosλ

Q₂＝sin u+2e_xsin u cosλ-2e_y cos u cosλ

P₃＝1+e_xcos u+e_ysin u

$Q_{11} = \frac{{2 a}^{2}}{\sqrt{μp}} [Q_{1} + \frac{e_{y} (e_{x} \sin u - e_{y} \cos u)}{1 - e_{x}^{2} - e_{y}^{2}}]$

$Q_{41} = - 2 \sqrt{\frac{p}{μ}} \cos u \cos λ - \sqrt{\frac{p}{μ}} \frac{e_{y}}{1 - e_{x}^{2} - e_{y}^{2}} \sin u$

$Q_{51} = - 2 \sqrt{\frac{p}{μ}} \sin u \cos λ + \sqrt{\frac{p}{μ}} \frac{e_{y}}{1 - e_{x}^{2} - e_{y}^{2}} \cos u$

$Q_{61} = \sqrt{\frac{p}{μ}} [- \frac{1}{2} Q_{2} - \frac{2 e_{y} P_{3} + 2 (1 - e_{x}^{2} - e_{y}^{2}) Q_{2}}{\sqrt{1 - e_{x}^{2} - e_{y}^{2}} {P_{3}}^{2}}] - \sqrt{\frac{p}{μ}} \frac{e_{y}}{1 - e_{x}^{2} - e_{y}^{2}} [- \frac{1}{2} (e_{x} \cos u + e_{y} \sin u) - \frac{2 r}{p} \sqrt{1 - e_{x}^{2} - e_{y}^{2}}]$

$Q_{12} = \frac{{2 a}^{2}}{\sqrt{μp}} [Q_{2} + \frac{e_{y} P_{3}}{1 - e_{x}^{2} - e_{y}^{2}}]$

$Q_{42} = \sqrt{\frac{p}{μ}} [2 \sin u \cos λ + \frac{2 \sin u \cos λ P_{3} - (\cos u + e_{x}) Q_{2}}{{P_{3}}^{2}}] - \sqrt{\frac{p}{μ}} \frac{e_{y}}{1 - e_{x}^{2} - e_{y}^{2}} [(1 + \frac{r}{p}) \cos u + \frac{e_{x} r}{p}]$

$Q_{52} = \sqrt{\frac{p}{μ}} [- 2 \cos u \cos λ + \frac{(1 - 2 \cos u \cos λ) P_{3} - (\sin u + e_{y}) Q_{2}}{{P_{3}}^{2}}] - \sqrt{\frac{p}{μ}} \frac{e_{y}}{1 - e_{x}^{2} - e_{y}^{2}} [(1 + \frac{r}{p}) \sin u + \frac{e_{y} r}{p}]$

$Q_{62} = \frac{1}{2} \sqrt{\frac{p}{μ}} [Q_{1} - \frac{e_{y} (e_{x} \sin u - e_{y} \cos u)}{1 - e_{x}^{2} - e_{y}^{2}}]$

$Q_{23} = \sqrt{\frac{p}{μ}} \frac{1}{\sin i} \frac{- 2 \cos u \cos λ P_{3} - \sin u Q_{2}}{{P_{3}}^{2}} - \sqrt{\frac{p}{μ}} \frac{1}{\sin i} \frac{e_{y}}{1 - e_{x}^{2} - e_{y}^{2}} \frac{\sin u}{P_{3}}$

$Q_{33} = \sqrt{\frac{p}{μ}} \frac{2 \sin u \cos λ P_{3} - \cos u Q_{2}}{{P_{3}}^{2}} - \sqrt{\frac{p}{μ}} \frac{e_{y}}{1 - e_{x}^{2} - e_{y}^{2}} \frac{\cos u}{P_{3}}$

$Q_{43} = \sqrt{\frac{p}{μ}} \frac{1}{\tan i} [\frac{(\sin u - 2 e_{y} \cos u \cos λ) P_{3} - e_{y} \sin u Q_{2}}{{P_{3}}^{2}} - \frac{e_{y}}{1 - e_{x}^{2} - e_{y}^{2}} \frac{e_{y} \sin u}{P_{3}}]$

$Q_{53} = e_{x} \sqrt{\frac{p}{μ}} \frac{1}{\tan i} [\frac{- 2 \cos u \cos λ P_{3} - \sin u Q_{2}}{{P_{3}}^{2}} + \frac{e_{y}}{1 - e_{x}^{2} - e_{y}^{2}} \frac{\sin u}{P_{3}}]$

$Q_{63} = - \sqrt{\frac{p}{μ}} \frac{1}{\tan i} \frac{- 2 \cos u \cos λ P_{3} - \sin u Q_{2}}{{P_{3}}^{2}} + \sqrt{\frac{p}{μ}} \frac{1}{\tan i} \frac{e_{y}}{1 - e_{x}^{2} - e_{y}^{2}} \frac{\sin u}{P_{3}}$

d.瞬时轨道根数对i的偏导数

$\frac{\partial}{\partial i} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) = (\begin{matrix} 0 & 0 & 0 \\ 0 & 0 & - \frac{r \sin u \cos i}{\sqrt{μp} \sin^{2} i} \\ 0 & 0 & 0 \\ 0 & 0 & - \sqrt{\frac{p}{μ}} \frac{e_{y} r \sin u}{\sin^{2} i} \\ 0 & 0 & \sqrt{\frac{p}{μ}} \frac{e_{x} r \sin u}{\sin^{2} i} \\ 0 & 0 & \sqrt{\frac{p}{μ}} \frac{r \sin u}{p \sin^{2} i} \end{matrix}) (\begin{matrix} f_{r} \\ f_{u} \\ f_{h} \end{matrix}) + Φ (\begin{matrix} \frac{9 J_{2} {Re}^{2} μ}{{2 r}^{4}} \sin 2 i \sin^{2} u \\ - \frac{3 J_{2} {Re}^{2} μ}{{2 r}^{4}} \sin 2 i \sin 2 u \\ - \frac{6 J_{2} {Re}^{2} μ}{{2 r}^{4}} \cos 2 i \sin u \end{matrix})$

e.瞬时轨道根数对Ω的偏导数

$\frac{\partial}{\partial Ω} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) = (\begin{matrix} 0 \\ 0 \\ 0 \\ 0 \\ 0 \\ 0 \end{matrix})$

f.瞬时轨道根数对λ的偏导数

$\frac{\partial}{\partial λ} (\begin{matrix} \overset{\cdot}{a} \\ \overset{\cdot}{Ω} \\ \overset{\cdot}{i} \\ {\overset{\cdot}{e}}_{x} \\ {\overset{\cdot}{e}}_{y} \\ \overset{\cdot}{λ} \end{matrix}) = (1 + 2 e_{x} \cos λ + 2 e_{y} \sin λ) (\begin{matrix} R_{11} & R_{12} & 0 \\ 0 & 0 & \frac{r (\cos u + e_{x})}{\sqrt{μp} \sin i (1 + e_{x} \cos u + e_{y} \sin u)} \\ 0 & 0 & - \frac{r (\cos u + e_{y})}{\sqrt{μp} (1 + e_{x} \cos u + e_{y} \sin u)} \\ \sqrt{\frac{p}{μ}} \cos u & R_{42} & \sqrt{\frac{p}{μ}} \frac{e_{y} r (\cos u + e_{x})}{p \tan i (1 + e_{x} \cos u + e_{y} \sin u)} \\ \sqrt{\frac{p}{μ}} \sin u & R_{52} & - \sqrt{\frac{p}{μ}} \frac{e_{x} r (\cos u + e_{x})}{p \tan i (1 + e_{x} \cos u + e_{y} \sin u)} \\ R_{61} & R_{62} & - \sqrt{\frac{p}{μ}} \frac{r (\cos u + e_{x})}{p \tan i (1 + e_{x} \cos u + e_{y} \sin u)} \end{matrix}) (\begin{matrix} f_{r} \\ f_{u} \\ f_{h} \end{matrix})$

$+ (1 + {2 e}_{x} \cos λ + {2 e}_{y} \sin λ) Φ (\begin{matrix} \frac{3 J_{2} {Re}^{2} μ}{{2 r}^{4}} [4 \frac{- e_{x} \sin u + e_{y} \cos u}{1 + e_{x} \cos u + e_{y} \sin u} (3 \sin^{2} i \sin^{2} u - 1) + 3 \sin 2 u \sin^{2} i] \\ - \frac{3 J_{2} {Re}^{2} μ}{{2 r}^{4}} [4 \frac{- e_{x} \sin u + e_{y} \cos u}{1 + e_{x} \cos u + e_{y} \sin u} \sin^{2} i \sin 2 u + 2 \cos 2 u \sin^{2} i] \\ - \frac{3 J_{2} {Re}^{2} μ}{{2 r}^{4}} [4 \frac{- e_{x} \sin u + e_{y} \cos u}{1 + e_{x} \cos u + e_{y} \sin u} \sin 2 i \sin u + \cos u \sin 2 i] \end{matrix})$

$R_{11} = \frac{{2 a}^{2}}{\sqrt{μp}} (e_{x} \cos >+e_{y}\sin>)$

$R_{61} = \sqrt{\frac{p}{μ}} (- e_{x} \sin >+e_{y}\cos>)[-\frac{1}{2}+\frac{2 r}{p}\frac{\sqrt{1 - e_{x}^{2} - e_{y}^{2}}}{1 + e_{x} \cos >+e_{y}\sin>}]$

$R_{12} = \frac{{2 a}^{2}}{\sqrt{μp}} ({- e}_{x} \sin u + e_{y} \cos u)$

$R_{42} = \sqrt{\frac{p}{μ}} (\frac{- \sin u - e_{y} - e_{x}^{2} \sin u + e_{x} e_{y} \cos u}{{(1 + e_{x} \cos u + e_{y} \sin u)}^{2}} - \sin u)$

$R_{52} = \sqrt{\frac{p}{μ}} (\frac{\cos u + e_{x} + e_{y}^{2} \cos u - e_{x} e_{y} \sin u}{{(1 + e_{x} \cos u + e_{y} \sin u)}^{2}} + \cos u)$

$R_{62} = \frac{1}{2} \sqrt{\frac{p}{μ}} (e_{x} \cos u + e_{y} \sin u)$

(6)同理，计算卫星b在初始时刻t₀瞬时轨道根数对星间距离计算值的状态转移矩阵：

$\frac{\partial L_{i}}{\partial σ_{b 0}} = \frac{\partial L_{i}}{\partial x_{bi}} \frac{\partial x_{bi}}{\partial σ_{bi}} \frac{\partial σ_{bi}}{\partial σ_{b 0}} + \frac{\partial L_{i}}{\partial y_{bi}} \frac{\partial y_{bi}}{\partial σ_{bi}} \frac{\partial σ_{bi}}{\partial σ_{b 0}} + \frac{\partial L_{i}}{\partial z_{bi}} \frac{\partial z_{bi}}{\partial σ_{bi}} \frac{\partial σ_{bi}}{\partial σ_{b 0}}$

其中，相对距离对各坐标分量的偏导数如下：

$\frac{\partial L_{i}}{\partial x_{bi}} = \frac{x_{bi} - x_{ai}}{L_{i}}$

$\frac{\partial L_{i}}{\partial y_{bi}} = \frac{y_{bi} - y_{ai}}{L_{i}}$

$\frac{\partial L_{i}}{\partial z_{bi}} = \frac{z_{bi} - z_{ai}}{L_{i}}$

坐标分量对轨道根数的偏导数以及瞬时轨道根数对初始轨道根数的偏导数的计算方法与上一步中针对卫星a的计算方法相同。

(7)根据星间相对距离的计算值与测量真值的偏差δL_i，计算卫星a和卫星b在初始时刻t₀的瞬时轨道根数估值的修正量：

将初始时刻瞬时轨道根数的微小改变量Δσ_a0将引起星间距离改变量以矩阵形式表示为：

$(\begin{matrix} Δ L_{0} \\ Δ L_{1} \\ . \\ . \\ . \\ Δ L_{n} \end{matrix}) = (\begin{matrix} \frac{\partial L_{0}}{\partial σ_{a 0}} & \frac{\partial L_{0}}{\partial σ_{b 0}} \\ \frac{\partial L_{1}}{\partial σ_{a 0}} & \frac{\partial L_{1}}{\partial σ_{b 0}} \\ . & . \\ . & . \\ . & . \\ \frac{\partial L_{n}}{\partial σ_{a 0}} & \frac{\partial L_{n}}{\partial σ_{b 0}} \end{matrix}) (\begin{matrix} {Δσ}_{a 0} \\ {Δσ}_{b 0} \end{matrix})$

将(n+1)×12维矩阵M记为

$M = (\begin{matrix} \frac{\partial L_{0}}{\partial σ_{a 0}} & \frac{\partial L_{0}}{\partial σ_{b 0}} \\ \frac{\partial L_{1}}{\partial σ_{a 0}} & \frac{\partial L_{1}}{\partial σ_{b 0}} \\ . & . \\ . & . \\ . & . \\ \frac{\partial L_{n}}{\partial σ_{a 0}} & \frac{\partial L_{n}}{\partial σ_{b 0}} \end{matrix})$

将星间距离改正为真值，即改正量为ΔL_i＝-δL_i；则初始时刻瞬时轨道根数的微小改变量Δσ_a0根据

$- (\begin{matrix} δ L_{0} \\ δ L_{1} \\ . \\ . \\ . \\ δ L_{n} \end{matrix}) = M \cdot (\begin{matrix} {Δσ}_{a 0} \\ {Δσ}_{b 0} \end{matrix})$

进行求解，得到

$(\begin{matrix} {Δσ}_{a 0} \\ {Δσ}_{b 0} \end{matrix}) = - {(M^{T} M)}^{- 1} M^{T} (\begin{matrix} δ L_{0} \\ δ L_{1} \\ . \\ . \\ . \\ δ L_{n} \end{matrix}) \cdot$

(8)将σ_a0+Δσ_a0和σ_b0+Δσ_b0作为卫星a和卫星b在初始时刻t₀瞬时轨道根数估计值，重复上述步骤(1)～(7)的计算过程进行迭代，直至半长轴修正量δa_a0和δa_b0小于定位精度，如0.1m；

(9)根据卫星a和卫星b在初始时刻t₀瞬时轨道根数的迭代结果σ_a0和σ_b0，按照第(2)步的计算方法可以得到任意时刻的瞬时轨道根数σ_a和σ_b，以及相对瞬时轨道根数Δσ＝σ_b-σ_a。

实施例

本实施例所选取卫星a和卫星b在初始时刻的瞬时轨道根数的真值为：

星间距离的真值如图2所示。本实施例的解算输入为以32s为采样间隔，均值为0且方差为3m的星间伪距。

应用本发明方法对4个轨道周期内的星间伪距进行解算，卫星a和卫星b 在初始时刻的瞬时轨道根数的估值如下：

经过10次迭代运算，得到解算值如下：

绝对轨道根数解算结果比较如下：

相对轨道根数解算结果比较如下：

项目真值解算结果偏差半长轴(m) 19.5 19.4 0.1 倾角(°) -0.001 -0.001 1×10^-4 偏心率矢量e_x 0.731×10^-4 0.744×10^-4 1×10^-6 偏心率矢量e_y 0.6404×10^-4 0.6902×10^-4 5×10^-6 升交点赤经(°) 0.895 0.895 1×10^-4 平纬度幅角(°) -0.895 -0.895 1×10^-4

本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。

去获取专利，查看全文>

相似文献

专利
中文文献
外文文献

1. 一种编队飞行卫星绝对和相对轨道确定方法 [P] . 中国专利： CN102322862B . 2013.07.24
2. 一种圆参考轨道下给定边界的卫星初始相对状态确定方法 [P] . 中国专利： CN104200030B . 2017.03.29
3. Absolute angular position determination method for determining an object's position relative to a reference object using an optically sensed information track on a data medium, said track incorporating 2D pixel information [P] . 德国专利： DE10233158A1 . 2004-02-12

机译：用于使用数据介质上的光学感测信息轨道确定物体相对于参考物体的位置的绝对角位置确定方法，所述轨道包含2D像素信息
4. SATELLITE ORBITAL MODELING PROPRIETY DETERMINATION METHOD, LONG TERM PREDICTION ORBITAL DATA SUPPLY METHOD, AND SATELLITE ORBITAL MODELING PROPRIETY DETERMINATION DEVICE [P] . 日本专利： JP2010060489A . 2010-03-18

机译：卫星轨道建模属性确定方法，长期预测轨道数据提供方法和卫星轨道建模属性确定设备
5. Component`s e.g. piston rod, absolute position determining method for motor vehicle, involves generating switching signal by movement of wheel side axle part with vehicle body when reference position of part with respect to body is attained [P] . 法国专利： FR2863353A1 . 2005-06-10

机译：组件的活塞杆，一种用于汽车的绝对位置确定方法，涉及当达到部件相对于车身的基准位置时，通过轮侧轴部件与车身的运动来产生切换信号。