GNSS多普勒测速

GNSS多普勒测速

某一时刻t测站i至卫星j的伪距观测方程为:
(1) ρ ~ i j ( t ) = ρ i j ( t ) + c d t i − c d t j + I i j + T i + ε ρ \begin{aligned} \tilde{\rho}_i^j(t) = \rho_i^j(t)+cdt_i-cdt^j+I_i^j+T_i+\varepsilon_{\rho} \tag{1} \end{aligned} ρ~ij(t)=ρij(t)+cdticdtj+Iij+Ti+ερ(1)式中, ρ ~ i j ( t ) \tilde{\rho}_i^j(t) ρ~ij(t) ρ i j ( t ) \rho_i^j(t) ρij(t)分别为卫星j与接收机i在t时刻的测码伪距和几何距离; d t i dt_i dti d t j dt^j dtj分别为接收机i和卫星j在t时刻的钟差; I i j I_i^j Iij T i T_i Ti分别为电离层和对流层延迟; ε ρ \varepsilon_{\rho} ερ为观测值噪声。对(1)式微分可得:
(2) ρ ~ i j ( t ) ˙ = ρ i j ( t ) ˙ + c d t i ˙ − c d t j ˙ + I ˙ i j + T ˙ i \begin{aligned} \dot{\tilde{\rho}_i^j(t)} = \dot{\rho_i^j(t)}+c\dot{dt_i}-c\dot{dt^j}+\dot{I}_i^j+\dot{T}_i \tag{2} \end{aligned} ρ~ij(t)˙=ρij(t)˙+cdti˙cdtj˙+I˙ij+T˙i(2)由此,可以得出多普勒测速观测方程,即:
(3) λ D i j = ρ i j ( t ) ˙ + c d t i ˙ − c d t j ˙ + I ˙ i j + T ˙ i + ε D \begin{aligned} \lambda D_i^j = \dot{\rho_i^j(t)}+c\dot{dt_i}-c\dot{dt^j}+\dot{I}_i^j+\dot{T}_i+\varepsilon_D \tag{3} \end{aligned} λDij=ρij(t)˙+cdti˙cdtj˙+I˙ij+T˙i+εD(3)其中,
ρ i j ( t ) ˙ = [ l i j ( t ) m i j ( t ) n i j ( t ) ] [ x ˙ j ( t ) − x ˙ i ( t ) y ˙ j ( t ) − y ˙ i ( t ) z ˙ j ( t ) − z ˙ i ( t ) ] \begin{aligned} \dot{\rho_i^j(t)} = \begin{bmatrix} l_i^j(t) & m_i^j(t) & n_i^j(t) \end{bmatrix}\begin{bmatrix} \dot{x}^j(t) - \dot{x}_i(t) \\ \dot{y}^j(t) - \dot{y}_i(t) \\ \dot{z}^j(t) - \dot{z}_i(t) \end{bmatrix} \end{aligned} ρij(t)˙=[lij(t)mij(t)nij(t)]x˙j(t)x˙i(t)y˙j(t)y˙i(t)z˙j(t)z˙i(t)式中,各“.”项为各相应变量的时间变化率; [ l i j ( t ) m i j ( t ) n i j ( t ) ] \begin{bmatrix}l_i^j(t) & m_i^j(t) & n_i^j(t) \end{bmatrix} [lij(t)mij(t)nij(t)]为测站i与卫星j在t时刻的方向余弦; ε D \varepsilon_D εD为观测噪声。

你可能感兴趣的:(GNSS,卫星定位)