一种基于伪距的陆基长波定位解算方法与流程

技术特征：

1.一种基于伪距的陆基长波定位解算方法，其特征在于包括下述步骤：

步骤1，接收三个发射台的授时电文并进行解析和跟踪，得到每个发射台的信号标志脉冲相对于系统时钟的发射时刻TOT_i以及每个发射台的信号标志脉冲相对于接收机时钟的到达时刻TOA_i，计算出接收点到各发射台的伪距r_i＝(TOA_i-TOT_i)*c，其中，c为电磁波传播速度，i＝1、2、3；

步骤2，设接收点P未知的球面经纬度为发射台的已知球面经纬度为接收点到各发射台的未知大地距离为ρ_i，接收机与发射台的钟差为δt，则ρ_i＝aσ_0i+aδ_si，其中，a为地球的椭球形模型长半轴的长度，aσ_0i为半径等于a的地球球形模型上接收点到发射台的球面距离，aδ_si为对球面距离的修正量；

建立基于椭球面大地距离的伪距方程组

步骤3，确定发射台几何关系，约定三个发射台中的一个为台2，另外两个发射台分别为台1和台3，则台1的基线围绕台2沿顺时针方向旋转至台3的基线，转角小于π；

步骤4，计算接收点的位置和钟差δt(σ₁*a-r₁)/c，其中，σ_i为接收点到各发射台的球面角距，γ₃为台3与正北方向的夹角，θ为台3与接收点的夹角；将作为位置初值δt作为钟差初值δt₀；

步骤5，对基于椭球面大地距离的伪距方程组进行泰勒展开，得到

忽略高次项并将常数项移项，写成矩阵形式A·X＝B；

其中，

求出矩阵A中的偏导数后得到其解向量X，用新值代替旧值，即用代替用Δλ+λ₀代替λ₀，用Δδt+δt₀代替δt₀，进行下一次迭代运算，直到且Δλ|<ε为止，ε是预先设定的收敛门限。