[发明专利]一种确定整周模糊度的方法及装置

权利要求书

1.一种确定整周模糊度的方法，其特征在于，所述方法包括：

建立载波相位双差观测方程，所述方程中包括待确定的模糊度；

根据最小二乘算法，确定所述载波相位双差观测方程中模糊度的浮点解矩阵以及所述浮点解矩阵对应的方差-协方差矩阵；

对所述方差-协方差矩阵进行矩阵分解，得到包含下三角矩阵和对角矩阵的分解结果，并对分解后的矩阵进行模糊度去相关处理，得到变换矩阵和经过所述变换矩阵变换后的模糊度整数矩阵；

根据所述模糊度整数矩阵，确定模糊度的第一组模糊度候选解，并对所述第一组模糊度候选解中的各个模糊度在当前搜索空间中进行搜索取整，得到第二组模糊度候选解；其中，所述当前搜索空间的数值为第一组模糊度候选解的范数值和初始搜索空间的数值中的较小者；

根据所述变换矩阵对范数值最小的两组候选解进行反变换，得到最终的整周模糊度。

2.根据权利要求1所述的确定整周模糊度的方法，其特征在于，按照下述公式建立所述载波相位双差观测方程：

其中，表示第u个流动站与第r个参考站之间的载波相位双差观测值，λ表示载波相位波长，表示第u个流动站与第r个参考站之间真距的双差值，表示第u个流动站与第r个参考站之间的双差模糊度值，表示第u个流动站与第r个参考站之间的双差噪声误差值。

3.根据权利要求1所述的确定整周模糊度的方法，其特征在于，对分解后的矩阵进行模糊度去相关处理包括：

对所述分解结果中的下三角矩阵进行初等列变换，以使得经过初等列变换的下三角矩阵中非对角元素中的数值减小至0.5以下。

4.根据权利要求1或3所述的确定整周模糊度的方法，其特征在于，按照下述公式确定模糊度整数矩阵：

其中，表示所述模糊度整数矩阵，Z^T表示所述变换矩阵Z的转置，a表示所述浮点解矩阵。

5.根据权利要求1所述的确定整周模糊度的方法，其特征在于，所述初始搜索空间的数值按照下述公式确定：

其中，χ²表示所述初始搜索空间的数值，fac表示预设放大系数，P表示模糊度候选解的组数，Q表示所述方差-协方差矩阵，n表示模糊度的维数。

6.根据权利要求5所述的确定整周模糊度的方法，其特征在于，按照下述公式确定所述第一组模糊度候选解的范数值：

其中，DIST表示所述第一组模糊度候选解的范数值，z_i表示所述第一组模糊度候选解中的第i个候选解，表示估算精度排在第i位的估算整周模糊度值，d_i表示所述对角矩阵中第i行第i列的元素。

7.根据权利要求6所述的确定整周模糊度的方法，其特征在于，对所述第一组模糊度候选解中的各个模糊度在当前搜索空间中进行搜索取整包括：

根据所述第一组模糊度候选解中的第i个候选解与估算精度排在第i位的估算整周模糊度值之间的差值，确定待增加的整数值，并将所述第i个候选解与所述待增加的整数值的和作为搜索取整后的候选解，其中，按照下述公式确定所述待增加的整数值：

其中，Δ(i)表示第i个候选解对应的待增加的整数值。

8.根据权利要求1所述的确定整周模糊度的方法，其特征在于，按照下述公式确定最终的整周模糊度：

z'＝(Z^T)^-1z

其中，z'表示最终的整周模糊度，Z表示所述变换矩阵，z表示范数值最小的两组候选解。