[发明专利]一种特殊测量体制下的多测速数据处理方法

摘要

本发明提供了一种特殊测量体制下的多测速数据处理方法，首先针对一主两副的测量体制，建立测量方程组，然后进行模型设计，求取目标速度，最后计算目标速度的精度。本发明不仅远远满足了轨道数据处理精度要求，而且计算结果与GNSS外测数据处理结果的精度量级相当，明显提高了轨道数据处理精度。

主权项

一种特殊测量体制下的多测速数据处理方法，其特征在于包括下述步骤：步骤一，设分别为主站发送、接收距离变化率，和分别为两个副站的接收距离变化率，为主站测速数据，和分别为两个副站的测速数据，则建立测量方程组$\left\{\begin{array}{c}{\stackrel{\·}{s}}_1={\stackrel{\·}{r}}_0+{\stackrel{\·}{r}}_1\\ {\stackrel{\·}{s}}_2={\stackrel{\·}{r}}_0+{\stackrel{\·}{r}}_2\\ {\stackrel{\·}{s}}_3={\stackrel{\·}{r}}_0+{\stackrel{\·}{r}}_3\end{array}\right.,$式中，${\stackrel{\·}{r}}_1=\frac{\stackrel{\·}{x}(x-x_i)+\stackrel{\·}{y}(y-y_i)+\stackrel{\·}{z}(z-z_i)}{r_i},$i＝0,1,2,3代表测站，$r_i=\sqrt{{(x-x_i)}^2+{(y-y_i)}^2+{(z-z_i)}^2},$x_i,y_i,z_i为测站站址；x,y,z为目标位置分量，为目标速度分量；步骤二，设主站发送、接收距离变化率相同，则目标速度$\left[\begin{array}{c}\stackrel{\·}{x}\\ \stackrel{\·}{y}\\ \stackrel{\·}{z}\end{array}\right]={\left[\begin{array}{ccc}x-x_1& y-y_1& z-z_1\\ x-x_2& y-y_2& z-z_2\\ x-x_3& y-y_2& z-z_3\end{array}\right]}^{-1}\left[\begin{array}{c}\frac{{\stackrel{\·}{s}}_1{r}_1}{2}\\ ({\stackrel{\·}{s}}_2-\frac{{\stackrel{\·}{s}}_1}{2})r_2\\ ({\stackrel{\·}{s}}_3-\frac{{\stackrel{\·}{s}}_1}{2})r_3\end{array}\right];$步骤三，将步骤二的公式写成$\left[\begin{array}{c}\stackrel{\·}{x}\\ \stackrel{\·}{y}\\ \stackrel{\·}{z}\end{array}\right]={\left[\begin{array}{ccc}a& b& c\\ d& e& f\\ g& h& i\end{array}\right]}^{-1}\left[\begin{array}{c}\frac{{\stackrel{\·}{s}}_1{r}_1}{2}\\ ({\stackrel{\·}{s}}_2-\frac{{\stackrel{\·}{s}}_1}{2})r_2\\ ({\stackrel{\·}{s}}_3-\frac{{\stackrel{\·}{s}}_1}{2})r_3\end{array}\right],$则目标速度的精度为$\left[\begin{array}{c}{\mathrm{\σ}}_{\stackrel{\·}{x}}\\ {\mathrm{\σ}}_{\stackrel{\·}{y}}\\ {\mathrm{\σ}}_{\stackrel{\·}{z}}\end{array}\right]=\left[\begin{array}{c}{[{(\frac{\∂\stackrel{\·}{x}}{\∂{\stackrel{\·}{s}}_1}\·{\mathrm{\σ}}_{{\stackrel{\·}{s}}_1})}^2+{(\frac{\∂\stackrel{\·}{x}}{{\∂\stackrel{\·}{s}}_2}\·{\mathrm{\σ}}_{{\stackrel{\·}{s}}_2})}^2+{(\frac{\∂\stackrel{\·}{x}}{\∂{\stackrel{\·}{s}}_3}\·{\mathrm{\σ}}_{{\stackrel{\·}{s}}_3})}^2]}^{\frac{1}{2}}\\ {[{(\frac{\∂\stackrel{\·}{y}}{\∂{\stackrel{\·}{s}}_1}\·{\mathrm{\σ}}_{{\stackrel{\·}{s}}_1})}^2+{(\frac{\∂\stackrel{\·}{y}}{\∂{\stackrel{\·}{s}}_2}\·{\mathrm{\σ}}_{{\stackrel{\·}{s}}_2})}^2+{(\frac{\∂\stackrel{\·}{y}}{\∂{\stackrel{\·}{s}}_3}\·{\mathrm{\σ}}_{{\stackrel{\·}{s}}_3})}^2]}^{\frac{1}{2}}\\ {[{(\frac{\∂\stackrel{\·}{z}}{\∂{\stackrel{\·}{s}}_1}\·{\mathrm{\σ}}_{{\stackrel{\·}{s}}_1})}^2+(\frac{\∂\stackrel{\·}{z}}{\∂{\stackrel{\·}{s}}_2}\·{\mathrm{\σ}}_{{\stackrel{\·}{s}}_2})+{(\frac{\∂\stackrel{\·}{z}}{\∂{\stackrel{\·}{s}}_3}\·{\mathrm{\σ}}_{{\stackrel{\·}{s}}_3})}^2]}^{\frac{1}{2}}\end{array}\right],$其中，分别为的精度。