一、計算流程
1.讀取N檔案（導航電文），將所有星歷放到一個串列（陣列）ephlst 中
2.讀取O檔案（觀測值檔案），先讀取一個歷元觀測值epoch
3.資料預處理和程式初始化
（1）根據epoch中衛星號和歷元時刻Tr在星歷中查找相應的衛星星歷
（2）確定測站的概略位置和接識訓鐘差初值
4.計算衛星在瞬時地球坐標系中的位置
（1）計算衛星運行的平均角速度n
（2）計算歸化時間
（3）計算觀測時刻的平近點角
（4）計算偏近點角
（5）計算真近點角
（6）計算升交點角距
； w—近地點角距
（7）計算攝動改正項
（8）計算經過攝動改正的升交角距uk，矢徑rk，軌道偏角ik，
（9）計算衛星在軌道平面坐標系中的坐標
（10）計算觀測時刻升交點經度
（11）計算衛星在瞬時地球坐標系中的坐標
5.計算衛星si信號發射的概略時刻
（1）選擇epoch中的一顆衛星觀測值，設其偽距為psi
（2）利用上面的衛星坐標和測站概略坐標計算幾何距離Rsi
（3）再根據偽距和幾何距離計算衛星信號的傳播時間，利用迭代法求出衛星信號的傳播時間tsi和發射時刻Tsi
（4）其中迭代條件為
6.求衛星si方向余弦
7求衛星si在觀測方程中的余數項
8.選擇epoch中下一顆衛星sj觀測值，設其偽距為psj
9.重復5—8步驟，計算每顆衛星的系數和余數項
10.組成誤差方程，以
為未知數進行求解，方程為：
11.求解法方程
求出定位結果
12.輸出該歷元定位結果
13.繼續讀取下一個歷元并進行上述所有步驟計算，直到全部歷元計算完畢
為止 ，且下一歷元坐標初始值=上一歷元的定位結果
GPS單點定位程式及演算法流程
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/394110.html
標籤：其他