[发明专利]基于北斗四频信号的中长基线模糊度解算方法

权利要求书

1.一种基于北斗四频信号的中长基线模糊度解算方法，其特征在于，包括以下步骤：

S1：利用四频组合伪距和四频组合载波计算出两个宽巷模糊度；

S2：利用两个宽巷载波相位的双差观测量及相应模糊度，计算消电离层伪距和B1频点的双差电离层延迟；

S3：根据消电离层伪距和B1频点的双差电离层延迟，计算与第四频率相关的两个模糊度；

S4：根据S1～S3计算出的四个模糊度计算北斗四频载波的独立模糊度。

2.如权利要求1所述的基于北斗四频信号的中长基线模糊度解算方法，其特征在于，所述步骤S1具体包括：

S1.1：构造四频组合伪距和四频组合载波；

北斗伪距、载波相位的观测方程分别为：

Pi=ρ(ts,tr)+C(dtr-dts)+dorb+dtrop+Kfi2+MPi+ϵPi---(1)]]>

式中，下标i表示载波B_i和S相应参数，i＝1、2、3；Pi为载波B_i和S相应的伪距观测量，单位为米；Φ_i，分别为载波B_i和S的相位观测量，单位分别为米、周；λ_i为载波B_i和S的波长，单位为米；N_i为载波B_i和S的整周模糊度，单位为周；dt_s、dt_r分别为卫星钟差、接收机钟差，单位为秒；C为光速，单位为米/秒；为载波B_i和S相应的电离层延迟，单位为米；f_i为载波B_i和S的频率，单位为赫兹；d_trop为对流层的延迟量，单位为米；M_Pi、M_Φi分别为B_i和S频点上伪距、载波相位的多路径效应，单位为米；ε_Pi、ε_Φi分别为伪距、载波相位的观测噪声，单位为米；ρ为卫星至接收机天线的几何距离，单位为米；

以周为单位，四频组合载波的一般形式为：

以米为单位，四频组合载波的一般形式为：

Φi,j,k,m=i·f1·Φ1+j·f2·Φ2+k·f3·Φ3+m·f4·Φ4i·f1+j·f2+k·f3+m·f4---(4)]]>

组合载波的模糊度、频率和波长分别为：

N_{i，j，k，m}＝iN₁+jN₂+kN₃+mN₄                                        (5)

f_{i，j，k，m}＝if₁+jf₂+kf₃+mf₄                                        (6)

λi,j,k,m=Cfi,j,k=λ1λ2λ3λ4i·λ2λ3λ4+j·λ1λ3λ4+kλ1λ2λ4+mλ1λ2λ3---(7)]]>

四频组合伪距的一般形式为：

Pa,b,c,d=a·f1·P1+b·f2·P2+c·f3·P3+d·f4·P4a·f1+b·f2+c·f3+d·f4---(8)]]>

四频组合伪距和四频组合载波的观测方程分别表示为：

Pa,b,c,d=ρ+C(dtr-dts)+dorb+dtrop+βa,b,c,d·Kf12+ϵPa,b,c,d---(9)]]>

式中，为B1频点电离层延迟，β_{a，b，c，d}和β_{i，j，k，m}分别为组合伪距、组合载波的电离层系数：

βa,b,c,d=(af1+bf2+cf3+df4)·f12fa,b,c,d---(11)]]>

βi,j,k,m=(if1+jf2+kf3+mf4)·f12fi,j,k,m---(12)]]>

和分别为组合伪距、组合载波的观测噪声；

根据伪距组合系数(a，b，c，d)和载波组合系数(i，j，k，m)，计算出四频组合伪距和四频组合载波；

S1.2：构造组合伪距和组合载波的双差观测量，方程如下：

▿ΔPa,b,c,d=▿Δρ(ts,tr)+▿Δdorb+▿Δdtrop+βa,b,c,d·▿ΔKf12+▿ΔMP+▿ΔϵPa,b,c,d---(13)]]>

S1.3：利用组合伪距和组合载波的双差观测量计算宽巷模糊度，计算公式为：

▿ΔNi,j,k,m=int[▿ΔΦi,j,k,m-▿ΔPa,b,c,dλi,j,k,m+βi,j,k,m+βa,b,c,dλi,j,k,m·▿ΔKf12-▿ΔϵΦi,j,k,m-▿ΔϵPa,b,c,dλi,j,k,m]---(15)]]>

式中，表示双差算符，int[·]表示四舍五入。