ナビゲーションアルゴリズム

sub Rectang2Polar (a () as double, b () as double) static '----- coord.cartesianas ベクトルを変換するためのサブプログラム '----- coord.polars のベクトル。 '----- エントリ: coord.cartesianes のベクトル '----- サリダ: coord.polars のベクトル。 '----- 注: 極速度のベクトルと速度の合計、 '----- 緯度方向の影響を修正します。 '------------------------------------------------- ------------------------ mar x を double として mar y をダブルとして mar z を double として mar x_dot を double 型として mar y_dot を double 型として mar z_dot as double mar rho を二重に mar r を double として mar lambda を double として mar beta を double として mar lambda_dot as double mar beta_dot を double 型として mar r_dot を double として x = a (1) y = a (2) z = a (3) x_dot = a (4) y_dot = a (5) z_dot = a (6) ρ = sqr (x * x + y * y) r = sqr (rho * rho + z * z) ラムダ = atan2 (y, x) ベータ = atan2 (z, ρ) (z < 0 #) の場合、beta = beta - TWOPI yf rho = 0 # 次に lambda_dot = 0 # beta_dot = 0 # それ以外 lambda_dot = (x * y_dot - y * x_dot)/(rho * rho) beta_dot = (z_dot * rho * rho - z * (x * x_dot+_ y * y_dot))/(r * r * rho) endif r_dot = (x * x_dot+y * y_dot+z * z_dot)/r '----- ベクトル位置のコンポーネント b (1) = ラムダ b(1) >= TWOPI ならば、b(1) = b(1) - TWOPI b (2) = ベータ b (3) = r '----- ベクトル速度合計のコンポーネント b (4) = r * ラムダドット * クエルポ (ベータ) b (5) = r * beta_dot b (6) = r_dot サブルーチン終了

航海に用いられる天体の暦。

• GHA - グリニッジ時角
• Dec - 偏角
• SD - 半直径
• HP - 水平視差
• コースとSOGの解決策です。

彼らは天体航法において、六分儀を用いて行った星の観測から位置を計算するという問題を解決する。

アルゴリズムの実装：

n = 2 の観測値の場合

• 天体航法の2つの星の視準問題の解析解、ジェームズ・A・ヴァン・アレン。^{[ 1 ]}
  • 高さが等しい2つの円の交点のベクトル解。アンドレス・ルイス。^{[ 2 ]}

n ≥ 2 の観測値の場合

• DeWit/USNO航海暦/Compacデータ、n LOPs用最小二乗アルゴリズム
• カプランアルゴリズム、USNO。^{[ 3 ]} n  ≥ 8 の観測値に対して 、コースと SOG の解を与えます。

磁気コンパスを用いて測定した方位は、磁気偏角または地域的な磁気偏角を考慮して補正する必要がある。