其坐标系的坐标原点 O n 可根据需要选在任意点, O n Y n 与过原点的地心矢径 r 的方向一致, O n X n 轴指向北极且与 O n Y n 轴垂直, O n Z n 轴指向东方且与 O n X n 、 O n Y n 成右手系(见图2.8-1)。
假定某时刻火箭位于空间点 K ,且以点 K 为原点作北天东坐标系 K-X n Y n Z n 。下面基于定义推导惯性坐标系 O-XYZ 至北天东坐标系 K- X n Y n Z n 的方向余弦矩阵 。
图2.8-1 北天东坐标系与惯性坐标系的关系
如图2.8-1所示,不妨记惯性坐标系原点 O 地心矢径为 R 0 、北天东坐标系原点 K 地心矢径为 r K 、 O 至点 K 矢径为 ρ ,则由式(2.4-5)可知点 O 地心矢径在惯性坐标系投影形式( R 0 x R 0 y R 0 z ) T 。
地心 O e 至点 K 矢量 r K = r 0 + ρ ,其分量形式为
式中, x K 、 y K 、 z K 为 K 点在 O-XYZ 中坐标。
设 e X 、 e Y 、 e Z 为坐标系 O-XYZ 三轴单位矢量, 、 、 为坐标系 K - x n y n z n 三轴单位矢量, e r 为 r K 单位矢量, e ω 为地球自转角速度 ω e 单位矢量,且 e r 、 e ω 可写为
由图2.5-2所示及式(2.5-8),式中 、 、 易知为
又由图2.8-1所示可知
且由矢量点积、叉积运算可知, ,考虑到 ϕ K ∈[-90°,90°]及式(2.8-2),则有
将式(2.8-2)代入式(2.8-6)得
式中
又由矢量积公式 ( a × b )× c = b ( a · c )- a ( b · c ) 、 e r · e r =1 、 e Ω · e r =sin ϕ K 知
考虑到 ϕ K ∈[-90°,90°]及式(2.8-9),有
则将式(2.8-3)、式(2.8-10)、式(2.8-11)代入式(2.8-7)得
式中
由式(2.3-8)可知
则将式(2.8-2)、式(2.8-8)、式(2.8-12)代入式(2.8-14)得发射惯性系至北天东坐标系的方向余弦矩阵:
发射坐标系 O f -X f Y f Z f 至北天东坐标系 K-X n Y n Z n 的方向余弦矩阵 形式上与 完全一致,只是式(2.8-1)中的 x K 、 y K 、 z K 需要使用 K 点地心矢径 r K 在发射坐标系下投影。