[发明专利]一种实时北斗精密相对定位的方法

权利要求书

1.一种实时北斗精密相对定位的方法，其特征在于，包括以下步骤：

步骤一，整周模糊度参数初始化：利用基准站和移动站的同步北斗观测数据初始化载波相位整周模糊度参数；

步骤二，周跳探测：对站星间的非同步双差载波相位进行周跳探测，对发生周跳的卫星进行标记；

步骤三，伪距单点定位：利用基准站、移动站的北斗伪距观测信息，分别进行基准站和移动站伪距单点定位，获得基准站和移动站的绝对位置，获得基准站测量时的北斗卫星信号发送时刻的卫星坐标及钟差和移动站测量时的北斗卫星信号发送时刻的卫星坐标及钟差；

步骤四，基准站与移动站的观测历元差判断：计算基准站的观测历元与移动站的观测历元之差，如果差值大于15s，则返回步骤一，否则进入步骤五；

步骤五，周跳修复：采用附加模糊度参数方法修复发生周跳的卫星的整周模糊度，即利用未发生周跳的卫星的整周模糊度计算发生周跳的卫星的整周模糊度；

步骤六，求解基准站与移动站的相对位置：利用基准站和移动站的非差非同步载波相位观测数据和卫星坐标以及钟差，建立非同步载波相位双差相对定位观测方程，计算基准站与移动站的相对位置。

2.如权利要求1所述的一种实时北斗精密相对定位的方法，其特征在于，所述步骤六中建立非同步载波相位双差相对定位观测方程，计算基准站与移动站的相对位置，具体过程为：

(S1)建立载波相位双差观测方程组：利用基准站和移动站的非差载波相位观测值及广播星历信息，建立双差载波相位相对定位观测方程；设基准站A在观测时刻t₀，测得i号卫星的非差载波相位观测值为和j号卫星的非差载波相位观测值为移动站B在观测时刻t₁，测得i号卫星的非差载波相位观测值为和j号卫星的非差载波相位观测值为i,j为整数；

定义：

(·)ABij(t0,t1)=((·)Bj(t1)-(·)Aj(t0))-((·)Bi(t1)-(·)Ai(t0))(·)ABij(T0,T1)=((·)Bj( T1j)-(·)Aj(T0j))-((·)Bi(T1 i)-(·)Ai(T0i))]]>

其中，上标i和j代表卫星号，i号卫星为参考星，下标A和B代表基准站和移动站；T代表卫星信号发送时刻，如表示在观测时刻t₀时，对应i号卫星的信号发送时刻，单位为秒；

则非同步载波相位双差观测方程为：

ΦABij(t0,t1)=ρABij(T0,T1)+λNABij(t0,t1)-cδtABij(T0,T1)+ϵABij(t0,t1)---(1)]]>

其中，为不同时刻的双差载波相位观测值，单位为米；为不同时刻的双差几何距离，单位为米；为不同时刻的双差整周模糊度，单位为周；λ为某个频点的对应波长，单位为米；为不同时刻的双差卫星钟差，单位为秒；c为光速，单位为米/秒；为不同时刻的双差噪声，单位为米；

式(1)中，基站与移动站的相对位置参数隐含于双差几何距离量中，其具体为表达式为：

ρABij(T0,T1)=(||rB(t1)-rj(T1 j)||-||rA(t0)-rj(T0j)||)-(||rB(t1)-ri( T1 i)||-||rA(t0)-ri(T0i)||)---(2)]]>

式(2)中，和为j号卫星在信号发送时刻T₁和T₀的三维位置坐标，单位为米；和为i号卫星在信号发送时刻T₁和T₀的三维位置坐标，单位为米；||.||表示卫星到测站之间的几何距离，单位为米；r_A(t₀)为基准站在观测时刻t₀的三维位置坐标；r_B(t₁)为移动站在观测时刻t₁的三维位置坐标，单位为米；

设基准站和移动站共同观测了m颗卫星，则有建立(m-1)个双差载波相位观测方程：

ΦABi1(t0,t1)=ρABi1(T0,T1)+λNABi1(t0,t1)-cδtABi1(T0,T1)+ϵABi1(t0,t1)ΦABi2(t0,t1)=ρABi2(T0,T1)+λNABi2(t0,t1)-cδtABi2(T0,T1)+ϵABi2(t0,t1)···ΦABik(t0,t1)=ρABik(T0,T1)+λNABik(t0,t1)-cδtABik(T0,T1)+ϵABik(t0,t1)---(3)]]>

式中，k表示卫星号，k＝1,2,…，j,…m-1；k≠i；k,m为整数；

(S2)相对定位观测方程组线性化：根据泰勒展开原理，对非线性相对定位观测方程组进行线性化；

在式(2)中，卫星的坐标以实时播发的广播星历参数计算；在基准站位置线性泰勒展开至一阶项，式(2)变为

ρABij(T0,T1)=ρ0j-(uBj-uBi)rABρ0j=(||rA(t0)-rj( T1j)||-||rA(t0)-rj(T0j)||)-(||rA(t0)-ri(T1 i)||-||rA(t0)-ri(T0i)||)---(4)]]>

其中，u代表单位矢量，表示移动站B到j号卫星的单位矢量；表示移动站B到i号卫星的单位矢量；表示j号卫星的双差几何距离泰勒展开的零阶项；r_AB为基准站到移动站的基线矢量；

将式(4)代入式(3)，则线性化的观测方程组为：

ΦABi1(t0,t1)-ρ01-λNABi1(t0,t1)+cδtABi1(T0,T1)=-(uB1-uBi)rAB+ϵABi1(t0,t1)ΦABi2(t0,t1)-ρ02-λNABi2(t0,t1)+cδtABi2(T0,T1)=-(uB2-uBi)rAB+ϵABi2(t0,t1)···ΦABik(t0,t1)-ρ0k-λNABik(t0,t1)+cδtABik(T0,T1)=-(uBk-uBi)rAB+ϵABik(t0,t1)---(5)]]>

式中，k表示卫星号，k＝1,2,…,j,…m-1；k≠i；k,m为整数；

式(5)中，卫星钟差双差根据实时广播星历参数计算；

(S3)线性观测方程组求解：根据最小二乘估计，即可求得基准站到移动站的基线矢量；

记：

Y=ΦABi1(t0,t1)-ρ01-λNABi1(t0,t1)+cδtABi1(T0,T1)ΦABi2(t0,t1)-ρ02-λNABi2(t0,t1)+cδtABi2(T0,T1)···ΦABik(t0,t1)-ρ0k-λNABik(t0,t1)+cδtABik(T0,T1),Q=-(uB1-uBi)-(uB2-uBi)···-(uBk-uBi)---(6)]]>

式中，k表示卫星号，k＝1,2,…,j,…m-1；k≠i；k,m为整数；

根据最小二乘估计原理，求解式(5)，得到：

r_AB＝(Q^TQ)^-1Q^TY          (7)

r_AB即可确定基准站与移动站之间的相对位置。