[发明专利]基于奇偶矢量法的双星故障识别方法

权利要求书

1.基于奇偶矢量法的双星故障识别方法，包括步骤：

(1)、从现有的n颗卫星中剔除卫星i，n为卫星定位系统中接收机可见卫星数量，n≥7；

(2)、对剩余的(n-1)颗卫星根据故障检测判决量和故障检测门限σ_T，(n-1)进行故障检测，若则认为剩余卫星中无故障卫星存在，检测结束；否则，用奇偶矢量法进行故障卫星识别，根据故障识别判决量r_i的最大值r_j＝max(r₁，r₂…,r_n)，将卫星j剔除；

(3)、对剩余的(n-2)颗卫星根据故障检测判决量和故障检测门限σ_T，(n-2)进行故障检测，若则确认卫星i无故障，令无故障卫星数目加一，进入步骤(4)；否则，记录剔除的卫星组合(i_k,j_k)并用最小二乘法计算接收机的位置和钟差x_k，转到步骤(5)；

(4)、判断当前所得的无故障卫星数目是否达到4颗，若达到4颗则解算出接收机的粗略位置和钟差并从剩余(n-4)颗卫星中识别出故障卫星，重新解算接收机位置和钟差，检测结束；否则，进入步骤(5)；

(5)、令i＝i+1，若i＞n，则无故障卫星数目小于4颗，此时，从剔除的卫星组合中以伪距残差统计量λ_k的最小值λ_min＝min(λ₁,λ₂,…,λ_m)所对应的卫星组合作为故障卫星，并计算出与该卫星组合相对应的接收机位置和钟差作为定位结果，检测结束；否则返回步骤(1)。

2.根据权利要求1所述的基于奇偶矢量法的双星故障识别方法，其特征在于，故障检测判决量由下式得到：

σ^n-1=SSE(n-1)-4]]>

其中，SSE为伪距残差平方和；

故障检测门限σ_T，(n-1)由下式得到：

σT,(n-1)=σ0Tn-1(n-1)-4Tn-1=icdf(χ2(1-α,(n-1)-4))]]>

其中，icdf为计算卡方分布的(1-α)分位数；

故障识别判决量r_i由下式得到：

ri=||pQp.i||||Qp.i||.]]>